现代智能雷达技术18——算法（4）

气象雷达通过双极化技术实现从"看雨强"到"识雨形"的飞跃，利用ZDR、KDP和ρHV等参量精准区分雨、雪、冰雹，并克服传统Z-R关系的局限。相控阵雷达则以电子扫描取代机械转动，将体扫时间由6分钟缩短至30秒，实现快速、自适应观测，显著提升对强对流天气的捕捉能力。尽管双极化与相控阵融合面临校准难题，但其在灾害预警中的潜力已使其成为下一代气象观测的核心基础设施。

气象雷达是雷达家族中一个极其特殊的成员。它不想看飞机，只想看水（云、雨、雪）。

它的目标不是一个点，而是充满了整个空间的无数个微小粒子。这要求我们在信号处理上，从关注"位置"转向关注"分布特性"。

双极化 (Dual-Polarization) 技术的引入，让气象雷达从"黑白电视"（只能看雨强）进化到了"彩色电视"（能区分雨、雪、雹）。

---

11.1 气象方程：从回波到雨量

气象雷达的核心任务是定量估测降水。

• 气象雷达方程 ：

  P_r = \\frac{C \\cdot P_t \\cdot G\^2 \\cdot \\theta \\cdot \\phi \\cdot \\tau}{R\^2} \\cdot \|K\|\^2 \\cdot Z

  • R2R2 衰减：注意这里是 R2R2 而不是点目标的 R4R4，因为波束越远，照射的体积越大。
  • ZZ (反射率因子) ：这是气象学的核心物理量。单位 mm6/m3mm6/m3。
    • Z=∑Di6Z=∑Di6（单位体积内所有粒子直径的 6 次方之和）。
    • 含义：水滴越大，ZZ 值指数级暴增。大雨点的回波比毛毛雨强百万倍。
  • dBZdBZ：工程上取对数，10log⁡Z10logZ。一般 10dBZ 是云，30dBZ 是小雨，50dBZ 是暴雨，60dBZ+ 可能是冰雹。

• Z-R 关系：

  • 如何从 ZZ 算出降雨率 RR (mm/h)？
  • 经验公式：Z=aRbZ=aRb（通常 a=200,b=1.6a=200,b=1.6）。
  • 痛点：这个关系很不准。同样的 ZZ 值，可能是一堆小雨滴，也可能是几个大雨滴，但降雨量完全不同。

---

在气象雷达中，双极化 (Dual-Polarization) 技术的引入是一场革命。它不仅让雷达看到了雨的"强度"，还让雷达看到了雨的"形状"。

通过同时发射和接收水平 (HH) 和垂直 (VV) 极化的电磁波，我们可以提取出一组被称为 双极化参量 的特征。这些参量就像粒子的"X光片"，揭示了云层内部的微观秘密。

---

11.2 双极化技术：粒子的"X光片"

1. 差分反射率 (ZDRZDR​) ------ 粒子的扁平度

• 定义 ：水平回波功率与垂直回波功率的比值（对数形式）。

  Z_{DR} = 10 \\log_{10} (\\frac{Z_H}{Z_V})

• 物理机制 ：
  • 小雨滴 ：表面张力主导，呈完美的球形。水平截面 = 垂直截面。
    • 特征：ZDR≈0 dBZDR≈0 dB。
  • 大雨滴 ：下落时受空气阻力挤压，变成扁平的"汉堡包"形状。水平截面 > 垂直截面。
    • 特征：ZDR>0ZDR>0（通常在 1~4 dB）。雨滴越大，ZDRZDR 越大。
  • 冰雹 ：虽然个体形状不规则，但在乱翻滚下，统计平均呈球形。
    • 特征：ZDR≈0ZDR≈0。但它的反射率 ZZ 极大（>60dBZ）。
• 应用 ：区分大雨和冰雹。如果 ZZ 很大（红区），怎么知道是暴雨还是冰雹？看 ZDRZDR。如果 ZDRZDR 大，是暴雨；如果 ZDRZDR 小，是冰雹。

2. 差分相位 (KDPKDP​) ------ 粒子的含水量

• 定义 ：水平波和垂直波穿过降水区时，相位差随距离的变化率（单位：度/公里）。

  K_{DP} = \\frac{\\Delta \\Phi_{DP}}{\\Delta R}

• 物理机制 ：
  • 电磁波在"水"中传播比在真空中慢。
  • 扁平的大雨滴在水平方向上的"水厚度"比垂直方向大，导致水平波比垂直波跑得更慢（相位滞后更多）。
  • 雨越大，水越多，这种相位差积累得越快。
• 优势 ：
  • 不受衰减影响：相位是时间量，信号衰减只会让幅度变小，不会改变相位。
  • 不受标定误差影响：雷达发射功率波动也不影响相位。
• 应用 ：暴雨定量估测。在极端暴雨中，传统的 Z-R 关系会因为信号衰减而严重低估雨量，而 KDPKDP 依然精准。

3. 相关系数 (ρHVρHV​) ------ 粒子的一致性

• 定义 ：水平回波序列和垂直回波序列在时间上的相关性（0到1之间）。

  \\rho_{HV} = \\frac{\|\\sum H \\cdot V\^\*\|}{\\sqrt{\\sum \|H\|\^2 \\sum \|V\|\^2}}

• 物理机制 ：描述一个分辨单元内，所有粒子形状的"整齐程度"。
  • 纯雨/纯雪 ：大家都是形状规则的水滴或雪花，动作整齐划一。
    • 特征：ρHV>0.97ρHV>0.97（极高）。
  • 雨夹雪/冰雹 ：有的圆，有的扁，有的乱翻滚。回波信号乱套了。
    • 特征：ρHVρHV 下降到 0.85 ~ 0.95。
  • 非气象回波 ：飞鸟、昆虫、地杂波。形状极其不规则，且方位随机。
    • 特征：ρHV<0.8ρHV<0.8（极低）。
• 应用 ：融化层识别与杂波剔除 。
  • 零度层亮带：冰晶融化成水包裹的冰，形状极不规则，ρHVρHV 会出现一个明显的低值环，这标志着 0℃ 层的高度。
  • 剔除假回波：凡是 ρHVρHV 很低的目标，直接标为"非气象"，不参与雨量计算。

总结 ：

双极化参量 (ZDR,KDP,ρHVZDR​,KDP​,ρHV​) 构成了气象雷达的多维特征空间。通过模糊逻辑 (Fuzzy Logic) 或 AI 分类算法，我们可以精准地画出空中的"相态分布图"，告诉气象员：这里在下雨，那里在下雪，中间夹着冰雹，底下还有一群鸟在飞。

---

传统的气象雷达是一个巨大的抛物面天线（大锅盖），它需要机械转动来扫描整个天空。完成一次体扫 (Volume Scan) （扫描 9-14 个不同仰角层），通常需要 5-6 分钟。

然而，对于那些极具破坏力的强对流天气（如龙卷风、下击暴流），它们的生命周期极短，几分钟内就可能生成并消亡。传统的"慢扫描"很容易漏掉关键的演变过程。

相控阵气象雷达 就是为了解决"快"的问题而生的。

---

11.3 相控阵气象雷达：追逐风暴的速度

1. 快速扫描：从分钟级到秒级

相控阵雷达利用电子扫描代替机械扫描（至少在俯仰维上），彻底改变了观测节奏。

• 一维相控阵 (1D-PAR) ：
  • 架构：水平方向机械旋转，垂直方向电子扫描。
  • 原理：天线转一圈的同时，利用电子波束在垂直方向上瞬间完成所有仰角的扫描（或者发射一个扇形宽波束覆盖所有仰角，接收端用多波束并行接收）。
  • 效果 ：原来要转 14 圈才能扫完一个体，现在转 1 圈就扫完了。
  • 速度提升 ：体扫时间从 6 分钟缩短到 30 秒。
• 全数字阵列 (Digital Array) ：
  • 利用 DBF 技术，同时形成几十个接收波束。
  • 不仅快，而且驻留时间更长，多普勒测速精度更高。

2. 自适应观测 (Adaptive Scanning)：聪明的眼睛

传统雷达很"傻"，不管晴天还是暴雨，都按固定的模式（VCP）傻扫。相控阵雷达具备了认知能力。

• 哪里有云扫哪里 ：
  • 雷达先用宽波束快速粗扫全空域。
  • 一旦发现某个区域有强回波（比如雷暴云团），立刻控制波束回头，对该区域进行高密度、高重频的精细扫描。
  • 对于晴空区域，则跳过或降低扫描频率。
• 多任务并发 ：
  • 它可以在追踪一个龙卷风的同时，分出波束去扫描机场跑道的风切变，还能顺便跟踪一架飞过的飞机。这是机械雷达绝对做不到的。

3. 校准挑战：双极化的难题

相控阵虽好，但把它和双极化结合（PAR + Dual-Pol）是世界级难题。

• 极化一致性 ：
  • 在抛物面天线中，波束中心是对称的，H/V 极化很容易匹配。
  • 在相控阵中，当波束偏离法线方向扫描时，极化状态会发生畸变 (Cross-Polarization)。H 波可能会串扰到 V 波里。
  • 后果：导致 ZDRZDR 测量误差巨大（气象要求误差 < 0.1dB）。
  • 对策 ：需要极其复杂的逐阵元校准 算法和极化补偿矩阵。这是目前各国气象雷达正在攻关的核心技术。

总结 ：

相控阵气象雷达是下一代天气预报的基础设施。它能像拍电影一样，捕捉到风暴生消的每一个瞬间，为防灾减灾争取宝贵的"提前量"。