不是订正的订正：使用 Python 计算混合层 CAPE (MLCAPE)

用户11172986

发布于 2026-04-24 19:16:13

860

文章被收录于专栏：气python风雨气python风雨

【进阶篇】拒绝“虚高”：使用 Python 计算混合层 CAPE (MLCAPE)

在上一篇文章中，我们尝试了“手动修改地面温度”来预估下午的对流潜势。虽然这种方法直观，但它依赖于主观的温度预报（猜测下午能升到多少度）。

如果地面因为局地因素（比如一片水泥地晒得很热）出现了一个很薄的超绝热层，计算出的 SBCAPE（基于地面的 CAPE）往往会虚高。气块一旦上升几百米脱离这个薄层，浮力可能就没那么大了。

为了消除这种“地面噪点”，更科学的方法是计算 **混合层 CAPE (Mixed Layer CAPE, MLCAPE)**。

1. 什么是混合层订正？

它的核心逻辑是：不要只看脚底下的气块，而要看整个“大气边界层”的平均状态。

操作方法：通常将地面以上 50hPa 或 100hPa 厚度内的空气进行平均（温度、露点混合），用这个混合后的平均状态作为气块的起点。

物理意义：这代表了边界层内充分混合后的空气团，它更能代表真实参与对流的大尺度气团性质。

对比：

**SBCAPE (Surface Based)**：激进派，往往最大，容易高估。
**MLCAPE (Mixed Layer)**：稳健派，通常比 SBCAPE 小，但如果它依然很大，说明爆发对流的把握非常大。

2. Python 代码实现

MetPy 极其贴心地内置了 mixed_layer_cape_cin ，让我们无需手写平均算法。

import matplotlib.pyplot as plt
from netCDF4 import Dataset
from wrf import getvar, ll_to_xy
import metpy.calc as mpcalc
from metpy.plots import SkewT
from metpy.units import units

# ================= 1. 数据准备 (同前文) =================
# 假设 ncfile 已经读取
ncfile = Dataset("/home/mw/input/typhoon9537/wrfout_d01_2019-08-08_19_00_00") 
lat0, lon0 = 25.0, 120.4
x, y = ll_to_xy(ncfile, lat0, lon0)

# 读取数据并赋予单位 (必须!)
p  = getvar(ncfile, "pressure")[:, y, x].values * units.hPa
tc = getvar(ncfile, "tc")[:, y, x].values       * units.degC
td = getvar(ncfile, "td")[:, y, x].values       * units.degC

# ================= 2. 计算两种 CAPE 进行对比 =================

# --- A. 计算原始的基于地面的 CAPE (SBCAPE) ---
prof_sb = mpcalc.parcel_profile(p, tc[0], td[0])
sb_cape, sb_cin = mpcalc.surface_based_cape_cin(p, tc, td)

# --- B. 计算混合层 CAPE (MLCAPE) ---
ml_cape, ml_cin = mpcalc.mixed_layer_cape_cin(p, tc, td, depth=50 * units.hPa)

# 分两步获取混合层气块路径 (用于画图)
# 1. 获取混合后的气块起始属性 (气压, 温度, 露点)
mp_mixed, mt_mixed, mtd_mixed = mpcalc.mixed_parcel(p, tc, td, depth=50 * units.hPa)

# 2. 将混合后的 T 和 Td 作为起点，计算上升路径
prof_ml = mpcalc.parcel_profile(p, mt_mixed, mtd_mixed)

print(f"激进派 (SBCAPE): {sb_cape:.2f}")
print(f"稳健派 (MLCAPE): {ml_cape:.2f}")

# ================= 3. 绘图分析 =================
fig = plt.figure(figsize=(9, 9))
skew = SkewT(fig, rotation=45)

# 1. 绘制环境场
skew.plot(p, tc, 'r', lw=2, label='Env Temp')
skew.plot(p, td, 'g', lw=2, label='Env Dewpoint')

# 2. 绘制 SBCAPE 路径 (虚线，作为参考)
skew.plot(p, prof_sb, 'gray', lw=1, ls='-', alpha=0.5, label='SB Parcel (Surface)')

# 3. 绘制 MLCAPE 路径 (实线，作为重点)
# 注意：混合层气块的起点通常不在地面，而在地面以上一点点，或者温度比地面低
skew.plot(p, prof_ml, 'blue', lw=2, ls='-', label='ML Parcel (Mixed 50hPa)')

# 4. 填充 MLCAPE 面积
# 只有当混合层气块比环境温度高时才填充
skew.shade_cape(p, tc, prof_ml, alpha=0.3, color='royalblue') 

# 5. 辅助线
skew.plot_dry_adiabats(alpha=0.25)
skew.plot_moist_adiabats(alpha=0.25)
skew.plot_mixing_lines(alpha=0.25)

# 设置范围
skew.ax.set_ylim(1050, 100)
skew.ax.set_xlim(-30, 40)

# 标题与统计
title_info = (f"Comparison @ ({lat0}, {lon0})\n"
              f"SBCAPE={sb_cape.m:.0f} (Unstable Surface)\n"
              f"MLCAPE={ml_cape.m:.0f} (Mixed Layer)")
plt.title(title_info, loc='left')
plt.legend(loc='upper right')

plt.show()

输出示例：

激进派 (SBCAPE): 1691.50 joule / kilogram
稳健派 (MLCAPE): 1182.78 joule / kilogram

image

3. 结果解读

运行这段代码后，你通常会看到以下现象：

SBCAPE 数值较大：灰色的虚线（SBCAPE路径）通常在最右侧，因为它是直接把地面那个最热的气块提上去的。
MLCAPE 数值较小：蓝色的实线（MLCAPE路径）通常在灰色虚线左侧。这是因为把地面以上 50hPa 的冷空气和地面热空气一平均，气块的初始温度就降低了。

真实性判断：

如果 MLCAPE 依然很大（例如 > 1500 J/kg），说明即使考虑了边界层的湍流混合，环境依然极度不稳定，这种强对流预报的可信度极高。
如果 SBCAPE 很大但 MLCAPE 很小，说明不稳定能量仅限于极浅的地面层，一旦有扰动混合一下，能量就消散了，这种情况下预报强对流要非常谨慎（容易空报）。