CN108227037A - 一种降低三分量磁测系统载体干扰的差分磁补偿方法 - Google Patents

一种降低三分量磁测系统载体干扰的差分磁补偿方法 Download PDF

Info

Publication number
CN108227037A
CN108227037A CN201810024784.XA CN201810024784A CN108227037A CN 108227037 A CN108227037 A CN 108227037A CN 201810024784 A CN201810024784 A CN 201810024784A CN 108227037 A CN108227037 A CN 108227037A
Authority
CN
China
Prior art keywords
magnetic
interference
compensation
magnetic field
difference
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.)
Granted
Application number
CN201810024784.XA
Other languages
English (en)
Other versions
CN108227037B (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201810024784.XA priority Critical patent/CN108227037B/zh
Publication of CN108227037A publication Critical patent/CN108227037A/zh
Application granted granted Critical
Publication of CN108227037B publication Critical patent/CN108227037B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/081Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/40Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Electromagnetism (AREA)
  • Manufacturing & Machinery (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明涉及一种降低三分量磁测系统载体干扰的差分磁补偿方法,建立差分磁补偿工作平台,建立笛卡尔坐标系,对磁原体初始位置参数进行修正,计算干扰磁场值,通过两路磁力仪测量值差分,便可计算出干扰磁场Hd,将计算出的ΔK、Hd带入便求出背景地磁场He。本发明与现有的椭圆拟合法和参数估计法相比,本发明采用磁源体和两路三分量磁力仪位置参数以及磁测数据求解出每一组磁测数据对应的干扰磁场强度,有效地解决了因剩磁变化导致补偿精度偏低问题;本发明引入差分思想,有效消除了外界环境干扰,求取的干扰磁场强度更加准确;补偿处理计算量小,能够对磁测数据进行实时高精度补偿,缩短了测量时间,提高了效率,降低了野外工作成本,具有较高的实用性。

Description

一种降低三分量磁测系统载体干扰的差分磁补偿方法
技术领域
本发明涉及一种降低三分量磁测系统载体干扰的差分磁补偿方法,尤其是对三分量磁测系统载体上的磁源体产生的固定干扰磁场和感应干扰磁场进行有效补偿。
技术背景
三分量磁测方法利用三分量磁力仪测量地磁场分量信息,与传统的总场测量相比,可获得更丰富的地磁场信息,有效减少反演中的多解性,有助于对磁性体定量解释,提高地下矿体探测分辨率和定位精度,在地质调查、矿产勘探和地球科学研究工作中发挥着重要作用。
三分量磁测系统一般由三分量磁力仪和载体组成,载体上的磁源体会产生磁干扰,因此需要对磁测数据进行补偿处理。Tolles和Lawson提出了Tolles-Lawson磁补偿模型补偿飞机产生的磁干扰;李季,张琦,潘孟春,等.《载体干扰磁场补偿方法》国防科技大学学报,2013,35(3),将参数估计方法与Tolles-Lawson模型有效结合,引入信赖域法估计模型中的12个补偿参数,但是该方法所使用的Tolles-Lawson模型过于理想化,无法全面、准确表达干扰磁场特性,且长时间磁测后磁源体的剩磁会发生变化,因此需要从新估计补偿参数。
张晓明等《基于椭圆约束的新型载体磁场标定及补偿技术》,仪器仪表学报,2009,30(11),提出了一种基于椭圆假设的补偿方法,该方法通过磁测数据拟合出一个椭圆,进而求解出各补偿参数;这种方法要求磁测数据必须确保能够拟合出一个椭圆,限制了其实用性。
杨云涛等《一种基于磁偶极子磁场分布理论的磁场干扰补偿方法》,兵工学报,2008,29(12),基于磁偶极子磁场分布理论提出了一种补偿方法;该方法首先建立了地磁测量模型,然后测量4组载体处于不同姿态下三分量磁力仪的测量值以及该姿态所对应的无外界干扰情况下的地磁场强度,最后利用这4组测量值和地磁场强度计算出磁源体的固定干扰磁场和载体三轴方向的磁化系数,进而求解出地磁测量模型中的12个补偿参数;该方法求解出的补偿参数精度较低,而且该方法未给出明确的磁源体位置参数获取方法,导致建立的地磁测量模型无法准确表达磁源体产生的干扰磁场特性,限制了其补偿精度。
现有的磁补偿方法主要有椭圆拟合法和基于Tolles-Lawson磁补偿模型的参数估计法;其中,椭球拟合法需要获取载体全空间姿态下磁测数据,进而确保磁测数据能够拟合出一个椭圆,该方法计算量大,补偿过程过于复杂,实用性大高;参数估计法所使用的Tolles-Lawson磁补偿模型过于理想化,无法全面、准确表达干扰磁场特性;杨云涛建立的地磁测量模型存在求解出的补偿参数精度低,无法准确表达磁源体产生的干扰磁场特性问题,限制了其补偿精度;此外,长时间磁测后磁源体的剩磁会发生变化,降低补偿精度,因此长时间磁测后上述方法均需要从新估计补偿参数。
发明内容
本发明的目的是针对上述现有技术的不足,尤其是三分量磁测过程中容易受到载体上的磁源体干扰的问题,提出一种降低三分量磁测系统载体干扰的差分磁补偿方法。
本发明的目的是通过以下技术方案实现的:
一种降低三分量磁测系统载体干扰的差分磁补偿方法,包括以下步骤:
a、建立差分磁补偿工作平台;选取一块平整的无磁平板作为工作平台,然后将磁源体和两路三分量磁力仪刚性固定在工作平台上,建立差分磁补偿工作平台;
b、建立笛卡尔坐标系,确定磁源体和两路三分量磁力仪的位置参数;以磁源体所在位置中心为坐标系原点构建笛卡尔坐标系,定义两路三分量磁力仪位置(x1,y1,z1)、(x2,y2,z2),设定磁源体几何中心点作为其在差分补偿模型中的初始位置参数(x0,y0,z0);
c、对磁原体初始位置参数进行修正;
d、计算干扰磁场值;将两路三分量磁力仪和磁源体位置参数带入
计算出矩阵ΔK=K1-K2
e、ΔH=H1-H2通过两路磁力仪测量值差分获取,将ΔK和ΔH带入公式(7)便计算出干扰磁场Hd
f、求解地磁场He,将步骤d计算出的ΔK、Hd带入
便求出背景地磁场He
步骤c所述的磁源体位置参数修正,包括以下步骤:
c1首先选取一块磁场变化平缓且周围无磁干扰物的区域,测量出该区域的地磁总场强度T;
c2然后改变差分磁补偿系统的姿态,记录不同姿态下两路三分量磁力仪的测量值c=1,2,3...N;以磁源体位置参数为待估计参数,测量值和磁源体的初始位置参数带入公式(8)计算出N组干扰磁场
c3最后将求取的干扰磁场带入式(5)计算出补偿后的地磁场He,进而利用非限制最小二乘法对磁源体位置参数进行修正。
有益效果:本发明与现有常用的椭圆拟合法和参数估计法两种磁补偿方法相比,本发明采用磁源体和两路三分量磁力仪位置参数以及磁测数据求解出每一组磁测数据对应的干扰磁场强度,有效地解决了因剩磁变化导致补偿精度偏低问题;本发明与现有最接近的“基于磁偶极子磁场分布理论的磁补偿方法”相比,本发明采用非线性最小二乘法对磁源体位置参数进行修正,建立的差分磁补偿模型能够更加准确地表达磁源体产生的干扰磁场特性,而且本发明引入差分思想,有效消除了外界环境干扰,求取的干扰磁场强度更加准确;本发明补偿处理计算量小,能够对磁测数据进行实时高精度补偿,缩短了测量时间,提高了效率,降低了野外工作成本,具有较高的实用性。
附图说明:
附图1为一种降低三分量磁测系统载体干扰的差分磁补偿方法流程图
附图2为磁源体与磁力仪在补偿平台相对位置示意图
附图3为补偿前后误差对比图
附图4为位置参数修正前后误差对比图
具体实施方式:
下面结合附图和实施例对本发明作进一步的详细说明:
一种降低三分量磁测系统载体干扰的差分磁补偿方法,包括以下步骤:
a、建立差分磁补偿工作平台;选取一块平整的无磁平板作为工作平台,然后将磁源体和两路三分量磁力仪刚性固定在工作平台上,建立差分磁补偿工作平台;
b、建立笛卡尔坐标系,确定磁源体和两路三分量磁力仪的位置参数;以磁源体所在位置中心为坐标系原点构建笛卡尔坐标系,定义两路三分量磁力仪位置(x1,y1,z1)、(x2,y2,z2),设定磁源体几何中心点作为其在差分补偿模型中的初始位置参数(x0,y0,z0);
c、对磁原体初始位置参数进行修正;
d、计算干扰磁场值;将两路三分量磁力仪和磁源体位置参数带入
计算出矩阵ΔK=K1-K2
e、ΔH=H1-H2通过两路磁力仪测量值差分获取,将ΔK和ΔH带入公式(7)便计算出干扰磁场Hd
f、求解地磁场He,将步骤d计算出的ΔK、Hd带入
便求出背景地磁场He
步骤c所述的磁源体位置参数修正,包括以下步骤:
c1首先选取一块磁场变化平缓且周围无磁干扰物的区域,测量出该区域的地磁总场强度T;
c2然后改变差分磁补偿系统的姿态,记录不同姿态下两路三分量磁力仪的测量值c=1,2,3...N;以磁源体位置参数为待估计参数,测量值和磁源体的初始位置参数带入公式(8)计算出N组干扰磁场
c3最后将求取的干扰磁场带入式(5)计算出补偿后的地磁场He,进而利用非限制最小二乘法对磁源体位置参数进行修正。
实施例1
建立差分磁补偿模型:
磁场解析模型;根据磁偶极子磁场解析模型,空间一点P(x′,y′,z′)处的感应矢量磁异常为:
其中,μ(Tm/A)为磁导率,空气的磁导率μ=4π×10-7Tm/A;B=Bxi+Byj+Bzk(T)为点P处的感应矢量磁异常;M=Mxi+Myj+Mzk(Agm2)为磁源矢量磁矩;r=rxi+ryj+rzk(m)为点P到磁源中心间矢径。
定义磁偶极子位置为(x,y,z),用矩阵的形式展开式(1):
定义矩阵K为:
则式(2)可化简为:
B=K gM (4)
a:建立差分磁补偿模型;根据式(4),两路三分量磁力仪测量值表达式为:
其中,Hp为固定干扰磁场;Hi为感应干扰磁场;为两路三分量磁力仪对应的背景地磁场;K为3×3系数矩阵。
由于地磁场各分量的梯度较小,对于两路三分量磁力仪来说,背景地磁场属于共模信号,即:根据以上分析,对式(5)进行差分处理,得到差分磁补偿模型:
H1-H2=(K1-K2)g(Hp+Hi) (6)
b:建立笛卡尔坐标系,确定磁源体和两路三分量磁力仪的位置参数;选取一块平整的无磁平板作为工作平台,然后将磁源体和两路三分量磁力仪刚性固定在工作平面上,以磁源体所在位置的中心点为坐标系原点构建笛卡尔坐标系,定义两路三分量磁力仪位置(x1,y1,z1)、(x2,y2,z2);设定磁源体几何中心点作为其在差分补偿模型中的位置参数(x0,y0,z0);则式(6)可以展开为:
其中,ΔH=H1-H2,ΔK=K1-K2,Hd=Hp+Hi.
式(7)中矩阵ΔK为常数矩阵,此外ΔH可通过两路磁力仪测量值差分获取,因此可利用式(7)计算出干扰磁场Hd
将计算出的Hd带入式(5)便可求取地磁场He
磁源体位置参数修正方法如下:
D、首先选取一块磁场变化平缓且周围无磁干扰物的区域,测量出该区域的地磁总场强度T;然后改变差分磁补偿平台的姿态,记录不同姿态下两路三分量磁力仪的测量值(c=1,2,3...N)。
E、以磁源体位置参数为待估计参数,设定磁源体几何中心点作为其位置参数初始值,将测量值和磁源体位置参数带入公式(8)计算出N组干扰磁场最后将求取的干扰磁场带入式(5)计算出补偿后的地磁场He
F、非线性最小二乘法的目标函数可设定为:
其中,Δε表示绝对误差和。
根据非线性最小二乘法准则,待估计参数应使绝对误差和最小,即:
Δε→min (11)
三分量磁测数据实时补偿:
G、将修正后的磁源体位置参数和两路三分量磁力仪位置参数带入式(3)计算出矩阵ΔK;
H、将ΔK和两路三分量磁力仪测量值带入式(8)求取出干扰磁场Hd
I、将矩阵ΔK和干扰磁场Hd带入式(5),计算出补偿后的地磁场He

Claims (2)

1.一种降低三分量磁测系统载体干扰的差分磁补偿方法,其特征在于,包括以下步骤:
a、建立差分磁补偿工作平台;选取一块平整的无磁平板作为工作平台,然后将磁源体和两路三分量磁力仪刚性固定在工作平台上,建立差分磁补偿工作平台;
b、建立笛卡尔坐标系,确定磁源体和两路三分量磁力仪的位置参数;以磁源体所在位置中心为坐标系原点构建笛卡尔坐标系,定义两路三分量磁力仪位置(x1,y1,z1)、(x2,y2,z2),设定磁源体几何中心点作为其在差分补偿模型中的初始位置参数(x0,y0,z0);
c、对磁原体初始位置参数进行修正;
d、计算干扰磁场值;将两路三分量磁力仪和磁源体位置参数带入
计算出矩阵ΔK=K1-K2
e、ΔH=H1-H2通过两路磁力仪测量值差分获取,将ΔK和ΔH带入公式(7)便计算出干扰磁场Hd
f、求解地磁场He,将步骤d计算出的ΔK、Hd带入
便求出背景地磁场He
2.按照权利要求1所述的一种降低三分量磁测系统载体干扰的差分磁补偿方法,其特征在于,步骤c所述的磁源体位置参数修正,包括以下步骤:
c1首先选取一块磁场变化平缓且周围无磁干扰物的区域,测量出该区域的地磁总场强度T;
c2然后改变差分磁补偿系统的姿态,记录不同姿态下两路三分量磁力仪的测量值以磁源体位置参数为待估计参数,测量值和磁源体的初始位置参数带入公式(8)计算出N组干扰磁场
c3最后将求取的干扰磁场带入式(5)计算出补偿后的地磁场He,进而利用非限制最小二乘法对磁源体位置参数进行修正。
CN201810024784.XA 2018-01-11 2018-01-11 一种降低三分量磁测系统载体干扰的差分磁补偿方法 Active CN108227037B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810024784.XA CN108227037B (zh) 2018-01-11 2018-01-11 一种降低三分量磁测系统载体干扰的差分磁补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810024784.XA CN108227037B (zh) 2018-01-11 2018-01-11 一种降低三分量磁测系统载体干扰的差分磁补偿方法

Publications (2)

Publication Number Publication Date
CN108227037A true CN108227037A (zh) 2018-06-29
CN108227037B CN108227037B (zh) 2019-11-29

Family

ID=62641744

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810024784.XA Active CN108227037B (zh) 2018-01-11 2018-01-11 一种降低三分量磁测系统载体干扰的差分磁补偿方法

Country Status (1)

Country Link
CN (1) CN108227037B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6617850B2 (en) * 2001-12-03 2003-09-09 Mayo Foundation For Medical Education And Research Motion correction of magnetic resonance images using phase difference of two orthogonal acquisitions
CN105388533A (zh) * 2015-10-23 2016-03-09 国家海洋局第二海洋研究所 安装于潜器中磁力仪磁干扰的陆上校正方法
CN105510849A (zh) * 2015-11-26 2016-04-20 哈尔滨工业大学 航磁干扰补偿方法
CN106353824A (zh) * 2016-09-29 2017-01-25 吉林大学 航空磁通门磁梯度张量仪的系统校正及磁干扰补偿融合方法
CN106405666A (zh) * 2016-08-29 2017-02-15 东南大学 基于改进的约束总体最小二乘的地磁测量误差补偿方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6617850B2 (en) * 2001-12-03 2003-09-09 Mayo Foundation For Medical Education And Research Motion correction of magnetic resonance images using phase difference of two orthogonal acquisitions
CN105388533A (zh) * 2015-10-23 2016-03-09 国家海洋局第二海洋研究所 安装于潜器中磁力仪磁干扰的陆上校正方法
CN105510849A (zh) * 2015-11-26 2016-04-20 哈尔滨工业大学 航磁干扰补偿方法
CN106405666A (zh) * 2016-08-29 2017-02-15 东南大学 基于改进的约束总体最小二乘的地磁测量误差补偿方法
CN106353824A (zh) * 2016-09-29 2017-01-25 吉林大学 航空磁通门磁梯度张量仪的系统校正及磁干扰补偿融合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄学功 等: ""地磁信号检测系统误差分析与补偿方法研究"", 《兵工学报》 *

Also Published As

Publication number Publication date
CN108227037B (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
CN107272069B (zh) 基于磁异常梯度的磁性目标追踪方法
Liu et al. Calibration and compensation of geomagnetic vector measurement system and improvement of magnetic anomaly detection
CN109856689B (zh) 一种超导航磁梯度张量数据抑噪处理方法和系统
CN107544042B (zh) 一种磁力计阵列校正方法
CN110146839A (zh) 一种移动平台磁梯度张量系统校正方法
CN110274586B (zh) 包含多光系原子磁力仪方向误差补偿的航空磁补偿方法
CN109814163B (zh) 一种基于扩展补偿模型的航磁张量数据抑噪方法及系统
KR101548667B1 (ko) 지자기교란요소 제거를 통한 해양자력탐사 총지자장 측정방법, 시스템 및 이를 위한 기록매체
CN109444774B (zh) 基于小信号法测量水下航行器磁场噪声系数的方法
CN109521384A (zh) 一种基于原子磁强计的矢量磁补偿方法
CN103926627A (zh) 水下载体地磁三分量测量方法
CN108919304B (zh) 一种基于参考平面的移动测量系统中pos误差补偿方法
CN109725361B (zh) 基于磁梯度张量不变量的一种磁性目标定位方法
CN110187393B (zh) 一种基于广义回归神经网络的航磁补偿方法
CN102954804B (zh) 一种石油钻探测量用传感器正交校准方法
CN105466458A (zh) 正六面体磁梯度张量系统误差直接校正方法
CN111399066A (zh) 一种基于正交基函数处理标量磁异常梯度信号的方法
CN112362048A (zh) 一种实用的磁梯度张量高精度单点定位方法
Gang et al. Integrated calibration of magnetic gradient tensor system
Pang et al. Integrated calibration of strap-down geomagnetic vector measurement system
CN108227037B (zh) 一种降低三分量磁测系统载体干扰的差分磁补偿方法
CN108508479B (zh) 一种空地井立体重磁数据协同目标位置反演方法
CN114740540B (zh) 一种基于方向约束的洋中脊区磁异常图构建方法及系统
CN108254796A (zh) 一种标量磁力仪阵列基线的优化方法
CN113985494A (zh) 一种基于无迹卡尔曼算法海底地震计中电子罗盘误差补偿方法

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