中國海陸1∶500萬莫霍面深度圖及其所反映的地質內涵

郝天珧，胡衛(wèi)劍，邢健，2，胡立天，2，徐亞，秦靜欣，3，劉少華，雷受旻

1中國科學院地質與地球物理研究所，中國科學院油氣資源研究重點實驗室，北京 100029

2中國科學院大學，北京 100049

3中國科學院遙感與數字地球研究所，北京 100101

4國土資源部青島海洋地質研究所，青島 266071

1 引言

自國際巖石圈研究計劃（ILP）開展以來，人們越來越多地認識到對莫霍面深度分布特征的研究是巖石圈研究計劃中的重要組成部分，同樣涉及到地球科學的一些本質問題，應給予充分的重視.近年來，有關地殼厚度分布的研究成果不斷發(fā)表，獲得了許多新認識，對于全球巖石圈構造特征及其深部過程的研究有重要的推動作用.

中國大陸及海域位于歐亞、印度與太平洋三大板塊的交匯部位，是地球動力學研究的重點地區(qū)，因此，中國海陸莫霍面深度分布特征的研究一直也是科研人員關注的熱點.近年來，有不少學者依據對中國地殼厚度的研究結果編制了各種的圖件（劉光鼎，1992，1993；曾融生等，1995；朱介壽，1996；Li and Mooney，1998；方劍，1999；滕吉文等，2002；高星等，2005；Li et al.，2006；黃建平等，2006），其中最有代表性的為劉光鼎等（1992，1993）第一次將我國海、陸結合在一起編繪的地球物理系列圖中，將海、陸置于統(tǒng)一構造演化體系中進行研究時所編繪的1∶1000萬的莫霍面深度分布圖.但美中不足的是該圖雖對海陸結合帶給予了足夠的重視，卻未能完全覆蓋整個中國大陸地區(qū).此外，近年來很多根據人工地震深部探測結果結合重力資料的反演結果也有發(fā)表，如朱介壽等（1996）編制的1∶1200萬莫霍面深度分布圖，以中國大陸地區(qū)的人工地震測深剖面為基礎，在全國范圍內取得400多個莫霍面深度值.曾融生等（1995）利用20世紀90年代以前的地震測深結果，編制中國大陸莫霍面深度圖，西部地區(qū)采用了面波和其他地球物理資料.為了容易與大地構造相比較還采用了分區(qū)的辦法.Li和 Mooney（1998）以及Li等（2006）總結了自1958年以來的由中國地震局地球物理勘探中心等單位完成的深地震探測剖面（DDS）約90條，根據地震測深數據繪制了中國大陸地區(qū)地殼結構示意圖.滕吉文等（2002）對東亞地區(qū)莫霍面深度分布開展了研究并編繪了1∶2500萬的地殼厚度圖，高星等（2005）則利用中國及鄰近地區(qū)數字地震臺網提供的三分量地震波形記錄，應用轉換函數及快速模擬退火算法對中國及鄰區(qū)61個地震臺站下的地殼橫波速度結構進行了反演.結合已發(fā)表的人工地震測深資料獲得了該地區(qū)地殼厚度分布.這些工作為我國地殼深部結構的研究打下了堅實的基礎.但由于這些研究大多局限在中國大陸內部或向海域延伸有限，對于探索中國大陸宏觀構造格架的延伸以及海陸相互作用關系等問題而言，仍存在一定的局限.因此，有必要在對我國50年來在中國海域與陸區(qū)的地球物理調查成果，特別是近年來一系列專項調查獲得的最新資料進行充分研究的基礎上，編制中國海陸1∶500萬莫霍面深度圖.并聚焦歐亞板塊、太平洋板塊以及印度洋板塊之間的相互作用及其在中國大陸與海區(qū)所引發(fā)的構造效應，以中國大陸板塊邊緣和板塊內部的構造演化及其宏觀格架為重點，為中國大陸邊緣和中國海動力學等方面的研究提供基礎圖件.

本圖的編圖范圍為：0°—55°N，（71.5°E，40°N）—（135°E，40°N），采用蘭伯特投影，大地基 準面：WGS84，原點經緯度：105°E，0°，標準緯度：15°N、40°N，比例尺：1∶500萬.

中國海陸1∶500萬莫霍面深度圖采用分區(qū)編繪的方法，以中國海陸及鄰區(qū)布格重力異常與空間重力異常數據為基礎，在收集、整理我國近50年來莫霍面深度調查成果以及近年來一系列專項調查獲得的最新資料基礎上，對重力異常數據進行合理分區(qū)，然后通過完全布格改正、格萊尼改正、區(qū)域場分離以及自適應模擬退火界面反演等方法，得到各分區(qū)莫霍面深度，再進行合理拼接，得到中國海陸莫霍面深度初部結果.同時，引入中國海陸及鄰區(qū)人工地震測深、地學斷面、海底地震儀地震測深等各類地球物理調查剖面120條，通過統(tǒng)一的數字化、投影、檢驗、校正等一系列處理后，得到了用于約束莫霍面深度反演的控制點共計8948個.在對初步結果進行精度評價和實測數據點檢驗后，經計算機修改，調色定稿.

2 重力數據

編圖中用于反演莫霍面深度的原始重力數據在中國陸區(qū)及鄰域為布格重力異常數據，中國海區(qū)及鄰域為空間重力異常數據，其數據分布情況見表1和表2.

表1 中國海陸莫霍面深度圖中采用的中國大陸地區(qū)重力異常數據信息Table 1 Information of gravity anomalies in the mainland of China used for preparation of the Moho depth map in the land and seas of China and adjacent areas

表2 本文收集的鄰國或地區(qū)重力資料網度一覽表Table 2 Grid spacing of gravity data collected from neighboring countries and regions

陸區(qū)重力資料是地礦部門多年開展區(qū)域重力調查所實測的數據，按照驗收時給出的布格重力異常總精度來看，其精度均優(yōu)于±2.0×10－5m·s－2.圖幅內海區(qū)采用美國Scripps研究所發(fā)布衛(wèi)星測高重力數據，經緯度網格數據2′×2′，數據精度±（4～6）×10－5m·s－2.

3 編圖方法與指導思想

本次編圖的指導思想是：以板塊構造活動論為指導，關注、聚焦歐亞板塊、太平洋板塊以及印度洋板塊之間的相互作用及其在中國大陸與海區(qū)所引發(fā)的構造效應，以中國大陸板塊邊緣和板塊內部的宏觀構造格架為重點，編制中國海陸莫霍面深度圖，為中國大陸邊緣和中國海的動力學演化等研究提供基礎圖件.

本次編圖采用了劉光鼎（1992）提出的方法，即以板塊構造活動論為指導，以巖石物性為連接地質與地球物理之間的紐帶，以地質模型為跨越地球物理場的定性解釋與定量計算之間的橋梁；采用地質與地球物理、正演與反演、定性解釋與定量計算相結合，在大量高精度數據作為先驗信息約束下，通過信息的多次反饋，反復修改地質模型，逐步加深對地質、地球物理信息的認識，有效地克服反演問題的多解性，使解釋結果漸趨合理、完善和精確.本圖編圖流程如圖1所示.

圖1 中國海陸莫霍面深度圖編圖流程Fig.1 Compilation flow of the Moho depth map in the land and seas of China and adjacent areas

4 莫霍面深度反演與偏差分析

4.1 莫霍面深度反演方法與分區(qū)

本圖范圍較大，從大陸到大洋，地形起伏劇烈，跨越了多個構造單元，構造面貌十分復雜.假若全圖均采用單一的界面反演方法和相同的反演系數，往往顧此失彼，很難獲得最佳的反演效果.因此，在本次莫霍面深度反演時，依據中國海陸及鄰區(qū)重力異常及構造單元特征，對圖幅區(qū)進行了分區(qū)，將陸區(qū)分為10個分區(qū)，海區(qū)分為6個分區(qū)（圖2）.對不同分區(qū)進行了方法優(yōu)化，選擇適合該區(qū)的重力場分離、界面反演方法和相應的反演系數（表3），最終得到中國海陸1∶500萬莫霍面深度分布圖（圖3）.

4.2 控制點數據

為了提高反演結果的精度，本次編圖引入了以人工地震等實測剖面數據為主的控制點數據8948個，對反演結果進行校正，以提高可靠性.這樣做的效果比較明顯.以南海地區(qū)為例，圖4是南海地區(qū)經過控制點校正前后的莫霍面深度結果的比較.由于在南海北部陸緣地區(qū)引入了多條OBS測深剖面數據的控制，校正后的莫霍面深度分布更為合理，與該區(qū)實際地震探測結果更為接近.

本次編圖共收集和整理了中國海陸及鄰區(qū)人工地震測深、地學斷面、海底地震儀測深等各種實測剖面120條，莫霍面深度與地殼厚度圖件18余幅.通過統(tǒng)一的數字化、投影、檢驗、校正等一系列處理，得到了以人工地震等實測剖面數據為主的控制點數據8948個，為提高中國海陸及鄰區(qū)莫霍面深度圖的反演精度提供了相應的數據.圖5為中國海陸1∶500萬莫霍面深度圖控制剖面位置和部分控制點信息（表4）.

圖2 中國海陸1∶500萬莫霍面深度反演分區(qū)Fig.2 Region division for inversion of the Moho depth（1∶5000000）in the land and seas of China and adjacent areas

圖3 中國海陸1∶500萬莫霍面深度分布Fig.3 Contours of the Moho depth（1∶5000000）in the land and seas of China and adjacent areas

圖4 南海地區(qū)莫霍面深度反演結果（a）控制點數據校正前；（b）控制點數據校正后.圖中黑色圓點為OBS實測剖面控制點位置，黑色粗實線所示為圖8中剖面II′的位置.Fig.4 Inversion results of the Moho depth in the South China Sea（a）Before control point correction；（b）After control point correction.Black dots are the locations of control points from OBS profiles.Black thick solid line is Profile II′shown in Fig.8.

表3 各分區(qū)重力場分離方法和莫霍面深度反演方法優(yōu)選Table 3 Gravity field separation and inversion method optimization in each sub－region

控制點的使用：首先使用地震測深等得到的部分控制點計算各分區(qū)（圖2）莫霍面深度的平均深度，這個平均深度值作為一個重要參數參與莫霍面深度反演計算.在獲得各分區(qū)莫霍面深度后，任意抽取各分區(qū)若干條剖面上莫霍面深度反演結果與地震實測控制點結果對比并進行偏差分析，對于偏差較大地區(qū)，采用小范圍重新進行位場分離和平均深度參數選取并重新反演來進行校正，使莫霍面深度分布更為合理，與該區(qū)實際地震探測結果更為接近.

4.3 偏差分析

為了評價反演結果的可靠性，本文選擇了若干條剖面對莫霍面深度反演結果和控制點資料進行偏差分析，以平均偏差和標準偏差作為衡量反演結果的可靠性和準確度的標準.

將莫霍面深度的反演結果與某一點的實際地震測深結果的差值定義為偏差，用Xi表示，即Xi＝X反演－X地震測深.

將所有偏差的算數平均值定義為偏差平均值，用表示.

表4 中國海陸1∶500萬莫霍面深度圖部分控制點信息Table 4 Part information of control points on the Moho depth map（1∶5000000）in the land and seas of China and adjacent areas

平均偏差是指偏差與偏差平均值之差的算數平均值.其表達式為：

標準偏差又稱為均方根偏差，以σ表示，其表達式為：

以長江三峽剖面為例進行偏差分析（表5），圖6給出沿該剖面任意抽取若干點的莫霍面深度反演結果（表5）與Zhang等（2009）發(fā)表的地震實測結果對比，偏差的平均偏差1.02km，標準偏差1.52km，表明莫霍面深度的反演結果比較接近真實地質情況.本文采用同樣的方法對圖幅內各分區(qū)莫霍面深度反演結果進行了偏差分析（參見表6），最大標準偏差2.99km，最大平均偏差1.86km，主要集中在中、西部地區(qū).

通過對比可知，反演偏差與莫霍面深度的實際起伏程度有關.一般而言，莫霍面深度起伏相對平緩的地區(qū)莫霍面深度平均偏差較小，而莫霍面深度陡變大的地區(qū)平均偏差較大.比如我國中、西部的莫霍面的實際起伏幅度較大，而重力反演結果的偏差也較大.特別是在西部地區(qū)，陸殼受到嚴重的擠壓和改造，反演結果的平均偏差值可達1.86km.對于海區(qū)，黃、渤海地區(qū)莫霍面分布總體的趨勢較為平緩，偏差也相對較小.而東海的溝－?。梵w系地殼異常區(qū)，南海、菲律賓海、安達曼海等地區(qū)均涉及到俯沖帶，莫霍梯度陡變，橫向變化劇烈，因此這幾個區(qū)域的偏差值相對較大.一般而言，利用重力數據反演得到的莫霍面是一種“趨勢面”，可以近似地用連續(xù)曲面表示.但板塊俯沖帶地區(qū)（特別是洋－陸強烈作用地區(qū)），地質結構與其他地區(qū)顯著不同，莫霍面不能近似用連續(xù)曲面表示，這正是在這些地區(qū)常規(guī)重力反演方法無法準確地反映莫霍面深度分布情況的原因.為此，本圖特在莫霍面深度分布圖上對俯沖帶地區(qū)標出了震源位置、震級、俯沖帶位置等信息，以示這些地區(qū)的特殊性并供使用者參考.

圖5 中國海陸1∶500萬莫霍面深度圖控制剖面位置Fig.5 Location of control profiles on the Moho depth map（1∶5000000）in the land and seas of China and adjacent areas

表5 長江三峽地區(qū)深部結構剖面偏差分析表Table 5 Errors of the profile of the deep structure in the area of Three Gorges，Yangtze

表6 各分區(qū)莫霍面深度反演結果偏差Table 6 Deviations of inversion results of the Moho depth in each region

5 莫霍面深度分布特征

圖6 （a）長江三峽剖面位置（黑色粗實線所示），（b）地震測深給出的速度結構（據Zhang et al.，2009），（c）長江三峽深部結構剖面中地震測深與重力反演得到的莫霍面深度對比Fig.6 （a）Location of Three Gorges，Yangtze.（b）Velocity structure from seismic sounding（after Zhang et al.，2009）.（c）Comparison of the Moho depths from seismic sounding and gravity inversion along the profile of deep structure in Three Gorges，Yangtze

以莫霍面深度分布特征和地殼性質為基礎，結合重力異常及塊體構造差異等，本次編圖共劃分了12個莫霍面深度梯級帶、六大分區(qū)和24個二級分區(qū)（表7、表8和圖7）.其中一級分區(qū)主要以地殼性質的不同劃分為厚型地殼分區(qū)、正常型地殼分區(qū)、減薄型地殼分區(qū)、溝－弧－盆體系地殼異常分區(qū)、過渡型地殼分區(qū)和大洋型地殼分區(qū)；二級分區(qū)更側重于莫霍面深度分布特征，并參考了重力異常及塊體構造差異特征等.

5.1 莫霍面深度梯級帶

（1）大興安嶺—太行山—武陵山梯級帶

該梯級帶長度約為4400km，呈NNE向.莫霍面深度等值線梯度變化較強，梯級帶從西到東約200km寬度內莫霍面深度變化為42～34km，抬升8km，是減薄型陸殼和正常陸殼的分界帶，該梯級帶在中國境外仍然向北延伸，向南與中南半島東緣梯級帶相連.

表7 中國海陸及鄰區(qū)莫霍面深度梯級帶Table 7 Moho depth gradient belts in the land and seas of China and adjacent areas

圖7 中國海陸莫霍面深度梯級帶及分區(qū)Fig.7 Moho depth gradient belts and sub－regions

（2）賀蘭山—龍門山梯級帶

該梯級帶長度約4000km，近SN向.莫霍面深度等值線梯度變化在龍門山段十分強烈，從西到東約260km寬度內莫霍面深度變化為58～44km，抬升14km.該梯級帶在中國境外向北一直延伸與貝加爾—色楞格梯級帶相接，但隨著向北延伸梯級帶的梯度呈變小趨勢，向南因為受到印度板塊與歐亞板塊碰撞構造效應的影響，而未能延伸很遠，形成東亞地區(qū)一條大型SN向構造帶.

（3）天山—阿爾金山—祁連山梯級帶

該梯級帶長度約4350km，近EW向，呈狹長狀地帶，該梯級帶構成我國西部南北不同地殼性質的界帶.莫霍面深度等值線梯度變化強烈，地殼厚度較周圍盆地區(qū)增厚約10km.阿爾金山附近，從北到南約130km寬度內，莫霍面深度變化為49～59km，地殼厚度增厚10km.該梯級帶是青藏高原周邊增厚亞區(qū)與北側正常地殼亞區(qū)的分界，將塔里木盆地與柴達木盆地分隔，形成沉積盆地區(qū)莫霍面上隆、造山帶區(qū)莫霍面深度下凹的特殊分布.

（4）哈薩克丘陵—阿爾泰山梯級帶

該梯級帶總長約3350km，莫霍面深度等值線梯度變化比較平緩，莫霍面深度變化為42～46km.

（5）喜馬拉雅梯級帶

該梯級帶總長約3760km，寬度較大，梯級帶走向由西段的NW向轉為東段的EW向.莫霍面深度等值線梯度變化強烈，從印度到青藏在垂直于梯級帶約500km距離內，莫霍面深度變化為38～66km，深度增加28km.該梯級帶是恒河平原—印度次大陸正常地殼與青藏高原增厚地殼的分界，是印度與歐亞兩大板塊碰撞、擠壓形成的過渡帶.

（6）印度半島東緣梯級帶

該梯級帶長度約3000km，主要沿著印度半島東緣海岸線，走向由斯里蘭卡東部的SN向—NE—孟加拉灣北部的近EW向，并且沿著該走向莫霍面深度等值線梯度變化由大逐漸變小，其中在斯里蘭卡東部從大陸到印度洋在垂直于梯級帶約140km距離內，莫霍面深度變化為35～25km，抬升10km.該梯級帶是印度半島正常地殼與印度洋大洋地殼的分界帶.

表8 中國海陸及鄰區(qū)莫霍面深度圖分區(qū)Table 8 Regional division of Moho depth in the land and seas of China and adjacent areas

（7）安達曼群島—尼亞斯島梯級帶

該梯級帶長度約3200km，呈近SN向，莫霍面深度等值線梯度變化強烈，在垂直于梯級帶約150km距離內，莫霍面深度變化為20～30km.

（8）中南半島西南梯級帶

該梯級帶長度約1890km，沿中南半島西南邊緣，從北到南，走向由SN變?yōu)镹NW向，莫霍面深度等值線梯度變化較平緩，呈弧狀分布，莫霍面深度變化為36～39km.

（9）中南半島東緣梯級帶

該梯級帶長度約1700km，走向沿中南半島東部邊緣呈近SN向，莫霍面深度等值線梯度變化較強，在垂直于梯級帶約200km距離內，莫霍面深度變化為33～40km.與中南半島西南梯級帶共同構成東南亞的正常大陸地殼分區(qū).

（10）南海過渡帶梯級帶

該梯級帶長度約4720km.莫霍面深度等值線梯度變化強烈，莫霍面深度變化為20～25km.該梯級帶是大陸型地殼與過渡型地殼的分界帶.

（11）日本?！獤|海梯級帶

該梯級帶長度約3150km，呈NNE向.莫霍面深度等值線梯度在日本海地區(qū)較大，在東海地區(qū)梯度相對較緩，莫霍面深度變化為25～30km.該梯級帶是東亞陸架地區(qū)的過渡帶，向南延伸至臺灣與南海過渡帶梯級帶相連.

（12）菲律賓海梯級帶

該梯級帶長度約3450km，梯級帶在臺灣島以北走向為NNE，以南為NNW向.莫霍面深度等值線梯度變化較強烈，莫霍面深度變化為15～24km，此梯級帶是菲律賓海板塊向歐亞大陸俯沖形成的莫霍面深度梯級帶，將菲律賓海與東海、南海隔開，是過渡型地殼與大洋型地殼的分界帶.

5.2 莫霍面深度分區(qū)

（Ⅰ）增厚型地殼區(qū)

增厚型地殼分區(qū)主體為青藏高原，以地殼明顯增厚為特征.包括羌塘盆地、措勤盆地和柴達木盆地，莫霍面深度變化為50～72km.根據深度變化可分為青藏高原地殼隆升亞區(qū)（Ⅰ1）和青藏高原周邊地殼增厚亞區(qū)（Ⅰ2）.前者為青藏高原主體，莫霍面深度為60～72km；后者包括了柴達木盆地及其他青藏高原周邊地區(qū)，莫霍面深度在50～60km之間分布.

（Ⅱ）正常型地殼區(qū)

正常型地殼區(qū)的東西邊界分別為大興安嶺—太行山—武陵山梯級帶和賀蘭山—龍門山—貝加爾—色楞格梯級帶.進一步可劃分為8個亞區(qū).

（1）哈薩克斯坦亞區(qū)（Ⅱ1）：主要由巴爾哈什湖和準噶爾盆地組成，區(qū)內莫霍面深度等值線近EW向，莫霍面深度變化范圍為38～45km.

（2）蒙古高原亞區(qū)（Ⅱ2）：主體為蒙古高原，莫霍面深度等值線變化平緩.莫霍面深度變化為45～50km，最深處位于埃恩赫塔伊萬山，為50km.

（3）中國西北亞區(qū)（Ⅱ3）：在蒙古高原亞區(qū)以南，主體為塔里木盆地，莫霍面深度等值線近EW向.莫霍面深度變化為45～50km，區(qū)內兩個莫霍面深度低值圈閉分別為塔里木盆地和吐魯番盆地，均為46km.

（4）貝加爾湖—蒙古東部亞區(qū)（Ⅱ4）：位于蒙古高原亞區(qū)以東，莫霍面深度等值線呈SN向，與大興安嶺—太行山—武陵山梯級帶方向基本一致，區(qū)內莫霍面從西到東逐漸抬升，深度變化為38～45km.

（5）鄂爾多斯盆地亞區(qū)（Ⅱ5）：北以陰山，南以秦嶺—大別山為界，由鄂爾多斯盆地和汾渭裂谷系構成，近NS向.莫霍面深度變化為38～46km.

（6）四川盆地及周邊亞區(qū)（Ⅱ6）：北以秦嶺—大別山，南以中南半島東緣梯級帶為界，由四川盆地和華南塊體西部構成，莫霍面深度變化為38～47km.從莫霍面深度等值線形態(tài)可以看出，青藏受印度板塊碰撞擠壓后向東溢出，遇到硬度較大的四川盆地基底受阻而變?yōu)闁|南向溢出形成云貴高原.

（7）中南半島北部亞區(qū)（Ⅱ7）：以賀蘭山—龍門山—貝加爾—色楞格梯級帶、中南半島西南梯級帶和中南半島東緣梯級帶為界，構成以中南半島為主體，主要包括東南亞一系列國家.其莫霍面深度變化為37～47km.

（8）恒河平原—印度半島正常地殼亞區(qū)（Ⅱ8）：位于喜馬拉雅造山帶以南，即恒河平原與印度半島地區(qū)，莫霍面深度變化為30～45km.

（Ⅲ）減薄型地殼區(qū)

減薄型地殼分區(qū)以地殼減薄為特征，東西以日本海—東?！獤|南沿海梯級帶和大興安嶺—太行山—武陵山梯級帶為界，進一步可劃分為5個亞區(qū).

（1）東北—鄂霍次克海亞區(qū)（Ⅲ1）：南以燕山為界，即東北塊體.莫霍面深度變化為32～38km，整體變化平緩.區(qū)內莫霍面深度等值線低值圈閉對應著松遼盆地，其西側等值線呈NNE向，北側近EW向，莫霍面深度為32km.

（2）華北—朝鮮半島亞區(qū)（Ⅲ2）：主體由華北塊體、朝鮮半島組成.莫霍面深度變化為30～38km.其中渤中地區(qū)莫霍面略有隆起.

（3）華南及沿海亞區(qū)（Ⅲ3）：由華南塊體和東海陸架和南海北部陸架構成.莫霍面深度等值線走向NNE，深度變化為25～38km.

（4）中南半島南部亞區(qū)（Ⅲ4）：區(qū)內包括安達曼海、馬來西亞、印尼部分地區(qū)，莫霍面深度為23～35km.西部安達曼海莫霍面深度等值線呈近SN向長軸橢圓狀隆升，莫霍面深度變化從約30km隆升至23km.

（5）加里曼丹及周邊亞區(qū)（Ⅲ5）：屬過渡型地殼區(qū)，由南海南部的加里曼丹島及周邊海域構成，莫霍面深度變化為25～28km.

（Ⅳ）溝弧盆體系地殼異常區(qū)

（1）東海溝弧盆亞區(qū)（Ⅳ1）：屬過渡型及大洋型地殼區(qū)，由日本海塊體與沖繩海槽組成，琉球海溝由菲律賓海板塊NW向俯沖歐亞大陸而形成，區(qū)內地震活動劇烈.莫霍面深度等值線呈NNE向，深度變化為16～30km.區(qū)內日本海、沖繩海槽有洋殼分布，其中沖繩海槽莫霍面深度最小約為16km.

（2）呂宋—棉蘭老島亞區(qū)（Ⅳ2）：由呂宋島、棉蘭老島等一系列群島組成，形成特有的溝－?。梵w系地殼異常區(qū)，夾于兩條俯沖帶之間，地震活動劇烈.莫霍面深度為20～26km.

（3）安達曼群島—尼亞斯島亞區(qū)（Ⅳ3）：由科科群島、安達曼群島、尼科巴群島、尼亞斯島等一系列群島組成，形成SN向，狹長狀的溝－?。梵w系地殼異常區(qū)，沿俯沖帶分布，地震活動劇烈.莫霍面深度為15～25km.

（Ⅴ）過渡型地殼區(qū)

南海周邊過渡型地殼亞區(qū)（Ⅴ1）：該區(qū)包含南海海盆、蘇祿海盆與蘇拉威西海盆周邊地殼陡變帶地區(qū)，該區(qū)莫霍面深度等值線走向均是沿海盆呈閉合狀，深度范圍在15～25km之間變化.

（Ⅵ）大洋型地殼區(qū)

該區(qū)主要位于南海、菲律賓海、印度洋等海域，莫霍面深度變化相對陡變帶地區(qū)較為平緩，屬于大洋型地殼，可進一步分為5個亞區(qū).除這5個較大區(qū)域的洋殼亞區(qū)外，在日本海、沖繩海槽、望加錫海峽、安達曼海也有洋殼分布.

（1）南海海盆洋殼亞區(qū)（Ⅵ1）：該地區(qū)莫霍面深度特征是，在中央海盆呈近SN向展布，西南次海盆則呈NE向，西北次海盆呈NEE向.海盆區(qū)受海底擴張方向的影響，莫霍面深度等值線走向也相應改變.區(qū)內莫霍面深度范圍在8～15km之間變化.總體上海盆區(qū)莫霍面深度淺，海山區(qū)莫霍面相對較深，如黃巖海山、憲南海山、珍貝海山、中南海山、龍北海山和龍南海山等形成的高值圈閉，這些海山區(qū)的莫霍面基本上比周邊海盆深約2km.

（2）蘇祿海洋殼亞區(qū)（Ⅵ2）：莫霍面深度等值線走向為NE向，深度范圍為10～15km，形態(tài)呈閉合狀.

（3）蘇拉威西海盆洋殼亞區(qū)（Ⅵ3）：莫霍面深度為8～15km，形態(tài)呈閉合狀，其西南側的望加錫海峽也有洋殼分布.深源地震密集分布于該區(qū)東部，最大震源深度達650km.

（4）菲律賓海洋殼亞區(qū)（Ⅵ4）：莫霍面深度等值線走向基本沿俯沖帶方向，深度范圍為7～20km.

（5）印度洋北部洋殼亞區(qū)（Ⅵ5）：該區(qū)莫霍面深度變化相對平緩，中央90°SN向隆起帶將該區(qū)分為兩部分.西部莫霍面深度變化較緩，在13～25km之間變化；而東部莫霍面深度變化相對較陡，近SN走向，莫霍面深度范圍11～17km.

6 幾點認識

本次編圖涉及多個海域，有些海域（如南海）海底地形變化大，起伏劇烈.傳統(tǒng)布格改正與海底地形呈負相關，海底地形的影響很大，為了得到精確的布格重力異常，并解決海陸拼接的問題，本次編圖采用對空間重力異常進行完全布格改正的方法，消除了海底地形的影響.格萊尼重力異常消除了全球地形質量的影響和166.7km以外補償質量的影響，可以認為它是166.7km以內的深部和淺部異常地質體的重力效應，為了更好地反演莫霍面起伏，本次編圖對編圖范圍內進行了格萊尼校正，得到了格萊尼重力異常，使反演結果更加準確，也使編圖工作更加規(guī)范.

從莫霍面深度圖可以得知中國海陸及鄰區(qū)莫霍面深度的總體特征呈現(xiàn)“東西分帶，南北分塊”的特點.在中國東部地區(qū)，莫霍面深度從西向東的變化幅度，遠遠大于從北向南（不包括海區(qū)）的變化幅度.在整體宏觀構造格架上，劉光鼎（2007）提出的“三橫兩豎”中國大地構造宏觀構造格架理論，不僅適用于中國大陸內部，也適用于整個東亞地區(qū).從本次編繪的莫霍面深度圖來看，“兩豎”（大興安嶺—太行山—武陵山梯級帶和賀蘭山—龍門山梯級帶）有清晰的表現(xiàn)，表明它們是深部構造的地球物理場響應.此外，中國海陸及鄰區(qū)莫霍面深度圖中太平洋第一島鏈梯級帶，同樣也是明顯的地殼厚度變化帶，揭示出了兩側地殼性質的改變.

從本文得到的莫霍面深度圖可以發(fā)現(xiàn)，大興安嶺－太行山－武陵山莫霍面深度梯級帶向北經呼倫貝爾延伸進入了蒙古國和俄羅斯，但梯級帶越向北越寬緩；而梯級帶向南進入中南塊體后，顯然被后期印度板塊向歐亞板塊俯沖所引發(fā)的地殼結構變形蹤跡所阻斷扭曲，梯度帶的方向發(fā)生了明顯的變化，應是多次構造運動疊加效應的反映.賀蘭山－龍門山莫霍面深度梯度帶向北進入蒙古國后方向就發(fā)生了改變，推斷是受到了興蒙塊體的影響；向南則基本未延伸出中國大陸境外，同樣與印度板塊向北的俯沖碰撞有密切的關系.以上分析表明，莫霍面深度圖中反映的信息，應是歐亞、印度和太平洋三大板塊相互作用下深部構造的反映，更多的是三大板塊作用加上板塊內部各塊體間相互作用的綜合表現(xiàn).

從莫霍面深度圖中可見，雅魯藏布江處的莫霍面深度要大于喜馬拉雅山處，但根據前人研究結果（李廷棟，2010），雅魯藏布江處的巖石圈厚度卻明顯薄于喜馬拉雅山處，形成較大的反差.由此可以推斷出，由南向北發(fā)生俯沖的板塊前緣為大洋型巖石圈性質，地質歷史上應屬于特提斯的一部分.這與劉光鼎（1992）所指出的“中生代時，中國大陸東西兩側受一系列SN向的轉換斷層和EW向延展的、連接特提斯?！《妊笈c太平洋的洋脊體系相互作用，形成目前表現(xiàn)為U字形的特提斯構造格局”的認識是一致的.由于SN向的持續(xù)擠壓作用，在中國西南部形成大規(guī)模特提斯－喜馬拉雅構造帶，也就是劉光鼎（1992）提出的三條西部鋒線.這三條西部鋒線的形成過程中，在SN向的擠壓應力作用下，伴有EW向的拉張作用，因而形成深達殼幔邊界的SN向“縱谷”.這些縱谷在莫霍面深度圖上表現(xiàn)得很清晰，進而也證實了青藏高原的地殼增厚是在板塊俯沖過程中，特提斯不斷關閉、拉開過程的產物.

從大陸地區(qū)向海區(qū)，莫霍面深度也反映了地殼性質的轉變.但同為洋殼的情況下，圖8中南海中央海盆處的莫霍面起伏形態(tài)及埋深情況與蘇祿海的情況仍有所不同，推其原因應是拉張機制與拉張速率存在著差異的反映.

本次編圖涉及到西太平洋俯沖帶和印度洋俯沖帶，存在著不同板塊在不同方向的俯沖、對沖等復雜的地質現(xiàn)象，在這樣的地區(qū)莫霍面不僅起伏劇烈，而且構造復雜.如圖幅內西太平洋俯沖帶東側俯沖方向由最北端（日本南部）的NW向南逐漸變?yōu)榉坡少e島的SWW向，至馬來西亞北部沙巴地區(qū)轉變?yōu)镹WW向，俯沖帶寬度和深度都急劇增加，俯沖深度可達700km；西側俯沖帶俯沖方向由中國臺灣島北部的E向轉為呂宋島南部的NEE向，至巴拉望地區(qū)又轉為E向，在馬來西亞北部沙巴地區(qū)表現(xiàn)為SEE向.對這個地區(qū)截取的震中分布剖面甚至可觀察到俯沖帶對沖現(xiàn)象1）郝天珧等，2013.中國海陸1∶500萬莫霍面深度圖驗收報告.中國科學院地質與地球物理研究所..西太平洋俯沖帶反映了菲律賓海板塊沿琉球海溝對歐亞板塊的俯沖，同時伴隨著15Ma以來南海洋殼向東的俯沖消亡.印度洋俯沖帶大致沿爪哇海溝向北延伸，俯沖方向大體上為E向，在11°N—13°N附近出現(xiàn)東側俯沖帶，俯沖方向向西，反映了印度—澳大利亞板塊沿印度洋東北緣的爪哇海溝向蘇門答臘、爪哇等島嶼之下的俯沖.因此，洋陸相互作用和洋－洋的相互作用是一個復雜的過程，而俯沖帶處莫霍面深度的計算與描述是未來地球物理研究面臨的重要挑戰(zhàn)之一，也是“俯沖工廠”研究的重要組成部分.

圖8 南海沿剖面II′的重力異常特征與莫霍面深度分布（據秦靜欣等，2011，有修改，剖面位置見圖4b）Fig.8 Characteristics of gravity anomalies and distribution of the Moho depths along the profile II′in the South China Sea（after Qin et al.，2011，Location of the profile is shown in Fig.4b）

致謝 本文得益于劉光鼎院士的學術思想和所提出的“一、二、三、多”的編圖方法論，在此表示衷心感謝！中國地質調查局青島海洋地質研究所張訓華研究員、溫珍河研究員、楊金玉研究員、韓波副研究員、郭興偉副研究員，王忠蕾助理研究員，中國地質調查局發(fā)展研究中心張明華教授、喬計花高級工程師，中國地質調查局航空物探遙感中心周堅鑫教授、尹航高級工程師以及中國地震局地殼應力研究所的黃忠賢研究員等都對本文給予了大量的幫助.此外，同濟大學王家林教授，中國地質調查局青島海洋地質研究所的郭振軒研究員、中國地質調查局廣州海洋地質調查局的陳邦彥教授、新疆石油學院李良辰教授等，在本文形成與編圖過程中提出大量的寶貴意見和提供了十分珍貴的數據資料，在此一并致謝！

Fang J.1999.Global crustal and lithospheric thickness inversed by using satellite gravity data.CrustalDefformationandEarthquake（in Chinese），19（1）：26－31.

Gao X，Wang W M，Yao Z X.2005.Crustal structure of China mainland and its adjacent regions.ChineseJ.Geophys.（in Chinese），48（3）：591－601.

Huang J P，F(xiàn)u R S，Xu P，et al.2006.Inversion of gravity and topography data for crust thickness of China and its adjacency.ActaSeismologicaSinica（in Chinese），28（3）：251－258.

Jia S X，Zhang X K.2005.Crustal structure and comparison of different tectonic blocks in North China.ChineseJ.Geophys.（in Chinese），48（3）：611－620.

Li S L，Mooney W D.1998.Crustal structure of China from deep seismic sounding profiles.Tectonophysics，288（1－4）：105－113.

Li S L，Mooney W D，F(xiàn)an J C.2006.Crustal structure of mainland China from deep seismic sounding data.Tectonophysics，420（1－2）：239－252.

Li T D.2010.Principal characteristics of the lithosphere of China.EarthScienceFrontiers（in Chinese），17（3）：1－13.

Liu F T，Xu P F，Liu J S，et al.2003.The crustal velocity structure of the continental deep subduction belt：study on the eastern Dabie orogen by seismic wide－angle reflection／refraction.ChineseJ.Geophys.（in Chinese），46（3）：366－372.

Liu G D.1992.The characteristic of geological and geophysical in China Seas and adjacent regions（in Chinese）.Beijing：Science Press.

Liu G D.1993.Map series of geology and geophysics of China Seas and adjacent regions（in Chinese）.Beijing：Science Press.

Liu G D.2007.Geodynamical evolution and tectonic framework of China.EarthScienceFrontiers（in Chinese），14（3）：39－46.

Lu Z X，Meng P Z，Zhao J M，et al.1989.Velocity structure characteristics of the crust and upper mantle in Luyang－Haicheng－Donggou DSS profile and the earthquake.Northeastern SeismologicalResearch（in Chinese），5（4）：1－11.

Mcintosh K D，Nakamura Y，Wang T K，et al.2005.Crustal－scale seismic profiles across Taiwan and western Philippine Sea.Tectonophysics，401（1－2）：23－54.

Qin J X，Hao T Y，Xu Y，et al.2011.The distribution characteristics and the relationship between the tectonic units of the Moho Depth in South China Sea and Adjacent Areas.ChineseJ.Geophys.（in Chinese），54（12）：3171－3183.

Qiu X L，Ye S Y，Wu S M，et al.2001.Crustal structure across the Xisha Trough，northwestern South China Sea.Tectonophysics，341（1－4）：179－193.

Qiu X L，Zhao M H，Ao W，et al.2011.OBS survey and crustal structure of the Southwest Sub－basin and Nansha Block，South China Sea.ChineseJ.Geophys.（in Chinese），54（12）：3117－3128.

Ruan A G，Niu X W，Qiu X L，et al.2011.A wide angle Ocean Bottom Seismometer profile across Liyue Bank，the southern margin of South China Sea.ChineseJ.Geophys.（in Chinese），54（12）：3139－3149.

Teng J W，Yan Y F，Zhang H，et al.2002.Depth distribution of Moho and tectonic framework in eastern Asian continent and its adjacent ocean areas.ScienceinChina（SeriesD：Earth Sciences）（in Chinese），32（2）：89－100.

Wang C Y，Lin Z Y，Chen X B.1995.Comprehensive study of geophysics on geoscience transect from Menyuan，Qinghai province，to Ningde，F(xiàn)ujian province，China.ChineseJ.Geophys.（ActaGeophysicaSinica）（in Chinese），38（5）：590－598.

Xiong S B，Jin D M，Sun K Z，et al.1991.Some characteristics of deep structure of the Zhangzhou geothermal field and it′s neighborhood in the Fujian province.ChineseJ.Geophys.（ActaGeophysicaSinica）（in Chinese），34（1）：55－63.

Yan P，Zhou D，Liu Z S.2001.A crustal structure profile across the northern continental margin of the South China Sea.Tectonophysics，338（1）：1－21.

Zeng R S，Sun W G，Mao T E，et al.1995.The map of Moho depth in China.ActaSeismologicaSinica（in Chinese），17（3）：322－327.

Zhang C K，Zhao J R，Ren Q F，et al.1994.Study on crust and upper mantle structure in north Henan and its surroundings.SeismologyandGeology（in Chinese），16（3）：243－253.

Zhang Z J，Bai Z M，Mooney W，et al.2009.Crustal structure across the Three Gorges area of the Yangtze platform，central China，from seismic refraction／wide－angle reflection data.Tectonophysics，475（3－4）：423－437.

Zhao J R，Zhang X K，Zhang C K，et al.1999.The crust mantle tectonics and velocity structure characteristics in Xianghe－Beijing－Zhulu and its adjacent areas.SeismologyandGeology（in Chinese），21（1）：29－36.

Zhao M H，Qiu X L，Ye C M，et al.2004.Analysis on deep crustal structure along the onshore－offshore seismic profile across the Binhai（Littoral）Fault Zone in northeastern South China Sea.ChineseJ.Geophys.（in Chinese），47（5）：845－852.

Zhu J S.1996.The map of Moho depth（1∶12000000）.／／Yuan X C.Geophysical Series Maps of China（in Chinese）.Beijing：Geological Publishing House.

附中文參考文獻

方劍.1999.利用衛(wèi)星重力資料反演地殼及巖石圈厚度.地殼形變與地震，19（1）：26－31.

高星，王衛(wèi)民，姚振興.2005.中國及鄰近地區(qū)地殼結構.地球物理學報，48（3）：591－601.

黃建平，傅容珊，許萍等.2006.利用重力和地形觀測反演中國及鄰區(qū)地殼厚度.地震學報，28（3）：251－258.

嘉世旭，張先康.2005.華北不同構造塊體地殼結構及其對比研究.地球物理學報，48（3）：611－620.

李廷棟.2010.中國巖石圈的基本特征.地學前緣，17（3）：1－13.

劉福田，徐佩芬，劉勁松等.2003.大陸深俯沖帶的地殼速度結構——東大別造山帶深地震寬角反射／折射研究.地球物理學報，46（3）：366－372.

劉光鼎.1993.中國海區(qū)及鄰域地質地球物理圖集.北京：科學出版社.

劉光鼎.1992.中國海區(qū)及鄰域地質地球物理特征.北京：科學出版社.

劉光鼎.2007.中國大陸構造格架的動力學演化.地學前緣，14（3）：39－46.

盧造勛，孟樸在，趙俊猛等.1989.閭陽－海城－東溝剖面地殼與上地幔速度結構特征與地震.東北地震研究，5（4）：1－11.

秦靜欣，郝天珧，徐亞等.2011.南海及鄰區(qū)莫霍面深度分布特征及其與各構造單元的關系.地球物理學報，54（12）：3171－3183.

丘學林，趙明輝，敖威等.2011.南海西南次海盆與南沙地塊的OBS探測和地殼結構.地球物理學報，54（12）：3117－3128.

阮愛國，牛雄偉，丘學林等.2011.穿越南沙禮樂灘的海底地震儀廣角地震試驗.地球物理學報，54（12）：3139－3149.

滕吉文，閆雅芬，張慧等.2002.東亞大陸及周邊海域Moho界面深度分布和基本構造格局.中國科學（D輯：地球科學），32（2）：89－100.

王椿鏞，林中洋，陳學波.1995.青海門源—福建寧德地學斷面綜合地球物理研究.地球物理學報，38（5）：590－598.

熊紹柏，金東敏，孫克忠等.1991.福建漳州地熱田及其鄰近地區(qū)的地殼深部構造特征.地球物理學報，34（1）：55－63.

曾融生，孫為國，毛桐恩等.1995.中國大陸莫霍界面深度圖.地震學報，17（3）：322－327.

張成科，趙金仁，任青芳等.1994.豫北及其外圍地區(qū)地殼上地幔結果研究.地震地質，16（3）：243－253.

趙金仁，張先康，張成科等.1999.香河—北京—涿鹿及其相鄰地區(qū)殼幔構造與速度結構特征.地震地質，21（1）：29－36.

趙明輝，丘學林，葉春明等.2004.南海東北部海陸深地震聯(lián)測與濱海斷裂帶兩側地殼結構分析.地球物理學報，47（5）：845－852.

朱介壽.1996.1∶1200萬莫霍界面深度圖.／／袁學誠.中國地球物理圖集.北京：地質出版社.