1. 海洋声速剖面基础
水声传播的一切计算,几乎都从一个量出发——海水中的声速 c。它决定了声波在海洋中走多远、往哪里弯、在何处汇聚或形成影区。理解声速剖面,是后续射线声学、简正波建模、声呐作用距离预报的共同起点。本节先把声速本身讲清楚,再说明它随深度变化的典型形态,以及这些形态如何改变声传播路径。
一、声速的物理定义与经验公式
声速是介质中微弱压缩扰动的传播速度。对海水这样的流体,声速由介质的绝热压缩系数与密度共同决定:
c = sqrt( K / ρ )
K:绝热体积弹性模量(Pa)
ρ:海水密度(kg/m³)
工程上很少直接测 K 和 ρ,而是用温度 T、盐度 S、压力(或深度 z)拟合的经验公式。目前水声建模中最常用的是 Mackenzie 九项公式(1981),适用范围为 T ∈ [−2, 30] ℃、S ∈ [25, 40]、z ∈ [0, 8000] m:
c(T,S,z) = 1448.96 + 4.591T − 5.304e-2·T² + 2.374e-4·T³
+ 1.340(S − 35) + 1.630e-2·z + 1.675e-7·z²
− 1.025e-2·T(S − 35) − 7.139e-13·T·z³
单位:c 为 m/s,T 为 ℃,S 为千分比,z 为 m(压力深度)
代入一组典型示例参数:T = 15 ℃、S = 35、z = 0,可得 c ≈ 1506.6 m/s;若保持 T、S 不变、z 增至 1000 m,则 c ≈ 1522.5 m/s。这组数值是公式代入的示例结果,用于说明量级,不代表某次实测。可以看到,仅压力一项在千米深度就贡献约 16 m/s 的增量,这正是深海声速在深处回升的原因。
二、三个影响因素的量级与敏感性
三个因素对声速的贡献并不对等,建模时必须分清主次:
- 温度:最敏感。约每升高 1 ℃,声速增加 3~5 m/s,且随温度升高灵敏度略降。上层海洋的季节与日变化主要通过温度体现。
- 盐度:约每增加 1 个千分比,声速增加约 1.2~1.3 m/s。除河口、融冰区和强蒸发区外,开阔大洋盐度变化小,通常不是主导项。
- 压力(深度):约每加深 100 m,声速增加约 1.6~1.7 m/s,近似线性但略有二次项。它在上层被温度压制,在深层成为主导。
三、典型声速剖面类型及其传播效应
把 c 随深度 z 的曲线画出来,就是声速剖面(SSP)。海洋中常见的形态有三类,它们对应完全不同的声线弯曲方式。
等温层(混合层):风浪搅拌使上层数十米温度、盐度近似均匀,声速随深度因压力缓慢增大,形成微弱正梯度。声线在此层内向上弯曲,贴近海面传播,利于浅层声呐与水面舰艇间的近程通信,但能量被限制在薄层内,向下泄漏少。
跃层(温跃层):温度随深度急剧下降,声速出现明显负梯度。声线强烈向下弯曲,在上方形成声影区,常规声呐在此区域探测距离骤降。跃层强度(dT/dz)越大,影区越明显,这也是夏季浅海声呐性能变差的主因。
深海声道(SOFAR 声道):上层负梯度与深层正梯度在某一深度交汇,声速取极小值,该深度称为声道轴,典型在 800~1200 m。声线被反复折射回声道轴附近,形成低损耗的远程传播通道,爆炸声源可传播数千公里。声道轴深度由温度与压力两项的竞争决定,因此随纬度、季节移动。
四、从剖面到传播:建模时的实际处理
实际建模中,SSP 通常以离散的(深度,声速)点对给出,例如每 1 m 或每 5 m 一个采样。射线追踪需要计算声速梯度 dc/dz,因此第一步是对剖面做平滑或分段线性化,避免实测噪声导致梯度剧烈跳变。一个可操作的检查流程是:
- 检查采样深度范围是否覆盖声源与接收器之间的全部深度,深海问题不能只给 0~300 m。
- 计算相邻点的 dc/dz,标出梯度符号变化处,那里就是声道轴或跃层顶底的候选位置。
- 用 Mackenzie 公式反算几个点的声速,与实测值比对,偏差超过约 1 m/s 时应检查盐度单位(千分比还是百分比)和深度是压力深度还是几何深度。
- 把剖面输入射线模型后,先看本征声线的弯曲方向是否与梯度符号一致,再判断会聚区、影区位置是否合理。
掌握声速剖面的意义在于:它把“海洋环境”翻译成了“声速随深度的分布”,而声线的每一次折射、每一个影区与会聚区,都是这条曲线梯度的直接结果。后续的射线声学与声场建模,都是在这条曲线之上展开的计算。