CN112344925B - 基于归一化磁矩的水下目标相关性检测定位方法 - Google Patents
基于归一化磁矩的水下目标相关性检测定位方法 Download PDFInfo
- Publication number
- CN112344925B CN112344925B CN202011081070.6A CN202011081070A CN112344925B CN 112344925 B CN112344925 B CN 112344925B CN 202011081070 A CN202011081070 A CN 202011081070A CN 112344925 B CN112344925 B CN 112344925B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- detection
- magnetic moment
- correlation
- moment
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005259 measurement Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 8
- 238000005314 correlation function Methods 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 2
- 230000000996 additive effect Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000005307 time correlation function Methods 0.000 claims 1
- 230000002159 abnormal effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
- G01C21/08—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于归一化磁矩的水下目标相关性检测定位方法，属于磁探技术领域和磁偶极子定位领域，该方法通过最简张量测量磁探仪动态获取归一化磁矩，对其进行相关性检测将有效能量从含噪信号中提取出来，进一步利用特征时间搜索确定过顶点距离。本发明提出的一次性搜索探测方法，能够有效规避目标磁矩方向的影响以及克服信噪比带来的无法定位问题，对特定场景下的磁异常目标定位应用具有重要的实际意义。
Description
技术领域
本发明属于磁探技术领域和磁偶极子定位领域，具体涉及一种基于归一化磁矩的水下目标相关性检测定位方法。
背景技术
由于磁场具有高穿透性、高抗扰性、信号被动接收、低功耗等优势，磁异常探测(MAD：magnetic anomaly detection)已经被广泛应用于遥感探测领域。在磁异常探测领域，一个远距离磁性目标通常可以视为磁偶极子。针对磁偶极子的探测定位问题，通过多传感器布局并采用非线性优化方法进行位置参数的拟合求解，存在系统复杂和计算量大的问题。
在静态磁探测领域，各种基于磁梯度测量的解析算法被广泛研究应用，主要以STAR方法(Wiegert,R.,and J.Oeschger,”Generalized magnetic gradient contractionbased method for detection,localization and discrimination of underwatermines and unexploded ordnance,”Oceans IEEE,vol.2,pp.1325-1332,Sept.2005.)和改进STAR方法等，其优点在于以解析方式求解目标位置等参数，计算量小，磁梯度测量能在一定程度上降低载体和背景磁场的干扰。但这些解析方法由于系统过于复杂，很难应用于动态磁探测领域。
基于正交基的能量检测方法能够提高磁异常探测过程中的信噪比，但传统的正交基方法都是基于背景磁场恒定的假设，与事实不符。
发明内容
针对传统正交基方法存在的无法适应实际动态磁探测需求的缺陷，本发明提供一种基于归一化磁矩的水下目标相关性检测定位方法，该方法从原理上消除了对背景磁场恒定假设的依赖，具备一次探测定位能力。
实现本发明目的的技术方案为：一种基于归一化磁矩的水下目标相关性检测定位方法，运用三轴矢量磁强计构建最简张量测量磁探仪进行磁目标动态定位，具体实现步骤如下：
步骤1、选取磁探仪构型，获取各磁强计三轴测量分量；
步骤2、利用磁强计三轴测量分量，求取磁张量矩阵；
步骤3、求解磁张量矩阵的三个特征值，得到动态磁探测测线上每一点的归一化磁矩；
步骤4、在航磁探测坐标系下，归一化磁矩表示成：其中，μ为相对磁导率，ρ为最小接近距离，/>tρ为磁测仪以速度V沿测量航迹到过靶点N的飞行时间；/>为归一化磁矩的特征函数；
步骤5、将归一化磁矩在设定的时间窗口内进行相关性检测得到检测函数，由检测函数的最大值时刻可获取航磁探测过目标顶点时刻；
步骤6、特征时间搜索获取过目标顶点距离。
与现有技术相比，本发明的有益效果在于：本发明通过将旋转不变量归一化磁矩应用在动态磁探测领域，解决了测量平台姿态不稳以及载体和环境磁干扰问题，基于正交基思想，通过相关性检测有效提高了被探信号的信噪比，并提出特征时间搜索方法获取过顶点距离，提高磁目标动态定位能力，具有实际意义。
附图说明
图1为基于归一化磁矩的水下目标相关性检测定位方法流程图。
图2为张量矩阵计算坐标系示意图。
图3为航磁探测坐标系示意图。
图4为最简张量测量磁探仪架构示意图。
图5为特征函数曲线图。
图6为含高斯白噪声的归一化磁矩测量及相关性检测曲线图。
图7为相关函数随特征时间比值变化关系示意图。
具体实施方式
如图1所示，本发明的一种基于归一化磁矩的水下目标相关性检测定位方法，运用三轴矢量磁强计构建最简张量测量磁探仪进行磁目标动态定位，包括以下步骤：
步骤1.选取磁探仪构型，获取各磁强计三轴测量分量，分别用Bix，Biy，Biz表示，i＝1,2,3；
步骤2.利用磁强计三轴测量分量Bix，Biy，Biz，求取磁张量矩阵：
其中，Txy＝(B3x-B2x)/L，Tyy＝(B3y-B2y)/L，Tzy＝(B3z-B2z)/L，Txz＝Tzx，Tyx＝Txy，Tyz＝Tzy，Tzz＝-Txx-Tyy。L为梯度仪基线即正三角形边长。
步骤3.求解磁张量矩阵的三个特征值分别为λ1，λ2和λ3，满足λ1≥λ2≥λ3，得到动态磁探测测线上每一点的归一化磁矩为：
步骤4.在航磁探测坐标系下，归一化磁矩可以表示成：
其中，μ为相对磁导率，ρ为最小接近距离tρ为磁测仪以速度V沿测量航迹到过靶点N的飞行时间，/>为归一化磁矩的特征函数；
步骤5.将归一化磁矩在[-2,2]的时间窗口内进行相关性检测得到检测函数为：
其中，m为检测时间序列，由检测函数的最大值时刻可获取航磁探测过目标顶点时刻。
步骤6.特征时间搜索获取过目标顶点距离。特征时间τ＝ρ/V，通过搜索获取真实的特征时间可以得到真实的过顶点距离。取搜索特征时间为τs，不同特征时间对应求解检测函数r(τs)，过顶点时刻检测函数最大值对应真实特征时间，从而得到真实的过顶点距离ρ。
在大部分磁探问题中，磁目标距磁探仪距离大于磁目标本体尺寸3倍以上，即磁目标可近似等效为磁偶极子。
下面结合附图对本发明做进一步的说明，为了更好理解本发明的技术方案，现将其原理及计算公式结合具体实施例推导详细叙述如下：
定义归一化磁矩为
如图2所示的张量矩阵计算坐标系其中/>为目标磁矩矢量，/>为测点相对目标位移矢量。则测点坐标为(r,0,0)，磁矩矢量为(Mcosφ,0,Msinφ)。
其中，φ为磁矩矢量和相对位置矢量的夹角，则张量矩阵在此坐标系下表示为：
特征方程为
求解特征方程得三个特征值分别为：
从上式可以得出归一化磁矩可以由磁张量矩阵的特征值表示成：
其中，λ1≥λ2≥λ3。
航磁探测坐标系如图3所示。设磁偶极子目标位于坐标原点，航磁测量航迹与X轴方向平行。定义N(0,yρ,zρ)为过顶点，即ON垂直于测量直线航迹，航迹方程为x＝V(t-tρ)，y＝yρ，z＝zρ。
式中tρ为磁测仪以速度V沿测量航迹到过靶点N的飞行时间，也称为过顶时间。任意时刻距离r可表示为r2＝ρ2+V2(t-tρ)2＝ρ2(1+Γ2)，式中ρ为最小接近距离结合归一化磁矩的定义式，可得在接近磁异常目标过程中，归一化磁矩的大小为：
如图4所示，利用最简三磁强计的架构，利用平面内磁梯度对磁导数的近似等效，仅采用平面内等边三角形架构即能近似获取关于目标的磁张量矩阵。根据静磁场的散度和旋度都为0，近似得到磁张量矩阵：
假设归一化磁矩的测量值为M′N(i)，测量值由真实值MN(i)和加性测量噪声n(i)共同组成，即M′N(i)＝MN(i)+n(i)。归一化磁矩将/>称之为归一化磁矩的特征函数，将其与测量值M′N做完整相关运算可得：
假设初始探测时刻测点离磁异常目标较远，且由于噪声与特征函数没有相关性，则上式可近似等效成
特征函数曲线如图5所示，可知其能量主要集中于[-2,2]的时间窗口内，为便于检测，将相关性检测集中在窗口内计算可得
假设目标磁矩大小M＝2×106A·m2，过顶点距离ρ＝300m，对目标归一化磁矩信号加入分布为信噪比为-15dB的高斯白噪声进行时间窗口内的相关性检测，如图6所示。可以看到在较低信噪比条件下仍然可以无偏差的获取过顶点时刻。
定义特征时间为τ＝ρ/V，则无量纲参数从而通过获取真实的特征时间可以得到真实的过顶点距离。假定真实过顶点距离ρ0，其对应的真实特征时间为τ0。通过一组特征时间搜索获取真实特征时间，假定搜索特征时间为τ′，其与真实特征时间比值定义为∈＝τ′/τ0，则用于相关性检测的归一化磁矩特征函数为
由式(11)可知做相关性检测时，m＝0即时间窗口过顶点时相关函数得到最大值，因此在特征时间搜索过程中，比较过顶点时刻不同特征时间对应的相关函数大小。定义新的特征函数则不同∈对应的过顶点时刻相关函数(考虑目标函数连续情况)为
令则相关函数r(∈)随∈的变化关系等同于rCh(∈)随∈的变化关系。图7为归一化rCh(∈)随∈的变化关系，可以看出rCh(∈)在∈＝1即搜索特征时间τ′等于真实特征时间τ0时取得最大值，且具有较好的分辨能力。
本发明对测量平台姿态稳定性不敏感，通过相关性检测有效提高了被探信号的信噪比，通过特征时间搜索方法获取过顶点距离可以实现对磁目标动态探测过程中的单点实时定位。
Claims (7)
1.一种基于归一化磁矩的水下目标相关性检测定位方法，运用三轴矢量磁强计构建最简张量测量磁探仪进行磁目标动态定位，其特征在于，具体实现步骤如下：
步骤1、选取磁探仪构型，获取各磁强计三轴测量分量；
步骤2、利用磁强计三轴测量分量，求取磁张量矩阵；
步骤3、求解磁张量矩阵的三个特征值，得到动态磁探测测线上每一点的归一化磁矩；
步骤4、在航磁探测坐标系下，归一化磁矩表示成：其中，μ为相对磁导率，ρ为最小接近距离，/>tρ为磁测仪以速度V沿测量航迹到过靶点N的飞行时间；为归一化磁矩的特征函数；
步骤5、将归一化磁矩在设定的时间窗口内进行相关性检测得到检测函数，由检测函数的最大值时刻可获取航磁探测过目标顶点时刻；
步骤6、定义特征时间为τ＝ρ/V，通过搜索获取真实的特征时间得到真实的过顶点距离；取搜索特征时间为τs，不同特征时间对应求解检测函数r(τs)，过顶点时刻检测函数最大值对应真实特征时间，从而得到真实的过顶点距离。
2.根据权利要求1所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，利用磁强计三轴测量分量Bix，Biy，Biz，i＝1,2,3，求取磁张量矩阵：
其中，Txy＝(B3x-L2x)/L，Tyy＝(B3y-B2y)/L，Tzy＝(B3z-B2z)/L，Txz＝Tzx，Tyx＝Txy，Tyz＝Tzy，Tzz＝-Txx-Tyy，L为梯度仪基线。
3.根据权利要求1所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，求解磁张量矩阵的三个特征值分别为λ1，λ2和λ3，满足λ1≥λ2≥λ3，得到动态磁探测测线上每一点的归一化磁矩为
4.根据权利要求1所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，在航磁探测坐标系下，设磁偶极子目标位于坐标原点，航磁测量航迹与X轴方向平行；定义N(0,yρ,zρ)为过顶点，即ON垂直于测量直线航迹，航迹方程为x＝V(t-tρ)，y＝yρ，z＝zρ。
5.根据权利要求1所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，步骤5将归一化磁矩在[-2,2]的时间窗口内进行相关性检测得到检测函数为：
其中，m为检测时间序列，由检测函数的最大值时刻可获取航磁探测过目标顶点时刻。
6.根据权利要求5所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，假设归一化磁矩的测量值为M′N(i)，测量值由真实值MN(i)和加性测量噪声n(i)共同组成，即M′N(i)＝MN(i)+n(i)；归一化磁矩将/>称为归一化磁矩的特征函数，将其与测量值M′N做完整相关运算可得：
将上式近似等效成
将相关性检测集中在[-2,2]的时间窗口内计算可得
7.根据权利要求1所述的基于归一化磁矩的水下目标相关性检测定位方法，其特征在于，假定真实过顶点距离ρ0，其对应的真实特征时间为τ0；通过一组特征时间搜索获取真实特征时间，假定搜索特征时间为τ′，其与真实特征时间比值定义为∈＝τ′/τ0，则用于相关性检测的归一化磁矩特征函数为
定义新的特征函数则不同∈对应的过顶点时刻相关函数为：
令则相关函数r(∈)随∈的变化关系等同于rCh(∈)随∈的变化关系。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011081070.6A CN112344925B (zh) | 2020-10-11 | 2020-10-11 | 基于归一化磁矩的水下目标相关性检测定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011081070.6A CN112344925B (zh) | 2020-10-11 | 2020-10-11 | 基于归一化磁矩的水下目标相关性检测定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112344925A CN112344925A (zh) | 2021-02-09 |
| CN112344925B true CN112344925B (zh) | 2023-11-21 |
Family
ID=74361691
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011081070.6A Active CN112344925B (zh) | 2020-10-11 | 2020-10-11 | 基于归一化磁矩的水下目标相关性检测定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112344925B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113238289A (zh) * | 2021-03-18 | 2021-08-10 | 北京军懋国兴科技股份有限公司 | 一种基于检测依据的磁探仪探测宽度预报模型及求解方法 |
| CN114440870B (zh) * | 2021-12-29 | 2023-06-27 | 宜昌测试技术研究所 | 一种基于自适应互补滤波的九轴磁罗盘姿态获取方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101796433A (zh) * | 2007-06-27 | 2010-08-04 | 普拉德研究及开发股份有限公司 | 在电磁成像勘探中用于除去导电套管以及井眼和地面多相性的影响的方法和系统 |
| CN108931241A (zh) * | 2018-03-22 | 2018-12-04 | 南京理工大学 | 基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法 |
| JP6739130B1 (ja) * | 2019-09-30 | 2020-08-12 | フジデノロ株式会社 | 磁気検出装置、検出方法、及び検出プログラム |
-
2020
- 2020-10-11 CN CN202011081070.6A patent/CN112344925B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101796433A (zh) * | 2007-06-27 | 2010-08-04 | 普拉德研究及开发股份有限公司 | 在电磁成像勘探中用于除去导电套管以及井眼和地面多相性的影响的方法和系统 |
| CN108931241A (zh) * | 2018-03-22 | 2018-12-04 | 南京理工大学 | 基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法 |
| JP6739130B1 (ja) * | 2019-09-30 | 2020-08-12 | フジデノロ株式会社 | 磁気検出装置、検出方法、及び検出プログラム |
Non-Patent Citations (4)
| Title |
|---|
| Magnetic Anomaly Detection and Localization Using Orthogonal Basis of Magnetic Tensor Contraction;Huanghuang Jin 等;《IEEE Transactions on Geoscience and Remote Sensing》;全文 * |
| 基于磁体模拟法的潜艇磁矩分离数学模型;张维 等;《舰船科学技术》(第5期);全文 * |
| 基于磁异常的目标体定位反演方法研究;吴国超;《中国优秀硕士论文全文数据库》(第8期);全文 * |
| 水中非合作运动磁性目标跟踪及参数估计;高俊吉;刘大明;周国华;;哈尔滨工程大学学报(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112344925A (zh) | 2021-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Clark | Corrigendum to: New methods for interpretation of magnetic vector and gradient tensor data I: eigenvector analysis and the normalised source strength | |
| CN112344925B (zh) | 基于归一化磁矩的水下目标相关性检测定位方法 | |
| CN108931241B (zh) | 基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法 | |
| Jin et al. | Magnetic anomaly detection and localization using orthogonal basis of magnetic tensor contraction | |
| Huang et al. | Characterization of UXO-like targets using broadband electromagnetic induction sensors | |
| CN106990440A (zh) | 一种基于两个探测位置处磁场空间梯度的潜艇定位方法 | |
| CN111399066B (zh) | 一种基于正交基函数处理标量磁异常梯度信号的方法 | |
| Gang et al. | Detection, localization and classification of multiple dipole-like magnetic sources using magnetic gradient tensor data | |
| CN109425906A (zh) | 一种磁异常探测矢量磁目标识别方法 | |
| Li et al. | Preferred configuration and detection limits estimation of magnetic gradient tensor system | |
| Liu et al. | Magneto-inductive magnetic gradient tensor system for detection of ferromagnetic objects | |
| Teixeira et al. | Magnetic navigation and tracking of underwater vehicles | |
| Wang et al. | Frequency characteristics analysis for magnetic anomaly detection | |
| Li et al. | Application of Helbig integrals to magnetic gradient tensor multi-target detection | |
| Gang et al. | Integrated calibration of magnetic gradient tensor system | |
| Zhou et al. | Magnetic anomaly detection via a combination approach of minimum entropy and gradient orthogonal functions | |
| Shen et al. | Geometry structure optimization of hexagonal pyramidal full tensor magnetic gradient probe | |
| Wang et al. | Joint use of scalar and vector magnetometers for underwater magnetic target localization | |
| Narkhov et al. | The perspective of full-gradient magnetometry system application in searching for magnetic objects | |
| Liu et al. | Detection and localization of improvised explosive devices based on 3-axis magnetic sensor array system | |
| Mersch et al. | Underwater magnetic target localization and characterization using a three-axis gradiometer | |
| CN111596369A (zh) | 面向埋地未爆弹的磁场梯度探测仪 | |
| Tingyu et al. | Noise suppression based on temporal coherent and it’s application to magnetic anomaly detection | |
| Wang et al. | A new data processing method for magnetic anomaly detection and localization based on 2-D orthonormal basis functions | |
| Clark | Interpretation of the magnetic gradient tensor and normalized source strength applied to the Tallawang magnetite skarn deposit, New South Wales, Australia |
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 |