CN106546235B - 一种基于载体补偿的磁性目标定位方法 - Google Patents

一种基于载体补偿的磁性目标定位方法 Download PDF

Info

Publication number
CN106546235B
CN106546235B CN201610942236.6A CN201610942236A CN106546235B CN 106546235 B CN106546235 B CN 106546235B CN 201610942236 A CN201610942236 A CN 201610942236A CN 106546235 B CN106546235 B CN 106546235B
Authority
CN
China
Prior art keywords
magnetic
magnetometer
target
measurement
array
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610942236.6A
Other languages
English (en)
Other versions
CN106546235A (zh
Inventor
康崇
郑权
樊黎明
周健
张晓峻
王明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201610942236.6A priority Critical patent/CN106546235B/zh
Publication of CN106546235A publication Critical patent/CN106546235A/zh
Application granted granted Critical
Publication of CN106546235B publication Critical patent/CN106546235B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于磁场测量领域,具体涉及一种基于载体补偿的磁性目标定位方法。本发明包括:在水面或水下利用4台磁传感器构建磁力仪阵列;测量在无磁性目标的条件下,记录每个磁力仪测量的随时间变换的磁场强度值,把磁场强度值和其中一个磁力仪记录的值进行自回归分析,得出测量值之间的线性关系;利用阵列中的传感器测量磁性目标产生的磁异常;构建磁偶极子模型,获得磁性目标在测量点处产生的磁场等。本发明所提出的基于载体补偿的磁性目标定位方法可以消除载体对磁力仪的影响,提高对磁性目标定位的精度,同时标量传感器测量的地磁信息是一种旋转不变量,从而传感器阵列的布放与方位无关,因此该探测方法实施简单,定位精度高,定位距离远。

Description

一种基于载体补偿的磁性目标定位方法
技术领域
本发明属于磁场测量领域,具体涉及一种基于载体补偿的磁性目标定位方法。
背景技术
地磁场是反映宇宙演变、地球演变、地质构造演变及地震活动等过程的重要物理量之一。地磁场研究成果在航海、航空、航天、能源、矿产、安全、考古等领域中有着广泛而重要的应用[2,3]
发明所涉及的是一种基于载体补偿的磁性目标定位方法。具体地说是构建由四个标量磁力仪组成的传感器阵列,对传感器测得的数据进行自回归分析,获得不同磁力仪测量值之间的线性关系,进而消除载体对磁力仪的影响,最后测得地磁总场数据,依据磁性目标磁偶极子的远场理论,将磁性目标的空间位置信息转化到对应的地磁总场信息中,再通过改进的粒子群算法计算出磁性目标的位置坐标,实现磁性目标的定位。由于消除了载体对磁力仪的影响,可以提高测量的精度。在水下能源矿藏勘测、水下管线维护监测、水下考古、沉船勘测、扫雷反潜等方面有重要应用。
地磁场是地球的一个天然的物理场,它有各种不同的起源,由不同变化规律的磁场成分叠加而成。按照场源位置划分,地磁场可以分为内源场和外源场。如果考虑地磁场随时间的变化特征,将随时间变化较快的地磁场成为地球的变化磁场,随时间变化较慢或者基本不变的地磁场成为地球的稳定磁场。
在各种应用领域中,精确确定目标物的位置是一项首要任务,是进行后续工作的前提。如军事上需要进行的沉没船只的货物抢救、排雷、海滩救援作业、港口船舶监测、反潜应用等,都需要对水下目标物进行准确而快速的定位。我国黄海平均海深50米,东海多为200米的大陆架,在这种环境下,海况和目标噪声是决定声呐探测距离的最大因素。而基于磁场探测则不用考虑这些因素。由于磁性目标的存在,其产生的感应磁场会导致空间地磁场分布的变化,从而在该空间中产生磁异常。因此磁测技术是非常有效的方法,人们可以通过对磁异常的反演,获得该目标物体的一些信息(如,几何参数,位置参数等)。通过消除载体的影响可以提高探测的精度,提高信息的可用度。
对磁性目标进行定位时,一般需要能够测量地磁分离的矢量传感器或者能够测量地磁总场的标量传感器中的一种。在应用矢量传感器进行测量过程中,传感器的安装很复杂,安装时姿态方位一定要严格校正。当传感器的角度误差为0.05°时,测量的地磁误差大概为50nT左右。因此在运动过程中仍要实时补偿姿态和方位变化的影响,校正姿态方位还要使用其他高精度定位***。同时由于地磁场随时间变化的影响,基于矢量传感器的方法的测量距离不能太长。
相比于矢量传感而言,探测地磁总场的标量传感器光泵磁力仪具有高可靠高精度的特点,测量的地磁总场值不会因为传感器的旋转而产生变化,同时最高分辨率可达fT量级。因此,在fT量级下的传感器阵列的极限探测距离超过10km。由于测量地磁总场,光泵磁力仪安装使用不需要姿态方位校准,非常方便。
发明内容
本发明的目的在于提供一种基于载体补偿的磁性目标定位方法。
本发明的目的是这样实现的:
(1)在水面或水下利用4台磁传感器构建磁力仪阵列;
(2)测量在无磁性目标的条件下,记录每个磁力仪测量的随时间变换的磁场强度值,把磁场强度值和其中一个磁力仪记录的值进行自回归分析,得出测量值之间的线性关系:
Nm1=α12Nm212
其中:Nm1是第一个磁力仪的输出值,Nm2是第二个磁力仪的输出值;其中k1,k2是感生磁场与地磁场的比例系数,即Hid1=k1He,Hid2=k2He,Hpd1和Hpd2分别恒定磁场的值;
(3)利用阵列中的传感器测量磁性目标产生的磁异常ΔB,磁力仪测量值为:
Hm1=He+Hpd1+Hid1+ΔΒ1
=He+Hpd1+k1He+ΔΒ1
=Nm1+ΔΒ1
Hm2=He+Hpd2+Hid2+ΔΒ2
=He+Hpd2+k2He+ΔΒ2
=Nm2+ΔΒ2
其中,Hm1是有目标时第一个磁力仪的输出值,Hm2是有目标时第二个磁力仪的输出值,得到:
ΔB112ΔB2=Hm112Hm212
(4)构建磁偶极子模型,获得磁性目标在测量点处产生的磁场Ba的表示形式:
其中:μ0为真空中的磁导率,Ba为地磁场矢量,m是磁偶极子磁矩;(0,0,0)表示磁偶极子的位置坐标,(x,y,z)表示测量点处的位置坐标,
(5)建立磁异常ΔB和磁性目标位置信息(x,y,z)的关系;
其中:G=[cos(I0)cos(D0) cos(I0)sin(D0) sin(I0)]
其中I0是地磁倾角,D0是地磁偏角;
得到:
ΔBiijΔBj=G(KiijKj)M;
(6)用矩阵变换分离物体的位置和磁矩;
其中,MT(MMT)-1GT是Mx,My,Mz,I0,D0的函数,是x,y,zD0,ΔB,I0的函数,由此构建磁力仪阵列,MT(MMT)-1GT对每个磁力仪探测器是一定的:
(7)利用粒子群算法求解磁性目标的位置信息,对目标的追踪和定位;
目标的位置由下列函数的最小值得到;
本发明的有益效果在于:
本发明所提出的基于载体补偿的磁性目标定位方法可以消除载体对磁力仪的影响,提高对磁性目标定位的精度,同时标量传感器测量的地磁信息是一种旋转不变量,从而传感器阵列的布放与方位无关,因此该探测方法实施简单,定位精度高,定位距离远。
附图说明
图1标量传感器阵列示意图;
图2传感器1测量值和其余传感器测量值之间自回归分析;
图3磁性目标的运动轨迹和预测的运动轨迹。
具体实施方式
下面结合附图对本发明做进一步描述。
一种基于载体补偿的磁性目标定位方法,构建如图1所示的标量传感器阵列,由四个标量传感器记录磁场强度的信息,获得了标量传感器之间在测量磁场强度时的函数关系,可以消除目标定位中载体的影响,获得对应的磁异常,从而实现对磁性目标的三维追踪定位。
本发明采用图1所示阵列中多个标量传感器输出值之函数补偿载体的干扰。由四个标量传感器记录磁场强度的信息,利用线性回归方法确定标量传感器之间在测量磁场强度时的函数关系,利用此函数关系整体补偿载体的干扰,具体步骤见发明内容步骤三。
本发明构建目标定位的函数模型。用粒子群算法构建出图1所示阵列的目标定位函数模型,实现对目标的磁场定位。
本发明的目的在于消除载体的影响,对磁性目标进行高精度的追踪和定位。先记录四个磁力仪测量的地磁强度,通过自回归分析这些数据,得到两个传感器测量磁场强度之间的关系,从而消除载体对于磁力仪的影响,再利用四个标量磁力仪构成的传感器阵列获得的地磁总场数据,依据磁性目标磁偶极子的远场理论,将磁性目标的空间位置信息转化到对应的地磁总场信息中,最后通过改进的粒子群算法计算出磁性目标的位置坐标,实现对目标的精确定位。
本发明是通过以下技术方案实现的:
基于载体补偿的磁性目标定位方法,包括以下步骤:
步骤一:在水面或水下利用4台磁传感器构建磁力仪阵列;
步骤二:测量在无磁性目标的条件下,记录每个磁力仪测量的随时间变换的磁场强度值。把这些值和其中一个磁力仪记录的值进行自回归分析,得出测量值之间的线性关系:
Nm1=α12Nm212 (1)
其中:Nm1是第一个磁力仪的输出值,Nm2是第二个磁力仪的输出值。其中k1,k2是感生磁场与地磁场的比例系数,即Hid1=k1He,Hid2=k2He,Hpd1和Hpd2分别恒定磁场的值。
步骤三:利用阵列中的传感器测量磁性目标产生的磁异常ΔB,磁力仪测量值可以表示为:
其中,Hm1是有目标时第一个磁力仪的输出值,Hm2是有目标时第二个磁力仪的输出值。由公式(1),(2),(3)可以得到:
ΔB112ΔB2=Hm112Hm212 (4)
由此,可以消除载体对磁力仪的影响。
步骤四:构建磁偶极子模型,获得磁性目标在测量点处产生的磁场Ba的表示形式:
其中:μ0为真空中的磁导率(μ0=4π10-7),Ba为地磁场矢量,m是磁偶极子磁矩。(0,0,0)表示磁偶极子的位置坐标,(x,y,z)表示测量点处的位置坐标,
步骤五:建立磁异常ΔB和磁性目标位置信息(x,y,z)的关系;
其中:G=[cos(I0)cos(D0) cos(I0)sin(D0) sin(I0)]
其中I0是地磁倾角,D0是地磁偏角。
由公式(4)和(6)可以得到:
ΔBiijΔBj=G(KiijKj)M (7)
步骤六:用矩阵变换分离物体的位置和磁矩;
其中,MT(MMT)-1GT是Mx,My,Mz,I0,D0的函数。是x,y,zD0,ΔB,I0的函数。由此构建磁力仪阵列,MT(MMT)-1GT对每个磁力仪探测器来说是一定的。得到:
步骤七:利用粒子群算法求解磁性目标的位置信息,实现对目标的追踪和定位。
当I0,D0,ΔB时已知时,目标的位置可由下列函数的最小值得到;
本发明是一种基于载体补偿的磁性目标定位方法,还可以包括:
1、阵列几何形状
如附图1,在水面或水下利用四台磁传感器排列成阵列,参照附图1。
2、传感器布设方向
采用光泵磁传感器时,传感器光轴取向与地磁场矢量T0方向夹角为[10,85]或者[95,170],同时传感器阵列中T1所在边对准地理北极方向。
3、确定传感器之间的间距L1、L2的原则:
首先在各种客观条件允许的情况下,相邻传感器之间的距离L越大,阵列的定位精度越高。
距离L不能超过传感器的有效测量范围,这是由构成阵列传感器的分辨率决定的。
距离L不要妨碍阵列的机动航行,这是由装载或拖曳阵列载体的载荷能力及航行速度决定的。
4、磁力仪探测器测量值之间进行自回归分析结果参照附图2。
5、通过获得的磁性目标的位置信息,可以估算出对应的磁矩M,从而可以初步判断目标的尺寸大小。
以图1结构搭建标量磁传感器阵列,阵列中的传感器间距L1=0.8m,L2=0.6m磁传感器采用灵敏度为0.6pT的CS-L光泵磁力仪。传感器阵列沿着平行与地磁北极的方向进行放置。获得传感器之间的线性关系如图2。磁性目标在水平面内移动,在x方向上的移动速度为0m/s,在Y方向上的移动速度为1m/s。采样间隔为1s。磁性目标沿着规划的航迹从(20,-20)点运动到(20,20)点,如图3所示。采用本方法的追踪定位的结果显示:X,Y和Z方向上的相对误差小于6%。磁性目标到传感器1的距离的平均相对误差为0.05%。可见该方案能够对磁性目标进行高精度的追踪和定位。同时计算出的磁性目标磁矩为P=145A·m2

Claims (1)

1.一种基于载体补偿的磁性目标定位方法,其特征在于,包括以下步骤:
(1)在水面或水下利用4台磁传感器构建磁力仪阵列;
(2)测量在无磁性目标的条件下,记录每个磁力仪测量的随时间变换的磁场强度值,把磁场强度值和其中一个磁力仪记录的值进行自回归分析,得出测量值之间的线性关系:
Nm1=α12Nm212
其中:Nm1是第一个磁力仪的输出值,Nm2是第二个磁力仪的输出值;其中k1,k2是感生磁场与地磁场的比例系数,即Hid1=k1He,Hid2=k2He,Hpd1和Hpd2分别恒定磁场的值;
(3)利用阵列中的传感器测量磁性目标产生的磁异常ΔB,磁力仪测量值为:
Hm1=He+Hpd1+Hid1+ΔB1
=He+Hpd1+k1He+ΔB1
=Nm1+ΔB1
Hm2=He+Hpd2+Hid2+ΔB2
=He+Hpd2+k2He+ΔB2
=Nm2+ΔB2
其中,Hm1是有目标时第一个磁力仪的输出值,Hm2是有目标时第二个磁力仪的输出值,得到:
ΔB112ΔB2=Hm112Hm212
(4)构建磁偶极子模型,获得磁性目标在测量点处产生的磁场Ba的表示形式:
其中:μ0为真空中的磁导率,Ba为地磁场矢量,m是磁偶极子磁矩;(0,0,0)表示磁偶极子的位置坐标,(x,y,z)表示测量点处的位置坐标,
(5)建立磁异常ΔB和磁性目标位置信息(x,y,z)的关系;
其中:G=[cos(I0)cos(D0) cos(I0)sin(D0) sin(I0)]
其中I0是地磁倾角,D0是地磁偏角;
得到:
ΔBiijΔBj=G(KiijKj)M;
(6)用矩阵变换分离物体的位置和磁矩;
其中,MT(MMT)-1GT是Mx,My,Mz,I0,D0的函数,是x,y,z D0,ΔB,I0的函数,由此构建磁力仪阵列,MT(MMT)-1GT对每个磁力仪探测器是一定的:
(7)利用粒子群算法求解磁性目标的位置信息,对目标的追踪和定位;
目标的位置由下列函数的最小值得到;
CN201610942236.6A 2016-11-02 2016-11-02 一种基于载体补偿的磁性目标定位方法 Active CN106546235B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610942236.6A CN106546235B (zh) 2016-11-02 2016-11-02 一种基于载体补偿的磁性目标定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610942236.6A CN106546235B (zh) 2016-11-02 2016-11-02 一种基于载体补偿的磁性目标定位方法

Publications (2)

Publication Number Publication Date
CN106546235A CN106546235A (zh) 2017-03-29
CN106546235B true CN106546235B (zh) 2019-07-16

Family

ID=58392333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610942236.6A Active CN106546235B (zh) 2016-11-02 2016-11-02 一种基于载体补偿的磁性目标定位方法

Country Status (1)

Country Link
CN (1) CN106546235B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108759643B (zh) * 2018-03-28 2019-08-09 华中科技大学 一种基于磁传感器阵列坐标系的磁性目标定位方法
CN109725360B (zh) * 2019-01-25 2020-06-09 中北大学 基于磁梯度张量不变量的单点定位方法
CN111273202A (zh) * 2020-02-25 2020-06-12 中国电子科技集团公司第二十九研究所 一种基于阵列的磁传感器补偿方法
CN111399066B (zh) * 2020-04-03 2022-05-17 西北工业大学青岛研究院 一种基于正交基函数处理标量磁异常梯度信号的方法
CN111913227B (zh) * 2020-08-11 2023-03-28 哈尔滨工程大学 一种磁异常信号频率特征的计算方法和***
CN113218393B (zh) * 2021-04-06 2023-06-13 青岛海月辉科技有限公司 基于磁异常总场匹配定位算法的水下目标分布式组网定位方法
CN114440871A (zh) * 2021-12-29 2022-05-06 宜昌测试技术研究所 一种基于自适应互补滤波的九轴磁罗盘数据融合方法
CN114440870B (zh) * 2021-12-29 2023-06-27 宜昌测试技术研究所 一种基于自适应互补滤波的九轴磁罗盘姿态获取方法
CN116975991B (zh) * 2022-04-24 2024-04-26 中国人民解放军海军工程大学 一种基于粒子群算法的磁目标三维寻优反演定位方法
CN115236746A (zh) * 2022-06-02 2022-10-25 吉林大学 一种潜航器搭载式水下多参量磁测***及磁测定位方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1388727A1 (en) * 2002-08-06 2004-02-11 Chan-Gook Park Apparatus and method for calculating azimuth angle
CN103885002A (zh) * 2014-03-06 2014-06-25 中国船舶重工集团公司第七一〇研究所 磁传感器阵列测量中的平行度误差补偿方法和***
CN105091880A (zh) * 2015-07-17 2015-11-25 哈尔滨工程大学 一种基于标量传感器阵列的追踪定位水下远距离磁性目标的方法
CN105783919A (zh) * 2016-03-09 2016-07-20 哈尔滨工程大学 H型标量传感器阵列对磁性目标的追踪定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1388727A1 (en) * 2002-08-06 2004-02-11 Chan-Gook Park Apparatus and method for calculating azimuth angle
CN103885002A (zh) * 2014-03-06 2014-06-25 中国船舶重工集团公司第七一〇研究所 磁传感器阵列测量中的平行度误差补偿方法和***
CN105091880A (zh) * 2015-07-17 2015-11-25 哈尔滨工程大学 一种基于标量传感器阵列的追踪定位水下远距离磁性目标的方法
CN105783919A (zh) * 2016-03-09 2016-07-20 哈尔滨工程大学 H型标量传感器阵列对磁性目标的追踪定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
运动平台的磁性目标定时定位方法;于振涛等;《哈尔滨工程大学学报》;20150531;第36卷(第5期);第606-610页

Also Published As

Publication number Publication date
CN106546235A (zh) 2017-03-29

Similar Documents

Publication Publication Date Title
CN106546235B (zh) 一种基于载体补偿的磁性目标定位方法
CN105783919B (zh) H型标量传感器阵列对磁性目标的追踪定位方法
CN103926625B (zh) 一种利用地磁总场对水下磁目标高精度远距离的定位方法
CN107272069B (zh) 基于磁异常梯度的磁性目标追踪方法
CN105091880B (zh) 一种基于标量传感器阵列的追踪定位水下远距离磁性目标的方法
Clark Corrigendum to: New methods for interpretation of magnetic vector and gradient tensor data I: eigenvector analysis and the normalised source strength
US6476610B1 (en) Magnetic anomaly sensing system and methods for maneuverable sensing platforms
Storms et al. Magnetic field navigation in an indoor environment
Li et al. Magnetic sensors for navigation applications: an overview
Liu et al. Calibration and compensation of geomagnetic vector measurement system and improvement of magnetic anomaly detection
CN112558006B (zh) 一种用于海底缆线的双三轴电磁探测及路由定位方法
Pang et al. Integrated compensation of magnetometer array magnetic distortion field and improvement of magnetic object localization
US7038458B1 (en) Magnetic anomaly homing system and method using rotationally invariant scalar contractions of magnetic gradient tensors
CN103926627A (zh) 水下载体地磁三分量测量方法
CN110308490A (zh) 基于磁传感器阵列的旋转运动磁异常目标探测装置及方法
CN101476860A (zh) 一种高背景磁场下的磁定位方法及装置
CN207440306U (zh) 一种深埋废弃井探测装置
CN109901205A (zh) 一种水下机器人多传感器融合与运动轨迹预测方法
CN107044854A (zh) 一种基于标量磁力仪阵列对远距离磁性目标定位的方法
CN109141407A (zh) 一种长输地下管道定位***的误差补偿***及方法
CN110082611A (zh) 一种电场测量装置的定位方法
CN102928884A (zh) 一种磁传感器定位方法
Liu et al. Compensation of geomagnetic vector measurement system with differential magnetic field method
Getscher et al. Magnetic gradient tensor framework for attitude-free position estimation
Pei et al. Magnetic gradiometer on an AUV for buried object detection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant