一种量化流域水文生态系统服务间权衡与协同关系的计算方法技术方案

本发明专利技术涉及一种量化流域水文生态系统服务间权衡与协同关系的计算方法，包括：采用流域水文模型将整个流域进行划分，得到水文响应单元空间分布图层并对其进行处理；统计各水文响应单元相应位置内的总产水量、水土保持和水质净化水文生态系统服务，并采用R软件和WinBUG软件通过蒙特卡罗采样建立各水文生态系统服务间的线性关系；对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果建立各水文生态系统服务间的贝叶斯线性拟合模型和绘制帕累托效益曲线。本发明专利技术通过分别建立的各水文生态系统服务间的线性拟合模型和绘制的帕累托效率曲线，可以为流域土地规划、生物多样性保护和生态补偿提供科学依据。

【技术实现步骤摘要】
一种量化流域水文生态系统服务间权衡与协同关系的计算方法
本专利技术涉及流域水文生态领域，尤其涉及一种量化流域水文生态系统服务间权衡与协同关系的计算方法。
技术介绍
生态系统服务是指生态系统所形成和维持的人类赖以生存和发展的环境条件与效用，为人类直接或间接从生态系统得到的所有效益。其中，流域生态系统提供的总产水量、水土保持和水质净化(营养物质保持量)水文生态系统服务对人类福祉产生重要影响，且这些水文生态系统服务与流域水文循环过程与营养物质循环过程密切相关。生态系统对人类社会福祉提供的服务是多重的，并且各种服务之间相互作用相互联系。所谓生态系统服务权衡是指某些类型生态系统服务的供给由于其他类型的生态系统服务使用的增加而减少的状况。协同是指两种或多种生态系统服务同时增强的情形。生态系统服务权衡与协同是进行区域生态系统服务管理与决策的基础。因此，近来如何构建量化生态系统服务间权衡与协同关系的计算方法成为流域生态学和水环境管理研究的热点。传统的量化流域水文生态系统服务间权衡与协同关系的计算方法多以土地利用与覆被变化为媒介定性或定量地分析生态系统服务，具体的方法有地图对比法、情景分析法和生态-经济综合模型法。以上方法未能从生态水文角度和地理化角度考虑生态系统服务之间的相互作用关系。另外，已有的方法未能以“格局-过程-功能-服务”为主线考虑生态水文计算单元和空间地理单元内各水文生态系统服务间的相互作用关系，且未能将水质净化服务进行细分，忽略了有机营养物质保持量和无机营养物质保持量的空间分布的异质性和权衡性。目前，我国量化流域水文生态系统服务间的权衡与协同关系的计算主要以流域和行政管理单元为计算单元计算各生态系统服务间的相关系数，该方法虽在理论上比较完备，但由于其计算单元的选择未考虑生态系统服务模拟过程的生态水文过程和生态系统服务的地理化，会在很大程度上制约研究结果的实际应用效果；因此，采用何种模型可以模拟基于“格局-过程-功能-服务”的流域水文生态系统服务；如何从生态水文过程角度量化各水文生态系统服务间的相互作用关系；如何从生态系统服务的空间地理化量化各水文生态系统服务间的相互作用关系，是现阶段需要解决的问题。
技术实现思路
本专利技术的目的在于克服现有技术的缺点，提供了一种量化流域水文生态系统服务间权衡与协同关系的计算方法，解决了现有计算管理方法存在的不足。本专利技术的目的通过以下技术方案来实现：一种量化流域水文生态系统服务间权衡与协同关系的计算方法，所述计算方法包括：采用流域水文模型ArcSWAT的流域划分工具将整个流域划分为多个子流域，并将每一个子流域划分为水文响应单元，得到全流域水文响应单元的空间分布矢量图；在每个水文响应单元内模拟水文循环过程和营养物质循环过程，进而模拟得到河道的流量、泥沙含量和营养物质含量，并将其与已知监测断面的流量、泥沙含量和营养物质含量观测值进行对比分析、率定和校正水文模型；在不同条件下模拟各水文响应单元内的水文循环过程和营养物质循环过程，得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量；将它两次得到的结果相减得到泥沙含量、有机营养物质和无机营养物质含量之差，并将总产水量、水土保持和水质净化作为流域水文生态系统服务，进而得到总产水量、水土保持和水质净化水文生态系统服务空间分布图层；统计各水文响应单元相应位置内的总产水量、水土保持和水质净化水文生态系统服务，并采用R软件和WinBUG软件通过蒙特卡罗采样建立各水文生态系统服务间的贝叶斯线性拟合关系；对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果绘制帕累托效益曲线。进一步地，所述在不同条件下计算模拟水文响应单元饿的模拟水文循环过程和营养物质循环过程，得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量包括：根据率定和校正的水文模型模拟各个水文响应单元的水文循环过程和营养物质循环过程，得到各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量；将土壤截留和植被吸收参数配置为0、植被管理和植被覆盖参数配置为1，再次运转模型，模拟水文循环过程和营养物质循环过程，再次得到各HRU内的总产水量、泥沙产生量、有机和无机营养物质产生量。进一步地，所述得到总产水量水文生态系统服务空间分布图层包括：根据率定和校正的模型得到各个水文响应单元内的总产水量数据，并将该数据转化为固定格式后通过ArcGIS软件将其与水文响应单元空间分布图层进行连接，得到总产水量水文生态系统服务的空间分布图层。进一步地，所述得到水土保持和水质净化水文生态系统服务空间分布图层包括：根据不同条件下得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量的差值，并将该差值数据转化为固定格式后通过ArcGIS软件将其与水文响应单元空间分布矢量图层进行连接，得到水土保持和水质净化水文生态系统服务的空间分布图层。进一步地，所述对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果建立各水文生态系统服务间的贝叶斯线性拟合模型包括：统计各水文响应单元相应位置内的总产水量、水土保持和水质净化水文生态系统服务，并采用R软件和WinBUG软件通过蒙特卡罗采样建立各水文生态系统服务间的贝叶斯线性拟合关系；进一步地，所述对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果绘制帕累托效益曲线包括：将基于水文响应单元的各水文生态系统服务空间分布图层转化为栅格格式，并将栅格作为基本地理计算单元；将各水文生态系统服务相除得到相应的比值栅格图层，图层中每个单元格为对应地理位置的两种水文生态系统服务的比值；根据两种水文生态系统服务的比值大学对比值图层每个单元格进行升序排列；根据升序排列的顺序依次分别对单元格对应地理位置的两种水文生态系统服务进行累计求和，并按照结果绘制帕累托效率曲线，根据该曲线的凹凸增减变化趋势判断各生态系统服务之间的相互作用关系。如果该曲线向下凹，且一种生态系统服务呈增加或递减趋势，而另一种生态系统服务呈递减或递增趋势，则该关系为权衡关系；如果该曲线向上凸，且一种生态系统服务呈增加或递减趋势，而另一种生态系统服务也呈增加或递减趋势，则该关系为协同关系。本专利技术具有以下优点：。(1)本专利技术在流域水文模型的基本计算单元(水文响应单元)尺度上分析水文生态系统服务间的关系，使其具有水文意义和生态学意义，能够有效地提高水文生态系统服务的关系效能，使得管理措施更易于实施。(2)本专利技术通过分别建立的各水文生态系统服务间的线性拟合模型和绘制的帕累托效率曲线，可以为流域土地规划、生物多样性保护和生态补偿提供科学依据。附图说明图1为本专利技术方法的流程示意图；图2为本专利技术各流域水文生态系统服务的空间分布图层；图3为本专利技术各流域水文生态系统服务间权衡与协同本文档来自技高网...

【技术保护点】
1.一种量化流域水文生态系统服务间权衡与协同关系的计算方法，其特征在于：所述计算方法包括：/n基于流域地形、气象、土地利用和土壤等空间分布图层基础信息，采用流域水文模型ArcSWAT的流域划分工具将整个流域进行划分，得到子流域分布图，再通过设定土地利用、土壤类型和地形坡度的分布面积阈值，进行空间叠加将每个子流域划分为单个或多个水文响应单元，进而得到全流域水文响应单元空间分布图层；/n在不同条件下模拟各水文响应单元内的水文循环过程和营养物质循环过程，得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量；/n将它两次得到的结果相减得到泥沙含量、有机营养物质和无机营养物质含量之差，并将总产水量、水土保持和水质净化作为流域水文生态系统服务，进而得到总产水量、水土保持和水质净化水文生态系统服务空间分布图层；/n统计各水文响应单元相应位置内的总产水量、水土保持和水质净化水文生态系统服务，并采用R软件和WinBUG软件通过蒙特卡罗采样建立各水文生态系统服务间的贝叶斯线性拟合关系；/n对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果绘制帕累托效益曲线。/n...

【技术特征摘要】
1.一种量化流域水文生态系统服务间权衡与协同关系的计算方法，其特征在于：所述计算方法包括：
基于流域地形、气象、土地利用和土壤等空间分布图层基础信息，采用流域水文模型ArcSWAT的流域划分工具将整个流域进行划分，得到子流域分布图，再通过设定土地利用、土壤类型和地形坡度的分布面积阈值，进行空间叠加将每个子流域划分为单个或多个水文响应单元，进而得到全流域水文响应单元空间分布图层；
在不同条件下模拟各水文响应单元内的水文循环过程和营养物质循环过程，得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量；
将它两次得到的结果相减得到泥沙含量、有机营养物质和无机营养物质含量之差，并将总产水量、水土保持和水质净化作为流域水文生态系统服务，进而得到总产水量、水土保持和水质净化水文生态系统服务空间分布图层；
统计各水文响应单元相应位置内的总产水量、水土保持和水质净化水文生态系统服务，并采用R软件和WinBUG软件通过蒙特卡罗采样建立各水文生态系统服务间的贝叶斯线性拟合关系；
对各水文响应单元内的总产水量、水土保持和水质净化水文生态系统服务进行处理，并根据处理结果绘制帕累托效益曲线。


2.根据权利要求1所述的一种量化流域水文生态系统服务间权衡与协同关系的计算方法，其特征在于：所述在不同条件下模拟各水文响应单元内的水文循环过程和营养物质循环过程，得到两次各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量包括：
通过机理性流域水文模型ArcSWAT的流域划分工具将整个研究流域划分为多个子流域，再通过设定土地利用类型、土壤类型和地形坡度的空间分布面积阈值，并通过空间叠加分析将每一个子流域划分为水文响应单元，进而得到全流域水文响应单元的空间分布矢量图；
模拟在每个水文响应单元上的流域水文循环过程和营养物质循环过程，进而模拟得到河道的流量、泥沙含量和营养物质含量，并将其与已知监测断面的流量、泥沙含量和营养物质含量观测值进行对比分析、率定和校正水文模型；
基于率定和校正的水文模型，模拟各个水文响应单元的水文循环过程和营养物质循环过程，得到各水文响应单元内的总产水量、泥沙产生量、有机和无机营养物质产生量；
将土壤截留和植被吸收参数配置为0、植被管理和植被覆盖参数配置为1，再次运转模型，模拟水文循环过程和营养物质循环过程，再次得到各HRU内的总产水量、泥沙产生量、有机和无机营养物质产生量。


3.根据...

【专利技术属性】
技术研发人员：樊敏，谌书，肖宇婷，赵丽，
申请(专利权)人：西南科技大学，
类型：发明
国别省市：四川;51