简介TVAPM.zip围绕水声学中的声波传播与仿真建模面向水声研究人员、工程师及具备MATLAB基础的学习者。压缩包共23个文件包含18个.m脚本、4个.dat数据文件和1个pdf说明文档整体仅406KB。MATLAB代码覆盖海面/海底边界处理、声速剖面、信道滤波、波导仿真及接收/发射位置评估等模块dat文件提供声速、海底地形、信号等样例数据pdf可作为模型使用手册。该资源将MATLAB数值计算与水声传播物理模型相结合尤其关注与Bellhop等模型的衔接应用。目前已有319人浏览学习。通过研读源码与运行示例可掌握水下声信道仿真流程、理解声速与边界条件对传播的影响并在此基础上扩展或改进自己的水声模型适用于科研验证与工程预研。1. 直接在便携设备上跑通的 TVAPM 水声模型先要搞清它算什么野外做海试时最怕的不是声速剖面测不全而是地形变化的地方没有合适的快速预报工具。TVAPM 就是那种能应对水平变化环境的时变水声传播模型它从抛物方程出发把声速场和海底地形随时间的变化也放进递推步进里。pocketi1i 是这套源码里针对小算力设备的裁剪分支解压 tvapm.zip 后可以用一条编译宏把整个水声模型压进一台便携单板机上。这篇文章写给准备在真实海深数据上用 TVAPM 算传播损失的人不聊大而全的声学理论只讲我怎么选择模型、怎么把参数设对、怎么让输出能用来做分析。2. 水声模型选型TVAPM 为什么适合水平变化环境算传播损失之前要回答两个问题频率多高海底地形在中距离上有没有明显变化。如果海底一直是平的声速剖面也稳定简正波模型往往最快一旦出现斜坡、海山或锋面能量在距离上就不是单一模态能描述的。TVAPM 的处理方式和经典抛物方程一样把声场沿距离一步步往前推每一步都能看到当前地形和水文条件这就把水平变化变成逐段近似问题。2.1 三种经典水声模型的分工TVAPM 处在哪个位置水声行业里最常用的是射线、简正波和抛物方程三类模型各有各的适用区间。射线模型把声线当成高速路适合高频和短距离但在地形剧烈变化的区域声线要不停重算多途到达结构也很难收敛。简正波模型在低频低地形变化时非常准但它的基本假设是环境沿距离不变对斜坡只能分段拼接。抛物方程模型则不同它把距离方向当作递推轴每一步都更新环境参数天然能处理水平变化。模型类型代表性方法典型频率对水平变化的处理能力射线模型声线追踪通常 1 kHz 以上逐段追踪拐点和焦散处容易失真简正波模型模态叠加低于 500 Hz弱水平变化斜坡需分段抛物方程模型逐步递推0.1–10 kHz天然支持水平变化TVAPM 属于这一类我会在斜坡、内波和锋面环境下优先选 TVAPM而不是普通射线模型因为 TVAPM 能输出带有慢时间维的声压场而射线模型每次只能给出一个频点下的本征声线。2.2 TVAPM 的时变项改进的是 Helmholtz 方程中的哪一项抛物方程最基本的改写思路是把亥姆霍兹方程中的声压拆成缓变包络和快速振荡项的乘积。假设沿距离方向的主传播项是 e^{ik0r}忽略包络对距离的二阶导数椭圆型方程就简化成抛物型初值问题。TVAPM 的改动是把波数 k 从只依赖空间改成依赖空间和慢时间 t也就是 k(r,z,t)用来描述内波、潮汐或声速剖面缓慢漂移的阶段。实际计算时不会真的把时间变量也做一次步进。常见做法是将整个观测时间切成窗口每个窗口内冻结声速剖面和海底地形跑一次水平步进窗口之间用时间轴上的三阶插值来衔接声压包络。这样最终得到的是多个时刻的传播损失切片而不是一个长时间平均值。pitfall 是窗口长度不能取得比内波周期还小否则插值出来的相位会出现不连续跳跃。2.3 pocketi1i 裁剪版本保留了哪些计算路径pocketi1i 是 tvapm.zip 源码里针对嵌入式设备的一个编译目标与其说它是一个新模型不如说是一套削减了资源占用的构建路径。解压后先看源码里哪些地方引用了这个宏直接 grep 就能搜到grep -rn POCKETI1I src/ config.mk上面的命令会列出宏在源码里的定义和所有引用点。常见的设计是当 POCKETI1I 被定义时编译器把默认的 double 浮点改成 float把内部 FFT 缓冲上限压到 4096并且关闭逐时刻的历史场输出。这样在 ARM 单板机上单个频点的运行时间能下降一个级别代价是远距离传播损失精度从 0.01 dB 量级放宽到 0.1 dB 量级。对海上现场判断够用但不适合拿来做高分辨反演。3. 在 pocketi1i 上跑通 TVAPM 的最小环境与命令很多人拿到 zip 包直接就开始改声速剖面这往往浪费一整天。TVAPM 这类抛物线方程代码对环境输入非常敏感先应该用一个简单斜坡场景把流程跑通再换真实数据。下面这组命令和文件就是我在调试时反复使用的最小集。3.1 解压 tvapm.zip 后先确认这些文件把压缩包解开之后先不要进 src 目录先看有没有输入文件。常见流程是unzip tvapm.zip -d tvapm cd tvapm ls -la解压后的根目录里通常会有 src、inp 和若干扩展名为.env或.prm的示例文件。.env文件是环境配置也就是模型的输入面。很多移植版本把示例放在 inp 目录里直接用它默认的slope.env跑一次能验证编译是否成功。如果手里没有现成示例就按照下一节的格式自己写一个最小环境文件。不要一上来就改 Makefile先确认源码在 x86 机器上能编译再交叉编译到 pocketi1i 目标板这种顺序能省掉大半环境问题。3.2 最小环境文件让模型跑一个 2 km 斜坡斜坡场景是最合适的最小测试因为它能暴露出网格和海底参数问题。下面是一份我按照 TVAPM 类 RAM 风格输入格式写的最小配置文件title slope_2km freq 800.0 source_depth 20.0 recv_max_depth 100.0 range_max 2000.0 azimuth_angle 0.0 dz 0.3 dr 10.0 nsources 1 bottom_model sediment sediment_cp 1700.0 sediment_cs 300.0 sediment_rho 1.8 sediment_alpha 1.2 bathy_start 100.0 bathy_slope -0.02 output_file slope_2km.tl这里的 freq 是中心频率 800 Hz海水声速按 1500 m/s 估算时波长约 1.875 mdz 取 0.3 m 保证每波长约 6 个格点。dr 取 10 m斜坡斜率 -0.02每个水平步内水深只变化 0.2 m不会在一步内出现地形跳变。sediment_cp 是海底纵波速度sediment_cs 是横波速度如果模型允许浅海区不要把 cs 设为 0否则横波能量无法从海床带走远距离传播损失会偏大。3.3 编译命令和后台运行技巧在 PC 上先不做裁剪直接用默认方式编译cd tvapm make clean make -j4 ./bin/tvapm inp/slope_2km.env run.log 21运行日志会记录每个水平步的完成情况。如果 run.log 里出现能量非单调上升通常是海底参数或网格问题不是程序 bug。确认 PC 版正常后再编译给 pocketi1i 用的版本make clean CFLAGS-DPOCKETI1I1 -O2 -ffast-math make这里的关键参数是-DPOCKETI1I1它触发源码里的单精度路径-ffast-math能让速度更快但会把浮点运算顺序重排。我一般只在需要快速看趋势时加这个开关做精确对比时会去掉因为差 0.1 dB 的传播损失在反演里可能被当成真实信号。3.4 用 Python 把二进制输出读成传播损失沿深度变化TVAPM 的输出大多数是按水平步长和深度网格排列的复声压通常以 float32 写入。下面这个函数是我在现场电脑上用来快速检查结果的做法import numpy as np def read_tvapm_tl(fname, nz, nr): 读取 TVAPM 输出的复声压返回 TL 矩阵。 nz: 垂直网格点数nr: 水平步数 raw np.fromfile(fname, dtypenp.float32) comp raw.view(np.complex64).reshape(nr, nz) p comp.T tl -20.0 * np.log10(np.abs(p) 1e-12) return tl tl read_tvapm_tl(slope_2km.tl, nz201, nr201)需要留意的是 reshape 的 nz 和 nr 顺序取决于输出文件实际存储方式。如果画出来的曲线出现规律的锯齿不要急着做平滑先把 nr 和 nz 对调再读一次。锯齿通常只是数组形状反了平滑反而会把真实的干涉条纹抹掉。4. 水声模型关键参数海底、网格和吸收层的调法TVAPM 画出来很容易难的是让传播损失曲线和实验数据对上。我见过大多数错误都集中在三个地方垂直网格太粗、海底声速和水体不匹配、吸收层不够长。下面把这三块的调参思路拆开说。4.1 网格步长要跟着频率走别照抄同行的配置抛物方程方法的稳定性依赖波形分解精度。垂直方向如果每个波长分配不到六个网格点相速度误差就会让声场出现虚假色散。水平步长则受海底地形起伏约束地形变化越剧烈dr 越小。常用的经验配置如下频率垂直网格 dz (m)水平步长 dr (m)FFT 点数参考200 Hz1.225256800 Hz0.3105122 kHz0.1251024注意例子中的 dz 是按波长 1500/f 计算确保每波长 6 个采样点。FFT 点数则主要按水深和 dz 的比值来取比如水深 100 m、dz 0.3 m 时垂直点数约 333取 512 比较妥当。如果 FFT 点数小于这个值波数域的周期延拓会把边界反射带进计算区表现出来就是传播损失曲线呈周期性的波浪纹。4.2 海底声速、密度和吸收系数的设定方法海底参数不能随便抄。TVAPM 支持多层海底时每一层的纵波速度、横波速度、密度和吸收系数都会改变透射能量。常见的海底沉积层参数范围如下海底类型纵波声速 (m/s)横波声速 (m/s)密度 (g/cm³)纵波吸收 (dB/λ)淤泥15502001.51.0粗砂17003001.91.2紧凑砂18003502.01.5如果实验区海底是粗砂就不要用淤泥的参数。区别在低掠射角时非常明显海底声速比水层高时声波会穿透进海底再折射回来形成较强的泄漏能量如果海底声速低于水体声速则容易出现假的全反射导致传播损失被低估。频率较低时横波影响显著在砂质海底上忽略横波会让 500 Hz 以下的传播损失误差达到几个 dB。4.3 吸收层和边界反射的必调参数浅海计算里计算域底部必须有一个人工吸收层否则声场在底部边界反射回来会叠加出假干涉条纹。常见做法是让海底往下延伸 10 到 20 个网格点作为吸收过渡带衰减系数从 0 平滑增加到 5 dB/λ 左右。在我的工作版本里环境文件中可以这样写absorb_layer_start 95.0 absorb_layer_delta 20.0 absorb_layer_max 5.0这里 absorb_layer_start 是吸收层起始深度通常取大于最大水深absorb_layer_delta 是过渡带厚度absorb_layer_max 是底部最大衰减。如果发现距离末端传播损失突然冲高优先把 absorb_layer_delta 翻倍而不是减小 dz。吸收层的参数对单板机影响不大但调错会让后续所有参数白调。4.4 发散和振荡时的三个排查点运行 TVAPM 时最让人头疼的是计算到某个距离声场突然爆炸或振荡。第一个排查点是垂直网格 dz看它是否满足每波长至少 6 个点频率提高后忘记缩小 dz 是最常见原因。第二个排查点是海底声速如果 sediment_cp 小于水体声速就会出现非物理的泄漏模需要让海底纵波速度至少比当前水层声速高 20 到 50 m/s。第三个排查点是吸收层厚度吸收层太薄时能量在底部形成驻波表现为某一距离之后 TL 曲线周期起伏。检查顺序就按这个来不要一开始就怀疑源码有问题。5. 验证 TVAPM 结果的两个快招交叉对比和局部加密跑通只是第一步真正影响判断的是结果是否可信。没有海试数据的时候我会先用手里的射线模型和网格加密做交叉验证这两个方法能快速定位参数错误也能省下在场景上反复试错的时间。5.1 用 Bellhop 射线模型交叉验证平坦海域的传播损失找一个平坦海底区间深度设为 100 m声速剖面用等声速 1500 m/s分别用 Bellhop 和 TVAPM 计算同一个接收深度上的传播损失。Bellhop 命令很简单bellhop flat_100m.env然后把输出文件读出来和 TVAPM 的结果画在同一个图上。平坦海域里两者差异通常不超过 2 dB如果一致说明 TVAPM 的海底和吸收层参数没有系统性偏差。如果只在超过 1 km 的距离上出现差异优先检查横波参数和海底吸收而不是怀疑网格。5.2 只加密源深附近网格降低 pocketi1i 单次运行耗时在 pocketi1i 这类小设备上FFT 长度往往是最大瓶颈。一个实用技巧是只在声源附近加密垂直网格远离源深的位置逐渐变稀。这样能保证近距离的干涉结构够精细同时让垂直网格点数从 800 降到 500 左右单次运行时间大约能减少三分之一。如果源码支持垂直网格拉伸我一般把源深处的 dz 取 0.1 倍波长海底以下取 0.3 倍波长中间用线性过渡。要注意加密区域不能只覆盖几个点至少包括源深度上下各一个波长否则近轴误差反而变大。5.3 从输出中提取 1/3 倍频程 TL 曲线实验分析经常需要看某个频带的平均传播损失不建议直接对复声压做算术平均而是先转成单独频点的 TL再按能量做平均freqs [500, 800, 1000, 1250, 1600] tl_slices [read_tvapm_tl(ftl_{f}hz.tl, nz201, nr401) for f in freqs] tl_third_octave -10.0 * np.log10(np.mean(10**(-np.asarray(tl_slices) / 10), axis0))这段代码把所有频点的 TL 从 dB 转回线性能量取平均后再转成 dB。这样做比直接平均 dB 值更符合人耳和声呐处理的能量叠加逻辑。使用这段代码时要保证所有频点使用同一个网格范围和输出尺寸否则不同频点的空间点无法对应平均结果会出现空间错位。本文还有配套的精品资源点击获取