基于磁梯度测量的磁校准方法

摘 要：      利用三轴磁传感器进行飞行器航向测量或者磁异常探测之前， 需要对磁传感器进行磁干扰校准。 为了避免由于地磁场短时波动造成校准精度下降， 本文基于近地空间地磁场梯度分布特性， 提出一种利用磁梯度测量的磁传感器干扰校准方法， 建立了磁梯度测量系统磁干扰的24参数干扰模型， 基于单传感器的椭球拟合法多组任意旋转构建多元非线性方程组， 随后利用具有较高区域搜索效率的遗传算法进行非线性方程组的求解。 仿真结果表明， 和传统的椭球拟合校准方法相比， 该方法克服了地磁场短时波动的影响， 提升了磁梯度测量系统的校准精度， 其磁梯度差分测量误差可校准至［1.848 2 1.845 3 1.915 0］T nT。 该研究为航磁系统的载体干扰补偿和高精度测量提供了可靠保障。

关键词：     磁梯度测量； 磁干扰校准； 遗传算法； 三轴磁强计； 椭球拟合; 航向测量;  飞行器

中图分类号：      TJ760

文献标识码：    A

文章编号：     1673-5048（2024）01-0117-05

DOI： 10.12132/ISSN.1673-5048.2023.0100

0 引  言

三轴矢量磁传感器作为重要的航向测量器件， 被广泛应用于无人机、 导弹、 炮弹等飞行器的组合导航中［1］， 同时搭载三轴矢量磁传感器进行的航空磁探测在矿产勘探、 未爆弹探测以及航空磁探潜［2-4］等领域也得到各国科研人员的关注。

由于安装误差、 零点漂移、 灵敏度差异［5-6］以及载体的铁磁性物质与电子线路［7-9］对磁传感器的测量造成干扰， 为了获取精确的航空磁传感器测量数据， 需要在测量前及测量过程中对磁传感器进行磁干扰校准和补偿。

本文重点研究三轴磁传感器的静态磁干扰校准， 即在航空磁传感器系统安装完成后进行测量之前必须进行的校准工作。 其目的在于减小由于磁传感器三轴正交误差、 三轴灵敏度差异、 载体静态硬磁干扰［10-12］等误差源对磁传感器测量精度的影响。

现有的磁传感器静态干扰校准方法主要包括椭球拟合法［13］和系数校正法［14］等。 这些方法均基于环境磁场恒定的假设， 该假设只在人为构建的磁屏蔽环境中才能实现。 在外场试验过程中， 发现以地磁场为主的环境磁场在短时间内存在高达上百纳特斯拉的波动， 因此上述两种方法在进行磁传感器静态干扰校准时存在较大的系统误差， 在高精度磁异常探测领域需要更高精度的磁校准方法。

本文提出的基于磁梯度测量的磁校准方法， 借鉴了目前磁异常定位的主要手段， 即利用近地空间地磁场梯度较小的特点［15-18］， 构建磁梯度测量系统的24参数干扰模型， 并采用任意自由度旋转和遗传算法进行干扰方程组的建立和求解， 从而准确获得系统干扰参数。

1 磁传感器静态干扰模型

首先建立单个磁传感器静态干扰模型。 对单个三轴矢量磁传感器来说， 其静态干扰源主要来自四个方面： 安装误差导致的三轴非正交误差、 三轴标定系数误差、 磁滞现象导致的零点漂移误差、 传感器测量噪声及载体的硬磁干扰， 再依次对上述干扰源进行建模。

三轴矢量磁传感器由于安装工艺导致其三个实际测量轴不可能完全正交。 如图1所示， 假设其理想正交测量坐标系为OXYZ， 实际测量坐标系为OX′Y′Z′， 其中OZ与OZ′重合， YOZ与Y′OZ′共面， 且OY与OY′之间的夹角为α， OX′与XOZ面的夹角为β， OX′在XOZ面上的投影与OX的夹角为γ。

则磁传感器测得的磁场矢量Bm与实际磁场矢量Br之间的关系， 可由两个坐标系之间的旋转变换得到：

式中： R为正交转换矩阵。 实际磁传感器的非正交误差一般不大， 以Mag-03MS系列高精度磁通门传感器为例， 其非正交误差小于0.1°， 因此可将正交转换矩阵R近似为

考虑到磁传感器的3个测量轴标定系数差异， 设其三轴标定系数分别为Kx， Ky和Kz （Mag-03MS系列磁通门传感器的标定误差小于±0.5%）， 则标定系数矩阵K可表示为

由于磁滞现象以及信号处理电路导致的零点漂移和载体的硬磁干扰对磁传感器来说可以统一建模成硬磁干扰， 硬磁干扰是磁传感器静态干扰源的最主要成分， 可表示为固定矢量磁场Bh。

因此， 磁传感器静态干扰模型可表示为

Bm=K·R·Br+Bh+Bn（4）

式中： Bn为传感器测量噪声， 一般为高斯白噪声［19］。

2 磁梯度测量系统静态干扰模型

现有的矢量磁传感器静态干扰校准大多是利用单传感器静态干扰模型， 并基于地磁场在短时间内恒定的假设， 对传感器进行三维旋转后拟合计算得到校准参数。 但是在绝大多数自然环境条件下， 地磁场在短时内存在较大波动， 如图2所示。 该曲线为野外空旷地带获取的地磁场大小实测数据， 可见地磁场大小在短时仍有高达200 nT的波动， 显然与单传感器校准的前提假设不符。

研究表明， 我国地表附近的地磁场磁感应强度梯度量级为10-12～10-11 T/m［20］， 在工程应用中相对于磁异常梯度场可以将其近似为0。 基于该特性， 在进行磁异常探测时， 可以利用磁场差分梯度将地磁信号抵消。 为了规避地磁场短时波动造成磁传感器校准精度下降的问题， 利用近地空间地磁场梯度极小的特性， 构建了磁梯度测量系统的静态干扰模型。

如图3所示， 对于两个磁传感器构成的最简磁梯度测量系统而言， 由于传感器自身的装配误差和系统安装误差， 导致在传感器理想测量正交坐标系间存在一定的对准误差， 在进行磁梯度测量系统的静态干扰模型时， 需将其考虑进去。

通过多次采样均值滤波可以有效减小测量高斯白噪声带来的干扰， 因此本文忽略测量噪声的影响， 由式（4）可得两个磁传感器理想正交坐标系下的磁场强度分别为

考虑两个磁传感器理想正交坐标系间的对准误差， 需要引入旋转变换矩阵Rg。 另外， 由于磁传感器自身安装精度造成测量坐标系与机械坐标系的角度偏差（Mag-03MS系列磁通门传感器的Z轴安装偏差小于0.1°）， 需要分别引入旋转变换矩阵R1a和R2a， 因此磁梯度系统的静态干扰模型可表示为

Gr=R1a·B1r－Rg·R2a·B2r=R1a·（K1·R1）－1·（B1m－B1h）－Rg·R2a·（K2·R2）－1·

（B2m－B2h） （6）

式中： R1a·（K1·R1）－1可简化为包含9个未知参数的矩阵M1; Rg·R2a·（K2·R2）－1可简化为包含9个未知参数的矩阵M2。 又由地磁场梯度特性可将Gr近似为0， 则式（6）可简化为

M1·（B1m－B1h）－M2·（B2m－B2h）=0（7）

3 磁梯度干扰校准算法

以图4所示的磁梯度测量系统为例， 测量系统在任意旋转角度下的两个磁传感器测量值分别为B1m（θ）和B2m（θ）， 则无磁异常目标下的磁梯度矢量为

M1·（B1m（θ）－B1h）－M2·（B2m（θ）－B2h）=0（8）

因此， 将该磁梯度测量系统任意旋转多个不同角度， 即可得到如下的非线性方程组：

则磁梯度测量系统的静态干扰校准问题就转化为对上述非线性方程组的求解问题。 通过求解未知矩阵和未知向量M1，  M2，  B1h， B2h， 即可得到磁传感器在同一坐标系下的真实测量值， 并保证磁梯度测量系统的测量精度。

B1r=M1·（B1m－B1h）

B2r=M2·（B2m－B2h）（10）

对于简单非线性方程组的求解， 牛顿迭代法在待求零点邻近区域内具有平方收敛的性能， 但在较大区域内存在多解的情况下， 其很难得到最优解。 本文基于单传感器椭球拟合方法获取的先验知识， 利用遗传算法（GA）对上述非线性方程组进行优化求解。

遗传算法是计算机科学人工智能领域用于解决最优化的一种搜索启发式算法［21］， 是进化算法的一种。 进化算法最初是借鉴进化生物学中的一些现象而发展起来的， 这些现象包括遗传、 突变、 自然选择以及杂交等。 遗传算法的基本运算过程如下。

初始化种群： 构建染色体个体， 在给定区间范围内随机生成初始种群；

个体评价： 计算种群中各个个体的适应度；

选择运算： 选择的过程是一种基于适应度的优胜劣汰的过程， 指以一定的选择算法从种群中选择若干个体的操作；

交叉（基因重组）运算： 交叉运算是按照一定的交叉概率结合来自父代交配种群中的信息产生新的个体；

变异运算： 交叉之后的子代进行变异， 即子代个体按小概率扰动产生的变化； 种群经过选择、 交叉、 变异运算后得到下一代种群。

终止条件判断： 若种群中的个体满足给定的适应度条件， 则终止进化过程。

遗传算法的算法流程图如图5所示。

根据式（10）， 需进行优化的校准参数共有24个， 若将24个校准参数作为单个遗传算法里的染色体基因， 会造成算法效率下降以及多解问题。 为确保算法效率及校准精度， 在选择染色体基因及适应度函数方面， 本文做了相应处理。

（1） 将式（8）分解成3个彼此独立的非线性方程， 即

mi1（b1mx（θ）－b1hx）+mi2（b1my（θ）－b1hy）+

mi3（b1mz（θ）－b1hz）－［m′i1（b2mx（θ）－b2hx）+

m′i2（b2my（θ）－b2hy）+m′i3（b2mz（θ）－b2hz）］=0

（11）

式中： i=1， 2， 3。 mi1， mi2， mi3和m′i1， m′i2， m′i3分别为校准矩阵M1和M2的分量， b1mx， b1my， b1mz和b2mx， b2my， b2mz分别为两个传感器测量矢量B1m和B2m的分量， b1hx， b1hy， b1hz和b2hx， b2hy， b2hz分别为两个传感器硬磁干扰矢量B1h和B2h的分量。 每个方程分别用独立的遗传算法进行优化求解， 以减少各遗传算法染色体基因数量。

（2） 仿真结果表明， 磁梯度测量精度受硬磁干扰相比校准矩阵影响更小。 为进一步减少基因数量， 将单传感器椭球拟合得到的硬磁干扰矢量作为已知量， 进行适应度函数的计算。

（3） 三个遗传算法分别进行校准矩阵三个行向量6分量的优化， 传感器测量值任取12组旋转值， 适应度函数为求取12项式（11）左侧梯度矢量的模， 在椭球拟合获取的校准矩阵基础上， 取适当邻域范围作为染色体初始种群。

4 仿真试验验证

4.1 仿真条件

使用仿真软件进行磁梯度干扰校准方法的仿真试验验证， 并将本文方法与单传感器的椭球拟合校准方法进行比较。 根据实测环境地磁场的短时波动特征， 给出地磁场大小仿真曲线如图6所示， 并以当地实际地磁场情况为参考， 取仿真地磁场磁倾角为50°， 磁偏角为-5°（偏东为正）。 传感器的误差参数取Mag-03MS系列磁通门传感器的典型值。 传感器1的三轴非正交角度取0.1°， 三轴标定系数误差取0.5%， 硬磁干扰B1h=［1001000］T nT， Z轴安装偏差取0.05°； 传感器2的三轴非正交角度取0.05°， 三轴标定系数误差取0.4%， 硬磁干扰B2h=［－2001000］T nT， Z轴安装偏差取0.05°。 两个传感器间的坐标系对准误差旋转矩阵（依次绕X轴、 Y轴旋转0.5°）为