物理先验引导的暴雨临近预报

确定性建模与基于扩散的时空细化

高楚林, 骆一川, 冷亮, 张彤

研究背景

暴雨的混沌性和快速演变特征使得临近预报充满挑战。近年来,研究逐渐从单一的确定性或概率性方法转向结合两者的混合框架。

然而,许多现有混合模型依然是纯数据驱动的,缺乏显式的物理先验约束,导致其物理一致性和泛化能力有限,尤其在极端降水情境下表现不足。

核心创新

物理先验表征

基于偏微分方程(PDE)构建门控卷积单元,将物理定律融入深度学习架构,确保预测的物理一致性

确定性预测模块

利用多层PDEGiU解码器生成物理一致的预测结果,捕捉大尺度降水演变特征

扩散校正网络

在物理感知的隐状态条件下建模残差不确定性,细化局地结构,提升预测精度

方法框架

STEP 1 暴雨先验表征网络

采用编码器-解码器结构的自编码器,通过PDE引导的门控卷积单元(PDEGiU)编码雷达观测与降水数据为多层时空特征表征

E₁:₃,ₗ = Fᵣ(Pₗ₋ₘ₊₁:ₗ, Rₗ₋ₘ₊₁:ₗ)

STEP 2 先验引导的确定性预测

通过时空相似性聚合历史事件先验,利用多层PDEGiU解码器递归生成未来降水场,并施加时空相干性损失确保物理一致性

P̂'ₗ₊₁:ₗ₊ₙ = Fₚᵣₑ(Ẽ₁:₃,ₗ:ₗ₊ₙ₋₁)

STEP 3 先验控制的扩散校正

针对每个预报时效训练独立的扩散模型,基于先验表征和确定性预测,通过去噪扩散概率模型逐帧校正残差误差

P̂"ₗ₊₁:ₗ₊ₙ = Fₐₘ(Ẽ₁:₃,ₗ₊ₙ₋₁, P̂'ₗ₊₁:ₗ₊ₙ)

实验结果

性能指标

本文方法 CSI 0.235
vs DiffCast
+5.4%
vs PhyDNet
+10.3%
vs DGMR
+67.9%

逐时效表现

第1小时 +8.3%
第2小时 +8.6%
第3小时 +6.2%

消融实验

模型配置 Total CSI 性能变化
完整模型 0.235 基准
去除扩散模块 0.212 -9.8%
去除物理先验 0.194 -17.4%
端到端训练 0.218 -7.2%

关键发现:物理先验表征是实现物理一致且高精度降水预报的核心要素,其移除导致性能下降最为显著(-17.4%)

103.58M
模型参数量
扩散模块占101.15M
4.4s
单样本推理时间
128×128网格 / A100 GPU
0.725
结构相似性指数
SSIM / 3小时预报
引入物理先验并结合扩散残差修正能显著提升深度学习临近预报的准确性和物理一致性,特别是在极端降水情境下表现突出

数据与实验设置

数据来源

山东省及邻近海域两年历史气象观测数据(每年4-10月),包含16个气象变量

时空分辨率

时间:1小时间隔
空间:1km分辨率,128×128像素网格

预测任务

输入:6小时历史雷达反射率与降水观测
输出:未来3小时降水场预测

评估指标

关键成功指数(CSI),强降水阈值:16 mm/h

方法优势

物理一致性强

通过PDE约束确保预测符合物理定律,避免不合理的降水形态

时空结构清晰

扩散校正模块有效缓解确定性模型的过度平滑问题

极端事件表现优

在强降水场景下相比基线模型优势明显,CSI提升显著

时效稳定性好

各预报时效均保持稳定领先,误差累积得到有效控制