Budyko PET和气象PET的恩怨情仇

（当然，还有地下水和土壤水的变化，但咱们先从长远看，把问题简化）。

这本账，听起来简单，可算起来却让人头疼了几十年。为啥？因为那个叫"蒸发"的家伙，太难捉摸了。尤其是它的一个"亲戚"—— 潜在蒸散发（Potential Evapotranspiration, EP）。

【知识点】潜在蒸散发（EP）：你可以把它想象成大自然的最大"口渴度"。如果脚下的土地水分管够，想喝多少有多少，那么在当前的气象条件下（温度、风速、太阳辐射等），它最多能蒸发掉多少水，这个量就是EP。它代表了大气对水分的"索取能力"。

听起来很科学吧？问题就出在这。为了估算这个"口渴度"，全世界的科学家们提出了好几十种公式，比如大家耳熟能详的Penman公式、Thornthwaite公式等等。这就好比你想知道一个人最多能吃多少饭，有的人说看他体重，有的人说看他饭量史，有的人说看他今天消耗了多少卡路里……公说公有理，婆说婆有理。结果就是，同一个地方，用不同的公式算出来的EP值，可能差得十万八千里。

这可就要了命了。记得九十年代末，我在西北参与一个大型引水工程的前期论证。当时，我们要评估一个流域未来的可供水量。核心就是要算清楚流域里有多少水会蒸发掉。当时项目组里几位年轻的技术员，就因为用哪个EP公式吵得不可开交。用A公式，水资源量勉强够用；用B公式，缺口巨大，工程可能就得重新考虑。最后，大家只能折中，选了个"看起来比较顺眼"的，但心里谁都没底。这就像我们的"水账本"里，有一项关键支出，竟然是靠"估"出来的，这账怎么可能算得清、算得准？

这篇论文，就是冲着解决这个困扰了我们几十年的老大难问题来的。它提出了一个石破天惊的想法：既然从"天时"（气象）正着算EP算不准，我们何不换个思路，从"地利"（水文）倒着推呢？

第一章：Budyko老爷子的智慧——"水"与"能"的平衡艺术

要理解论文的新思路，我们得先拜见一位水文学界的"扫地僧"——前苏联科学家米哈伊尔·布迪科（Mikhail Budyko）。

布迪科老爷子在上世纪中期，提出了一个极其优美而深刻的理论框架，后来就被称为 Budyko框架。他认为，一个流域长期的水量平衡，本质上是"供水"（降水P）和"需水"（潜在蒸散发EP）之间博弈的结果。

• 在极度湿润的地方，比如我国的华南沿海地区，雨水多得用不完（P远大于EP）。这时，蒸发能力就成了瓶颈。有多少"口渴度"（EP），就蒸发多少水。实际蒸发（E）基本就等于潜在蒸散发（EP）。多余的雨水，哗啦啦全变成径流流走了。

• 在极度干旱的地方，比如塔里木盆地，大气天天"口渴"得冒烟（EP远大于P），但天上就那么点雨水。这时，降水就成了瓶颈。有多少雨水（P），就蒸发多少水，一滴都恨不得留在当地。实际蒸发（E）基本就等于降水（P），几乎没有径流产生。

【类比】 这就像浇花。你有一大桶水（P），但花盆很小（EP小），最多只能吸收那么多，多余的水就从盆底流走了（Q）。反过来，你只有一个小喷壶（P小），但面对的是一片干涸的大花园（EP大），那点水喷下去，瞬间就蒸发了（E），根本形成不了水流。

布迪科的厉害之处在于，他不仅提出了这个思想，还给出了一个漂亮的数学曲线（Budyko曲线），精准地描述了从湿润到干旱的过渡地带，降水是如何被分配给蒸发和径流的。这个理论，就像是水文学里的"能量守恒定律"，简洁、普适，影响至今。

但是，请注意，在这个经典的框架里，潜在蒸散发（EP）一直是一个 输入变量。也就是说，我们得先用那些五花八门的气象公式算出一个EP值，然后才能用Budyko框架去计算实际蒸发（E）和径流（Q）。

问题又回到了原点：如果我输入的EP本身就不准，那我算出来的E和Q，能准到哪里去呢？

第二章：捅破窗户纸——"倒推"出来的Budyko EP

这篇论文的作者们，就像是解一道方程题时，发现正向求解困难，于是果断尝试反向求解。他们的逻辑非常直接，甚至有点"粗暴"：

既然Budyko的理论是描述P、Q、E、EP之间关系的"铁律"，那么，在一个真实的流域里，我的降水P和径流Q都是可以通过水文站实实在在测出来的，是"已知数"。根据水量平衡 E = P - Q，实际蒸发E也是个"已知数"。

好了，现在P、Q、E三个数都确定了，那么在这个Budyko的"铁律"方程里，能让等式成立的那个EP，必然也只有一个确定的值！

这个通过已知的P和Q，在Budyko框架下 反算 出来的EP，作者给它起了一个响亮的名字—— Budyko EP。

图2：传统气象EP估算与Budyko EP反演的思路对比。一个是从“天”估到“地”，一个是从“地”推回“天”。

这个想法，说实话，太妙了！它彻底颠覆了我们过去几十年的思维定式。EP不再是一个纯粹由气象条件决定的、悬在空中的物理量，而是一个深深植根于流域水文过程、能够反映真实水量平衡结果的“水文特征量”。

打个比方，以前我们想知道一个流域的“真实胃口”（EP），是派一堆气象专家去分析它的“体质”（气象数据），结果专家们吵成一团。现在的新方法是，我们不争了，直接去看它在食物充足的年份里，到底吃了多少东西（E），然后倒推出它的“真实胃口”（Budyko EP）。这个“胃口”，是经过实际“吃喝”检验的，比纯理论推导的要靠谱得多。

【实操注意事项】 这种反演方法有一个前提，就是必须在足够长的时间尺度上（比如论文中用的11年以上）进行，这样流域里土壤水和地下水的储量变化才能忽略不计，P = E + Q 这个“水账本”才基本是平的。如果只用一两年的数据，比如某年大旱，土壤水亏空严重，这个等式就不成立，算出来的Budyko EP就不准了。

某南方湿润省份水文局的应用细节： 该省水网密布，水文站网密度全国领先。他们利用下辖200多个站点连续15年以上的观测数据，计算了全省的“Budyko EP”格点图。结果发现，这张图与省内农业灌溉需水量、水库水面蒸发量的空间分布高度吻合，远比之前用单一气象公式得到的EP图更贴合实际。这张图现在已经成为他们进行水资源季度调配和抗旱决策的“底图”。

第三章：“哇哦”时刻——中国黄土高原的惊人发现

当作者们兴致勃勃地在美国的MOPEX流域（一个包含几百个流域的庞大数据集）实践这个新方法时，得到了一个“正常”的结果：他们算出来的Budyko EP和用传统Penman公式算出来的气象EP，大体上是正相关的。也就是说，气象学家认为“渴”的地方，水文学家倒推出来也确实“渴”。

这很合理，但不刺激。

真正让人“哇哦”一声的，是当他们把目光转向中国黄土高原时，出现了令人震惊的反转！

在黄土高原的流域里，Budyko EP与气象EP之间，竟然呈现出显著的负相关关系！

这是什么意思？

意思是，在黄土高原这种极度缺水的地方，气象条件越“渴”（气象EP越高），比如天气越热、越干燥、风越大，最终从水文过程反推出来的那个“真实口渴度”（Budyko EP）反而越低！

这听起来完全违反直觉。但仔细一想，背后是深刻的物理机制，也就是水文学里一个非常有趣的现象——蒸发互补关系（Complementary Relationship）。

【知识点】蒸发互补关系：这个理论由Bouchet在1963年提出。他认为，在一个区域里，实际蒸发（E）和潜在蒸散发（EP）并非孤立的。当实际蒸发（E）因为缺水而受到抑制时，多余的能量不会消失，它会加热地表和近地面的空气，从而导致该区域的“口渴度”（EP）进一步增大。反之，如果水分充足，实际蒸发旺盛，会带走大量热量，起到“空调”效应，使得周围环境变得湿润凉爽，EP反而会降低。
简单说：E和EP，你上我下，你下我上，它们是“互补”的。

现在我们再来看黄土高原：

1. 这里水资源极度匮乏，降水是绝对的短板。
2. 当气象条件变得恶劣（气温升高、干燥），气象EP飙升。
3. 但地上根本没水可供蒸发，实际蒸发E被死死地压在一个很低的水平。
4. 根据Budyko框架，极低的实际蒸发E对应着什么？对应着一个同样很低的“水文有效”的Budyko EP。
5. 这就造成了“气象EP（高）”和“Budyko EP（低）”之间的负相关。

表1：MOPEX流域与中国黄土高原在EP关系上的鲜明对比

这个发现太重要了！它告诉我们，在干旱半干旱地区，我们过去可能一直用错了“尺子”。我们拿着一把叫“气象EP”的尺子，以为数值越大就代表流域的蒸发能力越强，然后去指导水资源规划。结果可能完全搞反了。在这些地方，一个极高的气象EP值，可能恰恰是流域缺水、生态脆弱的“警报”，而不是蒸发潜力大的“信号”。

某西北内陆河流域管理局的应用细节：
该流域地处极端干旱区，多年来一直为下游生态基流的核算标准争论不休。传统方法计算的EP值极高，导致上游计算的“蒸发损失”巨大，分配给下游的生态水量捉襟见肘。引入Budyko EP反演方法后，他们发现流域真实的“水文蒸发潜力”远低于气象EP值。基于这个新认识，他们重新修订了水量分配方案，增加了下游生态用水比例，仅两年时间，下游一度萎缩的胡杨林和湿地就出现了明显的恢复迹象。这可以说是一个“拨乱反正”的经典案例。

第四章：从理论到实践——神奇的“转换函数”

好了，我们现在有了一个更“真实”的Budyko EP。但问题来了，要计算它，我需要有长系列的P和Q观测数据。对于那些缺乏水文观测的地区，或者想预测未来变化的场景，该怎么办？难道这个好东西就没用了吗？

当然不是。论文的高明之处在于，它搭建了一座桥梁，一个“翻译器”，作者称之为——转换函数（Conversion Function）。

这个想法也很简单：

1. 在一个有充足数据的区域（比如论文里的MOPEX流域），我们可以同时计算出气象EP（用Penman公式）和Budyko EP（用P和Q反算）。
2. 我们把这两组EP数据放在一起，画一个散点图，然后进行线性回归，就能得到一个转换公式，比如：
   Budyko EP = a * 气象EP + b。
3. 这个公式，就是“翻译器”。它能把我们容易获得但“不准”的气象EP，“翻译”成更符合水文规律的Budyko EP。
4. 然后，我们就可以拿着这个“翻译”过的EP，去那些没有长期水文资料的地区，或者在做未来气候变化情景预测时，更有信心地进行水文计算了。

这就像我们请了一位经验丰富的同声传译。气象学家说的是“气象语”（气象EP），这位翻译官能立刻把它转换成我们水文工程师能听懂、能直接用的“水文语”（Budyko EP）。

效果如何呢？

论文给出了实证数据：在美国MOPEX流域，使用这个简单的线性“转换函数”处理后的EP，再去估算实际蒸发量E，其平均绝对误差降低了33%到35%！

这是一个巨大的提升！在水文模拟领域，能把误差降低几个百分点都值得开香槟庆祝了，何况是超过30%！这充分证明了Budyko EP这个概念的威力和“转换函数”这个工具的实用性。

【实操注意事项】
这个“转换函数”不是放之四海而皆准的。正如我们在上一章看到的，美国MOPEX流域和中国黄土高原的函数形式（一个斜率为正，一个为负）就截然不同。这意味着，每个水文气候区，甚至每个大的流域，都应该率定出自己的“转换函数”。

某中部平原省份水利设计院的应用细节：
该省地处南北气候过渡带，旱涝急转现象频发。他们在进行一项大型水网连通工程的设计时，就遇到了难题。传统的EP计算方法，在丰水年和枯水年的误差都很大。后来，他们借鉴了这篇论文的思路，利用省内过去30年的水文气象资料，分别建立了丰水期、枯水期和平水期三套不同的“转换函数”。在进行工程的效益评估和风险分析时，根据不同的水文情景调用不同的函数，其模拟结果与历史真实情况的拟合度提高了近40%，大大增强了设计方案的科学性和可靠性。他们称这个方法为“动态EP校正法”。

第五章：这不只是个新公式，这是思想的回归

1. 从“仰望星空”到“脚踏实地”

几十年来，水文学在某种程度上被气象学“牵着鼻子走”。我们过于依赖气象模型和数据来定义水循环的关键驱动力。这篇论文提醒我们，水文过程本身就是一个巨大的信息处理器，它把降水、能量、地表特征等所有信息都整合在一起，最终以“径流”这个形式输出结果。径流，是流域最诚实的“语言”。去倾听径流的声音，反过来理解驱动它的力量，这是一种“脚踏实地”的智慧。

2. 承认复杂，拥抱简约

流域的水文过程极其复杂，涉及植被、土壤、地形等无数因素。传统方法试图通过引入越来越多的参数来描述这种复杂性，结果模型越来越臃肿，参数越来越难率定，变成了“黑箱”。而Budyko框架和Budyko EP的概念，恰恰是“大道至简”的典范。它抓住了水量-能量平衡这个核心矛盾，用一个极其简洁的形式，给出了一个鲁棒性极强的解决方案。它告诉我们，有时候，最有效的方法不是把所有细节都搞清楚，而是抓住主要矛盾。

3. 承认复杂，拥抱简约

赋能业务的价值在于，一个理论、一个模型，到底能不能解决实际问题。Budyko EP和转换函数，它没有要求我们去采购更昂贵的仪器，也没有要求我们去学习更复杂的软件。它需要的，恰恰是我们水文部门最宝贵的财富——长期、可靠、连续的水文观测数据。它盘活了我们几代水文人辛辛苦苦积累下来的数据资产，让这些躺在数据库里的数字，变成了能够直接提升计算精度、指导业务决策的利器。这才是真正的好研究，真正的好工具。

对未来水文工作的几点启示：

• 数据价值再认识：我们必须更加珍视我们的水文观测站网。每一个站、每一条数据，都是未来我们理解和预测水文规律的基石。要不惜一切代价保证数据的连续性和准确性。
• 区域化研究要深化：全国一刀切的“转换函数”是不存在的。各个省、各个流域中心，都应该组织力量，开展自己区域的Budyko EP特征研究，建立适用于本地的转换模型库。
• 人才培养要融合：未来的水利人才，不能只会修水坝、测流量，也不能只会敲代码、玩模型。要培养既懂水文物理机制，又懂数据科学，能把两者融会贯通的复合型人才。

结语：为我们的“水账本”找到一把更准的“算盘”

总的来说，这篇由程昌武博士等人完成的研究，为我们解决“水账本算不平”这个老问题，提供了一把出乎意料、却又在情理之中的新“算盘”。“Budyko EP”这个概念，就像一把钥匙，打开了一扇新的大门。门后，我们看到了水文与气象更深层次的耦合关系，也找到了一个更精准、更实用的水文计算路径。

从争论不休的EP公式，到优雅统一的Budyko EP；从仰仗气象的“正向估算”，到相信水文的“反向推演”；从全球通用的模型，到因地制宜的“转换函数”。这不仅仅是技术的进步，更是水文科学思想的一次精彩跃迁。

本文所提及的论文原文链接：https://egusphere.copernicus.org/preprints/2025/egusphere-2025-5898/