3. 水下导航传感器及其误差特性

组合导航的精度上限,首先由参与融合的传感器决定。EKF 只是把各传感器的信息按不确定度加权融合,它无法凭空创造观测信息。因此理解每类传感器的测量原理、输出物理量、误差来源与误差随时间的演化规律,是设计状态量、建立量测方程、设置噪声矩阵的前提。本节按“惯性传感器—深度计—DVL—磁力计—USBL”的顺序,逐一分析其特性,并给出可用于 EKF 建模的误差描述方式。

一、惯性传感器:加速度计与陀螺仪

惯性测量单元(IMU)由三轴加速度计与三轴陀螺仪组成,是水下组合导航的核心。加速度计测量比力,即载体所受非引力外力与质量的比值;陀螺仪测量载体相对惯性空间的角速度。二者输出经姿态矩阵投影与积分,才能得到速度、位置与姿态。这里必须明确坐标系约定:本文采用右前上(RFU)机体坐标系,导航系取北东地(NED),陀螺仪输出为机体系相对惯性系的角速度在机体系的投影,加速度计输出为机体系比力。

惯性传感器的误差分为两类。第一类是确定性误差,包括零偏、标度因数误差、安装失准角,可通过实验室标定补偿。第二类是随机误差,无法完全标定,只能用随机过程描述。工程上最常用的是 Allan 方差分析,从静态长时数据中辨识出角度随机游走(ARW)、速度随机游走(VRW)、零偏不稳定性与零偏随机游走等系数。在 EKF 中,通常把陀螺零偏与加速度计零偏作为状态量增广,用一阶马尔可夫过程或随机游走建模:

δb_g(k+1) = δb_g(k) + w_g,   w_g ~ N(0, Q_g·Δt)
δb_a(k+1) = δb_a(k) + w_a,   w_a ~ N(0, Q_a·Δt)

需要强调的是,惯性误差是随时间累积的:姿态误差会通过比力投影污染速度,速度误差再积分成位置误差,形成“姿态—速度—位置”的误差传播链。这正是纯惯导必须依靠辅助传感器定期校正的根本原因。教学示例中,若某 MEMS 陀螺零偏稳定性为 10°/h,则 1 分钟内由零偏引起的姿态漂移约 0.17°,这个量级足以说明为何水下任务不能长时间仅靠惯导。

二、深度计:最稳定的绝对观测量

深度计通常为压力式,通过测量静水压 p 反算深度:h ≈ (p − p₀)/(ρg),其中 p₀ 为水面大气压,ρ 为水密度。它的关键特性是误差不随时间累积,只与压力传感器精度、水密度假设和波浪引起的动态压力有关。因此深度几乎总是作为 EKF 的直接观测量,用于约束天向(NED 下的 D 轴)位置与垂向速度。建模时量测方程可写为 z_depth = h_true + b_d + v_d,其中 b_d 为缓慢变化的深度零偏,可作为状态量或在线估计。需要注意的是,若把深度直接当作“高度”使用,必须统一符号:NED 下深度为正向下,与高度符号相反。

三、DVL:速度观测的主力

多普勒计程仪(DVL)向海底发射声波,利用回波多普勒频移测量载体相对海底(底跟踪)或相对水团(水跟踪)的三维速度。底跟踪时输出的是相对大地(海底)的速度,是极有价值的速度观测;水跟踪时输出相对水流速度,会引入洋流偏差。DVL 的主要误差包括:标度因数误差(约 0.5%~2%)、安装偏角、波束倾角引起的误差放大,以及底跟踪失锁时的野值。在 EKF 中,DVL 速度通常作为量测,量测方程需通过姿态矩阵把机体系速度转换到导航系:

v_nav = C_b^n · v_dvl
z_dvl = v_nav + v_noise

实操中应设置野值剔除逻辑:当 DVL 输出的有效标志位为假、或速度跳变超过阈值(如 2 m/s)时,不更新该帧量测,避免污染滤波器。教学示例中,若 DVL 标度因数误差为 1%,在 1 m/s 航速下会带来 0.01 m/s 的速度偏差,短时影响小,但长时间会累积成航迹偏移,因此高精度任务中常把 DVL 标度因数也纳入标定或在线估计。

四、磁力计与 USBL:航向与位置的辅助约束

磁力计测量地磁场在机体系的分量,经倾斜补偿与磁偏角修正后可解算航向角,为惯导提供航向绝对基准。其误差主要来自硬磁与软磁干扰、载体铁磁材料以及局部磁异常,水下近底航行时尤其容易受干扰,因此使用前需做椭球拟合标定,并在 EKF 中给予较大的量测噪声。USBL(超短基线)通过测量声学应答器到基阵的时延与相位差,解算目标相对基阵的位置,再结合母船 GPS 得到水下目标的绝对位置。它的误差随距离增大而增大,且存在声速误差与多径影响,通常以较低频率提供位置观测,用于抑制惯导位置漂移。

水下组合导航传感器特性与误差传播关系 惯性传感器 IMU 加速度计:比力 f_b 陀螺仪:角速度 ω_b 零偏 / 标度因数 / 随机游走 误差随时间累积 输出:姿态、速度、位置(递推) 惯导误差传播链 姿态误差 → 比力投影偏差 → 速度误差 → 位置误差 EKF 增广状态:δb_g、δb_a 需辅助观测定期校正 辅助传感器 深度计:绝对深度(不累积) DVL:底跟踪速度 磁力计:航向基准 USBL:绝对位置(低频) 误差有界,可作量测 EKF 融合层 状态:位置、速度、姿态、零偏 量测:深度 / DVL 速度 / 航向 / USBL 位置 按噪声协方差加权,抑制惯导漂移 误差特性对比: IMU:随时间累积,需增广零偏状态 深度计:不累积,约束垂向 DVL:标度因数与野值需处理 磁力计:受磁干扰,噪声取大 USBL:误差随距离增大

五、误差建模的工程要点

把上述特性落到 EKF 实现上,有三点值得特别注意。第一,噪声矩阵要与传感器实际水平匹配:Q 阵过小会导致滤波器过度信任惯导、对辅助量测反应迟钝;Q 阵过大会使估计抖动。建议用 Allan 方差得到的 ARW、VRW 直接换算过程噪声,用静态数据的标准差设定量测噪声 R。第二,量测更新要有时序与有效性管理:深度、DVL、磁力计、USBL 的更新频率不同,应按到达时刻依次做量测更新,并对每类量测设置有效性标志与野值门限。第三,坐标系与符号必须全程一致:NED 下深度向下为正、航向以北为起点顺时针为正,DVL 机体系速度经 C_b^n 转换时姿态矩阵的定义要与状态量中的姿态误差定义对应,否则融合结果会出现系统性偏差。

实践建议:在正式下水前,先做一次静态与匀速直航的传感器数据采集,用静态段估计零偏与噪声,用直航段检查 DVL 标度因数与磁航向偏差,再把辨识结果写入 EKF 的初始协方差与噪声矩阵,可显著缩短现场调试时间。

综上,惯性传感器提供高频递推但误差累积,深度计、DVL、磁力计与 USBL 提供误差有界的绝对或相对观测。EKF 的价值正在于用后者的“有界”去约束前者的“累积”,而这一约束能否有效,取决于对每类传感器误差特性的准确刻画。