📚 课程目录第4章降雨衰落物理机制与统计特性

4. 降雨衰落物理机制与统计特性

卫星链路工作在10 GHz以上时，降雨是造成信号衰落最主要的因素。理解雨衰，不能只停留在“下雨会衰减信号”这一层，而要从雨滴的物理形态出发，逐级建立从微观雨滴到宏观链路衰落的定量关系。本节沿着“雨滴尺寸分布 → 降雨率 → 比衰减 → 路径缩减因子 → ITU-R雨衰模型”这条主线展开，帮助读者建立完整的物理与统计认知。

一、雨滴尺寸分布：衰落的微观起点

降雨并非由单一尺寸的雨滴组成，而是包含从零点几毫米到数毫米的雨滴群体。描述这一群体的核心工具是雨滴尺寸分布（DSD，Drop Size Distribution），即单位体积内不同等效直径雨滴的数量密度 N(D)，单位为 m⁻³·mm⁻¹。工程上最常用的是伽马分布和负指数分布，其中伽马分布形式为：

N(D) = N₀ · D^μ · exp(−Λ·D)

D   ：雨滴等效直径（mm）
N₀  ：截距参数
μ   ：形状参数（反映分布弯曲程度）
Λ   ：斜率参数（mm⁻¹），与降雨率强相关

雨滴尺寸分布之所以关键，是因为衰减与雨滴直径的高次方成正比。在瑞利近似下，单个雨滴的消光截面近似正比于 D⁶；即使进入米氏散射区，大直径雨滴的贡献仍远高于小滴。因此，当降雨率增大时，不仅雨滴数量增多，大直径雨滴占比也上升，衰减会以远快于线性增长的方式加剧。这也是为什么雨衰对降雨率的敏感度极高。

二、降雨率 R：统计特性的核心参数

降雨率 R 通常以 mm/h 表示，是描述降雨强度的宏观量。它与雨滴尺寸分布通过雨滴下落末速度 v(D) 建立联系：

R = 6π×10⁻⁴ · ∫ D³ · v(D) · N(D) dD

降雨率具有强烈的随机性和地域性。ITU-R P.837 建议书给出了全球降雨率统计模型，工程上最常用的是R₀.₀₁，即一年中 0.01% 时间被超过的降雨率。例如，我国华南地区 R₀.₀₁ 可达 60～90 mm/h，而西北干旱区可能仅 20～30 mm/h。这个参数直接决定了雨衰的严重程度，是链路预算中必须查表或按当地气象数据确定的关键输入。

三、比衰减 γ_R：单位路径的衰减能力

比衰减 γ_R（dB/km）表示电磁波在雨介质中每公里传播所损失的功率，是连接降雨率与链路总衰落的桥梁。ITU-R P.838 给出的经验公式为：

γ_R = k · R^α

k、α 为与频率、极化方式相关的系数
水平极化与垂直极化取值不同
圆极化时取两者组合

以 12 GHz 水平极化为例，k ≈ 0.0188，α ≈ 1.217；当 R = 50 mm/h 时，γ_R ≈ 0.0188 × 50^1.217 ≈ 2.2 dB/km。若频率升至 20 GHz，k 和 α 均增大，同样降雨率下比衰减可翻倍以上。这说明频率越高，雨衰越严重，Ku、Ka 频段链路必须认真对待雨衰问题。

四、路径缩减因子：从“均匀雨”到“真实雨”

如果降雨沿整条斜路径均匀分布，总雨衰就是 γ_R 乘以路径长度。但真实降雨具有强烈的空间不均匀性：强降雨核心区通常只有几公里，路径两端往往只有小雨甚至无雨。若直接用 γ_R × L 计算，会严重高估衰落。为此引入路径缩减因子 r（0 < r ≤ 1），把实际非均匀路径等效为一段较短的均匀强降雨路径：

A = γ_R · L_E · r

L_E：斜路径在雨区中的几何长度
r ：路径缩减因子，ITU-R P.530 给出经验式
    r = 1 / (1 + L_G / L_R)
    L_G：水平投影长度
    L_R = 35 · exp(−0.015 · R₀.₀₁)

以 R₀.₀₁ = 50 mm/h 为例，L_R = 35·e^(−0.75) ≈ 16.5 km。若地球站仰角 30°、雨区高度 4 km，则斜路径长度约 8 km，水平投影 L_G ≈ 6.9 km，r ≈ 1/(1+6.9/16.5) ≈ 0.70。可见缩减因子把 8 km 的几何路径折算为约 5.6 km 的等效强降雨路径，更符合实际观测。

五、ITU-R 雨衰模型：工程计算的完整链条

ITU-R P.618 建议书整合了上述各环节，形成标准化的长期雨衰预测流程。其核心步骤可概括为：

由当地降雨率统计确定 R₀.₀₁；
由频率、极化查 P.838 得到 k、α，计算 γ_R = k·R₀.₀₁^α；
由仰角 θ 和雨区高度 h_R 计算斜路径长度 L_S = (h_R − h_s)/sinθ，再得水平投影 L_G = L_S·cosθ；
计算缩减因子 r₀.₀₁；
得到 0.01% 时间 exceeded 的雨衰 A₀.₀₁ = γ_R·L_S·r₀.₀₁；
对其他时间百分比 p，用 A_p = A₀.₀₁·(p/0.01)^(−(0.655+0.033·ln p − 0.045·ln A₀.₀₁)) 外推。

需要特别说明的是，上述公式中的仰角 θ 以地平线为 0°、天顶为 90° 定义，全文及图示均采用这一约定。雨区高度 h_R 在温带地区常取 3～5 km，具体可按 ITU-R P.839 的全球雨顶高度图查取。

工程提示：雨衰计算中，R₀.₀₁ 的取值对结果影响最大。若当地缺乏长期雨量计数据，应优先采用 ITU-R P.837 数字地图，而不是简单套用邻近城市数值。此外，圆极化链路的 k、α 需按 P.838 的极化组合公式计算，不能直接取水平或垂直极化值。

六、统计特性与可用度评估的衔接

雨衰的统计特性体现在“时间百分比”这一维度上。A₀.₀₁ 表示一年中 0.01% 时间（约 52.6 分钟）雨衰不超过该值，对应链路可用度 99.99%。若系统要求 99.9% 可用度，则需关注 A₀.₁；若要求 99.99%，则必须满足 A₀.₀₁ 对应的衰落余量。链路预算中，雨衰余量通常按可用度目标反查 A_p，再与自由空间损耗、大气吸收、天线增益等合并计算。