《Sensors》：Real-Time Reconstruction of the Eddy-Current Magnetic Field in a Local Region of a Maneuvering UUV Platform
编辑推荐：
当具有铝合金壳体的无人水下航行器（UUV）机动时，其在地磁场中的运动产生涡流磁场，对搭载的磁强计形成干扰。研究人员提出一种方法，离线提取低维空间基，并利用两个磁强计在线重构局部区域内的涡流磁场。离线阶段，沿三个机体轴施加单位激励的瞬态有限元法（FEM）仿真为每
当具有铝合金壳体的无人水下航行器（UUV）机动时，其在地磁场中的运动产生涡流磁场，对搭载的磁强计形成干扰。研究人员提出一种方法，离线提取低维空间基，并利用两个磁强计在线重构局部区域内的涡流磁场。离线阶段，沿三个机体轴施加单位激励的瞬态有限元法（FEM）仿真为每个激励通道提供一个准静态模态和一个弛豫模态及其阶跃响应。在线阶段，主磁强计输出的时间导数驱动前馈重构通过阶跃响应进行，单个辅助磁强计的差分测量给出一个幅度因子，逐样本估计以校正离线模型的幅度偏差。该方法既不要求角速率测量，也不要求校准机动。六个空间模态在整个传感器候选区域内表示偏航、俯仰和滚转机动的磁场，误差低于0.8%；对于同时激励全部三个通道的组合换向机动，误差为1.2%。在注入五类磁强计误差和两类模型误差的情况下，有利安装条件下区域重构误差为2.7–3.2%，典型安装条件下为6.1–7.3%。主磁强计处的残差，即矢量补偿后剩余的涡流干扰，分别为8–9 nT和14–15 nT。开环前馈重构产生9.7–10.9%的误差，而在其他六次机动上校准的完整Tolles–Lawson模型产生9–10 nT和13–14 nT的残差。在进一步测试中，一个独立的FEM模型以低10%的电导率、附加导体和新网格生成数据，而该方法保持原始离线模型。幅度因子收敛至0.90，即两个模型磁场之间的真实比值，残差保持在4–10 nT以内。在主磁强计重新定位后，所提方法无需重新校准的残差为10 nT，而在原始位置校准的Tolles–Lawson模型为25–40 nT。
地磁导航与磁测任务要求水下平台搭载高精度磁强计，而平台自身产生的磁干扰是限制测量精度的首要因素。机动平台的磁干扰通常分为三部分，即永久磁场、感应磁场和涡流磁场。前两部分的矢量校准流程已较为成熟。涡流磁场是平台在地磁场中运动时在导电结构中感应产生的，与机体系下地磁矢量的变化率成正比，仅在机动期间出现，通常被认为较小。然而在舰船尺度上，涡流磁场是仅次于铁磁效应的第二大特征源，铝和不锈钢等非磁性导电材料同样产生明显的涡流特征。对于具有铝合金壳体的无人水下航行器（UUV），壳体本身为非磁性，永久磁场和感应磁场主要来源于壳内设备，可通过矢量校准去除。但壳体在机动期间仍产生显著的涡流磁场，成为校准后的主导残余干扰，且该残余无法通过滤波去除。低通滤波曾被应用于假设涡流磁场由高频成分组成的场合，然而涡流磁场跟随地磁矢量的变化率，因而在0.1–1 Hz的机动频率范围内变化，而这正是机动期间被测地磁信号变化的频带，因此去除干扰的滤波器同时也会去除该信号。
现有涡流磁场处理方式遵循三种途径。第一种途径中，部分水下磁矢量补偿研究以低速时地磁矢量变化率低为由在模型中忽略涡流磁场，或将之视为测量噪声并通过低通滤波去除，这两个前提在姿态机动期间均不成立。第二种途径中，航磁补偿采用Tolles–Lawson（T-L）模型描述涡流磁场，以磁强计输出的时间导数作为回归量。该模型从八个涡流系数发展到含附加干扰项的扩展模型，再到具有在线更新的3×3系数矩阵。然而系数的在线更新依赖于地磁矢量的参考，该参考由姿态信息或标量总场约束提供。T-L模型的系数通过校准机动辨识，仅描述磁强计安装点处的测量畸变。因此，当在设计阶段比较候选安装位置时，或在任务之间随载荷一起更换位置时，每个位置都必须重新校准。第三种途径中，舰船特征领域进行高保真离线正向计算和涡流设施的测量，并探索了从船体参数统计预测特征的方法。这些提供了涡流磁场的空间分布，但单次有限元法（FEM）求解需要数小时，且模型与真实结构在电导率、焊缝和几何细节方面的差异直接作为误差进入计算结果。
上述三种途径均未解决的问题是：涡流磁场的空间分布和幅度无法同时在线获取。涡流磁场由分布在整个导电结构上的源产生，任意位置的场由结构的瞬态响应和位置本身共同决定。T-L模型仅在校准点成立且不含空间结构。离线正向计算含有空间结构但无法实时更新，且其幅度受模型与真实结构之间差异的支配。因此，在局部区域内在线重构涡流磁场的时空分布，同时使其幅度独立于离线模型的绝对精度，仍是一个未解决的问题。
为解决该问题，研究人员提出一种涡流磁场实时重构方法。从沿三个机体轴施加单位激励的瞬态FEM仿真中离线提取低维空间基及其模态的阶跃响应。在线阶段，主磁强计输出的变化率作为前馈重构的激励，单个辅助磁强计的差分测量提供幅度因子的在线估计，以校正离线模型与真实平台之间的幅度偏差。
该方法的两要素均有先例。基于快照的模型降阶为瞬态涡流问题提供低维空间基，但这些基被用于加速正向计算而非从磁强计数据在线重构。在流体动力学中，gappy本征正交分解通过求解预计算基的系数从少量传感器重构非定常场，此类重构的传感器位置从基本身选取。同样，闭环消磁系统磁强计的部署区域已从计算的特征变化中选取。成对磁强计的差分测量已被用于地磁矢量测量中抵消平台干扰和小卫星上。在后一应用中，两个传感器处干扰的比值在干扰与环境场不相关的假设下从数据中估计。姿态信息也被用于抑制干扰成分。这些方案估计的是干扰本身而非离线空间模型的幅度。将降阶模型的预计算空间基与来自单次差分测量的在线幅度校正相结合，尚未在涡流磁场中尝试。
本研究的主要贡献如下：（1）研究人员使用沿机体轴三个单位激励和四次独立机动（包括一次三个激励通道同时作用的机动）的瞬态FEM数据，验证了平台在传感器候选区域内的涡流磁场由六个与机动无关的空间模态张成，并验证了各模态系数对任意激励的响应完全由其阶跃响应确定。（2）研究人员建立了以主磁强计为参考的差分测量模型。以机动前的稳态姿态区间作为偏移参考，从单个辅助磁强计的差分测量中逐样本估计一个幅度因子。该方法既不要求角速率测量也不要求校准机动，磁强计的系统误差和离线模型的幅度误差均仅通过幅度因子进入结果。（3）研究人员从材料手册和器件数据表给出的参数构建误差预算。在磁强计误差和模型误差同时注入的情况下，以区域重构误差和主点残差两个指标，将所提方法与T-L模型、开环前馈重构和全系数估计进行比较。
在技术方法方面，研究采用三维瞬态FEM（COMSOL Multiphysics 6.1，约7.0 m长铝合金UUV简化实体模型，6063-T83铝合金，电导率约3.03×10? S/m，291万二阶矢量单元，1840万自由度）生成单位激励集（沿三个机体轴的斜坡-保持背景场）和机动集（偏航、俯仰、滚转及组合换向机动），用于模态提取和阶跃响应拟合；阶跃响应以多指数函数在固定时间常数网格上拟合，在线计算归结为一组一阶线性滤波器的递推；幅度因子通过堆叠差分测量方程以噪声协方差加权逐样本求解；并以一个独立FEM模型（电导率降低10%、附加铝圆柱、重新划分网格）作为地面真值进行验证。样本数据来源于上述FEM仿真，不涉及实测队列。
研究结果如下：
4.1 仿真设计与低维模态库。通过单位激励解确定各激励方向的模态数，两个模态将单方向重构误差降至1.4%以内，三个模态降至0.2%以内；三个方向子空间之间的主角为82.6°–90.0°，表明不同方向的模态相互独立。模态库通过三项独立检验：从斜坡阶段提取的模态重构未用于提取的衰减阶段，误差不超过0.8%；仅从轴线数据预测偏移线上的场，误差为0.08–0.49%；将公共基投影到四次未用于提取的机动数据集，六个模态在候选区域体积网格上的重构误差为偏航0.60%、俯仰0.78%、滚转0.61%，增至九个模态仅降至0.57–0.64%；组合换向机动的六模态误差为体积网格1.24%、两条线1.00%，恒速段为0.7–0.8%，在0.1 s斜坡和换向期间升至3–7%。阶跃响应的七项指数拟合对準静态模态的相对RMS误差为0.3–2.7%，对弛豫模态为1.4–7.1%；以精确激励代入重构未用于拟合的机动，主点残差偏航4.7 nT、俯仰5.2 nT。
4.2 传感器配置、误差注入与参考方法。主磁强计置于候选区域前端轴线点(?2.70, 0, 0) m；按第3.3节准则，在2244个可安装网格点中，涡流场与主点之差的RMS最大值650 nT出现在候选区域后部隔舱接口框架附近的径向点(?0.90, ?0.15, 0) m。注入的误差包括模型幅度缩放误差（最大±20%）和形状误差（通过对单位激励快照施加随机扰动重新提取模态和阶跃响应生成）。参考方法包括：在主点采用完整矢量形式的T-L模型（18个系数，由六次其他机动最小二乘辨识）；开环前馈（去除幅度因子校正的消融参考）；全系数估计（去除前馈结构的消融参考，逐样本求解全部六个系数）；以及成对磁强计自适应干扰抵消。比较在理想情况、有利安装和典型安装三种主要情况下进行，每种情况200次蒙特卡罗运行。
4.3 在线重构结果。在理想情况下，开环前馈的残差略低于所提方法，因后者在启动阶段乘以尚未收敛的幅度因子；完整T-L模型的残差为8.1–9.5 nT。注入误差后方法排序改变：有利安装下，所提方法偏航/俯仰/滚转的区域重构误差为2.7%/3.2%/2.9%，开环前馈为10.5%/10.7%/9.7%，全系数估计为12.0%/8.4%/12.1%；主点残差所提方法为8.0/9.3/8.4 nT，T-L模型为10.3/9.3/10.4 nT。典型安装下，所提方法为6.1%/6.3%/7.3%和13.6/15.0/15.4 nT，开环前馈为10.6%/10.9%/9.9%，全系数估计为45.9%/28.1%/46.3%，T-L模型为13.1/13.7/13.1 nT。组合换向机动下，所提方法在有利安装下为5.0%和11.6 nT，典型安装下为5.8%和14.1 nT。自适应抵消在所有情况下产生约5 μT（单轴机动）和1.2 μT（组合机动）的残差，因其比值估计所依据的读数由机动期间地磁矢量的变化主导，违反了干扰与环境场不相关的假设。所提方法的误差仅通过幅度因子进入，典型安装下干扰仍被抑制92–93%。在理想情况下，所提方法残差集中在机动开始后前0.3 s内；偏航0.40 s后、俯仰和滚转0.45 s后，200次运行的幅度因子均值保持在1±0.01以内；组合换向机动中幅度因子在启动和角速率换向之后两次偏离并各在0.55 s后恢复。噪声从1 nT增至10 nT使幅度因子标准差增大十倍，但残差仅增加1.1 nT。在主磁强计重新定位至候选区域中部轴线点(?1.70, 0, 0) m且辅助位置和离线模型不变时，有利安装下残差为9.6/10.0 nT，而在原始位置校准的T-L模型直接转移至新位置产生39.6/24.7 nT。标量补偿残差在主点为0.6 nT，补偿前总场涡流干扰RMS为10.1/20.3 nT。
4.4 误差预算与鲁棒性。1°轴失准和0.5%比例因子误差各使残差增加约3 nT，20 mm位置误差改变小于0.1 nT，2 ms传感器间时间偏移增加约4 nT；对于19°旋转角和10 nT的单因素残差限值，容差约为1°轴失准、0.5%比例因子和2 ms时间偏移。偏移去除使失准和比例因子误差仅通过地磁矢量的变化泄漏。幅度缩放误差ε_a为±20%时，所提方法残差为6.9–7.4/7.4–9.1 nT，开环前馈为19.3–20.7/23.5–24.0 nT；形状误差ε_s单独作用时所提方法为9.9/9.6和16.4/14.7 nT，开环前馈为8.1/8.1和13.9/13.2 nT；所提方法优于开环前馈的条件是幅度缩放误差超过形状误差，交叉点约在两者均为5%处。独立FEM模型测试中，其涡流场在兩個磁强计位置均为标称值的0.90倍，幅度因子收敛至0.899，主点残差在1 nT噪声下为4.2 nT，在有利、典型和校准典型安装的磁强计误差下为5.1/9.8/6.6 nT；区域误差在全部2845个空气点为15.8–16.8%，但远离圆柱处为2.5/3.0/6.3%。外部异常影响通过叠加量化：只要异常差值的RMS不超过50 nT，残差增加小于3 nT（俯仰为5 nT）；无异常时差分测量中方程(10)未解释部分的RMS为1–4 nT，而残差增加超过3 nT的每种情况下均超过20 nT。角速率增至更高值时，涡流干扰升至578/638 nT，所提方法在理想、有利和典型情况下的区域误差为1.6%/2.2%、2.0%/2.6%和3.7%/4.0%，抑制比为95–98%。
讨论部分总结如下：所提方法的优势在于区域重构、无需重新校准的重新定位以及无需校准机动，而非校准点处的精度。T-L模型在其校准点与所提方法相当，但其感应项吸收了比例因子和失准泄漏，且其感应项与涡流项的分离需要在每个安装位置进行绕全部三个轴、多于一个角速率并含角速率换向的校准机动；在原位置校准的模型用于新位置时残差升至25–40 nT。所提方法的误差仅通过幅度因子进入，五个误差来源的单项增加之和略大于组合增加，表明误差源部分独立；将比例因子误差从0.5%校准至0.1%是收益最大的单项措施。启动阶段的差距由涡流建立期间辅助位置处的模态失配而非噪声引起，先验正则化、阈值切换和累积形式切换均未能消除该差距。研究结论为：对于低频机动和线性材料，平台涡流磁场由六个与机动无关的模态张成，即每个机体轴激励通道一个准静态模态和一个弛豫模态，这些模态以低于0.8%的重构误差表示未用于模态提取的偏航、俯仰和滚转机动，并以1.2%的误差表示三通道同时作用的组合换向机动；各模态系数的时间响应以多指数形式的阶跃响应表示，在精确激励下前馈重构误差约为2%。主磁强计输出的变化率可作为激励，离线模型的幅度偏差通过由单个辅助磁强计差分测量获得的幅度因子逐样本校正；在注入五类磁强计误差和两类模型误差的情况下，有利安装下区域重构误差为2.7–3.2%、主点残差为8–9 nT，典型安装下为6.1–7.3%和14–15 nT，组合换向机动为5.0%和12 nT与5.8%和14 nT；相同条件下开环前馈的区域误差为9.7–10.9%，全系数估计为8.4–46%；主磁强计移至另一位置而无需重新校准时残差为10 nT；在独立FEM模型测试中幅度因子收敛至0.90，主点残差保持在4–10 nT，区域误差仅在附加导体0.125 m范围内上升。磁强计误差仅通过幅度因子进入结果，对于不超过10 nT的单因素残差，容差约为1°轴失准、0.5%比例因子和2 ms传感器间时间偏移；离线模型±20%的幅度误差使残差变化小于3 nT，10%的形状误差使其升至15–16 nT。上述结论针对线性材料、低频机动和单一平台的FEM模型获得，形状类模型误差、绕时变轴的机动、涡流建立期间辅助位置的模态失配、失准漂移、近距铁磁物体以及真实磁强计更高采样率下的在线行为等六方面仍有待解决。
