CN112241850A - 一种网络化调水的生态补偿标准测算方法 - Google Patents
一种网络化调水的生态补偿标准测算方法 Download PDFInfo
- Publication number
- CN112241850A CN112241850A CN202011171605.9A CN202011171605A CN112241850A CN 112241850 A CN112241850 A CN 112241850A CN 202011171605 A CN202011171605 A CN 202011171605A CN 112241850 A CN112241850 A CN 112241850A
- Authority
- CN
- China
- Prior art keywords
- water
- compensation
- area
- project
- transfer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 667
- 238000012546 transfer Methods 0.000 title claims abstract description 171
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000011161 development Methods 0.000 claims abstract description 29
- 230000008901 benefit Effects 0.000 claims abstract description 28
- 238000011160 research Methods 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000011144 upstream manufacturing Methods 0.000 claims description 26
- 238000005192 partition Methods 0.000 claims description 17
- 230000001105 regulatory effect Effects 0.000 claims description 17
- 238000013461 design Methods 0.000 claims description 11
- 238000010248 power generation Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 8
- 230000002265 prevention Effects 0.000 claims description 8
- 239000002689 soil Substances 0.000 claims description 8
- 238000003911 water pollution Methods 0.000 claims description 8
- 239000010865 sewage Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 4
- 239000003344 environmental pollutant Substances 0.000 claims description 4
- 231100000719 pollutant Toxicity 0.000 claims description 4
- 238000013316 zoning Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 21
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000018109 developmental process Effects 0.000 description 25
- 238000011156 evaluation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000000638 solvent extraction Methods 0.000 description 3
- 239000000356 contaminant Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000007261 regionalization Effects 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009182 swimming Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Geometry (AREA)
- Strategic Management (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Development Economics (AREA)
- Pure & Applied Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Sewage (AREA)
Abstract
本发明涉及一种网络化调水的生态补偿标准测算方法,包括:研究分区;水源区保护成本与发展机会损失测算;调水工程造成损失测算;水源区水量贡献系数和调水对下游影响贡献系数;补偿量计算;补偿主体分担补偿量的计算。本发明通过建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质的变化情况等将技术手段引入工程经济的运算中,使经济预测更加符合现实的自然规律。本发明突破网络化调水格局下的生态补偿标准计算技术瓶颈,有利于明晰利益相关地区的补偿关系,科学制定补偿标准。可用于计算复杂调水格局中分解各项调水工程补偿量和水源流域各地区所获补偿量。
Description
技术领域
本发明涉及一种网络化调水的生态补偿标准测算方法,是一种水利计算方法。
背景技术
我国水资源短缺且时空分布不均,与经济社会发展布局不匹配,建设调水工程是破解这一难题的重要途径。随着水资源开发利用程度的提高和供需矛盾的加剧,部分地区调水工程呈现逐渐增多及加快的趋势,流域水量调入调出关系日趋复杂,跨流域调水呈现网络化发展趋势。为解决水资源开发利用与保护的外部效应问题,国家高度重视实施水生态补偿机制。
补偿标准的计算是建立生态补偿机制的关键技术和难点所在。流域内水源地生态补偿标准的常用计算方法包括保护成本法、机会损失法、水资源或水生态服务价值法、水量水质参数法、支付意愿法、生态足迹法、能值分析法、博弈论方法等。其中,服务价值尚未有公认的评价指标和方法,不同评估方法结果差异较大,且难以从中分离出水源保护的边际效益,难以直接用于补偿实践;以水源保护成本和发展机会损失为依据的计算方法得到较为广泛的认可和应用。
现有跨流域调水工程生态补偿标准的计算多通过应用或改进流域内生态补偿标准计算方法来实现。学者们应用保护成本与机会损失、支付意愿、排污权价值损失、生态足迹、博弈论等方法,计算了南水北调中线水源流域部分地市的补偿标准。对水源区补偿量的分摊方法包括考虑用水量、水资源价值、支付意愿和生态系统服务价值等因素的各区域占比来进行计算。
在调水工程日渐增多,呈现网络化调水格局的发展趋势下,流域水量调入调出关系日趋复杂,现有的流域内生态补偿或单项调水工程生态补偿研究方法具有一定的局限性,难以支撑补偿标准的科学测算和机制的有效建立。
①在研究范围上,缺乏以水源流域为整体的系统研究。调水生态补偿与流域内生态补偿的一个重大区别,就是受影响范围不仅限于水源地,还包括对取水点下游影响河段的补偿。当前大多数调水生态补偿标准测算的补偿对象范围是水源区的部分地区,且缺乏对下游损失的补偿量测算。
②在补偿关系上,网络化调水格局下补偿主体与对象关系更为复杂。在水源流域同时具有多项调水工程情况下,会出现不同工程的水源区相互包含、对下游产生叠加影响等情况,为补偿主体与对象对应关系的确定带来困难,不能简单按照水源流域保护成本和损失的总量来计算补偿标准。
③在计算方法上,网络化调水生态补偿标准研究不能直接沿用流域生态补偿标准的研究思路,应深入研究网络化调水格局下水量平衡、水资源价值运移传递的特点和规律,研发适用的标准计算方法。
现有的流域水量调入调出关系日趋复杂,不同调水工程水源区相互重叠、对下游产生叠加影响等情况,仅从经济角度考虑解决网络化调水格局下生态补偿标准的问题无法满足要求。在同时存在多项调水工程情况下,由于水流运动复杂性,不能简单按照水源流域保护成本和损失的总量来计算补偿标准,需要根据调水工程布局进行逐级分析,考虑各项调水工程与上游保护和下游影响的不同作用关系,同时考虑受水区分享的效益,研发能够从多线路调水格局中分离出单项调水补偿标准的测算方法,并量化水源流域不同地区获得的补偿量。
发明内容
为了克服现有技术的问题,本发明提出了一种网络化调水的生态补偿标准测算方法。所述的方法基于水量平衡和水资源价值运移传递规律,系统考虑水源流域上游保护成本与下游损失,研发适用于网络化调水生态补偿的标准测算方法,为调水生态补偿机制的研究与实践提供科学基础和技术支撑。
本发明的目的是这样实现的:一种网络化调水的生态补偿标准测算方法,所述方法的步骤如下:
步骤1,研究分区:按照调水工程的取水节点,将水源区进行一级分区,在一级分区基础上,按照流域内市域区界进一步做二级分区;设有n个水量调出工程,根据取水位置从上游到下游分别编号为D1,D2,……,Dj,……,Dn,根据n个水量调出工程的取水口节点断面,将水源流域分为m个集水区域,Bj为调水工程Dj与Dj-1取水点区间流域集水范围;
步骤2,水源区保护成本与发展机会损失测算:统计研究区域数年的水源区涵养林建设、水土流失治理、生态移民、水污染防治的主要保护措施投入成本,考虑资金的时间价值,将其核算到现状年平均成本,计算源区单位水量保护成本和发展机会损失,计算如下:
在水源区水质达标情况下,n个水量调出工程中,工程Dj(j=1,2,3,……,n)取水断面上游流域各区域B1,B2,……,Bj水源涵养林建设、水土流失治理、生态移民、水污染防治的保护措施投入,则Bj区域的实际保护成本为CPj,则单位水量保护成本Cj为:
Bj区域发展损失为CRj,单位水量相应的发展机会损失Cxj为:
其中:Wn为Bn区域产生的地表径流量;WIn为该区域调入水量;
步骤3,调水工程造成损失测算:采用情景对比法,分析测算调水工程对下游造成的损失;
情景工况1:仅考虑水源流域已有水利工程,不考虑调水工程;
情景工况2:既考虑已有水利工程,又考虑调水工程;
建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质和环境容量的变化情况,水动力模型:
式中:x为流程;Q为流量;ql为旁侧入流量;ul为旁侧入流主流流速分量;U为主流流速;R为水力半径;Z为水位;A为过水断面面积;t为时间;n为糙率系数;α为动力修正系数;g为重力加速度;
水环境容量模型:
M=31.5[Cd-C0exp(-kL/u)]exp(kX/u)Qr
式中:M为环境容量;Cd为水功能区下断面水质目标控制因子浓度;C0为水功能区上断面水质目标控制因子浓度;k为污染物综合衰减系数;L为水功能区长度;X为排污口下游断面距水功能区下断面流程;u为设计流量下的平均流速;Qr为设计流量;
采用替代成本法计算受影响河段河道内取水成本增加、环境容量下降、航运费用增加的经济损失;调水造成下游水库发电效益损失采用单方水损失法,即利用逐月调水量、单方水发电量及上网电价计算调水工程造成下游Bi区域的取水、治污、发电、航运的损失合计为Li;
步骤4,水源区水量贡献系数和调水对下游影响贡献系数:将水源区不同分区对调水工程的水量贡献程度作为分摊系数,用于计算调水工程所需分摊的源区保护成本及机会损失;
上游第i个区域Bi对调水工程Dj的水量贡献系数βij公式为:
式中:WQi为Bi区域出流径流量;WQj为Bj区域出流径流量;WDj为调水工程Dj的调出水量;
Dj对下游某区域Bi的影响贡献系数γij为:
流域的n项水量调出工程及其划分的m个集水区域,调水工程源区水量贡献系数和对下游影响贡献系数m×n阶矩阵Am×n:
步骤5,补偿量计算:
利用计算得到的水源区单位水量保护成本Cj和机会损失Cxj,与矩阵Am×n的β值及调出水量相乘,得到各调水工程分摊的水源区保护成本Ej和机会损失补偿量Exj:
Ej=Cj×Am×n|β×WDj
Exj=Cxj×Am×n|β×WDj
将调水工程造成下游地区损失Li与矩阵Am×n中的γ值相乘,得到各调水工程需要分摊的下游损失补偿量BLi:
ELi=Li×Am×n|γ;
得到调水工程Dj分摊到区域Bj的补偿标准为CTij,n项调水工程的补偿量形成如下m×n阶矩阵Fm×n:
纵向列的合计值为调水工程Dj对水源流域的补偿量CTj:
横向行的合计值为水源流域某区域应获得的补偿量CTi:
步骤6,补偿主体分担补偿量的计算:采用效益分摊系数法,测算调水工程Dj受水区分享的水资源效益合计为BNj,得到调水工程Dj的受水区域Sij分担的补偿量CSij:
式BNij中为第j项调水工程在受水区Sij的水资源效益,受水区分为r个,则得第j项调水工程的r个受水区分担的补偿量矩阵Gl为:
本发明的优点和有益效果是:本发明通过建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质的变化情况等,将技术手段引入工程经济的运算中,使经济预测更加符合现实的自然规律。本发明拟突破网络化调水格局下的生态补偿标准计算技术瓶颈,有利于明晰利益相关地区的补偿关系,科学制定补偿标准。可用于计算复杂调水格局中分解各项调水工程补偿量和水源流域各地区所获补偿量。为网络化调水格局下生态补偿的有效实施提供技术基础。对于跨流域调水生态补偿标准的计算具有较好的参考价值和应用前景。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明实施例所述测算方法的原理框图;
图2是本发明实施例所述方法的流域调水工程分布示意图;
图3是本发明的应用实例所研究的调水工程分布图;
图4是本发明所述方法的流程图。
具体实施方式
实施例:
本实施例是一种网络化调水的生态补偿标准测算方法,所述方法的通过深入进行多学科理论分析,建立资源与环境经济模型、水动力与水质模型和地理信息系统模型,研究网络化调水生态补偿的基本原理和标准计算方法,原理框图如图1所示。
(1)基本原理:
应用水量平衡、水资源价值运移传递规律和外部效应内部化理论,研究调水生态补偿的基本原理。对于某流域来讲,若存在n个水量调出工程(n为正整数,指水量调出工程的个数),根据取水位置从上游到下游分别编号为D1,D2,……,Dj,……,Dn,如图2所示。根据n个水量调出工程的取水口节点断面,将水源流域分为m个集水区域。Bj为调水工程Dj与Dj-1取水点区间流域集水范围。
水文循环和水量平衡是水文水资源学的基本原理。水量平衡方程式可由水量的收支情况来制定。利用水量平衡规律进行分析,按研究对象系统的空间尺度,大可到全球,小至一个区域或一个河段;也可从大气层到地下水的任何层次,均可写出不同的水量平衡方程式。本发明针对某河段,其输入水量与输出水量之差就是该河段内的水量变化量,如河道槽蓄变量或水库蓄水变化量。在年周期内,调水工程Dn取水点处河段的水量平衡方程为:
Wn+WIn+WQ(n-1)=WQn+WDn+ΔWRn (1)
式中:Wn为Bn区域产生的地表径流量;若Bn区域有水量调入工程,WIn为该区域调入水量;WQ(n-1)为Bn区域上游入流径流量。WQn为Bn区域出流径流量;WDn为调水工程Dn的调出水量;ΔWRn为河道槽蓄量(包括水库蓄水变化量),在多年平均的年尺度上,可认为ΔWRn为零。
水资源价值流具有连续传递特性,即水资源循环过程中入流处获得的价值与出流处传递上来的价值处于平衡状态,出流的总价值等于入流的间接价值。水在向下游流动过程中会产生机会效益价值,调水工程会改变这种机会效益价值的实现区域,带来水资源价值分布格局的变化,在受水区产生价值增值,在水源流域调水断面的下游因水量衰减而造成价值衰减。调水工程Dn取水点处河流断面的水资源价值平衡方程为:
VTn=VQn+VDn (2)
式中:VTn为Dn取水点处河流断面的水资源总价值;VQn为断面下泄水量在下游区域产生的价值;VDn为调出水量在受水区产生的价值。
由于调水工程仅会得到取水点以上源区产水量,并造成下游水量衰减,多项调水工程对应不同的源区和下游影响区,不能简单按照水源流域保护成本和损失的总量来计算补偿标准,需要根据调水工程布局进行逐级分析。图1所示流域的n项水量调出工程D1,D2,......,Dn,对水源流域m个区域的水资源价值影响关系如下表所示。其中Vmn为调水工程Dn对区域Bm水资源价值的影响,对上游区域带来价值增值,给取水断面以下流域带来价值衰减,因而表中虚线右上方的各项V值代表对上游水源区的补偿量,左下方的各项V值代表对下游影响区的补偿量,由此可以形成调水工程对水源流域水资源价值影响的关系矩阵。
表1调水工程对水源流域各区域水资源价值影响关系
应用上述矩阵关系,纵向列可用于分析某项调水工程对水源流域的补偿量,横向行可用于分析测算水源流域某区域应获得的补偿量。如某调水工程Dj的补偿量包括对其取水点上游水源区B1,B2,……,Bj的补偿量V1j,V2j,……,Vjj,还包括对其取水点下游受影响区域Bj+1,……,Bm的补偿量V(j+1)j,V(j+2)j,……,Vmj。
根据上述水量平衡和水资源价值的运移传递规律,在调水工程水源流域水质达标的情况下,受水区为补偿主体区域,水源流域为补偿对象区域,可以摸清各调水工程与取水点上游源区和下游影响区的补偿关系,应用外部成本内部化的经济学原理,测算补偿标准。
(2)补偿标准计算模型:
水源保护成本和损失是确定补偿标准的基础,以此为依据测算的补偿标准具有较好的科学性和可行性。测算某调水工程Dj水源区各项保护成本和机会损失,及其给下游地区带来的损失。按照调水工程的取水节点,将研究对象流域进行一级分区。在一级分区基础上,按照流域内省市域区界可进一步做二级分区。根据该调水工程水源分区的水量贡献程度,及该工程对下游损失的影响程度,测算源区水量贡献系数和调水对下游区域影响贡献系数,进而计算该调水工程所需分摊的补偿量。
①水源区保护成本与发展机会损失:
在水源区水质达标情况下,重点考虑水量调出工程Dj取水断面上游流域各区域B1,B2,……,Bj水源涵养林建设、水土流失治理、生态移民、水污染防治等主要保护措施投入。以Bj区域为例,实际保护成本为CPj,则单位水量保护成本为:
水源区发展机会损失采用区域对比法进行计算,选取与水源区自然条件相近但未因涵养水源而限制经济发展的地区作为参照对象,如水源区周边或下游邻近地区,比较人均财政收入和工资报酬的差异,并通过对关、停、并、转企业等产值、税收等实证调查数据进行校核,评价限制发展的经济损失。如Bj区域发展损失CRj,单位水量相应的发展机会损失为:
②调水工程给下游地区带来的损失:
调水工程可对水源流域取水点下游区域产生水文、水质、水资源开发利用、航运等方面影响。采用情景对比法,分析测算调水工程对下游造成的损失。
情景工况1:仅考虑水源流域已有水利工程,不考虑调水工程。
情景工况2:既考虑已有水利工程,又考虑调水工程。
建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质等变化情况。利用河道内高程信息,概化河道断面,采用一维水流数学模型进行水动力模拟。基于垂向积分的质量和动量守恒方程即圣维南方程组建立的,方程组如下:
式中:x为流程(m);Q为流量(m3/s);ql为旁侧入流量(m3/s);ul为旁侧入流主流流速分量(m/s);U为主流流速(m/s);R为水力半径(m);Z为水位(m);A为过水断面面积(m2);t为时间(s);n为糙率系数;α为动力修正系数;g为重力加速度。
方程组用隐式有限差分法离散,用追赶法求解。计算网格布置为交叉网格方式(交替水位点和流量点)。有断面数据的点设置为水位点,流量点自动布置在相邻水位点的中间和水工建筑物点。
考虑污染物沿河长衰减及稀释倍数变化,采用一维水质模型进行水功能区环境容量计算:
M=31.5[Cd-C0exp(-kL/u)]exp(kX/u)Qr (6)
式中:M—环境容量(t/a);
Cd—水功能区下断面水质目标控制因子浓度(mg/L);
C0—水功能区上断面水质目标控制因子浓度(mg/L);
k—污染物综合衰减系数(1/s);
L—水功能区长度(m);
X—排污口下游断面距水功能区下断面纵向距离(m);
u—设计流量下的平均流速(m/s);
Qr—设计流量(m3/s)。
采用替代成本法计算受影响河段河道内取水成本增加、环境容量下降、航运费用增加等经济损失。替代成本法,亦称替代措施法,是指当水资源供需矛盾日趋尖锐,水量衰减或水质恶化,严重制约水资源开发利用时,可采取诸如节水、增加污水处理和利用成本、增加取水成本等措施来满足用水,可用替代措施的边际成本作为经济损失的间接估算值。调水造成下游水库发电效益损失采用单方水损失法,即利用逐月调水量、单方水发电量及上网电价进行计算。调水工程造成下游Bi区域的取水、治污、发电、航运等方面损失合计为Li,其中,i=1,2,3,……,m,m为正整数。
③源区水量贡献系数和调水对下游影响贡献系数:
源区水量贡献系数:将源区不同分区对调水工程的水量贡献程度作为分摊系数,用于计算调水工程所需分摊的源区保护成本及机会损失。调水工程Dj的调出水量来源于取水点断面上游流域,包括B1,B2,……,Bj范围产生的地表径流量及调入水量。因而,根据水量平衡规律公式(1),可以得到上游各区域产水量对Dj调水量的水量贡献系数,如第i个区域Bi对调水工程Dj的水量贡献系数公式为:
式中βij为水量贡献系数;WQi为Bi区域出流径流量;WQj为Bj区域出流径流量;WDj为调水工程Dj的调出水量。
调水对下游影响贡献系数:该系数用于计算调水工程所需分摊的下游损失。调水工程Dj会影响取水点断面下游的区域Bj+1,Bj+2,……,Bm。下游区域受到上游水量调出与调入工程的叠加影响,因而Dj对下游某区域Bi的影响贡献系数为:
对于图2所示流域的n项水量调出工程及其划分的m个集水区域,由表1和公式(5)、(6)可得调水工程源区水量贡献系数和对下游影响贡献系数m×n阶矩阵如下:
④补偿量计算:
利用公式(3)和公式(4)计算得到的源区单位水量保护成本Cj和机会损失Cxj,与矩阵(9)中的β值及调出水量相乘,得到各调水工程分摊的源区保护成本和机会损失补偿量。
将调水工程造成下游地区损失Li与矩阵(9)中的γ值相乘,得到各调水工程需要分摊的下游损失补偿量。
基于表1调水工程对水源流域各区域水资源价值影响关系,对图1所示流域可得到调水工程Dj分摊到区域Bj的补偿标准为CTij,n项调水工程的补偿量形成如下m×n阶矩阵:
纵向列的合计值为调水工程Dj对水源流域的补偿量CTj:
横向行的合计值为水源流域某区域应获得的补偿量。
⑤补偿主体分担补偿量测算方法:
采用效益分摊系数法,测算调水工程Dj受水区分享的水资源效益合计为BNj。调水工程受水区各地按照所分享国民经济效益和生态效益所占比例,分担水源流域的补偿量,可得到调水工程Dj的受水区域Sij分担的补偿量CSij。
(3)本实施例所述方法的应用实例和具体步骤如下:
应用实例数据:
汉江发源于秦岭南麓,流经陕豫鄂三省,于武汉市汇入长江,全长1577km,流域面积为15.9万km2,多年平均水资源总量573亿m3。已建成南水北调中线一期、引江济汉工程,在建引汉济渭工程、鄂北水资源配置工程,如图3所示。
计算流程如图4所示,流程如下:
步骤1,研究分区:按照调水工程的取水节点,将水源区进行一级分区,在一级分区基础上,按照流域内市域区界进一步做二级分区;设有n个水量调出工程,根据取水位置从上游到下游分别编号为D1,D2,......,Dn,Bn为调水工程Dn与Dn-1取水点区间流域集水范围;
应用实例:考虑汉江流域已建的南水北调一期工程和引汉济渭工程和鄂北水资源配置工程,按照调水工程的取水节点,将汉江流域进行一级分区。在一级分区基础上,按照流域内市域区界进一步做二级分区。
步骤2,水源区保护成本与发展机会损失测算:统计研究区域数年的水源区涵养林建设、水土流失治理、生态移民、水污染防治的主要保护措施投入成本,考虑资金的时间价值,将其核算到现状年平均成本,计算源区单位水量保护成本和发展机会损失,计算如下:
在水源区水质达标情况下,水量调出工程Dj取水断面上游流域各区域B1,B2,……,Bj水源涵养林建设、水土流失治理、生态移民、水污染防治的保护措施投入,则Bj区域的实际保护成本为CPj,则单位水量保护成本Cj为:
水源区发展机会损失采用区域对比法进行计算,选取与水源区自然条件相近但未因涵养水源而限制经济发展的地区作为参照对象,如水源区周边或下游邻近地区,比较人均财政收入和工资报酬的差异,并通过对关、停、并、转企业等产值、税收等实证调查数据进行校核,评价限制发展的经济损失。Bj区域发展损失为CRj,单位水量相应的发展机会损失Cxj为:
其中:Wn为Bn区域产生的地表径流量;WIn为该区域调入水量。
接上述实例:统计2015年以来水源区涵养林建设、水土流失治理、生态移民、水污染防治等主要保护措施投入成本,考虑资金的时间价值,将其核算到现状年平均成本,按公式计算源区单位水量保护成本和发展机会损失。
步骤3,调水工程造成损失测算:调水工程可对水源流域取水点下游区域产生水文、水质、水资源开发利用、航运等方面影响。采用情景对比法,分析测算调水工程对下游造成的损失。
情景工况1:仅考虑水源流域已有水利工程,不考虑调水工程。
情景工况2:既考虑已有水利工程,又考虑调水工程。
建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质的变化情况,水动力模型:
式中:x为流程(m);Q为流量(m3/s);ql为旁侧入流量(m3/s);ul为旁侧入流主流流速分量(m/s);U为主流流速(m/s);R为水力半径(m);Z为水位(m);A为过水断面面积(m2);t为时间(s);n为糙率系数;α为动力修正系数;g为重力加速度。
水环境容量模型:
M=31.5[Cd-C0exp(-kL/u)]exp(kX/u)Qr
式中:M—环境容量(t/a);
Cd—水功能区下断面水质目标控制因子浓度(mg/L);
C0—水功能区上断面水质目标控制因子浓度(mg/L);
k—污染物综合衰减系数(1/s);
L—水功能区长度(m);
X—排污口下游断面距水功能区下断面纵向距离(m);
u—设计流量下的平均流速(m/s);
Qr—设计流量(m3/s)。
采用替代成本法计算受影响河段河道内取水成本增加、环境容量下降、航运费用增加的经济损失。调水造成下游水库发电效益损失采用单方水损失法,即利用逐月调水量、单方水发电量及上网电价计算调水工程造成下游Bi区域的取水、治污、发电、航运的损失合计为Li。
接上述应用实例:
情景工况1:考虑汉江流域已有蓄水工程,包括丹江口水库加高工程,不考虑调水工程。
情景工况2:既考虑已有水利工程,又考虑南水北调一期工程、引江济汉工程、引汉济渭工程和鄂北调水工程。
采用情景对比法进行测算。建立引汉济渭工程取水断面下游至丹江口水库、丹江口水库下游干流(黄家港水文站以下)一维水动力与水质模型,利用多年平均情况下不同工况来水条件,模拟比较逐月流量、水位、水质等参数在有无调水工况下的变化情况。水动力模型:
式中:x为流程;Q为流量;ql为旁侧入流量;ul为旁侧入流主流流速分量;U为主流流速;R为水力半径;Z为水位;A为过水断面面积;t为时间;n为糙率系数;α为动力修正系数;g为重力加速度。
水环境容量模型:
M=31.5[Cd-C0exp(-kL/u)]exp(kX/u)Qr
式中:M为环境容量;Cd为水功能区下断面水质目标控制因子浓度;C0为水功能区上断面水质目标控制因子浓度;k为污染物综合衰减系数;L为水功能区长度;X为排污口下游断面距水功能区下断面纵向距离;u为设计流量下的平均流速;Qr为设计流量。
由于调水对下游地区社会影响和生态影响难以量化,主要采用单方水发电损失逐月调算法计算水电站(石泉、喜河、蜀河等)发电损失,采用替代成本法计算干流取水、治污、航运等损失。
步骤4,水源区水量贡献系数和调水对下游影响贡献系数:将水源区不同分区对调水工程的水量贡献程度作为分摊系数,用于计算调水工程所需分摊的源区保护成本及机会损失;
上游第i个区域Bi对调水工程Dj的水量贡献系数βij公式为:
式中:WQi为Bi区域出流径流量;WQj为Bj区域出流径流量;WDj为调水工程Dj的调出水量。
Dj对下游某区域Bi的影响贡献系数γij为:
流域的n项水量调出工程及其划分的m个集水区域,调水工程源区水量贡献系数和对下游影响贡献系数m×n阶矩阵Am×n:
针对引汉济渭工程、南水北调一期工程、鄂北调水工程及其对应的水源区和下游影响区,应用公式(7)~(9),考虑多年平均水文条件下的水量平衡和设计调水规模,可得水源流域调水生态补偿标准分摊系数m×n阶矩阵Am×n,包括源区水量贡献系数β和调水对下游影响贡献系数γ。
步骤5,补偿量计算:
利用计算得到的水源区单位水量保护成本Cj和机会损失Cxj,与矩阵Am×n的β值及调出水量相乘,得到各调水工程分摊的水源区保护成本Ej和机会损失补偿量Exj:
Ej=Cj×Am×n|β×WDj
Exj=Cxj×Am×n|β×WDj
将调水工程造成下游地区损失Li与矩阵Am×n中的γ值相乘,得到各调水工程需要分摊的下游损失补偿量BLi:
ELi=Li×Am×n|γ。
得到调水工程Dj分摊到区域Bj的补偿标准为CTij,n项调水工程的补偿量形成如下m×n阶矩阵Fm×n:
纵向列的合计值为调水工程Dj对水源流域的补偿量CTj:
横向行的合计值为水源流域某区域应获得的补偿量CTi:
应用实例计算:利用计算得到的汉江水源区单位水量保护成本和机会损失,与矩阵(9)中β值及调出水量相乘,得到各调水工程分摊的源区补偿量。将调水工程造成下游地区损失与矩阵(9)中γ值相乘,得到各调水工程分摊的下游损失补偿量。
得到形如矩阵(10)的各项调水工程补偿量矩阵。按照公式(11)计算矩阵的列合计值,作为某单项调水工程对水源流域的补偿量,各行的合计值为各区域应获得补偿量。
步骤6,补偿主体分担补偿量的计算:采用效益分摊系数法,测算调水工程Dj受水区分享的水资源效益合计为BNj,得到调水工程Dj的受水区域Sij分担的补偿量CSij:
式BNij中为第j项调水工程在受水区Sij的水资源效益,受水区分为r个,则得第j项调水工程的r个受水区分担的补偿量矩阵Gr为:
与现状年的保护成本和损失测算相协调,利用现状年经济社会发展数据和多年平均条件下的调水规模,采用效益分摊系数法,测算引汉济渭工程、南水北调一期工程、鄂北调水工程受水区效益,包括国民经济用水和生态用水效益,得到受水区各地分享效益所占比例。按照公式(12),测算受水区域分担水源流域的补偿量。
最后应说明的是,以上仅用以说明本发明的技术方案而非限制,尽管参照较佳布置方案对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案(比如各种模型的形式、各种公式的运用、步骤的先后顺序等)进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (1)
1.一种网络化调水的生态补偿标准测算方法,其特征在于,所述方法的步骤如下:
步骤1,研究分区:按照调水工程的取水节点,将水源区进行一级分区,在一级分区基础上,按照流域内市域区界进一步做二级分区;设有n个水量调出工程,根据取水位置从上游到下游分别编号为D1,D2,……,Dj,……,Dn,根据n个水量调出工程的取水口节点断面,将水源流域分为m个集水区域,Bj为调水工程Dj与Dj-1取水点区间流域集水范围;
步骤2,水源区保护成本与发展机会损失测算:统计研究区域数年的水源区涵养林建设、水土流失治理、生态移民、水污染防治的主要保护措施投入成本,考虑资金的时间价值,将其核算到现状年平均成本,计算源区单位水量保护成本和发展机会损失,计算如下:
在水源区水质达标情况下,水量调出工程Dj取水断面上游流域各区域B1,B2,……,Bj水源涵养林建设、水土流失治理、生态移民、水污染防治的保护措施投入,则Bj区域的实际保护成本为CPj,则单位水量保护成本Cj为:
Bj区域发展损失为CRj,单位水量相应的发展机会损失Cxj为:
其中:Wn为Bn区域产生的地表径流量;WIn为该区域调入水量;
步骤3,调水工程造成损失测算:采用情景对比法,分析测算调水工程对下游造成的损失;
情景工况1:仅考虑水源流域已有水利工程,不考虑调水工程;
情景工况2:既考虑已有水利工程,又考虑调水工程;
建立调水工程取水断面下游河段水动力与水质模型,模拟不同工况条件下水利工程下游河段流量、水位、水质的变化情况,水动力模型:
式中:x为流程;Q为流量;ql为旁侧入流量;ul为旁侧入流主流流速分量;U为主流流速;R为水力半径;Z为水位;A为过水断面面积;t为时间;n为糙率系数;α为动力修正系数;g为重力加速度;
水环境容量模型:
M=31.5[Cd-C0exp(-kL/u)]exp(kX/u)Qr
式中:M为环境容量;Cd为水功能区下断面水质目标控制因子浓度;C0为水功能区上断面水质目标控制因子浓度;k为污染物综合衰减系数;L为水功能区长度;X为排污口下游断面距水功能区下断面纵向距离;u为设计流量下的平均流速;Qr为设计流量;
采用替代成本法计算受影响河段河道内取水成本增加、环境容量下降、航运费用增加的经济损失;调水造成下游水库发电效益损失采用单方水损失法,即利用逐月调水量、单方水发电量及上网电价计算调水工程造成下游Bi区域的取水、治污、发电、航运的损失合计为Li;
步骤4,水源区水量贡献系数和调水对下游影响贡献系数:将水源区不同分区对调水工程的水量贡献程度作为分摊系数,用于计算调水工程所需分摊的源区保护成本及机会损失;
上游第i个区域Bi对调水工程Dj的水量贡献系数βij公式为:
式中:WQi为Bi区域出流径流量;WQj为Bj区域出流径流量;WDj为调水工程Dj的调出水量;
Dj对下游某区域Bi的影响贡献系数γij为:
流域的n项水量调出工程及其划分的m个集水区域,调水工程源区水量贡献系数和对下游影响贡献系数m×n阶矩阵Am×n:
步骤5,补偿量计算:
利用计算得到的水源区单位水量保护成本Cj和机会损失Cxj,与矩阵Am×n的β值及调出水量相乘,得到各调水工程分摊的水源区保护成本Ej和机会损失补偿量Exj:
Ej=Cj×Am×n|β×WDj
Exj=Cxj×Am×n|β×WDj
将调水工程造成下游地区损失Li与矩阵Am×n中的γ值相乘,得到各调水工程需要分摊的下游损失补偿量BLi:
ELi=Li×Am×n|γ;
得到调水工程Dj分摊到区域Bj的补偿标准为CTij,n项调水工程的补偿量形成如下m×n阶矩阵Fm×n:
纵向列的合计值为调水工程Dj对水源流域的补偿量CTj:
横向行的合计值为水源流域某区域应获得的补偿量CTi:
步骤6,补偿主体分担补偿量的计算:采用效益分摊系数法,测算调水工程Dj受水区分享的水资源效益合计为BNj,得到调水工程Dj的受水区域Sij分担的补偿量CSij:
式BNij中为第j项调水工程在受水区Sij的水资源效益,受水区分为r个,则得第j项调水工程的r个受水区分担的补偿量矩阵Gr为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171605.9A CN112241850A (zh) | 2020-10-28 | 2020-10-28 | 一种网络化调水的生态补偿标准测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171605.9A CN112241850A (zh) | 2020-10-28 | 2020-10-28 | 一种网络化调水的生态补偿标准测算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112241850A true CN112241850A (zh) | 2021-01-19 |
Family
ID=74170270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011171605.9A Pending CN112241850A (zh) | 2020-10-28 | 2020-10-28 | 一种网络化调水的生态补偿标准测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112241850A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111667145A (zh) * | 2020-05-06 | 2020-09-15 | 武汉大学 | 一种基于非合作博弈的河道内外用水冲突协商方法 |
CN113420512A (zh) * | 2021-06-11 | 2021-09-21 | 中国水利水电科学研究院 | 一种计算跨界河流水功能区水质达标率的方法 |
CN114117848A (zh) * | 2021-11-10 | 2022-03-01 | 北京师范大学 | 一种基于多模型耦合的流域水环境模拟预测的方法及装置 |
CN114462120A (zh) * | 2022-01-12 | 2022-05-10 | 中国长江三峡集团有限公司 | 一种基于一维水动力模型的河道槽蓄量计算方法 |
CN118094323A (zh) * | 2024-02-27 | 2024-05-28 | 中国水利水电科学研究院 | 一种基于svr耦合机理模型的调水工程水质预测方法 |
-
2020
- 2020-10-28 CN CN202011171605.9A patent/CN112241850A/zh active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111667145A (zh) * | 2020-05-06 | 2020-09-15 | 武汉大学 | 一种基于非合作博弈的河道内外用水冲突协商方法 |
CN111667145B (zh) * | 2020-05-06 | 2023-04-28 | 武汉大学 | 一种基于非合作博弈的河道内外用水冲突协商方法 |
CN113420512A (zh) * | 2021-06-11 | 2021-09-21 | 中国水利水电科学研究院 | 一种计算跨界河流水功能区水质达标率的方法 |
CN113420512B (zh) * | 2021-06-11 | 2022-04-01 | 中国水利水电科学研究院 | 一种计算跨界河流水功能区水质达标率的方法 |
CN114117848A (zh) * | 2021-11-10 | 2022-03-01 | 北京师范大学 | 一种基于多模型耦合的流域水环境模拟预测的方法及装置 |
CN114117848B (zh) * | 2021-11-10 | 2022-11-11 | 北京师范大学 | 一种基于多模型耦合的流域水环境模拟预测的方法及装置 |
CN114462120A (zh) * | 2022-01-12 | 2022-05-10 | 中国长江三峡集团有限公司 | 一种基于一维水动力模型的河道槽蓄量计算方法 |
CN114462120B (zh) * | 2022-01-12 | 2024-05-28 | 中国长江三峡集团有限公司 | 一种基于一维水动力模型的河道槽蓄量计算方法 |
CN118094323A (zh) * | 2024-02-27 | 2024-05-28 | 中国水利水电科学研究院 | 一种基于svr耦合机理模型的调水工程水质预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112241850A (zh) | 一种网络化调水的生态补偿标准测算方法 | |
Zhou et al. | Industrial structural upgrading and spatial optimization based on water environment carrying capacity | |
Wang et al. | Application of genetic algorithm to land use optimization for non-point source pollution control based on CLUE-S and SWAT | |
Wang et al. | Optimization of storage tank locations in an urban stormwater drainage system using a two-stage approach | |
Dou et al. | Measurement and assessment of water resources carrying capacity in Henan Province, China | |
Xie | An environmentally extended social accounting matrix | |
Saadatpour et al. | Meta-model assisted 2D hydrodynamic and thermal simulation model (CE-QUAL-W2) in deriving optimal reservoir operational strategy in selective withdrawal scheme | |
CN105913146B (zh) | 南方湿润地区的水资源优化配置系统 | |
CN111259530B (zh) | 多入河口的中小流域水质达标动态补水量预测方法 | |
Liu et al. | Optimal allocation of water quantity and waste load in the Northwest Pearl River Delta, China | |
CN105297827A (zh) | 一种考虑多用户需水和多水源供水的水资源配置方法 | |
CN109598376B (zh) | 一种优化生态补水量的方法 | |
CN105389743A (zh) | 一种流域水量-水质联合调控优化方法 | |
Kang et al. | Optimal allocation of water resources considering water quality and the absorbing pollution capacity of water | |
Deng et al. | Modelling the coupling evolution of the water environment and social economic system using PSO-SVM in the Yangtze River Economic Belt, China | |
Zhang et al. | Simulation study on the impact of ecological water replenishment on reservoir water environment based on MIKE21——taking Baiguishan Reservoir as an example | |
Nozari et al. | Simulation and optimization of control system operation and surface water allocation based on system dynamics modeling | |
Tang et al. | A numerical study of hydrodynamic processes and flood mitigation in a large river-lake system | |
Yang et al. | Developing a comprehensive evaluation method for Interconnected River System Network assessment: A case study in Tangxun Lake group | |
Cheng et al. | Application of HEC-HMS parameter regionalization in small watershed of hilly area | |
Li et al. | Comprehensive evaluation and scenario simulation of water resources carrying capacity: a case study in Xiong’an New Area, China | |
Ji et al. | Predictive analysis of water resource carrying capacity based on system dynamics and improved fuzzy comprehensive evaluation method in Henan Province | |
Wang et al. | A hydraulics-based analytical method for artificial water replenishment in wetlands by reservoir operation | |
Wang et al. | Calculation of city total maximum allocated load of total nitrogen for jurisdictions in Qingdao, China: a water quality-based modeling approach | |
Jieying et al. | Spatial and temporal variation and driving factors of ecological carrying capacity in the Pan-Pearl River Basin, China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210119 |
|
RJ01 | Rejection of invention patent application after publication |