CN108919368A - 一种用于消除微小卫星剩磁干扰的系统及方法 - Google Patents
一种用于消除微小卫星剩磁干扰的系统及方法 Download PDFInfo
- Publication number
- CN108919368A CN108919368A CN201810291494.1A CN201810291494A CN108919368A CN 108919368 A CN108919368 A CN 108919368A CN 201810291494 A CN201810291494 A CN 201810291494A CN 108919368 A CN108919368 A CN 108919368A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- vector
- satellite platform
- microsatellite
- probe
- 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
Links
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/40—Electric 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
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明提供了一种用于消除微小卫星剩磁干扰的系统及方法,将微小卫星平台剩磁建模为磁偶极子,在卫星表面两对称伸杆的末端处各布设一个矢量磁力仪探头,利用磁场梯度积分与磁场值之间的内在联系,估算卫星平台的剩磁干扰,从而给出卫星本体位置处的地磁场值。本发明采用磁场梯度积分的方法给出卫星平台剩磁估计方法,可达到消除卫星平台剩磁干扰,提高空间地磁场测量准确性的目的。
Description
技术领域
本发明属于地球物理卫星磁法勘探领域,涉及一种基于磁场梯度测量消除微小卫星平台磁干扰的系统及方法。
背景技术
卫星磁测作为一种重要的地磁场探测方法在地球物理卫星磁法勘探领域获得广泛的应用。现有卫星磁测法将磁力仪固定在卫星伸杆上,通过增加磁力仪传感部分与卫星平台本体的距离来减小卫星本体剩磁对磁力仪测量结果的影响,随着磁力仪指标的提高,仅依靠这种方法已经无法有效消除卫星平台本体剩磁的干扰,因此,采用有效方法计算并消除卫星平台造成的磁干扰具有重要的意义。
传统的卫星平台磁干扰计算方法将两个磁力仪探头布置在同一伸杆上,其中一个探头设置在伸杆末端,另一个磁力仪探头设置在伸杆末端和卫星本体之间某一位置处,通过双探头差分测量方法计算卫星平台的剩磁。但微小卫星相对于普通卫星来说磁矩更小,采用双探头差分测量方法估计微小卫星平台剩磁的精度无法满足需求,而且传统方法采用的伸杆的长度较长,磁力仪探头的安装方式不适合微小卫星。
发明内容
本发明解决的技术问题是:提供了一种基于磁场梯度测量消除微小卫星平台磁干扰的系统及方法,在卫星上布设两个伸杆,保证卫星中心与两伸杆共线,在两伸杆末端各布设一个矢量磁力仪探头,结合消除磁干扰的算法解决微小卫星平台带来的磁干扰问题。
本发明的技术方案是:
一种用于消除微小卫星剩磁干扰的系统,包括第一矢量磁力仪探头、第二矢量磁力仪探头、磁力仪电子箱、第一线缆、第二线缆、第一伸杆、第二伸杆和卫星平台;
第一伸杆和第二伸杆分别安装在卫星平台的两侧,且卫星平台的中心与第一伸杆和第二伸杆共线,第一矢量磁力仪探头与第二矢量磁力仪探头分别固定在第一伸杆末端以及第二伸杆的末端,磁力仪电子箱设置在卫星平台中心处,第一线缆和第二线缆用于实现第一矢量磁力仪探头和第二矢量磁力仪探头与磁力仪电子箱之间的信号传输;
磁力仪电子箱对第一矢量磁力仪探头和第二矢量磁力仪探头测量得到的包含磁干扰的磁场值进行处理,从而消除微小卫星的剩磁干扰,得到真实地磁场值。
所述第一伸杆和第二伸杆关于卫星平台的中心对称放置,且第一伸杆和第二伸杆长度相同。
第一矢量磁力仪探头和第二矢量磁力仪探头与卫星平台中心之间的距离相同。
第一矢量磁力仪探头、第二矢量磁力仪探头是磁敏感单元。
卫星平台为微小卫星。
一种消除微小卫星剩磁干扰的方法,步骤如下:
(1)记录第一矢量磁力仪探头(101)和第二矢量磁力仪探头(105)测量得到的磁场值,分别为和
(2)根据所述磁场值和计算与卫星平台(103)剩磁相关的三个磁矩参数k1、k2和k3;
其中L1是第一矢量磁力仪探头(101)在卫星平台坐标系X轴的投影,其为负,L2是第二矢量磁力仪探头(105)在卫星平台坐标系X轴的投影,其为正;
为磁场的三个分量,即
为磁场的三个分量,即
所述卫星平台坐标系的原点为卫星平台(103)中心,以伸杆的轴向为X轴,第二伸杆所在方向为X轴正方向,卫星平台(103)指向地心的方向为Z轴正向,Y轴由右手定则确定;
(3)利用磁矩参数k1、k2和k3计算第一矢量磁力仪探头(101)处的卫星平台(103)的磁干扰
(4)利用磁矩参数k1、k2和k3计算第二矢量磁力仪探头(105)处的卫星平台(103)的磁干扰
(5)将第一矢量磁力仪探头(101)的磁场值与卫星平台(103)在第一矢量磁力仪探头(101)处的磁干扰作差,计算第一矢量磁力仪探头(101)处地磁场
(6)将第二矢量磁力仪探头(105)的磁场值与卫星平台(103)在第二矢量磁力仪探头(105)处的磁干扰作差,计算第二矢量磁力仪探头(105)处地磁场
(7)计算第一矢量磁力仪探头(101)处的背景磁场与第二矢量磁力仪探头(105)处的背景磁场的平均值,得到空间地磁场值从而完成微小卫星剩磁干扰的消除;
本发明与现有技术相比的有益效果:
(1)本发明只需在微小卫星两侧安装矢量磁力仪探头,从而消除卫星本体引起的磁场干扰。
(2)本发明采用磁场梯度积分的方法计算卫星本体剩磁,能够提高磁场测量精度。
(3)本发明的系统结构简单,可用于微小卫星平台。
(4)本发明采用矢量磁力仪,可得到空间磁场矢量信息。
(5)本发明算法计算量小,占用资源少,能够实现地磁场在轨实时校准。
附图说明
图1是本发明的消除微小卫星剩磁干扰系统结构示意图;
图2是本发明中L1、L2、Ln参数的说明示意图;
图3是本发明的微小卫星平台磁干扰消除流程图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步说明,但不作为对本发明的限定。
如图1所示,本发明给出的消除微小卫星剩磁干扰系统结构示意图,包括第一矢量磁力仪探头101、第二矢量磁力仪探头105、磁力仪电子箱106、第一线缆107、第二线缆108、第一伸杆102、第二伸杆104和卫星平台103;
第一伸杆102和第二伸杆104分别安装在卫星平台103的两侧,且卫星平台103的中心与第一伸杆102和第二伸杆104共线,如图2所示,第一矢量磁力仪探头101与第二矢量磁力仪探头105分别固定在第一伸杆102末端以及第二伸杆104的末端,磁力仪电子箱106设置在卫星平台103中心处,第一线缆107和第二线缆108用于实现第一矢量磁力仪探头101和第二矢量磁力仪探头105与磁力仪电子箱106之间的信号传输;
磁力仪电子箱106对第一矢量磁力仪探头101和第二矢量磁力仪探头105测量得到的包含磁干扰的磁场值进行处理,从而消除微小卫星的剩磁干扰,得到真实地磁场值。
所述第一伸杆102和第二伸杆104关于卫星平台103的中心对称放置,且第一伸杆102和第二伸杆104长度相同,本发明设置为5m。
第一矢量磁力仪探头101和第二矢量磁力仪探头105与卫星平台103中心之间的距离相同。第一矢量磁力仪探头101、第二矢量磁力仪探头105是磁敏感单元。本发明中的卫星平台103指微小卫星。
本发明基于上述系统,还提供了一种用于消除微小卫星剩磁干扰的方法,在卫星平台上对称布设两伸杆101和102,在两伸杆末端分别固定矢量磁力仪探头101和105,将微小卫星平台103剩磁建模为磁偶极子,利用梯度积分为磁场的原理计算卫星平台103剩磁,达到准确求解卫星平台本体处地磁场、计算两伸杆处地磁梯度的目的。
所依据原理为:(1)磁场梯度是磁场强度在磁场方向上单位距离的变化量且磁场强度与距离的三次方成反比,故磁场梯度与距离的四次方成反比。(2)无穷远处微小卫星平台产生的磁干扰为零。(3)在轨道高度处,磁力仪单次测量过程中,直径10米区域内地磁场值为均一、恒定值。
本发明提出的消除剩磁干扰实施例如图3所示。
步骤301:磁测系统各部件安装好以后,开始进入卫星平台磁干扰计算阶段;
步骤302:读取第一矢量磁力仪探头101数据
第一矢量磁力仪探头101所在位置处矢量磁场表达式为
式(1)中为第一矢量磁力仪探头测量得到的地磁场值,卫星平台103在第一矢量磁力仪探头101处产生干扰可表达为式(2)。
其中gx(r)、gy(r)、gz(r)为第一矢量磁力仪探头101所在位置处矢量磁场所在三个轴的磁场梯度,k1、k2和k3为与卫星平台103剩磁有关的参数,以卫星平台103中心为原点建立卫星平台坐标系,L1是第一矢量磁力仪探头101在卫星平台坐标系X轴的投影,其为负,Ln则表示正无穷远点+∞,这里注意到磁场梯度与距离的四次方成反比。
所述卫星平台坐标系的原点为卫星平台103中心,以伸杆的轴向为X轴,第二伸杆所在方向为X轴正方向,卫星平台103指向地心的方向为Z轴正向,Y轴由右手定则确定;
步骤303:读取第二矢量磁力仪探头105数据
第二矢量磁力仪探头105所在位置处矢量磁场表达式为
式(3)中卫星平台103在第二矢量磁力仪探头105处产生干扰可表达为式(4)。
其中L2是第二矢量磁力仪探头105在卫星平台坐标系X轴的投影,其为正;
步骤304:第一矢量磁力仪探头101与第二矢量磁力仪探头105的测量结果作差,表达为式(5)。
式(5)利用积分区间可加性,同时等式左侧的测量值已知,可直接求解与卫星剩磁相关的三个参数:
其中为磁场的三个分量,即
为磁场的三个分量,即
步骤305:利用计算所得参数k1、k2和k3,求出卫星平台103在第一矢量磁力仪探头101所在位置处磁干扰
步骤306:计算得到第一矢量磁力仪探头测得的地磁场值为
步骤307:利用计算所得参数k1、k2和k3,求出卫星平台103在第二矢量磁力仪探头105所在位置处磁干扰
步骤308:计算得到第二矢量磁力仪探头测得的地磁场值为
步骤309:根据两个矢量磁力仪的测量结果计算平均值,即为地磁场值,故地磁场表达为式(10)。
步骤310:结束一次磁场测量过程。
可以理解的是,以上所述的实施例只是本发明较优选的具体实施方式,对本领域技术人员来说,对本发明的技术方案及发明构思加以等同替换或改变都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种用于消除微小卫星剩磁干扰的系统,其特征在于:包括第一矢量磁力仪探头(101)、第二矢量磁力仪探头(105)、磁力仪电子箱(106)、第一线缆(107)、第二线缆(108)、第一伸杆(102)、第二伸杆(104)和卫星平台(103);
第一伸杆(102)和第二伸杆(104)分别安装在卫星平台(103)的两侧,且卫星平台(103)的中心与第一伸杆(102)和第二伸杆(104)共线,第一矢量磁力仪探头(101)与第二矢量磁力仪探头(105)分别固定在第一伸杆(102)末端以及第二伸杆(104)的末端,磁力仪电子箱(106)设置在卫星平台(103)中心处,第一线缆(107)和第二线缆(108)用于实现第一矢量磁力仪探头(101)和第二矢量磁力仪探头(105)与磁力仪电子箱(106)之间的信号传输;
磁力仪电子箱(106)对第一矢量磁力仪探头(101)和第二矢量磁力仪探头(105)测量得到的包含磁干扰的磁场值进行处理,从而消除微小卫星的剩磁干扰,得到真实地磁场值。
2.根据权利要求1所述的一种用于消除微小卫星剩磁干扰的系统,其特征在于:所述第一伸杆(102)和第二伸杆(104)关于卫星平台(103)的中心对称放置,且第一伸杆(102)和第二伸杆(104)长度相同。
3.根据权利要求1所述的一种用于消除微小卫星剩磁干扰的系统,其特征在于:第一矢量磁力仪探头(101)和第二矢量磁力仪探头(105)与卫星平台(103)中心之间的距离相同。
4.根据权利要求1所述的一种用于消除微小卫星剩磁干扰的系统,其特征在于:第一矢量磁力仪探头(101)、第二矢量磁力仪探头(105)是磁敏感单元。
5.根据权利要求1~4中任一项所述的一种用于消除微小卫星剩磁干扰的系统,其特征在于:卫星平台(103)为微小卫星。
6.一种基于所述消除微小卫星剩磁干扰系统实现的消除微小卫星剩磁干扰的方法,其特征在于步骤如下:
(1)记录第一矢量磁力仪探头(101)和第二矢量磁力仪探头(105)测量得到的磁场值,分别为和
(2)根据所述磁场值和计算与卫星平台(103)剩磁相关的三个磁矩参数k1、k2和k3;
(3)利用磁矩参数k1、k2和k3计算第一矢量磁力仪探头(101)处的卫星平台(103)的磁干扰
(4)利用磁矩参数k1、k2和k3计算第二矢量磁力仪探头(105)处的卫星平台(103)的磁干扰
(5)将第一矢量磁力仪探头(101)的磁场值与卫星平台(103)在第一矢量磁力仪探头(101)处的磁干扰作差,计算第一矢量磁力仪探头(101)处地磁场
(6)将第二矢量磁力仪探头(105)的磁场值与卫星平台(103)在第二矢量磁力仪探头(105)处的磁干扰作差,计算第二矢量磁力仪探头(105)处地磁场
(7)计算第一矢量磁力仪探头(101)处的背景磁场与第二矢量磁力仪探头(105)处的背景磁场的平均值,得到空间地磁场值从而完成微小卫星剩磁干扰的消除。
7.根据权利要求6所述的一种消除微小卫星剩磁干扰的方法,其特征在于:所述步骤2中磁矩参数k1、k2和k3,具体通过如下方式计算得到:
其中L1是第一矢量磁力仪探头(101)在卫星平台坐标系X轴的投影,其为负,L2是第二矢量磁力仪探头(105)在卫星平台坐标系X轴的投影,其为正;
为磁场的三个分量,即
为磁场的三个分量,即
所述卫星平台坐标系的原点为卫星平台(103)中心,以伸杆的轴向为X轴,第二伸杆所在方向为X轴正方向,卫星平台(103)指向地心的方向为Z轴正向,Y轴由右手定则确定。
8.根据权利要求6所述的一种消除微小卫星剩磁干扰的方法,其特征在于:所述步骤(3)第一矢量磁力仪探头(101)处的卫星平台(103)的磁干扰具体为:
第二矢量磁力仪探头(105)处的卫星平台(103)的磁干扰具体为:
9.根据权利要求6所述的一种消除微小卫星剩磁干扰的方法,其特征在于:第一矢量磁力仪探头(101)处地磁场具体为:
第二矢量磁力仪探头(105)处地磁场具体为:
10.根据权利要求6所述的一种消除微小卫星剩磁干扰的方法,其特征在于:所述步骤(7)空间地磁场值具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810291494.1A CN108919368B (zh) | 2018-04-03 | 2018-04-03 | 一种用于消除微小卫星剩磁干扰的系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810291494.1A CN108919368B (zh) | 2018-04-03 | 2018-04-03 | 一种用于消除微小卫星剩磁干扰的系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919368A true CN108919368A (zh) | 2018-11-30 |
CN108919368B CN108919368B (zh) | 2020-06-09 |
Family
ID=64403487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810291494.1A Active CN108919368B (zh) | 2018-04-03 | 2018-04-03 | 一种用于消除微小卫星剩磁干扰的系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919368B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113985327A (zh) * | 2021-10-29 | 2022-01-28 | 北京卫星环境工程研究所 | 一种基于多点磁场传感的在轨磁扰动信号监测系统及方法 |
CN115561681A (zh) * | 2022-09-28 | 2023-01-03 | 北京大学 | 一种卫星磁场测量扰动修正方法及系统、存储介质及终端 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5047945A (en) * | 1990-08-09 | 1991-09-10 | General Electric Company | Residual magnetic dipole cancellation system for satellites |
CN101452063A (zh) * | 2007-12-04 | 2009-06-10 | 北京卫星环境工程研究所 | 地磁场中卫星剩磁矩的测量系统及测量方法 |
CN104090250A (zh) * | 2014-07-25 | 2014-10-08 | 中国科学院空间科学与应用研究中心 | 零磁空间中测量卫星的剩磁和感磁的装置及方法 |
CN104335716B (zh) * | 2005-07-06 | 2015-04-29 | 上海卫星装备研究所 | 地磁环境磁场中卫星磁矩测量方法 |
-
2018
- 2018-04-03 CN CN201810291494.1A patent/CN108919368B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5047945A (en) * | 1990-08-09 | 1991-09-10 | General Electric Company | Residual magnetic dipole cancellation system for satellites |
CN104335716B (zh) * | 2005-07-06 | 2015-04-29 | 上海卫星装备研究所 | 地磁环境磁场中卫星磁矩测量方法 |
CN101452063A (zh) * | 2007-12-04 | 2009-06-10 | 北京卫星环境工程研究所 | 地磁场中卫星剩磁矩的测量系统及测量方法 |
CN104090250A (zh) * | 2014-07-25 | 2014-10-08 | 中国科学院空间科学与应用研究中心 | 零磁空间中测量卫星的剩磁和感磁的装置及方法 |
Non-Patent Citations (4)
Title |
---|
周斌,等: "卫星磁部件分布对梯度法消除剩磁的影响分析", 《中国空间科学技术》 * |
张艺腾,等: "磁场梯度张量测量法消除卫星磁干扰", 《北京航空航天大学学报》 * |
易忠,等: "地磁场中测量卫星磁矩的一种新的方法", 《环模技术》 * |
翁成翰,等: "一种卫星有效载荷零部件剩磁分析方法", 《第二十六届全国空间探测学术研讨会》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113985327A (zh) * | 2021-10-29 | 2022-01-28 | 北京卫星环境工程研究所 | 一种基于多点磁场传感的在轨磁扰动信号监测系统及方法 |
CN113985327B (zh) * | 2021-10-29 | 2024-01-26 | 北京卫星环境工程研究所 | 一种基于多点磁场传感的在轨磁扰动信号监测系统及方法 |
CN115561681A (zh) * | 2022-09-28 | 2023-01-03 | 北京大学 | 一种卫星磁场测量扰动修正方法及系统、存储介质及终端 |
CN115561681B (zh) * | 2022-09-28 | 2023-05-23 | 北京大学 | 一种卫星磁场测量扰动修正方法及系统、存储介质及终端 |
Also Published As
Publication number | Publication date |
---|---|
CN108919368B (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487770B (zh) | 一种用于测量空间磁场梯度全张量信息的无源探测装置 | |
CN103837900B (zh) | 一种基于矢量磁场探测的地下电缆定位方法及装置 | |
CN109856689B (zh) | 一种超导航磁梯度张量数据抑噪处理方法和系统 | |
CN106959471B (zh) | 基于非线性航磁总场梯度补偿模型的航磁补偿方法 | |
CN109856690B (zh) | 基于混合范数拟合的航磁梯度张量数据抑噪方法及系统 | |
CN110274586A (zh) | 包含多光系原子磁力仪方向误差补偿的航空磁补偿方法 | |
CN104535062A (zh) | 基于磁梯度张量和地磁矢量测量的运动式定位方法 | |
CN108520112A (zh) | 一种基于吉洪诺夫正则化的飞机干扰磁场补偿方法 | |
CN107121655A (zh) | 非屏蔽serf原子磁力仪磁场抵消线圈非正交角测量装置及测量方法 | |
CN105785288A (zh) | 一种基于低温超导 squid 的航空磁测量装置 | |
CN108919368A (zh) | 一种用于消除微小卫星剩磁干扰的系统及方法 | |
CN102928713B (zh) | 一种磁场天线的本底噪声测量方法 | |
Chwala et al. | Full tensor SQUID gradiometer for airborne exploration | |
CN106289206A (zh) | 一种提供稳定地磁场环境的装置和方法 | |
Wan et al. | Improved component compensation for geomagnetic field vector measurement using Lagrange multiplier method | |
Mester et al. | Development and drift-analysis of a modular electromagnetic induction system for shallow ground conductivity measurements | |
CN206146426U (zh) | 一种提供稳定地磁场环境的装置 | |
CN206546425U (zh) | 一种实现标量矢量同时测量的氦光泵磁力仪探头 | |
CN102944901A (zh) | 一种大地电磁阻抗估计方法 | |
CN108254796B (zh) | 一种标量磁力仪阵列基线的优化方法 | |
RU2628448C1 (ru) | Способ измерения магнитных моментов объекта | |
Noltimier | Use of the spinner magnetometer for anisotropy measurements | |
Ripka et al. | A 3-phase current transducer based on microfluxgate sensors | |
CN112284377A (zh) | 一种应用于飞行器的地磁场测量系统及方法 | |
Felch et al. | Air-borne magnetometers for search and survey |
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 |