青藏高原東緣及四川盆地的殼幔導(dǎo)電性結(jié)構(gòu)研究

張樂天，金 勝，3＊，魏文博，3，葉高峰，3，段書新，董 浩，張 帆，謝成良

1 中國(guó)地質(zhì)大學(xué)（北京）地球物理與信息技術(shù)學(xué)院，北京 100083

2 地下信息探測(cè)技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室，北京 100083

3 地質(zhì)過程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100083

4 中核集團(tuán)核工業(yè)北京地質(zhì)研究院，北京 100092

1 引 言

始于大約55Ma前［1］的印度板塊與亞歐板塊之間的陸－陸碰撞過程造就了當(dāng)今青藏高原雄偉壯觀的地表特征.這一自新生代以來持續(xù)進(jìn)行的造山運(yùn)動(dòng)在漫長(zhǎng)的地質(zhì)歷史時(shí)期中致使印度板塊與亞歐板塊之間的距離縮短量高達(dá)1360km左右［2］.Tapponnier等［3］通過分析力學(xué)模擬實(shí)驗(yàn)的結(jié)果認(rèn)為這一碰撞過程主要是通過高原內(nèi)部物質(zhì)的側(cè)向擠出來吸收的，即在印度板塊的擠壓作用下，青藏高原東緣的物質(zhì)沿著各個(gè)次級(jí)塊體之間的走滑斷裂不斷向青藏高原東南部及云南西北部地區(qū)逃逸.Clark與Royden［4］通過分析青藏高原及其周邊地區(qū)的地形高程分布情況認(rèn)為，在重力差異的作用下青藏高原內(nèi)部應(yīng)當(dāng)廣泛存在著殼內(nèi)物質(zhì)流，這些管流通道在高原內(nèi)部物質(zhì)運(yùn)移及地表隆升的過程中發(fā)揮了重要作用.Clark與Royden還進(jìn)一步推斷了這些物質(zhì)運(yùn)移通道的分布情況，認(rèn)為在青藏高原東緣，由于受到剛性的四川盆地的阻擋作用，物質(zhì)轉(zhuǎn)而向高原的東南緣地區(qū)逃逸.大地電磁測(cè)深方法利用具有豐富頻譜成分的天然電磁場(chǎng)作為場(chǎng)源，可以實(shí)現(xiàn)對(duì)地下物質(zhì)的電性結(jié)構(gòu)（電阻率或電導(dǎo)率）進(jìn)行巖石圈尺度的探測(cè)［5］.天然大地電磁場(chǎng)信號(hào)可以輕易穿透高阻區(qū)域，并對(duì)低阻體的存在十分敏感，可以有效地反映深部流體及高溫區(qū)域的分布情況，對(duì)巖石圈內(nèi)部的流變結(jié)構(gòu)具有重要的約束作用.近年來在青藏高原東緣地區(qū)開展的大地電磁測(cè)深工作在中下地殼范圍內(nèi)發(fā)現(xiàn)了分布較為廣泛的殼內(nèi)低阻層［6－8］，并被認(rèn)為是殼內(nèi)物質(zhì)流存在的電性反映，從而為上述論點(diǎn)提供了有力的證據(jù).綜合以上觀點(diǎn)可以看出，青藏高原東緣的構(gòu)造模式與高原內(nèi)部的物質(zhì)運(yùn)移以及高原的隆升機(jī)制之間存在著密切的聯(lián)系，高原東緣在吸收印度板塊與亞歐板塊之間的陸－陸碰撞過程中發(fā)揮著重要的作用.

位于青藏高原東緣的龍門山斷裂帶標(biāo)志著青藏高原與四川盆地之間的邊界.在這一區(qū)域，自西北向東南依次分布有松潘—甘孜地塊、四川盆地、川東滑脫褶皺帶等地質(zhì)構(gòu)造單元以及龍門山逆沖推覆斷裂帶、龍泉山斷裂、華鎣山斷裂等斷裂構(gòu)造體系.由于松潘—甘孜地塊的側(cè)向擠出運(yùn)動(dòng)受到剛性的四川盆地的阻擋作用，導(dǎo)致在龍門山斷裂帶區(qū)間形成了一系列呈疊瓦狀分布的逆沖推覆構(gòu)造體，并造成了非常明顯的地表高程差異.據(jù)估算，在龍門山斷裂帶區(qū)間東西寬約100km的范圍內(nèi)，地表高程落差可達(dá)6000m左右［9］.這種復(fù)雜多變的地形地貌分布預(yù)示著地下深部同樣存在著極其特殊的地質(zhì)構(gòu)造情況，在松潘—甘孜地塊與四川盆地之間不斷積累的應(yīng)力作用下，導(dǎo)致了該地區(qū)地震及各類地質(zhì)災(zāi)害頻發(fā)的現(xiàn)狀.發(fā)生于2008年的MS8.0級(jí)汶川大地震即為龍門山斷裂帶中的映秀—北川斷裂突發(fā)錯(cuò)動(dòng)的結(jié)果［10］.此外，四川盆地中部的龍泉山斷裂以及東部的華鎣山斷裂均為活動(dòng)斷裂，盆地內(nèi)部還有一些隱伏的盲斷層，上述這些斷裂均具有逆沖推覆或右旋斜沖推覆的性質(zhì)，存在著一定的地震危害性［11］.例如發(fā)生于2009年的四川瀘州MS4.0級(jí)地震就是由盆地東部華鎣山斷裂的活動(dòng)造成的.但是，前人的研究多集中于青藏高原東部及四川盆地的西部，主要是龍門山斷裂帶地區(qū)，對(duì)四川盆地東部的研究目前還較少.而且目前已有的研究對(duì)龍門山斷裂帶地區(qū)的構(gòu)造模式及地表的隆升機(jī)制等問題還存有爭(zhēng)議［12］，需要進(jìn)一步的深入研究.2011年，在國(guó)家深部探測(cè)技術(shù)與實(shí)驗(yàn)研究專項(xiàng)SinoProbe項(xiàng)目［13－14］的資助下，完成了一條跨越青藏高原東緣及整個(gè)四川盆地的大地電磁長(zhǎng)剖面［15］.該剖面自西北始于青藏高原內(nèi)部的松潘—甘孜地塊，向東南延伸穿過龍門山斷裂帶、四川盆地內(nèi)部及四川盆地東部的華鎣山斷裂，最終止于重慶東南的川東滑脫褶皺帶附近，旨在從電學(xué)性質(zhì)的角度研究青藏高原東緣及四川盆地區(qū)域地下深部的殼幔電性結(jié)構(gòu)，并由此進(jìn)一步認(rèn)識(shí)和探索這一區(qū)域內(nèi)的殼幔構(gòu)造模式及其大陸動(dòng)力學(xué)意義，為尚存有爭(zhēng)議的科學(xué)問題從電學(xué)性質(zhì)的角度提供約束與依據(jù)，以上也正是本文的研究意義所在.

2 大地電磁測(cè)深數(shù)據(jù)采集與處理

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

測(cè)區(qū)點(diǎn)位分布如圖1所示，大地電磁測(cè)深剖面整體呈北西—南東走向，基本沿與龍門山斷裂帶地表構(gòu)造走向垂直的方向布設(shè).剖面最西北端起自阿尼瑪卿縫合帶附近，在松潘—甘孜地塊內(nèi)部經(jīng)過若爾蓋與松潘等地，并以近垂直的方式穿過位于北川與汶川之間的龍門山斷裂帶中段；在四川盆地內(nèi)部，剖面于德陽附近穿過龍泉山斷裂北端，經(jīng)遂寧并穿過華鎣山斷裂，最后經(jīng)過重慶并止于川東滑脫褶皺帶以東南的地區(qū)，基本覆蓋了松潘—甘孜地塊東部、龍門山斷裂帶中段以及整個(gè)四川盆地的大部地區(qū).

野外的數(shù)據(jù)采集工作基本于2011年年內(nèi)完成，整條測(cè)線共包括68個(gè)寬頻大地電磁測(cè)深點(diǎn).采集儀器使用加拿大鳳凰公司生產(chǎn)的MTU－5寬頻大地電磁儀，平均點(diǎn)距約為10km左右，平均采集時(shí)間約為20h.原始時(shí)間序列經(jīng)過快速傅氏變換由時(shí)間域轉(zhuǎn)化為頻率域，并進(jìn)一步通過遠(yuǎn)參考［17］、Robust估計(jì)［18］等處理技術(shù)轉(zhuǎn)化為與頻率有關(guān)的阻抗張量信息.經(jīng)過處理后得到的視電阻率與相位曲線表明數(shù)據(jù)質(zhì)量普遍較好，可用頻點(diǎn)的平均周期范圍為0.003～2000s左右，部分?jǐn)?shù)據(jù)質(zhì)量較好的測(cè)點(diǎn)可用信號(hào)的最長(zhǎng)周期可達(dá)5000s以上.剖面所經(jīng)各個(gè)塊體內(nèi)部典型測(cè)點(diǎn)的視電阻率與相位曲線如圖2所示，其中100號(hào)測(cè)點(diǎn)位于松潘—甘孜地塊內(nèi)部若爾蓋附近，356號(hào)測(cè)點(diǎn)位于四川盆地內(nèi)部遂寧附近，432號(hào)測(cè)點(diǎn)位于重慶以東的川東滑脫褶皺帶地區(qū)，以上測(cè)點(diǎn)的具體位置已在圖1中標(biāo)出.由三個(gè)具代表性典型測(cè)點(diǎn)的視電阻率曲線可以看出，各個(gè)塊體內(nèi)部的電性結(jié)構(gòu)各不相同.松潘—甘孜地塊（100號(hào)測(cè)點(diǎn)）的電性結(jié)構(gòu)由淺部到深部表現(xiàn)為“低阻－高阻－低阻”的分布特征；四川盆地內(nèi)部（356號(hào)測(cè)點(diǎn)）則表現(xiàn)為淺部低阻、深部高阻的特點(diǎn)，并且TE（圖中xy）與TM（圖中yx）曲線基本重合，說明盆地內(nèi)部的電性結(jié)構(gòu)整體一維性較強(qiáng)；川東滑脫褶皺帶附近（432號(hào)測(cè)點(diǎn)）的電性結(jié)構(gòu)由淺部到深部則表現(xiàn)出由高阻到低阻的變化趨勢(shì).

2.2 維性分析

圖1 測(cè)區(qū)點(diǎn)位分布圖及構(gòu)造情況簡(jiǎn)圖紅色圓點(diǎn)為MT測(cè)點(diǎn)，構(gòu)造信息據(jù)Styron等［16］.Fig.1 Topography map showing MT stations layout and tectonic structures in the survey areaRed dots are MT stations，tectonic structures after Styron et al［16］.

圖2 各個(gè)地質(zhì)單元內(nèi)部典型測(cè)點(diǎn)的視電阻率與相位曲線Fig.2 Apparent resistivity and phase curves for typical stations inside of each geological unit

在獲得與頻率相關(guān)的阻抗張量信息后，需要對(duì)所得數(shù)據(jù)進(jìn)行維性分析，從而進(jìn)一步確認(rèn)適合的數(shù)據(jù)處理及反演方法（一維／二維／三維）.理想的一維情況下電性結(jié)構(gòu)僅在深度方向上有變化，而在水平方向上沒有變化，即一維層狀模型；理想的二維情況是指電性結(jié)構(gòu)在構(gòu)造走向方向上沒有變化，即假設(shè)電性結(jié)構(gòu)沿構(gòu)造走向方向是無限延伸的；三維情況則是指電性結(jié)構(gòu)在兩個(gè)正交的水平方向及深度方向上都有變化.可見，實(shí)際的地質(zhì)情況都是三維的，只有在滿足一些特殊條件的前提下，可以近似視為二維或一維情況處理.通常的二維大地電磁測(cè)深剖面都會(huì)盡量垂直于構(gòu)造走向方向布設(shè)就是為了盡可能的滿足二維性的要求.Swift［19］與 Bahr［20］各自提出了一種用以判別大地電磁測(cè)深數(shù)據(jù)是否滿足二維性要求的參數(shù)，即二維偏離度.其中，Swift二維偏離度較容易受各類畸變效應(yīng)的影響，從而影響判別的準(zhǔn)確性；相對(duì)來說Bahr二維偏離度受畸變效應(yīng)影響更小一些，判別結(jié)果也更為穩(wěn)定可靠.圖3中給出了這兩種二維偏離度沿剖面方向分布的擬斷面圖.一般來說，二維偏離度的值越小（藍(lán)色），大地電磁測(cè)深數(shù)據(jù)越接近二維情況，當(dāng)二維偏離度為0時(shí)則表明數(shù)據(jù)符合理想的二維條件；反之二維偏離度越大（紅色），則大地電磁測(cè)深數(shù)據(jù)的三維性更強(qiáng).通常認(rèn)為當(dāng)二維偏離度的值小于0.3時(shí)，可以近似視為二維情況.從圖3中可以看出，沿剖面大部區(qū)域基本滿足二維性假設(shè)（藍(lán)色區(qū)域），只有在剖面西北端的阿尼瑪卿縫合帶以及剖面中部的龍門山斷裂帶附近表現(xiàn)得具有三維性（紅色區(qū)域），這可能表明這些地區(qū)的構(gòu)造情況比較復(fù)雜.松潘—甘孜地塊、四川盆地以及川東滑脫褶皺帶各個(gè)塊體內(nèi)部的數(shù)據(jù)均表現(xiàn)出較好的二維性，特別是四川盆地內(nèi)部的很多測(cè)點(diǎn)接近于滿足理想的二維（或一維）情況.另外，相對(duì)來說，反映淺部的高頻段測(cè)深數(shù)據(jù)比反映深部的低頻段測(cè)深數(shù)據(jù)二維性要更好，說明深部可能存在三維效應(yīng)，或者是由于長(zhǎng)周期信號(hào)受干擾影響所致.總體說來，剖面整體表現(xiàn)基本滿足二維性假設(shè)，適于進(jìn)行二維反演.

2.3 構(gòu)造走向分析

在對(duì)大地電磁測(cè)深數(shù)據(jù)進(jìn)行二維反演之前，首先需要確定剖面所經(jīng)地區(qū)的區(qū)域構(gòu)造走向，并將MT數(shù)據(jù)旋轉(zhuǎn)至構(gòu)造走向方向上.此時(shí)才能將MT數(shù)據(jù)分解為兩組相互獨(dú)立的極化模式（TE／TM）并通過二維反演獲得比較準(zhǔn)確可靠的地下電性結(jié)構(gòu)信息.本文中利用相位張量［21－22］與阻抗張量 Groom－Bailey（GB）分解［23－24］兩種方法來綜合分析并確定區(qū)域構(gòu)造走向的方向.

相位張量被定義為阻抗張量實(shí)部矩陣的逆矩陣與虛部矩陣的乘積.其圖示可以用一橢圓來表示，橢圓長(zhǎng)軸對(duì)應(yīng)張量元素不變量的最大值，短軸對(duì)應(yīng)張量元素不變量的最小值，橢圓長(zhǎng)軸方向?qū)?yīng)電性主軸方向.一般情況下，橢圓長(zhǎng)軸與短軸越接近表明一維性越好，理想一維情況下，橢圓則退化為圓形.此外，每個(gè)相位張量橢圓還對(duì)應(yīng)一二維偏離角度值β，β的角度值越大表明MT數(shù)據(jù)的三維性越強(qiáng).相位張量工具在用于MT數(shù)據(jù)維性分析及區(qū)域構(gòu)造走向判斷上具有不易受三維畸變效應(yīng)影響的優(yōu)點(diǎn).圖4中給出了在0.1s、10s、100s、1000s四個(gè)周期值下相位張量的分布情況，其中每個(gè)橢圓對(duì)應(yīng)的二維偏離角β的值由不同顏色示出.由圖中可以看出，沿剖面絕大部分測(cè)點(diǎn)的二維偏離角β均小于5°（藍(lán)色），表明剖面數(shù)據(jù)二維性較好.特別是在淺部中、高頻段（0.1s與10s），位于四川盆地及松潘—甘孜地塊內(nèi)部的部分測(cè)點(diǎn)甚至比較接近一維情況（橢圓接近退化為圓形）.只有在龍門山斷裂帶附近出現(xiàn)二維偏離角β值較大的情況，說明在這一區(qū)域構(gòu)造情況比較復(fù)雜，可能存在三維效應(yīng).此外β的值在較低的頻段（100s及1000s）有逐漸變大的趨勢(shì)，這說明在研究區(qū)域深部可能存在三維結(jié)構(gòu)，或者是由于長(zhǎng)周期MT數(shù)據(jù)質(zhì)量相對(duì)來說較差造成的.以上由相位張量得出的維性分析結(jié)果與Swift及Bahr二維偏離度分析的結(jié)果也是基本一致的.由相位張量橢圓的長(zhǎng)軸方向分布可以看出，剖面上絕大部分測(cè)點(diǎn)數(shù)據(jù)或者表現(xiàn)出近似一維的特性，或者電性主軸方向大體平行于龍門山斷裂帶，而只有在較低的頻段（100s及1000s），四川盆地內(nèi)部的部分測(cè)點(diǎn)主軸方向逐漸轉(zhuǎn)變?yōu)榻媳毕?

在此基礎(chǔ)上，進(jìn)一步對(duì)MT數(shù)據(jù)進(jìn)行了多點(diǎn)、多頻段的阻抗張量GB分解［24］，并將各個(gè)頻段的統(tǒng)計(jì)結(jié)果以玫瑰圖的形式示于圖4中.可以看出，在高頻段（0.01～0.1s）電性主軸方向的分布比較均勻，玫瑰圖中沒有顯示出具有明顯優(yōu)勢(shì)的主軸方向，這應(yīng)當(dāng)是由于淺部數(shù)據(jù)的一維性較好造成的.在中、低頻段（1～10s及10～100s），玫瑰圖中指示出了比較明顯的電性主軸方向，結(jié)合地表地質(zhì)構(gòu)造的分布情況判斷區(qū)域構(gòu)造走向應(yīng)當(dāng)大致平行于龍門山斷裂帶的地表走向，即為北東35°左右.在低頻段（100～1000s），構(gòu)造走向方向逐漸轉(zhuǎn)變?yōu)榻媳钡姆较?以上結(jié)果與相位張量的分析結(jié)果也是一致的，從而進(jìn)一步證明了構(gòu)造走向分析結(jié)果的可靠性.據(jù)此，將剖面MT數(shù)據(jù)全部旋轉(zhuǎn)至區(qū)域構(gòu)造走向方向上（即北東35°）并進(jìn)行二維反演.

3 大地電磁測(cè)深數(shù)據(jù)反演

圖3 Swift與Bahr二維偏離度擬斷面圖Fig.3 Pseudosections of Swift and Bahr skewness along the profile

圖4 不同頻率的相位張量分布圖及對(duì)應(yīng)各頻段構(gòu)造走向分析結(jié)果玫瑰圖Fig.4 Phase tensor maps for 4different frequencies and rose diagrams showing strike analysis results for each corresponding frequency band

利用二維非線性共軛梯度算法（NLCG）［25］在不同反演參數(shù)條件下對(duì)不同模式的MT剖面數(shù)據(jù)（TE、TM以及TE＋TM）進(jìn)行了大量的反演.在反演過程中發(fā)現(xiàn)TE模式的數(shù)據(jù)受三維畸變效應(yīng)的影響比較嚴(yán)重，從而導(dǎo)致TE模式以及TE＋TM模式的二維反演結(jié)果RMS擬合差較大.TE模式的MT數(shù)據(jù)在反演過程中存在的問題主要表現(xiàn)在以下兩個(gè)方面：（1）長(zhǎng)周期的TE數(shù)據(jù)無法擬合；（2）部分測(cè)點(diǎn)TE模式的相位數(shù)據(jù)存在超象限的問題.蔡軍濤等［26］通過正演研究發(fā)現(xiàn)在對(duì)三維模型生成的MT數(shù)據(jù)進(jìn)行二維反演時(shí)，TE模式對(duì)模型的二維的近似程度要求遠(yuǎn)高于TM模式；當(dāng)三維結(jié)構(gòu)影響較明顯時(shí)，利用TM模式數(shù)據(jù)進(jìn)行二維反演比其它模式反演更合理，反演結(jié)果中的虛假異常也更少.同時(shí)，還建議對(duì)于實(shí)測(cè)數(shù)據(jù)的二維反演應(yīng)優(yōu)先考慮采用TM模式數(shù)據(jù)進(jìn)行.一般地說，TM模式數(shù)據(jù)能夠更好地反映與二維區(qū)域構(gòu)造走向方向近似平行的三維異常體［27］，并且更不容易受到由三維低阻異常體所產(chǎn)生的畸變效應(yīng)的影響［28］，因此單獨(dú)TM模式的反演應(yīng)當(dāng)能夠比較準(zhǔn)確地反映出地下介質(zhì)的電性結(jié)構(gòu)分布形態(tài).綜上，最終我們選取了擬合情況最好的TM模式反演結(jié)果，如圖5a所示.反演中使用的參數(shù)設(shè)置為：TM視電阻率誤差級(jí)數(shù)10%，相位誤差級(jí)數(shù)5%，正則化因子τ＝1.反演過程中使用了0.01～10000 s共六個(gè)數(shù)量級(jí)的MT數(shù)據(jù)，初始模型為100Ωm均勻半空間，經(jīng)過200次迭代，最終RMS反演擬合差為1.482，說明反演擬合情況較好.此外，在圖5中還給出了TM視電阻率與相位的實(shí)測(cè)數(shù)據(jù)（圖5b）以及模型響應(yīng)數(shù)據(jù)（圖5c）的擬斷面圖.通過對(duì)比可以看出，實(shí)測(cè)數(shù)據(jù)與由反演模型生成的響應(yīng)數(shù)據(jù)非常接近，更進(jìn)一步說明了反演的擬合情況比較理想以及反演模型的可靠性.

4 電性結(jié)構(gòu)模型分析

圖5 （a）TM模式二維反演電性結(jié)構(gòu)模型及TM視電阻率與相位的（b）實(shí)測(cè)數(shù)據(jù)及（c）響應(yīng)數(shù)據(jù)擬斷面圖Fig.5 （a）2DTM inversion model and pseudosections of（b）observed and（c）modeled TM data of apparent resistivity and phase

圖6 二維反演電性結(jié)構(gòu)模型構(gòu)造解釋圖莫霍面深度據(jù)文獻(xiàn)［29］［30］，并根據(jù)電性結(jié)構(gòu)模型略作修改，C為高導(dǎo)體，R為高阻體，CL為高導(dǎo)層，F(xiàn)為斷層.Fig.6 Interpretation of the electrical structure model obtained from 2DinversionMoho depth after references［29］［30］and modified based on the electrical structure.Abbreviations are：C-Conductor，R-Resistor，CL-Conductive Layer，F(xiàn)-Fault.

由圖6可以看出，二維剖面的電性結(jié)構(gòu)反演模型總體來說具有“橫向分塊，縱向分層”的特點(diǎn).剖面所經(jīng)各個(gè)構(gòu)造單元內(nèi)部的電性結(jié)構(gòu)各不相同，各主要斷裂帶在電性結(jié)構(gòu)上均表現(xiàn)為明顯的電性分界面或電性梯度帶.下文分別對(duì)各個(gè)構(gòu)造單元內(nèi)部的電性結(jié)構(gòu)進(jìn)行詳細(xì)描述和分析.

4.1 松潘—甘孜地塊

松潘—甘孜地塊位于北部的東昆侖—柴達(dá)木地塊與南部的羌塘地塊之間，其北部邊界為阿尼瑪卿縫合帶，南部則以金沙江縫合帶為界.松潘—甘孜地塊的西部在青藏高原內(nèi)部呈狹長(zhǎng)的條帶狀分布，而在青藏高原東緣則以龍門山斷裂帶為界，從而構(gòu)成了形如倒三角形的構(gòu)造單元（見圖1）.松潘—甘孜地塊最明顯的地表地質(zhì)特征為廣泛分布的巨厚的三疊紀(jì)復(fù)理石雜巖沉積［2］，一般厚度至少為幾公里，局部地區(qū)厚度可達(dá)到10km以上.并且在青藏高原東緣的逆沖作用下，這一雜巖地層產(chǎn)生了強(qiáng)烈的褶皺變形.從電性結(jié)構(gòu)來看，松潘甘孜地塊總體上表現(xiàn)為高阻的上地殼以及低阻的中下地殼.上地殼高阻層應(yīng)當(dāng)是三疊紀(jì)復(fù)理石雜巖的電性體現(xiàn)，高阻層底界面的起伏不平則反映了雜巖地層的褶皺變形作用.在中下地殼范圍內(nèi)分布有高導(dǎo)區(qū)域（圖6中CL1及C1），反映了青藏高原內(nèi)部地殼較弱的特點(diǎn).MT剖面在松潘—甘孜地塊內(nèi)部經(jīng)過了昆侖斷裂帶的最東端以及龍日壩斷裂帶的北端.高導(dǎo)層CL1在受到昆侖斷裂帶左行走滑的剪切作用下呈上凸的弧形分布，其高導(dǎo)電性的成因應(yīng)當(dāng)是在剪切作用下形成的破碎帶內(nèi)的流體造成的.Wang等人［31］通過在東昆侖斷裂帶開展的高精度反射地震剖面研究發(fā)現(xiàn)在這一地區(qū)昆侖斷裂帶表現(xiàn)為一殼內(nèi)尺度近垂直的左行走滑斷裂，其下方被一中地殼范圍內(nèi)近水平的剪切帶所截?cái)?，剪切帶下方的下地殼與巖石圈地幔耦合并被一系列的沖斷構(gòu)造所切割.從本文的電性結(jié)構(gòu)模型來看，昆侖斷裂帶附近的確沒有發(fā)現(xiàn)穿透整個(gè)巖石圈的電性結(jié)構(gòu)特征，因此很可能只是殼內(nèi)尺度的.而高導(dǎo)層CL1則很可能表征了中地殼內(nèi)受水平剪切作用所形成的滑脫面.此外，從高導(dǎo)層CL1的分布形態(tài)來看，其起伏變化較大且分布不是十分連續(xù)，因此可能表明這種中地殼范圍內(nèi)的剪切作用只是局部的，而這也與Wang等人［31］從地震剖面上所觀測(cè)到的結(jié)果是一致的.龍日壩斷裂帶是近年來通過GPS觀測(cè)發(fā)現(xiàn)并經(jīng)過衛(wèi)星影像及野外考察證實(shí)的一條北東向右行走滑斷裂帶，屬于晚第四紀(jì)新生的活動(dòng)斷裂帶［32］.在電性結(jié)構(gòu)模型中，龍日壩斷裂帶表現(xiàn)為一延伸深度較淺的電性分界面，同時(shí)還表現(xiàn)出具有一定的逆沖分量.在昆侖斷裂帶左行走滑以及龍日壩斷裂帶右行走滑的共同作用下，位于這兩條斷裂帶之間的物質(zhì)應(yīng)當(dāng)表現(xiàn)出整體向北東方向運(yùn)移的趨勢(shì).而位于這一區(qū)間中下地殼范圍內(nèi)的大規(guī)模高導(dǎo)體C1則很可能為這種運(yùn)動(dòng)提供了深部的動(dòng)力學(xué)條件.C1所在區(qū)域的電阻率值普遍不超過10Ωm，最小的電阻率值甚至不到1Ωm.這一低阻特征很可能是由局部熔融物質(zhì)或含鹽水流體形成的［7］，因此就有可能在C1所在區(qū)域內(nèi)形成殼內(nèi)物質(zhì)流動(dòng)的通道，從而在地下物質(zhì)發(fā)生大規(guī)模走滑運(yùn)動(dòng)的過程中起到引導(dǎo)的作用.值得一提的是，這一構(gòu)造特征與Bai等人在青藏高原東南緣鮮水河—小江斷裂帶以及班公—怒江縫合帶地區(qū)發(fā)現(xiàn)的殼內(nèi)物質(zhì)流動(dòng)通道［6］具有很高的相似性.在Bai等人的研究結(jié)果中所發(fā)現(xiàn)的兩組物質(zhì)流動(dòng)通道同樣表現(xiàn)為中下地殼范圍內(nèi)的大規(guī)模高導(dǎo)體，與本文模型中的高導(dǎo)體C1相比，無論在空間的分布形態(tài)與規(guī)模上還是在電阻率的數(shù)值上都具有很高的相似性.這些中下地殼的大規(guī)模高導(dǎo)體普遍分布于25km深度以下的中下地殼范圍內(nèi)，并可延伸至100km附近的上地幔區(qū)域；在水平方向上這些高導(dǎo)體的分布范圍也十分廣泛，由幾十至上百公里不等.這些高導(dǎo)區(qū)域的的電阻率值分布也非常相似，其核心區(qū)域的電阻率值普遍小于10Ωm.此外，這些中下地殼的大規(guī)模高導(dǎo)體大多都分布于一些大型走滑斷裂的附近.綜上，我們認(rèn)為這些高導(dǎo)體都應(yīng)當(dāng)具有類似的成因，并且反映了青藏高原東緣在殼內(nèi)物質(zhì)流動(dòng)與大型走滑斷裂的綜合作用下形成的復(fù)雜多變的構(gòu)造變形樣式.

4.2 龍門山斷裂帶

龍門山逆沖推覆斷裂帶標(biāo)志著青藏高原東部與四川盆地之間的邊界.該區(qū)地表地質(zhì)特征主要是廣泛分布的前寒武變質(zhì)巖帶.自東北向西南依次分布有彭灌雜巖、寶興雜巖及康定雜巖三個(gè)主要雜巖帶［12］.該斷裂帶主要由三條主干斷裂組成，自西向東分別為：汶川—茂縣斷裂、映秀—北川斷裂以及灌縣—江油斷裂［10］.這些斷裂均具逆沖推覆性質(zhì)，呈疊瓦狀分布于龍門山區(qū)間，并上覆于揚(yáng)子塊體西緣之上［33］.位于龍門山斷裂帶下方中下地殼范圍內(nèi)的高阻體R1很好地體現(xiàn)了揚(yáng)子地塊所應(yīng)具有的較為穩(wěn)定的高阻特征.同時(shí)，R1與C1之間形成了非常強(qiáng)烈的電性對(duì)比.R1體現(xiàn)了揚(yáng)子地塊所具有的較“冷”、較“硬”的特征，由于殼內(nèi)流體含量少，因此表現(xiàn)為高電阻率；而C1則代表了較“熱”、較“軟”的青藏高原內(nèi)部物質(zhì)，并且由于局部熔融體與含鹽水流體的存在使之形成了大規(guī)模的高導(dǎo)體.位于R1與C1之間的電性梯度帶則標(biāo)志著這兩類電性結(jié)構(gòu)之間的過渡與轉(zhuǎn)換.從R1的分布形態(tài)來看，揚(yáng)子塊體的西緣應(yīng)當(dāng)楔入了高原內(nèi)部的中下地殼，從而導(dǎo)致了青藏高原東緣地殼的增厚.另外，從電性結(jié)構(gòu)模型來看，龍門山斷裂帶地區(qū)的三條主干斷裂應(yīng)當(dāng)以較緩的傾角在20km左右的深度匯聚從而形成一向西北方向延伸的高導(dǎo)層CL2.CL2表征著位于推覆體下方的拆離面，其高導(dǎo)電性的成因應(yīng)當(dāng)是由破碎帶內(nèi)的流體造成的.CL2所代表的拆離面與該區(qū)內(nèi)的地震活動(dòng)具有密切的聯(lián)系.Wang等人［9］通過分析高精度測(cè)地?cái)?shù)據(jù)發(fā)現(xiàn)汶川地震所引起的滑動(dòng)位移主要沿一殼內(nèi)滑脫面發(fā)生，震源位置也位于這一滑脫面上.該面深度約為20km，呈近水平方向展布，在龍門山下方水平方向上延伸60km左右.其分布形態(tài)與CL2具有很高的一致性，據(jù)此推斷CL2應(yīng)當(dāng)是這一滑脫面的電性指示.王緒本等人［34］也曾在青藏高原東緣龍門山斷裂帶地區(qū)開展過大地電磁測(cè)深研究，所得出的龍門山斷裂帶電性結(jié)構(gòu)模型與本文的模型具有較好的一致性.龍門山斷裂帶在電性結(jié)構(gòu)上均表現(xiàn)為北西傾向連續(xù)展布的殼內(nèi)低阻層，其下方則是高阻的結(jié)晶基底，反映出了逆沖推覆構(gòu)造的特征.Bai等人［6］研究成果中的P3剖面盡管是從龍門山斷裂帶的南端穿過，但其電性結(jié)構(gòu)模型也依然反映出龍門山斷裂帶區(qū)間的電性結(jié)構(gòu)由淺部的殼內(nèi)低阻層以及其下的高阻基底構(gòu)成，具有與本文模型相似的逆沖推覆特征.綜上，在龍門山地區(qū)以CL2為界，其上部為以逆沖推覆為主要特征的薄皮構(gòu)造，而在深部則表現(xiàn)為楔入構(gòu)造.

4.3 四川盆地

四川中－新生代前陸盆地屬揚(yáng)子塊體的一部分，為較穩(wěn)定的地塊.盆地內(nèi)部地形平緩，平均海拔在300～500m左右，廣泛分布的沉積巖層呈近水平分布.以盆地內(nèi)部的龍泉山斷裂為界，可將四川盆地進(jìn)一步分為西部的成都平原以及位于龍泉山斷裂及華鎣山斷裂之間的盆中丘陵地區(qū).四川盆地在電性結(jié)構(gòu)模型上從淺部到深部表現(xiàn)出“低阻－高阻－低阻”的層狀分布的特點(diǎn).淺部的沉積層表現(xiàn)為低阻層，并且具有西厚東薄的特點(diǎn)，這應(yīng)當(dāng)與青藏高原東緣的逆沖推覆作用有關(guān).淺部低阻層下的高阻基底則具有西薄東厚的特點(diǎn)，且電阻率值向東不斷增大，從而在四川盆地東部的深部形成一大規(guī)模的高阻區(qū)域R2.龍泉山斷裂在電性結(jié)構(gòu)上表現(xiàn)為一北西傾向的低阻體，推測(cè)是逆沖推覆構(gòu)造的電性反映［35］.其規(guī)模與龍門山逆沖推覆體相比，在延伸深度及水平延伸距離上都要小很多.需要注意的是，由于本文中的MT剖面僅從龍泉山斷裂的北端穿過（如圖1），因此并未能反映出作為其主體的斷裂中段的電性特征.另外，從電性結(jié)構(gòu)模型推斷在龍門山斷裂帶東緣的山前地區(qū)深部應(yīng)存在一隱伏斷裂（圖6中F1）將龍門山逆沖推覆體下方的高阻基底R(shí)1與四川盆地的高阻基底分離開來.

4.4 華鎣山斷裂及川東滑脫褶皺帶

從地貌特征來看，位于四川盆地東緣的華鎣山是一明顯的地表地形分界線，標(biāo)志著由寬緩平坦的四川盆地向川東山嶺地帶的轉(zhuǎn)換.華鎣山以東的川東滑脫褶皺帶由多條近NE—SW走向的條狀背斜山地與向斜寬谷組成.在電性結(jié)構(gòu)模型上，分布于這一地區(qū)淺部的起伏不平的低阻層應(yīng)當(dāng)是這種滑脫褶皺構(gòu)造的電性體現(xiàn).在華鎣山斷裂下方存在一南東傾向的低阻體C2，延伸深度至少為30km左右，可能表明華鎣山斷裂具有向北西方向推覆的性質(zhì).同時(shí)，在華鎣山斷裂以東的滑脫褶皺帶內(nèi)還存在一組類似的具推覆特點(diǎn)的電性結(jié)構(gòu).據(jù)此推斷川東滑脫褶皺帶的淺部構(gòu)造可能具有整體向四川盆地內(nèi)部推覆的特征，而這一北西向逆沖推覆構(gòu)造的前緣應(yīng)當(dāng)就是華鎣山斷裂.川東滑脫褶皺帶的深部電性結(jié)構(gòu)整體表現(xiàn)非常高阻（R3），符合古老、穩(wěn)定的揚(yáng)子克拉通基底的電性特征.與淺部具滑脫褶皺特征的低阻層結(jié)合，很好地體現(xiàn)了這一地區(qū)薄皮構(gòu)造的特點(diǎn).此外，表征華鎣山斷裂的低阻體C2很可能在深部切割了高阻的揚(yáng)子基底，并構(gòu)成了四川盆地的東部邊界.

5 構(gòu)造涵義討論

以上詳細(xì)分析了剖面所經(jīng)各個(gè)主要構(gòu)造單元的電性結(jié)構(gòu).總體來看，青藏高原東緣及四川盆地地區(qū)的構(gòu)造模式可以分為兩個(gè)層次：在淺部中－上地殼主要是以逆沖推覆為特征的薄皮構(gòu)造；而在深部中－下地殼至上地幔范圍內(nèi)則表現(xiàn)為剛性的揚(yáng)子塊體嵌入柔性青藏塊體內(nèi)部的楔入式構(gòu)造.萬天豐［36］通過總結(jié)認(rèn)為這種鋸齒狀的嵌入構(gòu)造在陸－陸碰撞造山帶地區(qū)是一種普遍存在的構(gòu)造模式.而這種構(gòu)造模式的成因應(yīng)當(dāng)與由淺部到深部的物性變化有關(guān).從本文的電性結(jié)構(gòu)模型來看，淺部的電性變化劇烈，電性分塊表現(xiàn)較為破碎，反映了脆性物質(zhì)的逆沖推覆過程.在推覆體下方，如青藏高原東緣及川東滑脫褶皺帶地區(qū)的中－上地殼范圍內(nèi)存在著一組連續(xù)性較好的水平電性分界面（如CL2），標(biāo)志著推覆體下方的滑脫界面.在這組界面以下，電性結(jié)構(gòu)則轉(zhuǎn)換為大規(guī)模電性均勻分布的塊體（如C1、R1、R2及R3），標(biāo)志著深部以韌性變形為主的構(gòu)造特征.以滑脫界面為界，淺部的脆性物質(zhì)與深部的韌性物質(zhì)發(fā)生解耦，從而形成了這種構(gòu)造模式.該區(qū)域的殼幔S波速度結(jié)構(gòu)研究［37］也表明，在青藏高原東向擠壓與四川盆地強(qiáng)硬地殼阻擋的長(zhǎng)期作用下導(dǎo)致松潘—甘孜地塊內(nèi)部蓄積很大的應(yīng)變能量以及上、下地殼在殼內(nèi)低速層頂部邊界的解耦，從而在龍門山斷裂帶附近形成上地殼的鏟形逆沖推覆結(jié)構(gòu).龍門山斷裂帶地區(qū)的重力異常分布特征［38］也體現(xiàn)出了同樣的趨勢(shì)，由重力異常推測(cè)該區(qū)域密度分界面分為兩段，深部較陡的分界面為巖石圈塊體的邊界，淺部較緩的分界面則認(rèn)為與逆沖推覆構(gòu)造有關(guān).淺部較緩的分界面由地表大地構(gòu)造單元的界線往深部密度界面緩傾并與深部巖石圈塊體的邊界相連接，反映了中上地殼脆性變形層的地質(zhì)體由巖石圈塊體界線沿緩傾的密度界面推覆至地表大地構(gòu)造單元界線處的逆沖推覆過程.由此可見，波速結(jié)構(gòu)與重力異常場(chǎng)的分布特征從另一個(gè)側(cè)面印證了以上電性結(jié)構(gòu)模型解釋的合理性，降低了地球物理的多解性.

王緒本等人［34］通過大地電磁測(cè)深獲得的電性結(jié)構(gòu)模型推斷可能存在一個(gè)西延至若爾蓋地塊的泛揚(yáng)子陸塊.但是，從本文的電性結(jié)構(gòu)模型來看，揚(yáng)子塊體并沒有俯沖至松潘—甘孜塊體下方，而是呈楔入體插入青藏高原東部的中地殼下方.許志琴等人［39］根據(jù)地質(zhì)資料與地震剖面等地球物理資料的綜合研究也得出了相似的結(jié)論，并指出是否存在揚(yáng)子地塊往西運(yùn)動(dòng)及大陸俯沖作用尚待探究.地震P波三維速度結(jié)構(gòu)研究［40］表明，在龍門山斷裂帶30km深度之下的揚(yáng)子地塊具有明顯的高速特征，其前緣隨深度增加向青藏高原方向擴(kuò)展，在下地殼和上地幔頂部已達(dá)到龍門山斷裂帶以西.這一特征與本文電性結(jié)構(gòu)模型中的高阻體R1具有很好的一致性.此外，S波速度結(jié)構(gòu)的研究［37］也同樣表明并不存在四川盆地向西側(cè)的俯沖.Zhao等人［41］根據(jù)地震S波接收函數(shù)的研究提出亞洲巖石圈俯沖于青藏巖石圈之下分三個(gè)階段進(jìn)行的構(gòu)造模型.這三個(gè)階段分別為：首先是青藏巖石圈的增厚，第二階段是青藏巖石圈下部的剝離，最后是整個(gè)青藏巖石圈覆蓋于亞洲巖石圈之上.從本文得出的青藏高原東緣的構(gòu)造模式來看，應(yīng)當(dāng)處于第二階段的早期.我們認(rèn)為，在揚(yáng)子塊體向青藏高原中下地殼的楔入作用下，使得高原殼內(nèi)物質(zhì)流的東向運(yùn)移受到了阻擋，這些物質(zhì)在東緣的堆積導(dǎo)致了高原東部的地殼增厚.松潘—甘孜塊體內(nèi)部中下地殼范圍內(nèi)的大規(guī)模高導(dǎo)體C1應(yīng)當(dāng)就是在這種機(jī)制下形成的.龍門山斷裂帶中上地殼的P波高速異常表明介質(zhì)具有相對(duì)較高的強(qiáng)度［40］，在青藏高原物質(zhì)向東擠出過程中起到了較強(qiáng)的阻擋作用，有利于深部能量積累.S波速度結(jié)構(gòu)研究［37］也同樣表明在龍門山斷裂帶附近地殼平均泊松比僅為0.2，顯示了堅(jiān)硬地殼的特征，而松潘—甘孜地塊的地殼則相對(duì)軟弱.這些波速結(jié)構(gòu)研究的結(jié)果都進(jìn)一步印證了以上關(guān)于高導(dǎo)體C1形成機(jī)制的認(rèn)識(shí).同時(shí)，高導(dǎo)體C1的存在對(duì)高原東緣的物質(zhì)運(yùn)移、地表隆升等過程也都具有重要的意義.

本文的電性結(jié)構(gòu)模型對(duì)于青藏高原東緣造山帶的隆升機(jī)制也具有一定的啟示意義.長(zhǎng)期以來，針對(duì)青藏高原東緣的隆升機(jī)制主要存在以下兩類觀點(diǎn)：一種認(rèn)為隆升過程是通過逆沖推覆等一系列的殼內(nèi)脆性變形過程來實(shí)現(xiàn)的，龍門山斷裂帶區(qū)間頻發(fā)的地震活動(dòng)是這種觀點(diǎn)的最好例證［9，42］；另一種觀點(diǎn)則認(rèn)為隆升過程是由于高原內(nèi)部中－下地殼物質(zhì)流導(dǎo)致的地殼增厚引起的，殼內(nèi)物質(zhì)流的發(fā)現(xiàn)則是這種觀點(diǎn)的有力證據(jù)［6－8］.我們認(rèn)為，實(shí)際上這兩種機(jī)制可能在青藏高原東緣的隆升過程中同時(shí)發(fā)揮著作用，但是每種機(jī)制占主導(dǎo)地位的空間分布范圍又各不相同.電性結(jié)構(gòu)模型反映出的高阻脆性物質(zhì)特征主要表現(xiàn)在淺部中－上地殼以及位于最東邊緣的龍門山斷裂帶區(qū)間.在這些地區(qū)的隆升過程中，由逆沖推覆等殼內(nèi)脆性變形導(dǎo)致的地殼縮短起到了主要的作用.而殼內(nèi)物質(zhì)流導(dǎo)致的地殼增厚區(qū)域（圖6中C1）則集中在松潘—甘孜地塊內(nèi)部的中－下地殼范圍內(nèi)，并對(duì)高原內(nèi)部的隆升起到了主要的推動(dòng)作用.在這兩種機(jī)制的共同作用下，造成了青藏高原東部地區(qū)的整體隆升.當(dāng)然，需要指出的是，本文的大地電磁測(cè)深剖面僅僅從松潘—甘孜地塊的最東北端穿過，因此電性結(jié)構(gòu)模型也只能反映出這一局部地區(qū)的殼幔結(jié)構(gòu)特征，要對(duì)青藏高原東緣地區(qū)整體的殼幔結(jié)構(gòu)組成、物質(zhì)運(yùn)移形式以及地表隆升過程等問題加以全面的理解，還需要更進(jìn)一步全面和深入的研究.

6 結(jié) 論

6.1 在SinoProbe項(xiàng)目的資助下，完成了一條跨越青藏高原東緣及整個(gè)四川盆地的大地電磁測(cè)深剖面.經(jīng)過一系列常規(guī)處理流程，獲得了質(zhì)量較好的數(shù)據(jù)體.維性分析結(jié)果表明剖面數(shù)據(jù)整體二維性較好.在分析得出區(qū)域構(gòu)造走向后，通過二維反演方法得到了整條剖面的電性結(jié)構(gòu)模型.

6.2 電性結(jié)構(gòu)模型表明剖面整體具有“橫向分塊，縱向分層”的特點(diǎn).沿剖面可分為三個(gè)主要的電性結(jié)構(gòu)單元，分別為：淺部高阻、中下地殼存在大規(guī)模低阻體（C1）的松潘—甘孜地塊，淺部低阻、中下地殼相對(duì)高阻的四川盆地，以及華鎣山以東整體為高阻特征的揚(yáng)子克拉通地塊（R3）.四川盆地西部的龍門山區(qū)間存在著一組表征南東向逆沖推覆的電性結(jié)構(gòu).其下方的低阻層（CL2）說明這一系列推覆體在深部20km左右匯聚并形成一滑脫面.在滑脫面CL2的下方是高阻的揚(yáng)子塊體與低阻的青藏塊體在深部的過渡帶.四川盆地東部的華鎣山斷裂及川東滑脫褶皺帶附近同樣存在一組北西向逆沖推覆的電性結(jié)構(gòu)，其下方的大規(guī)模高阻體R3則表征了古老、穩(wěn)定的揚(yáng)子克拉通地塊.

6.3 青藏高原東緣及四川盆地地區(qū)的整體構(gòu)造模式可以分為兩個(gè)層次：在淺部中－上地殼主要是以逆沖推覆為特征的薄皮構(gòu)造；而在深部中－下地殼至上地幔范圍內(nèi)則表現(xiàn)為剛性的揚(yáng)子塊體嵌入柔性青藏塊體內(nèi)部的楔入式構(gòu)造.這一構(gòu)造模式與研究區(qū)域內(nèi)殼內(nèi)物質(zhì)流的分布、青藏塊體與揚(yáng)子塊體在深部的關(guān)系以及青藏高原東緣地區(qū)的隆升機(jī)制等問題都具有緊密的內(nèi)在聯(lián)系.

（References）

［1］Rowley D B.Age of initiation of collision between India and Asia：A review of stratigraphic data.Earth and Planetary Science Letters，1996，145（1）：1－13.

［2］Yin A，Harrison T M.Geologic Evolution of the Himalayan－Tibetan Orogen.Annual Review of Earth and Planetary Sciences，2000，28：211－280.

［3］Tapponnier P，Peltzer G，LeDain A Y，et al.Propagating extrusion tectonics in Asia： New insights from simple experiments with plasticine.Geology，1982，10（12）：611－616.

［4］Clark M K，Royden L H.Topographic ooze：Building the eastern margin of Tibet by lower crustal flow.Geology，2000，28（8）：703－706.

［5］Jones A G.Imaging the continental upper mantle using electromagnetic methods.Lithos，1999，48（1－4）：57－80.

［6］Bai D H，Unsworth M J，Meju M A，et al.Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging.Nature Geoscience，2010，3（5）：358－362.

［7］趙國(guó)澤，陳小斌，王立鳳等.青藏高原東邊緣地殼“管流”層的電磁探測(cè)證據(jù).科學(xué)通報(bào)，2008，53（3）：345－350.Zhao G Z，Chen X B，Wang L F，et al.Evidence of crustal＇channel flow＇in the eastern margin of Tibetan Plateau from MT measurements.Chinese Science Bulletin，2008，53（12）：1887－1893.

［8］萬戰(zhàn)生，趙國(guó)澤，湯吉等.青藏高原東邊緣冕寧—宜賓剖面電性結(jié)構(gòu)及其構(gòu)造意義.地球物理學(xué)報(bào)，2010，53（3）：585－594.Wan Z S，Zhao G Z，Tang J，et al.The electrical structure of the crust along Mianning－Yibin profile in the eastern edge of Tibetan plateau and its tectonic implications.Chinese J.Geophys.（in Chinese），2010，53（3）：585－594.

［9］Wang Q，Qiao X J，Lan Q G，et al.Rupture of deep faults in the 2008Wenchuan earthquake and uplift of the Longmen Shan.Nature Geoscience，2011，4（9）：634－640.

［10］張培震，徐錫偉，聞學(xué)澤等.2008年汶川8.0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因.地球物理學(xué)報(bào)，2008，51（4）：1066－1073.Zhang P Z，Xu X W，Wen X Z，et al.Slip rates and recurrence intervals of the Longmen Shan active fault zone，and tectonic implications for the mechanism of the May 12 Wenchuan earthquake，2008，Sichuan，China.Chinese J.Geophys.（in Chinese），2008，51（4）：1066－1073.

［11］嵇少丞，王茜，孫圣思等.亞洲大陸逃逸構(gòu)造與現(xiàn)今中國(guó)地震活動(dòng).地質(zhì)學(xué)報(bào)，2008，82（12）：1644－1667.Ji S C，Wang Q，Sun S S，et al.Continental extrusion and seismicity in China.Acta Geologica Sinica （in Chinese），2008，82（12）：1644－1667.

［12］Xu Z Q，Ji S C，Li H B，et al.Uplift of the Longmen Shan range and the Wenchuan earthquake.Episodes，2008，31（3）：291－301.

［13］Dong S W，Li T D，Gao R，et al.A multidisciplinary Earth science research program in China.Eos Trans.AGU，2011，92（38）：313－314.

［14］董樹文，李廷棟.SinoProbe——中國(guó)深部探測(cè)實(shí)驗(yàn).地質(zhì)學(xué)報(bào)，2009，83（7）：895－909.Dong S W，Li T D.SinoProbe：the exploration of the deep interior beneath the Chinese continent.Acta Geologica Sinica（in Chinese），2009，83（7）：895－909.

［15］金勝，張樂天，魏文博等.中國(guó)大陸深探測(cè)的大地電磁測(cè)深研究.地質(zhì)學(xué)報(bào)，2010，84（6）：808－817.Jin S，Zhang L T，Wei W B，et al.Magnetotelluric method for deep detection of Chinese continent.Acta Geologica Sinica （in Chinese），2010，84（6）：808－817.

［16］Styron R，Taylor M，Okoronkwo K.Database of active structures from the Indo－Asian collision.Eos Trans.AGU，2010，91（20）：181－182.

［17］Gamble T，Goubau W，Clarke J.Magnetotellurics with a remote magnetic reference.Geophysics，1979，44（1）：53－68.

［18］Egbert G D，Booker J R.Robust estimation of geomagnetic transfer functions.Geophysical Journal of the Royal Astronomical Society，1986，87（1）：173－194.

［19］Swift C M.A magnetotelluric investigation of an electrical conductivity anomaly in the southwestern United States［Ph.D.thesis］.Cambridge：Massachusetts Institute of Technology，1967.

［20］Bahr K. Geological noise in magnetotelluric data： a classification of distortion types.Physics of the Earth and Planetary Interiors，1991，66（1－2）：24－38.

［21］Caldwell T G，Bibby H M，Brown C.The magnetotelluric phase tensor.Geophysical Journal International，2004，158（2）：457－469.

［22］Moorkamp M.Comment on ‘The magnetotelluric phase tensor’by T.Grant Caldwell，Hugh M.Bibby and Colin Brown.Geophysical Journal International，2007，171（2）：565－566.

［23］Groom R W，Bailey R C.Decomposition of magnetotelluric impedance tensors in the presence of Local three－Dimensional galvanic distortion.Journal of Geophysical Research，1989，94（B2）：1913－1925.

［24］McNeice G W，Jones A G.Multisite，multifrequency tensor decomposition of magnetotelluric data.Geophysics，2001，66（1）：158－173.

［25］Rodi W， Mackie R L. Nonlinear conjugate gradients algorithm for 2－D magnetotelluric inversion.Geophysics，2001，66（1）：174－187.

［26］蔡軍濤，陳小斌.大地電磁資料精細(xì)處理和二維反演解釋技術(shù)研究（二）——反演數(shù)據(jù)極化模式選擇.地球物理學(xué)報(bào)，2010，53（11）：2703－2714.Cai J T，Chen X B.Refined techniques for data processing and two－dimensional inversion in magnetotelluric II：Which data polarization mode should be used in 2Dinversion.Chinese J.Geophys.（in Chinese），2010，53（11）：2703－2714.

［27］Ledo J.2－D versus 3－D magnetotelluric data interpretation.Surveys in Geophysics，2005，26（5）：511－543.

［28］Berdichevsky M N.Marginal Notes On Magnetotellurics.Surveys in Geophysics，1999，20（3－4）：341－375.

［29］Zhang Z，Wang Y，Chen Y，et al.Crustal structure across Longmenshan fault belt from passive source seismic profiling.Geophys.Res.Lett.，2009，36（17）：L17310.

［30］Robert A，Pubellier M，de Sigoyer J，et al.Structural and thermal characters of the Longmen Shan （Sichuan，China）.Tectonophysics，2010，491（1－4）：165－173.

［31］Wang C，Gao R，Yin A，et al.A mid－crustal strain－transfer model for continental deformation：A new perspective from high－resolution deep seismic－reflection profiling across NE Tibet.Earth and Planetary Science Letters，2011，306（3－4）：279－288.

［32］徐錫偉，聞學(xué)澤，陳桂華等.巴顏喀拉地塊東部龍日壩斷裂帶的發(fā)現(xiàn)及其大地構(gòu)造意義.中國(guó)科學(xué)（D輯：地球科學(xué)），2008，38（5）：529－542.Xu X W，Wen X Z，Chen G H，et al.Discovery of the Longriba faults，eastern part of the Bayankela tectonic block and its geodynamic implications.Sciences in China （Series D：Earth Sciences）（in Chinese），2008，38（5）：529－542.

［33］Burchfiel B C，Chen Z L，Yupinc L，et al.Tectonics of the Longmen Shan and Adjacent Regions，Central China.International Geology Review，1995，37（8）：661－735.

［34］王緒本，朱迎堂，趙錫奎等.青藏高原東緣龍門山逆沖構(gòu)造深部電性結(jié)構(gòu)特征.地球物理學(xué)報(bào)，2009，52（2）：564－571.Wang X B，Zhu Y T，Zhao X K，et al.Deep conductivity characteristics of the Longmen Shan，Eastern Qinghai－Tibet Plateau.Chinese J.Geophys.（in Chinese），2009，52（2）：564－571.

［35］王偉濤，賈東，李傳友等.四川龍泉山斷裂帶變形特征及其活動(dòng)性初步研究.地震地質(zhì)，2008，30（4）：968－979.Wang W T，Jia D，Li C Y，et al.Preliminary investigation on deformation characteristics and activity of Longquanshan fault belt in Sichuan.Seismology and Geology （in Chinese），2008，30（4）：968－979.

［36］萬天豐.中國(guó)大地構(gòu)造學(xué).北京：地質(zhì)出版社，2011：258－262.Wan T F.The Tectonics of China：Data， Maps and Evolution.New York：Springer，2011：300－304.

［37］劉啟元，李昱，陳九輝等.汶川MS8.0地震：地殼上地幔S波速度結(jié)構(gòu)的初步研究.地球物理學(xué)報(bào)，2009，52（2）：309－319.Liu Q Y，Li Y，Chen J H，et al.Wenchuan MS8.0 earthquake： preliminary study of the S－wave velocity structure of the crust and upper mantle.Chinese J.Geophys.（in Chinese），2009，52（2）：309－319.

［38］張季生，高銳，曾令森等.龍門山及鄰區(qū)重、磁異常特征及與地震關(guān)系的研究.地球物理學(xué)報(bào)，2009，52（2）：572－578.Zhang J S，Gao R，Zeng L S，et al.Relationship between characteristics of gravity and magnetic anomalies and the earthquakes in Longmenshan range and adjacent areas.Chinese J.Geophys.（in Chinese），2009，52（2）：572－578.

［39］許志琴，楊經(jīng)綏，姜枚等.大陸俯沖作用及青藏高原周緣造山帶的崛起.地學(xué)前緣，1999，6（3）：139－151.Xu Z Q，Yang J S，Jiang M，et al.Continental subduction and uplifting of the orogenic belts at the margin of the Qinghai－Tibet Plateau.Earth Science Frontiers （in Chinese），1999，6（3）：139－151.

［40］吳建平，黃媛，張?zhí)熘械?汶川MS8.0級(jí)地震余震分布及周邊區(qū)域P波三維速度結(jié)構(gòu)研究.地球物理學(xué)報(bào)，2009，52（2）：320－328.Wu J P，Huang Y，Zhang T Z，et al.Aftershock distribution of the MS8.0Wenchuan earthquake and three dimensional P－wave velocity structure in and around source region.Chinese J.Geophys.（in Chinese），2009，52（2）：320－328.

［41］Zhao W J，Kumar P， Mechie J，et al.Tibetan plate overriding the Asian plate in central and northern Tibet.Nature Geoscience，2011，4（12）：870－873.

［42］Hubbard J，Shaw J H.Uplift of the Longmen Shan and Tibetan plateau，and the 2008Wenchuan（M＝7.9）earthquake.Nature，2009，458（7235）：194－197.