CN113742636A - 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 - Google Patents
一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 Download PDFInfo
- Publication number
- CN113742636A CN113742636A CN202110731130.2A CN202110731130A CN113742636A CN 113742636 A CN113742636 A CN 113742636A CN 202110731130 A CN202110731130 A CN 202110731130A CN 113742636 A CN113742636 A CN 113742636A
- Authority
- CN
- China
- Prior art keywords
- evapotranspiration
- actual
- water
- evaporation capacity
- vegetation index
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000008020 evaporation Effects 0.000 title claims abstract description 35
- 238000001704 evaporation Methods 0.000 title claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 61
- 239000002689 soil Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000003673 groundwater Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 10
- 230000003993 interaction Effects 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 241000235935 Hilaria belangeri Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000012223 aqueous fraction Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- GMTYREVWZXJPLF-AFHUBHILSA-N butorphanol D-tartrate Chemical compound OC(=O)[C@@H](O)[C@H](O)C(O)=O.N1([C@@H]2CC3=CC=C(C=C3[C@@]3([C@]2(CCCC3)O)CC1)O)CC1CCC1 GMTYREVWZXJPLF-AFHUBHILSA-N 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明公开了一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,步骤1,计算潜在蒸散发;步骤2,通过水分系数和植被指数建立实际蒸散发和潜在蒸散发之间的函数关系;步骤3,计算多年平均蒸散发值;步骤4,率定Kw,Kv表达式,步骤4:率定出水分系数和植被指数中的参数;步骤5,利用步骤4得到的表达式,计算单个区域的实际蒸散发值。本发明解决了现有技术中通过水文模拟计算实际蒸散发的不确定性问题。相比于利用模型模拟的方法,能够更加精确获得水分及植被变化下蒸散发的变化,去除了模型不确定性的影响,为实际蒸散发的估算提供了有力支撑。
Description
技术领域
本发明属于实际蒸散发估算技术领域,具体涉及一种基于气象观测数据和植被指数的流域实际蒸发量估算方法。
背景技术
蒸腾与蒸发量,简称蒸散发(evapotranspiration,ET)作为区域水量平衡和能量平衡的主要组分,不仅在水循环和能量循环过程中起着极其重要的作用,也是连接生态与水文过程的重要纽带,更是评价农业用水效率的重要基础和关键环节。它决定着地球系统中“地圈—生物圈—大气圈”的相互作用和反馈。目前,人们对自由水面、饱和土壤表面和充分供水条件下的稠密植被冠层等下垫面蒸散发理解较为深入,将其称为潜在蒸散发。实际蒸散发量一般根据潜在蒸散发进行估算,主要理论方法有两种:Penman方法和互补原理。Penman方法从田间尺度蒸散发过程出发,将潜在蒸散发作为大气对陆面蒸散发过程的外部驱动力,根据实测或模拟的陆面水分状况估算实际蒸散发。互补原理则直接从区域尺度蒸散发过程出发,将潜在蒸散发作为内部变量,根据潜在蒸散发对实际蒸散发过程的反馈来估算实际蒸散发。传统的测定方法上具有较高的精度,但主要针对局地尺度和点状数据,难以实现区域大尺度的扩展。而水文模型对于实际蒸散发的模拟又存在很大不确定性的影响。
发明内容
本发明的目的是提供一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,解决了水文模型对于实际蒸散发的模拟存在较大不确定性的问题。
本发明所采用的技术方案是,
一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,具体按照以下步骤实施:
步骤1:利用彭曼公式计算潜在蒸散发值;
步骤2:基于潜在蒸散发值来计算实际蒸散发值;
步骤3:计算多年平均实际蒸散发;
步骤4:率定出水分系数和植被指数中的参数;基于水量平衡原理得到实际蒸散发;与步骤3中的实际蒸散发值作比较,得到水分系数和植被指数的表达式;
步骤5:利用步骤4得到的表达式,计算单个区域的实际蒸散发值。
本发明的特点还在于,
步骤1中,彭曼公式计算日潜在蒸散发值是基于逐日最高气温、最低气温、平均气温、平均水汽压、平均风速和日照时数的资料:
其中:Rn为作物表面净辐射量,MJ/(m2·d);G为土壤热通量,MJ/(m2·d);γ为湿度计常数,kPa/℃;Δ为饱和水汽压与温度关系曲线的斜率kPa/℃;T 为空气平均温度℃;u2为在地面以上2m处的风速m/s;es为空气饱和水汽压 kPa;ea为空气实际水压kPa。
步骤2具体为,依据水分供应不充分时,实际蒸散发与潜在蒸散发成正比的关系,即受到土壤供水能力和植被生长情况的影响,得到如下的实际蒸散发的计算公式:
其中,EA为年实际蒸散发值,Epmd为日潜在蒸散发之,Kw水分系数,是降水量p的函数,即:
Kw=f(p,psummer) (3),
Kv为植被指数,是指被覆盖指数Vc的函数,即
Kv=f(Vc) (4)。
步骤4中,水量平衡原理的理论基础是基于质量守恒原理,其是全面研究单个地区在一定时间段内地下水的补给量、储存量和消耗量之间的数量转化关系的平衡计算;由此可得,多年实际蒸散发值为:
根据由水量平衡原理和潜在蒸散发计算出的多年实际蒸散发值应相等,即可得到:
根据公式(8)可率定出水分系数Kw和植被指数Kv的函数表达式。
步骤5中,根据率定出的水分系数和植被指数表达式,可直接依据公式 (2)计算实际蒸散发值。
本发明的有益效果是:本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,基于水量平衡原理,提出了一种利用水分系数及植被指数关系推求实际蒸散发的方法。在实际蒸散发的估算研究中,相比于利用模型模拟的方法,能够更加精确获得水分及植被变化下蒸散发的变化,去除了模型不确定性的影响,为实际蒸散发的估算提供了有力支撑。
附图说明
图1是本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法的流程示意图。
具体实施方式
下面结合附图和具体实施方式对本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法进行进一步详细说明。
如图1所示,一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,根据彭曼公式计算潜在蒸散发,再利用水量平衡原理,率定出水分系数以及植被指数与降水量等因子的关系,从而直接计算某一流域的实际蒸散发,具体按照以下步骤实施:
步骤1,基于彭曼公式计算出逐日潜在蒸散发Epmd;
其中,Kw为水分系数,Kv为植被指数;
步骤4.2,水分系数kw=f(P,Psummer,…)为年降水量等变量的函数,例如 Kw=a0+a1P+a2Psummer或其他形式;Kv是植被覆盖指数Vc的函数,例如 Kv=b0+b1Vc+b2Vc或其他形式;根据由水量平衡原理和潜在蒸散发计算出的多年实际蒸散发值应相等即可得到Kw、Kv的函数表达式。
进一步地,步骤1具体为,计算逐日的潜在蒸散发值Epmd。
Penman(1948)第一次将能量平衡与空气动力学方法结合,并巧妙利用饱和水汽压-温度曲线斜率的概念,提出了适用于宽阔自由水面、湿润裸土和短草的蒸发解析式,之后Monteith(1965)在Penman公式的基础上,考虑了植被生理特征,引入冠层阻力和空气动力学阻力,仍然保持Penman公式的形式,建立了Penman-Montheith公式。该公式的自变量包含了辐射、温度、风速、湿度、和气压等气象数据,可分为辐射项和空气动力项两部分。综合了空气动力学的湍流传输与能量平衡,较为全面地考虑了影响潜在蒸散的气象因素,是目前较为普遍的计算潜在蒸散发的一种方法。
Penman-Montheith公式为:
进一步地,步骤2具体为,在蒸散发的研究中,实际蒸散发(EA)与潜在蒸散发(Epmd)之间的关系对于分析气候与水文相互作用关系具有非常重要的作用。当水分供应不充分时,陆面实际蒸散发与潜在蒸散发成正比,其大小取决于水分的有效性,即将实际蒸散发表征为潜在蒸散发的函数,采用土壤水分含量、植被生长状况等作为限制因子。Penman正比假设理论可以表示为EA=KwEpmd,Kw为水分系数,表示为土壤供水能力函数,在有植被或作物生长情况下,EA=KwKvEpmd,Kv为植被指数,表示植被生长函数(作物系数)。因此,建立基于水分系数和植被指数的实际蒸散发与潜在蒸散发的函数关系。
进一步地,步骤3具体为,根据步骤2中计算出的年潜在蒸散发EA,计算多年平均的潜在蒸散发:
进一步地,步骤4具体为,水量平衡原理是全面研究某一地区(或均衡区)在一定时间段内地下水的补给量、储存量和消耗量之间的数量转化关系的平衡计算,其理论基础是质量守恒原理。假定长期水量平衡后Δω=0 由于水文模拟的Δω不准确,所以无法直接用水量平衡计算实际蒸发EA,而由水量平衡可得:
进而得到:
既而得到水分系数Kw以及植被指数Kv的函数表达式;
下面通过具体的实施例对本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法进行进一步详细说明。
步骤1,计算潜在蒸散发:
基于彭曼公式计算潜在蒸散发,如图1所示。
步骤2:基于潜在蒸散发计算实际蒸散发值:
水分供应不充分时,陆面实际蒸散发与潜在蒸散发成正比,其大小取决于水分的有效性,即将实际蒸散发表征为潜在蒸散发的函数,采用土壤水分含量、植被生长状况等作为限制因子。即:
其中,公式(3)Kw=f(p,psummer) (3)的形式从以下公式中选择拟合效果最好的公式:
公式(4)Kv=f(Vc)的形式从以下公式中选择拟合效果最好的公式:
步骤4,确定Kw、Kv的函数形式:
步骤4.1,根据水量平衡原理计算实际蒸散发值;
步骤4.2,对比步骤3中计算的实际蒸散发结果,
由此得到,水分系数Kw以及植被指数Kv的函数表达式;
步骤5,根据水分系数Kw以及植被指数Kv表达式由公式(2)可直接计算实际蒸散发。
本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,解决了现有技术中通过水文模拟计算实际蒸散发的不确定性问题。相比于利用模型模拟的方法,能够更加精确获得水分及植被变化下蒸散发的变化,去除了模型不确定性的影响,为实际蒸散发的估算提供了有力支撑。
Claims (6)
1.一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,具体按照以下步骤实施:
步骤1:利用彭曼公式计算潜在蒸散发值;
步骤2:基于潜在蒸散发值来计算实际蒸散发值;
步骤3:计算多年平均实际蒸散发;
步骤4:率定出水分系数和植被指数中的参数;基于水量平衡原理得到实际蒸散发;与步骤3中的实际蒸散发值作比较,得到水分系数和植被指数的表达式;
步骤5:利用步骤4得到的表达式,计算单个区域的实际蒸散发值。
6.根据权利要求3所述的一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,根据步骤4中率定出的水分系数和植被指数表达式,可由公式(2)直接计算实际蒸散发值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110731130.2A CN113742636B (zh) | 2021-06-29 | 2021-06-29 | 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110731130.2A CN113742636B (zh) | 2021-06-29 | 2021-06-29 | 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113742636A true CN113742636A (zh) | 2021-12-03 |
CN113742636B CN113742636B (zh) | 2024-04-19 |
Family
ID=78728561
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110731130.2A Active CN113742636B (zh) | 2021-06-29 | 2021-06-29 | 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742636B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009025122A (ja) * | 2007-07-19 | 2009-02-05 | Tokyo Electric Power Co Inc:The | 洞道内の潜熱効果を考慮したケーブル導体温度推定方法、ケーブル導体温度推定システム、及びケーブル導体温度推定プログラム |
CN101551459A (zh) * | 2008-10-15 | 2009-10-07 | 北京天宏金睛信息技术有限公司 | 基于遥感的区域蒸散量监测方法 |
CN102097056A (zh) * | 2009-12-11 | 2011-06-15 | 宏达国际电子股份有限公司 | 电子装置及其控制背光亮度的方法 |
CN107084709A (zh) * | 2017-04-17 | 2017-08-22 | 河海大学 | 一种多弹性对径流变化驱动因素的定量分割方法 |
CN107258494A (zh) * | 2017-07-06 | 2017-10-20 | 郑州大学 | 一种植被生态高效用水智能化精准调配系统 |
CN107818238A (zh) * | 2017-09-28 | 2018-03-20 | 河海大学 | 一种确定蒸散发变化主因及判别因素间耦合作用的方法 |
CN108984867A (zh) * | 2018-06-28 | 2018-12-11 | 中国科学院地理科学与资源研究所 | 基于gf-4和modis结合的蒸散发遥感反演方法和系统 |
CN110727900A (zh) * | 2019-09-20 | 2020-01-24 | 中国科学院遥感与数字地球研究所 | 一种流域植被干旱发生遥感预警与缺水估算方法 |
CN112649370A (zh) * | 2019-11-13 | 2021-04-13 | 四川大学 | 一种基于遥感的区域蒸散发计算方法 |
CN112800636A (zh) * | 2021-04-08 | 2021-05-14 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种估算无资料地区流域地表水资源量的方法及系统 |
-
2021
- 2021-06-29 CN CN202110731130.2A patent/CN113742636B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009025122A (ja) * | 2007-07-19 | 2009-02-05 | Tokyo Electric Power Co Inc:The | 洞道内の潜熱効果を考慮したケーブル導体温度推定方法、ケーブル導体温度推定システム、及びケーブル導体温度推定プログラム |
CN101551459A (zh) * | 2008-10-15 | 2009-10-07 | 北京天宏金睛信息技术有限公司 | 基于遥感的区域蒸散量监测方法 |
CN102097056A (zh) * | 2009-12-11 | 2011-06-15 | 宏达国际电子股份有限公司 | 电子装置及其控制背光亮度的方法 |
CN107084709A (zh) * | 2017-04-17 | 2017-08-22 | 河海大学 | 一种多弹性对径流变化驱动因素的定量分割方法 |
CN107258494A (zh) * | 2017-07-06 | 2017-10-20 | 郑州大学 | 一种植被生态高效用水智能化精准调配系统 |
CN107818238A (zh) * | 2017-09-28 | 2018-03-20 | 河海大学 | 一种确定蒸散发变化主因及判别因素间耦合作用的方法 |
CN108984867A (zh) * | 2018-06-28 | 2018-12-11 | 中国科学院地理科学与资源研究所 | 基于gf-4和modis结合的蒸散发遥感反演方法和系统 |
CN110727900A (zh) * | 2019-09-20 | 2020-01-24 | 中国科学院遥感与数字地球研究所 | 一种流域植被干旱发生遥感预警与缺水估算方法 |
CN112649370A (zh) * | 2019-11-13 | 2021-04-13 | 四川大学 | 一种基于遥感的区域蒸散发计算方法 |
CN112800636A (zh) * | 2021-04-08 | 2021-05-14 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种估算无资料地区流域地表水资源量的方法及系统 |
Non-Patent Citations (4)
Title |
---|
李彤彤1,刘登峰,明广辉,韩松俊,孟宪萌,黄强: "基于多种方法的汉江和渭河流域实际蒸发量对比研究", 《人民珠江》, 5 January 2024 (2024-01-05) * |
牛文全,李 靖: "陕西省多年月潜在蒸发量 和降水量的随机特性", 《干旱地区农业研究》, 31 July 2001 (2001-07-31) * |
赵楠: "海流兔河流域地下水蒸散发量估算研究", 《中国优秀硕士学位论文全文数据库》, 15 January 2023 (2023-01-15) * |
马岚: "气象干旱向水文干旱传播的如影杻汲其驱动力研究", 《中国优秀硕士论文全文数据库》, 15 August 2019 (2019-08-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN113742636B (zh) | 2024-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hao et al. | Seasonal and interannual variation in water vapor and energy exchange over a typical steppe in Inner Mongolia, China | |
Arora et al. | A representation of variable root distribution in dynamic vegetation models | |
CN106771073B (zh) | 一种基于端元信息模型估算土壤和植被蒸散发的方法 | |
Mamadou et al. | Energy fluxes and surface characteristics over a cultivated area in Benin: daily and seasonal dynamics | |
Lindroth | Aerodynamic and canopy resistance of short-rotation forest in relation to leaf area index and climate | |
Tourula et al. | Modelling evapotranspiration from a barley field over the growing season | |
Raddatz | Evidence for the influence of agriculture on weather and climate through the transformation and management of vegetation: Illustrated by examples from the Canadian Prairies | |
Shuttleworth | Towards one-step estimation of crop water requirements | |
CN112733385B (zh) | 基于时变增益模型的水、能、碳耦合数值模拟方法及系统 | |
CN112699621A (zh) | 一种考虑下垫面条件变化的潜在蒸散发估算方法 | |
Arain et al. | Testing of vegetation parameter aggregation rules applicable to the Biosphere-Atmosphere Transfer Scheme (BATS) and the FIFE site | |
Ritchie et al. | Calculating evaporation from native grassland watersheds | |
Bailey et al. | The effect of uncertainty in aerodynamic resistance on evaporation estimates from the combination model | |
CN109188465B (zh) | 基于参考像元信息的区域地表蒸散发遥感估算方法 | |
Shang et al. | Simulation of water dynamics and irrigation scheduling for winter wheat and maize in seasonal frost areas | |
Liu et al. | An improved evapotranspiration model for an apple orchard in northwestern China | |
CN109858186B (zh) | 联合光学和热学数据的农田干旱评估方法 | |
Harazono et al. | Characteristics of energy and water budgets over wet sedge and tussock tundra ecosystems at North Slope in Alaska | |
CN113742636A (zh) | 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 | |
Pitacco et al. | Evaluation of actual evapotranspiration of a Quercus ilex L. stand by the Bowen ratio-energy budget method | |
QIU et al. | Theoretical analysis and experimental verification of a remotely measurable plant transpiration transfer coefficient | |
Duffková | Difference in canopy and air temperature as an indicator of grassland water stress | |
CN110083954B (zh) | 基于边界层阻力测量技术的温室作物需水量计算方法 | |
Ling et al. | Energy balance and evapotranspiration characteristics of rubber tree (Hevea brasiliensis) plantations in Xishuangbanna, Southwest of China. | |
Tesfuhuney et al. | Comparison of energy available for evapotranspiration under in‐field rainwater harvesting with wide and narrow runoff strips |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |