2. 水声信道与混响物理基础
混响抑制与回波处理的所有算法,最终都要在真实的水声信道中工作。若不了解声波在海水中如何衰减、如何被多途展宽、散射体如何分布,就无法判断一段接收信号里哪些成分是目标回波、哪些是混响、哪些只是信道畸变。本节从传播损失、多途扩展、散射体分布三个物理层面建立信道模型,再对比混响与回波在时频域的差异,最后给出浅海与深海两类典型混响环境的量级认识。
一、传播损失:能量随距离的衰减规律
传播损失 TL 描述声源级到接收点之间的能量衰减,工程上常用球面扩展加吸收的形式估算:
TL = 20·log10(r) + α·r×10⁻³ (r 单位 m,α 单位 dB/km)
例:r = 5 km,α = 0.5 dB/km(约 3 kHz 浅海)
TL ≈ 20·log10(5000) + 0.5×5 ≈ 74 + 2.5 ≈ 76.5 dB
球面扩展项 20log10(r) 反映波前面积随距离平方增长;吸收系数 α 随频率近似平方增长,是高频声呐作用距离受限的主因。浅海还存在柱面扩展(10log10(r))区间,因声能被海面与海底多次反射约束在水平方向传播。实际计算时,若已知声速剖面,可用射线追踪或简正波模型得到更准确的 TL 曲线;但作为混响建模的输入,球面扩展加吸收已能给出合理量级。需要注意:混响强度与 TL 的往返路径有关,散射体在距离 r 处产生的混响,其双程损失约为 2·TL(r),这决定了混响随距离的衰减速度远快于单程信号。
二、多途扩展:信道的时间展宽
多途扩展指同一发射脉冲经不同路径到达接收点,形成时间上展宽的到达结构。浅海中的主要路径包括:直达波、海面反射、海底反射、以及海面—海底多次反射。设海深 H、声源与接收器深度分别为 zs、zr,水平距离 r,则第 n 次海面—海底反射路径的到达时延可近似为:
τ_n ≈ (√(r² + (2nH + zs − zr)²) − √(r² + (zs − zr)²)) / c
例:H = 50 m,r = 3 km,c = 1500 m/s,zs = zr = 25 m
n=1 时延约 2×50/1500 ≈ 67 ms,n=2 约 133 ms
多途扩展的直接后果是:一个短脉冲被展宽为持续数十至数百毫秒的到达簇。对回波处理而言,这意味着目标回波不再是干净的单峰,而是被信道卷积后的波形;对混响而言,海面、海底界面的随机起伏使各条路径的时延和幅度随机变化,形成统计意义上的散射叠加。工程上常用多途扩展时延 τm 与信号带宽 B 的乘积判断信道是否“频率选择性”:当 B·τm ≫ 1 时,信道在频域出现明显梳状起伏,匹配滤波输出会出现多峰分裂,需要采用均衡或宽带波形设计。
三、散射体分布特性:混响的来源
混响由介质中大量随机散射体对入射声的散射叠加而成,按散射体位置分为海面混响、海底混响和体积混响。海面混响主要来自风浪引起的海面粗糙度与气泡层,散射强度随风速增大而增强;海底混响来自海底沉积物的粗糙界面与次表层不均匀体,其强度与底质类型(泥、砂、砾石)和入射掠射角密切相关;体积混响来自水体中的鱼群、浮游生物、气泡和温度微结构,深海还存在著名的深海散射层(DSL),在数百米深度形成强散射带。
散射体分布可用单位体积(或单位面积)散射截面描述,混响强度随时间的衰减规律近似为:
RL(t) = SL − 2·TL(c·t/2) + S_b + 10·log10(A) (面混响)
其中 S_b 为散射强度(dB),A 为照射面积,随 t 增大而增大
关键认识是:混响强度随时间的变化由“双程传播损失下降”与“照射面积增大”两个相反因素共同决定,前者使混响快速衰减,后者使其缓慢回升,二者竞争形成混响强度曲线的典型形状。这也是混响抑制算法必须处理非平稳背景的根本原因。
四、混响与回波的时频域差异
从处理角度看,混响与回波最本质的区别在于“相干性”与“随机性”。目标回波是发射波形经目标散射后的确定性副本,其时域包络窄而突出,频域能量集中在发射带宽内,相位与发射信号保持较强相关;混响是大量随机散射体回波的叠加,按中心极限定理趋于高斯分布,时域包络宽而起伏随机,频域谱形平缓且覆盖更宽,相位随机。这一差异正是匹配滤波、谱减、预白化等抑制手段的物理依据:匹配滤波对相干回波产生尖锐峰值,对随机混响只产生较低的旁瓣响应;预白化则利用混响谱的平缓性估计背景谱并加以抑制。
需要强调的是,混响并非完全“白”的:海面混响在低掠射角下具有方向性,海底混响的谱形与底质有关,体积混响在深海散射层处呈现局部增强。因此实际处理中常先做分频段、分方位的混响统计估计,再设计抑制滤波器,而不是简单假设混响为均匀白噪声。
五、典型浅海与深海混响环境
浅海(如黄海、东海大陆架,水深 20~100 m):海面与海底反射强烈,多途扩展显著,混响以海面混响和海底混响为主,且随距离呈现明显的“近程海面主导、远程海底主导”过渡。由于声能被界面多次反射约束,传播损失接近柱面扩展,混响强度高、持续时间长,对低速小目标检测极为不利。典型参数示例:风速 5~8 m/s 时海面散射强度可比平静海面高 5~10 dB;砂质海底的底散射强度在中等掠射角下约为 −25~−20 dB。
深海(如南海海盆,水深 3000~4000 m):声传播以深海声道(SOFAR)为主,声线在声道轴附近汇聚,界面反射路径少,多途扩展相对较小。混响以体积混响和深海散射层为主,散射层在 300~600 m 深度形成强散射带,夜间上浮、白天下沉,导致混响强度具有明显的昼夜变化。深海混响的绝对强度通常低于浅海,但散射层的局部增强和声道汇聚效应会使混响在特定距离上出现“亮带”,处理时需结合声速剖面和散射层深度进行建模。