4. 海面粗糙度对镜面反射相干分量的影响计算

问题拆解与计算链

海面镜面反射的相干分量取决于海面相对于声波波长的“粗糙程度”。工程上常用瑞利粗糙度参数判断：当海面均方根波高 σ_h 增大，反射波相位随机化加剧，镜面反射系数被相干损失因子 ρ 衰减。本例以 Pierson-Moskowitz（P-M）充分成长海浪谱为输入，从风速反演有效波高与均方根波高，再结合入射掠射角计算 ρ，最后对比不同风速下海面路径的等效反射幅度。

P-M 谱与均方根波高

P-M 谱以风速 U₁₀（海面以上 10 m 处风速）为唯一参数，其有效波高与风速的经验关系为：

H_s ≈ 0.21 · U10² / g        （g = 9.81 m/s²）
σ_h = H_s / 4                 （窄带高斯海面近似，H_s ≈ 4σ_h）

该式给出的是充分成长海况下的统计量，属于工程近似，实际海况受风区、风时限制会偏离。本例取 U₁₀ = 4、8、12 m/s 三档，分别对应轻浪、中浪、大浪量级。

瑞利粗糙度判据与相干损失因子

对掠射角 θ（相对海面）的入射声波，海面两点间相位差的标准差由瑞利参数描述：

Γ = 2 · k · σ_h · sinθ        k = 2π/λ
ρ = exp( -Γ² / 2 )            （高斯相位起伏下的相干分量衰减）

当 Γ ≪ 1 时 ρ→1，海面近似镜面；Γ 增大到约 1 时 ρ 降至 0.6 左右，镜面分量显著削弱。注意 θ 为掠射角，θ 越小（越接近水平入射）sinθ 越小，同样海况下相干性反而更好——这是海面多径建模中容易被忽略的方向性。

计算示例（f = 1 kHz，λ = 1.5 m，θ = 10°）

表中 Γ 按 k = 4.19 rad/m、sin10° ≈ 0.174 计算。可见风速从 4 m/s 增至 12 m/s，海面镜面反射路径幅度衰减超过 20 dB，多径模型中若仍按理想镜面反射系数处理，会严重高估海面反射路径强度。

结果判读与建模建议

• 轻风（U≤4 m/s）：ρ 接近 1，海面可近似为镜面，两径/三径模型直接使用理想反射系数即可，误差在 0.3 dB 量级。
• 中风（U≈8 m/s）：ρ 约 0.6，海面路径幅度下降 4 dB 以上，若忽略粗糙度会高估海面反射抽头增益，导致多径干涉深度预测偏大。
• 大风（U≥12 m/s）：ρ 趋近 0，镜面反射路径基本消失，能量转入漫散射背景，此时应改用统计散射模型而非确定性镜面抽头。
• 掠射角敏感性：同一风速下 θ 从 10° 增至 30°，sinθ 增大近 3 倍，Γ 随之增大，ρ 显著下降。建模时须与实际传播几何一致地取掠射角，不能统一按垂直入射估算。