熊 維 譚 凱 喬學(xué)軍 王 偉

1 中國地震局地震研究所（地震大地測量重點實驗室），武漢市洪山側(cè)路40號，430071

2 中國地震局地殼應(yīng)力研究所武漢科技創(chuàng)新基地，武漢市洪山側(cè)路40號，430071

李強等［1］處理了1999～2007網(wǎng)絡(luò)工程GPS數(shù)據(jù)和2009、2011年陸態(tài)網(wǎng)絡(luò)GPS數(shù)據(jù)，獲得較新的中國大陸地殼運動速度場。陸態(tài)網(wǎng)絡(luò)在2013年又觀測了1期，增強了速度場分析的精度和可靠性。然而，2008年龍門山斷裂中段發(fā)生汶川Mw7.9地震，龍門山斷裂的解鎖和地震釋放的巨大應(yīng)力勢必對龍門山周邊區(qū)域構(gòu)造運動產(chǎn)生重要影響。

本文基于1999～2013陸態(tài)網(wǎng)絡(luò)工程GPS觀測資料，結(jié)合地質(zhì)學(xué)和地球物理學(xué)資料，通過對汶川地震前后研究區(qū)域地塊形變模式和斷層活動性狀的研究，分析了研究區(qū)現(xiàn)今構(gòu)造運動模式與應(yīng)變分配及其與區(qū)域地震活動性的關(guān)系，并討論了汶川地震對上述過程的影響。

1 活動塊體劃分與活動性狀

1.1 活動塊體劃分

龍門山地區(qū)斷層十分發(fā)育，地震頻發(fā)，可分為巴顏喀拉地塊（包括東部的龍門山塊體）、川西北地塊（分為南北兩部分）、滇中地塊、華南地塊和柴達木地塊等數(shù)個次級塊體（圖1）。

1.2 活動塊體運動學(xué)特征

1.2.1 塊體整體旋轉(zhuǎn)與均勻應(yīng)變模型

本文參考趙靜等［2］的做法，利用塊體整體旋轉(zhuǎn)與均勻應(yīng)變模型描述各活動地塊。塊體旋轉(zhuǎn)與均勻應(yīng)變模型［3］認為，塊體運動包括剛性旋轉(zhuǎn)和塊體內(nèi)部變形，且塊體內(nèi)部變形是均勻的。觀測點的運動速度為由相鄰塊體的剛性旋轉(zhuǎn)引起的速度與塊體內(nèi)部均勻變形引起的速度之和。其公式如下：

其中，

式中，λ、φ分別為觀測點經(jīng)、緯度，Ve、Vn為東方向速度和北方向速度，ωx、ωy、ωz為塊體的歐拉矢量，r為塊體的平均曲率半徑，λ0、φ0為塊體的幾何中心經(jīng)度和幾何中心緯度，εe、εne、εen、εn為區(qū)域內(nèi)應(yīng)變率張量。由于塊體應(yīng)變均勻，因此εe、εne、εen、εn為常數(shù)且εen＝εne。式（1）內(nèi)有ωx、ωy、ωz、εe、εne、εn6個未知量，因此需要3個以上的GPS站點數(shù)據(jù)進行解算。

1.2.2 塊體邊界差異運動

塊體邊界斷層的長期滑移速率——塊體邊界差異運動速率可以通過兩個塊體的旋轉(zhuǎn)在其邊界處的速度差得到，而這兩個塊體的運動都可以通過球坐標下的歐拉矢量表示。兩個相鄰塊體i和j邊界斷層上長期滑移速率計算公式如下：

式中，Ωi、Ωj分別為i、j板塊的歐拉矢量，Ωij為兩個板塊的相對運動。利用式（3）可由兩個板塊的絕對運動求得這兩個板塊的相對運動。

在板塊運動模型已知的情況下，可以求得兩個相鄰板塊邊界上的相對運動：

式中，r為測站的位置矢量，為r處兩個板塊的邊界差異運動速率。由活動地塊的歐拉運動矢量和式（3）、（4），計算活動地塊邊界相對運動（表1）。

圖1 研究區(qū)主要活動斷裂、塊體劃分與速度場Fig.1 Distribution map of main active faults，blocksand GPS－derived horizontal velocity relative to Eurasia plate in the study region

表1 研究區(qū)活動塊體邊界差異運動Tab.1 The differential movement on block boundary

2 研究區(qū)主要斷裂活動性狀

斷層滑動速率和斷層閉鎖深度是描述斷裂活動性狀的重要依據(jù)，同時也是地震危險性分析的重要基礎(chǔ)數(shù)據(jù)。本文基于1999～2007和2009～2013GPS速度場資料，以主要塊體間差異運動結(jié)果為約束，利用基于Okada 模型的Rngchn 程序［6］，反演研究區(qū)主要活動斷裂在汶川地震前后的滑動速率與斷層閉鎖深度，部分反演模擬值與實際觀測值的比較見圖2。下面簡要分析主要斷裂的反演結(jié)果（表2，圖3）。

表2 研究區(qū)主要活動斷裂滑動速率與斷層閉鎖深度Tab.2 Slip rate and locking depth of main active faults

圖2 反演模擬值與觀測值的比較（藍色箭頭表示速度場觀測值，紅色箭頭表示速度場模擬值）Fig.2 Comparison between analog value and observed value（The blue arrows represent observations while the red arrow indicates the value of the simulation）

圖3 汶川地震前后斷裂滑動速率（藍色箭頭表示傾滑分量，紅色箭頭表示走滑分量）Fig.3 Slip rate of main active faults before and after Wenchuan earthquake（The blue arrows represent dip－slip rate while the red arrow indicates strike－slip rate）

從反演結(jié)果來看，塊體整體旋轉(zhuǎn)和內(nèi)部均勻應(yīng)變模型與塊體運動和邊界斷層彈性位錯模型所得斷裂活動速率大部分相差不大，差異較大的區(qū)域集中在昆侖斷裂西部，除了模型差異原因，主要是數(shù)據(jù)量較少所致。塊體運動和內(nèi)部均勻變形模型的優(yōu)點是根據(jù)大尺度的宏觀場給出結(jié)果，結(jié)果穩(wěn)定性好；而塊體運動和邊界彈性位錯模型考慮了邊界閉鎖深度和傾滑運動的影響，理論比較嚴密。下面主要根據(jù)后者的結(jié)果進行討論分析。

2.1 東昆侖斷裂

東昆侖斷裂是柴達木地塊與巴顏喀拉地塊的分界，以左旋走滑為主，近代地震活動性很強?；谛l(wèi)星影像解譯和野外考察測量，李陳俠等［9］認為東昆侖斷裂西段（托索湖斷裂）平均水平速率為10.8±1mm／a，中段（瑪沁斷裂）平均水平滑動速率為9.3±2mm／a，東段（瑪曲斷裂）平均水平滑動速率為4.9±1.3mm／a。本文結(jié)果較之略大，這可能是由于不同的研究方法引起的（李陳俠等利用的是地質(zhì)學(xué)方法）。本文結(jié)果與文獻結(jié)果均揭示，東昆侖斷裂水平滑動速率從西至東呈梯度下降，但是并沒有出現(xiàn)Kirby所認為的150km 內(nèi)滑動速率迅速衰減的現(xiàn)象［10］。

2.2 鮮水河－安寧河－則木河－小江斷裂

鮮水河－安寧河－則木河－小江斷裂是研究區(qū)活動性最強的左旋走滑斷裂帶，由北到南依次為鮮水河斷裂、安寧河斷裂、則木河斷裂和小江斷裂。王閻昭等［8］利用連接斷層元模型給出鮮水河斷裂、安寧河斷裂、則木河斷裂和小江斷裂左旋走滑速率依次為8.9～17.1mm／a、5.1±2.5mm／a、2.8±2.3mm／a和9.4±1.2mm／a；丁開華等［7］利用斷層速度剖面法給出的以上斷層左旋走滑速率分別為10.0±1.0 mm／a、11.4±0.4 mm／a、7.1±0.5 mm／a和7.8±0.8 mm／a。本文結(jié)果與上述文獻結(jié)果基本吻合。

魏文薪等［11］認為，鮮水河南段、安寧河、則木河斷裂閉鎖程度相對于川滇地區(qū)其他活動斷裂較高。由本文結(jié)果可看出，鮮水河瀘定－石棉段、安寧河和則木河斷裂震后斷層閉鎖深度均大于20km。

2.3 龍門山斷裂

趙靜等［2］給出了龍門山斷裂在汶川地震前后的滑動速率與閉鎖深度，其結(jié)果與本文十分吻合。龍門山斷裂南段震后滑動速率與斷層閉鎖深度均明顯上升，說明汶川地震發(fā)生后，龍門山斷裂北段發(fā)生破裂并解鎖，其同震擾動和震后蠕滑作用加快了龍門山斷裂南段的應(yīng)力積累，并最終促使了2013年蘆山地震的發(fā)生。

3 研究區(qū)構(gòu)造運動與應(yīng)變分配

3.1 巴顏喀拉塊體

本文結(jié)果與文獻結(jié)果均揭示東昆侖斷裂水平滑動速率呈現(xiàn)出西強東弱的特征，由北西托索湖段的13mm／a左右逐漸下降至南東瑪曲段的6.7 mm／a左右，說明青藏高原物質(zhì)向東擠出作用由西向東連續(xù)而漸變地減弱。由于受到華南塊體阻擋，巴顏喀拉塊體東部在龍門山地區(qū)形成巨大的推覆構(gòu)造。在此作用下，東昆侖斷裂走向在其東端塔藏段迅速偏轉(zhuǎn)為近南北向，并與岷江斷裂、龍日壩斷裂、虎牙斷裂交匯，在此，沿東昆侖斷裂帶向東逸出的青藏高原物質(zhì)部分被地殼縮短所吸收，驅(qū)使岷山－龍門山地區(qū)山體抬升。

東昆侖斷裂閉鎖深度由西向東逐漸升高。西段在經(jīng)過了瑪尼地震、昆侖山口西地震后，應(yīng)力基本釋放完畢，處于震間蠕滑階段，應(yīng)力并未充分積累，滑動速率較快且閉鎖深度較低。東昆侖斷裂東部瑪曲段是該斷裂著名和唯一的地震空區(qū)，缺乏歷史地震使得該地區(qū)應(yīng)力長期積累，這種狀況在該地區(qū)較低的滑動速率和較高的閉鎖深度上有所體現(xiàn)。

受汶川地震的影響，東昆侖斷裂的滑動速率在震后有所上升，且東段增加量大于中段和西段，說明汶川地震的發(fā)生加速了巴顏喀拉塊體向南東方向的運動，且距離汶川地震震中越近，這種加速作用越大。

3.2 川滇菱形塊體

3.2.1 川西北塊體

鮮水河斷裂滑動速率結(jié)果呈現(xiàn)出與東昆侖斷裂相似的西強東弱的特征，其走向在接近華南地塊部分顯著向南北向偏轉(zhuǎn)，說明沿羌塘、川西北塊體逸出的高原物質(zhì)運動方向由西向東逐漸向南偏轉(zhuǎn)。川西北塊體受到來自羌塘塊體的向東推擠和華南塊體的強烈阻擋作用，其運移方向被迫向南偏轉(zhuǎn)，并在川滇菱形塊體中部遇到雜巖體的阻礙，形成帶有較大逆沖分量的麗江－小金河斷裂。

龍門山斷裂與鮮水河斷裂交匯于瀘定附近，鮮水河斷裂的爐霍－瀘定段處于巴顏喀拉塊體與川西北塊體的交界處，位于龍門山斷裂與麗江－小金河斷裂的西側(cè)，瀘定－石棉段位于龍門山斷裂以東、華南塊體與川西北塊體的交界處。比較汶川地震前后的結(jié)果發(fā)現(xiàn)，鮮水河爐霍－瀘定段滑動速率在震后明顯下降，而瀘定－石棉段在震后上升。此外，麗江－小金河斷裂滑動速率逆沖分量在汶川震后亦顯著上升。這說明在汶川地震后，巴顏喀拉塊體相對于羌塘塊體和川西北塊體加速向東運動，而川西北塊體亦受到汶川地震影響，相對于穩(wěn)定的華南塊體加速向南東方向運動。

3.2.2 滇中塊體

川西北塊體向南東方向的滑移量經(jīng)過麗江－小金河斷裂的吸收與調(diào)整，傳遞給位于其南部的滇中塊體。安寧河－則木河－小江斷裂系是滇中塊體的東邊界，其走向由南東逐漸偏轉(zhuǎn)向近南北，證明了滇中塊體向南南東向的運移和自身的順時針旋轉(zhuǎn)。滇中塊體的自旋作用使塊體邊界存在區(qū)域應(yīng)變響應(yīng)，形成小江斷裂系南部張性傾滑分量和北部雁列的多條北東向右旋走滑斷裂。值得一提的是，其中的包谷垴－小河斷裂即是2014－08－03云南魯?shù)榈卣鸬脑姓饠鄬?。此外，安寧河－則木河－小江斷裂系的傾滑分量由北段的擠壓逐漸變化至南段的拉張。

比對汶川地震庫侖應(yīng)力結(jié)果［12］可以看出，同震釋放的應(yīng)力對滇中地區(qū)影響已非常微弱。然而安寧河、則木河和小江斷裂在汶川震后滑動速率均有所上升，說明在羌塘和川西北塊體的加速推擠下，滇中塊體的運移速率有所增加。

4 結(jié) 語

本文基于1999～2013陸態(tài)網(wǎng)絡(luò)工程GPS觀測資料，分析了汶川地震前后龍門山周邊地區(qū)活動地塊運動學(xué)特征和主要斷裂活動性狀。結(jié)果顯示，巴顏喀拉、川西北和滇中塊體在向南東運移的同時，自身作順時針旋轉(zhuǎn)；東昆侖、鮮水河等大型走滑斷裂左旋速率自西向東逐漸降低，在其南東段走向向近南北向偏轉(zhuǎn)；安寧河－則木河－小江斷裂系的傾滑分量由北段的擠壓逐漸變化至南段的拉張。這些現(xiàn)象說明，在印度板塊向歐亞大陸的北北東向俯沖和青藏高原東西向的重力滑塌作用下，青藏高原物質(zhì)的東向擠出遇到華南地塊強硬阻擋后，其運移方向發(fā)生向南偏轉(zhuǎn)，造就了東昆侖、鮮水河等大型走滑斷裂和龍門山斷裂、麗江－小金河斷裂等逆沖構(gòu)造，并造成活動地塊邊界的多處應(yīng)變響應(yīng)；在多條逆沖斷層的吸收和調(diào)整作用下，高原物質(zhì)運動幅度由西向東、由北至南逐漸降低，并驅(qū)使逆沖帶區(qū)域地殼增厚與山體抬升。

2008汶川地震使龍門山斷裂解鎖，龍門山斷裂、麗江－小金河斷裂的逆沖速率和安寧河－則木河－小江斷裂走滑速率顯著上升，巴顏喀拉、羌塘、川西北塊體加速向南東運移，并推動滇中塊體的運動。這種現(xiàn)象表明，研究區(qū)地殼構(gòu)造運動是一個完整有機的系統(tǒng)，不同地塊間通過邊界斷裂的錯動進行相互作用與調(diào)整。

［1］李強，游新兆，楊少敏，等.中國大陸構(gòu)造變形高精度大密度GPS監(jiān)測－現(xiàn)今速度場［J］.中國科學(xué)：地球科學(xué)，2012，42（5）：629－632（Li Qiang，You Xinzhao，Yang Shaomin，et al.A Precise Velocity Field of Tectonic Deformation in China as Inferred from Intensive GPS Observations［J］.Science China Earth Science，2012，42（5）：629－632）

［2］趙靜，江在森，武艷強，等.汶川地震前龍門山斷裂帶閉鎖程度和滑動虧損分布研究［J］.地球物理學(xué)報，2012，55（9）：2 963－2 972（Zhao Jing，Jiang Zaisen，Wu Yanqiang，et al.Study on Fault Locking and Fault Slip Deficit of the Longmenshan Fault Zone before the Wenchuan Earthquake［J］.Chinese J Geophys，2012，55（9）：2 963－2 972）

［3］李延興，黃珹，胡新康，等.板內(nèi)塊體的剛性彈塑性運動模型與中國大陸主要塊體的應(yīng)變狀態(tài)［J］.地震學(xué)報，2001，24（6）：565－572（Li Yanxing，Huang Cheng，Hu Xinkang，et al.The Rigidand Elastic－Plastic Model of the Blocks in Intro－Plate and Strain Status of Principal Blocks in the Continent of China［J］.Actaseismologica Sinica，2001，24（6）：565－572）

［4］Yang S M，Wang W，Wang Q.Crustal Block Rotations in Chinese Mainland Revealed by GPS Measurements［J］.Earthquake Science，2009，22（6）：639－649

［5］喬學(xué)軍，王琪，杜瑞林，等.川滇地區(qū)活動地塊現(xiàn)今地殼形變特征［J］.地球物理學(xué)報，2004，47（5）：805－811（Qiao Xuejun，Wang Qi，Du Ruilin.Characteristics of Current Crustal Deformation of Active Blocks in the Sichuan－Yunnan Region［J］.Chinese J Geophys，2004，47（5）：805－811）

［6］Feigl K L，DupréE.RNGCHN：A Program to Calculate Displacement Components from Dislocations in an Elastic Half－Space with Applications for Modeling Geodetic Measurements of Crustal Deformation［J］.Computers ＆Geosciences，1999，25（6）：695－704

［7］丁開華，許才軍，鄒蓉，等.利用GPS分析川滇地區(qū)活動地塊運動與應(yīng)變模型［J］.武漢大學(xué)學(xué)報：信息科學(xué)版，2013，38（7）：822－827（Ding Kaihua，Xu Caijun，Zhou Rong，et al.Crustal Movement and Strain Model of Active Blocks Analyzed by GPS in Sichuan－Yunnan Region［J］.Geomatics and Information Science of Wuhan University，2013，38（7）：822－827）

［8］王閻昭，王恩寧，沈正康，等.基于GPS資料約束反演川滇地區(qū)主要斷裂活動速率［J］.中國科學(xué)D 輯：地球科學(xué)，2008，38（5）：582－597（Wang Yanzhao，Wang Enning，Shen Zhengkang，et al.GPS－Constrained Inversion of Present－Day Slip Rates along Major Faults of the Sichuan－Yunnan Region，China［J］.Science in China（Series D）：Earth Science，2008，38（5）：582－597）

［9］Li C X，Xu X W，Wen X Z，et al.Rupture Segmentation and Slip Partitioning of the Mid－Eastern Part of the Kunlun Fault，North Tibetan Plateau［J］.Science China Earth Science，2011，54（11）：1 730－1 745

［10］Kirby E，Harkins N，Wang E Q，et al.Slip Rate Gradients Along the Eastern Kunlun Fault［J］.Tectonics，2007，26（2）：485－493

［11］魏文薪，江在森，武艷強，等.利用GPS數(shù)據(jù)研究川滇塊體東邊界主要斷裂帶運動特性［J］.武漢大學(xué)學(xué)報：信息科學(xué)版，2012，37（9）：1 041－1 044（Wei Wenxin，Jiang Zaisen，Wu Yanqiang，et al.Motion Characteristics of Major Faults in East Boundary of Sichuan－Yunnan Block Obtained with GPS Data［J］.Geomatics and Information Science of Wuhan University，2012，37（9）：1 041－1 044）

［12］單斌，熊熊，鄭勇，等.2008年5月12日Mw7.9汶川地震導(dǎo)致的周邊斷層應(yīng)力變化［J］.中國科學(xué)D 輯：地球科學(xué)，2009，39（5）：537－545（Shan Bin，Xiong Xiong，Zheng Yong，et al.Stress Changes on Major Faults Caused by Mw7.9 Wenchuan Earthquake，May 12，2008［J］.Science in China（Series D）：Earth Science，2009，39（5）：537－545）