徐超群，易忠，王斌，劉超波

?

基于磁場梯度張量的航天器內(nèi)部多磁場源探測技術(shù)

徐超群1, 2，易忠1, 2，王斌2，劉超波2

（1. 北京衛(wèi)星環(huán)境工程研究所可靠性與環(huán)境工程技術(shù)重點實驗室；2. 北京衛(wèi)星環(huán)境工程研究所：北京 100094）

基于磁場梯度張量的探測技術(shù)具有高分辨率、高精度等優(yōu)點，是航天器多磁源分辨的有效方法。文章建立了航天器磁場源模型，提出了航天器多磁場源模型擬合方法，利用磁場梯度張量的主不變量的極值確定磁場源個數(shù)和水平位置，然后聯(lián)合歐拉方程計算磁場源深度，完成航天器內(nèi)部多磁場源仿真與計算，分辨率最高可達0.012m，滿足工程需求。該研究開辟了航天器內(nèi)部多磁場源目標探測的新途徑。

航天器；多磁場源；擬合方法；磁梯度張量；主不變量

0 引言

航天器在研制過程中，為使其磁性控制達到要求，需對航天器整體磁矩進行測試和補償。目前航天器磁測試大多將航天器整體視作一個磁場源[1]，故測試結(jié)果是航天器內(nèi)部多個磁性組件的綜合信息，不能真實反映航天器內(nèi)部的磁場源狀況。隨著空間技術(shù)的快速發(fā)展，新的磁測試要求不僅需要對航天器整體磁矩做出評估，還需要分析航天器內(nèi)部主要磁性分布狀況，因此開展航天器內(nèi)部磁場源的分辨研究十分重要。

國外采用等效源方法對航天器內(nèi)多磁場源進行擬合[2]，但此方法不能真實反映星內(nèi)磁場源分布情況。目前我國幾乎沒有針對小尺度范圍（航天器尺度）內(nèi)的多磁場源探測技術(shù)的研究應(yīng)用，有人提出用歐拉方法[3]求解航天器內(nèi)部磁場源，但該方法漏解、多解情況嚴重，用于衛(wèi)星的磁測主要集中在消除磁干擾[4-5]。類似的磁探測技術(shù)主要用于海洋、陸地磁測[6-7]等領(lǐng)域。

本文首先建立航天器多磁場源模型，然后給出模型的逼近擬合方法，利用磁場梯度張量的主不變量的極大值對多磁場源目標進行預(yù)判，通過聯(lián)合歐拉反演方法，完成航天器內(nèi)部多磁場源測量與計算。

1 航天器磁場源建模

對于單個理想磁場源（即只有磁矩和位置信息的質(zhì)點），可以對磁場公式進行勒讓德函數(shù)多項式展開，通過對展開式各階項的計算來分析磁場細部特征。在多數(shù)情況下，只保留展開式中的磁偶極子項，而忽略后面的高階項，因為它們隨距離衰減得很快。

一個磁偶極子項可表示為

式中：為磁場源的磁矩；為由磁場源到探測點的位移矢量；0為真空磁導(dǎo)率，H/m。

在實際工況中，當探測點至磁場源的距離為源自身尺寸的2.5倍或以上時，可采用磁偶極子近似法[8]進行磁場源分析。航天器內(nèi)部磁場源情況復(fù)雜，磁性部件的磁矩、位置和形狀各不相同，而磁偶極子近似法是處理多目標磁場源問題的重要手段。在初步的航天器磁性建模中，可忽略磁性實體的形狀大小，只考慮其磁性位置和磁矩參數(shù)，但要注意兩點：一是忽略磁矩相對較小的磁場源；二是把十分靠近的磁場源當做一個等效源來處理。最終，航天器可被等效為一個內(nèi)部包含多個理想磁場源的立方體，如圖1所示，紅色點表示位置，箭頭表示磁矩。

圖1 航天器磁偶極子模型

使用磁偶極子近似法得到的航天器總磁場為

Σ。 (2)

把航天器的一個表面作為探測平面，則多磁偶極子在這個探測平面的磁通密度分布如圖2所示。

圖2 模型上方磁通密度分布

2 航天器磁性擬合技術(shù)

建立航天器多磁場源模型后，常常利用逼近擬合的方法來推斷各磁場源的位置和磁矩信息。其方法如下：f

1）測量計算

①設(shè)置航天器的參考坐標系和測量范圍；

②設(shè)置航天器磁場的測量平面（曲面）和測量方式；

③在測量平面布置磁場探測器；

④計算反演航天器內(nèi)部多個等效磁場源的位置、磁矩信息。

2）擬合優(yōu)化

①設(shè)定磁場源數(shù)目、位置和磁矩參數(shù)的初始值，建立模擬模型；

②計算模擬模型在測量平面的磁場數(shù)據(jù)信息；

③對比模擬數(shù)據(jù)與真實測量數(shù)據(jù)，求出兩者的差值；

④根據(jù)一定的算法程序修改模擬模型中的各項參數(shù)；

⑤重復(fù)步驟②、③，直到模擬值與測量值的差值滿足一定的允差條件。

在逼近擬合方法中，K Mehelm等[2]提出了一種簡單、有效的算法。近些年來，很多研究把人工智能算法（如遺傳算法、神經(jīng)網(wǎng)絡(luò)算法等）引入多磁場源定位問題，但是，這些方法在應(yīng)用過程中都遇到一個難題，即初始值的選?。汉芸赡軐?dǎo)致冗繁的計算量或是離真實情況偏差大。因此，本文引入磁場梯度張量來改善擬合初始值的選取。

3 磁場梯度張量

磁場梯度張量可表示為

根據(jù)電磁場麥克斯韋方程組，對于磁場中的無源區(qū)，有下列關(guān)系成立：

。(4)

這樣，是1個對稱張量，實際只有5個獨立分量，這就為實際測量提供了便利。在實際工況中一般通過測量磁感應(yīng)強度，再利用差分計算[9]來獲得。基于上述測量原理的自制磁場梯度張量測量結(jié)構(gòu)[10]如圖3所示。

根據(jù)實對稱張量的主不變量，即

可在探測平面內(nèi)繪制出如圖4所示的“斑點”圖。

圖4 多源磁場梯度張量主不變量I2平面分布

圖中每個“斑點”代表對應(yīng)位置下方存在一個磁場源[11]。這種通過探測主不變量極大值來尋找目標磁場源水平位置的方法，經(jīng)常應(yīng)用于地質(zhì)領(lǐng)域的大面積探測[12]。

4 數(shù)據(jù)采集

計算所需的磁場矢量、磁場梯度張量數(shù)據(jù)一般采用網(wǎng)格法探測，即：把測量平面分為個網(wǎng)格，采集每個結(jié)點的磁感應(yīng)強度。這種方法可以并行采集結(jié)點數(shù)據(jù)，還能靈活調(diào)整網(wǎng)格邊距（基線），關(guān)鍵是可以利用差分法消除背景噪聲干擾，得到磁場梯度張量。理想的采集系統(tǒng)如圖5所示。

圖5 試驗數(shù)據(jù)采集系統(tǒng)示意

5 仿真計算

基本的歐拉齊次方程為

式中：(0,0,0)是磁場源所在位置；探測點(,,)處測得總場為，背景場為，它可被差分法消除；齊次方程系數(shù)稱作結(jié)構(gòu)因子，單個磁偶極子點源的磁場隨距離成三次方衰減，其結(jié)構(gòu)因子為3。根據(jù)磁偶極子的表達式，都可以用B、B、B替換而滿足上述方程。因此，在每一個測量點，都可以得到如下的一個方程組：

。 (7)

對于探測區(qū)域只有單個磁場源的情形，將測量的磁場矢量、磁場梯度張量數(shù)據(jù)代入式(7)即可準確得到磁場源的位置(0,0,0)；而當探測區(qū)域中同時存在多個磁場源時，所有測點的計算值在探測區(qū)域（航天器范圍）內(nèi)形成一定的分布，而分布密度最大的地方，即最大似然值可以作為當下目標的真實位置。將此位置數(shù)據(jù)代入式(1)，便可求得磁場源磁矩參數(shù)。

這里利用磁場梯度張量的主不變量（2）極大值來確定目標磁場源個數(shù)和水平位置，然后把水平位置代入式(7)，即可得到磁偶極子的埋深。

設(shè)置航天器模型尺寸為3m×3m×3m，預(yù)設(shè)內(nèi)部磁場源4個，反演結(jié)果如圖6所示，其中，紅點表示磁場源的位置，箭頭表示磁矩；藍點表示磁場梯度張量反演計算位置。該模型磁場源的原始值（數(shù)字）與反演值（字母）如表1所示。

圖6 埋深3.0 m內(nèi)的多磁場源反演

表1 數(shù)值模擬值與預(yù)設(shè)值比較

通過對比，可以看出利用磁場梯度張量的主不變量探測得到的磁源水平位置具有較高的準確度，空間距離上的最小分辨率可以達到0.012m。

6 結(jié)束語

多磁場源分辨一直是磁測領(lǐng)域的難題。細致的分辨主要依靠逼近擬合的“猜測”式算法。而利用磁場梯度張量算法，可以快速判定磁場源數(shù)目、位置、磁矩等信息，可以為逼近擬合方法的初始值選擇提供參考，加快計算的收斂速度，同時提高計算結(jié)果的可靠性。

本文利用磁場梯度張量的主不變量極大值來確定目標磁場源個數(shù)和水平位置。為了提高準確性，可測量2～3平面并進行統(tǒng)計，然后利用歐拉方程反演磁偶極子埋深，可同時準確測量和反演航天器內(nèi)部多個磁場源的位置信息。對于小范圍多磁場源相互干擾的情況，磁場梯度張量反演計算可以獲得較為準確、可靠的結(jié)果。

這種多目標分辨的方法可對航天器內(nèi)部的磁場分布進行判斷，精確探測和反演多個磁場源位置，提高了航天器磁測試和磁補償?shù)尼槍π院蜏蚀_度，是未來航天器磁測試的理想方法。

[1] 易忠, 孟立飛. “探測”雙星磁性仿真分析計算[J]. 裝備環(huán)境工程, 2006, 3(3): 37-42

YI Z, MENG L F. Calculation and analysis of magnetic simulation for the Double Star[J]. Equipment Environmental Engineering, 2006, 3(3): 37-42

[2] MEHLEM K, WIEGAND A. Magnetostatic cleanliness of spacecraft[C]//Electromagnetic Compatibility (APEMC), 2010 Asia-Pacific Symposium. Beijing, 2010: 936-944

[3] 陳俊杰, 易忠, 孟立飛, 等. 基于歐拉方法的多磁偶極子分辨技術(shù)[J]. 航天器環(huán)境工程, 2013, 30(4): 401-406

CHEN J J, YI Z, MENG L F, et al. Multi-dipole discrimination technology based on Euler inverse method[J]. Spacecraft Environment Engineering, 2013, 30(4): 401-406

[4] YIN F. Mathematic approaches for the calibration of the CHAMP satellite magnetic field measurements[D]. Potsdam: The University of Potsdam, 2010

[5] 周斌, 王勁東. 衛(wèi)星磁部件分布對梯度法消除剩磁的影響分析[J]. 中國空間科學(xué)技術(shù), 2013, 33(5): 29-34

ZHOU B, WANG J D. Influence of magnetic component distribution of satellite on eliminating remainant magnetic field by gradient method[J]. Chinese Space Science and Technology, 2013, 33(5): 29-34

[6] PEI Y H, YEO H G, KANG X Y, et al. Magnetic gradiometer on an AUV for buried object detection[C]//Oceans 2010. Seattle, 2010: 1-8. DI: 10.1109/OCEANS.2010.5664272

[7] 柳建新, 孫歡樂, 陳波, 等. 重磁方法在國內(nèi)外金屬礦中的研究進展[J]. 地球物理學(xué)進展, 2016, 31(2): 0713-0722

LIU J X, SUN H L, CHEN B, et al. Review of the gravity and magnetic methods in the exploration of metal deposits[J]. Progress in Geophysics, 2016, 31(2): 0713-0722

[8] 張朝陽, 肖昌漢, 高俊吉, 等. 磁性物體磁偶極子模型適用性的試驗研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報, 2010, 18(5): 862-867

ZHANG Z Y, XIAO C H, GAO J J, et al. Experiment research of magnetic dipole model applicability for a magnetic object[J]. Journal of Basic Science and Engineering, 2010, 18(5): 862-867

[9] 劉超波, 王斌, 陳金剛, 等. 磁場全張量測量計算方法與誤差分析[J]. 航天器環(huán)境工程, 2015, 32(1): 63-66

LIU C B, WANG B, CHEN J G, et al.Measurement method and error analysis of magnetic field full tensor[J]. Spacecraft Environment Engineering, 2015, 32(1): 63-66

[10] 徐超群, 易忠, 劉超波, 等. 十字形磁梯度張量結(jié)構(gòu)在多磁偶極子測量中的誤差與修正方法[J]. 航天器環(huán)境工程, 2017, 34(3): 312-316

XU C Q, YI Z, LIU C B, et al. The error and correction methods of the cross-shaped magnetic gradient tensor in the multi-magnetic dipole measurement[J]. Spacecraft Environment Engineering, 2017, 34(3): 312-316

[11] PEDERSEN L B, RASMUSSEN T M. The gradient tensor of potential field anomalies: some implications on data collection and data processing of maps[J]. Geophysics, 1990, 55(12): 1558-1566

[12] SCHMIDT P, CLARK D, LESLIE K, et al. GETMAG: A SQUID magnetic tensor gradiometer for mineral and oil exploration[J]. Exploration Geophysics, 2004, 35(4): 297-305

（編輯：許京媛）

Detection technology of multi-magnetic source in spacecraft based on magnetic field gradient tensor

XU Chaoqun1,2, YI Zhong1,2, WANG Bin2, LIU Chaobo2

(1. Science and Technology on Reliability and Environmental Engineering Laboratory, Beijing Institute of Spacecraft Environment Engineering; 2. Beijing Institute of Spacecraft Environment Engineering: Beijing 100094, China)

The multi-magnetic source resolution in spacecraft is a difficult problem in the field of magnetic survey. The technology based on the magnetic gradient tensor can be used for high resolution and high precision detection. In this paper, a spacecraft magnetic source model is established, and a multi-magnetic source model fitting method for spacecraft is proposed. The principal invariants of the magnetic field gradient tensor are used to determine the number and the horizontal location of the source, and the Euler equations are combined to complete the multi-magnetic source depth calculation, with best resolution up to 0.012m, meeting the engineering requirements. The study opens up a new way of the multi-magnetic source target detection in spacecraft.

spacecraft; multi-magnetic source; fitting method; magnetic gradient tensor; principal invariants

V416.5

A

1673-1379(2017)04-0398-05

10.3969/j.issn.1673-1379.2017.04.010

徐超群（1985—），男，博士研究生，研究方向為航天器磁環(huán)境工程；E-mail: xucq111@163.com。

易忠（1968—），男，博士生導(dǎo)師，研究員，主要從事航天器磁環(huán)境效應(yīng)與測試技術(shù)研究；E-mail: yizhong6808@sina.com。

2017-02-22；

2017-07-13

國家自然科學(xué)基金青年科學(xué)基金項目“基于磁場梯度張量的多磁源目標反演方法研究”（編號：51207011）

XU C Q,YI Z, WANG B, et al. Detection technology of multi-magnetic source in spacecraft based on magnetic field gradient tensor[J].Spacecraft Environment Engineering, 2017, 34(4): 398-402