基于MF6，EMM2010和CGRF2010模型的中國大陸地殼磁異常特征

馮麗麗 王 粲 陳 斌 袁潔浩

1）中國北京100081中國地震局地球物理研究所

2）中國西寧810001青海省地震局

引言

地磁場由主磁場、巖石圈磁場、變化磁場和感應(yīng)磁場等4部分組成.從地磁場中減去主磁場，再消去變化磁場和感應(yīng)磁場，剩余部分稱為巖石圈磁場，又稱為地殼磁場或地磁異常場.巖石圈磁場的異常分布是由巖石圈物質(zhì)磁化強度的不均勻分布所引起的，具有空間分布極其復(fù)雜、而時間上卻異常穩(wěn)定的特點.在地震過程中巖體受力產(chǎn)生形變直至破裂，巖石磁性會伴隨這一過程發(fā)生變化，并表現(xiàn)為巖石圈磁場的變化.若這種變化引起的量級能夠被地面儀器觀測到，則可以利用這種效應(yīng)（地震磁效應(yīng)）獲知震源巖體的受力狀態(tài)和運動過程，進而可以根據(jù)巖體應(yīng)力、應(yīng)變狀態(tài)并結(jié)合地體條件預(yù)測地震（丁鑒海等，2011）；巖石圈磁場能夠反映地殼和上地幔的結(jié)構(gòu)特征和物質(zhì)組成，可為地球物理構(gòu)造研究提供依據(jù)（車自成等，2002；吳招才等，2011）；根據(jù)磁異常并結(jié)合地質(zhì)及地球物理資料可以圈定巖體和劃分不同巖性區(qū)，推斷斷裂、破碎帶及褶皺位置及產(chǎn)狀，并圈定和劃分成礦區(qū)（熊光楚，1994；管志寧，1997；李才明，李軍，2013）.

巖石圈磁場一直是地磁學(xué)中的重要研究對象.得益于美國、德國、前蘇聯(lián)和丹麥等國家陸續(xù)發(fā)射的多顆磁測衛(wèi)星，對全球范圍的地磁場實現(xiàn)了高密度連續(xù)觀測.現(xiàn)有大部分全球地磁模型都是基于衛(wèi)星資料建立的，如Pomme模型（Maus et al，2010）和綜合磁場（magnetic field，簡寫為MF）模型（Maus et al，2008）等.美國和英國的科研人員將衛(wèi)星、航空、海洋和地面磁測數(shù)據(jù)綜合匯編，建立了分辨率達到56km的NGDC-720模型及增強模型（enhanced magnetic model 2010，簡寫為EMM2010）（Maus，2010）.該模型是利用球諧分析法描述巖石圈磁場精度最高的模型之一.利用上述全球模型提供的球諧系數(shù)可以計算中國地區(qū)的地磁場強度.

基于實際需要，我國地磁科研人員開展了大量關(guān)于中國區(qū)域地磁場模型的研究工作.1950—2000年期間，每10年繪制出版一次中國地磁圖，2000年以后每5年繪制出版一次中國地磁圖.目前最新一代中國地磁圖為《2010.0年代中國地磁圖》，該地磁圖及與其相配套的2010.0年代中國地磁參考場（Chinese geomagnetic reference field on 2010epoch，簡寫為CGRF2010）模型描述的是2010.0年代地球主磁場和巖石圈磁場在中國地區(qū)的分布狀態(tài).

在建立地磁場模型過程中，由于地磁場觀測數(shù)據(jù)來源不同，以及建立地磁場模型所選用方法不同，可能導(dǎo)致不同地磁模型存在較為顯著的差異.因此，在開展研究工作時，根據(jù)實際需要選擇適當(dāng)?shù)牡卮艌瞿Ｐ褪种匾?本文首先分別計算了第六代綜合模型MF6及EMM2010模型在中國大陸地區(qū)內(nèi)（18°N—54°N，73°E—136°E）內(nèi)巖石圈磁場各要素，研究了垂直分量ΔZ的強度及其分布特征；然后對MF6模型、EMM2010模型及CGRF2010模型的巖石圈磁場分布特征進行了對比，指出3個模型的差異，并分析造成差異的可能原因；最后分析了各模型的優(yōu)點和不足，以便根據(jù)研究的實際需要選擇適當(dāng)?shù)哪Ｐ瓦M行計算.

1 建模理論

1.1 高斯球諧分析

高斯在1837年將球諧分析法應(yīng)用于對地磁場的數(shù)學(xué)描述中.地磁場B是一個位場，其位函數(shù)可以用一個標(biāo)量函數(shù)U來描述，即

在地表與電離層底部之間的自由空間，地磁位U滿足拉普拉斯方程，即

在球坐標(biāo)系下展開拉普拉斯方程：

則式（3）的解可寫為

式中中括號內(nèi)有兩項，第1項表示內(nèi)源場，第2項表示外源場.地球表面地磁場的絕大部分來源于地球內(nèi)部，地磁場中內(nèi)源場部分為基本磁場，其相應(yīng)的位函數(shù)為

將式（5）代入式（1）得

式中：Br＝－Z，Bθ＝－X，Bλ＝Y(jié)；X，Y和Z分別為球坐標(biāo)系下地磁場北向分量、東向分量和垂直分量；λ和θ分別為球坐標(biāo)系中經(jīng)度和余緯；a為地球參考半徑（a＝6 371.2km）；r為離開地心的徑向距離；Pmn（cosθ）為締合勒讓德函數(shù)；n為球諧級數(shù)的階，m為球諧級數(shù)的次；gmn和hmn為球諧系數(shù).

1.2 曲面樣條擬合

在磁異常擬合中，常用曲面樣條函數(shù)擬合，即

式中：x和y分別為待擬合點的空間橫坐標(biāo)和縱坐標(biāo)；xi和yi分別為第i個原始測點的橫坐標(biāo)和縱坐標(biāo)；a，b，c和Fi為曲面樣條擬合系數(shù)；為待擬合點（x，y）與第i個原始測點的距離；N為原始數(shù)據(jù)個數(shù)；ε為曲率因子，控制曲面擬合數(shù)值與原始數(shù)值的最大差值.Harder和Desmarais（1972）首次將該方法用于地表離散數(shù)據(jù)的插值，安振昌和徐元芳（1981）將該方法引入到中國地區(qū)地磁場模型計算中.

2 模型簡介

2.1 全球模型

近年來日益發(fā)展的近地表磁測（地面、海洋和航空磁測）以及衛(wèi)星磁測積累了大量的地磁場觀測數(shù)據(jù).其中，衛(wèi)星磁測手段對地磁場研究作出了突出貢獻.衛(wèi)星磁場測量的顯著優(yōu)點是能夠在地表以上數(shù)百千米的高度對地磁場進行高密度掃描式的測量，可獲得國境線以外地區(qū)及自然條件惡劣地區(qū)的地磁測量數(shù)據(jù)，彌補了其它地磁觀測手段所不能覆蓋的監(jiān)測空白.以高斯理論為基礎(chǔ)并結(jié)合大量地磁場實測資料，已建立了多種全球地磁場模型.

2.1.1 MF系列模型

MF模型是利用CHAMP衛(wèi)星數(shù)據(jù)構(gòu)建的系列模型.第1版MF1模型于2002年發(fā)布，用15—80階球諧函數(shù)代表巖石圈磁場的可見部分（Maus et al，2008）.隨著CHAMP衛(wèi)星磁測數(shù)據(jù)的不斷積累，陸續(xù)建立了新的MF系列模型.其中，MF6模型使用2004—2007年的觀測數(shù)據(jù)建立，最大截止階數(shù)提高到133階，所對應(yīng)的最小空間波長為300km.該模型是第一個基于衛(wèi)星解決海底磁條帶方向的巖石圈磁場模型，可以揭示海洋地殼的年齡結(jié)構(gòu).由MF6模型計算得到的中國大陸及其鄰區(qū)巖石圈磁場如圖1所示.可以看出，其正負(fù)異常條帶似蠕蟲狀相間分布，條帶寬度約為150—300km.

2.1.2 NGDC-720模型和EMM2010模型

NGDC-720模型是由地面、航空、海洋和衛(wèi)星磁測數(shù)據(jù)綜合匯編而成.該模型建立過程中使用了經(jīng)掃描數(shù)字化后的中國大陸地區(qū)航磁圖數(shù)據(jù).由于衛(wèi)星數(shù)據(jù)決定的長波長場更為可靠，故波長大于330km的成分由CHAMP衛(wèi)星得到的地殼磁場所代替.最后，由最小二乘法反演求出地磁位的球諧系數(shù).NGDC-720模型提供了16—720階球諧級數(shù)展開來表示地殼場，對應(yīng)于2 500—56km的空間波長范圍.720階的截斷階數(shù)相當(dāng)于角波長為30弧分，即模型分辨率為15弧分.從此意義上看，該模型具有極高的空間分辨率.

作為EMM模型的前身，NGDC-720模型將不再進行更新.以后的版本則稱為EMM2010模型.該模型包含1—720階的球諧級數(shù)，時間跨度為2000—2015年.其中2000—2010年的主磁場（1—15階）是由POMME7模型的前15階所代替.巖石圈磁場部分（16—720階）則來源于NGDC-720模型.由EMM2010模型計算得到的中國大陸及其鄰區(qū)巖石圈磁場分布如圖2所示.理論上該模型分辨率可達到56km，因此其所顯示的磁異常細(xì)節(jié)特征十分顯著，碎屑狀磁異常散布于整個區(qū)域.在大規(guī)模盆地及山脈地區(qū)，一般可出現(xiàn)較大尺度的片狀磁異常區(qū).

圖1 由MF6模型得到的中國大陸及其鄰區(qū)巖石圈磁場北向分量ΔX（a），東向分量ΔY （b）和垂直分量ΔZ（c）A—E分別為西域巖石圈、青藏巖石圈、松遼巖石圈、華北巖石圈和華南巖石圈地塊，下同F(xiàn)ig.1 The north componentΔX （a），east componentΔY （b）and vertical componentΔZ （c）of lithospheric magnetic field of Chinese mainland and its adjacent areas obtained by the model MF6 The capital letters A，B，C，D，Erepresent the Xiyu，Qinghai-Xizang（Tibet），Songliao，North China and South China lithosphere block，respectively，the same below

圖2 由EMM2010模型得到的中國大陸及其鄰區(qū)巖石圈磁場北向分量ΔX（a），東向分量ΔY （b）和垂直分量ΔZ（c）Fig.2 The north componentΔX （a），east componentΔY （b）and vertical componentΔZ （c）of lithospheric magnetic field of Chinese mainland and its adjacent areas obtained by the model EMM2010

2.2 2010.0年代中國區(qū)域地磁場曲面樣條模型

2008—2009年，中國地震局地球物理研究所承擔(dān)了中國大陸地區(qū)共851個野外地磁三分量測量點的測量工作（陳斌，2013）.這些測點的觀測數(shù)據(jù)是《2010.0年代中國地磁圖》的主要數(shù)據(jù)來源.除此之外，2002—2004年431個測點的觀測數(shù)據(jù)也參與到《2010.0年代中國地磁圖》的編制過程中.另外，《2010.0年代中國地磁圖》還使用了鄰近海島及周邊國家190個測點的歷史地磁觀測數(shù)據(jù)，以及32個國內(nèi)地磁臺和11個周邊國家地磁臺的觀測數(shù)據(jù).將全部測點的地磁觀測數(shù)據(jù)進行日變化改正和長期變化改正后，統(tǒng)一歸算到2010年1月1日午夜北京時00—03時的地磁場均值中.為降低邊界效應(yīng)，還在計算區(qū)域邊緣均勻補充了36個測點的國際地磁參考場計算值.

圖3 2010.0年代中國地磁參考場模型野外地磁測點和臺站分布Fig.3 Distribution of stations and observatories of the model CGRF2010

對上述數(shù)據(jù)直接進行曲面樣條擬合，計算了中國大陸地區(qū)地磁場，得到了2010.0年代中國地區(qū)地磁場曲面樣條模型.該數(shù)據(jù)已經(jīng)進行了日變化改正，扣除了外源場的影響.由于測量過程中避開了快速變化的外源磁場活動時段，降低了感應(yīng)磁場的影響，故建模過程中可忽略這部分能量，可以認(rèn)為擬合得到的模型中包含了主磁場和巖石圈磁場.其中主磁場部分根據(jù)2010.0年代中國地區(qū)地磁場球冠諧模型（陳斌，2013）計算得到.該模型利用8階球冠諧函數(shù)擬合中國地區(qū)851個測點的流動地磁觀測數(shù)據(jù)，得到主磁場；將其從2010.0年代中國地區(qū)地磁場曲面樣條模型中消除，即可得到2010.0年代中國地區(qū)巖石圈磁場.

圖4為由CGRF2010模型得到的中國大陸巖石圈磁場.可以看出，該區(qū)域巖石圈磁場以串珠狀磁異常為主，東部與西部存在較大差異.東部由于測點密集，顯示出大量的磁異常細(xì)節(jié)，而西部尤其是青藏高原地區(qū)由于測點稀少分辨率則較低.在CGRF2010模型中中國大陸巖石圈磁場存在若干大尺度的異常區(qū)，下文將作進一步討論.

3 中國大陸地區(qū)巖石圈磁場

圖4 由CGRF2010模型得到的中國大陸巖石圈磁場北向分量ΔX（a），東向分量ΔY （b）和垂直分量ΔZ（c）Fig.4 The north componentΔX （a），east componentΔY （b）and vertical componentΔZ （c）of lithospheric magnetic field of Chinese mainland obtained by the model CGRF2010

中國大陸位于歐亞、印度－澳大利亞和太平洋（菲律賓）三大板塊交匯部位，經(jīng)歷了漫長的地質(zhì)演化歷史，地質(zhì)構(gòu)造十分復(fù)雜，巖石圈結(jié)構(gòu)頗具特色.根據(jù)李廷棟（2006）研究可以將中國大陸劃分為5個巖石圈地塊，分別為西域巖石圈、青藏巖石圈、松遼巖石圈、華北巖石圈和華南巖石圈地塊.下面分別討論由3個模型計算得出的巖石圈磁場垂直分量ΔZ在5個巖石圈地塊內(nèi)的分布特征.由于EMM2010模型正負(fù)磁異常區(qū)分布概貌與MF6模型十分類似，因此將二者作為全球模型共同討論.

3.1 西域巖石圈地塊

西域巖石圈地塊以康西瓦斷裂帶—阿爾金斷裂帶—祁連北緣斷裂帶與青藏巖石圈地塊為分界，東部為南北構(gòu)造帶，北部和西部延伸至境外.地殼表層構(gòu)造為近東西向造山帶與構(gòu)造盆地相間排列的構(gòu)造格局，形成“三山二盆”結(jié)構(gòu)，從北到南依次為阿爾泰山、準(zhǔn)噶爾盆地、天山、塔里木盆地和昆侖山.各模型計算得出的西域巖石圈地塊ΔZ分量分布如圖5所示.可以看出，全球模型ΔZ分量山脈沿線以負(fù)磁異常為主，盆地地區(qū)多為正磁異常.具體表現(xiàn)為：

1）阿爾泰山負(fù)異常.由MF6和EMM2010模型得出的巖石圈磁場在該地區(qū)均表現(xiàn)為負(fù)磁異常，最大異常強度均在－110nT左右.異常區(qū)邊緣磁場水平梯度大，向內(nèi)逐漸趨于平緩.

2）南北天山負(fù)異常.沿南北天山出現(xiàn)呈V字型分布的負(fù)異常.在南天山南部負(fù)異常區(qū)范圍向南擴散到塔里木盆地北部，覆蓋了40°N以北的大部分地區(qū).由MF6模型得到的該地區(qū)巖石圈磁場最大負(fù)異常為－150nT（85.4°E，41.8°N）；由EMM2010模型得到的該地區(qū)巖石圈磁場最大負(fù)異常為－223nT（88°E，41°N）.

圖5 由MF6（a），EMM2010（b）和CGRF2010（c）模型得到的中國大陸西部地區(qū)巖石圈磁場垂直分量ΔZFig.5 The vertical componentΔZof lithospheric magnetic field of western Chinese mianland obtained by the models MF6（a），EMM2010（b）and CGRF2010（c）

3）準(zhǔn)噶爾盆地正異常.準(zhǔn)噶爾盆地除小部分邊緣地區(qū)以外，巖石圈磁場均表現(xiàn)為正磁異常.由MF6和EMM2010模型得到的該地區(qū)巖石圈磁場最大強度分別為87nT（84.7°E，46.4°N）和207nT（85.5°E，45.2°N）.EMM2010模型相比于 MF6模型的一個顯著特征表現(xiàn)為：在準(zhǔn)噶爾盆地中部，存在一條貫穿東西、略向南凸出的狹窄負(fù)磁異常帶，強度約為－110nT.該特征與趙俊猛等（2008）提供的航磁異常圖相符合.

4）塔里木盆地中南部正異常.塔里木盆地在40°N以南為正異常區(qū)，該區(qū)域范圍與盆地邊緣基本吻合.由MF6和EMM2010模型得到的該區(qū)域巖石圈磁場最大強度分別為171nT（87.3°E，39.9°N）和389nT（82.4°E，39.6°N）.EMM2010模型得到的巖石圈磁場在塔里木盆地南部可見一條東西向延展的狹窄負(fù)磁異常帶，幾乎貫穿了整個盆地，但其異常強度并不大，約為－20—－30nT.該結(jié)果與楊文采等（2012）的結(jié)論一致.

根據(jù)地磁場球諧理論，EMM2010模型所描述的巖石圈磁場比MF6模型包含了更多的短波成分，表明EMM2010模型所顯示的小尺度異常與淺表的磁化率分布有關(guān).CGRF2010模型ΔZ分量與全球模型在該地區(qū)表現(xiàn)出顯著的不同，具體表現(xiàn)為：

1）準(zhǔn)噶爾盆地核心區(qū)正異常.準(zhǔn)噶爾盆地的巖石圈磁場分布較為不規(guī)則，與地表構(gòu)造的關(guān)系似乎并不明確，其最大正異常值為259nT（86.4°E，45.7°N）.

2）南北天山及塔里木盆地北部以負(fù)異常為主.最大負(fù)異常為－228nT（83.6°E，41.8°N），其延展優(yōu)勢方向為東西向，在東緣被一弧形正異常斬斷，而在北部似乎向北延伸出兩條南北向負(fù)異常帶.

3）塔里木盆地中南部正異常.其分布范圍比全球模型略向北收攏，最大正異常為234nT（88.8°E，40.4°N）.雖然該區(qū)域幾個核心位置的磁異常正負(fù)性并未見大面積反轉(zhuǎn)，但其異常區(qū)分布范圍和形態(tài)均明顯有別于全球模型.尤其在北部地區(qū)，全球模型異常條帶基本沿東西向延展，而CGRF2010模型似乎以南北向條帶居多.

由于CGRF2010模型測點僅分布于境內(nèi)，阿爾泰山絕大部分地區(qū)由于缺乏磁場觀測資料，故在此不作討論.

3.2 青藏巖石圈地塊

青藏巖石圈地塊包括整個青藏高原，其北部康西瓦斷裂帶—阿爾金斷裂帶—祁連北緣斷裂帶與西域巖石圈地塊相鄰，南部喜馬拉雅主邊界斷裂帶與印巴巖石圈地塊相隔，東部龍門山造山帶與華南巖石圈地塊相接.各模型計算得到的青藏巖石圈地塊ΔZ分量分布如圖5所示.全球模型ΔZ分量在青藏巖石圈地塊的分布特征康國發(fā)等（2011）曾作過詳細(xì)討論，在此簡要說明以下兩點：

1）在MF6模型和EMM2010模型磁異常圖上，青藏高原地區(qū)巖石圈磁場總體表現(xiàn)為低值異常區(qū)，與其周邊地區(qū)（北部塔里木盆地強正異常區(qū)、西南部印度板塊團狀高值異常區(qū)）形成鮮明對比.駱遙（2013）研究表明，NGDC-720模型在青藏高原為一巨大資料空白區(qū)，導(dǎo)致其不能反映高原內(nèi)各地塊分界，并表現(xiàn)出低強度且均勻的負(fù)異常特征.

2）在青藏高原低弱的負(fù)磁異常背景下，柴達木盆地顯示出較大范圍的微弱正異常.在MF6和EMM2010模型巖石圈磁場圖上，該地區(qū)磁異常強度不超過50nT，且分布范圍僅限于柴達木盆地東部.

CGRF2010模型ΔZ分量與全球模型存在較大差異.由于可可西里無人區(qū)自然條件過于惡劣，因此缺乏有效的地面觀測數(shù)據(jù).除該地區(qū)外，西部地區(qū)地磁測點平均密度約為150km.盡管地磁測點的空間分布密度遠(yuǎn)低于衛(wèi)星磁測，但將地面磁測數(shù)據(jù)用于中國地區(qū)地磁場建模能夠真實有效地反映地表磁場.基于地表的磁場觀測相比于350—400km高度的衛(wèi)星觀測包含了更多的磁場信息，尤其是巖石圈磁場（短波長的巖石圈磁場成分會隨著高度增加而迅速衰減）.數(shù)據(jù)源的差異使得CGRF2010模型相比于衛(wèi)星模型得到的巖石圈磁場分布形態(tài)及異常強度均存在較大差異.該模型得到的青藏巖石圈地塊磁異常特征為：

1）青藏高原東北部正異常.在青藏高原東北部存在一個上千千米尺度的正異常區(qū)，其范圍覆蓋了整個柴達木盆地且向南突出，貫穿巴顏喀拉地塊中段并延伸到羌塘地塊中部.該異常強度達到231nT（37°N，91°E）.沿著該區(qū)域向東、西、南3個方向輻射有若干個異常條帶，但強度并不大.

2）青藏高原南部南北向磁異常帶.與全球模型不同，CGRF2010模型得到的巖石圈磁場在青藏高原南部存在明顯的正負(fù)異常分區(qū).該異常條帶優(yōu)勢方向為南北向，與青藏高原南部已探明的南北向構(gòu)造走向（薛典軍等，2006）一致.與此同時，CGRF2010模型得到的巖石圈磁異常強度普遍比全球模型高，其最大正異常強度達到208nT（82.6°E，32.7°N），而最大負(fù)異常也達到－178nT（80.6°E，32.3°N）.

應(yīng)指出的是，由于測點分布稀疏，監(jiān)測空區(qū)（可可西里地區(qū)）磁場模型誤差較大.而在青藏高原南部分布有若干地磁測點，從而保證了CGRF2010模型在該地區(qū)可信度較高.

3.3 松遼巖石圈地塊

松遼巖石圈地塊包括東北東部、松遼平原和大興安嶺地區(qū)南部.該地塊以中朝地臺北緣斷裂帶為界與華北巖石圈地塊相接；向西向北可能包括了蒙古東部及俄羅斯外興安嶺以南廣大地區(qū)；向東連接錫霍特地區(qū).衛(wèi)星磁異常條帶走向以ENE向為主，與控制該地區(qū)的主要構(gòu)造走向一致（圖6）.全球模型ΔZ分量在該地區(qū)的主要分布特征如下：

1）大興安嶺正異常.MF6模型巖石圈磁場圖上，從西伯利亞東南部到該區(qū)域分布一系列ENE向近似平行延伸的蠕蟲狀正負(fù)異常條帶.在中國境內(nèi)，沿大興安嶺北段存在一條寬約500km的正異常帶，強度達到150nT.在45°N以南靠近陰山—燕山構(gòu)造帶北部地區(qū)異常條帶轉(zhuǎn)為NE向，并向西延伸且逐漸向陰山—燕山構(gòu)造帶靠攏.EMM2010模型得到的巖石圈磁場在整體上符合上述特征，局部強度比MF6模型得到的異常強度要大.

2）松遼盆地正負(fù)異常交錯分布.在大興安嶺、小興安嶺及長白山環(huán)繞的松遼盆地，MF6模型巖石圈磁場顯示出西北負(fù)異常、東南正異常的格局，異常強度不大.而EMM2010模型磁異常分布十分復(fù)雜，呈碎屑狀正負(fù)異常交錯分布，意味著該地區(qū)巖石圈淺部磁化率分布不均勻程度高.值得注意的是，在正負(fù)異常錯雜分布的大背景下，在松遼盆地東北角，小興安嶺與長白山交匯處所圍的三角形地區(qū)存在范圍約為200km的正異常區(qū)，該異常區(qū)可能與該地區(qū)較深的居里面有關(guān)（胡旭芝等，2006）.

圖6 由MF6（a），EMM2010（b）和CGRF2010（c）模型得到的中國大陸東部地區(qū)巖石圈磁場垂直分量ΔZFig.6 The vertical componentΔZof lithospheric magnetic field of eastern Chinese mainland obtained by the models MF6（a），EMM2010（b）and CGRF2010（c）

3）長白山低值異常.在松遼盆地東南部的長白山地區(qū)，MF6模型得到的磁異常強度絕對值在50nT以內(nèi)，從東北到西南呈正負(fù)異常相間排列.而EMM2010模型得到的巖石圈磁場則以從東北部的三江平原地區(qū)延伸的負(fù)異常為背景，在該背景下呈現(xiàn)小尺度的正異常零星分布.

CGRF2010模型ΔZ分量在42°E以北地區(qū)，地磁測點空間分布較為均勻，平均間隔約為110km.該地區(qū)巖石圈磁異常分布呈現(xiàn)十分明顯的團狀不均勻性，磁異常強度大且區(qū)域內(nèi)變化劇烈，與同區(qū)域的衛(wèi)星磁異常特征有著顯著差別.具體分布特征為：

1）大興安嶺串珠狀負(fù)異常.在我國境內(nèi)大興安嶺南段，串珠狀負(fù)異常沿北東向分布并向西南方向延伸，最大強度約為－340nT（120°E，45.5°N），其兩側(cè)伴隨有團狀正異常，沿西北方向磁異常變化劇烈；在大興安嶺西南段，負(fù)異常面積增大、強度減弱，并逐漸轉(zhuǎn)向東北偏東方向.

2）松遼盆地正異常.松遼盆地整體表現(xiàn)為數(shù)個正異常團，北部正異常團強度最大，超過500nT，且該正異常范圍在盆地東北部向長白山北緣突出，這一細(xì)節(jié)特征與衛(wèi)星模型磁異常圖是一致的.

3）長白山負(fù)異常.除北部小區(qū)域外，長白山脈沿線以負(fù)異常為主，強度約為150nT.

3.4 華北巖石圈地塊

華北巖石圈地塊北以中朝地臺北緣斷裂帶為界，與松遼巖石圈地塊相鄰；南以商丹斷裂帶—磨子潭斷裂帶—舟山斷裂帶為界，與華南巖石圈地塊相接；西界為南北構(gòu)造帶；向東包括黃海.在大地構(gòu)造上主要屬于中朝地臺范疇.該區(qū)域地殼表層構(gòu)造呈NNE向隆起與坳陷相間排列的構(gòu)造格局.全球模型ΔZ分量顯示華北巖石圈地塊磁異常與其東部松遼巖石圈地塊存在較為明顯差異（圖6）.MF6模型得到的巖石圈磁異常尺度較小，正負(fù)異常形態(tài)多呈團狀，條帶狀不明顯；而EMM2010模型得到的巖石圈磁場則可以分辨出較為清晰的NE向正負(fù)異常條帶，但異常條帶并不沿構(gòu)造邊界延展.具體區(qū)域特征如下：

1）鄂爾多斯NE向正負(fù)異常條帶.在MF6模型得到的巖石圈磁場圖上該特征并不明顯，但在EMM2010模型得到的巖石圈磁場圖上卻顯示了4條清晰的NE向正負(fù)異常條帶，其分布方向與構(gòu)造走向存在一個夾角.

2）華北細(xì)碎的巖石圈磁場.在太行山以西，相對于鄂爾多斯，MF6模型得到的磁異常弱且平緩，焦點異常值約為70nT（116.1°E，34.8°N）.而EMM2010模型得到的巖石圈磁場則顯示了非常細(xì)碎的正負(fù)異常細(xì)節(jié)，除南部邊界附近的狹窄異常帶以外，并未見其它顯著的規(guī)律性分布，焦點異常值約為370nT（117.8°E，37.8°N）.

3）南黃海陸上負(fù)異常.在MF6模型得到的巖石圈磁場圖上，郯廬斷裂帶南段（華北東部）與磁異常零值線基本重合.其西側(cè)磁異常為正異常，東側(cè)南黃海陸地部分為負(fù)異常，最大異常值約為－75nT（119.7°E，32.5°N）.EMM2010模型在負(fù)異常背景下顯示了一些細(xì)碎的正異常，似乎存在兩條分別為NE向和ENE向負(fù)異常帶從陸地向黃海延伸出去.陸上最大負(fù)異常值約為－190nT（120.4°E，32°N）.

CGRF2010模型ΔZ分量在該地區(qū)測點準(zhǔn)均勻分布，點間距為70—100km.從該模型得到的巖石圈磁場圖來看，具有以下幾個顯著的區(qū)域特征：

1）鄂爾多斯NE向磁異常帶.該地區(qū)CGRF2010模型與EMM2010模型得到的巖石圈磁場類似.所不同的是，CGRF2010模型得到的該地區(qū)正異常條帶寬于負(fù)異常條帶，且NE向分布的正異常帶有跨過構(gòu)造邊界向東延伸的趨勢.該地區(qū)最大異常值約為300nT（111.2°E，38.2°N）.

2）華北團狀磁異常.該地區(qū)CGRF2010模型得到的巖石圈磁場分布與衛(wèi)星模型差異很大.泰山附近負(fù)異常顯著，異常值達到－347nT（117°E，36.2°N）；北京及河北地區(qū)主要為正異常，最大異常值為373nT（116.2°E，39.4°N）；南部未見明顯沿構(gòu)造邊界分布的磁異常條帶.

3）南黃海陸上團狀磁異常.該地區(qū)范圍不大，但CGRF2010模型顯示的巖石圈磁場分布卻存在異常強度大、形態(tài)復(fù)雜的特點.可能是由于海上無實測數(shù)據(jù)，導(dǎo)致海陸交界地區(qū)模型結(jié)果準(zhǔn)確性降低.該模型與衛(wèi)星模型得到的巖石圈磁場的顯著區(qū)別在于該地區(qū)存在明顯的正異常團，例如：南京附近200km尺度的正異常團，其最大強度為265nT（118.8°E，31.8°N）；山東省境內(nèi)郯蘆斷裂帶東側(cè)的正異常團.

3.5 華南巖石圈地塊

華南巖石圈地塊北以商丹斷裂帶—磨子潭斷裂帶—舟山斷裂帶為界，與華北巖石圈地塊相隔；南以濱海斷裂帶為界，與南海巖石圈地塊相鄰；西抵川滇南北向構(gòu)造帶；向東包括東海和琉球島弧隆起帶.全球模型ΔZ分量沿北部邊界可見狹窄的磁異常條帶，邊界兩側(cè)磁異常特征迥異.與華北巖石圈地塊的中小尺度團狀磁異常不同，華南巖石圈地塊除四川盆地外，其余地區(qū)磁異常強度低，變化平緩，且以負(fù)異常為主（圖6）.具體區(qū)域特征為：

1）四川盆地正異常區(qū)及其北部負(fù)異常區(qū).四川盆地及其北部地區(qū)在中國地區(qū)衛(wèi)星磁異常圖上十分顯眼，其原因在于該地區(qū)出現(xiàn)一對尺度約為800km×300km、近東西向延展的正負(fù)異常對.其中四川盆地以強正異常為主，由MF6模型得到的巖石圈磁場最大磁異常為189nT，而EMM2010模型得到的最大磁異常達到457nT（106.4°E，30.8°N）；四川盆地北部巖石圈磁場強度向北迅速衰減，零值線沿盆地北緣展布；繼續(xù)向北沿大巴山脈以及漢中盆地存在一個與該正異常區(qū)尺度相當(dāng)?shù)呢?fù)磁異常區(qū)，由MF6和EMM2010模型得到的巖石圈磁場強度最大分別為－90nT（107°E，32.7°N）和－205nT（107.4°E，32.6°N）.

2）華南巖石圈地塊其它區(qū)域巖石圈磁場低弱平緩.由MF6模型得到的巖石圈磁場顯示華南腹地為廣泛寬緩的負(fù)值區(qū)，而一條低值正異常帶沿華南地塊北部邊界、東海、臺灣海峽及西部邊界綿延不斷將整個負(fù)異常區(qū)包圍；EMM2010模型與MF6模型得到的巖石圈磁場的區(qū)別在于，在大的區(qū)域異常背景之下，有大量細(xì)小的正負(fù)異常碎片散落各處.除西北部的四川盆地、大巴山地區(qū)以外，華南巖石圈地塊其它區(qū)域顯示了低弱的巖石圈磁場值，絕大部分地區(qū)MF6模型得到的巖石圈磁場強度在50nT以下，EMM2010模型得到的巖石圈磁場強度在100nT以下.

CGRF2010模型ΔZ分量在四川盆地為正異常區(qū)，其北部為負(fù)磁異常區(qū)，該形態(tài)特征與MF6及EMM2010模型得到的巖石圈磁場相近.但北部負(fù)異常區(qū)南北略微收緊呈帶狀，最大磁異常為－367nT（103.6°E，31.8°N）；南部正異常區(qū)幾乎完整地覆蓋了整個四川盆地，磁場強度為395nT（106°E，30.8°N）.在華夏復(fù)合陸內(nèi)造山區(qū)（張國偉等，2013）則主要表現(xiàn)為正異常，除浙江麗水附近的強正異常團最大值達到567nT（120.1°E，28°N）以外，其余地區(qū)均在150nT以下.尤其值得指出的是：在東南沿?；鸪蓭r帶為正異常；在正異常背景下有小范圍負(fù)異常區(qū)，如杭州南部及武夷山沿線；江南造山帶為斷續(xù)低弱的負(fù)異常區(qū)；揚子地塊南部以負(fù)異常為主，強度在－110nT以下.

4 討論與結(jié)論

本文通過對MF6模型、EMM2010模型及CGRF2010模型分別計算的中國大陸地區(qū)巖石圈磁場ΔZ分量分布特征進行分析，得到以下結(jié)論：

1）3個模型ΔZ分量在準(zhǔn)噶爾盆地、塔里木盆地、柴達木盆地及四川盆地等地區(qū)均為正值，而在天山、大巴山等地區(qū)則為負(fù)值.異常強度CGRF2010模型最大，EMM2010次之，MF6最小，說明利用不同數(shù)據(jù)源計算得到的模型所包含的巖石圈磁場成分不盡相同.CGRF2010模型主要使用地面磁測數(shù)據(jù)，最大限度地保留了地表以下巖石所產(chǎn)生的磁場信息；EMM2010模型對航空磁場數(shù)據(jù)的應(yīng)用以及高達720階的截斷階數(shù)，使得其相對于MF6模型包含更多的磁場信息.3種模型得到的巖石圈磁場異常區(qū)范圍大致吻合，僅細(xì)節(jié)略有不同.巖石圈磁場ΔZ分量分布在華北巖石圈地塊與華南巖石圈地塊吻合程度高；在青藏高原南部、西域巖石圈地塊北部及松遼巖石圈地塊，全球模型與CGRF2010模型存在顯著差別.

2）MF6模型能夠顯示出大尺度構(gòu)造單元所反映的巖石圈磁場輪廓，如塔里木盆地和四川盆地正異常，大興安嶺—陰山、南北天山和大巴山負(fù)異常等.由于該模型是由衛(wèi)星數(shù)據(jù)建立，因此利用該模型能夠計算得到中國及周邊地區(qū)的巖石圈磁場強度，對于研究跨國界的構(gòu)造及地震活動性具有重要意義.但由于衛(wèi)星飛行高度為地表以上數(shù)百千米，不能對地磁場的短波成分形成有效觀測，因此僅僅基于衛(wèi)星數(shù)據(jù)建立的MF6模型其分辨率較低，不能滿足對中小尺度構(gòu)造的研究.

3）EMM2010模型顯示的磁異常背景與MF6模型相似，但在細(xì)節(jié)上卻不盡相同.例如塔里木盆地和四川盆地正異常，南北天山、大巴山和大興安嶺—陰山負(fù)異常等大尺度異常區(qū)，在主體異常背景之下可見小范圍細(xì)碎的相反性質(zhì)的磁異常分布，體現(xiàn)了地表淺層磁性物質(zhì)分布的復(fù)雜性.尤其在松遼巖石圈地塊和華北巖石圈地塊，該模型巖石圈磁場細(xì)節(jié)特征顯著，可能代表著該地區(qū)淺表磁性物質(zhì)分布相當(dāng)復(fù)雜.

由于整合了全球的地面、航空、海洋及衛(wèi)星磁測數(shù)據(jù)并將球諧函數(shù)的截斷階數(shù)提高到720階，故EMM2010模型能夠顯示全球范圍的巖石圈磁場細(xì)節(jié).高達56km的空間分辨率主要依托稀疏的地面臺站觀測數(shù)據(jù)及航空磁測數(shù)據(jù).但在青藏高原等航空磁測空區(qū)EMM2010模型的準(zhǔn)確性還有待進一步論證.

4）CGRF2010模型得到的巖石圈磁場與前兩者有較大區(qū)別，尤其在中國西部地區(qū)更為明顯.青藏高原地區(qū)存在較強磁異常區(qū)，例如柴達木盆地，以及在青藏高原南部存在若干個南北向展布的正異常區(qū)，與已探明的南北向構(gòu)造特征相符.而在西域巖石圈地塊北部地區(qū)，CGRF2010模型地磁異常則表現(xiàn)出強度強、梯度大的特點，且有南北向展布的趨勢.東部地區(qū)由于地磁測點分布密集，保證了該模型的分辨率及可靠性.相鄰巖石圈地塊顯示出了獨特的磁異常特征：松遼巖石圈地塊的團狀強磁異常，華北巖石圈地塊細(xì)碎的正負(fù)異常，華南巖石圈地塊一對大規(guī)模磁異常區(qū)以及周邊相對低弱而均勻的磁異常區(qū).

CGRF2010模型是基于中國大陸地面磁場觀測數(shù)據(jù)建立的，能在最大程度上準(zhǔn)確真實地反映了中國大陸地區(qū)地磁場的空間分布特征，尤其是地殼淺層磁性物質(zhì)產(chǎn)生的磁場成分保存得比較完整.但由于受地理環(huán)境及自然條件制約，地面磁場觀測很難覆蓋所有地區(qū)，如可可西里無人區(qū)、國境線以外等地區(qū).由于缺乏有效的數(shù)據(jù)控制，地磁場模型在邊界附近可能存在畸變，且不能擬合境外的巖石圈磁場.在研究跨國界的地震活動和大地構(gòu)造時，需要能夠描述邊境附近巖石圈磁場的地磁模型.為滿足上述要求，值得探索的方法一是補充監(jiān)測空白區(qū)的地面地磁測點；二是整合多源數(shù)據(jù)，如衛(wèi)星磁測數(shù)據(jù)與地面磁測數(shù)據(jù)的整合，能夠從一定程度上增強對地面監(jiān)測空區(qū)及國界周邊地區(qū)的監(jiān)測能力.

安振昌，徐元芳.1981.局部地區(qū)高空磁場的計算方法［J］.空間科學(xué)學(xué)報，1（1）：68-73.

An Z C，Xu Y F.1981.Methods of computation of geomagnetic field at greater altitude in a local region［J］.Chinese Journal of Space Sciences，1（1）：68-73（in Chinese）.

車自成，羅金海，劉良.2002.中國及其鄰區(qū)區(qū)域大地構(gòu)造學(xué)［M］.北京：科學(xué)出版社：5-6.

Che Z C，Luo J H，Liu L.2002.Geotectonics of China and Its Adjacent Areas［M］.Beijing：Science Press：5-6（in Chinese）.

陳斌.2013.區(qū)域地磁場時空變化特征及模型研究［D］.北京：中國地震局地球物理研究所：52-60.

Chen B.2013.Regional Geomagnetic Field Temporal-Spatial Variation Characteristics and Models Study［D］.Beijing：Institute of Geophysics，China Earthquake Administration：52-60（in Chinese）.

丁鑒海，盧振業(yè)，余素榮.2011.地震地磁學(xué)概論［M］.合肥：中國科學(xué)技術(shù)大學(xué)出版社：226-370.

Ding J H，Lu Z Y，Yu S R.2011.A Brief Treatise on Seismomagnetism［M］.Hefei：University of Science and Technology of China Press：226-370（in Chinese）.

管志寧.1997.我國磁法勘探的研究與進展［J］.地球物理學(xué)報，40（增刊）：299-307.

Guan Z N.1997.Researches and progresses of magnetic prospecting in China［J］.Acta Geophysica Sinica，40（Suppl）：299-307（in Chinese）.

胡旭芝，徐鳴潔，謝曉安，王良書，張慶龍，劉紹文，解國愛，馮昌閣.2006.中國東北地區(qū)航磁特征及居里面分析［J］.地球物理學(xué)報，49（6）：1674-1681.

Hu X Z，Xu M J，Xie X A，Wang L S，Zhang Q L，Liu S W，Xie G A，F(xiàn)eng C G.2006.A characteristic analysis of aeromagnetic anomalies and Curie point isotherms in northeast China［J］.Chinese Journal of Geophysics，49（6）：1674-1681（in Chinese）.

康國發(fā)，高國明，白春華，韶丹，馮麗麗.2011.青藏高原及鄰區(qū)的地殼磁異常特征與區(qū)域構(gòu)造［J］.中國科學(xué)：地球科學(xué)，41（11）：1577-1585.

Kang G F，Gao G M，Bai C H，Shao D，F(xiàn)eng L L.2012.Characteristics of the crustal magnetic anomaly and regional tectonics in the Qinghai-Tibet Plateau and the adjacent areas［J］.Science China：Earth Sciences，55（6）：1028-1036.

李才明，李軍.2013.重磁勘探原理與方法［M］.北京：科學(xué)出版社：256-285.

Li C M，Li J.2013.The Principle and Method of Gravity and Magnetic Exploration［M］.Beijing：Science Press：256-285（in Chinese）.

李廷棟.2006.中國巖石圈構(gòu)造單元［J］.中國地質(zhì)，33（4）：700-710.

Li T D.2006.Lithospheric tectonic units of China［J］.Geology in China，33（4）：700-710（in Chinese）.

駱遙.2013.青藏高原在巖石圈磁場模型中的巨大空白：同“青藏高原及鄰區(qū)的地殼磁異常特征與區(qū)域構(gòu)造”作者商榷［J］.中國科學(xué)：地球科學(xué)，43（8）：1376-1378.

Luo Y.2013.A huge gap of the lithospheric magnetic field model in the Qinghai-Tibet Plateau：Comment on“characteristics of the crustal magnetic anomaly and regional tectonics in the Qinghai-Tibet Plateau and the adjacent areas”［J］.Science China：Earth Sciences，56（9）：1623-1626.

吳招才，高金耀，李家彪，張濤，沈中延，楊春國.2011.南海北部磁異常特征及對前新生代構(gòu)造的指示［J］.地球物理學(xué)報，54（12）：3292-3302.

Wu Z C，Gao J Y，Li J B，Zhang T，Shen Z Y，Yang C G.2011.The characteristics of magnetic anomalies：Implications for Pre-Cenozoic tectonics of the northern South China Sea［J］.Chinese Journal of Geophysics，54（12）：3292-3302（in Chinese）.

熊光楚.1994.金屬礦區(qū)磁法勘探的進展與展望［J］.地球物理學(xué)報，37（1）：437-443.

Xiong G C.1994.Development of magnetic exploration searching for metallic ore deposits in China［J］.Chinese Journal of Geophysics，37（1）：437-443（in Chinese）.

薛典軍，姜枚，吳良士，李慶慶，汪興旺.2006.青藏高原區(qū)域重磁異常的東西向分區(qū)及其構(gòu)造地質(zhì)特征［J］.中國地質(zhì)，33（4）：912-919.

Xue D J，Jiang M，Wu L S，Li Q Q，Wang X W.2006.East-west division of regional gravity and magnetic anomalies on the Qinghai-Tibet Plateau and its tectonic features［J］.Geology in China，33（4）：912-919（in Chinese）.

楊文采，王家林，鐘慧智，陳冰.2012.塔里木盆地航磁場分析與磁源體結(jié)構(gòu)［J］.地球物理學(xué)報，55（4）：1278-1287.

Yang W C，Wang J L，Zhong H Z，Chen B.2012.Analysis of regional magnetic field and source structure in Tarim Basin［J］.Chinese Journal of Geophysics，55（4）：1278-1287（in Chinese）.

張國偉，郭安林，王岳軍，李三忠，董云鵬，劉少峰，何登發(fā)，程順有，魯如魁，姚安平.2013.中國華南大陸構(gòu)造與問題［J］.中國科學(xué)：地球科學(xué)，43（10）：1553-1582.

Zhang G W，Guo A L，Wang Y J，Li S Z，Dong Y P，Liu S F，He D F，Cheng S Y，Lu R K，Yao A P.2013.Tectonics of South China continent and its implications［J］.Science China：Earth Sciences，56（11）：1804-1828.

趙俊猛，盧造勛，姚長利，黎益仕，劉占坡.2008.準(zhǔn)噶爾盆地基底斷裂的重磁學(xué)研究［J］.地震地質(zhì)，30（1）：132-143.

Zhao J M，Lu Z X，Yao C L，Li Y S，Liu Z P.2008.Crystalline basement faults in Junggar Basin determined by gravity and geomagnetism［J］.Seismology and Geology，30（1）：132-143（in Chinese）.

Harder R L，Desmarais R N.1972.Interpolation using surface splines［J］.J Aircr，9（2）：189-191.

Maus S，Yin F，Lühr H，Manoj C，Rother M，Rauberg J，Michaelis I，Stolle C，Müller R D.2008.Resolution of direction of oceanic magnetic lineations by the sixth-generation lithospheric magnetic field model from CHAMP satellite magnetic measurements［J］.Geochem Geophys Geosyst，9（7）：Q07021.doi：10.1029／2008GC001949.

Maus S.2010.An ellipsoidal harmonic representation of Earth’s lithospheric magnetic field to degree and order 720［J］.Geochem Geophys Geosyst，11（6）：Q06015.doi：10.1029／2010GC003026.

Maus S，Manoj C，Rauberg J，Michaelis I，Lühr H.2010.NOAA／NGDC candidate models for the 11th generation International Geomagnetic Reference Field and the concurrent release of the 6th generation Pomme magnetic model［J］.Earth Planets Space，62（10）：729-735.