4. 海底边界建模与声学参数
海水并不是无限延伸的，声波最终会与海底相遇。海底对声场的影响，往往比海面更复杂：海面近似为理想压力释放边界，反射系数接近 −1；而海底是沉积物与基岩构成的多层介质，声波会部分反射、部分透射，并在沉积层内被吸收和折射。因此，海底边界建模的核心任务，是用一组可量化的声学参数描述“声波在这里如何被反射、透射和衰减”，并据此选择合适的地球声学模型。
一、海底的分层结构与关键参数
从声学角度看，典型海底自上而下可简化为：海水—沉积物层（可含多层）—基岩半空间。最上层是松散或半固结的砂、粉砂、黏土，孔隙率高、声速低；向下压实程度增加，声速与密度上升；再往下是砂岩、石灰岩等基岩，声速可能显著高于海水。对建模而言，每一层需要给出的基本参数是：
- 声速 c（m/s）：决定折射与临界角，沉积物中常见约 1450～1800 m/s，基岩可达 2500～6000 m/s；
- 密度 ρ（g/cm³ 或 kg/m³）：与声速共同决定特性阻抗 Z = ρc，控制反射强弱；
- 衰减系数 α（dB/λ 或 dB/m）：描述沉积物吸收，通常随频率近似线性或幂律增长，砂质衰减较小，泥质衰减较大；
- 层厚 h（m）：决定层内干涉与透射的频段特征。
反射系数由阻抗失配决定。对垂直入射，界面反射系数为 R = (Z₂ − Z₁)/(Z₂ + Z₁)。例如海水 ρ₁ = 1.025 g/cm³、c₁ = 1500 m/s，得 Z₁ ≈ 1.54×10⁶ kg/(m²·s)；若下层砂质沉积物 ρ₂ = 1.9 g/cm³、c₂ = 1700 m/s，则 Z₂ ≈ 3.23×10⁶，R ≈ 0.35，反射能量约 12%。这说明砂质海底是“较硬”的边界，反射较强；而高声速基岩界面 R 可超过 0.5。斜入射时还需用 Snell 定律处理，并注意是否存在临界角与全内反射。
二、三类建模方法及其适用条件
均匀半空间模型是最简形式：把海底视为单一均匀介质，只给 c、ρ、α 三个参数。它计算量小，适合中远距离、低频、海底分层不明显或只需一阶反射估计的场景。缺点是忽略了层内多次反射与干涉，在浅海低频段容易低估或高估传播损失。
分层介质模型把海底离散为若干水平层，逐层给出厚度与声学参数，用传播矩阵或递推方法计算总反射系数。它能刻画层内驻波、频率选择性透射和“高声速层下掩低速层”等实际结构。建模时要注意层厚与波长之比：当层厚远小于波长时，该层近似“薄层”，其影响可用等效参数合并；当层厚与波长可比时，必须保留分层。
粗糙界面模型处理海底并非理想平面的情况。常用方法是引入界面均方根粗糙度 σ 和相关长度 L，用 Kirchhoff 近似或微扰法计算相干反射损失，工程上也常用经验散射系数。粗糙度会削弱镜面反射、增强散射，在较高频率（波长与粗糙度可比）时影响明显。实际建模中，若缺乏测量数据，可按底质类型给出合理假设，并在结果中标注不确定性。
三、参数如何影响传播损失：一个示例推演
下面用同一套浅海条件对比三种海底假设，观察传播损失差异。设海深 60 m，海水 c = 1500 m/s、ρ = 1.025 g/cm³，声源与接收器均位于水下 30 m，频率 500 Hz，距离 5 km。三种海底分别为：A 均匀软泥（c = 1480 m/s，ρ = 1.4，α = 0.5 dB/λ）；B 均匀砂质（c = 1700 m/s，ρ = 1.9，α = 0.2 dB/λ）；C 两层结构（上层 5 m 砂质，下层基岩 c = 3000 m/s、ρ = 2.4、α = 0.1 dB/λ）。
按简正波或射线理论估算，A 的阻抗接近海水，反射弱、透射强，能量大量进入沉积层被吸收，5 km 处传播损失相对最大；B 阻抗失配明显，反射增强，损失相对减小；C 中薄砂层与基岩形成复合反射，在 500 Hz 附近层厚约 1.7 个波长，出现干涉起伏，损失介于两者之间且随距离呈非单调变化。这个示例说明：海底参数不是“背景常数”，而是直接改变声场结构的控制量；仅用均匀半空间可能把分层干涉造成的数十 dB 起伏抹平。
四、建模实施要点与结果判读
实际建模时，建议按以下顺序推进。第一步确定几何与频段：海深、声源与接收深度、工作频率决定波长，进而决定哪些层必须保留。第二步按底质类型赋初值：泥质取低声速、较高衰减，砂质取中高声速、低衰减，基岩取高声速、低衰减，并记录参数来源。第三步选择模型：若只需平均损失，用均匀半空间快速估算；若关注低频干涉或分层透射，用分层模型；若频率较高且海底起伏明显，叠加粗糙界面修正。第四步做敏感性检查：把声速、密度、衰减各扰动 ±5%～10%，观察传播损失曲线变化幅度，变化剧烈的参数就是后续最需要实测校准的对象。
结果判读时不要只看单点损失值，应关注损失随距离的斜率、干涉起伏的周期与深度、以及不同海底假设之间的分叉位置。若两条曲线在近场重合、远场分离，说明差异主要来自海底反射与层内衰减的累积效应；若全程存在周期性起伏，则提示分层干涉或粗糙散射在起作用。把这些特征与参数表对应起来，才能判断模型是否“物理上说得通”，而不是只追求数值吻合。
综上，海底边界建模的实质是“用分层结构 + 声学参数 + 界面特性”三要素逼近真实海底。均匀半空间给出基线，分层介质刻画结构，粗糙界面补充散射，三者按需组合，再通过敏感性分析锁定关键参数，才能让传播损失预测既有物理依据又可复现。