CN103810401B - 一种人类活动影响分离的二维径流还原方法 - Google Patents
一种人类活动影响分离的二维径流还原方法 Download PDFInfo
- Publication number
- CN103810401B CN103810401B CN201410092691.2A CN201410092691A CN103810401B CN 103810401 B CN103810401 B CN 103810401B CN 201410092691 A CN201410092691 A CN 201410092691A CN 103810401 B CN103810401 B CN 103810401B
- Authority
- CN
- China
- Prior art keywords
- activity
- runoff
- mankind
- state
- watershed
- 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.)
- Expired - Fee Related
Links
- 230000000694 effects Effects 0.000 title claims abstract description 199
- 238000000034 method Methods 0.000 title claims abstract description 87
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 63
- 238000004162 soil erosion Methods 0.000 claims description 29
- 238000001556 precipitation Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 206010016807 Fluid retention Diseases 0.000 claims description 4
- 238000001704 evaporation Methods 0.000 claims description 4
- 230000008020 evaporation Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 235000012149 noodles Nutrition 0.000 claims description 4
- 239000003673 groundwater Substances 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 description 22
- 230000008569 process Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 239000002352 surface water Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002262 irrigation Effects 0.000 description 2
- 238000003973 irrigation Methods 0.000 description 2
- 238000002715 modification method Methods 0.000 description 2
- 230000036299 sexual function Effects 0.000 description 2
- 241000196171 Hydrodictyon reticulatum Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 244000144980 herd Species 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000003621 irrigation water Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000005413 snowmelt Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Volume Flow (AREA)
Abstract
一种人类活动影响分离的二维径流还原方法,包括:步骤一:构建界壳‑泛系二元自动机模型;步骤二:确定人类活动‑流域下垫面系统的输入、输出和内部状态要素;步骤三:计算人类活动对流域下垫面条件的影响量;步骤四:获得人类活动对径流的间接影响量;步骤五:确定满足下垫面条件一致性要求的人类活动‑径流系统的输入、输出和内部状态要素;步骤六:基于已构建的人类活动对径流直接影响分离的系统分析法,求出人类活动对径流的直接影响量;以及步骤七:根据计算出来的人类活动对径流的间接影响量和人类活动对径流的直接影响量,计算出在满足流域下垫面条件一致性的天然径流量。该方法简单易操作,计算速度快、精度高,计算成果更全面、更科学。
Description
技术领域
本发明涉及径流还原方法领域,尤其涉及一种人类活动影响分离的二维径流还原方法。
背景技术
实测径流的还原计算,是目前国内外进行水资源评价和工程设计时首先要做的工作,还原的目的,是为了使河川径流计算成果能够基本上反映天然情况,使资料系列具有一致性,以便采用数理统计方法对其进行频率分析。水量还原工作,在很长时间里,均是通过对实测水文要素中隐含的人类活动影响扣除掉,“还原”到流域水资源的天然“本底”状态,即“实测-还原”的一元静态模式。但传统的还原方法仅考虑人类活动对径流的直接影响量,没有考虑由于人类活动对流域下垫面条件的改变而对径流的间接影响,导致还原成果精度不高、可靠性较低。为应对日益强烈的人类活动对流域的直接和间接影响,保护水资源安全和促进流域水资源的可持续开发利用,开发一种耦合人类活动对径流的直接和间接影响的二维径流还原方法是十分必要的。
综合目前的径流还原方法,从评价的理论和考虑的因素来看主要有分项还原法、气象因子关系法(径流双曲线法、流域蒸发差值法、降雨径流关系法)、综合修正法、向后还原法、水文模拟法、投影寻踪回归分析法、二元动态水资源演化法、BP神经网络法。
(1)分项还原法:分项还原法是生产中最常用的方法。该方法通过对农业灌溉耗水量、工业用水耗损量、生活用水耗损量、水库渗漏量、分洪量以及跨流域调水量等各项水量进行计算,依据水量平衡原理,把控制断面实测径流量与控制断面以上各项还原水量之和认为是控制断面的天然年径流量。分项还原法其需要对社会经济、农业灌溉、工业用水、城镇生活用水、水利工程等方面进行调查,调查的年份较长,需要的资料太多,工作量巨大,调查的数据和涉及的参数繁多,且这些参数大多是各地实验和经验值,在不同的时间或空间尺度上是变化的,在还原时会引起累积误差,还原精度不高。同时该方法只能用于人类活动影响较小的流域,在人类活动影响比较强烈的流域,该法计算得到的天然水量,不能代表水资源的演化、发展趋势,不能用于实践。
(2)气象因子关系法:该法主要通过蒸发量、降水量等因子建立与径流量关系,基于该关系进行天然径流计算。该法计算精度较差,且要进行外延时,外延部分的精度不易控制等等。
(3)综合修正法:该法对操作人员要求很高,既要求熟悉研究区域的水文情势,又要具有丰富的水文工作经验,故该法在实际中的应用并不多。
(4)向后还原法:该法将区域内各代表站的天然年径流量系列折算(换算)为现状条件下的年径流系列。但由于实际当中“现状条件”的不断变化,“向后还原”方法在实际应用中较难操作。
(5)水文模拟法:该法需要确定一定的模型结构和参数,这些结构与参数代表了流域的水文规律与水文特性,通过对这些结构与参数进行率定,实现对水文过程的模拟。但该方法需要确定模型的结构,率定与检验各种参数,模型结构与参数的确定受制于一定的经验,同时,模型存在着一定的结构误差与参数误差,使得该方法的应用受到一定的限制。
(6)投影寻踪回归分析法:该方法是把影响产流量的因素如降水(也包括前期降水)、蒸发、气温、融雪等与径流量之间的非线性关系用投影寻踪法进行降维(降到1-3维空间上),找出径流量与这些影响因素之间的内在联系,确定径流序列的突变点,获得还原水量。但该法受制于投影方向、岭函数和突变点,随意性较大,计算成果精度较低。
(7)二元动态水资源演化法:该法基于二元动态水资源演化模式的“实测-分离-耦合-建模-调控”思想和天然状态下的降雨P—径流R函数关系式R=f(P),利用经验统计方法,找出影响径流形成的主要因子,利用实测径流量与人工耗水量之和建立与这些影响因子之间的经验关系,从而获得受人类活动影响后的流域下垫面条件的降水径流函数关系f'。但该方法考虑的人工耗水量与水库蓄水量的变化量相当于分项还原法中的分项还原量,因此,该法本质上讲属于分项还原法范畴,仍无法克服分项还原法的不足。
(8)BP神经网络法:应用BP网络,以影响因子如农业灌溉地表水净耗水量、工业引用地表水净耗水量、林牧渔业引用地表水净耗水量、水库蓄水变量、水库附加蒸发损失量、跨流域调水量、城镇生活净耗水量、农村人畜吃水净耗水量和水库渗漏损失量、降雨量、下渗系数、年平均温度等作为输入层,中间层为隐含层,输出层为天然径流量或还原水量,通过对历史数据的训练和拟合,当网络满足精度要求后,可用于预测天然径流量或还原水量。但该法需要考虑大量的影响因子和历史数据序列,这些值的确定,需要耗费大量的工作量,具有较大的难度和不确定性,当引入BP神经网络后,可能会加重还原计算的工作量和不确定性。
上述还原方法各有特色,总的来说,其主要还是离不开分项还原法所要求的对人类消耗水和流域蓄水变化、引水、分洪等的还原和对影响径流的气象因素的分析。但在强烈人类活动影响的流域,由于人类活动的影响改变了流域下垫面条件,在这种下垫面条件改变前后测得的实测径流值,其形成机理是不一致的,但上述方法均未考虑人类活动对径流的这种间接影响。因此,在强烈人类活动影响的流域,为有效、准确计算天然径流量和评价区域水资源可利用量,有必要寻求一种有效的强烈人类活动影响下的径流还原方法。
发明内容
针对上述问题,本发明的目的是提供一种人类活动影响分离的二维径流还原方法,它以界壳理论-泛系方法为基础,巧妙结合系统分析法与水保措施计算法,从人类活动对径流的间接影响和直接影响角度,提出一种新的人类活动影响分离的二维径流还原方法。本发明首先建立界壳-泛系二元自动机模型,界定人类活动-流域下垫面系统的输入、输出和内部状态要素,然后提出人类活动对流域下垫面条件影响分离的系统分析法,计算人类活动对流域下垫面条件的影响量;进而提出水保措施计算法来量化该下垫面条件的变化量对径流的影响,确定人类活动对径流的间接影响量;同时界定下垫面条件满足一致性要求的人类活动-径流系统的输入、输出和内部状态要素,基于构建的人类活动对径流直接影响分离的系统分析法,求出人类活动对径流的直接影响量;最后结合人类活动对径流的间接影响量和直接影响量,得到人类活动对径流的总影响量,获得还原后的天然径流量。
为解决上述问题,本发明采取以下技术方案:
一种人类活动影响分离的二维径流还原方法,它适用于分析和计算人类活动影响下的人类活动对径流的间接影响和直接影响,该方法具体步骤如下:
步骤一:构建界壳-泛系二元自动机模型。基于界壳理论和泛系方法,构建水资源系统的界壳-泛系二元自动机模型,如下所示:
式中,S为水资源系统;A为水资源系统周界;B为水资源系统的系里;X为水资源系统的输入;K为水资源系统的内部状态;Y为水资源系统的输出;f为水资源系统内部状态与输入之间的表里、因果关系,体现的是系里的可控性功能;g为水资源系统内部状态的外部表现,体现的是系里的一种观测效果;h为水资源系统输入输出的转换机制,即系统周界作用。
步骤二:确定人类活动-流域下垫面系统的输入、输出和内部状态要素。当不存在人类活动对流域的下垫面条件影响时,流域下垫面系统在系统输入及其自身的作用下发生着演化。当受人类活动影响时,形成人类活动-流域下垫面系统。此时,系统内部状态取为流域的梯田、林地、草地和坝地四种水保措施面积,输入取为流域年平均降水量、气温和日照小时数,系统输出取为下一年度的流域的梯田、林地、草地和坝地四种水保措施面积。
步骤三:基于建立的人类活动对流域下垫面条件影响分离的系统分析法,计算人类活动对流域下垫面条件的影响量。设没有人类活动影响的流域下垫面系统内部状态的时间序列为{K1,K2,…,Ki,Kj,…},由系统周界A的关系式有:
Kj=f(X,Ki) (2)
式中,Ki是i时段初的状态量;Kj是i时段末(即j时段初)的状态量;X是系统输入。输入、状态和输出均是时间t的函数。
设人类活动-流域下垫面系统中的人类活动在i时段对系统内部状态的影响作用为Mi,由于在当前条件下人类很难对系统输入施加影响,故系统内部状态在天然的系统输入和人类活动的影响下发生着演化,故有:
Kj=f(X,MiKi) (3)
令U=MiKi,对式(3)两边求时间t的导数,有:
令分别称为输入X对系统输出Kj的弹性系数、包含有人类活动影响作用Mi的状态Ki对系统输出Kj的弹性系数。当计算的时间步长相同时,上式经整理可化简为:
从而由上式可求出i时段的人类活动对状态的影响量令其为hi,即把hi从1时段开始累积,即可得到以时间序列{K1,K2,…,Ki,Kj,…}的初始状态K1为参考点的i时段人类活动对状态的影响量Hi(即人类活动对流域下垫面条件的影响量):
从而可确定分离人类活动影响后的系统内部状态量Kni为:
Kni=Ki-Hi (7)
式中,Kni表示以研究序列的初始状态K1为参照点且没有人类活动影响的天然系统自身演化时状态量,即为分离人类活动影响后的流域天然下垫面条件变量值。
步骤四:利用提出的水保措施计算法(如式(8)所示)计算由于人类活动对流域下垫面的影响而导致的径流变化量,从而获得人类活动对径流的间接影响量。采用梯田、林地、草地和坝地四种水保措施面积变化量来表征人类活动对流域下垫面条件的影响(由式(6)可确定),从而由下式可求出这几种水保措施对径流的影响量,确定i时段人类活动对径流的间接影响量Wki。
式中,j=1,2,3,4,分别表示梯田、林地、草地和坝地四种水保措施;θ为地下水补给系数;di为i时段的天然径流模数,m3·km-2;Hji为i时段第j种水保措施面积,km2;uj为水保措施的利用率,%;ηj为水保措施的保水率,%。
步骤五:在分离出人类活动对流域下垫面条件的影响量之后,确定满足下垫面条件一致性要求的人类活动-径流系统的输入、输出和内部状态要素。由式(7)求出获得分离人类活动影响后的天然下垫面条件Kni后,系统内部状态量就在自然状态下演化,如果此时没有人类活动对径流进行直接影响时,系统的输入在这种天然状态值情况下会有一个输出,设为Wn。当有人类活动对径流施加直接影响时,形成人类活动-径流影响系统。此时,系统内部状态取为剔除人类活动影响后的梯田面积、林地面积、草地面积和坝地面积;系统输入取为流域的年平均降水量、年平均蒸发量和年平均气温;系统输出取为流域出口水文站的实测年径流量。
步骤六:基于已构建的人类活动对径流直接影响分离的系统分析法,求出人类活动对径流的直接影响量。设人类活动-径流系统中的人类活动对系统输出(径流)的直接影响系数为Md,系统的输出取为实测径流(用Wm表示),系统的输入取为降水量、气温和蒸发气象要素(用Xr表示),系统的内部状态取为分离人类活动影响后的梯田面积、林地面积、草地面积和坝地面积(用Kn表示),则有Wm=MdWn(Xr,Kn),经分析可有:
式中,为系统输入对系统输出的弹性系数;为系统状态对系统输出的弹性系数。
因此,由式(9)可求出i时段人类活动对径流的直接影响系数Mdi,故i时段人类活动对径流的直接影响量Wdi为:
Wdi=(1-Mdi)Wmi/Mdi (10)
步骤七:根据计算出来的人类活动对径流的间接影响量Wk和人类活动对径流的直接影响量Wd,从间接影响和直接影响二维角度获得在满足流域下垫面条件一致性的天然径流量W:
W=Wm+Wd+Wk (11)
本发明由于采取以上技术方案,其具有以下优点:
1.该方法需要的资料较少,且所需资料大都可以免费获取,大大降低了计算成本。
2.该方法简单易操作,计算速度快,计算精度高,计算成果更全面、更科学,既考虑了人类活动对径流的直接影响,也考虑了人类活动对径流的间接影响。
3.方法具有更好的适用性,既适用于人类活动影响强度较低的流域,也适用于人类活动影响强烈的流域。
附图说明
图1为本发明方法的流程框图。
图2为人类活动对流域下垫面四种水保措施(梯田、林地、草地和坝地)的面积影响量图。
图3为人类活动对径流的间接、直接和总影响量图。
图4为还原所得的天然径流结果与实测径流对比图。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步详细说明。
如图1所示,本发明一种人类活动影响分离的二维径流还原方法,包括界壳-泛系二元自动机模型的构建、人类活动-流域下垫面系统要素的确定、人类活动对流域下垫面条件的影响计算、人类活动对径流的间接影响量计算、人类活动-径流系统要素的确定、人类活动对径流直接影响量的计算以及天然径流量的确定共七部分。
以已建立的人类活动对流域下垫面条件影响分离的系统分析法和人类活动对径流直接影响分离的系统分析法为例,本发明的具体实施按照以下步骤进行:
步骤一:构建界壳-泛系二元自动机模型。基于界壳理论和泛系方法,构建水资源系统的界壳-泛系二元自动机模型;如下所示:
式中,S为水资源系统;A为水资源系统周界;B为水资源系统的系里;X为水资源系统的输入;K为水资源系统的内部状态;Y为水资源系统的输出;f为水资源系统内部状态与输入之间的表里、因果关系,体现的是系里的可控性功能;g为水资源系统内部状态的外部表现,体现的是系里的一种观测效果;h为水资源系统输入输出的转换机制,即系统周界作用。
步骤二:根据构建的界壳-泛系二元自动机模型,确定人类活动-流域下垫面系统的输入、输出和内部状态要素。系统内部状态取为流域的梯田、林地、草地和坝地四种水保措施面积(分别以ft、fl、fc和fb表示);系统输入取为流域年平均降水量P、气温T和日照小时数S;系统输出取为下一年度的流域的梯田、林地、草地和坝地四种水保措施面积。
步骤三:根据已构建的人类活动对流域下垫面条件影响分离的系统分析法,计算人类活动-流域下垫面系统中的人类活动对流域下垫面条件的影响量。
(1)根据步骤一中确定的系统输入、输出和内部状态数据,采用“势平衡法”计算弹性系数α和β:
式中,r1=P1/P01,r2=P2/P02,P1=Kj/Xj,P2=Kj/Kj-1,P01和P02分别为P1和P2的前一年值。
(2)求出α和β后,把其代入到下式中,求出第i年人类活动对下垫面条件的影响强度
从而确定第i年的人类活动对状态的影响量把其从1时段开始累积,即可得到以初始状态K1为参考点的i时段人类活动对状态的影响量Hi(即人类活动对流域下垫面条件的影响量):
进而确定出人类活动对梯田、林地、草地和坝地四种水保措施的面积影响值(fhti、fhli、fhci和fhbi),且由式Fn=f-fh可获得分离人类活动影响后的这四种水保措施面积Fnti、Fnli、Fnci和Fnbi。
步骤四:利用提出的水保措施计算法(下式(5)所示)计算出四种水保措施面积fht、fhl、fhc和fhb对径流的影响量,获得人类活动对径流的间接影响量Wk:
式中,j=1,2,3,4,分别表示梯田、林地、草地和坝地四种水保措施;θ为地下水补给系数;di为i时段的天然径流模数,m3·km-2;fhji为i时段第j种水保措施面积,km2;uj为水保措施的利用率,%;ηj为水保措施的保水率,%。
步骤五:在分离出人类活动对流域下垫面条件的影响量之后,确定下垫面条件满足一致性要求的人类活动-径流系统的输入、输出和内部状态要素。系统内部状态取为剔除人类活动影响后的梯田面积Knti、林地面积Knli、草地面积Knci和坝地面积Knbi;系统输入取为流域的年平均降水量P、年平均蒸发量E和年平均气温T;系统输出取为流域出口断面的实测年径流量Wm。
步骤六:由已构建的人类活动对径流直接影响分离的系统分析法,计算人类活动-径流系统中的人类活动对径流的直接影响量。
(1)由下式计算出人类活动对径流的直接影响系数Md:
(2)求出人类活动对径流的直接影响系数Md后,代入下式计算出人类活动对径流的直接影响量Wd:
Wd=(1-Md)Wm/Md (7)
步骤七:根据计算出来的人类活动对径流的间接影响量Wk和人类活动对径流的直接影响量Wd,由下式计算出在满足流域下垫面条件一致性的天然径流量W:
W=Wm+Wd+Wk (8)
实施案例
按照上述方法分析流域人类活动对径流的间接影响和直接影响,获得还原后的天然径流量。本发明以泾河张家山水文站以上流域47年(1957-2003年)的径流还原为例,通过本发明提出的新方法分析人类活动对流域下垫面条件的影响,获得人类活动对径流的间接影响量,并通过构建下垫面条件满足一致性要求的人类活动-径流系统,确定人类活动对径流的直接影响,从而计算出人类活动对径流的总影响量,获得还原后的天然径流量。
其结果,分别见图2、图3、图4。
图2为人类活动对流域下垫面四种水保措施(梯田、林地、草地和坝地)的面积影响量图。
图3为人类活动对径流的间接、直接和总影响量图。
图4为还原所得的天然径流结果与实测径流对比图。
从上述实例可以看出,本发明提供的一种新的人类活动影响分离的二维径流还原方法,可以分析人类活动对径流的直接和间接影响,对于分析人类活动的水文效应、解决人类活动影响下复杂水文系统的预报和模拟等问题具有重大的现实意义。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (6)
1.一种人类活动影响分离的二维径流还原方法,包括以下步骤:
步骤一:构建界壳-泛系二元自动机模型;
步骤二:确定人类活动-流域下垫面系统的输入、输出和内部状态要素;
步骤三:基于建立的人类活动对流域下垫面条件影响分离的系统分析法,计算人类活动对流域下垫面条件的影响量;设没有人类活动影响的流域下垫面系统内部状态的时间序列为{K1,K2,…,Ki,Kj,…},由系统周界A的关系式有:
Kj=f(X,Ki) (2)
式中,Ki是i时段初的状态量;Kj是即j时段初的状态量;X是系统输入;输入、状态和输出均是时间t的函数;
设人类活动-流域下垫面系统中的人类活动在i时段对系统内部状态的影响作用为Mi,由于在当前条件下人类很难对系统输入施加影响,故系统内部状态在天然的系统输入和人类活动的影响下发生着演化,故有:
Kj=f(X,MiKi) (3)
令U=MiKi,对式(3)两边求时间t的导数,有:
令分别称为输入X对系统输出Kj的弹性系数、包含有人类活动影响作用Mi的状态Ki对系统输出Kj的弹性系数;当计算的时间步长相同时,上式经整理可化简为:
从而由上式可求出i时段的人类活动对状态的影响量令其为hi,即把hi从1时段开始累积,即可得到以时间序列{K1,K2,…,Ki,Kj,…}的初始状态K1为参考点的i时段人类活动对状态的影响量Hi(即人类活动对流域下垫面条件的影响量):
从而可确定分离人类活动影响后的系统内部状态量Kni为:
Kni=Ki-Hi (7);
步骤四:利用提出的水保措施计算法计算由于人类活动对流域下垫面的影响而导致的径流变化量,从而获得人类活动对径流的间接影响量;
步骤五:在分离出人类活动对流域下垫面条件的影响量之后,确定满足下垫面条件一致性要求的人类活动-径流系统的输入、输出和内部状态要素;
步骤六:基于已构建的人类活动对径流直接影响分离的系统分析法,求出人类活动对径流的直接影响量;设人类活动-径流系统中的人类活动对系统输出的直接影响系数为Md,系统的输出取为实测径流Wm,系统的输入取为降水量、气温和蒸发气象要素Xr,系统的内部状态取为分离人类活动影响后的梯田面积、林地面积、草地面积和坝地面积Kn,则有Wm=MdWn(Xr,Kn),经分析可有:
式中,为系统输入对系统输出的弹性系数;为系统状态对系统输出的弹性系数;
因此,由式(9)可求出i时段人类活动对径流的直接影响系数Mdi,故i时段人类活动对径流的直接影响量Wdi为:
Wdi=(1-Mdi)Wmi/Mdi (10);
以及
步骤七:根据计算出来的人类活动对径流的间接影响量和人类活动对径流的直接影响量,计算出在满足流域下垫面条件一致性的天然径流量。
2.根据权利要求1所述的人类活动影响分离的二维径流还原方法,其特征在于步骤一中:基于界壳理论和泛系方法,构建水资源系统的界壳-泛系二元自动机模型,如下所示:
式中,S为水资源系统;A为水资源系统周界;B为水资源系统的系里;X为水资源系统的输入;K为水资源系统的内部状态;Y为水资源系统的输出;h为水资源系统周界对输入输出的转换机制。
3.根据权利要求2所述的人类活动影响分离的二维径流还原方法,其特征在于步骤二中:当不存在人类活动对流域的下垫面条件影响时,流域下垫面系统在系统输入及其自身的作用下发生着演化;当受人类活动影响时,形成人类活动-流域下垫面系统,系统内部状态取为流域的梯田、林地、草地和坝地四种水保措施面积,输入取为流域年平均降水量、气温和日照小时数,系统输出取为下一年度的流域的梯田、林地、草地和坝地四种水保措施面积。
4.根据权利要求3所述的人类活动影响分离的二维径流还原方法,其特征在于步骤四中:采用梯田、林地、草地和坝地四种水保措施面积变化量来表征人类活动对流域下垫面条件的影响,由上述式(6)可确定,从而由下式可求出这几种水保措施对径流的影响量,确定i时段人类活动对径流的间接影响量Wki;
式中,j=1,2,3,4,分别表示梯田、林地、草地和坝地四种水保措施;θ为地下水补给系数;di为i时段的天然径流模数,m3·km-2;Hji为i时段第j种水保措施面积,km2;uj为水保措施的利用率,%;ηj为水保措施的保水率,%。
5.根据权利要求4所述的人类活动影响分离的二维径流还原方法,其特征在于步骤五中:由上述式(7)求出获得分离人类活动影响后的天然下垫面条件变量值Kni后,系统内部状态量就在自然状态下演化,如果此时没有人类活动对径流进行直接影响时,系统的输入在这种天然状态值情况下会有一个输出,设为Wn;当有人类活动对径流施加直接影响时,形成人类活动-径流影响系统;此时,系统内部状态取为剔除人类活动影响后的梯田面积、林地面积、草地面积和坝地面积;系统输入取为流域的年平均降水量、年平均蒸发量和年平均气温;系统输出取为流域出口水文站的实测年径流量。
6.根据权利要求5所述的人类活动影响分离的二维径流还原方法,其特征在于步骤七中:据计算出来的人类活动对径流的间接影响量Wk和人类活动对径流的直接影响量Wd,从间接影响和直接影响二维角度获得在满足流域下垫面条件一致性的天然径流量W:
W=Wm+Wd+Wk (11)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410092691.2A CN103810401B (zh) | 2014-03-13 | 2014-03-13 | 一种人类活动影响分离的二维径流还原方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410092691.2A CN103810401B (zh) | 2014-03-13 | 2014-03-13 | 一种人类活动影响分离的二维径流还原方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103810401A CN103810401A (zh) | 2014-05-21 |
CN103810401B true CN103810401B (zh) | 2017-02-15 |
Family
ID=50707161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410092691.2A Expired - Fee Related CN103810401B (zh) | 2014-03-13 | 2014-03-13 | 一种人类活动影响分离的二维径流还原方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103810401B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110619111A (zh) * | 2019-09-19 | 2019-12-27 | 中国水利水电科学研究院 | 一种天然径流系列一致性修正方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105893785B (zh) * | 2014-09-04 | 2021-03-05 | 中国环境科学研究院 | 一种分离气候变化和人类活动因素对生物多样性影响的检测方法 |
US9652712B2 (en) * | 2015-07-27 | 2017-05-16 | Google Inc. | Analyzing health events using recurrent neural networks |
CN106846475B (zh) * | 2017-02-08 | 2018-03-23 | 长江水利委员会长江科学院 | 一种河流历史形态数字复原方法 |
CN110110464B (zh) * | 2019-05-15 | 2021-03-23 | 中国科学院地理科学与资源研究所 | 一种定量求解径流的流域水储量变化弹性系数的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103088783A (zh) * | 2013-01-31 | 2013-05-08 | 武汉大学 | 一种面向生态的水库调度方案生成方法 |
CN103268522A (zh) * | 2013-05-14 | 2013-08-28 | 长沙理工大学 | 径流算法 |
-
2014
- 2014-03-13 CN CN201410092691.2A patent/CN103810401B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103088783A (zh) * | 2013-01-31 | 2013-05-08 | 武汉大学 | 一种面向生态的水库调度方案生成方法 |
CN103268522A (zh) * | 2013-05-14 | 2013-08-28 | 长沙理工大学 | 径流算法 |
Non-Patent Citations (5)
Title |
---|
《基于系统周界模型的典型坝控流域人类活动影响》;李勋贵,等.;《兰州大学学报(自然科学版)》;20100228;全文 * |
Soil erosion analysis of human influence on the controlled basin system of check dams in small watersheds of the Loess Plateau;Xun-Gui Li, et al.;《Expert Systems with Applications》;20110430;全文 * |
人工神经网络在河道水量还原计算中的应用;吴建华,等.;《太原理工大学学报》;20050930;第36卷(第5期);全文 * |
水利水保措施对潮河流域年径流量的影响——基于经验统计模型的评估;李子君,等.;《地理学报》;20080930;第63卷(第9期);全文 * |
系统周界观控模型的人类活动影响定量分析;李勋贵,等.;《系统工程学报》;20121031;第27卷(第5期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110619111A (zh) * | 2019-09-19 | 2019-12-27 | 中国水利水电科学研究院 | 一种天然径流系列一致性修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103810401A (zh) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zuo et al. | Simulating spatiotemporal variability of blue and green water resources availability with uncertainty analysis | |
CN103810401B (zh) | 一种人类活动影响分离的二维径流还原方法 | |
Ye et al. | A modeling study of catchment discharge to Poyang Lake under future climate in China | |
CN104318077B (zh) | 气候变化和人类活动对河川径流变化定量分析方法 | |
Guan et al. | Study on spatiotemporal distribution characteristics of flood and drought disaster impacts on agriculture in China | |
CN106845771A (zh) | 一种基于前期雨量优选参数的洪水预报方法 | |
CN108614915B (zh) | 基于情景驱动的水文模型自由组建策略方法 | |
Min et al. | Quantitative simulation of the relationships between cultivated land-use patterns and non-point source pollutant loads at a township scale in Chaohu Lake Basin, China | |
Li et al. | Human and climatic drivers of land and water use from 1997 to 2019 in Tarim River basin, China | |
CN111199298B (zh) | 基于神经网络的洪水预报方法与系统 | |
Wang et al. | Long-term variations in water balance by increase in percent imperviousness of urban regions | |
Hao et al. | Coupling the Xinanjiang model with geomorphologic instantaneous unit hydrograph for flood forecasting in northeast China | |
Liu et al. | Impacts of land use/cover change on water balance by using the SWAT model in a typical loess hilly watershed of China | |
CN109979172A (zh) | 一种基于新安江模型的动态山洪临界雨量预报方法 | |
Long et al. | Improving streamflow simulation in Dongting Lake Basin by coupling hydrological and hydrodynamic models and considering water yields in data-scarce areas | |
Chen et al. | Identify runoff generation patterns of check dams and terraces and the effects on runoff: a case study | |
CN106446359A (zh) | 基于流式大数据处理方式的降雨径流预测计算方法 | |
Yang et al. | Research on intelligent prediction and zonation of basin-scale flood risk based on LSTM method | |
Lin et al. | Improving calibration of two key parameters in Hydrologic Engineering Center hydrologic modelling system, and analysing the influence of initial loss on flood peak flows | |
Zhou et al. | Future runoff forecast in Hanjiang River Basin based on Wetspa model and CMIP6 model | |
Teng et al. | Early warning index of flash flood disaster: a case study of Shuyuan watershed in Qufu City | |
Zhu et al. | Flood forecasting methods for a semi‐arid and semi‐humid area in Northern China | |
Li et al. | Long-lead daily streamflow forecasting using Long Short-Term Memory model with different predictors | |
Fu et al. | Water Resource Availability Assessment Through Hydrological Simulation Under Climate Change in the Huangshui Watershed of the Qinghai–Tibet Plateau | |
CN104794539A (zh) | 一种洪涝损失模型构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 |