CN102928794A - 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 - Google Patents
一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 Download PDFInfo
- Publication number
- CN102928794A CN102928794A CN2012104043740A CN201210404374A CN102928794A CN 102928794 A CN102928794 A CN 102928794A CN 2012104043740 A CN2012104043740 A CN 2012104043740A CN 201210404374 A CN201210404374 A CN 201210404374A CN 102928794 A CN102928794 A CN 102928794A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- rightarrow
- target
- sensors
- gradient tensor
- 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.)
- Pending
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于三轴矢量磁传感器阵的磁性目标磁矩的测量方法,该方法利用由五个三轴矢量磁传感器组成磁梯度张量测量阵列来测量磁性目标的磁矩,能够进行单点实时测量,测量精度高。具体方案为:首先利用磁梯度张量测量阵列测量磁性目标周围在探测点处的磁梯度张量;然后利用磁梯度张量对磁性目标与探测点的相对距离及相对坐标进行解算;最后根据所解算出的相对距离及相对位置坐标,结合磁矩求解公式计算出磁性目标在当前坐标系下的三轴磁矩。该测量方法只需要测量阵列在某一个探测点的测量数据即可计算出磁性目标的磁矩,可以达到较高的测量精度。
Description
技术领域
本发明涉及一种高精度的磁性目标磁矩测量方法,具体涉及一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法,属于磁测量领域。
背景技术
磁性目标的磁矩是一个重要物理量,广泛应用于各种理论推导和实际工程中,例如通过推导磁性体的磁矩可以进行磁性目标体的反演、通过磁矩可以对舰船磁场进行外推从而进行舰船消磁或者舰船磁隐身,而在磁探测定位应用领域磁矩这一物理量更是显得尤为重要。
目前还没有一种能够实时测量磁性目标磁矩的方法。现有的磁性目标磁矩测量方法都是一种事后处理分析方法,并且精度不高,其中最常用的就是根据所建立的复杂的磁性目标磁场模型推导出相应的磁矩求解方程组,然后再根据所测量的大量的磁场数据进行反演,这一过程,涉及到求解非线性方程组,常见的求解方法有神经网络法,逐次回归法,遗传算法搜索求解,POWELL法,遗传算法和单纯形法联合求解等。
上述这种传统的目标磁矩确定方法的存在数据量大,计算复杂,求解结果不稳定,实时性差,测量误差大等缺点。
发明内容
有鉴于此,本发明提出一种基于三轴矢量磁传感器阵的磁性目标磁矩的测量方法,该方法利用由五个三轴矢量磁传感器组成磁梯度张量测量阵列来测量磁性目标的磁矩,能够进行单点实时测量,测量精度高。
采用该方法进行磁性目标磁矩的测量的步骤为:
步骤一:布置磁梯度张量测量阵列
磁梯度张量测量阵列由五个三轴矢量磁传感器组成,布置磁梯度张量测量阵列时,位于探测点处的三轴矢量磁传感器为基准传感器。其它三轴矢量磁传感器分布在基准传感器周围,为参考传感器。设基准传感器的三个敏感轴方向分别为x向、y向和z向,布置磁梯度张量测量阵列时应保证所有三轴矢量磁传感器对应的三个敏感轴均相互平行。在x向的正向和负向各布置一个参考传感器,在y向的正向和负向各布置一个参考传感器,四个参考传感器与基准传感器之间的距离均相等。
步骤二:利用步骤一所建立的磁梯度张量测量阵列测量磁性目标在探测点处的磁梯度张量G和磁感应强度B
设基准传感器为1号传感器。位于x向正向上的参考传感器为2号传感器,位于x向负向上的参考传感器为4号传感器,位于y向正向上的参考传感器为3号传感器,位于y向负向上的参考传感器为5号传感器。四个参考传感器与基准传感器之间的距离均为d。
五个三轴矢量磁传感器输出的数据分别为:(B1x,B1y,B1z),(B2x,B2y,B2z)、(B3x,B3y,B3z)、(B4x,B4y,B4z)、(B5x,B5y,B5z)。
其中Bab表示第a个三轴矢量磁传感器测得的在b方向上的磁感应强度,a=1、2、3、4、5,b=x、y、z。
则磁性目标在探测点处的磁梯度张量G为:
所述磁性目标在探测点处的磁感应强度B为:(B1x,B1y,B1z) (2)
式(1)中Bij表示i方向上的磁感应强度在j方向上的导数,i=x、y、z,j=x、y、z。
步骤三:利用磁梯度张量G和磁感应强度B计算磁梯度张量测量阵列与磁性目标的相对位置关系
将步骤二得到的磁性目标在探测点处磁感应强度B和磁梯度张量G代入公式(3),则探测点在以磁偶极子为原点的坐标系中的位置坐标(x,y,z)为:
步骤四:依据步骤三计算的磁梯度张量测量阵列和磁性目标的相对位置计算磁性目标在当前坐标系下的三轴磁矩
有益效果
本发明利用由五个三分量矢量磁传感器组成测量阵列进行磁性目标的磁矩测量,能够进行单点实时测量,测量精度高,且计算过程简单。
附图说明
图1五个三轴矢量磁传感器组成的磁测阵列的结构示意图;
图2为该方法的计算流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本实施例提供一种基于三轴矢量磁传感器阵的磁性目标磁矩的测量方法,该方法利用五个三轴矢量磁传感器组成的磁梯度张量测量阵列来测量磁性目标的磁矩。
磁梯度张量测量阵列的布置方式如图1所示,设五个三轴矢量磁传感器的编号分别为1号~5号。将布置在探测点处的1号三轴矢量磁传感器作为基准传感器,其它三轴矢量磁传感器作为参考传感器。设基准传感器的三个敏感轴方向分别为x向、y向和z向,建立坐标系xyz,所有三轴矢量磁传感器对应的三个敏感轴均相互平行。四个参考传感器中2号和4号三轴矢量磁传感器分别布置在x向的正向和负向,3号和5号三轴矢量磁传感器分别布置在y向的正向和负向;四个参考传感器与基准传感器之间的距离均为d。
采用该磁梯度张量测量阵列进行磁性目标磁矩测量的具体步骤如图2所示:
步骤一:按照上述要求布置磁梯度张量测量阵列;
步骤二:利用磁梯度张量测量阵列测量磁性目标在探测点处的磁梯度张量G和磁感应强度B
磁场是矢量场,在磁性目标周围存在磁场,磁场的二阶张量称为磁梯度张量,记为G。磁梯度张量G的表达式为:
其中Bij表示i方向上的磁感应强度在j方向上的导数(i=x、y、z,j=x、y、z)。
通过电磁场理论,磁性物体所产生的磁场为无源无旋场,因此,磁梯度张量G的矩阵具有对称性,即当i≠j时,Bij=Bji,且有Bzz-Bxx-Byy=0
则有:
图1中基准传感器用于测量磁性目标的磁感应强度B,参考传感器用于测量磁性目标的磁梯度张量G。设1号~5号三轴矢量磁传感器输出的数据分别为:(B1x,B1y,B1z),(B2x,B2y,B2z)、(B3x,B3y,B3z)、(B4x,B4y,B4z)、(B5x,B5y,B5z),利用差分运算代替式(1)中的微分运算,则磁性目标在探测点处的磁梯度张量G为:
磁性目标在探测点处的磁感应强度B为:B=(B1x,B1y,B1z) (3)
其中Bab表示第a个三轴矢量磁传感器测得的在b方向上的磁感应强度(a=1、2、3、4、5,b=x、y、z)
步骤三:利用磁梯度张量G计算磁性目标与磁梯度张量测量阵列间的相对坐标
由于任何磁性目标在一定的距离外均可以等效为磁偶极子模型,因此,将待测量的磁性目标等效为磁偶极子,以磁偶极子为坐标原点,设探测点与磁偶极子间的距离为,则探测点处的磁感应强度B为:
结合式(1)和式(4)有:
式中G-1为磁梯度张量G的逆矩阵。
将步骤一中的式(2)和式(3)代入式(5),则探测点在以磁偶极子为原点的坐标系中的位置坐标(x,y,z)为:
步骤四:计算磁性目标在当前坐标系下的三轴磁矩
将式(5)、式(6)、式(7)写成矩阵形式,即:
令 则上式(8)可以写成:
式(11)中C为常数,磁偶极子在探测点所产生的三分量磁感应强度值B由基准传感器测得;该探测点在以磁偶极子为中心的三轴坐标系内的坐标值已由式(6)得到,则该磁偶极子(即磁性目标)在当前坐标系下三轴磁矩的计算公式为:
将上述磁矩测量方法利用五传感器磁测阵列对某磁性目标进行磁矩测量,其测量值与真实值的对比如表1所示。
表1磁性目标的三轴磁矩
三轴磁矩 | 真实值 | 测量值 |
mx | 150A·m2 | 152.9535A·m2 |
my | 0 | 0.2865A·m2 |
mz | 0 | 0.1307A·m2 |
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于三轴矢量磁传感器阵的磁性目标磁矩的测量方法,其特征在于,
步骤一:布置磁梯度张量测量阵列
所述磁梯度张量测量阵列由五个三轴矢量磁传感器组成,其中位于探测点处的三轴矢量磁传感器为基准传感器;其它三轴矢量磁传感器分布在基准传感器周围,为参考传感器;设基准传感器的三个敏感轴方向分别为x向、y向和z向,布置磁梯度张量测量阵列时应保证所有三轴矢量磁传感器对应的三个敏感轴均相互平行;在x向的正向和负向各布置一个参考传感器,在y向的正向和负向各布置一个参考传感器,四个参考传感器与基准传感器之间的距离均相等;
步骤二:利用步骤一所建立的磁梯度张量测量阵列测量磁性目标在探测点处的磁梯度张量G和磁感应强度B
设基准传感器为1号传感器;位于x向正向上的参考传感器为2号传感器,位于x向负向上的参考传感器为4号传感器,位于y向正向上的参考传感器为3号传感器,位于y向负向上的参考传感器为5号传感器;四个参考传感器与基准传感器之间的距离均为d;
五个三轴矢量磁传感器输出的数据分别为:(B1x,B1y,B1z),(b2x,B2y,B2z)、(B3x,B3y,B3z)、(B4x,B4y,B4z)、(B5x,B5y,B5z);
其中Bab表示第a个三轴矢量磁传感器测得的在b方向上的磁感应强度,a=1、2、3、4、5,b=x、y、z;
则磁性目标在探测点处的磁梯度张量G为:
所述磁性目标在探测点处的磁感应强度B为:(B1x,B1y,B1z) (2)
式(1)中Bij表示i方向上的磁感应强度在j方向上的导数,i=x、y、z,j=x、y、z;
步骤三:利用磁梯度张量G和磁感应强度B计算磁梯度张量测量阵列与磁性目标的相对位置关系
将步骤二得到的磁性目标在探测点处磁感应强度B和磁梯度张量G代入公式(3),则探测点在以磁偶极子为原点的坐标系中的位置坐标(x,y,z)为:
步骤四:依据步骤三计算的磁梯度张量测量阵列和磁性目标的相对位置计算磁性目标在当前坐标系下的三轴磁矩
距离磁偶极子处的磁感应强度为B为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104043740A CN102928794A (zh) | 2012-10-22 | 2012-10-22 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104043740A CN102928794A (zh) | 2012-10-22 | 2012-10-22 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102928794A true CN102928794A (zh) | 2013-02-13 |
Family
ID=47643638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012104043740A Pending CN102928794A (zh) | 2012-10-22 | 2012-10-22 | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102928794A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559376A (zh) * | 2013-11-25 | 2014-02-05 | 中国舰船研究设计中心 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
CN103885002A (zh) * | 2014-03-06 | 2014-06-25 | 中国船舶重工集团公司第七一〇研究所 | 磁传感器阵列测量中的平行度误差补偿方法和*** |
CN103995239A (zh) * | 2014-05-09 | 2014-08-20 | 北京航空航天大学 | 一种磁场梯度参数测量新方法 |
CN105467000A (zh) * | 2015-12-21 | 2016-04-06 | 中国石油大学(北京) | 埋地管道管体缺陷非开挖检测方法及装置 |
CN106404891A (zh) * | 2016-08-29 | 2017-02-15 | 中国石油大学(北京) | 一种埋地管道管体缺陷非开挖检测方法及装置 |
CN106405658A (zh) * | 2016-08-30 | 2017-02-15 | 中国人民解放军海军工程大学 | 一种基于矢量磁梯度计的运动式磁性目标定位方法 |
CN107607999A (zh) * | 2017-08-22 | 2018-01-19 | 哈尔滨工程大学 | 一种对铁磁目标远场磁矩矢量的测量方法 |
CN107817457A (zh) * | 2017-10-13 | 2018-03-20 | 北京工业大学 | 一种地磁梯度张量测量阵列的设计方法 |
CN109613458A (zh) * | 2019-01-14 | 2019-04-12 | 北京卫星环境工程研究所 | 航天器磁异常定位方法 |
CN109781093A (zh) * | 2018-10-09 | 2019-05-21 | 武昌工学院 | 一种水下目标磁性定位跟踪***与方法 |
CN110146839A (zh) * | 2019-05-30 | 2019-08-20 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110568387A (zh) * | 2019-09-16 | 2019-12-13 | 吉林大学 | 一种基于磁梯度张量的航天器磁矩测试方法 |
CN112050800A (zh) * | 2020-08-19 | 2020-12-08 | 哈尔滨工程大学 | 一种基于日字型三轴磁强计对称配置平面阵列的磁梯度张量定位方法 |
CN112327230A (zh) * | 2020-10-28 | 2021-02-05 | 吉林大学 | 一种基于磁梯度张量反演磁化率张量的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7342399B1 (en) * | 2006-04-17 | 2008-03-11 | The United States Of America As Represented By The Secretary Of The Navy | Magnetic anomaly sensing-based system for tracking a moving magnetic target |
US20090006009A1 (en) * | 2007-06-26 | 2009-01-01 | Peter Victor Czipott | Method and system for improving target localization and characterization |
CN101361660A (zh) * | 2008-05-16 | 2009-02-11 | 深圳先进技术研究院 | 一种多磁性目标的定位方法及定位*** |
CN102540110A (zh) * | 2011-12-19 | 2012-07-04 | 北京卫星环境工程研究所 | 磁矩测量方法与测量线圈组 |
-
2012
- 2012-10-22 CN CN2012104043740A patent/CN102928794A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7342399B1 (en) * | 2006-04-17 | 2008-03-11 | The United States Of America As Represented By The Secretary Of The Navy | Magnetic anomaly sensing-based system for tracking a moving magnetic target |
US20090006009A1 (en) * | 2007-06-26 | 2009-01-01 | Peter Victor Czipott | Method and system for improving target localization and characterization |
CN101361660A (zh) * | 2008-05-16 | 2009-02-11 | 深圳先进技术研究院 | 一种多磁性目标的定位方法及定位*** |
CN102540110A (zh) * | 2011-12-19 | 2012-07-04 | 北京卫星环境工程研究所 | 磁矩测量方法与测量线圈组 |
Non-Patent Citations (3)
Title |
---|
吴招才等: "磁力梯度张量测量及应用", 《地质科技情报》 * |
张昌达: "航空磁力梯度张量测量——航空磁测技术的最新进展", 《工程地球物理学报》 * |
张朝阳等: "磁性目标的单点磁梯度张量定位方法", 《探测与控制学报》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559376A (zh) * | 2013-11-25 | 2014-02-05 | 中国舰船研究设计中心 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
CN103559376B (zh) * | 2013-11-25 | 2016-08-17 | 中国舰船研究设计中心 | 基于等效磁矩的潜航平台磁异特征场表征方法 |
CN103885002A (zh) * | 2014-03-06 | 2014-06-25 | 中国船舶重工集团公司第七一〇研究所 | 磁传感器阵列测量中的平行度误差补偿方法和*** |
CN103885002B (zh) * | 2014-03-06 | 2016-06-15 | 中国船舶重工集团公司第七一〇研究所 | 磁传感器阵列测量中的平行度误差补偿方法和*** |
CN103995239A (zh) * | 2014-05-09 | 2014-08-20 | 北京航空航天大学 | 一种磁场梯度参数测量新方法 |
CN103995239B (zh) * | 2014-05-09 | 2016-10-05 | 北京航空航天大学 | 一种磁场梯度参数测量新方法 |
CN105467000A (zh) * | 2015-12-21 | 2016-04-06 | 中国石油大学(北京) | 埋地管道管体缺陷非开挖检测方法及装置 |
CN105467000B (zh) * | 2015-12-21 | 2019-05-21 | 中国石油大学(北京) | 埋地管道管体缺陷非开挖检测方法及装置 |
WO2017107913A1 (zh) * | 2015-12-21 | 2017-06-29 | 中国石油大学(北京) | 埋地管道管体缺陷非开挖检测方法及装置 |
CN106404891A (zh) * | 2016-08-29 | 2017-02-15 | 中国石油大学(北京) | 一种埋地管道管体缺陷非开挖检测方法及装置 |
CN106405658B (zh) * | 2016-08-30 | 2018-03-27 | 中国人民解放军海军工程大学 | 一种基于矢量磁梯度计的运动式磁性目标定位方法 |
CN106405658A (zh) * | 2016-08-30 | 2017-02-15 | 中国人民解放军海军工程大学 | 一种基于矢量磁梯度计的运动式磁性目标定位方法 |
CN107607999A (zh) * | 2017-08-22 | 2018-01-19 | 哈尔滨工程大学 | 一种对铁磁目标远场磁矩矢量的测量方法 |
CN107817457A (zh) * | 2017-10-13 | 2018-03-20 | 北京工业大学 | 一种地磁梯度张量测量阵列的设计方法 |
CN107817457B (zh) * | 2017-10-13 | 2020-03-17 | 北京工业大学 | 一种地磁梯度张量测量阵列的设计方法 |
CN109781093B (zh) * | 2018-10-09 | 2021-07-27 | 杨露雪 | 一种水下目标磁性定位跟踪***与方法 |
CN109781093A (zh) * | 2018-10-09 | 2019-05-21 | 武昌工学院 | 一种水下目标磁性定位跟踪***与方法 |
CN109781093B8 (zh) * | 2018-10-09 | 2021-08-13 | 武昌首义学院 | 一种水下目标磁性定位跟踪***与方法 |
CN109613458A (zh) * | 2019-01-14 | 2019-04-12 | 北京卫星环境工程研究所 | 航天器磁异常定位方法 |
CN109613458B (zh) * | 2019-01-14 | 2020-12-01 | 北京卫星环境工程研究所 | 航天器磁异常定位方法 |
CN110146839A (zh) * | 2019-05-30 | 2019-08-20 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110146839B (zh) * | 2019-05-30 | 2021-11-16 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110568387A (zh) * | 2019-09-16 | 2019-12-13 | 吉林大学 | 一种基于磁梯度张量的航天器磁矩测试方法 |
CN112050800A (zh) * | 2020-08-19 | 2020-12-08 | 哈尔滨工程大学 | 一种基于日字型三轴磁强计对称配置平面阵列的磁梯度张量定位方法 |
CN112327230A (zh) * | 2020-10-28 | 2021-02-05 | 吉林大学 | 一种基于磁梯度张量反演磁化率张量的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102928794A (zh) | 一种基于三轴矢量磁传感器阵的磁性目标磁矩测量方法 | |
CN202815202U (zh) | 基于九个三轴矢量磁传感器的磁性目标磁矩测量装置 | |
CN102927981A (zh) | 一种基于三轴矢量磁传感器阵的磁性目标定位方法 | |
CN202975350U (zh) | 基于五个三轴矢量磁传感器的磁性目标定位装置 | |
CN110333536B (zh) | 一种测距线性定位算法 | |
CN102865862B (zh) | 全站仪免置平自由设站的测量方法 | |
CN104374385B (zh) | 一种海底磁传感器阵列目标定位的新方法 | |
CN105091880B (zh) | 一种基于标量传感器阵列的追踪定位水下远距离磁性目标的方法 | |
Li et al. | An efficient closed-form solution for acoustic emission source location in three-dimensional structures | |
CN103792513B (zh) | 一种雷声定位***及方法 | |
CN102540110B (zh) | 磁矩测量线圈组 | |
CN106990440A (zh) | 一种基于两个探测位置处磁场空间梯度的潜艇定位方法 | |
CN105589108A (zh) | 基于不同约束条件的瞬变电磁快速三维反演方法 | |
CN103852795A (zh) | 一种水下小目标的磁异常信号提取方法 | |
CN111551996A (zh) | 基于磁张量的合作磁目标实时定位方法 | |
CN109597125B (zh) | 一种基于p波到时与最大振幅波形的微震源定位方法 | |
CN103955002B (zh) | 基于磁异常一阶导数的磁偶极子目标位置测量方法 | |
CN112749493A (zh) | 基于全磁力梯度张量特征值的地质体边界检测方法及*** | |
CN104776827A (zh) | Gps高程异常数据的粗差探测方法 | |
CN104363649A (zh) | 带有约束条件的ukf的wsn节点定位方法 | |
CN104182648A (zh) | 反演航天器内部多磁源分布的方法 | |
JP5866893B2 (ja) | 地磁気測定装置 | |
CN111399066A (zh) | 一种基于正交基函数处理标量磁异常梯度信号的方法 | |
CN105046046A (zh) | 一种集合卡尔曼滤波局地化方法 | |
CN112362048A (zh) | 一种实用的磁梯度张量高精度单点定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C05 | Deemed withdrawal (patent law before 1993) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130213 |