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

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

0 引言

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

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

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

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

1 基本原理

地磁場(chǎng)的磁位勢(shì)滿足高斯理論，磁位勢(shì)U的表達(dá)式為

(1)

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

(2)

(3)

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

1.1 地殼磁場(chǎng)的反演

1.2 ML地磁能譜

(4)

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

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

1.3 線性擬合趨勢(shì)項(xiàng)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

表1 MF系列模型簡(jiǎn)介Table 1 Introduction to series of MF models

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

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

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

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