2. 声呐信道与混响特性分析

声呐虚警控制的第一步,不是调检测门限,而是弄清楚“背景”长什么样。声呐接收到的信号从来不是目标回波加白噪声这么简单:声波在海洋信道中传播会衰减、扩展、被界面反射折射;海洋环境本身有持续的环境噪声;主动声呐发射的脉冲还会被海面、海底和水中散射体“反弹”回来,形成混响。混响是主动声呐虚警的主要来源,理解它的统计特性和时空频变化规律,才能设计出合理的检测门限和虚警控制策略。

海洋声传播损失模型

传播损失(Transmission Loss, TL)描述声源到接收点之间声强衰减的程度,通常用分贝表示。工程上最常用的是球面扩展加吸收的组合模型:

TL = 20·log10(r) + α·r        (球面扩展,r 单位 m,α 单位 dB/m)
TL = 15·log10(r) + α·r        (柱面/混合扩展,浅海远距离常用)

其中 20log10(r) 项来自球面波能量随距离平方衰减,α 是海水吸收系数,与频率强相关。以 30 kHz 为例,典型海水吸收系数约 6~8 dB/km;而 3 kHz 时仅约 0.2~0.5 dB/km。这意味着高频声呐作用距离受吸收限制明显,低频声呐则受扩展和海底反射损失主导。实际建模时还要考虑声速剖面的折射效应:负梯度(声速随深度减小)使声线向下弯曲,形成表面声道或深海声道轴,传播损失可能显著偏离简单扩展模型。对虚警控制而言,传播损失模型决定了目标回波与混响随距离的相对变化趋势,是设定距离相关门限的基础。

环境噪声与混响的统计特性

环境噪声主要来自航运、风浪、生物和热噪声,其功率谱随频率变化:低频段航运噪声占优,高频段风浪噪声和热噪声占优。噪声经窄带滤波和包络检波后,包络近似服从瑞利分布,这是平方律检波器输出在纯噪声下的经典结果。瑞利分布的特点是尾部较“薄”,超过门限的虚警概率可以解析计算:

Pfa = exp(-T²/(2σ²))

其中 T 为门限,σ 为噪声包络的尺度参数。给定 Pfa 要求,就能反解门限。但混响的统计特性完全不同:混响是大量随机散射体回波叠加的结果,当散射体数量足够多且相互独立时,中心极限定理使其趋于复高斯分布,包络同样近似瑞利;然而实际海面/海底混响存在散射体分布不均匀、界面粗糙度起伏、气泡云团等“斑点”结构,包络分布往往比瑞利更“厚尾”,用 K 分布、Weibull 分布或对数正态分布描述更贴切。厚尾意味着同样门限下虚警率更高,若仍按瑞利分布设门限,实际虚警会明显超标。

工程判断:若实测混响包络的尾部明显高于瑞利分布预测,应改用 K 分布或 Weibull 分布拟合,并据此重新计算门限,而不是简单提高几 dB 了事。

混响强度随时间与频率的变化规律

主动声呐发射脉冲后,混响强度随时间(等效于随距离)单调下降。在均匀散射体假设下,混响强度近似按 10log10(t) 或更快的速率衰减,具体取决于波束照射面积随距离的增长方式:球面扩展加脉冲长度限制时,混响强度约按 t⁻³ 到 t⁻⁴ 衰减。这形成一条“混响衰减曲线”,是距离相关门限(TVG 之外的检测门限曲线)的直接依据。

频率方面,混响强度与散射强度、吸收和波束几何都有关。高频吸收大,混响随距离衰减更快,但近程混响绝对强度也更高;低频混响衰减慢,影响距离更远。此外,海面混响在风浪大时显著增强,海底混响则与底质类型(泥、砂、岩)强相关。实际分析时,通常按距离单元统计混响包络的均值与方差,绘制“混响强度—距离”曲线,并与噪声本底比较,找出混响主导区与噪声主导区的分界点。

混响强度随距离衰减与噪声本底的关系 (示意图,参数为教学示例) 距离 r(对数刻度,m) 强度(dB,相对值) 10 100 1k 10k 100k 高频混响(30 kHz) 低频混响(3 kHz) 环境噪声本底 混响/噪声分界 混响主导区:门限需随距离调整 噪声主导区:门限趋于恒定

混响对虚警率的影响机制

混响对虚警率的影响可以从三个层面理解。第一是均值抬升:混响使背景功率高于纯噪声本底,若门限仍按噪声设定,虚警率会成倍上升。第二是统计厚尾:混响包络的尾部比瑞利分布更重,即使均值相同,超过门限的概率也更大,等效于“有效虚警率”被放大。第三是非平稳性:混响强度随距离快速变化,海面风浪、平台运动、海底底质变化都会让混响统计量在时间和空间上起伏,固定门限在某些距离段必然偏松或偏紧。

因此,虚警控制的核心思路是“让门限跟着背景走”。工程上常用恒虚警率(CFAR)处理:在待检测单元周围取参考单元估计背景功率,按预设 Pfa 乘以估计值得到门限。在混响主导区,参考单元应沿距离方向选取,并注意混响衰减带来的背景不均匀性,必要时采用分段 CFAR 或有序统计 CFAR(OS-CFAR)抑制强散射亮点对背景估计的污染。一个可操作的检查流程是:先按距离分段统计混响包络的均值与方差,拟合分布类型,计算各段在目标 Pfa 下所需门限,再与实际检测结果对比,若某段虚警明显偏高,优先检查该段是否存在未建模的强散射体或分布厚尾。

注意:混响是非平稳背景,直接用长时间平均估计背景功率会同时“抹平”混响衰减和真实目标,导致近程门限偏松、远程门限偏紧。参考单元长度应与混响相关半径匹配,通常取几个到十几个距离单元。