CN104182648B - 反演航天器内部多磁源分布的方法 - Google Patents
反演航天器内部多磁源分布的方法 Download PDFInfo
- Publication number
- CN104182648B CN104182648B CN201410446670.6A CN201410446670A CN104182648B CN 104182648 B CN104182648 B CN 104182648B CN 201410446670 A CN201410446670 A CN 201410446670A CN 104182648 B CN104182648 B CN 104182648B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- spacecraft
- magnetic
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种反演航天器内部多磁源分布的方法，该方法根据航天器产品的外包络尺寸，确定磁场梯度扫描测量平面，在扫描测量平面上划分若干等距的测试点，利用磁场梯度计测量每一个测量点位置的三方向磁感应强度以及磁场梯度值，根据各测量点磁感应强度分量计算合量并找到平面上磁感应强度的每个极值点并计算出磁源位置信息，从而得到航天器内部各磁源的位置信息。相比目前采取的将航天器简化为单偶极子的方法，无论磁场测量精度还是内部精细磁结构的分辨能力都得到了显著提高。
Description
技术领域
本发明属于航天器磁性测量技术领域，具体涉及航天器产品内部的磁性测量。
背景技术
通过外部磁场测量对物体内部磁源状况进行反演，可以获得物体内部的磁结构信息和磁场分布情况，是了解物体磁性的重要方法。目前，该类方法主要应用在矿藏探测、未爆弹探测等大空间尺度范围探测上，且多为单个反演目标。
在航天器磁性测量领域，目前还只能将航天器产品整体作为单个磁偶极子进行测量和反演，尚无有效方法对其内部磁结构进行测量。为精确掌握航天器产品内部磁性状态，北京卫星环境工程研究所研发了一种可应用在航天器这样小空间尺度内，并能实现多磁源定位的反演方法。该方法主要基于航天器产品外部磁场梯度张量测量数据，利用欧拉反摺积算法结合区域磁场极值扫描技术来实现内部多磁源扫描的目的。在本方法中，申请人首先将欧拉反演方法进行了改进，使其适应小尺度空间内的应用。其次为了解决多磁源相互干扰对目标反演的影响问题，申请人发明了平面极值搜索方法，并利用局域磁场极值反演实现了多目标的定位和反演。在该方法发明过程中，申请人设计了大量验证性实验，对方法的可行性进行了验证。
该方法应用，大大提升了航天器磁性测试能力，解决了原测试方法将整个航天器简化单磁偶极子，不能检测其内部磁结构的问题。
发明内容
本发明旨在提供一种能够通过测量航天器外部磁场状况，精确反演其内部磁结构和磁场分布的方法。
为了实现上述目的，本发明采用了如下的技术方案：
反演航天器内部多磁源分布的方法，包括如下步骤：
1)根据航天器产品的外包络尺寸，确定磁场梯度扫描测量平面，扫描测量平面为包含航天器产品的最小正六面体的外表面(六个面)；
2)在扫描测量平面上划分若干等距的测试点，利用磁场梯度计测量每一个测量点位置的三方向磁感应强度Bx,By,Bz，以及9个磁场梯度值
3)根据各测量点磁感应强度分量计算合量
并找到平面上磁感应强度的每个极值点；
4)按照航天器的三方向坐标系，将每个极值点处测到的9个磁场梯度值带入欧拉方程：
其中x0,y0,z0是各测量点坐标，计算出磁源位置信息x,y,z；
5)将上一步每个极值点处计算出的磁源位置(X1,Y1,Z1),(X2,Y2,Z2)……等记录，并剔除重复的位置点，即可得到航天器内部各磁源的位置信息。
其中，精度要求越高测试点越密集。
利用该方法，成功完成了航天器内部8个及以下磁偶极子的精确位置反演。相比目前采取的将航天器简化为单偶极子的方法，无论磁场测量精度还是内部精细磁结构的分辨能力都得到了显著提高。在航天器研制工程实践中，该方法可有效指导航天器磁性设计，并可用来诊断和定位由于航天器内部发生的磁、电异常状况。
附图说明
图1为本发明的反演航天器内部多磁源分布的方法实施的流程图。
具体实施方式
以下介绍的是作为本发明所述内容的具体实施方式，下面通过具体实施方式对本发明的所述内容作进一步的阐明。当然，描述下列具体实施方式只为示例本发明的不同方面的内容，而不应理解为限制本发明范围。
参见图1，图1为本发明的反演航天器内部多磁源分布的方法实施的流程图。以下以立方体型(3m×3m×3m)某型号航天器为例介绍该方法。
第一步，考虑航天器表面存在的突起物会于测量仪器发生干涉，扫描平面选择为距离航天器平面0.2m处，与航天器表面平行的正方形(3.2m×3.2m)平面区域。由于航天器有6个表面，扫描平面也为6个。
第二步，根据航天器测量需求，确定其内部磁源的位置反演精度为0.1m，并据此将扫描平面化分为上下、左右间距均为0.1m的测量点网格。
第三步，以航天器产品地面中心点为原点(0，0，0)，以通过原点垂直于表面的直线为X轴，建立空间直角坐标系，并确定网格内测量点坐标。
第四步，利用磁场梯度测量仪对每个测量点的三分量磁场值、磁场合量值以及9个梯度张量值进行测量并记录。
第五步，在每个扫描平面内，寻找磁场合量为极值的测量点，并记录这些点的位置和9个梯度张量信息。
第六步，将极值点处的磁场和张量测量数据带入欧拉方程，计算磁源位置信息。
第七步，记录反演出的磁源位置信息。
第八步，判断是否完成所有极值点反演。判断是否完成了所有极值点的反演计算，如没有完成择对下一个极值点重复第六到八步操作直至反演完所有的极值点。
第十步，在所有反演位置记录中查找重复点并从中去除。
第十一步，剩下的位置坐标即为被测产品内部磁源的位置坐标。
尽管上文对本发明的具体实施方式进行了详细的描述和说明，但应该指明的是，我们可以对上述实施方式进行各种改变和修改，但这些都不脱离本发明的精神和所附的权利要求所记载的范围。
Claims (2)
1.反演航天器内部多磁源分布的方法，包括如下步骤：
1)根据航天器产品的外包络尺寸，确定磁场梯度扫描测量平面，考虑航天器表面存在的突起物会与测量仪器发生干涉，扫描测量平面选择为距离航天器平面0.2m处，与航天器表面平行的正方形平面区域，其面积为3.2m×3.2m；
2)根据航天器测量需求，确定其内部磁源的位置反演精度为0.1m，并据此将扫描平面化分为上下、左右间距均为0.1m的测量点网格，以航天器产品地面中心点为原点(0，0，0)，以通过原点垂直于表面的直线为X轴，建立空间直角坐标系，并确定网格内测量点坐标，利用磁场梯度计测量每一个测量点位置的三方向磁感应强度Bx,By,Bz，以及9个磁场梯度值
3)根据各测量点磁感应强度分量计算合量
<mrow>
<mi>B</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>B</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>B</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>B</mi>
<mi>z</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
在每个扫描平面内，寻找磁场合量为极值的测量点，并记录这些点的位置和9个梯度张量信息；
4)按照航天器的三方向坐标系，将每个极值点处测到的9个磁场梯度值带入欧拉方程：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>B</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中x0,y0,z0是各测量点坐标，计算出磁源位置信息x,y,z；
5)将上一步每个极值点处计算出的磁源位置进行记录，判断是否完成了所有极值点的反演计算，如没有完成则对下一个极值点重复操作直至反演完所有的极值点，剔除重复的位置点，即得到航天器内部各磁源的位置信息。
2.如权利要求1所述的方法，其中，精度要求越高测试点越密集。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410446670.6A CN104182648B (zh) | 2014-09-03 | 2014-09-03 | 反演航天器内部多磁源分布的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410446670.6A CN104182648B (zh) | 2014-09-03 | 2014-09-03 | 反演航天器内部多磁源分布的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104182648A CN104182648A (zh) | 2014-12-03 |
| CN104182648B true CN104182648B (zh) | 2017-09-12 |
Family
ID=51963683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410446670.6A Expired - Fee Related CN104182648B (zh) | 2014-09-03 | 2014-09-03 | 反演航天器内部多磁源分布的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104182648B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105486943B (zh) * | 2015-11-19 | 2019-03-29 | 江西洪都航空工业集团有限责任公司 | 一种复杂环境下电子部件干扰磁场的测量方法 |
| JP7209176B2 (ja) * | 2018-10-26 | 2023-01-20 | スミダコーポレーション株式会社 | 磁場発生源検出装置および磁場発生源検出方法 |
| CN109613458B (zh) * | 2019-01-14 | 2020-12-01 | 北京卫星环境工程研究所 | 航天器磁异常定位方法 |
| CN109633494B (zh) * | 2019-01-14 | 2020-10-30 | 北京卫星环境工程研究所 | 航天器磁场分布信息成像方法 |
| CN111766549B (zh) * | 2020-07-07 | 2023-03-31 | 北京卫星环境工程研究所 | 一种可穿戴式磁场梯度探测仪及探测方法 |
| CN113124882B (zh) * | 2021-06-17 | 2021-09-28 | 天津大学 | 一种背景磁场未知情形下的多磁偶极子磁源反演定位方法 |
| CN113124881B (zh) * | 2021-06-17 | 2021-10-08 | 天津大学 | 一种基于磁信标的同步定位与构图系统的故障恢复方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6982520B1 (en) * | 2001-09-10 | 2006-01-03 | Aerojet-General Corporation | Hall effect thruster with anode having magnetic field barrier |
| CN102446242A (zh) * | 2011-12-09 | 2012-05-09 | 北京卫星环境工程研究所 | 航天器磁性的仿真方法 |
-
2014
- 2014-09-03 CN CN201410446670.6A patent/CN104182648B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6982520B1 (en) * | 2001-09-10 | 2006-01-03 | Aerojet-General Corporation | Hall effect thruster with anode having magnetic field barrier |
| CN102446242A (zh) * | 2011-12-09 | 2012-05-09 | 北京卫星环境工程研究所 | 航天器磁性的仿真方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于欧拉方法的多磁偶极子分辨技术;陈俊杰 等;《航天器环境工程》;20130815;第30卷(第4期);摘要,第401-403页第1-3小节,图1-4 * |
| 磁力梯度张量及在磁偶极子定位中的应用;吴招才 等;《2009年国家安全地球物理专题研讨会论文集》;20130709;第81-87页 * |
| 磁性物体磁偶极子模型适用性的试验研究;张朝阳 等;《应用基础与工程科学学报》;20101015;第18卷(第5期);第862-868页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104182648A (zh) | 2014-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104182648B (zh) | 反演航天器内部多磁源分布的方法 | |
| CN104535062B (zh) | 基于磁梯度张量和地磁矢量测量的运动式定位方法 | |
| JP5620475B2 (ja) | 磁気によって位置を求めるための方法及び装置 | |
| CN110333536B (zh) | 一种测距线性定位算法 | |
| CA2990468C (en) | A modified halbach nmr array sensor for analyzing core or fluid samples from a subsurface formation | |
| CN108197389A (zh) | 二维强磁性体磁场的快速、高精度数值模拟方法 | |
| CN109254327B (zh) | 三维强磁性体的勘探方法及勘探系统 | |
| CN103486984A (zh) | 一种风洞内型面同轴度的检测方法 | |
| CN105354850A (zh) | 基于电场性质的复杂曲面零件尺寸三维匹配检测方法 | |
| CN106940343A (zh) | 一种基于阵列电磁传感器的材料微损伤检测方法及系统 | |
| Tarnawski et al. | Improving the quality of magnetic signature reproduction by increasing flexibility of multi-dipole model structure and enriching measurement information | |
| Vuillermet et al. | Scalar potential formulation and inverse problem applied to thin magnetic sheets | |
| Guo et al. | Modeling on ground magnetic anomaly detection of underground ferromagnetic metal pipeline | |
| CN103955002A (zh) | 基于磁异常一阶导数的磁偶极子目标位置测量方法 | |
| CN107092001B (zh) | 基于面剖分与测点间磁场映射的船舶感应磁场测量方法 | |
| CN104391260A (zh) | 航天器及其产品磁矩自适应动态测试方法 | |
| Shur et al. | On the inverse problem of magnetostatics | |
| CN109633494B (zh) | 航天器磁场分布信息成像方法 | |
| CN107748834A (zh) | 一种计算起伏观测面磁场的快速、高精度数值模拟方法 | |
| CN115728829A (zh) | 一种基于磁梯度全张量的磁性目标定位方法 | |
| CN116736199A (zh) | 一种空间静态磁场分布测量系统及方法 | |
| Nowicki et al. | Determination of the location and magnetic moment of ferromagnetic objects based on the analysis of magnetovision measurements | |
| Nowicki et al. | Modelling of the magnetovision image with the finite element method | |
| Shubitidze et al. | Comparison of the physically complete model with a simple dipole model for UXO detection and discrimination | |
| Rosu et al. | Assessment of the ellipsoidal shell model for ship magnetic signature |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170912 Termination date: 20210903 |