蔣鋒云 朱良玉 李玉江

摘要：在考慮川滇地區(qū)地殼介質(zhì)的橫向分塊、縱向分層特征，及汶川地震同震斷層面靜態(tài)位錯(cuò)上下盤(pán)非對(duì)稱(chēng)性的基礎(chǔ)上，構(gòu)建川滇地區(qū)地殼三維粘彈性有限元模型，研究了汶川地震對(duì)川滇地區(qū)的震后影響。結(jié)果表明：（1）震后l0 a由于介質(zhì)粘彈性松弛效應(yīng)產(chǎn)生的遠(yuǎn)場(chǎng)水平變形在川滇菱形地塊內(nèi)部約0～20mm，對(duì)其北部影響較大，南部較刁、，其產(chǎn)生的遠(yuǎn)場(chǎng)垂直變形較小，在川滇菱形地塊大部分區(qū)域表現(xiàn)為0～4mm的上升：（2）離發(fā)震斷層較近的鮮水河斷裂、東昆侖斷裂，震后水平和垂直形變呈現(xiàn)高梯度帶：（3）為了和汶川地震之后川滇地區(qū)強(qiáng)震活動(dòng)進(jìn)行對(duì)比，計(jì)算了主要塊體邊界活動(dòng)斷裂水平應(yīng)力狀態(tài)和庫(kù)侖破裂應(yīng)力變化。從同震及震后位移場(chǎng)空間分布特征、應(yīng)力狀態(tài)反映的斷層活動(dòng)特征及庫(kù)侖破裂應(yīng)力變化反映的斷層面應(yīng)力加卸載情況來(lái)看，與該區(qū)域強(qiáng)震活動(dòng)分布在空間上存在一定的相關(guān)性。

關(guān)鍵詞：汶川地震；有限元；粘彈性松弛；川滇地區(qū)

中圖分類(lèi)號(hào)：P315.725 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-0666（2018）02-0233-11

0 引言

自前，普遍的認(rèn)識(shí)是一次大地震震后短期發(fā)震斷層近場(chǎng)的地殼變形主要由發(fā)生在斷層面及其延伸部分的滑移（震后余滑）和孔隙介質(zhì)反彈控制（Jouanne et al，2011；Jonsson et al，2003），而震后長(zhǎng)期的遠(yuǎn)場(chǎng)變形則主要受脆性上地殼下部的粘性下地殼和上地幔內(nèi)部同震應(yīng)力變化導(dǎo)致的粘彈性松弛效應(yīng)控制（Rice，Gu，1983）。而有關(guān)研究表明，由于介質(zhì)粘彈性松弛效應(yīng)引起的斷層面應(yīng)力變化（庫(kù)侖破裂應(yīng)力變化）甚至比同震影響更為顯著，可以導(dǎo)致周?chē)鷶鄬邮Х€(wěn)，從而觸發(fā)強(qiáng)震。如沈正康等（2003）模擬了東昆侖斷裂帶1937年以來(lái)的5次7級(jí)以上強(qiáng)震在黏彈性成層介質(zhì)中地震斷層錯(cuò)動(dòng)產(chǎn)生的應(yīng)力演化過(guò)程，并計(jì)算了在后續(xù)地震破裂面上產(chǎn)生的庫(kù)侖破裂應(yīng)力變化，發(fā)現(xiàn)前4次地震均使2001年昆侖山口西8.1級(jí)大震發(fā)震斷層面庫(kù)侖破裂應(yīng)力增加，且中地殼和下地殼的黏彈性松弛效應(yīng)使得庫(kù)侖破裂應(yīng)力場(chǎng)隨著時(shí)間的推移而逐漸加強(qiáng)，前3個(gè)地震震后松弛效應(yīng)相比同震影響更為顯著。萬(wàn)永革等（2007）采用多層彈粘性介質(zhì)中的位錯(cuò)模型，以M≥7.0大地震和GPS數(shù)據(jù)得出的長(zhǎng)期構(gòu)造加載為形變?cè)矗?jì)算得到了地震產(chǎn)生的應(yīng)力變化，給出了青藏高原東北部1920年以來(lái)積累庫(kù)侖破裂應(yīng)力演化，結(jié)果表明：17次大地震均發(fā)生在庫(kù)侖破裂應(yīng)力變化為正的區(qū)域，觸發(fā)率達(dá)85%。

2008年5月12日汶川8.0級(jí)特大地震發(fā)生至今已接近10年，有關(guān)該地震引起的震后遠(yuǎn)場(chǎng)變形的研究不多，少數(shù)的研究（談洪波等，2009；李強(qiáng)等，2013）局限于利用彈性位錯(cuò)模型，僅考慮介質(zhì)的縱向分層而沒(méi)有考慮研究區(qū)復(fù)雜的介質(zhì)橫向非均勻性，而川滇地區(qū)相比華南地塊普遍存在流變性更強(qiáng)的中下地殼，對(duì)震后形變場(chǎng)產(chǎn)生重要的影響（季靈運(yùn)等，2015）。相比位錯(cuò)模型，有限元技術(shù)由于能充分考慮地球介質(zhì)的非均勻性、調(diào)整斷層面上下盤(pán)同震靜態(tài)位移而在數(shù)值模擬中得到廣泛的應(yīng)用。因此，本文借助有限元技術(shù)，綜合考慮活動(dòng)地塊、地球物理、大地測(cè)量、地震活動(dòng)等多學(xué)科資料，構(gòu)建考慮地殼介質(zhì)縱向分層、橫向分塊的三維粘彈性有限元模型，以USGS給出的汶川地震同震破裂（Ji，Hayes，2008）作為加載，研究汶川地震對(duì)川滇地區(qū)由介質(zhì)粘彈性松弛效應(yīng)引起的震后形變場(chǎng)、應(yīng)力場(chǎng)的時(shí)空分配關(guān)系。

汶川地震之后，川滇地區(qū)中強(qiáng)地震活動(dòng)相比地震前出現(xiàn)了一些顯著的變化，如鮮水河斷裂（3次MS≥5.0地震）、大涼山地塊附近（除蘆山地震外，發(fā)生7次MS≥5.0地震）強(qiáng)震活動(dòng)較為活躍。本文利用有限元模擬結(jié)果計(jì)算了汶川地震同震及震后不同時(shí)段研究區(qū)域主要塊體邊界活動(dòng)斷裂應(yīng)力狀態(tài)和庫(kù)侖應(yīng)力變化，并結(jié)合震后形變應(yīng)變場(chǎng)時(shí)空演化規(guī)律，進(jìn)一步分析了汶川地震對(duì)川滇地區(qū)的震后影響和地震活動(dòng)之間的相關(guān)性。

川滇地區(qū)三維粘彈性有限元模型

根據(jù)川滇地區(qū)全新世活動(dòng)斷裂、地殼上地幔三維波速結(jié)構(gòu)等資料，充分考慮對(duì)地質(zhì)構(gòu)造運(yùn)動(dòng)和地震活動(dòng)起決定作用的活動(dòng)斷裂和活動(dòng)斷塊，建立川滇地區(qū)三維地質(zhì)構(gòu)造模型。

1.1 介質(zhì)分區(qū)及物性參數(shù)

川滇地區(qū)的活動(dòng)斷裂控制了主要構(gòu)造活動(dòng)的發(fā)生和發(fā)展，同時(shí)成為活動(dòng)地塊的邊界（張培震等，2003），依據(jù)川滇地區(qū)的構(gòu)造特征模型劃分出10個(gè)塊體，具體見(jiàn)圖1。

構(gòu)建有限元幾何模型時(shí)，模型橫向分塊主要依照前述10個(gè)塊體劃分結(jié)果（圖1）設(shè)定，縱向分層及各塊體介質(zhì)參數(shù)主要根據(jù)川滇地區(qū)地震波速資料揭示的地殼上地幔結(jié)構(gòu)和P波、S波速度結(jié)果設(shè)定（王椿墉等，2015；吳建平等，2006，2009），具體設(shè)定為：（1）吳建平等（2006，2009）研究結(jié)果表明川滇地區(qū)中下地殼普遍存在低速層，而華南地塊并不明顯。因此，在介質(zhì)參數(shù)設(shè)定時(shí)，考慮了中下地殼低速體的存在及其橫向差異性，將巴顏喀拉塊體、川西北塊體、滇中塊體、大涼山塊體、保山塊體、景谷塊體6個(gè)塊體中下地殼粘滯系數(shù)設(shè)定為1.0×1018Pa·s，而穩(wěn)定的華南地塊設(shè)定為1.0×1021Pa·s，其他區(qū)域?yàn)?.0×1020Pa·s；（2）已有研究（王達(dá)鏞等，2015；吳建平等，2006，2009）表明，南北地震帶中南段上地殼統(tǒng)一為強(qiáng)度較高的介質(zhì)，橫向和構(gòu)造存在顯著的相關(guān)性。因此，上地殼統(tǒng)一設(shè)定為較硬且粘滯系數(shù)較高的脆性層，但在彈性模量的設(shè)置上體現(xiàn)了華南地塊相比其他區(qū)域介質(zhì)相對(duì)較硬的特性：（3）整個(gè)區(qū)域上地幔具有相同的介質(zhì)參數(shù)，王輝等（2007）根據(jù)地震波速進(jìn)行統(tǒng)計(jì)計(jì)算得，粘滯系數(shù)約為1.0×1020Pa·s：（4）由于需要在龍門(mén)山斷裂上加載汶川地震的同震位移，而同震滑動(dòng)可以看作無(wú)摩擦的瞬間滑動(dòng)，因此龍門(mén)山發(fā)震斷層采用無(wú)摩擦的接觸面來(lái)設(shè)計(jì)，具體的斷層深度、傾角的設(shè)定和滑動(dòng)分布的加載依據(jù)Ji和Hayes（2008）給出的有限斷層面解來(lái)進(jìn)行：（5）為簡(jiǎn)化模型，研究區(qū)其他斷裂均采用傾角為90°的直立斷層，參考王輝等（2007）的做法，斷層的設(shè)定用粘滯系數(shù)較小（1.0×1020Pa·s）的軟弱體代替，具體模型分塊及分層參數(shù)見(jiàn)表1。為了盡可能地減少模型邊界效應(yīng)對(duì)計(jì)算結(jié)果的影響，我們將研究區(qū)域向外圍擴(kuò)展了200km。在上述地質(zhì)模型的基礎(chǔ)上，借助ANSYS有限元軟件建立川滇地區(qū)三維粘彈性有限元模型（圖2）。模型研究范圍為（94°～110°E，21°～34°N），垂直區(qū)域?yàn)榈乇碇辽系伥犿敳?00km深處，采用SOLID185和SOLID186這2種實(shí)體單元進(jìn)行劃分，整個(gè)有限元模型共劃分555330個(gè)單元，737989個(gè)節(jié)點(diǎn)。

1.2 汶川地震同震發(fā)震斷層模型

有關(guān)汶川地震的同震斷層破裂模型已發(fā)表的研究成果有很多，Ji和Hayes（2008）的模型僅考慮中遠(yuǎn)場(chǎng)地震波，斷層僅為一個(gè)面，比較簡(jiǎn)單。陳運(yùn)泰等（2008）的研究結(jié)果雖然也只有一個(gè)斷層面，但在資料的使用上，增加了近場(chǎng)的地震波資料約束。而之后的地震破裂模型由于考慮了大地測(cè)量、地震波、InSAR等多種觀測(cè)（畢研磊等，2017；Shen et al，2009；Wang et al，2011；Field-ing et al，2013；Wan et al，2017），斷層分了很多段，且考慮了斷層深淺部?jī)A角的差異，相對(duì)比較復(fù)雜，也更為精確。談洪波等（2009）對(duì)陳運(yùn)泰等（2008）以及ii和Hayes（2008）研究中的不同模型進(jìn)行比較，發(fā)現(xiàn)不同模型同震及震后影響差異主要集中在斷層附近，中遠(yuǎn)場(chǎng)差異較小?？紤]到本文模型空間尺度范圍較大，且主要考查中遠(yuǎn)場(chǎng)同震及震后形變，為了使有限元模型建立過(guò)程中同震斷層面處理更方便，本文采用Ji和Hayes（2008）的簡(jiǎn)單模型（圖3），這種模型在斷層附近和實(shí)際結(jié)果存在一定的差異，但對(duì)于中遠(yuǎn)場(chǎng)的分析而言，差異可以忽略。圖3中箭頭表示下盤(pán)相對(duì)上盤(pán)的運(yùn)動(dòng)方向，最優(yōu)選取的發(fā)震斷層參數(shù)為走向229°，傾角33°，分別沿走向和傾向均勻的分成21塊和8塊子斷層，每個(gè)子斷層的尺度為15km×5km，并給出了每個(gè)子斷層面的滑動(dòng)角和滑動(dòng)量。

1.3 模型邊界約束與同震加載方式

為了得到汶川地震對(duì)川滇地區(qū)真實(shí)的震后形變場(chǎng)，必須先得到和實(shí)際觀測(cè)接近的同震位移場(chǎng)。Wang等（2011）為我們提供了可靠的、高密度的、大范圍的汶川地震的同震觀測(cè)水平位移場(chǎng)（圖4a），其結(jié)果表明：（1）汶川地震同震水平形變存在明顯的非對(duì)稱(chēng)性，即上盤(pán)（青藏高原一側(cè)）形變強(qiáng)度和范圍明顯大于下盤(pán)（華南地塊一側(cè)），這種現(xiàn)象在INSAR同震近場(chǎng)形變觀測(cè)中也可以看到（Shen et al，2009）：（2）汶川地震同震影響在川滇菱形地塊東邊界以西從南至北，具有和背景構(gòu)造水平運(yùn)動(dòng)場(chǎng)接近的順時(shí)針旋轉(zhuǎn)特征，在滇西北地塊表現(xiàn)為近南北向運(yùn)動(dòng)，而在滇中地塊表現(xiàn)出明顯的南西向整體運(yùn)動(dòng)，遠(yuǎn)至滇西北和滇中塊體南部仍然有將近10mm的同震形變量。

相關(guān)研究（Wang et al，2007；Funning et al，2007；邵志剛等，2008；李寧等，2017；季靈運(yùn)等，2017，焦佳爽等，2017）表明，同震形變場(chǎng)的非對(duì)稱(chēng)性成因機(jī)理主要有以下幾種：（1）發(fā)震斷層兩盤(pán)介質(zhì)性質(zhì)存在顯著差異：（2）震后粘彈性松弛引起的有效粘度發(fā)生橫向變化：（3）斷層存在多條近平行的剪切帶：（4）斷層產(chǎn)狀，即斷層并非直立斷層，存在傾角。對(duì)于汶川地震的發(fā)震斷層龍門(mén)山斷裂而言，斷層兩側(cè)青藏高原相對(duì)華南地塊明顯的介質(zhì)性質(zhì)差異和斷層產(chǎn)狀都可造成同震形變場(chǎng)的非對(duì)稱(chēng)性。而第2種情況，相對(duì)同震短時(shí)間而言，影響可忽略不計(jì)。由于龍門(mén)山斷裂相對(duì)比較集中，第3種情況也可忽略不計(jì)。我們?cè)谟邢拊Ｐ蜆?gòu)建時(shí)就考慮發(fā)震斷層330的北西傾向，也考慮了發(fā)震斷層兩盤(pán)介質(zhì)性質(zhì)的差異。在此基礎(chǔ)上，在邊界條件的設(shè)定上，參照李玉江等（2013）將模型四周和底部做法向約束，其他方向可自由移動(dòng)。目前大多數(shù)同震靜態(tài)位移斷層面加載方式為上下兩盤(pán)斷層面對(duì)稱(chēng)加載（李玉江等，2013），即上下兩盤(pán)同震靜態(tài)位移各占一半。通過(guò)試算，我們發(fā)現(xiàn)這種加載方式對(duì)近場(chǎng)的擬合效果較好，而對(duì)中遠(yuǎn)場(chǎng)擬合效果較差，無(wú)法完全體現(xiàn)出上下兩盤(pán)形變的非對(duì)稱(chēng)性。通過(guò)分析，我們認(rèn)為除了上述4種情況，同震形變的非對(duì)稱(chēng)性可能與斷層兩盤(pán)本身的運(yùn)動(dòng)狀態(tài)有關(guān)，即斷層的主動(dòng)盤(pán)與被動(dòng)盤(pán)。青藏高原受到印度板塊北東向推擠，在外圍構(gòu)造環(huán)境約束下，龍門(mén)山以西表現(xiàn)為垂直龍門(mén)山斷裂近東西向的運(yùn)動(dòng)特征，而華南地塊，特別是穩(wěn)定的四川盆地，相對(duì)青藏高原而言是穩(wěn)定不動(dòng)的。這表明龍門(mén)山斷裂上盤(pán)青藏高原一側(cè)是主動(dòng)盤(pán)，而華南地塊一側(cè)是相對(duì)穩(wěn)定的被動(dòng)盤(pán)。換句話(huà)說(shuō)，華南地塊和青藏高原本身的構(gòu)造邊界約束條件是不一致的，這種不一致性有可能造成上盤(pán)斷層面同震滑移要比下盤(pán)斷層面大，從而使得同震形變?cè)诟蠓秶鷥?nèi)表現(xiàn)出上盤(pán)比下盤(pán)顯著的特征。因此，我們擬通過(guò)調(diào)整斷層上下兩盤(pán)同震靜態(tài)滑移比例來(lái)實(shí)現(xiàn)模擬同震水平位移（主要是中遠(yuǎn)場(chǎng)）和實(shí)際觀測(cè)的最佳匹配，在實(shí)際給定比例時(shí)，要保證上盤(pán)一直大于下盤(pán)，這樣才符合物理機(jī)理。試算結(jié)果表明上下兩盤(pán)同震滑移比例為7.5：2.5時(shí)獲得的同震水平移場(chǎng)（圖4b）和實(shí)際觀測(cè)結(jié)果（圖4a）較為接近。

需要強(qiáng)調(diào)的是雖然受到有限元模型邊界條件、地殼介質(zhì)的簡(jiǎn)化以及發(fā)震斷層產(chǎn)狀與同震滑動(dòng)分布簡(jiǎn)略等影響，模擬結(jié)果和觀測(cè)結(jié)果仍然存在顯著的差異，從模擬值與觀測(cè)值矢量差可以更為清晰地看出（圖4c），差異主要表現(xiàn)為：（1）在斷層附近模擬值比實(shí)際觀測(cè)值要小，這可能是由于使用的同震破裂模型是用遠(yuǎn)場(chǎng)地震波資料反演得到，而沒(méi)有近場(chǎng)的約束引起的：（2）整個(gè)川滇地區(qū)和四川盆地遠(yuǎn)離發(fā)震斷裂地區(qū)，模擬得到速度場(chǎng)比觀測(cè)速度場(chǎng)量值大，且存在系統(tǒng)性偏差。其原因一方面可能是有限元模擬中邊界條件僅做法向約束，和實(shí)際情況可能存在差異，另一方面可能是上下盤(pán)斷層面同震靜態(tài)位移比例調(diào)整，僅用一個(gè)比例系數(shù)，過(guò)于簡(jiǎn)單：（3）整個(gè)區(qū)域除了發(fā)震斷層附近上下兩盤(pán)一定范圍內(nèi)，大部分地區(qū)模擬得到的速度場(chǎng)和實(shí)際觀測(cè)存在10mm以下的殘差。

雖然模擬結(jié)果與觀測(cè)結(jié)果存在上述明顯的差異，但其反映的整個(gè)同震形變場(chǎng)的整體運(yùn)動(dòng)形態(tài)，特別是中遠(yuǎn)場(chǎng)形態(tài)和實(shí)際觀測(cè)結(jié)果（圖4a）具有較好的相似性。因此，利用上述得到的上下盤(pán)同震加載比例及邊界約束方式進(jìn)行汶川地震同震及震后形變場(chǎng)的研究具有一定合理性。

2 汶川地震對(duì)川滇地區(qū)震后影響

采用上述邊界條件及斷層面靜態(tài)同震破裂加載，計(jì)算了汶川地震同震及震后1a、5a、10a這3個(gè)時(shí)間點(diǎn)由于介質(zhì)粘彈性松弛效應(yīng)引起的川滇地區(qū)位移場(chǎng)以及主要塊體邊界活動(dòng)斷裂水平主應(yīng)力狀態(tài)和庫(kù)侖應(yīng)力的變化，并結(jié)合川滇地區(qū)近年來(lái)強(qiáng)震活動(dòng)特征，對(duì)這些物理量的時(shí)空演化特征進(jìn)行分析驗(yàn)證。

從模擬同震水平位移等值線(xiàn)圖（圖5a）可以看出，和實(shí)際觀測(cè)結(jié)果類(lèi)似，汶川地震的同震影響主要沿垂直發(fā)震斷裂兩側(cè)方向延伸，上盤(pán)形變強(qiáng)度和范圍明顯大于下盤(pán)，上下兩盤(pán)存在明顯的非對(duì)稱(chēng)性。同時(shí)平行斷裂方向也存在一定分量，表現(xiàn)出同震形變主要以斷裂逆沖為主，同時(shí)兼具右旋走滑的特征。需要說(shuō)明的是模擬得到的汶川地震對(duì)川滇交界東部地區(qū)同震影響比實(shí)際觀測(cè)結(jié)果量值偏大。

從震后水平形變場(chǎng)來(lái)看，震后1a（圖5b）形變主要集中在發(fā)震斷裂附近，且上盤(pán)變形強(qiáng)度明顯大于下盤(pán)。但就影響范圍來(lái)看，下盤(pán)影響范圍比上盤(pán)要大。震后5a、10a形變（圖5c、d）則在短期形變的基礎(chǔ)上隨時(shí)間推移向外圍擴(kuò)展。其中從川滇交界東部地區(qū)向滇西南、滇西北方向擴(kuò)展特征更為明顯，這和汶川地震之后川滇地區(qū)中強(qiáng)地震活動(dòng)主要呈北東向條帶分布一致。另外，在形變擴(kuò)展過(guò)程中，沿著巴顏喀拉地塊南北兩側(cè)邊界帶出現(xiàn)形變高梯度帶，且隨時(shí)間變化逐漸加強(qiáng)。在這一過(guò)程中，先后在其南邊界鮮水河斷裂上發(fā)生了2010年道孚5.4級(jí)、2011年?duì)t霍5.7級(jí)、2014康定6.3級(jí)強(qiáng)震，近期又在其北邊界東北角發(fā)生了2017年8月8日九寨溝7.0級(jí)強(qiáng)震，可能與汶川地震對(duì)該區(qū)域的震后影響有關(guān)。大涼山地塊震后形變隨著時(shí)間變化增強(qiáng)顯著，水平變形量震后10a達(dá)到20mm，該區(qū)域也是汶川震后強(qiáng)震活躍顯著的區(qū)域。

從垂直形變特征（圖6）來(lái)看，和水平同震形變類(lèi)似，同震垂直形變主要沿著龍門(mén)山斷裂東西兩側(cè)垂直斷層方向延伸，靠近斷層附近下盤(pán)表現(xiàn)為下沉，上盤(pán)表現(xiàn)為上升。距斷層一定范圍內(nèi)（上盤(pán)約300km，下盤(pán)約150km）表現(xiàn)為上升，之后向兩側(cè)延伸又表現(xiàn)為下沉。隨著離開(kāi)斷層距離的增加，同震垂直形變衰減較快。斷裂西南部表現(xiàn)為顯著的上升，而斷裂的東北部表現(xiàn)為上升-下沉交替變化特征。震后1a尺度垂直形變（扣除同震）則主要集中在斷層附近較小的區(qū)域，靠近斷層附近上盤(pán)和同震形變一致繼續(xù)表現(xiàn)為上升，下盤(pán)仍然表現(xiàn)為下沉。相對(duì)遠(yuǎn)離斷層區(qū)域和同震形變相反，表現(xiàn)為下沉區(qū)。隨著時(shí)間的推移，震后垂直形變不斷向外擴(kuò)展，和水平形變類(lèi)似，擴(kuò)展方向除了沿著斷裂兩側(cè)垂直斷裂方向外，還具有從川滇交界東部向滇西北、滇西南方向擴(kuò)展的特征。

從同震及震后主要塊體邊界水平向主應(yīng)力變化特征（圖7）來(lái)看：（1）甘孜—玉樹(shù)斷裂南段、鮮水河斷裂爐霍—道孚—乾寧段以及金沙江斷裂北段同震影響以垂直斷裂方向張應(yīng)力為主，震后隨時(shí)間推移逐漸加強(qiáng)。尤其是鮮水河斷裂，震后10a的張應(yīng)力累積變化量已超過(guò)同震變化量。而該斷裂上相繼發(fā)生了2010年道孚5.4級(jí)、2011年?duì)t霍5.7級(jí)、2014康定6.3級(jí)強(qiáng)震，隨時(shí)間的推移震級(jí)呈增強(qiáng)特征，和汶川地震震后應(yīng)力場(chǎng)變化具有顯著相關(guān)性；（2）安寧河—?jiǎng)t木河—小江斷裂、馬邊斷裂、紅河斷裂、金沙江斷裂南段同震及震后影響均表現(xiàn)為與斷層走向呈一定夾角的莊應(yīng)力變化特征。安寧河—?jiǎng)t木河—小江斷裂主壓應(yīng)力方向（北北東）和背景構(gòu)造主壓應(yīng)力方向（北北西向）相反，使斷層表現(xiàn)為右旋，具有負(fù)影響，汶川震后強(qiáng)震活動(dòng)較弱。紅河斷裂應(yīng)力所反映的斷層走滑活動(dòng)特征和背景特征一致，表現(xiàn)為右旋走滑，但量值相對(duì)較小。金沙江斷裂中南段以垂直斷裂的壓應(yīng)力特征為主：（3）總體上，汶川地震同震及震后對(duì)鮮水河斷裂影響最大，其次為金沙江斷裂北段、安寧河—?jiǎng)t木河斷裂、馬邊斷裂、魯?shù)椤淹〝嗔?、金沙江斷裂中南段、紅河斷裂、麗江—小金河斷裂，對(duì)小江斷裂影響最小。從影響特征和影響強(qiáng)度2方面考慮，汶川震后川滇地區(qū)強(qiáng)震活動(dòng)的主體區(qū)域和其具有較好的一致性。

為了進(jìn)一步分析汶川地震對(duì)川滇地區(qū)主要塊體邊界活動(dòng)斷裂的影響，將汶川地震產(chǎn)生的應(yīng)力變化投影到斷裂滑動(dòng)方向上，用庫(kù)侖破裂應(yīng)力變化來(lái)表示該地震對(duì)周?chē)邮諗嗔旬a(chǎn)生的影響。具體計(jì)算公式參照J(rèn)aeger等（1969）的研究，在計(jì)算中有效摩擦系數(shù)參照Stein等（1992）和King等（1994）的研究，取值為0.4。有關(guān)接收斷裂活動(dòng)性質(zhì)幾何參數(shù)參考鄧起東等（2002）結(jié)果給出。由于研究區(qū)域的地震震源深度絕大多數(shù)在0～25km范圍內(nèi)，因此我們研究汶川地震在10km深度上產(chǎn)生的庫(kù)侖破裂應(yīng)力變化，我們也計(jì)算了其他深度上產(chǎn)生的庫(kù)侖破裂應(yīng)力變化，發(fā)現(xiàn)模式和量值變化不大。

從庫(kù)侖破裂應(yīng)力變化結(jié)果（圖8）來(lái)看：（1）同震產(chǎn)生的庫(kù)侖破裂應(yīng)力顯著增加斷裂主要包括龍門(mén)山斷裂南端和北端、甘孜—玉樹(shù)斷裂、鮮水河斷裂、白玉斷裂、東昆侖斷裂東段部分段落，這些部位庫(kù)侖破裂應(yīng)力增加值均超過(guò)0.01MPa，達(dá)到強(qiáng)震觸發(fā)的閾值（King et al，1994；Harris，Simgson，1998）。震后隨時(shí)間的變化，這些斷裂庫(kù)侖破裂應(yīng)力大部分在同震的基礎(chǔ)上有所增加，但增加量值不太顯著。汶川地震之后，這些斷裂相繼有5.0級(jí)以上強(qiáng)地震發(fā)生，其中包括在甘孜一玉樹(shù)斷裂發(fā)生的2010年玉樹(shù)7.1級(jí)、龍門(mén)山斷裂南段發(fā)生的2013年蘆山7.0級(jí)和鮮水河斷裂發(fā)生的2014年康定6.3級(jí)地震。（2）同震庫(kù)侖應(yīng)力有所增加，但增加不是很顯著的斷裂，主要包括大涼山地塊東邊界馬邊斷裂、大涼山地塊南邊界魯?shù)橐徽淹〝嗔?、金沙江斷裂、麗江小金河斷裂的南段、南汀河斷裂，?kù)侖應(yīng)力增加值小于0.01MPa。魯?shù)橐徽淹〝嗔褞?kù)侖破裂應(yīng)力震后隨時(shí)間變化逐漸增強(qiáng)，震后10a達(dá)到0.01MPa，該斷裂附近先后發(fā)生了幾次5.0級(jí)以上強(qiáng)震，其中包括2014年8月3日魯?shù)?.5級(jí)地震。麗江一小金河斷裂震后隨時(shí)間推移庫(kù)侖破裂應(yīng)力也逐漸增強(qiáng)，其南段震后10a庫(kù)侖應(yīng)力接近0.01MPa，在其附近發(fā)生了2012年寧蒗5.7級(jí)地震。金沙江斷裂震后隨時(shí)間變化庫(kù)侖破裂應(yīng)力有所增強(qiáng)，南汀河斷裂震后隨時(shí)間變化庫(kù)侖應(yīng)力有所增加，但并不顯著。（3）同震庫(kù)侖破裂應(yīng)力減少的斷裂包括龍門(mén)山斷裂中段發(fā)震斷裂、安寧河—?jiǎng)t木河—小江斷裂、麗江—小金河斷裂北段、紅河斷裂及東昆侖斷裂東段部分段落。除了麗江—小金河斷裂北段在震后隨時(shí)間推移，庫(kù)侖破裂應(yīng)力由同震的卸載轉(zhuǎn)為震后的加載外，其他斷裂震后隨時(shí)間變化情況和同震效應(yīng)一致，均存在不同程度的卸載增強(qiáng)的特征?？傊鸺罢鸷髷鄬訋?kù)侖破裂應(yīng)力反映的加卸載情況和強(qiáng)震活動(dòng)空間分布具有一定的相關(guān)性。需要說(shuō)明的是由于使用的模型及方法的不同，本文計(jì)算獲得的同震庫(kù)侖破裂應(yīng)力變化和萬(wàn)永革等（2009）采用位錯(cuò)方法得到的結(jié)果存在明顯的差異。

3 討論與結(jié)論

本文借助三維粘彈性有限元技術(shù)，以Wang等（2011）給出的大范圍汶川地震同震水平位移場(chǎng)作為約束，模擬了汶川地震對(duì)川滇地區(qū)的同震及震后形變場(chǎng)、應(yīng)力的時(shí)空演化特征，在模擬過(guò)程中，不但考慮了川滇地區(qū)地球介質(zhì)的橫向和縱向非均性，而且對(duì)同震上下盤(pán)斷層面靜態(tài)位移加載比例進(jìn)行了調(diào)整，盡可能使我們的模擬結(jié)果在中遠(yuǎn)場(chǎng)和實(shí)際觀測(cè)結(jié)果在形態(tài)上具有一致性。

模擬得到的川滇地區(qū)同震和震后形變場(chǎng)、應(yīng)力狀態(tài)反映的斷層活動(dòng)特征及庫(kù)侖破裂應(yīng)力變化反映的斷層面應(yīng)力加卸載情況，在空間上和強(qiáng)震活動(dòng)具有較好的相關(guān)性。對(duì)于理解和解釋汶川地震之后川滇地區(qū)出現(xiàn)的一些觀測(cè)異常具有一定的參考意義。但本文模擬結(jié)果相比實(shí)際觀測(cè)結(jié)果還存在許多不一致的地方，在模擬過(guò)程中仍然存在許多有待商榷和解決的問(wèn)題。主要包括以下幾個(gè)方面：（1）有限元邊界約束條件對(duì)結(jié)果有著重要的影響，位錯(cuò)模型認(rèn)為在離開(kāi)斷層無(wú)限遠(yuǎn)處和無(wú)限深處同震位移為0，即固定不動(dòng)。在有限元模型中，由于模型范圍總是有限的，我們無(wú)法做到這一點(diǎn)。本文采用大多數(shù)學(xué)者的做法，即模型邊界法向約束，切向自由移動(dòng)，這和實(shí)際情況存在明顯的差異，也可能是本文模擬結(jié)果相比觀測(cè)結(jié)果中遠(yuǎn)場(chǎng)同震形變相對(duì)較大的主要原因。因此探索更為有效、合理的邊界條件約束是在模擬過(guò)程需要進(jìn)一步解決的問(wèn)題：（2）本文簡(jiǎn)單地通過(guò)采用上下盤(pán)斷層面同震錯(cuò)動(dòng)位移比例，來(lái)擬合測(cè)同震形變場(chǎng)，僅僅采用一個(gè)比例系數(shù)，過(guò)于簡(jiǎn)單，也可能是擬合效果存在差異的原因之一：（3）本文模型沒(méi)有考慮重力、地殼背景構(gòu)造應(yīng)力場(chǎng)作用，重力、背景構(gòu)造應(yīng)力對(duì)同震及震后變形也存在一定的影響。另外，地殼介質(zhì)本身具有復(fù)雜性（地形、斷層的非連續(xù)性）等，這對(duì)模擬結(jié)果都會(huì)產(chǎn)生一定影響。在后續(xù)的數(shù)值模擬工作中需要進(jìn)一步探索和研究。

參考文獻(xiàn)：

畢研磊，張永志，曹海坤，等.2017.利用斷層自剖分技術(shù)反演汶川MW7.9地震滑動(dòng)分布[J].40（2）：211-215.

陳運(yùn)泰，許力生，張勇，等.2008.2008年5月12日汶川特大地震震源特性分析報(bào)告[EB/OL].（2017-09-20）[2008-05].ht-tp：//www.csi.ac.cn/sichuan/.

鄧起東，張培震，冉勇康，等.2002.中國(guó)活動(dòng)構(gòu)造基本特征[J].中國(guó)科學(xué)：地球科學(xué)，32（12）：1020-1030.

季靈運(yùn)，劉傳金，徐晶，等.2017.九寨溝MS7.0級(jí)地震的InSAR觀測(cè)及發(fā)震構(gòu)造分析[J].地球物理學(xué)報(bào)，60（10）：4069-4082.

季靈運(yùn)，劉立煒，郝明.2015.利用InSAR技術(shù)研究滇西南鎮(zhèn)康永德地區(qū)現(xiàn)今地殼形變特征[J].地震研究，38（1）：84-89.

焦佳爽，張永志，張凱南，等.2017.基于同震GPS數(shù)據(jù)的汶川大震同震應(yīng)變特征研究[J].地震工程學(xué)報(bào)，39（3）：227-233.

李寧，趙強(qiáng)，李金.2017.基于InSAR形變觀測(cè)反演2015年皮山MS6.5地震滑動(dòng)分布[J].地震，37（2）：67-76.

李強(qiáng)，江在森，武艷強(qiáng)，等.2013.汶川地震震后短期青藏高原東北緣地殼變形特征成因討論[J].地震，33（3）：124-132.

李玉江，陳連旺，楊樹(shù)新.2013.基于應(yīng)變能變化的蘆山強(qiáng)震同震效應(yīng)的數(shù)值模擬[J].震災(zāi)防御技術(shù)，8（4）：361-369.

邵志剛，傅容珊，薛霆虓，等.2008.昆侖山MS8.1級(jí)地震震后變形場(chǎng)數(shù)值模擬與成因機(jī)理探討[J].地球物理學(xué)報(bào)，51（3）：805-816.

沈正康，萬(wàn)永革，甘衛(wèi)軍，等.2003.東昆侖活動(dòng)斷裂帶大地震之間的黏彈性應(yīng)力觸發(fā)研究[J].地球物理學(xué)報(bào)，46（6）：786-795.

談洪波，申重陽(yáng)，李輝，等.2009.汶川大地震震后重力變化和形變的黏彈分層模擬[J].地震學(xué)報(bào)，31（5）：491-505.

萬(wàn)永革，沈正康，盛書(shū)中，等.2009.2008年汶川大地震對(duì)周?chē)鷶鄬拥挠绊慬J]地震學(xué)報(bào)，31（2）：129-139.

萬(wàn)永革，沈正康，曾躍華，等.2007.青藏高原東北部的庫(kù)侖應(yīng)力積累演化對(duì)大地震發(fā)生的影響[J].地震學(xué)報(bào)，29（2）115-129.

主椿鏞，楊文采，吳建平，等.2015.南北構(gòu)造帶巖石圈結(jié)構(gòu)與地震的研究[J].地球物理學(xué)報(bào)，58（11）：3867-3901.

王輝，曹建玲，張懷，等，2007.川滇地區(qū)下地殼流動(dòng)對(duì)上地殼運(yùn)動(dòng)變形影響的數(shù)值模擬[J].地震學(xué)報(bào)，29（6）：581-591.

吳建平，黃媛，張?zhí)熘校龋?009.汶川MS8.0級(jí)地震余震分布及周邊區(qū)域P波三維速度結(jié)構(gòu)研究[J].地球物理學(xué)報(bào)，52（2）：320-328.

吳建平，明躍紅，王椿鏞.2006.川滇地區(qū)速度結(jié)構(gòu)的區(qū)域地震波形反演研究[J].地球物理學(xué)報(bào)，49（5）：1369-1376.

張培震，鄧起東，張國(guó)民，等.2003.中國(guó)大陸的強(qiáng)震活動(dòng)與活動(dòng)地塊[J].中國(guó)科學(xué)：地球科學(xué)，33（增刊1）.12-20.

Fielding E J，Sladen A，Li Z，et aL 2013.Kipemalic fault slip evolutionsource models of he 2008 M7.9 Wenchuan.earthquake in Chinafrom SAIL interferometry，GPS and teleseismic analysis and implies-Lionsfor Lougmen Shan tectonics[J].Geophysical journal interns-tional，194（2）：1138-1166.

Funning G J，Parsons 13，Wright T J.2007.Fault slip in the 1997Manyi，Tibet earthquake from linear elastic modelling of InSAR displace-ments[J].Geophysical Journal International，169（3）：988-1008.

Ham's R A，Simpson R W.1998.Supp cession of large earquakes by stressshadow：A Conpaiison of Coulomb and rate-and-state failure[J].Journal of Geophysical Research solid Earth，103（B10）：24439-24451.

Jaeger J C，Cook N G W，Zimmerman R.1969.Fundamentals of rock me-chanics[M].Chapman and Hall'Landon.

Ji C，Hayes G.2008.NEIC，Preliminary result of theMay 12，2008MW7.9 eastern Sichuan，China earthquake[EB/OL].（2017-09-20）[2008-05].http；//earthquake.nags.gov/eqcenter/eqinthe-news/us2008ryan/.

Jonssons，Segall P，Pedersen R，et al.2003.Post-earthquake groundmovements correlated to pore-pressure transients[J].Nature，424（4945）：179-183：

Jouanne F，Awan A，Madji A，et al.2011.Postseismic deformation in Pa-kistan after the 8 October 2005 earthquake：evidence of afterslip a-long a flat north of the Balakot-Bagh thrust[J].Journal of Geo-physical Research：Solid Earth，116（B07401）：1-22.

King G C P，Stein Rs，Lin J.1994.Static stress changes and the trigge-ring of earthquakes[J].Bulletin of the Seismological Society of A-merica，84（3）：935-953.

Rice J R，Go J.1983.Earthquake after effects and triggered seismic phe-nomena[J].Pure and Applied Geophysics，121（2）：187-219.

Shen Z K，Sun J，Zhang P，et al.2009.Slip maxima at fault junctions andrupturing of barriers dining the 2008 Wenchuan earthquake[J].Na-ture geosoence，2009，2（10）：718-724.

Stein R S，King G C P，Lin J.1992.Change in failure stress on the south-ern San Andreas fault system caused by the 1992 magnitude=7.4Lenders earthquake[J].Science，258（5086）：1328-1332.

Wan Y，Shen Z K，Burgmann R，et at.2017.Fault geometry and slip dis-tribution of the 2008 MW7.9 Wenchuan，China earthquake，inferredfrom GPS and InSAR measurements[J].Geophysical Journal tnter-nationel，208（2）：748-766.

Wang，Xu C，Ge L.2007.Coseismic deformation and slip distribution ofthe 1997 MW7.5 Manyi，Tibet，earthquake from InSAR measure-ments[J].Journal of Geodynamics，44（3）：200-212.

Wang Q，Qiao X J，Len Q G，et al.2011.Rupture of deep faults in the2008 Wenchuan earthquake and uplift of the Longmen Shan[J].Na-Lure Ceascience，4（9）：634-640.