4. 水下地形测量与数字高程模型基础

地形辅助导航的核心思想，是把实时测得的局部海底地形与预先存储的数字高程模型（DEM）进行匹配，从而反推潜器的位置。因此，DEM 的质量直接决定了匹配定位的上限：如果参考图本身存在系统性偏差或分辨率不足，再精巧的匹配算法也无法给出可靠结果。本节围绕“DEM 从哪里来、怎样组织、如何评价”这条主线，讲清多波束测深原理、网格构建、插值方法与精度指标四个关键环节。

一、多波束测深：把一条线扩展成一个扇面

单波束测深仪每次只测船正下方一点，效率低且难以覆盖成片区域。多波束测深系统（MBES）在船底沿垂直于航向的方向布置一个发射换能器阵和一个接收阵，一次发射即可形成垂直于航向的扇形声波束，通常包含几十到几百个窄波束，覆盖宽度可达水深的 2～8 倍。船沿测线航行时，每个发射周期获得一条“测深条带”，相邻条带按一定重叠率拼接，就形成了面状覆盖。

测深值的计算依赖声速剖面：声波在水中的传播速度随温度、盐度和深度变化，系统用声速剖面仪（SVP）实测声速，将双程旅行时换算为斜距，再结合波束指向角解算海底点的三维坐标。同时，姿态传感器（横摇、纵摇、艏向）和 GNSS 提供换能器的瞬时位置与姿态，用于把每个波束点归算到统一的大地坐标系。这里必须明确坐标约定：本文统一采用“北-东-深”右手坐标系，深度以海平面为基准向下为正，后续 DEM 网格的行方向对应北向、列方向对应东向。

实际作业中，测线布设要保证相邻条带重叠 10%～20%，避免出现漏测的“数据空洞”；在坡度较大的区域，重叠率还需提高，否则边缘波束的测深精度下降会形成条带间的系统性台阶。此外，边缘波束的入射角大、 footprint 宽，精度明显低于中央波束，这也是后续精度评价需要分区考虑的原因。

二、从散点到规则网格：DEM 的构建

多波束原始输出是海量、不规则分布的离散点（点云），而地形辅助导航的匹配算法通常要求规则网格化的参考图，因此需要把散点重采样为按固定间距排列的格网，每个格网单元存储一个高程值。网格构建包含三个决定性的参数：

• 网格间距：决定 DEM 的空间分辨率。间距过大会丢失地形细节，降低匹配的可区分性；间距过小则会在数据稀疏区产生大量插值伪影，且存储与计算开销上升。
• 网格范围与原点：需与导航区域的坐标框架对齐，保证参考图与实时测量数据在同一坐标系下可直接比较。
• 坐标基准与深度基准：所有点必须归算到同一水平基准和同一垂直基准，否则会引入整体偏移。

一个常用的经验是：网格间距取测线间距的 1/3～1/2，或取平均水深的 1%～2%。例如在平均水深 100 m、测线间距 20 m 的测区，取 5 m 网格间距较为合理，既保留了主要地形起伏，又不至于让插值过度外推。

三、插值方法：不同地形用不同策略

散点到网格的转换本质是插值。常用方法各有适用场景：

对地形辅助导航而言，插值不仅要“填满网格”，更要避免抹平对匹配有利的地形特征。陡坡、海山、沟槽等区域若用强平滑的 IDW 或大范围克里金，会削弱地形可区分性，建议优先采用 TIN 线性插值或带地形约束的克里金。同时，对无数据覆盖的格网单元应明确标记为无效值，而不是用插值强行填充，否则会在参考图中制造虚假地形，诱导匹配算法产生错误定位。

四、分辨率与精度：两个不能混为一谈的概念

分辨率回答“能分辨多小的地形”，精度回答“高程值有多准”，二者相互独立。一个 5 m 网格的 DEM 可能分辨率很高，但如果声速改正错误或基准不一致，整体高程仍可能偏差数米。评价时通常关注以下指标：

• 垂直 RMSE：用独立检查点（checkpoint）与 DEM 插值高程比较，计算均方根误差，反映整体高程精度。
• 系统偏差 bias：检查点残差的均值，若显著非零，说明存在基准或声速系统性误差，需要整体改正。
• 水平分辨率：由网格间距、波束 footprint 和测线间距共同决定，可用地形功率谱或半变异函数分析其有效截止波长。
• 条带间一致性：相邻测线重叠区的测深差异，是发现姿态或声速问题的敏感指标。

以示例参数说明：某测区平均水深 100 m，网格间距 5 m，抽取 200 个独立检查点与 DEM 比较，若计算得 RMSE 约 0.8 m、bias 约 0.1 m，则说明高程精度较好且无明显系统偏差；若 bias 达到 1.5 m，则应优先排查深度基准或声速剖面是否存在整体性错误，而不是急于更换插值方法。需要强调，上述数值为教学示例，实际工程应以实测检查点统计为准。

综上，多波束测深提供原始观测，网格构建与插值把散点组织成可用参考图，而分辨率与精度指标则回答“这张图能不能用、能用到什么程度”。理解这四者的关系，才能在后续的地形匹配定位中正确设定参考图与实时测量的匹配尺度，并对定位结果的可靠性作出合理判断。