基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法技术

本发明专利技术公开了一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，该方法通过立方体磁探仪求取磁总场梯度的模量以及磁张量模量梯度的模量，进一步联立求取磁目标的方位和磁矩大小的解析解。本发明专利技术提出的磁目标单点定位方法，基于磁总场梯度和磁张量模量梯度的旋转不变性，能够有效规避磁目标和磁探仪相对旋转和相对运动带来的定位不准确问题，对特定场景下的磁异常目标定位应用具有重要的实际意义。

Single-point localization of magnetic targets based on total magnetic field gradient and magnetic tensor modulus gradient

The invention discloses a single point positioning method for magnetic target based on magnetic total field gradient and magnetic tensor modulus gradient. The method calculates the modulus of magnetic total field gradient and magnetic tensor modulus gradient by cubic magnetometer, and further solves the analytical solution of magnetic target orientation and magnetic moment. Based on the rotation invariance of the magnetic total field gradient and the magnetic tensor modulus gradient, the proposed single point positioning method can effectively avoid the inaccurate positioning problems caused by the relative rotation and relative motion of the magnetic target and the magnetic detector, and has important practical significance for the application of magnetic anomaly target positioning in a specific scene.

【技术实现步骤摘要】
基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法
本专利技术涉及磁偶极子定位技术，具体涉及一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法。
技术介绍
一个远距离磁性目标通常可以视为磁偶极子。基于磁偶极子的磁异常方程，对目标进行定位和磁矩参数估计至少需要求解6个方程，因此，通常采用多个磁传感器对目标进行测量进行参数估计。基于目标异常场总场或分量的定位方法，需要求解非线性方程，难以获得参数的解析解。基于多测量点数据的非线性方法和粒子算法被研究，存在定位算法复杂，计算量大和实时性差等不足，而且对测量平台或目标运动状态提出了限制。同时，磁异常信号相对地磁场是一个非常弱的信号，无论是标量测量或是向量测量，从地磁场分离出磁异常信号通常是困难的。在磁场方程的反演研究中，一直在试图利用磁场差分和模量对目标进行检测、定位和识别(DLC：detection,localizationandClassification)。前者可以规避航空磁探测的2个关键问题：探测设备与环境噪声的排除和补偿问题。后者利用模量的旋转不变性降低对运动平台的敏感性。基于梯度张量模量旋转不变量的定位方法得到了研究，并表现出比较好的性能。最早开始该方向研究的RoyWiegert提出STAR方法(Wiegert,R.,andJ.Oeschger,”Generalizedmagneticgradientcontractionbasedmethodfordetection,localizationanddiscriminationofunderwaterminesandunexplodedordnance,”OceansIEEE,vol.2,pp.1325-1332,Sept.2005.)。STAR方法基于梯度张量模量的旋转不变性解决了测量平台姿态不平稳以及磁目标和测量平台间存在相对运动等问题，但其系统固有误差目前还未有解析方法进行消除。
技术实现思路
本专利技术的目的在于提供一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，从原理上消除系统定位误差，提高磁目标定位精度。实现本专利技术目的的技术方案为：一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，运用三轴向量磁强计为传感器，具体步骤如下：步骤1，选取磁探仪正立方体构型，获取各磁强计三轴磁分量；步骤2，利用磁强计三轴磁分量求取各磁强计测量磁场的模量；步骤3，针对磁探仪正立方体构型，求取磁探仪中心处磁场模量的梯度；步骤4，针对磁探仪正立方体构型，求取磁探仪各面中心处及磁探仪中心处磁场的梯度张量；步骤5，利用步骤4求得的磁探仪各面中心处磁场的梯度张量，求取梯度张量的模量；步骤6，利用步骤5求得的张量模量，对磁探仪中心处张量模量的梯度进行工程化近似获取；步骤7，对步骤3和步骤6求取的磁探仪中心处梯度和张量模量梯度分别进行归一化处理；步骤8，求取磁目标位置以及磁目标磁矩大小，完成磁目标定位。与现有技术相比，本专利技术的有益效果在于：本专利技术结合总场梯度模量和磁张量模量梯度的模量这两种旋转不变量，解决了测量平台姿态不稳及测量平台和磁目标之间存在相对运动问题，同时基于这两种旋转不变量的归一化处理，给出了对磁目标相对位置以及磁目标磁矩大小的解析解，有效解决了以往方法存在的系统固有误差问题，提高磁目标定位精度，具有实际意义。附图说明图1为磁目标单点定位方法流程图。图2为磁偶极子模型示意图。图3为正立方体阵列磁探仪结构示意图。图4为磁张量模量梯度及总场梯度计算示意图。图5为短基线，无噪声情况下三种方法定位结果对比图。图6为长基线，有噪声情况下三种方法定位结果对比图。具体实施方式结合图1，一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位解析方法，运用三轴向量磁强计构建立方体磁探仪进行磁目标定位，包括以下步骤：步骤1，选取磁探仪正立方体构型，获取各磁强计三轴磁分量，分别用Bnx，Bny，Bnz表示。步骤2，利用磁强计三轴磁分量Bnx，Bny，Bnz，求取各磁强计测量磁场的模量：步骤3，针对磁探仪正立方体构型，求取磁探仪中心处磁场模量的梯度，模量梯度的表达式为：G≈iGx+jGy+kGz；利用各磁强计模量对模量梯度进行工程化近似获取，模量梯度的三轴分量可以表示为：式中，分别为正x方向、正y方向、正z方向的磁场模量，的值对应的分别是正立方体正x面、正y面、正z面四个顶点磁强计所测磁场模量的平均值；分别为负x方向、负y方向、负z方向的磁场模量，的值对应的分别是正立方体负x面、负y面、负z面四个顶点磁强计所测磁场模量的平均值。步骤4，针对磁探仪正立方体构型，求取磁探仪各面中心处及磁探仪中心处磁场的梯度张量，梯度张量的表达式为：其中Gxy＝Gyx,Gxz＝Gzx,Gyz＝Gzy,Gxx+Gyy+Gzz＝0。依据上述梯度张量特性，利用磁强计测得三轴磁分量，对梯度张量进行工程化近似获取，Gab为磁场a轴分量在b轴上的梯度，即a,b轴为x，y，z中任意轴。步骤5，利用步骤4求得的磁探仪各面中心处磁场的梯度张量，求取梯度张量的模量：步骤6，利用步骤5求得的张量模量，对磁探仪中心处张量模量的梯度进行工程化近似获取，张量模量梯度的三轴分量可以表示为：式中，分别为正x面、正y面、正z面中心处的梯度张量，分别为负x面、负y面、负z面中心处的梯度张量。步骤7，对步骤3和步骤6求取的磁探仪中心处梯度G和张量模量梯度D分别进行归一化处理，即G0＝G/|G|和D0＝D/|D|。步骤8，利用以下方程组联立求取磁目标位置以及磁目标磁矩大小，完成磁目标定位：η＝(|T|/|G|)2(5)M＝4πr4|T|/3μf(8)式中各向量均为在磁探仪本体坐标系所求值，r0为目标方位单位向量，m0为磁目标磁矩单位向量，为磁目标磁矩单位向量和目标方位单位向量的夹角，r为目标距磁探仪中心处距离大小，M为目标磁矩大小，μ为真空磁导率；式中各向量均以磁探仪本体坐标系为参考坐标系。在磁探问题中，磁目标距磁探仪距离为磁目标本体尺寸3倍以上，即磁目标可近似等效为磁偶极子。本专利技术对磁目标与磁探仪之间的相对运动不敏感，可以实现对磁目标的单点实时定位。下面结合附图对本专利技术做进一步的说明，为了更好理解本专利技术的技术方案，现将其原理及计算公式结合具体实施例推导详细叙述如下：如图2所示的磁偶极子参考坐标系，磁偶极子M在原点O即测量点处产生的磁场表示为式中μ为真空磁导率，r为目标方位向量，r为目标距磁探仪中心处距离大小，M为目标磁矩向量；B的模量表示为式中，M为目标磁矩大小；根据(11)式可得到|B|的梯度，也称为总场梯度，表达式为：式中i,j,k为直角坐标系单位向量。根据(12)式可以得到梯度的模|G|为根据(12)式和(13)式，可得模量|B|梯度的归一化表达式(TFG定位)由(10)式可得目标磁场B的张量表示为I为3×3的单位矩阵；磁场张量的模量可表示为式中根据(16)式，可得到张量模量的梯度由上式可得到D的模量|D|为根据(16)式和(17)式，可得(TMG定位)由(13)式和(16)式可得求解(20)式可得根据(21)式，容易验证具有唯一解将p2代入(16)式及(18)式，可以得到距离计算公式将上式代入(16)式，可以得到磁距计算公式由(14)式和(19)式，可得到r0的表达式(T本文档来自技高网...

【技术保护点】
1.一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，运用三轴向量磁强计为传感器，其特征在于，具体步骤如下：步骤1，选取磁探仪正立方体构型，获取各磁强计三轴磁分量；步骤2，利用磁强计三轴磁分量求取各磁强计测量磁场的模量；步骤3，针对磁探仪正立方体构型，求取磁探仪中心处磁场模量的梯度；步骤4，针对磁探仪正立方体构型，求取磁探仪各面中心处及磁探仪中心处磁场的梯度张量；步骤5，利用步骤4求得的磁探仪各面中心处磁场的梯度张量，求取梯度张量的模量；步骤6，利用步骤5求得的张量模量，对磁探仪中心处张量模量的梯度进行工程化近似获取；步骤7，对步骤3和步骤6求取的磁探仪中心处梯度和张量模量梯度分别进行归一化处理；步骤8，求取磁目标位置以及磁目标磁矩大小，完成磁目标定位。

【技术特征摘要】
1.一种基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，运用三轴向量磁强计为传感器，其特征在于，具体步骤如下：步骤1，选取磁探仪正立方体构型，获取各磁强计三轴磁分量；步骤2，利用磁强计三轴磁分量求取各磁强计测量磁场的模量；步骤3，针对磁探仪正立方体构型，求取磁探仪中心处磁场模量的梯度；步骤4，针对磁探仪正立方体构型，求取磁探仪各面中心处及磁探仪中心处磁场的梯度张量；步骤5，利用步骤4求得的磁探仪各面中心处磁场的梯度张量，求取梯度张量的模量；步骤6，利用步骤5求得的张量模量，对磁探仪中心处张量模量的梯度进行工程化近似获取；步骤7，对步骤3和步骤6求取的磁探仪中心处梯度和张量模量梯度分别进行归一化处理；步骤8，求取磁目标位置以及磁目标磁矩大小，完成磁目标定位。2.根据权利要求1所述的基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，其特征在于，步骤2利用磁强计三轴磁分量Bnx，Bny，Bnz，求取各磁强计测量磁场的模量：3.根据权利要求2所述的基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，其特征在于，步骤3具体为：针对磁探仪正立方体构型，求取磁探仪中心处磁场模量的梯度，模量梯度的表达式为：G≈iGx+jGy+kGz；利用各磁强计模量对模量梯度进行工程化近似获取，模量梯度的三轴分量表示为：式中，分别为正x方向、正y方向、正z方向的磁场模量，的值对应的分别是正立方体正x面、正y面、正z面四个顶点磁强计所测磁场模量的平均值；分别为负x方向、负y方向、负z方向的磁场模量，的值对应的分别是正立方体负x面、负y面、负z面四个顶点磁强计所测磁场模量的平均值。4.根据权利要求3所述的基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法，其特征...

【专利技术属性】
技术研发人员：庄志洪，金煌煌，王宏波，王容川，付梦印，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32