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

?

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

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

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

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

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

0 引言

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

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

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

1 航天器磁場源建模

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

一個磁偶極子項(xiàng)可表示為

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

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

圖1 航天器磁偶極子模型

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

Σ。 (2)

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

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

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

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

1）測量計(jì)算

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

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

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

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

2）擬合優(yōu)化

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

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

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

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

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

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

3 磁場梯度張量

磁場梯度張量可表示為

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

。(4)

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

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

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

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

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

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

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

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

5 仿真計(jì)算

基本的歐拉齊次方程為

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

。 (7)

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

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

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

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

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

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

6 結(jié)束語

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

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

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

[1] 易忠, 孟立飛. “探測”雙星磁性仿真分析計(jì)算[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ìn)展[J]. 地球物理學(xué)進(jìn)展, 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] 張朝陽, 肖昌漢, 高俊吉, 等. 磁性物體磁偶極子模型適用性的試驗(yàn)研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 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ì)算方法與誤差分析[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—），男，博士研究生，研究方向?yàn)楹教炱鞔怒h(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é)基金項(xiàng)目“基于磁場梯度張量的多磁源目標(biāo)反演方法研究”（編號：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