CN108061922B - 一种分布式磁异常探测***动态磁补偿方法 - Google Patents
一种分布式磁异常探测***动态磁补偿方法 Download PDFInfo
- Publication number
- CN108061922B CN108061922B CN201610975122.1A CN201610975122A CN108061922B CN 108061922 B CN108061922 B CN 108061922B CN 201610975122 A CN201610975122 A CN 201610975122A CN 108061922 B CN108061922 B CN 108061922B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- magnetic
- carrier
- compensation
- estimation
- 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 36
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims description 7
- 230000006698 induction Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 2
- 230000009471 action Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Manufacturing & Machinery (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明属于磁异常探测领域,具体公开了一种分布式磁异常探测***动态磁补偿方法,建立如下优化的载体磁场补偿模型,采用二级估计法确定x的最终参数估计值和独立同分布随机子样g(t)的非参数估计值
Description
技术领域
本发明属于磁异常探测领域,具体涉及一种分布式磁异常探测***动态磁补偿方法。
背景技术
磁异常探测***可实现对潜艇的探测、识别、跟踪与定位,具有潜艇定位精度高、全天候应用等突出特点,在复杂海况或复杂水声环境下进行磁异常探测、精确定位磁异常目标等方面相比传统声呐探潜具有不可替代的优势,已经成为世界各军事强国目前大力发展的探潜关键技术。围绕磁异常探潜的重大战略需求,提出了基于新型原子磁强计的分布式磁异常探测***,形成无人机编队、水面组网、水下阵列等分布式磁异常探测的新型探测模式,构建空中、水面、水下三位一体的全天时、全天候智能磁异常探测网络,大幅提升在复杂对抗环境下对磁异常目标的探测、识别、跟踪与定位能力,代表了新一代磁异常探测技术的发展方向。
分布式动态磁补偿技术是分布式磁异常探测***关键部分,影响***的综合探测效能。分布式磁异常探测***实现远距离、高效率的磁异常目标探测,对分布式动态磁补偿技术提出了更高的要求。磁干扰源包括与载体有关的干扰和与载体无关的干扰,本文研究的是与载体有关的磁干扰补偿方法。传统的载体干扰磁场模型主要有恒定干扰磁场、感应干扰磁场和涡流干扰磁场。但载体实际飞行过程中,除了上述载体干扰磁场之外,还会存在其他因素可能对磁强计测量产生影响,因此传统载体干扰磁场模型不能全面描述真实的物理过程,导致测量结果存在较大的误差,进而影响磁异常探测***的性能。
发明内容
本发明的目的是提供一种分布式磁异常探测***动态磁补偿方法,其能够提高载体磁补偿的精度,改善补偿效率。
本发明的技术方案如下:
一种分布式磁异常探测***动态磁补偿方法,该方法包括如下步骤:
步骤一、建立如下优化的载体磁场补偿模型,以确定载体磁补偿后的磁场强度HL
HL=f(x)+g(t)+Δ
Hm=Ho+Xp+XiHo+XeHo
其中,x=(xp,xi,…)为载体恒定、感应和涡流干扰磁场模型参数;g(t)为非参数独立同分布随机子样,t∈[0,1];Δ为随机误差;Hm为磁强计测量值矢量,Ho为传统模型中的地磁场真实值矢量,Xp为载体恒定干扰磁场矢量,Xi为载体感应干扰磁场参数矩阵,Xe为载体涡流干扰磁场参数矩阵;
步骤二、采用二级估计法确定x的最终参数估计值和独立同分布随机子样g(t)的非参数估计值具体为
1)对载体恒定、感应和涡流干扰磁场模型参数x进行一级估计,得到一级估计结果x*;
2)利用独立同分布随机子样g(t)进行一级估计,得到一级估计结果g*(t);
3)对上述的一级估计结果x*,进行二级估计,得到最终参数估计值
4)对上述的一级估计结果g*(t),进行二级估计,得到最终非参数估计值
5)带回优化的载体磁场补偿模型解算得出载体磁补偿后的磁场强度HL;
步骤三、对载体磁干扰进行动态补偿,具体为
1)确定动态载体的飞行轨迹
2)确定空间梯度引起的磁场变化
3)空间梯度引起的磁场变化对步骤一中的模型解算出的磁强计测量值矢量Hm进行补偿,得到补偿后的磁强计测量值矢量H′m。
在上述的一种分布式磁异常探测***动态磁补偿方法中:所述的步骤三,利用下式确定空间梯度引起的磁场变化
其中,kx,ky,kz为北向磁场、东向磁场以及垂向磁场的梯度常数,Δx,Δy,Δz为纬度、经度以及高度的变化量。
在上述的一种分布式磁异常探测***动态磁补偿方法中:所述的步骤三中,利用下式得到补偿后的磁强计测量值矢量H′m
H′m=Hm+ΔHg
然后利用下式解算载体磁补偿后的磁场强度HL
在上述的一种分布式磁异常探测***动态磁补偿方法中:所述的步骤二中对载体恒定、感应和涡流干扰磁场模型参数x进行一级估计和二级估计均采用最小二乘法。
本发明的显著效果如下:基于传统的磁场干扰模型设计了的优化的载体磁场补偿模型,以便未来求解载体磁场补偿后的磁场强度HL,通过二级估计法进行参数估计和非参数估计,带回到模型中解算得出载体磁干扰补偿后的磁场强度HL,之后利用空间梯度引起的磁场变化,对磁强计测量值矢量进行动态补偿,利用非参数项来表示测量中可能存在的影响因素,对于测量值得描述更加充分,一方面使载体干扰磁场模型与实际更加符合、提高了模型的准确度;另外一方面也将提高参数估计的精度、减小测量的误差;同时,实时高精度动态补偿提高了分布式磁异常探潜***的性能。
具体实施方式
下面结合具体实施方式对本发明作进一步说明。
步骤一、建立优化的载体磁场补偿模型。
传统的磁场干扰模型为:
Hm=Ho+Xp+XiHo+XeHo (1)
其中,Hm为磁强计测量值矢量,Ho为地磁场真实值矢量,Xp为载体恒定干扰磁场矢量,Xi为载体感应干扰磁场参数矩阵,Xe为载体涡流干扰磁场参数矩阵。
真实磁场为:
令
其中,x=(xp,xi,…)为载体恒定、感应和涡流干扰磁场模型参数。
考虑到实际测量中存在的非模型参数误差,引入半参数模型表示这一部分误差,优化的载体磁场补偿模型为:
HL=f(x)+g(t)+Δ (4)
其中,g(t)为非参数独立同分布(independent identically distributed简称idd)随机子样,t∈[0,1]。Δ为随机误差。
步骤二、采用二级估计法对式(4)进行参数估计和非参数估计。
首先对模型参数x和非参数g(t)进行一级估计,然后采用一级估计结果x*和g*(t)进行二级估计,获得最终的参数估计值和非参数估计值
模型参数x和g(t)的一级估计如下:
设α=E(g(ti)),E(g2(ti))<∞,记ei=g(ti)-α-Δi,i=1,2,…,
则E(ei)=0,E(ei 2)=δε 2=var(g(ti))+δ2<∞
则式(4)可写成:
HL=f(x)+αJ+e (5)
采用非线性最小二乘法(递推法)可得到x和α的一级估计值x*和α*。
其中,e=(e1,……eN)′,J=(1,……1),N为测量值序列的长度。
将x*代入式(4)可得
HLi-fi(x*)=gi(t)+Δi (6)
定义权函数
其中,为核函数,B=0.01N1/5为窗宽。
模型参数x和g(t)的二级估计如下:
将式(7)代入式(4),求下式的极小值:
上式求得的x的估计值为参数量的估计值,非参数量的估计为:
步骤三、动态载体干扰补偿
载体干扰磁场补偿分为两步,第一步进行补偿学习飞行,以获取相关参数;第二步在磁异常探测飞行中进行动态补偿。在磁干扰补偿学习飞行过程中,获得磁场测量值后,首先进行磁场梯度补偿,然后由步骤二获取载体干扰的相关系数。获取载体干扰参数成功后,实施动态补偿,完成磁异常探测任务。
补偿学习飞行在四个航向上进行,分别为东、南、西、北航向,且在每个航向上进行幅度为±10°的俯仰、摇摆、偏航三种机动动作。每个方向上三种动作的顺序要一致。每个动作做3个~4个周期,时间约30s,并且幅值和周期大小要相同。
地磁场的磁感应强度Hg随纬度、经度以及高度的变化而变化。在进行磁干扰补偿系数求解时,导致最后求解得到的值包含了无人机本体磁干扰以及空间的磁场变化,从而影响补偿系数的精度。分布式磁异常探测***的精度要求高,空间梯度带来的补偿系数求解误差需要补偿。在补偿飞行阶段,假设北向磁场、东向磁场以及垂向磁场的梯度kx,ky,kz为常数,若假设其纬度、经度以及高度的变化量为Δx,Δy,Δz,则空间梯度引起的磁场变化:
磁强计的测量值首先进行磁场梯度补偿,即载体磁场干扰为:
H′m=Hm+ΔHg (11)
则优化的载体磁场补偿模型修正为:
HL=f′(x)+g(t)+Δ (12)
由补偿飞行获得测量值后,根据二级估计法对参数进行估计。
磁异常探测飞行中,根据磁强计的测量值及补偿飞行估计的参数值实时解算式(12)获取真实磁场,进而完成实时完成磁异常探测任务。
Claims (2)
1.一种分布式磁异常探测***动态磁补偿方法,其特征在于,该方法包括如下步骤:
步骤一、建立如下优化的载体磁场补偿模型,以确定载体磁补偿后的磁场强度HL
HL=f(x)+g(t)+Δ
Hm=Ho+Xp+XiHo+XeHo
其中,x=(xp,xi,…)为载体恒定、感应和涡流干扰磁场模型参数;g(t)为非参数独立同分布随机子样,t∈[0,1];Δ为随机误差;Hm为磁强计测量值矢量,Ho为传统模型中的地磁场真实值矢量,Xp为载体恒定干扰磁场矢量,Xi为载体感应干扰磁场参数矩阵,Xe为载体涡流干扰磁场参数矩阵;
步骤二、采用二级估计法确定x的最终参数估计值和独立同分布随机子样g(t)的非参数估计值具体为
1)对载体恒定、感应和涡流干扰磁场模型参数x进行一级估计,得到一级估计结果x*;
2)利用独立同分布随机子样g(t)进行一级估计,得到一级估计结果g*(t);
3)对上述的一级估计结果x*,进行二级估计,得到最终参数估计值
4)对上述的一级估计结果g*(t),进行二级估计,得到最终非参数估计值
5)带回优化的载体磁场补偿模型解算得出载体磁补偿后的磁场强度HL;
步骤三、对载体磁干扰进行动态补偿,具体为
1)确定动态载体的飞行轨迹
2)确定空间梯度引起的磁场变化
在补偿飞行阶段,假设北向磁场、东向磁场以及垂向磁场的梯度kx,ky,kz为常数,若假设其纬度、经度以及高度的变化量为Δx,Δy,Δz,则空间梯度引起的磁场变化为
3)空间梯度引起的磁场变化对步骤一中的模型解算出的磁强计测量值矢量Hm进行补偿,得到补偿后的磁强计测量值矢量H′m;
磁强计的测量值首先进行磁场梯度补偿,即载体磁场干扰为:
Hm′=Hm+ΔHg
优化的载体磁场补偿模型修正为:
HL=f′(x)+g(t)+Δ。
2.如权利要求1所述的一种分布式磁异常探测***动态磁补偿方法,其特征在于:所述的步骤二中对载体恒定、感应和涡流干扰磁场模型参数x进行一级估计和二级估计均采用最小二乘法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610975122.1A CN108061922B (zh) | 2016-11-07 | 2016-11-07 | 一种分布式磁异常探测***动态磁补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610975122.1A CN108061922B (zh) | 2016-11-07 | 2016-11-07 | 一种分布式磁异常探测***动态磁补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108061922A CN108061922A (zh) | 2018-05-22 |
CN108061922B true CN108061922B (zh) | 2019-06-11 |
Family
ID=62136520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610975122.1A Active CN108061922B (zh) | 2016-11-07 | 2016-11-07 | 一种分布式磁异常探测***动态磁补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108061922B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109283585B (zh) * | 2018-10-15 | 2024-06-07 | 南京理工大学 | 一种采用滑轨的磁异常目标探测仪 |
CN111060980A (zh) * | 2018-10-17 | 2020-04-24 | 北京自动化控制设备研究所 | 一种分布式海浪磁干扰补偿方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236108A (zh) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | 一种磁性地表三维地形改正方法 |
CN102252689A (zh) * | 2010-05-19 | 2011-11-23 | 北京国浩传感器技术研究院(普通合伙) | 一种基于磁传感器的电子罗盘校准方法 |
CN102927984A (zh) * | 2012-10-26 | 2013-02-13 | 哈尔滨工程大学 | 消除载体磁化磁场对地磁测量影响的方法 |
CN103513284A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | 一种双磁性层磁异常的剥离方法 |
CN103577607A (zh) * | 2013-11-20 | 2014-02-12 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
CN103592691A (zh) * | 2013-10-30 | 2014-02-19 | 中国石油天然气集团公司 | 一种反向剩余磁性磁异常识别方法及装置 |
CN103837845A (zh) * | 2014-01-22 | 2014-06-04 | 哈尔滨工程大学 | 一种飞行器磁干扰场模型参数解算方法 |
WO2015013499A1 (en) * | 2013-07-24 | 2015-01-29 | Schlumberger Canada Limited | Method to predict local geomagnetic disturbance field and its practical application |
CN105509737A (zh) * | 2015-11-26 | 2016-04-20 | 哈尔滨工业大学 | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 |
-
2016
- 2016-11-07 CN CN201610975122.1A patent/CN108061922B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236108A (zh) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | 一种磁性地表三维地形改正方法 |
CN102252689A (zh) * | 2010-05-19 | 2011-11-23 | 北京国浩传感器技术研究院(普通合伙) | 一种基于磁传感器的电子罗盘校准方法 |
CN103513284A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | 一种双磁性层磁异常的剥离方法 |
CN102927984A (zh) * | 2012-10-26 | 2013-02-13 | 哈尔滨工程大学 | 消除载体磁化磁场对地磁测量影响的方法 |
WO2015013499A1 (en) * | 2013-07-24 | 2015-01-29 | Schlumberger Canada Limited | Method to predict local geomagnetic disturbance field and its practical application |
CN103592691A (zh) * | 2013-10-30 | 2014-02-19 | 中国石油天然气集团公司 | 一种反向剩余磁性磁异常识别方法及装置 |
CN103577607A (zh) * | 2013-11-20 | 2014-02-12 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
CN103837845A (zh) * | 2014-01-22 | 2014-06-04 | 哈尔滨工程大学 | 一种飞行器磁干扰场模型参数解算方法 |
CN105509737A (zh) * | 2015-11-26 | 2016-04-20 | 哈尔滨工业大学 | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 |
Non-Patent Citations (1)
Title |
---|
"基于磁异常反演的磁航向误差实时补偿方法";尹刚 等;《武汉大学学报(信息科学版)》;20160731;第41卷(第7期);第978-982页 |
Also Published As
Publication number | Publication date |
---|---|
CN108061922A (zh) | 2018-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107272069B (zh) | 基于磁异常梯度的磁性目标追踪方法 | |
CN103630137B (zh) | 一种用于导航***的姿态及航向角的校正方法 | |
CN107894235B (zh) | 一种超高速飞行器自主导航***的模型误差补偿方法 | |
Li et al. | Magnetic sensors for navigation applications: an overview | |
CN108844536B (zh) | 一种基于量测噪声协方差矩阵估计的地磁导航方法 | |
CN104406610B (zh) | 一种磁力计实时校正装置及方法 | |
CN102445201A (zh) | 用于水下载体的地磁异常特征点匹配导航方法 | |
CN108827310A (zh) | 一种船用星敏感器辅助陀螺仪在线标定方法 | |
CN109556631A (zh) | 一种基于最小二乘的ins/gnss/偏振/地磁组合导航***对准方法 | |
Chen et al. | A new geomagnetic matching navigation method based on multidimensional vector elements of earth’s magnetic field | |
CN103954303B (zh) | 一种用于磁力计导航***航向角动态计算及校正方法 | |
CN104535080B (zh) | 大方位失准角下基于误差四元数的传递对准方法 | |
CN109725360A (zh) | 基于磁梯度张量不变量的单点定位方法 | |
CN202209953U (zh) | 用于水下载体的地磁辅助惯性导航*** | |
CN107132542A (zh) | 一种基于光学和多普勒雷达的小天体软着陆自主导航方法 | |
CN109521384A (zh) | 一种基于原子磁强计的矢量磁补偿方法 | |
CN107270898A (zh) | 基于mems传感器和vlc定位融合的双粒子滤波导航装置和方法 | |
CN110736484B (zh) | 基于陀螺仪及磁传感器融合的背景磁场标定方法 | |
CN108061922B (zh) | 一种分布式磁异常探测***动态磁补偿方法 | |
Huang et al. | Method of separating dipole magnetic anomaly from geomagnetic field and application in underwater vehicle localization | |
CN106918350B (zh) | 一种应用于地磁导航中的地磁场模型误差补偿方法 | |
Chen et al. | An improved geomagnetic navigation method based on two-component gradient weighting | |
CN111220932A (zh) | 无人机磁干扰标定方法及分布式磁异常探测*** | |
Ning et al. | Improved MEMS magnetometer adaptive filter noise reduction and compensation method | |
Wang et al. | An adaptive Kalman filtering algorithm based on maximum likelihood estimation |
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 |