基于两点磁梯度全张量的线性定位方法

摘要

本发明公开了一种基于两点磁梯度全张量的线性定位方法，包括如下步骤：S1，以空间中任意一点为原点建立空间直角坐标系，确定第一观测点与第二观测点的坐标，得到待观测目标点的定位公式；S2，分别测量两个观测点的磁梯度全张量；S3，将步骤S2中测得的数据代入步骤S1中，计算得到待定位目标点的定位数据。本发明提供的基于两点磁梯度全张量的线性定位方法，可以根据两点磁梯度全张量以及两点的相对位置实现对目标点的定位，期间不需要测量地磁场值，减小了地磁场噪声所产生的定位误差，并且采用线性方法完成解算，求解过程简便，可以直接得出解析解。

说明书

技术领域

本发明涉及磁探测领域，尤其是涉及一种基于两点磁梯度全张量的线性定位方法。

背景技术

磁探测是研究磁性目标相关信息的重要技术手段，特别是在识别定位深海或地下磁性掩埋物方面，磁异常探测技术具有更好的主动性和隐蔽性。采用单个观测点实现的单点磁梯度全张量磁性目标定位方法存在受地磁影响大的问题。现有技术中也存在一些采用两个或多个观测点实现定位的方法，现有的两点或多点的定位方法虽然减小了地磁场的影响，却基本都面临着非线性方程组求解的问题，求解过程复杂、无法得到解析最优解。

发明内容

为解决上述背景技术中提出的问题，本发明的目的在于提供一种基于两点磁梯度全张量的线性定位方法。

为实现上述目的，本发明采取的技术方案为：

一种基于两点磁梯度全张量的线性定位方法，包括如下步骤：

S1，以空间中任意一点为原点建立空间直角坐标系，确定第一观测点与第二观测点的坐标，得到待定位目标点的定位公式：

其中，

式中，

S2，分别测量两个观测点的磁梯度全张量G

S3，将步骤S2中测得的数据代入步骤S1中，计算得到待定位目标点的定位数据。

在一些实施例中，步骤S1中，得到待定位目标点的定位公式的具体步骤如下：

对于任意点(x,y,z)的磁场矢量

在磁偶极子空间任意一点的磁场

其中，μ

设在位置矢量

根据公式(4)和(5)可得

在公式(6)中

将公式(7)代入公式(5)可得

在公式(8)中

因此

根据公式(8)和(10)可得：

将目标点到第一观测点的位置矢量表示为

由公式(11)可得

由公式(12)可得

根据公式(12)、(13)、(14)、(15)即可得到式(2)：

进一步便可得到式(1)：

在一些实施例中，步骤S2中，采用四个对称设置的磁通门传感器构成十字形磁梯度全张量测量系统，利用该系统测量磁场数据，并根据测得的数据分别计算得出第一观测点和第二观测点的磁梯度全张量G

在一些实施例中，步骤S2中，具体根据如下步骤测量两个观测点的磁梯度全张量：

根据麦克斯韦方程，有

由公式(16)和(17)可得

因此公式(3)可表示为：

将待测点放置在十字形磁梯度全张量测量系统的几何中心处，则待测点的磁梯度全张量可表示为：

其中，d表示每两个相对磁通门的基线距离，B

分别将第一观测点和第二观测点放置在十字形磁梯度全张量测量系统的几何中心处，并利用公式(20)即可分别测得第一观测点和第二观测点的磁梯度全张量G

与现有技术相比，本发明的有益效果是：

本发明提供的基于两点磁梯度全张量的线性定位方法，可以根据两点磁梯度全张量以及两点的相对位置实现对目标点的定位，期间不需要测量地磁场值，减小了地磁场噪声所产生的定位误差，并且采用线性方法完成解算，求解过程简便，可以直接得出解析解。

附图说明

图1为利用本发明提供的基于两点磁梯度全张量的线性定位方法时建立的空间直角坐标系的示意图；

图2为十字形磁梯度全张量测量系统模型示意图。

具体实施方式

为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合附图和具体实施方式，进一步阐述本发明是如何实施的。

本发明提供了一种基于两点磁梯度全张量的线性定位方法，包括如下步骤：

S1，如图1所示，以空间中任意一点为原点建立空间直角坐标系，确定第一观测点与第二观测点的坐标，得到待定位目标点的定位公式：

其中，

式中，

S2，利用图2所示的十字形磁梯度全张量测量系统分别测量两个观测点的磁梯度全张量G

S3，将步骤S2中测得的数据代入步骤S1中，计算得到待定位目标点的定位数据。

进一步地，步骤S1中，得到待定位目标点的定位公式的具体步骤如下：

对于任意点(x,y,z)的磁场矢量

磁性目标与探测点的距离大于2.5倍的目标长度时，磁性目标可以视为磁偶极子，在磁偶极子空间任意一点的磁场

其中，μ

设在位置矢量

根据公式(4)和(5)可得

在公式(6)中

将公式(7)代入公式(5)可得

在公式(8)中

因此

根据公式(8)和(10)可得：

将目标点到第一观测点的位置矢量表示为

由公式(11)可得

由公式(12)可得

根据公式(12)、(13)、(14)、(15)即可得到式(2)：

进一步便可得到式(1)：

进一步参照图2所示，步骤S2中，采用四个对称设置的磁通门传感器构成十字形磁梯度全张量测量系统，利用该系统测量磁场数据，并根据测得的数据分别计算得出第一观测点和第二观测点的磁梯度全张量G

具体地，步骤S2中，具体根据如下步骤测量两个观测点的磁梯度全张量：

根据麦克斯韦方程，有

由公式(16)和(17)可得

因此公式(3)可表示为：

将待测点放置在十字形磁梯度全张量测量系统的几何中心处，采用差分的方法即可获取张量各分量，则待测点的磁梯度全张量可表示为：

其中，d表示每两个相对磁通门的基线距离，B

分别将第一观测点和第二观测点放置在十字形磁梯度全张量测量系统的几何中心处，并利用公式(20)即可分别测得第一观测点和第二观测点的磁梯度全张量G

综上，本发明提供的基于两点磁梯度全张量的线性定位方法，可以根据两点磁梯度全张量以及两点的相对位置实现对目标点的定位，期间不需要测量地磁场值，减小了地磁场噪声所产生的定位误差，并且采用线性方法完成解算，求解过程简便，可以直接得出解析解。

最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围中。