简介:本资源为彭曼公式完整推导过程的docx文档,面向水文学、气象学及农业灌溉领域的学习者与研究者,帮助读者系统理解开放水面蒸发量估算公式的数学来源与物理意义。文档从能量守恒出发,逐步推导短波净辐射、长波净辐射、对流热交换与潜热蒸发之间的平衡关系,并引入饱和水蒸气压随温度变化的斜率Δ,结合Magnus经验公式与质量传输方法,最终得到可直接用气象参数估算蒸发量的表达式。压缩包内仅含1个docx文件,约24KB,内容为完整的公式推导与参数说明,适合作为课堂笔记或科研参考。目前已有3740人学习下载,读者可借此掌握彭曼公式的推导逻辑、关键参数含义及其在作物需水、水库蒸发损失评估等场景中的应用思路。
1. 彭曼公式到底是什么:一个“算蒸发”的万能公式
你有没有想过,太阳晒着一盆水,风一吹,水就慢慢变少了。这少掉的水到底去哪儿了?答案就是蒸发。可如果我问你:今天这盆水蒸发了多少毫米?你可能就答不上来了。彭曼公式,就是用来回答这个问题的。
它是一套基于能量守恒的数学推导,专门估算开放水面(比如水库、湖泊、农田里的水盆)每天蒸发掉多少水。通俗地说,它把“太阳晒进来的热量”和“风带走的水汽”这两件事放在一个天平上,让两边平衡,从而算出蒸发量。搞懂它,你就能理解农业灌溉为什么要看天气、水库为什么夏天水位下降得特别快。它不只是一个公式,更是一套把气象数据翻译成“水账本”的思维方式。
2. 为什么需要彭曼公式:从一盆水到一片农田的蒸发难题
2.1 直接测量蒸发有多难
要算蒸发,最直接的办法是拿一个标准蒸发皿,装上水,每天量水位下降了多少。这个方法确实准,但它有个致命问题:只能代表蒸发皿那一小块地方。农田里的作物、水库宽阔的水面、森林的树冠,它们的蒸发条件跟一个铁盆完全不一样。你不可能在每个田块、每片湖面都摆一个蒸发皿,更不可能在几十年前的气象记录里找到蒸发皿数据。
所以,科学家需要一种“间接算法”:用常规气象站就能测到的数据——温度、风速、日照、湿度——来推算蒸发量。彭曼公式就是干这个的。它把蒸发过程拆成两个驱动力:一个是能量,也就是太阳辐射提供的热量;另一个是空气动力学,也就是风和水汽压差把水汽带走的能力。这两股力量合在一起,决定了水蒸发的速度。
2.2 能量平衡:蒸发本质是一场“热量搬家”
彭曼公式的起点是一个能量守恒方程。想象一个开放水面,太阳的短波辐射(K)和大气长波辐射(L)不断给水体输入热量。这些热量有三个去处:一部分用来把水变成水蒸气,也就是蒸发消耗的潜热(E·lv·rw);一部分通过湍流对流交换给大气(H);还有一小部分储存在水体里或者跟周围环境交换,但彭曼在推导时直接忽略了这部分。
于是能量平衡写成:
K + L - H = E·lv·rw
这个式子的意思是:净辐射收入减去对流热损失,剩下的能量全部用来蒸发。这里E是蒸发量(厘米/天),lv是汽化潜热(卡/克),rw是水的密度(克/立方厘米)。注意单位,这是早期水文气象学常用的单位制,现在国际单位制下会换成瓦每平方米和毫米每天,但物理意义完全一样。
2.3 为什么不能直接解:水面温度是个“黑匣子”
能量平衡方程里,对流热交换H的表达式是:
H = KH · va · (Ts - Ta)
其中KH是对流热交换系数,va是风速,Ts是水面温度,Ta是大气温度。问题来了:Ts通常没有观测数据。气象站测的是百叶箱里的气温,不是水面温度。水面温度受辐射、蒸发、水体混合等影响,跟气温差别很大。没有Ts,H就算不出来,整个方程就卡住了。
彭曼的聪明之处在于,他绕开了Ts。他引入了一个关键参数Δ——饱和水汽压随温度变化的斜率。这个Δ把“温度差”和“水汽压差”联系了起来。因为水面饱和水汽压es和大气饱和水汽压ea的差值,除以温度差Ts-Ta,近似等于Δ。于是:
Ts - Ta ≈ (es - ea) / Δ
把这个关系代入H的表达式,Ts就被消掉了,H变成了用es、ea和Δ表达的式子。这一步是整个推导的转折点,也是彭曼公式最核心的技巧。
3. 推导过程拆解:从能量平衡到彭曼公式的几步关键跳跃
3.1 引入质量传输:蒸发量本身也能表达成水汽压差
消掉Ts之后,H的表达式里出现了es和ea。但这两个量还是未知的,因为e*s是水面饱和水汽压,需要知道水面温度才能算。彭曼再次用了一个近似:根据质量传输理论,蒸发量E本身可以写成:
E = KE · va · (e*s - ea)
这里KE是质量传输系数,ea是大气实际水汽压。这个式子的物理含义很直白:风越大、水面饱和水汽压与大气水汽压差越大,蒸发越快。于是(e*s - ea)可以用E除以KE·va来表示。把这个关系代回H的表达式,H就变成了包含E的式子。
3.2 代回能量平衡:得到一个关于E的方程
现在能量平衡方程里,H被替换成了包含E的表达式。整理之后,方程左边是E乘以一个系数,右边是净辐射(K+L)加上一个与(e*a - ea)相关的项。经过移项、提取公因子、乘以Δ等代数操作,最终得到:
E = [Δ(K+L) + KH·va·(e*a - ea)] / [lv·rw·(Δ + γ)]
这就是彭曼公式的经典形式。其中γ是湿度计常数,定义为γ = ca·P / (0.622·lv),ca是空气热容量,P是大气压。γ的物理意义是:它衡量了空气“容纳水汽”和“储存热量”之间的比例关系。
3.3 从水面到作物:Monteith的修正
彭曼公式原本是算开放水面蒸发的。但农田里的蒸发不只是水面蒸发,还有作物叶片的蒸腾。叶片气孔对水汽有阻力,这个阻力跟开放水面完全不同。Monteith在彭曼的基础上,引入了冠层导度Ccan和大气导度Cat,把水汽传输分成两步:先从气孔到冠层表面,再从冠层表面到大气。总导度Ct = Cat / (1 + Cat/Ccan)。
修正后的彭曼-蒙特斯公式变成:
ET = [Δ(K+L) + ca·ra·Cat·(1-Wa)·e*a] / [lv·rw·(Δ + γ(1 + Cat/Ccan))]
其中Wa是相对湿度。这个公式把作物本身的生理特性(气孔开闭)也纳入了计算,使得蒸发估算从“水面”扩展到了“农田”。这也是为什么现在农业气象里用的基本都是彭曼-蒙特斯公式,而不是原始彭曼公式。
3.4 一个简化计算示例
假设某天测得:净辐射(K+L) = 500 cal/cm²/day,气温Ta = 25°C,相对湿度Wa = 60%,风速va = 2 m/s。先算饱和水汽压e*a。用Magnus公式:
import math def sat_vapor_pressure(T_celsius): # Magnus公式,返回单位hPa return 6.112 * math.exp(17.67 * T_celsius / (T_celsius + 243.5)) Ta = 25.0 e_a_star = sat_vapor_pressure(Ta) # 饱和水汽压 Wa = 0.60 e_a = Wa * e_a_star # 实际水汽压 delta = 2508.3 * math.exp(17.3 * Ta) / ((Ta + 237.3) ** 2) # 斜率Δ print(f"e*a = {e_a_star:.2f} hPa, ea = {e_a:.2f} hPa, Δ = {delta:.2f} hPa/°C")这段代码先算饱和水汽压,再算实际水汽压和斜率Δ。Δ的公式来自原文中的经验式,注意温度单位是摄氏度。得到Δ之后,再结合风速、净辐射和湿度计常数,就能算出E。实际业务中,这些计算会被封装成函数,输入日常气象数据即可输出参考蒸发量。
4. 易混概念辨析:彭曼公式、彭曼-蒙特斯与蒸发皿系数
4.1 彭曼公式和彭曼-蒙特斯公式是一回事吗
易混点:很多人把两者混着叫,甚至以为蒙特斯只是改了个名字。
为什么混:两者形式相似,都包含Δ、γ、净辐射和空气动力学项,而且蒙特斯确实是在彭曼基础上改的。
怎么区分:原始彭曼公式算的是开放水面蒸发,没有考虑作物气孔阻力。彭曼-蒙特斯公式引入了冠层导度和大气导度,算的是作物蒸腾加土壤蒸发,也就是参考作物蒸发量。如果你要算水库蒸发,用原始彭曼;如果你要算农田灌溉需水,用彭曼-蒙特斯。现在FAO推荐的参考作物蒸发量公式就是彭曼-蒙特斯。
4.2 参考蒸发量和实际蒸发量有什么区别
易混点:看到“蒸发量”三个字就以为是一回事。
为什么混:都叫蒸发量,单位都是毫米每天。
怎么区分:参考蒸发量是在标准条件下(假设作物高度12厘米、表面阻力70秒/米、反射率0.23)算出来的,代表大气蒸发需求,跟实际种什么作物无关。实际蒸发量还要乘以作物系数Kc,Kc随作物种类和生长阶段变化。所以参考蒸发量是“天气有多想蒸发”,实际蒸发量是“地里到底蒸发了多少”。
4.3 湿度计常数γ和斜率Δ是一对容易搞反的参数
易混点:两个都是希腊字母,都跟温度、气压有关,经常被弄混。
为什么混:公式里Δ和γ总是成对出现,而且都影响分母。
怎么区分:Δ是饱和水汽压曲线的斜率,单位是hPa/°C,只跟温度有关,温度越高Δ越大。γ是湿度计常数,单位也是hPa/°C,跟气压和温度有关,但主要反映空气的物理属性。简单记:Δ描述“水汽压随温度变多快”,γ描述“空气容纳水汽的能力”。在彭曼公式里,Δ在分子,γ在分母,两者共同决定辐射项和空气动力学项的权重。
4.4 净辐射K+L到底包括什么
易混点:以为K就是太阳总辐射,L就是大气长波辐射。
为什么混:原文里K叫“短波净辐射输入”,L叫“长波净辐射输入”,但“净”字容易被忽略。
怎么区分:K是短波净辐射,等于入射太阳短波辐射减去反射部分。L是长波净辐射,等于大气向下长波辐射减去水面向上长波辐射。两者都是“净”值,不是单向通量。如果只拿太阳总辐射代入,算出来的蒸发量会偏大,因为没扣掉反射和长波损失。
4.5 为什么忽略水体储热项
易混点:能量平衡方程里没有储热项,有人以为推导错了。
为什么混:真实水体确实会吸热和放热,尤其深水湖。
怎么区分:彭曼推导针对的是“开放水面”且时间尺度为天。对于浅水或日尺度,水体储热变化相对净辐射很小,忽略它是合理简化。但对于深水库或季节尺度,储热项不能忽略,否则夏季蒸发会被高估。所以彭曼公式用的时候要看清适用条件,不是万能钥匙。
5. 什么时候用、什么时候别用:一份蒸发估算的边界清单
彭曼公式和它的修正版,在农业灌溉、水库调度、干旱监测里用得非常多。比如某高校实验室做区域水分循环模拟,输入逐日气温、风速、湿度和日照时数,用彭曼-蒙特斯算出参考蒸发量,再结合土壤水分模型判断旱情。又比如某公司做灌溉决策系统,把气象站数据接入公式,每天给出“今天该浇多少水”的建议。这些场景的共同点是:有常规气象数据,需要连续、大范围的蒸发估算,而且能接受一定误差。
但有些情况它就不灵了。第一,极端干旱地区,空气动力学项占主导,彭曼公式对风速和湿度误差非常敏感,算出来可能偏差百分之几十。第二,冠层结构复杂的森林或高秆作物,冠层导度不是常数,彭曼-蒙特斯需要参数本地化,直接套用FAO推荐值会翻车。第三,时间尺度小于小时,比如研究午后雷阵雨前的蒸发脉冲,日均公式会平滑掉关键波动。第四,水面温度与气温差异极大的场景,比如温泉湖或工业冷却池,Δ近似不再成立。
注意:彭曼公式算的是“潜在蒸发”或“参考蒸发”,不是实际蒸发。实际蒸发还受土壤水分限制。土壤干了,蒸发就断了,公式再准也没用。
我自己的习惯是:拿到一个蒸发估算需求,先问三个问题——时间尺度是小时还是天?下垫面是水面还是作物?有没有土壤水分数据?三个答案决定用原始彭曼、彭曼-蒙特斯还是更复杂的SVAT模型。从那以后我每次用彭曼公式之前,都强制走一遍这个清单,省掉了很多返工。希望帮到你。
本文还有配套的精品资源,点击获取