不同衛(wèi)星軌道高度對地殼磁場反演的影響研究

種洋， 常宜峰， 柴洪洲， 郭云飛， 敬子健， 劉毅

1 信息工程大學(xué)地理空間信息學(xué)院， 鄭州 450001 2 軍事科學(xué)院戰(zhàn)爭研究院， 北京 100091 3 自然資源部海洋環(huán)境探測技術(shù)與應(yīng)用重點實驗室， 廣州 510300 4 95948部隊， 蘭州 732750 5 93920部隊， 西安 710061 6 92292部隊， 青島 266000

0 引言

地磁場同重力場一樣，也是地球的基本物理場之一，在地質(zhì)與板塊研究、礦產(chǎn)資源勘探、地磁導(dǎo)航等領(lǐng)域發(fā)揮著重要的應(yīng)用價值.隨著國際上對衛(wèi)星磁力測量研究的逐步深入以及更多高精度衛(wèi)星磁測數(shù)據(jù)的積累，如何有效利用和挖掘衛(wèi)星磁測數(shù)據(jù)中的地磁場信息，成為亟待進一步研究的重要課題(徐文耀，2009；Hulot et al.，2010; 常宜峰，2015).相比于海洋磁測、航空磁測和地面磁測，衛(wèi)星磁測是唯一可以快速獲得全球均勻分布磁測數(shù)據(jù)的有效方法.利用衛(wèi)星磁測數(shù)據(jù)不僅可以構(gòu)建全球主磁場模型，同時也可研究全球范圍內(nèi)地磁異常場的分布情況，以及空間電離層、磁層等外源磁場的結(jié)構(gòu)及特點(Langlais et al.，2010；Alken et al.，2015；Sabaka et al.，2015；Vigneron et al.，2015).

國內(nèi)外學(xué)者都對衛(wèi)星磁測數(shù)據(jù)處理及其相關(guān)應(yīng)用方面進行了研究.安振昌等(1998)使用冠諧分析法對歐洲及其相鄰地區(qū)MAGSAT衛(wèi)星黎明資料處理，獲得1°×1°衛(wèi)星磁異常網(wǎng)格值.綜合IGRF模型數(shù)據(jù)和MAGSAT衛(wèi)星實測地磁數(shù)據(jù)，左文輯和宋福香等(2000)利用擴展卡爾曼濾波算法來確定衛(wèi)星的位置和速度，實現(xiàn)微小近地衛(wèi)星自主導(dǎo)航.Hulot等(2002)利用MAGSAT和Oersted衛(wèi)星磁測數(shù)據(jù)進行了地球磁場電流衰退研究.Macmillan和Thomson等(2003)經(jīng)過對MAGSAT、Oersted以及地磁臺站觀測數(shù)據(jù)處理，研究了地殼磁場變化和時間的相關(guān)性，分析得出磁層磁場和電離層可能是地殼磁場產(chǎn)生變化的主要原因.Maus等(2009)聯(lián)合利用臺站、CHAMP衛(wèi)星、船載和航空等磁測數(shù)據(jù)，構(gòu)建了全球地磁異常網(wǎng)格圖(Earth Magnetic Anomaly Grid 2，EMAG2)模型，精度約2弧分.針對球諧函數(shù)僅能對衛(wèi)星磁測數(shù)據(jù)低頻率球面場進行擬合的問題，王慧琳(2010)提出多尺度混合球面小波方法來融合衛(wèi)星、航空、地面以及海洋等多源磁測數(shù)據(jù)，完成低中高三種頻率的空間局部地磁數(shù)據(jù)擬合.康國發(fā)等(2010)利用CHAMP衛(wèi)星磁測資料構(gòu)建了高階地磁場模型POMME-4.2S，并針對中國及其鄰近地區(qū)400 km高度的衛(wèi)星磁異常和垂直梯度進行了計算分析，得出巖石圈構(gòu)造與磁異常存在相關(guān)性.白春華等(2011)系統(tǒng)總結(jié)了CHAMP衛(wèi)星磁測數(shù)據(jù)處理方法對地磁場建模的影響，重點研究了如何剔除來自電離層和磁層等外源場的噪聲干擾、海洋潮汐引起的地磁信號和模型估計需要顧及高度偏差等方面.區(qū)家明等(2012)結(jié)合CHAMP衛(wèi)星與地磁臺站的磁測數(shù)據(jù)構(gòu)建了中國及鄰近地區(qū)的巖石圈磁場模型，得到的磁場模型空間分辨率達到150 km.羅開奇和張有明(2012)研究了衛(wèi)星磁測數(shù)據(jù)的通化處理技術(shù)，構(gòu)建了川滇地區(qū)的泰勒多項式模型，為衛(wèi)星磁測數(shù)據(jù)的廣泛應(yīng)用提供了新思路.馮麗麗(2015)綜合利用衛(wèi)星和地面的磁測數(shù)據(jù)資料，構(gòu)建了中國及鄰近地區(qū)的諧和樣條地磁模型.杜勁松和陳超(2015) 指出衛(wèi)星磁測數(shù)據(jù)在全球巖石圈磁場建模方面具有獨特優(yōu)勢，詳細研究等效源方法在衛(wèi)星磁測數(shù)據(jù)處理方面的進展，為衛(wèi)星磁測數(shù)據(jù)的標(biāo)定、數(shù)據(jù)篩選及主磁場改正等技術(shù)的改進提供重要參考.Kotsiaros等(2015)利用近14個月的Swarm衛(wèi)星磁測數(shù)據(jù)恢復(fù)的巖石圈磁場只達到70階，但由反演結(jié)果可以得出在內(nèi)源場的建模中，差分地磁梯度數(shù)據(jù)發(fā)揮了重要作用.姜乙等(2015)結(jié)合CHAMP的衛(wèi)星磁測數(shù)據(jù)和地面臺站數(shù)據(jù)，通過球諧模型建立了全球的大尺度地磁場模式，并計劃開展模型磁場數(shù)據(jù)的實時3D顯示等工作.在假設(shè)CHAMP衛(wèi)星498個網(wǎng)格測點位于同一高度(307.69 km)的基礎(chǔ)上，基于該衛(wèi)星的矢量數(shù)據(jù)，F(xiàn)eng等(2016)結(jié)合最新的IGRF12模型構(gòu)建了2010.0年中國地區(qū)的衛(wèi)星磁異常球冠諧(SCH2010)模型，可反映地磁場更多細節(jié)信息.利用CHAMP衛(wèi)星地磁總強度F的觀測數(shù)據(jù)，王粲等(2017)通過相關(guān)參數(shù)選取地磁活動相對平靜期間的衛(wèi)星磁測數(shù)據(jù)，剔除主磁場和外源場部分信息得到中國及鄰區(qū)的地磁異常分布.馮彥等(2018)等基于三個高度的數(shù)據(jù)，包括地面實測數(shù)據(jù)以及子午工程測點數(shù)據(jù)，CHAMP衛(wèi)星高度實測數(shù)據(jù)，并結(jié)合高空的IGRF12(International Geomagnetic Reference Field 12)數(shù)據(jù)點，首次建立了中國地區(qū)地磁場3D Spline模型，結(jié)果表明該模型在不同高度的模擬結(jié)果取決于相應(yīng)高度的地磁實測值數(shù)量和精度.

從前期的研究可以看出，很多學(xué)者在衛(wèi)星磁測數(shù)據(jù)的處理、衛(wèi)星磁測數(shù)據(jù)高精度地磁場模型的構(gòu)建等方面做了大量工作.為了反演出更加精細的地殼磁場模型，就需要獲得具有一定強度、包含較寬頻率范圍的磁測衛(wèi)星觀測信息，這和衛(wèi)星軌道的高度設(shè)計緊密相關(guān).然而磁測衛(wèi)星在不同的軌道高度下得到的觀測信息是不同的.通常，磁測衛(wèi)星的軌道高度選擇范圍為距地球表面約1000 km以內(nèi)的近地圓形極軌模式.當(dāng)衛(wèi)星軌道較低時，磁測數(shù)據(jù)可以敏感到更加精細的高頻地磁場信息成分，但由于空氣阻力等因素的存在導(dǎo)致磁測衛(wèi)星壽命縮短，不能進行對地長期觀測.而當(dāng)衛(wèi)星軌道較高時，一方面，磁測數(shù)據(jù)會對地殼磁場的敏感性大大降低；另一方面，由于太陽風(fēng)磁場的影響，磁測數(shù)據(jù)質(zhì)量會受到嚴重影響.

因此，建模目的不同，對于衛(wèi)星軌道的高度選擇也會有所不同，就需要進一步研究分析衛(wèi)星軌道高度對地殼場反演的影響.研究如何設(shè)計更加合理的衛(wèi)星軌道高度顯得十分必要，本文通過仿真和實測實驗來研究衛(wèi)星軌道高度對地殼磁場反演的影響，為我國磁測衛(wèi)星系統(tǒng)總體設(shè)計提供一定的參考.

1 基本原理

地磁場的磁位勢滿足高斯理論，磁位勢U的表達式為

(1)

對磁位勢U在X、Y、Z三分量方向上分別進行求導(dǎo)，就可以得到磁位勢U在X、Y、Z方向上三分量的表達式為

(2)

(3)

根據(jù)式(3)計算，即可得全球地磁場總強度F的模型值.

1.1 地殼磁場的反演

1.2 ML地磁能譜

(4)

其中，r是研究球?qū)拥母叨龋琑n即為與球諧系數(shù)階數(shù)n相關(guān)的ML能譜.

由式(4)可得，位系數(shù)的平方和會隨階數(shù)的增大而減小.同時，系數(shù)因子(n+1)(R/r)2n+4呈現(xiàn)出指數(shù)變化的趨勢，通常有r≥R.因此ML地磁能譜是隨著高度增大而逐漸減小的.

1.3 線性擬合趨勢項

綜上，衛(wèi)星軌道高度與地殼磁場反演階數(shù)間的關(guān)系分析步驟總結(jié)如下：

①由衛(wèi)星觀測得到每點X、Y、Z三分量的磁測數(shù)據(jù)，根據(jù)式(2)列立方程組，將高斯系數(shù)作為未知參數(shù)改寫成標(biāo)準(zhǔn)的觀測方程.

②通過最小二乘原理，對觀測方程進行求解，可得高斯系數(shù).

③基于不同階數(shù)的高斯系數(shù)，根據(jù)式(4)求得不同階數(shù)對應(yīng)的ML能譜值.

④根據(jù)所求n與Rn之間的變化關(guān)系，對散點進行線性擬合，得到ML能譜的變化趨勢.

2 軌道高度仿真數(shù)據(jù)分析

首先介紹實驗所用衛(wèi)星磁測數(shù)據(jù)的仿真方法及特點，其次分析了不同軌道高度所測得F的變化特征，最后利用本文所提出的方法研究了不同軌道高度的衛(wèi)星磁測數(shù)據(jù)反演得到的地殼磁場最大階數(shù)問題.

2.1 仿真數(shù)據(jù)的生成

基于衛(wèi)星地磁場模型POMME-11，選取球諧系數(shù)中16～133階的位系數(shù)，根據(jù)式(2)仿真計算出全球地殼磁場X、Y、Z三分量.參考POGO、MAGSAT、Oersted、SAC-C、CHAMP和Swarm等磁測衛(wèi)星的基本軌道參數(shù)，分別選取250 km、350 km、450 km、550 km、650 km、750 km、1050 km等高度的1°×1°均勻格網(wǎng)點來模擬不同磁測衛(wèi)星的軌道測點，并結(jié)合由式(2)仿真得到的全球地殼磁場X、Y、Z三分量，即得衛(wèi)星磁測的仿真數(shù)據(jù).根據(jù)衛(wèi)星磁測數(shù)據(jù)預(yù)處理的數(shù)據(jù)篩選方法(杜勁松，2015)，對磁測數(shù)據(jù)進行時間格式和坐標(biāo)系統(tǒng)的統(tǒng)一，顧及緯度條件選取緯度范圍±86°以內(nèi)的衛(wèi)星磁測數(shù)據(jù)，得到不同高度下的衛(wèi)星磁測仿真數(shù)據(jù).

2.2 軌道高度與地磁場總強度F的關(guān)系

選取0°經(jīng)線上不同緯度處的特征點，并以20°為步長計算緯度范圍為-74°～86°的9個特征點的F值來對比不同高度和緯度處F值的變化特征，結(jié)果如圖1.

圖1 不同高度下F的變化特征Fig.1 The variation characteristics of F at different altitudes

不同的顏色曲線表示不同緯度下特征點的F值隨高度的變化情況.由圖1可得，在相同緯度情況下，隨著軌道高度的增加，F(xiàn)值均隨著軌道高度增加呈現(xiàn)出減小的趨勢；在不同緯度情況下，與中低緯度地區(qū)的F值相比，高緯度地區(qū)的F值波動范圍更大；在相同軌道高度情況下，兩極方向的高緯度地區(qū)F值要大于對應(yīng)的低緯度地區(qū)F值.

2.3 不同高度仿真磁測數(shù)據(jù)反演結(jié)果比較

選取軌道高度分別為250 km、350 km、450 km和750 km的衛(wèi)星磁測數(shù)據(jù)，分別反演出不同軌道高度下相應(yīng)地殼磁場的高斯系數(shù).通過ML能譜分析法計算得到不同軌道高度下對應(yīng)的ML能譜變化曲線，并利用線性擬合法給出不同曲線的變化趨勢項特征，結(jié)果如圖2.

圖2 不同軌道高度下的地磁ML能譜分布Fig.2 The distribution of geomagnetic ML energy spectrum at different orbital altitudes

圖2下半部分的4條實線分別表示在不同軌道高度下衛(wèi)星磁測反演得到的高斯系數(shù)ML能譜分布，4條實線相對應(yīng)的虛線部分能夠有效刻畫出往高階方向反演得到的變化趨勢.由于衛(wèi)星和地面之間高度差的存在，導(dǎo)致觀測到的地殼磁場信息有所衰減，進而其ML能譜值要比實際的地殼磁場的ML能譜值小，并隨著軌道高度的增加呈現(xiàn)出下降的趨勢.只有當(dāng)ML能譜標(biāo)準(zhǔn)大于0.1 nT2時，才能探測到各軌道高度處反演能力的趨勢變化特征(馮彥等，2014).將ML能譜標(biāo)準(zhǔn)的下限取為0.1 nT2時，軌道高度250 km、350 km、450 km和750 km對應(yīng)的最高反演階數(shù)分別為163階、108階、86階和42階.結(jié)果反映了不同軌道高度下，通過衛(wèi)星磁測數(shù)據(jù)反演地殼磁場最大階數(shù)的能力.

對不同軌道高度與地殼磁場反演最高階數(shù)間的數(shù)值關(guān)系進行曲線擬合分析，得出衛(wèi)星軌道高度與地殼磁場最高反演階數(shù)間的變化關(guān)系曲線，結(jié)果如圖3.

圖3 軌道高度與反演階數(shù)間的變化關(guān)系曲線Fig.3 The relationship between orbital altitude and inversion order

很明顯，由圖3可得，隨著衛(wèi)星的軌道高度降低，地殼磁場的最高反演階數(shù)會逐漸增大.由軌道高度與反演階數(shù)間的變化關(guān)系曲線，可以發(fā)現(xiàn)反演階數(shù)的增長幅度變化存在一個臨界值.根據(jù)擬合曲線的斜率變化情況，通過比較分析，當(dāng)軌道高度大于450 km時，擬合曲線的斜率較大，即隨著軌道高度的降低，反演階數(shù)的增長幅度相對平緩；而當(dāng)在軌道高度小于450 km時，擬合曲線的斜率明顯縮小，即隨著軌道高度的降低，反演階數(shù)的增長幅度逐漸變大.因此，可以將軌道高度450 km作為衛(wèi)星磁力測量系統(tǒng)軌道高度設(shè)計的參考高度，同時這也與當(dāng)前主流的CHAMP和Swarm磁測衛(wèi)星的軌道高度保持一致.

3 軌道高度實測數(shù)據(jù)分析

MF系列地磁場模型是使用了不同時間段內(nèi)的CHAMP衛(wèi)星磁測數(shù)據(jù)進行反演得到的地殼磁場，隨著后期衛(wèi)星軌道高度的逐漸降低，MF系列模型所對應(yīng)的衛(wèi)星平均軌道高度也是逐漸減小的.隨著衛(wèi)星軌道高度的降低，使得CHAMP衛(wèi)星對高頻的地磁場信息更加敏感，相應(yīng)地衛(wèi)星磁測數(shù)據(jù)中也就包含了更多的短波長信息，進而為反演高階地殼磁場模型提供了可能.將MF系列模型的階數(shù)、數(shù)據(jù)來源以及基本特點進行了總結(jié)，如表1所示.

表1 MF系列模型簡介Table 1 Introduction to series of MF models

由表1可得，MF系列模型一般是3～5年的衛(wèi)星磁測數(shù)據(jù)進行反演獲得，模型階數(shù)選擇的是16階以上的高階部分(地殼場部分)，反演前需要對矢量或標(biāo)量數(shù)據(jù)進行合理的篩選和磁測數(shù)據(jù)預(yù)處理.從MF4模型的反演最大階數(shù)90階增加到MF7模型的反演最大階數(shù)133，也意味著衛(wèi)星軌道高度的逐漸降低和衛(wèi)星磁測數(shù)據(jù)質(zhì)量的提高.

如圖4所示，3個虛線框內(nèi)標(biāo)記出了MF系列模型所使用的CHAMP衛(wèi)星磁測數(shù)據(jù)區(qū)間，同時用實線標(biāo)出了相應(yīng)模型的數(shù)據(jù)平均軌道高度以及發(fā)布時間.由左到右的3個虛線框內(nèi)分別代表MF4、MF6和MF7地殼磁場模型所使用的磁測數(shù)據(jù)，對應(yīng)的平均軌道高度分別約為430 km、370 km和320 km，對應(yīng)的發(fā)布時間分別為2005年6月、2007年12月和2010年8月.

圖4 CHAMP衛(wèi)星高度隨時間變化Fig.4 Altitude varies with time from CHAMP satellite

通過使用不同軌道高度下的CHAMP衛(wèi)星磁測數(shù)據(jù)反演得到MF系列模型，求得相應(yīng)模型的ML能譜分布，與POMME11地殼磁場模型的ML能譜分布對比結(jié)果如圖5.

從圖5可以看出，POMME11和MF7模型的ML能譜分布具有較強的一致性，兩者能夠在整個地殼磁場的階數(shù)范圍內(nèi)幾乎完全重合，ML能譜分布范圍基本在10～100 nT2.當(dāng)階數(shù)小于60階時， MF系列模型同POMME11模型具有較高的重合度，ML能譜分布特征相近；當(dāng)階數(shù)范圍為60～120階時，MF系列模型的ML能譜分布保持平穩(wěn)變化，僅MF4模型在80階之后出現(xiàn)了較為明顯的下降趨勢；當(dāng)階數(shù)大于120階時，MF6模型的ML能譜值急劇下降.分析出現(xiàn)上述變化趨勢的原因：MF系列模型的構(gòu)建是聯(lián)合使用了衛(wèi)星、航空和地面臺站等多源磁測數(shù)據(jù)，使其在地殼磁場的反演過程中(尤其是低階部分)仍保持了較高的ML能譜，隨著階數(shù)的增加，ML能譜曲線會出現(xiàn)急劇下降的趨勢，這也進一步反映了衛(wèi)星軌道高度和地殼磁場反演階數(shù)之間具有一定的相關(guān)性.

圖5 MF系列模型與POMME11模型對比Fig.5 Comparison between MF series model and POMME11 model

圖6給出了POMME11和MF7模型分別反演得到的距離地表1 km高度的全球地殼磁場分布情況來直觀描述衛(wèi)星軌道高度和地殼磁場反演階數(shù)的相關(guān)性.其中，左圖代表POMME11模型計算得到的全球地殼磁場分布，右圖代表MF7模型的分布結(jié)果.

從圖6可以看出，與MF7模型相比，POMME11模型能夠刻畫出更為精細的全球地殼磁場分布情況.MF7模型的局部精細地磁異常雖然已經(jīng)被忽略掉，但如非洲中部的班吉異常、歐洲中東部庫爾斯克異常、北美洲中部的肯德基異常以及南極威爾克斯異常等全球較為明顯的地磁異常分布基本得到呈現(xiàn).圖6的結(jié)果更加直觀地印證了衛(wèi)星軌道高度和地殼磁場反演結(jié)果具有相關(guān)性.通過衛(wèi)星軌道高度的降低來構(gòu)建更為精細的地磁異常分布，同時使得地殼磁場反演階數(shù)有所提高，為構(gòu)建更符合實際的地殼磁場模型提供了重要的指導(dǎo)和參考.

圖6 POMME11模型和MF7模型的全球地殼磁場分布情況對比Fig.6 Distribution comparison in the global crust magnetic field between POMME11 model and MF7 model

4 結(jié)論

通過仿真和實測數(shù)據(jù)實驗分析得出，衛(wèi)星軌道高度和磁測數(shù)據(jù)中包含的信息量有關(guān).隨著衛(wèi)星軌道高度逐漸降低，對應(yīng)磁測數(shù)據(jù)中將包含更多短波長的地殼磁場特征信息.基于已有實驗結(jié)果，得出以下結(jié)論：

(1)衛(wèi)星軌道高度和地殼磁場反演的最高階數(shù)之間具有較強的相關(guān)性.通常高度在300～500 km的衛(wèi)星軌道，對應(yīng)最高反演階數(shù)大致范圍為130～80階，在反演階數(shù)固定的條件下，衛(wèi)星軌道高度也會對應(yīng)一個最小高度值.顧及衛(wèi)星軌道高度和反演階數(shù)的相關(guān)關(guān)系，能夠為衛(wèi)星磁測數(shù)據(jù)處理以及我國自主研制的磁測衛(wèi)星提供一定參考.

(2)隨著衛(wèi)星的軌道高度降低，地殼磁場的最高反演階數(shù)會逐漸增大，進而反映更多的地磁場信息，構(gòu)建更為精細的地磁異常分布，軌道高度450 km可以作為磁測衛(wèi)星系統(tǒng)軌道高度設(shè)計的參考高度.通過比較分析，當(dāng)軌道高度大于450 km時，擬合曲線的斜率較大，即隨著軌道高度的降低，反演階數(shù)增長幅度相對平緩；而當(dāng)在軌道高度小于450 km時，擬合曲線的斜率明顯縮小，即隨著軌道高度的降低，反演階數(shù)的增長幅度逐漸變大.

(3)衛(wèi)星磁測是一個龐大的系統(tǒng)工程，磁測數(shù)據(jù)的處理應(yīng)用主要由地磁場的復(fù)雜性、磁測載荷儀器的精度以及磁測衛(wèi)星的組網(wǎng)方案設(shè)計等共同決定.隨著磁力測量儀器精度的不斷提高和多星協(xié)同探測優(yōu)化設(shè)計的開展，地殼磁場反演階數(shù)也會逐步提升.

致謝感謝GFZ提供的部分實驗數(shù)據(jù)和地磁場模型以及國家自然科學(xué)基金、地理信息工程國家重點實驗室開放基金對本文提供的資助.