1. 课程导入与声学基础

水声传播研究的是声波在海洋这一复杂介质中如何产生、如何传播、如何被接收,以及环境条件如何改变这一过程。它与雷达、光学探测最大的不同在于:电磁波在海水中衰减极快,而声波是目前唯一能远距离穿透海水的信息载体。因此,从潜艇探测、鱼雷制导、水声通信,到海洋资源勘探、渔业声学评估、海底地震监测,几乎所有水下感知任务都建立在声波传播规律之上。本课程的目标,是帮助学习者建立"声源—信道—接收"的完整物理图像,并掌握用环境参数预测传播损失的基本方法。

一、从声波方程到两种描述体系

声场的基本出发点是流体介质中的小振幅声波方程。在无黏、绝热、小扰动的假设下,由连续性方程、动量方程和状态方程线性化,可得到均匀介质中的三维波动方程:

∇²p − (1/c²)·∂²p/∂t² = 0

其中 p 为声压,c 为声速。若介质不均匀,密度 ρ 与声速 c 随空间变化,方程需改写为更一般的形式,并引入声源项 f(r,t):

ρ∇·[(1/ρ)∇p] − (1/c²)·∂²p/∂t² = −f(r,t)

这一方程是所有传播模型的共同起点,但直接求解三维时变问题代价极高。工程上通常假设声源为简谐振动,取 p(r,t)=P(r)e^(−iωt),得到亥姆霍兹方程 ∇²P + k²P = 0(k=ω/c 为波数)。围绕如何求解亥姆霍兹方程,形成了两条经典路线。

波动声学保留声波的相位与干涉信息,通过简正波、抛物方程、波数积分等方法求解,能够描述多途干涉、焦散、衍射和影区等精细现象,适合低频、复杂距离相关环境。其代价是计算量大,且对频率和环境的采样要求高。射线声学则在波长远小于环境尺度(高频近似)时,把声能看作沿射线传播,用斯涅尔定律和程函方程追踪声线,计算快、物理直观,适合高频、大范围、强折射环境,但无法给出焦散区和影区边界的正确声强,也丢失干涉结构。

两者的联系在于:射线是波动解在高频极限下的渐近结果,简正波的各阶模态在传播角度上对应不同射线族。实际建模中常按频率分工——低频用波动方法,高频用射线方法,中间频段则需权衡精度与效率。

水声传播预测的总体框架 环境输入 声速剖面 c(z) 跃层 / 声道轴 / 混合层 海底与海面 底质参数 / 粗糙度 声源与接收几何 深度 / 距离 / 频率 吸收与散射 体积衰减 / 内波 不确定性来源 剖面误差 / 底质失配 声场求解 亥姆霍兹方程 ∇²P + k²P = 0 波动声学 简正波 抛物方程 波数积分 保留相位与干涉 射线声学 斯涅尔定律 程函方程 声线追踪 高频近似、计算快 高频→射线 低频→波动 按频率与精度需求选择 传播损失 TL(r,z) 多途结构 / 会聚区 / 影区 结果判读与应用 传播损失曲线 球面扩展→柱面扩展 声呐作用距离 优质因数 FOM 通信误码评估 信噪比 / 多途时延 模型校验 与实测比对、调参 反馈修正环境 迭代逼近真实声场 实测比对 → 修正声速剖面与底质参数 → 重新计算

二、环境建模为什么是传播预测的核心

上图的框架说明了一个关键事实:声场求解器本身是成熟的,真正决定预测成败的是输入的环境是否可信。海洋声速剖面受温度、盐度和压力共同控制,典型深海剖面自上而下分为混合层、主跃层和深海等温层,声速极小值处形成声道轴,声线被不断折射回轴附近,形成远距离低损耗传播。若把声道轴深度估错 100 m,会聚区位置和影区范围就会明显偏移。海底则通过声速、密度、衰减系数和分层结构影响反射与透射,泥质海底与砂质海底的反射损失可相差数分贝以上。

因此,环境建模的任务不是"填几个参数",而是把现场可测的物理量组织成模型可用的输入:用 CTD 或声速仪获取剖面并做垂向插值,用底质取样或反演确定海底参数,用海况等级估计海面粗糙度,再明确声源与接收的深度、距离和频率。建模精度应与预测目的匹配——用于声呐作用距离粗估时,平滑后的剖面已足够;用于会聚区精细预报时,则需保留跃层的细结构。

三、一个可动手的示例:等温层中的传播损失估算

为建立数量级概念,考虑一个合理假设的示例:海水声速取常数 c=1500 m/s,声源与接收均位于同一深度,忽略海底与海面反射,只计球面扩展和海水吸收。吸收系数可用经验式估算,在频率 f(kHz)下近似取 α≈0.036·f^1.5 dB/km(淡水、常温的简化估计,仅作教学示例)。传播损失定义为 TL=20·lg(r)+α·r/1000,r 以米计。

取 f=10 kHz,则 α≈0.036×31.6≈1.14 dB/km。计算几个距离点:r=1 km 时 TL≈60+1.1≈61.1 dB;r=5 km 时 TL≈74+5.7≈79.7 dB;r=10 km 时 TL≈80+11.4≈91.4 dB。可以看到,近距离以扩展损失为主,远距离吸收项逐渐占优,且频率越高这一转折越早出现。这正是高频声呐作用距离受限的物理原因。

若进一步引入深海声道,声线不再自由扩展而近似柱面扩展,TL 变为 10·lg(r)+α·r/1000,同样 10 km 处仅约 40+11.4≈51.4 dB,比球面扩展少约 40 dB。这一对比说明:环境结构对传播损失的影响远大于许多设备参数的微调,也解释了为何本课程把大量篇幅放在环境建模上。

学习建议:每学一种传播模型,都回到"它如何处理声速剖面、如何处理边界、在什么频段和距离上可信"这三个问题,避免把模型当成黑箱。后续章节将围绕声速剖面测量与建模、海底声学参数、简正波与射线模型的具体实现、以及模型与实测的比对校验逐步展开。

本课程的整体脉络是:先建立声波方程与两种描述体系的物理基础,再深入海洋环境各要素的测量与参数化,随后学习典型传播模型的计算流程与适用条件,最后通过案例完成"环境输入—声场计算—结果判读—反馈修正"的闭环训练。掌握这一闭环,才能在真实海洋条件下对传播损失作出有依据、可检验的预测。