李 毅，刘忠乐，陈 俊，胡家文，文无敌

(1. 海军工程大学 兵器工程学院，湖北 武汉 430033；2. 华中科技大学 光学与电子信息学院，湖北 武汉 430030)

0 引 言

铁磁性目标物在地磁背景场中会被磁化，进而改变背景场原有的分布，产生磁异常现象。通过对磁场的测量，将测量的磁场信息进行一定的处理来判断目标物的有无、实现目标物的定位。磁异常探测技术中测量到的磁场信号含有丰富的噪声。目前磁异常探测常用手段主要是传感器和被探测物相对运动直到彼此错开，再利用生成的完备数据进行建模求解。然而现在的磁异常探测更需要能够及时预判被探测物的具体定位。

如何在有噪声干扰且信噪比极低的情况下，实现对磁性目标的预判定位是磁异常探测技术领域中一个比较棘手的问题。

正交基函数是磁异常探测技术中检测磁性目标物的经典算法。Frumkis 等[1]基于正交基函数的信号能量评估，提出了OBF 算法。通过该算法可以提高目标信号的信噪比，从而判断目标信号的有无。Qin 等[2]利用地磁坐标系3 个轴向半径的梯度信息，进一步对OBF算法进行改进。使其对目标磁异常检测有了更好的表征，显著提高信噪比。Ginzburg 等[3]利用2 个磁传感器测量总场梯度信息，通过改进的正交基函数对磁场信号进行分解，以提高信噪比和信号检测的特性。OBF算法有着噪声与正交基函数不相关的特性，因而有着对噪声不敏感，有效抑制噪声和提高信噪比的优点。但目前对OBF 算法的研究大部分用于判断磁异常信号的有无。利用OBF 强大的抗噪声性来实现磁目标的定位还不是很普遍。现阶段，利用磁场梯度信息对磁性目标物进行定位的研究越来越多。因其能有效规避背景场的干扰，而广泛应用到磁异常探测定位中。Sui 等[4]通过从旋转磁盘上的单轴磁传感器中提取二阶和三阶梯度张量，以此来实现磁偶极子的高效定位。Huang等[5]将磁梯度张量和吃水深度结合起来，实现水下航行器的有效定位。Fan 等[6]提出一种基于总磁场梯度定位目标的快速线性算法，实现了静态磁目标的高速定位。大多数利用梯度信息进行目标的定位方法，虽然有着不易受到背景场干扰，能较好描述磁异常场的优点。但是这些方法都是在背景噪声不明显，信噪比较高的情况下才能实现磁性目标的准确定位。

Birsan[7]通过传感器阵列，将正交基函数与单一梯度信息结合，将普遍用于判断磁异常有无的OBF 算法运用到磁异常定位中。目前这种利用OBF 实现磁异常探测定位的算法，针对的是信噪比大于0 的情况。对于信噪比小于0，甚至更低的复杂噪声背景下的定位效果还没有研究。

本文利用声呐系统提供测量的大致方位和距离，结合正交基函数天然的规避噪声的特性，将OBF 算法结合L-M 算法改进，成为可以运用到磁异常探测定位领域中的一种基于声磁复合的磁探测高效算法，有效提升声呐定位精度。

1 OBF 算法

在传统的MAD 系统中, 如图1 搭载有2 个标量磁强计的探测平台沿着笛卡尔坐标系的X轴方向移动，路径的方程见式（1）。其中x1为探测平台在X轴的分量初始值，且为负数，探测平台沿着平行于X轴的直线匀速以v向X轴正方向行进，经过声呐的初步探测定位点，把目标暂时假定为坐标系原点[8-13]。

图1 航空磁异常探测Fig. 1 Airborne magnetic anomaly detection

假定铁磁目标静止不动。当传感器距离铁磁目标大于目标尺寸3 倍时,铁磁目标产生的异常场可以视为偶极子场，并且可以描述为：

式中，M为磁偶极矩，r为偶极子到传感器的位置矢量， µ0是真空的磁导率。磁强计在轨道上以匀速v平行于x轴的直线上运动:，在t=t0时刻，传感器处于离目标最近的接近点R0(也即为最短横距CPA 值)。将代入式（2）化简，将的每个分量都可以写成3 个Anderson 函数的线性组合

其中， τ为归一化时间，t0为传感器和被测物最近的时候（最短横距为CPA）。

将B化简为：

其中：

此时S可以由3 个线性无关的基函数的加权和表示。

可以使用朗斯基行列式证明这3 个函数之间是线性无关的。

代入测量的数据，可以解出y0、R0和ai(i=1,2,3)这5 个未知数，而且上述过程假定地磁背景场是均匀的，使用2 个标量传感器取均值来计算更加准确。

2 探测定位流程框图

本文设计探测流程如图2 所示。由2 个标量磁性传感器测得的数据加入不同噪声，再经过数据预处理，结合声呐给定初始位置，进行L-M 算法的迭代处理，从而得出磁异常的具体参数，并根据结果分析验证准确性。

图2 探测定位流程框图Fig. 2 Probe the locating process block diagram

磁异常目标位置坐标为（x0,y0,z0），2 个磁传感器相距1 m，且S1和S2平行于水平面移动，磁强计的排布如图3 所示。将移动的方向定义为X轴正方向，则2 个磁强计采集的磁异常信号可以表示为：

图3 磁强计的排布Fig. 3 Magnetic sensor layout

对于测量轨迹上的偶极子异常信号S，其在X方向上的梯度为：

其中：

可以使用朗斯基行列式证明这4 个函数之间是线性无关的[14-17]。

沿X轴方向的梯度可以由2 个传感器的差分近似表达为：

至此根据x方向上的梯度得到了如上方程，未知数仍是x0，y0，z0，R0，ai(i=1,2,3)，因而可以根据最小均方误差的标准，使用Levenberg-Marquardt 算法求解[18-20]。但由于z0只存在于(zs-z0)2这一项中，即z0=-|z0|或者z0=|z0|都满足最优解，在此不直接用LM 算法求出z0，而是使用得到的R0、x0以及z0＜zs这一条件唯一确定z0。

3 实验验证

在数值验证部分将地磁场视为背景场，且认为信号中没有噪声，2 个磁强计的数据由传感器测量得到，如图4 所示。

图4 两个磁强计磁场数据Fig. 4 Magnetic data of two magnetic sensors

可以看出由于间距比较小，2 个磁强计的数据差别非常小。

计算过程中zs=100 m，xs=100 m，偶极子位置为(0,0,0)，地磁场的磁倾角I=70°，磁偏角D=81°。传感器间距ΔY=1 m。使用Levenberg-Marquardt 算法求解，但由于z0只存在于(zs-z0)2这一项中，即z0=-|z0|或者z0=|z0|都满足最优解，在此不直接用LM 算法求出z0，而是使用得到的R0、x0以及z0＜zs这一条件确定唯 -z0。

可以看出在声呐系统辅助下，将初始值设定为0.1，传感器和被测物垂直方向距离定位110 m，在运算过程中将b0，b1，b2都定初值为1 时（b0，b1，b2为L-M 算法中自主定义的3 个参数），Levenberg-Mar quardt 算法经过501 次迭代，拟合效果很好，如图5所示。

图5 测量磁场和拟合磁场对比Fig. 5 The measurement of magnetic field and the parallel field contrast

接下来进行磁场梯度的验证，按照L-M 算法计算的结果(197.90,296.38,-52.19)，估计的误差较小。使用求解出的系数代入式(12)、式(8)与式(9)、式(14)，得到的梯度与原始数据对比如图6 所示，拟合的结果与真实值非常接近。

图6 X 轴方向上的梯度Fig. 6 The gradient of the X-axis

再考虑偶极子朝向的不同，角度的变化也会导致磁场的结果随之发生改变。随机设置32 组偶极子朝向(θd,βd)组合，并生成仿真数据，使用算法反演磁异常位置结果如图7 所示。

图7 随机偶极子朝向时的预估位置和实际位置对比Fig. 7 The estimated position of the random dipole is compared to the actual position

可知，在R0=269.26 m的情况下，对于不同的(θd,βd)组合，x0的求解误差在±25 m，y0和z0的求解误差在 ±5 m，说明本文方法具有较好的求解精度。

而在实际探测到的信号一般包含各种噪声。为了检验本文方法在低信噪比情况下的表现，现在对仿真生成的偶极子信号人为添加不同强度的噪声，以模拟不同信噪比下采集到的信号。所加的噪声是标准的高斯白噪声，此处信噪比定义为：

由于磁强计之间的距离比较近，差分求梯度能在一定程度上减弱噪声的影响。图8 为信噪比在-10～10 dB的范围内，式(6)中R0的计算结果。

图8 不同信噪比下 R0的计算结果Fig. 8 The results of the different noise ratio of R0

可以看出，本文的解算方法对噪声具有一定的适应能力。在信噪比为-10 dB的情况下R0的求解相对误差小于5%。说明算法对于低信噪比的环境具有很好的适应能力[21-25]

4 结 语

本文提出一种信号处理技术，用于定位和识别一个模拟为磁偶极子的目标。该算法同时利用了在配备2 个磁强计的飞机上同时测量的总磁场和标量梯度。对于不同的实验情况和噪声水平，这2 个量都可以通过一组安德森型函数很好地建模。从大量的计算机模拟得到的结果表明，对目标位置的估计和磁矩精度的预测是可以接受的，即使在相当低的信噪比-10 dB的情况下，从2 个磁力计的信号依然可以较好反演出实际的磁异常信号位置。