作物蒸发蒸腾量计算公式
一、采用彭曼—蒙蒂斯(Penman—Monteith)法计算参考作物蒸发蒸腾量(ET0)
1、彭曼—蒙蒂斯(Penman—Monteith)公式
彭曼—蒙蒂斯(Penman—Monteith)公式是联合国粮农组织(FAO,1998)提出的最新修正彭曼公式,并已被广泛应用且已证实具有较高精度及可使用性。P-M公式对参考作物的蒸发蒸腾量定义如下:参考作物的蒸发蒸腾量为一种假想的参考作物冠层的蒸发蒸腾速率,假想作物的高度为0.12m,固定的叶面阻力为70s/m,反射率为0.23,非常类似于表面开阔、高度一致、生长旺盛、完全覆盖地面且不缺水的绿色草地蒸发蒸腾量。
Penman——Monteith公式:
ET00.408(RnG)900U2(eaed)T273(10.34U2) (1)
式中 ET0——参考作物蒸发蒸腾量,mm/d;
——温度~饱和水汽压关系曲线在T处的切线斜率,kPa∙℃-1;
4098ea(T237.3)2 (2)
T——平均气温,℃
ea——饱和水汽压,kpa;
.27Tea0.611expT17237.3 (3)
Rn——净辐射,MJ/(m2∙d);
Rns——净短波辐射,MJ/(m2∙d);Rnl——净长波辐射,MJ/(m2∙d);n——实际日照时数,h;
N——最大可能日照时数,h;
Ws——日照时数角,rad;
ψ——地理纬度,rad;
RnRnsRnl Rns0.77(0.250.5n/N)Ra N7.64Ws Wsarccos(tantan) 4)
5)
6)
7)
(
( ( (
δ——日倾角,rad;
0.409sin(0.0172J1.39) (8)
J——日序数(元月1日为1,逐日累加);
Ra——大气边缘太阳辐射,MJ/(m2∙d);
Ra37.6dr(WssinsincoscossinWs)dr——日地相对距离;
d0.033cos(2πr1365J) R2.45109(0.9n/N0.1)(0.340.14e44nld)(TkxTkn) ed——实际水汽压,kpa;
eTmin)ed(Tmax)1RHmaxded(22e(T1RHminamin)1002ea(Tmax)100 RHmax——日最大相对湿度,%;
Tmin——日最低气温;℃
ea(Tmin)——Tmin时饱和水汽压,kpa,可将Tmin代入(3)式求得;
(9)
10)
11)
12)
( ( (
ed(Tmin)——Tmin时实际水汽压,kpa;
RHmin——日最小相对湿度,%;
Tmax——日最高气温,℃
ea(Tmax)——Tmax时饱和水汽压,kpa,可将Tmax代入(3)式求得;
ed(Tmax)——Tmax时实际水汽压,kpa;
若资料不符合(12)式要求或计算较长时段ET0,也可采用下式计算ed,即e5050dRHmean/ea(Tmin)ea(Tmax) RHmean——平均相对湿度,%;
RHmeanRHmaxRHmin2 在最低气温等于或十分接近露点温度时,也可采用下式计算ed,即
e0.611exp17.27TmindTmin237.3 Tks——最高绝对温度,K;
Tkn——最低绝对温度,K;
13)14)15)
(
(
(
TksTmax273 (16)
TknTmin273 (17)
G——土壤热通量,MJ/(m2∙d);
对于逐日估算ET0,则第d日土壤热通量为
G0.38(TdTd1) 对于分月估算ET0,则第m月土壤热通量为:
G0.14(TmTm1) Td、Td-1——分别为第d、d-1日气温,℃;
Tm、Tm-1——分别为第m、m-1日气温,℃;
γ——湿度表常数,kpa∙℃-1;
0.00163P/ P——气压,kpa;
P101.3(2930.0065Z5.26293) 18)
19)
20)
21)
( ( ( (
Z——计算地点海拔高程,m;
λ——潜热,MJ∙kg-1;
2.501(2.361103)T (35)
u2——2m高处风速,m/s;
u24.87uh/ln(67.8h5.42) (36)
h——风标高度,m;
uh——实际风速,m/s。
2、实际作物需水计算,根据试验测定的蒸发蒸腾量和计算的参考作物蒸发蒸腾量,分析确定作物系数(Kc)和土壤水分修正系数(Ks)后,可计算实际作物需水量 。
ETcKcET0 ETaKsET0
Kc——作物系数; Ks——土壤水分修正系数; ET0——参考作物腾发量;ETc——作
物需水量;ETa——任意土壤水分条件下的作物腾发量。
二、波文比-能量平衡法
1926年Bowen从能量平衡方程出发,提出了计算水面蒸发的波文比-能量平衡模型。该方法的两大理论支柱是能量平衡原理和边界层扩散理论。假定植物和土壤是一个蒸发界
面,水分子可以从此界面逸出而进入大气,那么,对于这个面的垂直方向上的能量收支平衡可用下式描述:
RnGHET (1)
式中Rn—太阳净辐射;
G—土壤热通量;
H—感热通量;
ET—潜热通量,—水汽化潜热,ET—植物蒸发蒸腾量。
波文比定义为
HET (2)
综合式(1)和(2)可得:
ETRnG1 (3)
式(3)即为用波文比-能量平衡法估算植物蒸发蒸腾量的公式,其关键在于波文比的确定。
根据经验关系,感热通量、潜热通量可表示为:
HaCpkhTz (4)
ETaCpkvez (5)
式中a—空气密度;
CP—空气定压比热;
kh—感热交换系数;
kv—潜热交换系数;
—湿度计常数,约为0.66。
根据雷诺相似原理,假定感热和潜热的交换系数相等,即kh=kv,合并式(2)、式(4)和式(5)可得:
T/zTe/ze (6)
式中ΔT—上下空气温度差(℃);
Δe—上下饱和气压差。
ETRnGT(1)e (7)
利用波文比系统测得Rn,G,T和e后,就能够计算出该区域的潜热通量和相应的植物蒸发蒸腾量。
波文比-能量平衡法素以物理概念明确、计算方法简单而著称,且对大气层没有特别的要求和限制。该法只需要两个高度的要素观测值,不用求湍流交换系数,而且精度较高,可作为其他蒸发蒸腾量测定方法的准判别标准。但是,使用波文比系统观测的区域要具有开阔、均一的下垫面,且天气平稳少变,辐射和风速都没有过于剧烈的变化。该模型长期以来得到了较好的应用,但在下垫面极为潮湿或平流逆温条件下,计算结果偏低,精度下降。
三、彭曼综合法公式
Penman修正公式根据能量平衡原理、水汽扩散原理和空气导热定律等提出并经多次修正的,1979年联合国粮农组织(FAO)推荐的修正公式为:
P0Rn0.26(esea)(1Cu2)PET0P01.0P
式中P0—海平面平均气压,hPa;
P—计算点平均气压,hPa;
△为饱和水汽压—温度曲线上的斜率(mbar/℃);
—湿度计常数,约为0.66;
ea —空气中实际大气压(mbar);
es —饱和水汽压(mbar);
u2 —2 m处的风速,若用气象站常规的观测高度的风速则需乘以0.75的风速修正系数,m/s;
C—与最高气温和最低气温有关的风速修正系数。
四、布莱尼—克雷多公式
ET0c[p(0.46t8)]
式中ET0—月平均参考作物蒸发蒸腾量(mm/d);
t—平均气温(℃);
p—月内日平均可能日照时数占全年可能日照时数的百分比; c—根据最低相对湿度、日照时数、白天风速确定的修正系数。
五、水汽扩散公式
利用近地面大气层中的湍流交换规律为基础,其形式为:
(u2u1)(q2q1)[ln(z2/z1)]2ET0k
式中u1、u2—分别是z1、z2两个高度上的风速;
q1、q2—两个高度上的比湿; ρ—空气密度; k—卡曼常数。
六、水量平衡法
水量平衡法的基本原理是根据计算区域内水量的收人和支出的差额来推算植物蒸发蒸腾量,属于一种间接的测定方法。水量平衡方程式如下:
PeIWETDW
PePRP
式中 Pe——有效降雨量;
I——灌溉水量;
W——地下水补给量;
ET——植物蒸发蒸腾量;
D——深层渗漏量;
△W——阶段内土壤含水量的变化;
P——实际降雨量;
R——地表径流损失量;
α——降雨人渗系数,其值与一次降雨量、降雨强度、降雨历时、土壤性质、地面覆盖及地形等因素有关。
七、称重式测筒测定作物蒸发蒸腾量
蒸发蒸腾量计算公式:
G1G2GmGpGcSET12
式中 ETl-2——阶段蒸发蒸腾量,mm;
G1——时段开始时的土壤容器总重量,kg;
G2——时段末时的土壤容器总重量,kg;
Gm——时段内向土壤容器内的灌水量,kg;
Gp——时段内落入土壤容器内的降水量,kg;
Gc——时段内土壤容器中的土表及底层排水量之和,kg;
S——土壤容器内的水平截面积,m2。
八、α值法
大量灌溉试验资料表明,水面蒸发与作物需水量之间存在一定程度的相关关系。因此,可以用水面蒸发量这一参数来计算作物需水量。其计算公式如下:
ETE0
或ETE0b
式中ET—某时段作物需水量,以水层深度mm计;
E0—ET同时段的水面蒸发(E601型蒸发皿或80cm口径蒸发皿测定值),以水层深度mm计;
α—作物需水系数,即作物需水量与水面蒸发量比值,随作物生育阶段而改变,由实测资料确定,一般条件,水稻0.8~1.57,小麦0.3~0.9,棉花0.34~0.9,玉米0.33~1.0,谷糜0.5~0.72;
b—经验系数。
九、气温推算法
计算公式为
ETabT
式中ET——作物全生育期蒸发蒸腾量,mm;
T——作物全生育期日平均气温累积值,℃; b——以温度为指标的系数; a——经验常数。
十、用气温和水面蒸发推算蒸发量
其中主要有竺士林提出的方法,此法的计算公式为
ETii(ti50)E0i
式中ETi——作物第i阶段的蒸发蒸腾量,mm;
E0i ——第i 阶段水面蒸发量,mm;
ti——作物第i 阶段的日平均气温,℃;
i——阶段需水系数。
因篇幅问题不能全部显示,请点此查看更多更全内容