CN107817457B - 一种地磁梯度张量测量阵列的设计方法 - Google Patents
一种地磁梯度张量测量阵列的设计方法 Download PDFInfo
- Publication number
- CN107817457B CN107817457B CN201710955114.5A CN201710955114A CN107817457B CN 107817457 B CN107817457 B CN 107817457B CN 201710955114 A CN201710955114 A CN 201710955114A CN 107817457 B CN107817457 B CN 107817457B
- Authority
- CN
- China
- Prior art keywords
- array
- order
- gradient tensor
- geomagnetic
- axis
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种地磁梯度张量测量阵列的设计方法,属于地磁异常反演领域。阵列由9个探头在平面形成菱形排列,每个探头能够测量三轴磁场信号。本发明公开了一阶及二阶地磁梯度张量计算方法,公开了基线距离与一阶、二阶地磁梯度张量测量误差的关系。本发明能够提高地磁梯度张量的测量精度,进一步保证地磁异常反演可靠性。
Description
技术领域
本发明涉及一种地磁梯度张量测量阵列的设计方法,能够实现高精度的一阶及二阶地磁梯度张量测量,属于地磁异常反演领域。
背景技术
地磁场作为地球基本属性,具有相对稳定和均匀的特性。地磁场具有不需要发射源、全天候存在、全地域覆盖等优点。在地磁环境条件下的铁磁物质,能够感应出磁场,从而使得所在位置的磁场信号异于周围磁场。地磁异常信号的强弱和方向,与铁磁物质本身的特性、构造和外形尺寸等相关。通过地磁异常信号反演能得到目标体的形状和姿态特征,从而识别出目标物质的基本属性。
一阶地磁梯度张量由磁场分量在不同方向上的空间变化率组成,二阶地磁梯度张量是一阶地磁梯度张量进一步偏导。相比于单分量磁场、总磁场,一阶和二阶地磁梯度张量不受总场测量的限制,被广泛应用在磁目标定位、磁源几何参数反演等领域。
由于外界干扰、仪器零漂变化、地磁不稳定等随机噪声,给地磁异常信号反演造成了很大困难,对检测探头的灵敏度、准确性和鲁棒性提出了很高的要求。现有的地磁检测多采用单探头、十字形和六面体阵列型式,测量精度易受到噪声影响。
发明内容
本发明目的是提供一种地磁梯度张量测量阵列的设计方法,能够实现目标磁场的高精度测量,提高地磁反演的灵敏度、分辨率和鲁棒性。
为解决上述目的,本发明所采用的技术方案为
一种地磁梯度张量测量阵列的设计方法,测量阵列由9个探头在x-y平面形成菱形排列,每个探头能够测量三轴磁场信号,其中,每4个探头组成一个十字型阵列结构,基线距离d相等,则每个十字型阵列的中心磁场由周围4个探头磁场平均值求得,表示为
式中,i=1、2、3、4,代表4个十字型阵列,j=x、y、z,代表磁场三轴方向。x±、y±分别代表沿x轴和y轴正、反方向。
4个十字型阵列的中点可组成新的十字型阵列,根据公式(1)及张量对称性,阵列中心的一阶地磁梯度张量G表示为
进一步,分别求得4个十字型阵列的一阶地磁梯度张量Gk,表达式为
其中,k=1、2、3、4。G3、G1分别代表沿x轴正方向、反方向一阶地磁梯度张量,G2、G4分别代表沿y轴正方向、反方向一阶地磁梯度张量。
其中,C和D为
进一步,基线距离越小,一阶及二阶地磁梯度张量的测量精度越高。
附图说明
图1是本专利阵列结构示意图
图2是十字形阵列结构示意图
图3是六面体阵列结构示意图
图4是一阶地磁梯度张量误差对比
图5是十字形阵列和本专利阵列的二阶地磁梯度张量误差对比
图6是六面体阵列和本专利阵列的二阶地磁梯度张量误差对比
图7是一阶地磁梯度张量误差与基线距离的关系
图8是二阶地磁梯度张量误差与基线距离的关系
具体实施内容
如图1所示,地磁阵列结构由9个探头组成,每个探头测量三轴磁场。探头1、2、3、5组成十字形阵列;探头2、4、5、7组成十字形阵列;探头3、5、6、8组成十字形阵列;探头5、7、8、9组成十字形阵列,其基线距离d相等。每个十字形阵列中点位置磁场信号表示为
本方法的阵列分别求得4个十字型阵列的一阶地磁梯度张量Gk,k=1、2、3、4,沿x轴正方向的G3表达式为
进一步,阵列同样能够测量一阶地磁梯度张量在不同方向上的偏导数,即二阶地磁梯度张量,表达式为
设置磁源位置坐标为原点,测量平面为1m×1m的正方形区域,每隔50mm作为一个测量点,设置基线距离d=0.01m,提离高度z=1m。磁矩大小为M=(-12500,12500,10000)A·m2。为更好进行对比,引入十字形阵列及六面体阵列,其结构分别如图2和3所示。十字形阵列一阶地磁梯度张量为
十字形阵列无法完整计算二阶磁梯度张量,则部分二阶磁梯度张量的计算公式为
六面体阵列的磁梯度张量G3表示为
六面体阵列也只能计算部分二阶磁梯度张量,其计算公式为
利用式(13)计算一阶地磁梯度张量分量与理论值的误差,即
式中,Gij-a为一阶地磁梯度张量实际值,Gij-t为一阶地磁梯度张量理论值。
由于存在探测盲点,利用三次样条插值完成盲点近似计算。设置地磁场磁感应磁强度B=(27961.9,-2938.9,46792.8)nT。将地磁场背景下的磁场依次叠加噪声强度为0、5、10、15、20、25、30dBw的高斯白噪声,以所有测量点误差的平均值作为表征值,十字形阵列、六面体阵列、本专利阵列的一阶磁梯度张量测量误差如图5所示。可见,随着噪声强度的增加,三种阵列测量误差逐渐增加。在添加任何强度噪声时,本专利阵列测量误差都要比其他两个阵列小,测量精度更高。
二阶磁梯度张量分量与理论值间误差的计算公式为
式中,G′ijk-a为二阶地磁梯度张量测量值,G′ijk-t为二阶地磁梯度张量理论值。由于十字形和六面体只能测量部分二阶地磁梯度张量,所以将本专利阵列和十字形阵列和六面体阵列分别进行对比。将地磁场背景下的磁场依次叠加噪声强度为0、5、10、15、20、25、30dBW的高斯白噪声,以所有测量点误差的平均值作为表征值,十字形阵列和本专利阵列的二阶地磁梯度张量的测量误差如图5所示(n=6)。可见,在任何添加强度噪声时,本专利阵列测量误差都要比十字形阵列小。六面体阵列和本专利阵列的二阶地磁梯度张量的测量误差如图6所示(n=9)。可见,本专利阵列与六面体阵列测量误差接近,但考虑到二阶地磁梯度张量的完整性,本专利阵列更适合二阶地磁梯度张量测量。
设置磁源位置坐标为原点,阵列中心坐标为x=1,y=1,z=1,设置基线距离d变化范围为0.01m~1m,每隔0.01m作为一个距离点。如图7为一阶地磁梯度张量测量误差与基线距离d的关系。图8为二阶地磁梯度张量测量误差与基线距离d的关系。可知,基线距离d越小,一阶及二阶地磁梯度张量的测量精度越高。
Claims (2)
1.一种地磁梯度张量测量阵列的设计方法,其特征在于:阵列由9个探头在x-y平面形成菱形排列,每个探头能够测量三轴磁场信号,其中每4个探头组成一个十字型阵列结构,基线距离d相等,则每个十字型阵列的中心磁场由周围4个探头磁场平均值求得,表示为
式中,i=1、2、3、4,代表4个十字型阵列;j=x、y、z,代表磁场三轴方向;x±、y±分别代表沿x轴和y轴正、反方向;
4个十字型阵列的中点组成新的十字型阵列,根据公式(1)及张量对称性,则阵列中心的一阶地磁梯度张量G可表示为
分别求得4个十字型阵列的一阶地磁梯度张量Gk,表达式为
其中,k=1、2、3、4;G3、G1分别代表沿x轴正方向、反方向一阶地磁梯度张量,G2、G4分别代表沿y轴正方向、反方向一阶地磁梯度张量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710955114.5A CN107817457B (zh) | 2017-10-13 | 2017-10-13 | 一种地磁梯度张量测量阵列的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710955114.5A CN107817457B (zh) | 2017-10-13 | 2017-10-13 | 一种地磁梯度张量测量阵列的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107817457A CN107817457A (zh) | 2018-03-20 |
CN107817457B true CN107817457B (zh) | 2020-03-17 |
Family
ID=61607225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710955114.5A Active CN107817457B (zh) | 2017-10-13 | 2017-10-13 | 一种地磁梯度张量测量阵列的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107817457B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108873086B (zh) * | 2018-06-05 | 2019-09-27 | 哈尔滨工程大学 | 一种采用地磁总场梯度阵列对磁性目标定位的方法 |
CN109188315B (zh) * | 2018-08-31 | 2020-10-13 | 哈尔滨工业大学 | 一种传感器阵列基线距离可调的磁梯度仪 |
CN109307847B (zh) * | 2018-11-06 | 2021-01-05 | 哈尔滨工程大学 | 一种磁体的二阶欧拉反演法 |
CN110095738B (zh) * | 2019-04-15 | 2021-03-12 | 哈尔滨工业大学 | 磁梯度仪参数设计方法 |
CN111596369A (zh) * | 2020-06-01 | 2020-08-28 | 北京卫星环境工程研究所 | 面向埋地未爆弹的磁场梯度探测仪 |
CN111653104B (zh) * | 2020-07-13 | 2022-04-19 | 深圳华强技术有限公司 | 磁场基线自维护方法、装置、计算机设备及存储介质 |
CN112050799B (zh) * | 2020-08-19 | 2022-11-18 | 哈尔滨工程大学 | 一种基于磁梯度张量缩并量之比的测距定位方法 |
CN112050800B (zh) * | 2020-08-19 | 2022-11-18 | 哈尔滨工程大学 | 一种基于日字型三轴磁强计对称配置平面阵列的磁梯度张量定位方法 |
CN113253162B (zh) * | 2021-06-18 | 2022-04-26 | 上海交通大学 | 一种微机电系统磁通门地磁张量传感芯片 |
CN113720904B (zh) * | 2021-09-08 | 2024-05-10 | 中国人民解放军国防科技大学 | 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 |
CN114545489B (zh) * | 2022-03-09 | 2024-05-10 | 中国地质大学(武汉) | 一种基于高阶磁异常导数的地震短临预测方法 |
CN114413884B (zh) * | 2022-03-29 | 2022-07-15 | 西南民族大学 | 基于磁梯度仪阵列的定位方法、装置及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2171523A (en) * | 1985-02-23 | 1986-08-28 | Stc Plc | Magnetic gradient detection |
CN102590765A (zh) * | 2012-02-21 | 2012-07-18 | 大连理工大学 | 一种全张量磁力梯度仪 |
CN102928794A (zh) * | 2012-10-22 | 2013-02-13 | 中国船舶重工集团公司第七一〇研究所 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
CN202815202U (zh) * | 2012-10-22 | 2013-03-20 | 中国船舶重工集团公司第七一〇研究所 | 基于九个三轴矢量磁传感器的磁性目标磁矩测量装置 |
CN104215919A (zh) * | 2014-09-30 | 2014-12-17 | 北京卫星环境工程研究所 | 用于测量磁场梯度张量的装置及其测量方法 |
CN105022003A (zh) * | 2015-06-12 | 2015-11-04 | 于润桥 | 磁场梯度测量装置 |
-
2017
- 2017-10-13 CN CN201710955114.5A patent/CN107817457B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2171523A (en) * | 1985-02-23 | 1986-08-28 | Stc Plc | Magnetic gradient detection |
CN102590765A (zh) * | 2012-02-21 | 2012-07-18 | 大连理工大学 | 一种全张量磁力梯度仪 |
CN102928794A (zh) * | 2012-10-22 | 2013-02-13 | 中国船舶重工集团公司第七一〇研究所 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
CN202815202U (zh) * | 2012-10-22 | 2013-03-20 | 中国船舶重工集团公司第七一〇研究所 | 基于九个三轴矢量磁传感器的磁性目标磁矩测量装置 |
CN104215919A (zh) * | 2014-09-30 | 2014-12-17 | 北京卫星环境工程研究所 | 用于测量磁场梯度张量的装置及其测量方法 |
CN105022003A (zh) * | 2015-06-12 | 2015-11-04 | 于润桥 | 磁场梯度测量装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107817457A (zh) | 2018-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107817457B (zh) | 一种地磁梯度张量测量阵列的设计方法 | |
CN110095738B (zh) | 磁梯度仪参数设计方法 | |
CN111624389B (zh) | 一种用于倾斜多导体系统的电流测量系统及方法 | |
CN102313543B (zh) | 基于巨磁阻传感器的地磁方位角测量系统、测量方法及正交补偿方法 | |
CN110146839A (zh) | 一种移动平台磁梯度张量系统校正方法 | |
CN109725360B (zh) | 基于磁梯度张量不变量的单点定位方法 | |
CN104535062A (zh) | 基于磁梯度张量和地磁矢量测量的运动式定位方法 | |
CN110793439B (zh) | 一种多传感器测量机坐标统一的标准器和坐标统一方法 | |
CN106643792A (zh) | 惯性测量单元和地磁传感器整体标定装置及标定方法 | |
CN111190230A (zh) | 一种基于磁梯度张量的探测方法 | |
CN108344361B (zh) | 一种基于激光跟踪仪的平面法向量测量方法 | |
CN111551996A (zh) | 基于磁张量的合作磁目标实时定位方法 | |
CN107390155B (zh) | 一种磁传感器校准装置和方法 | |
CN109725361B (zh) | 基于磁梯度张量不变量的一种磁性目标定位方法 | |
CN112684511B (zh) | 基于两点磁梯度全张量的线性定位方法 | |
CN111190229B (zh) | 一种磁目标探测方法 | |
CN108398124B (zh) | 一种电子罗盘校准方法 | |
CN103353612B (zh) | 一种地下目标物体的测量定位设备及测量定位方法 | |
CN103307968A (zh) | 机器人承载平台姿态的检测方法 | |
Pang et al. | Misalignment error suppression between host frame and magnetic sensor array | |
CN109633540B (zh) | 一种磁源的实时定位系统及实时定位方法 | |
CN113866688B (zh) | 一种小姿态角条件下的三轴磁传感器误差校准方法 | |
CN109633541B (zh) | 一种磁源定位装置及磁源定位方法 | |
CN110596619B (zh) | 一种全张量磁梯度测量组件及其优化方法 | |
CN112050800A (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 |