本发明属于地磁探测技术领域，具体涉及一种采用地磁总场梯度阵列对磁性目标定位的方法。
背景技术
地磁场是地球的一个天然的物理场，它有各种不同的起源，由不同变化规律的磁场成分叠加而成。按照场源位置划分，地磁场可以分为内源场和外源场。如果考虑地磁场随时间的变化特征，随时间变化较快的地磁场成为地球的变化磁场，随时间变化较慢或者基本不变的地磁场成为地球的稳定磁场。
一个带有磁性的目标产生的磁场会导致空间地磁场分布的变化，从而能够产生磁异常。当观测点与目标距离大于目标尺度2-3倍及以上时，磁性目标在观测点产生的磁场一般被看作磁偶极子远场。
对磁性目标定位时，可采用能够测量地磁三分量的矢量磁传感器或者能够测量地磁总场的标量磁传感器。矢量磁传感器安装和使用相对复杂，初始姿态要严格校正，工作姿态要实时测量。矢量磁传感器的角度误差为0.05°时，测量的地磁误差大概为50nt。因此，实时补偿矢量磁传感器姿态变化对测量的影响是必需的，也是有一定难度的。一般情况下，市场上常见的矢量磁传感器(如普通磁通门磁力仪)分辨率相对偏低(相对于标量光泵磁力仪)，因此基于矢量测量的目标定位距离不能太长，目标磁矩不能太小。
综上所述，现有技术中存在分辨率低，安装使用复杂，探测距离近，被探测物体要求较高，适用范围小等问题。
技术实现要素：
本发明的目的在于对磁性目标进行准确快速简便的追踪和定位。首先构建由七个标量磁传感器(光泵磁力仪)组成的传感器阵列，然后基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系。采用矩阵变换，把求解参数消元缩减到三个。建立适应度函数，采用粒子群算法对磁性目标定位。
一种采用地磁总场梯度阵列对磁性目标定位的方法，包括以下步骤：
(1)构建由七个标量磁传感器组成的传感器阵列；
(2)根据所述传感器阵列，基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系；
(3)根据所述地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系，采用矩阵变换，把求解参数消元缩减到三个；
(4)建立适应度函数，用粒子群算法计算上述适应度函数的最小值，求解目标位置。
所述构建由七个标量磁传感器组成的传感器阵列，包括：
阵列中心为坐标原点o，七个传感器位于原点及三个坐标轴上，阵列中各传感器的坐标为t0(0,0,0)，
其中，l为每一坐标轴上两个远端传感器的距离。
所述根据所述传感器阵列，基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系，包括：
每台标量磁传感器测量值ti都是地磁正常场be和目标异常场bi相叠加的总场模，即
ti＝|be+bi|
u＝[cos(i0)cos(d0)cos(i0)sin(d0)sin(i0)]＝[abc]
远场条件表述为：
3la≤r、l＜r、bi＜be
在远场条件下：
ti＝|be+bi|＝be+u·bi
其中，u为be的方向向量，i0为地磁倾角，d0为地磁偏角，be是地磁正常场，bi是目标异常场，a＝cos(i0)cos(d0)、b＝cos(i0)sin(d0)、c＝sin(i0)，la为目标的最大几何尺度，r为观测点与目标的距离，l为每一坐标轴上两个远端传感器的距离，be为be的模，bi为bi的模。
所述根据所述传感器阵列，基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系，包括：
令(xa,ya,za)为目标位置，(xi,yi,zi)是传感器i位置，得：
令
ti＝be+ω·upim
ti＝be+ω(qimx+simy+himz)
其中，bi是目标异常场，bix、biy、biz是bi的直角坐标分量，mx、my、mz是目标磁矩矢量m的直角坐标分量，
qi＝a·fi11+b·fi21+c·fi31,si＝a·fi12+b·fi22+c·fi32,hi＝a·fi13+b·fi23+c·fi33。
所述根据所述地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系，采用矩阵变换，把求解参数消元缩减到三个，包括：
在七个传感器中，选出三组传感器，每组由两个磁传感器i和j组成，i、j＝0,1,2......6。每组传感器测量值的差值δtij＝ti-tj，且这三个差值线性无关；
将磁矩分量用磁性目标的位置坐标表示：
地磁总场测量值ti的梯度gi：
原点o的地磁总场梯度值g0为：
其中，
所述建立适应度函数，用粒子群算法计算上述适应度函数的最小值，求解目标位置，包括：
当l＜＜r的远场情况下，各传感器同步测量得到δt12、δt34、δt56，可得原点地磁总场梯度测量值
显然：
建立粒子群算法求解必需的适应度函数f的表达式为：
其中，l为每一坐标轴上两个远端传感器的距离，r为观测点与目标的距离，
本发明的有益效果在于：
分辨率相对较高，可以对弱磁目标定位，定位距离远、范围大。通过专门的算法设计，可以快速求解定位方程测量时只要避开光泵磁力仪死区方向，无需实时测量和补偿每个传感器姿态，简便快速。
附图说明
图1是地磁总场梯度测量阵列结构图。
图2是传感器测量差值δt20与目标位置xa的关系曲线。
图3是传感器测量差值δt40与目标位置xa的关系曲线。
图4是传感器测量差值δt60与目标位置xa的关系曲线。
图5是运动目标的磁法定位结果。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明的目的在于对磁性目标进行准确快速简便的追踪和定位。首先构建由七个标量磁传感器(光泵磁力仪)组成的传感器阵列，然后基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系。采用矩阵变换，把求解参数消元缩减到三个。建立适应度函数，采用粒子群算法对磁性目标定位。
本发明提供了一种对磁性目标定位的方法，采用多个标量磁传感器构成空间阵列，测量地磁总场梯度，通过专门设计的算法对磁性目标定位。该算法通过对地磁总场梯度矩阵进行变换，分离变量并消除目标的磁矩矢量，用目标空间位置坐标表示目标磁矩矢量三分量，将磁场梯度方程中六个未知数消元减少到三个。建立关于磁场梯度的适应度函数，采用粒子群算法对目标位置求解。本发明所提出的对磁性目标定位的方法，测量精度高，探测距离远，定位快速，操作简便可靠。为磁法目标定位提供了一种新的方法。为地下和水下考古、管线检测、能源矿藏勘测、扫雷反潜等磁性目标探测定位及识别提供参考，有一定应用意义。
本方法用七个标量磁传感器构建了地磁总场梯度测量阵列，设计了目标定位算法，提出了一种新的求解方法，可对磁性目标准确快速定位。其特征是：
(1)按照附图1所示，用七个标量磁传感器构成阵列。阵列中心为坐标原点o，七个传感器位于原点及三个坐标轴上，对称排布。每一坐标轴上两个远端传感器相距为l。阵列中各传感器的坐标为t0(0,0,0)，
(2)建立定位算法方程组。其过程是：基于上述阵列，建立传感器测量值与磁性目标位置矢量和磁矩矢量关系，如式(3)所示。采用矩阵变换，以目标位置坐标表示磁矩矢量的三分量，如(6)式所示。把(3)式中关于目标磁矩和位置的六个未知数消元缩减为三个，得到可求解的总场梯度方程组，如式(12)所示。建立关于地磁总场梯度的适应度函数，如式(13)所示。(3)应用粒子群算法求解地磁总场梯度测量阵列方程组，关键环节是应用粒子群算法计算上述适应度函数的最小值，求解目标位置实现目标定位。2、根据权利要求1所述的对磁性目标定位的方法，阵列中每一坐标轴上远端的两个传感器间距l存在最小值限制。其特征是：l存在最小值，由预设的目标磁矩、定位范围和定位精度代入适应度函数(13)式，求解l最小值。且存在另一必要条件，即采用的传感器工作时存在微量电磁辐射，为防止相互干扰，任意两个传感器之间的距离不小于0.5m。
3、根据权利要求1所述的对磁性目标定位的方法，其特征是：采用本发明选择的传感器构建阵列和使用阵列时，仅需避开光泵磁传感器死区方向，无需实时测量和补偿每个磁传感器姿态。因此该目标定位方法简便快速。
在地下及水下矿藏勘测、管线监测、考古、沉船勘测、扫雷反潜等方面磁法探测具有重要应用意义。因为上述目标使本区域地磁场产生异常，因而用地磁匹配技术可以实现目标定位和识别。在石油能源勘测领域也经常采用地磁探测技术作为辅助勘测手段。在卫星、飞机、舰船、潜艇、车辆等导航应用中，地磁辅助导航提供了一种非常重要的方法和途径。在火山、地震、海啸等自然灾害发生前后，磁暴、地磁异常现象经常伴随发生，这对灾害中心定位和预报提供了可能的技术途径。因此磁法定位技术具有广阔应用领域。本发明提出的一种采用地磁总场梯度对磁性目标定位的方法可应用于上述技术领域。
地磁场是地球的一个天然的物理场，它有各种不同的起源，由不同变化规律的磁场成分叠加而成。按照场源位置划分，地磁场可以分为内源场和外源场。如果考虑地磁场随时间的变化特征，随时间变化较快的地磁场成为地球的变化磁场，随时间变化较慢或者基本不变的地磁场成为地球的稳定磁场。
一个带有磁性的目标产生的磁场会导致空间地磁场分布的变化，从而能够产生磁异常。当观测点与目标距离大于目标尺度2-3倍及以上时，磁性目标在观测点产生的磁场一般被看作磁偶极子远场。
对磁性目标定位时，可采用能够测量地磁三分量的矢量磁传感器或者能够测量地磁总场的标量磁传感器。矢量磁传感器安装和使用相对复杂，初始姿态要严格校正，工作姿态要实时测量。矢量磁传感器的角度误差为0.05°时，测量的地磁误差大概为50nt。因此，实时补偿矢量磁传感器姿态变化对测量的影响是必需的，也是有一定难度的。一般情况下，市场上常见的矢量磁传感器(如普通磁通门磁力仪)分辨率相对偏低(相对于标量光泵磁力仪)，因此基于矢量测量的目标定位距离不能太长，目标磁矩不能太小。
本发明所涉及的是一种采用地磁总场梯度阵列对磁性目标定位的方法。采用的标量磁传感器(光泵磁力仪)分辨率相对较高，可以对弱磁目标定位，定位距离远、范围大。通过专门的算法设计，可以快速求解定位方程。测量时只要避开光泵磁力仪死区方向，无需实时测量和补偿每个传感器姿态，简便快速。这在能源矿藏勘测、地下及水下管线监测维护、考古、坠机沉船搜救、扫雷反潜等方面有一定的应用意义。
本发明的目的在于对磁性目标进行准确快速简便的追踪和定位。首先构建由七个标量磁传感器(光泵磁力仪)组成的传感器阵列，然后基于磁偶极子远场模型，建立地磁总场梯度与磁性目标位置坐标和磁矩矢量的关系。采用矩阵变换，把求解参数消元缩减到三个。建立适应度函数，采用粒子群算法对磁性目标定位。
本发明的技术方案是通过以下步骤实现的：
步骤一：用七个标量磁传感器ti(i＝0,1,2......6)构建如附图1所示的磁探阵列。阵列坐标系的原点为o，x轴正向指向地理北极，y轴正向指东，z轴正向向下，为右手螺旋系。阵列中各磁传感器的直角坐标为t0(0,0,0)，
步骤二：各台磁传感器同步测量，每个测量值ti都是地磁正常场be和目标异常场bi相叠加的总场模，即ti＝|be+bi|。令be的模为be，be是随时间变化的，这种变化在局部空间内，可认为是同步等幅的。令u为be的方向向量，i0和d0为地磁倾角和偏角，再令a＝cos(i0)cos(d0)，b＝cos(i0)sin(d0)，c＝sin(i0)，
得u＝[cos(i0)cos(d0)cos(i0)sin(d0)sin(i0)]＝[abc]。bix、biy、biz是bi的直角坐标分量，bi为bi的模。目标的最大几何尺度为la，观测点与目标的距离为r。远场条件可以表述为3la≤r、l＜＜r、bi＜＜be等。因此
ti＝|be+bi|≈be+u·bi(1)
以磁偶极子远场模型表示磁性目标在传感器i处产生的磁场bi：
式(2)中：μ0＝4π×10-7h/m为真空中的磁导率。令(xa,ya,za)为目标位置，(xi,yi,zi)是传感器i位置，
ti＝be+ω·upim(3)
qi＝a·fi11+b·fi21+c·fi31,si＝a·fi12+b·fi22+c·fi32,hi＝a·fi13+b·fi23+c·fi33。
式(3)可改写为：ti＝be+ω(qimx+simy+himz)(4)
步骤三：在七个传感器中，选出三组传感器，每组由两个磁传感器i和j组成，i、j＝0,1,2......6。每组传感器测量值的差值δtij＝ti-tj。选择原则是保证δtij、δt12、δt34、δt56这六个差值线性无关。根据(3)式，举例选择如下三组差值：
可见，δtij中消除了地磁正常场be，be是时间的函数，所以δtij不受地磁场随时间变化的影响。变换(5)式，把磁矩分量用磁性目标的位置坐标表示为：
在传感器i处，地磁总场测量值ti的梯度gi按(7)式定义，其中gix、giy、giz是gi的直角坐标分量，i＝0,1,2......6。
将(4)式代入(7)式：
根据(9)式，原点o的地磁总场梯度值g0为：
步骤四：当l＜＜r的远场情况下，各传感器同步测量得到δt12、δt34、δt56。可得原点地磁总场梯度测量值
联立(10)和(11)两式，显然
矢量方程(12)式是以目标位置(xa,ya,za)为未知数的高阶三元方程组，含有三个标量方程，约束条件完备，在三维空间内有定解，解的个数有限，根据实际情况可消除部分伪解。建立粒子群算法求解必需的适应度函数f的表达式为：
采用粒子群算法求解适应度函数f的最小值，实现对目标定位。进而代入(6)式求解目标磁矩，可对目标进行初步识别。
图1地磁总场梯度测量阵列结构图
图2传感器测量差值δt20与目标位置xa的关系曲线
图3传感器测量差值δt40与目标位置xa的关系曲线
图4传感器测量差值δt60与目标位置xa的关系曲线
图5运动目标的磁法定位结果
在哈尔滨市松花江沿岸某郊区开展目标定位实验，测量环境磁场数据，以实验过程对本发明实施方式举例说明。
用cs-l光泵磁传感器构建如图1所示的地磁总场梯度测量阵列，传感器分辨率0.6pt，最高采样率10hz，测量范围15000nt-105000nt。阵列参数l＝1m，阵列所在坐标系x轴正向指向地理北极，y轴正向指东，z轴正向向下。当地地磁倾角1.10rad，地磁偏角-0.18rad。磁目标(永久磁铁)磁矩大小为152am2，磁矩倾角1.13rad，磁矩偏角为-0.18rad，目标几何尺度143mm×120mm×40mm。目标运动过程中磁矩方向保持不变。目标初始位置(3.4m，4.21m，0m)，沿平行于x轴方向做直线运动，终止位置(-2.1m，4.21m，0m)。
目标在运动过程中，阵列各传感器同步测量地磁总场，传感器采样率选择5hz，计算原点传感器分别与三个坐标轴负半轴上传感器的测量值之差δt20，δt40，δt60。三个差值随目标位置xa的变化分别如图2、图3、图4所示。其中蓝色虚线表示磁法理论差值，经由公式(5)计算得到。红色圆圈点表示经传感器实际测量得到的差值。可见理论曲线与实际测量点基本吻合。
当目标运动至某一位置(xa,ya,za)，根据传感器同步输出数据计算出δt20、δt40、δt60、δt12、δt34、δt56。分别代入式(10)和(11)，得到g0x、g0y、g0z的表达式以及
图5是该方法对运动目标定位的实验结果。蓝色圆圈点是目标实际位置，红色圆圈点是磁法定位位置。从起点(3.4m，4.21m，0m)至终点(-2.1m，4.21m，0m)直线运动过程中，目标位置xa从2.4m至-2.1m，每间隔0.5m进行一次磁法定位，共计10次。定位结果表明：该方法对大小为143mm×120mm×40mm、磁矩为152am2的磁性运动目标，在4.21m的垂向距离上对目标定位的水平圆概率偏差cep等于0.287m，定位方法可行。