CN113742636B - 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 - Google Patents
一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 Download PDFInfo
- Publication number
- CN113742636B CN113742636B CN202110731130.2A CN202110731130A CN113742636B CN 113742636 B CN113742636 B CN 113742636B CN 202110731130 A CN202110731130 A CN 202110731130A CN 113742636 B CN113742636 B CN 113742636B
- Authority
- CN
- China
- Prior art keywords
- evaporation
- actual
- actual evaporation
- vegetation
- potential
- 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.)
- Active
Links
- 230000008020 evaporation Effects 0.000 title claims abstract description 98
- 238000001704 evaporation Methods 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 56
- 230000014509 gene expression Effects 0.000 claims abstract description 14
- 239000002689 soil Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000003673 groundwater Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 7
- 230000008859 change Effects 0.000 abstract description 6
- 230000008569 process Effects 0.000 description 6
- 230000007774 longterm Effects 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000005068 transpiration Effects 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 230000002354 daily effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 241000235935 Hilaria belangeri Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
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
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)。
步骤3具体为,根据步骤2中计算出的年潜在蒸散发值,在选定序列长度m的条件下,计算多年平均实际蒸散发值如下公式:
其中,m为序列长度,EAi为第i年实际蒸散发值,为m年平均实际散发值。
步骤4中,水量平衡原理的理论基础是基于质量守恒原理,其是全面研究单个地区在一定时间段内地下水的补给量、储存量和消耗量之间的数量转化关系的平衡计算;由此可得,多年实际蒸散发值为:
根据由水量平衡原理和潜在蒸散发计算出的多年实际蒸散发值应相等,即可得到:
根据公式(8)可率定出水分系数Kw和植被指数Kv的函数表达式。
步骤5中,根据率定出的水分系数和植被指数表达式,可直接依据公式 (2)计算实际蒸散发值。
本发明的有益效果是:本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,基于水量平衡原理,提出了一种利用水分系数及植被指数关系推求实际蒸散发的方法。在实际蒸散发的估算研究中,相比于利用模型模拟的方法,能够更加精确获得水分及植被变化下蒸散发的变化,去除了模型不确定性的影响,为实际蒸散发的估算提供了有力支撑。
附图说明
图1是本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法的流程示意图。
具体实施方式
下面结合附图和具体实施方式对本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法进行进一步详细说明。
如图1所示,一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,根据彭曼公式计算潜在蒸散发,再利用水量平衡原理,率定出水分系数以及植被指数与降水量等因子的关系,从而直接计算某一流域的实际蒸散发,具体按照以下步骤实施:
步骤1,基于彭曼公式计算出逐日潜在蒸散发Epmd;
步骤2,假设年实际蒸发
其中,Kw为水分系数,Kv为植被指数;
步骤3,多年平均实际蒸散发
步骤4,基于水量平衡得根据/>率定出Kw, Kv中的参数;
步骤4.1,假定长期水量平衡后Δω=0由于水文模拟的Δω不准确,所以不直接用水量平衡计算实际蒸发EA,而水量平衡得/>
步骤4.2,水分系数kw=f(P,Psummer,…)为年降水量等变量的函数,例如 Kw=a0+a1P+a2Psummer或其他形式;Kv是植被覆盖指数Vc的函数,例如 Kv=b0+b1Vc+b2Vc或其他形式;根据由水量平衡原理和潜在蒸散发计算出的多年实际蒸散发值应相等即可得到Kw、Kv的函数表达式。
步骤5,利用公式直接计算某一区域或者流域的实际蒸散发。
进一步地,步骤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的函数表达式;
步骤5,即为利用公式(2)计算实际蒸散发。
下面通过具体的实施例对本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法进行进一步详细说明。
步骤1,计算潜在蒸散发:
基于彭曼公式计算潜在蒸散发,如图1所示。
步骤2:基于潜在蒸散发计算实际蒸散发值:
水分供应不充分时,陆面实际蒸散发与潜在蒸散发成正比,其大小取决于水分的有效性,即将实际蒸散发表征为潜在蒸散发的函数,采用土壤水分含量、植被生长状况等作为限制因子。即:
其中,公式(3)Kw=f(p,psummer) (3)的形式从以下公式中选择拟合效果最好的公式:
公式(4)Kv=f(Vc)的形式从以下公式中选择拟合效果最好的公式:
步骤3:计算多年平均的潜在蒸散发
步骤4,确定Kw、Kv的函数形式:
步骤4.1,根据水量平衡原理计算实际蒸散发值;
假定长期水量平衡后Δω=0由于水文模拟的Δω不准确,所以不直接用水量平衡计算实际蒸发EA,而水量平衡得到/>
步骤4.2,对比步骤3中计算的实际蒸散发结果,
由此得到,水分系数Kw以及植被指数Kv的函数表达式;
步骤5,根据水分系数Kw以及植被指数Kv表达式由公式(2)可直接计算实际蒸散发。
本发明一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,解决了现有技术中通过水文模拟计算实际蒸散发的不确定性问题。相比于利用模型模拟的方法,能够更加精确获得水分及植被变化下蒸散发的变化,去除了模型不确定性的影响,为实际蒸散发的估算提供了有力支撑。
Claims (6)
1.一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,具体按照以下步骤实施:
步骤1:利用彭曼公式计算潜在蒸散发值;
步骤2:基于潜在蒸散发值来计算实际蒸散发值;
步骤3:计算多年平均实际蒸散发;
步骤4:率定出水分系数和植被指数中的参数;基于水量平衡原理得到实际蒸散发;与步骤3中的实际蒸散发值作比较,得到水分系数和植被指数的表达式;
步骤5:利用步骤4得到的表达式,计算单个区域的实际蒸散发值。
2.根据权利要求1所述的一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,步骤1中,所述彭曼公式计算日潜在蒸散发值是基于逐日最高气温、最低气温、平均气温、平均水汽压、平均风速和日照时数的资料,具体如下公式(1):
3.根据权利要求1所述的一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,步骤2具体为:依据水分供应不充分时,实际蒸散发与潜在蒸散发成正比的关系,即受到土壤供水能力和植被生长情况的影响,得到如下公式(2)的实际蒸散发的计算公式:
其中,Kw为水分系数,Kw为植被指数。
4.根据权利要求1所述的一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,所述步骤3具体为,根据步骤2中计算出的年潜在蒸散发值,在选定序列长度m的条件下,计算多年平均潜在蒸散发值如下公式:
其中,m为序列长度,EAi为第i年实际蒸散发值,为m年平均实际散发值。
5.根据权利要求1所述的一种基于气象观测数据和植被指数的流域实际蒸发量估算方法,其特征在于,步骤4中,水量平衡原理的理论基础是基于质量守恒原理,其是全面研究单个地区在一定时间段内地下水的补给量、储存量和消耗量之间的数量转化关系的平衡计算,因此,可由水量平衡计算多年实际蒸散发:
根据由水量平衡原理和潜在蒸散发计算出的多年实际蒸散发值应相等,即可得到:
进而即可率定出水分系数Kw和植被指数Kv的函数表达式。
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 CN113742636A (zh) | 2021-12-03 |
CN113742636B true 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,刘登峰,明广辉,韩松俊,孟宪萌,黄强;《人民珠江》;20240105;全文 * |
气象干旱向水文干旱传播的如影杻汲其驱动力研究;马岚;《中国优秀硕士论文全文数据库》;20190815;全文 * |
海流兔河流域地下水蒸散发量估算研究;赵楠;《中国优秀硕士学位论文全文数据库》;20230115;全文 * |
陕西省多年月潜在蒸发量 和降水量的随机特性;牛文全,李 靖;《干旱地区农业研究》;20010731;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113742636A (zh) | 2021-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jones et al. | Future climate impact on the productivity of sugar beet (Beta vulgaris L.) in Europe | |
Li et al. | Carbon dioxide exchange and the mechanism of environmental control in a farmland ecosystem in North China Plain | |
Zhao et al. | Canopy water use efficiency of winter wheat in the North China Plain | |
Adhikari et al. | Simulated effects of winter wheat cover crop on cotton production systems of the Texas Rolling Plains | |
Liu et al. | Optimizing irrigation rates for cotton production in an extremely arid area using RZWQM2-simulated water stress | |
CN110533287A (zh) | 夏玉米干旱灾损定量测评模型的构建与应用 | |
Chatterjee et al. | Actual evapotranspiration and crop coefficients for tropical lowland rice (Oryza sativa L.) in eastern India | |
Chuanyan et al. | Estimating water needs of maize (Zea mays L.) using the dual crop coefficient method in the arid region of northwestern China | |
CN113742636B (zh) | 一种基于气象观测数据和植被指数的流域实际蒸发量估算方法 | |
Sarma et al. | Determination of crop-coefficients and estimation of evapotranspiration of rapeseed using lysimeter and different reference evapotranspiration models | |
Meshram et al. | Water requirement of pomegranate (Punica granatum L.) for Solapur district of Maharashtra State | |
Huntington et al. | Evapotranspiration and net irrigation water requirements for Nevada | |
Treder et al. | An hourly reference evapotranspiration model as a tool for estimating plant water requirements | |
CN103218532B (zh) | 绿水流估算方法及装置 | |
Radersma et al. | Water use of tree lines: importance of leaf area and micrometeorology in sub-humid Kenya | |
Chaowiwat et al. | Prediction of Future Agriculture Water Demand in Thailand Using Multi Bias Corrected Climate Models | |
Vishnoi et al. | Study on water requirement of potato (Solanum tuberosum L.) using CROPWAT MODEL for Tarai region of Uttarakhand | |
Hu et al. | A study on crop coefficients of jujube under drip-irrigation in Loess Plateau of China | |
Bello et al. | Evaluation of some empirical methods of estimating potential evapotranspiration for determination of length of growing season in a tropical wet and dry climate | |
Chatterjee et al. | Determination of Actual Evapotranspiration and Crop Coefficients of Tropical Indian Lowland Rice (Oryza sativa) Using Eddy Covariance Approach | |
Olejnik et al. | Mitigation of radiation and heat balance structure by plant cover structure | |
CN117291047A (zh) | 一种基于动态冠层消光系数的作物蒸腾耗水量模拟方法 | |
CN106034973B (zh) | 一种杉木人工林防灾控制方法 | |
Mei et al. | A Study on Variation of Priestley-Taylor Model Parameter in Rain Fed Spring Maize Field | |
CN117474704A (zh) | 一种玉米昼夜生长差异和产量的估算方法 |
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 |