4. 水声信道仿真器实现

水声OFDM系统仿真中，信道仿真器决定了接收信号是否“像真实海洋”。一个可用的水声信道仿真器需要同时给出三件事：多径到达结构（幅度、时延、相位）、随时间变化的调制（多普勒与相位漂移）、以及叠加在接收端的海洋环境噪声。工程上通常有两条实现路线：以BELLHOP为代表的确定性射线追踪模型，和以统计分布描述多径幅度与时延的随机模型。本节围绕这两条路线，说明如何把它们组织成可复用的时变多径信道生成器，并给出多普勒与噪声的叠加方法。

一、两条实现路线的定位与选择

BELLHOP基于射线声学，输入是声速剖面、海深、海底/海面声学参数、收发位置与频率，输出是本征声线的到达角、时延、幅度与相位。它适合描述特定水文条件下的确定性多径结构，例如浅海负梯度下的表面反射路径与海底反射路径的分离。统计模型（如基于Rician、Rayleigh或Nakagami分布的多径幅度模型，配合指数或均匀时延分布）不依赖声场计算，参数少、生成快，适合大量蒙特卡洛仿真或链路级性能评估。实际工程中常见的做法是：用BELLHOP生成一组“基准多径剖面”，再用统计扰动在其上叠加时变性与随机性，兼顾物理合理性与仿真效率。

二、BELLHOP 多径剖面的生成与参数化

以浅海典型场景为例：海深 80 m，声速剖面为负梯度（表层 1520 m/s，底层 1490 m/s），收发水平距离 1 km，声源与接收器均位于水下 30 m，频率 12 kHz。BELLHOP 的 .env 环境文件需要给出 NSD（声速剖面点数）、NRD（接收深度数）、NR（接收距离数）以及海面、海底边界条件。运行后从 .arr 或 .ray 输出中提取每条本征声线的到达时延与复幅度，按最强径归一化，得到类似 τ=[0, 2.1, 4.7, 7.9] ms、a=[1.0, 0.62, 0.35, 0.18] 的多径表（示例值，实际由剖面与边界参数决定）。

需要注意两点工程细节。第一，BELLHOP 输出的是单频或窄带条件下的复幅度，用于宽带 OFDM 时，应把每条径视为带内平坦的复增益，并检查系统带宽内是否存在明显的频率选择性——若带宽达到数 kHz，需在多个频点重复计算或采用宽带射线模型。第二，多径表必须与 OFDM 采样率对齐：若采样率为 48 kHz，则 2.1 ms 对应约 101 个采样点，时延需量化为整数采样或采用分数延迟滤波器，否则会引入额外误差。

三、时变多径的生成：多普勒与增益扰动

静态多径表只是起点。真实水声信道中，海面起伏、平台运动与内波都会让每条径的时延和相位缓慢变化。工程上常用的简化是：对第 l 条径，保持时延 τ_l 在一个相干时间内近似不变，用复增益的时变描述其变化：

h_l(t) = a_l(t) · exp(j·2π·f_d,l·t + j·φ_l)
f_d,l  = (v/c) · f_c · cos(θ_l)      % 单程多普勒，θ_l 为到达角
a_l(t) = ρ·a_l(t-Δt) + sqrt(1-ρ²)·w(t)   % 一阶AR扰动，ρ 由相干时间决定

其中 v 为收发相对径向速度，c 为声速（约 1500 m/s），f_c 为载波频率。以 v=1 m/s、f_c=12 kHz 为例，最大单程多普勒约 8 Hz；若收发各运动，双程可达 16 Hz。对 OFDM 而言，多普勒频移会造成子载波间干扰（ICI），而各径多普勒不同则造成多普勒扩展，二者都应在仿真中体现。AR 系数 ρ 可由相干时间 T_c 估算：ρ≈exp(-Δt/T_c)，浅海典型 T_c 在几十到几百毫秒量级，需按场景设定并明确标注为假设值。

四、海洋环境噪声的叠加方法

噪声不能简单用高斯白噪声代替。水声环境噪声的功率谱级随频率变化，通常由湍流噪声、航运噪声、风关噪声和热噪声叠加而成。工程仿真中可按 Wenz 谱的简化分段模型生成噪声谱密度，再经逆傅里叶变换得到时域噪声，最后按目标信噪比定标。具体做法是：先计算信号段平均功率 P_s，再根据设定 SNR 求噪声功率 P_n=P_s/10^(SNR/10)，对生成的噪声序列做功率归一化后乘以 sqrt(P_n)。这样得到的接收信号 r(t)=Σ h_l(t)·x(t-τ_l)+n(t) 才具有可比性。

五、实现检查与结果判读

一个可信的信道仿真器应通过三项自检：其一，关闭多普勒与噪声时，接收信号应与静态卷积结果一致；其二，仅开启多普勒时，单径信道的相位应随时间线性旋转，旋转速率与设定 f_d 吻合；其三，噪声功率实测值应与目标 SNR 偏差在 0.5 dB 以内。若偏差过大，优先检查时延量化、归一化方式和功率计算窗口是否一致。完成这些检查后，再把时变多径信道接入 OFDM 收发链路，观察均衡前后的星座图与误码率变化，才能判断信道仿真器是否真正服务于系统性能评估。