曾祥方，羅 艷，韓立波，石耀霖＊

1 中國科學院計算地球動力學重點實驗室，中國科學院大學，北京 100049

2 中國地震局地震預測研究所，北京 100036

3 中國地震局地球物理研究所，北京 100081

1 引 言

據中國地震臺網測定，北京時間2013 年4 月20日8時02分46秒四川省雅安市蘆山縣（30．3°N，103．0°E）發(fā)生MS7．0 地震（圖1，以下簡稱蘆山地震）．震源機制解研究是認識地震發(fā)震斷層的重要手段，也是理解深部構造應力和地震發(fā)震機理的重要依據．地震發(fā)生后國內外研究機構，通過各種地震學方法，初步研究了震源機制解、破裂過程等震源參數［1－2］，結果表明蘆山地震是發(fā)生在龍門山斷裂帶南段上的一次逆沖地震，但是不同機構所給出的地震斷層面傾角存在一定差異．在震源參數研究中，遠震數據對逆沖地震的傾角具有較好的約束，而近震可以較好地確定斷層走向．本文利用近遠震數字波形記錄，分別利用P 波初動和波形聯(lián)合反演方法，探討本次地震的震源機制解，通過分析不同參數對斷層面傾角的影響，并結合短期余震機制解，確認了蘆山地震是一個高角度逆沖地震．

2 近遠震P波初動聯(lián)合反演

利用P波初動反演震源機制解的基本原理是通過對比不同機制解預測的各臺站P 波初動符號與實際觀測結果的差異，求得矛盾符號最小的一個或一組解作為地震機制解結果［4］．該方法較為簡單，但是結果可靠性較為依賴所用臺站分布和離源矢量誤差．臺站分布則決定了離源矢量在震源球上的采樣，不同震中距的采樣空間不同，遠震臺站的離源矢量主要集中在震源球中心區(qū)域，近震臺站除了少數震中距較小的臺站可以采樣中心區(qū)域外，大部分近震臺站集中在外圍區(qū)域．蘆山地震的初步研究結果表明其為一高角度逆沖地震，近震臺站初動符號對斷層的傾角和滑動角的約束不夠理想．而遠震臺站可以密集采樣震源球中心位置，對近震臺站數據形成了有效補充．我們從IRIS（Incorporated Research Institutions for Seismology）數據中心下載了蘆山地震在全球地震臺網（Global Seismographic Network，GSN）的寬頻帶數字波形記錄，選取震中距在35°～85°，清晰直達P波信號的36個臺站的記錄（圖1b）讀取遠震P波初動符號，而近震初動符號則在國家數字測震臺網數據備份中心［5］提供的波形數據中，選取了200km 以內初動最清楚的18 個臺站，進行P波初動拾?。?/p>

圖1 （a）近震臺站分布圖；（b）遠震臺站分布圖紅色五角星代表蘆山地震震中，藍色三角代表臺站，黑色實線代表活動斷層［3］．Fig．1 （a）Regional station distribution map；（b）Teleseismic station distribution mapStar denotes epicenter while triangles indicate stations．black lines denotes fault．

離源矢量計算中的誤差主要來自震源深度誤差和地震波速度模型誤差．本文采用中國地震臺網快報目錄給出的蘆山地震震中位置，深度為13km，將深度誤差設為1km．近震離源角計算中所用速度模型為前人汶川地震余震研究中所用速度模型［6］基礎上修改得到（表1，WC），而遠震離源角計算采用IASP91模型．Hardebeck 和Shearer提出速度模型和震源深度的誤差可以根據離源角、方位角、深度誤差，進行隨機擾動，生成新的離源矢量，然后進行網格搜索并得到最可信解［7］．由于地殼中地震波傳播較為復雜，將近震出射角的誤差設為5°，而遠震地震波出射角誤差一般較小，設為2°，方位角誤差均為1°．震源深度的誤差為1km，進行50次嘗試后，得到斷層面參數為212°/44°/92°（圖2），該解給出的理論初動符號與實際觀測結果相符，其走向與龍門山斷裂帶走向接近．

圖2 P波初動反演結果＋表示初動為正，○表示初動為負，紅色區(qū)域為壓縮區(qū)域，理論初動為正，白色區(qū)域為拉張區(qū)域，理論初動為負．Fig．2 P wave first－motion focal mechanismPlus denotes positive first motion，while circle means negative one．

表1 蘆山地震震源區(qū)地殼地震波速度模型Table 1 Seismic velocity model of epicenter regionof the Lushan earthquake

3 近遠震波形聯(lián)合反演

P波初動反演可以快速給出地震機制解，但是缺少其他震源參數，而利用波形反演方法，可以充分利用地震記錄中各種震相到時、振幅、相位等信息，對震源機制提供有效約束，不僅可以提供斷層面參數，還可以提供標量地震矩、矩心深度、破裂持續(xù)時間等震源參數［8－9］．傳統(tǒng)的地震波形反演方法一般只單獨采用近震或遠震波形數據，陳偉文等［9］發(fā)展了改進 的“剪切—粘貼”方法（Cut－And－Paste joint，CAPjoint），利用近遠震波形數據進行聯(lián)合反演，表明近遠震數據可以相互補充，更好地約束地震震源參數，尤其是對于逆沖地震，遠震對傾角的約束具有一定優(yōu)勢［9］，因此本文采用該方法進行波形反演．

遠震波形數據可包括豎直向P 波和切向SH波，其格林函數計算過程分為震源區(qū)地殼、地幔、臺站區(qū)地殼三步．由于不同區(qū)域的地殼地震波速度結構差異，需要選取特定震源區(qū)地殼模型，利用傳播矩陣方法計算震源區(qū)和臺站區(qū)地殼部分傳播效應，地幔部分則采用初步地球參考模型（Primarily Reference Earth Model，PREM），計算幾何擴散效應和非彈性衰減效應，最后將三者卷積得到格林函數［10］．近震波形數據包括體波和面波部分，格林函數計算采用頻率波數域積分方法（F－K）［11］．反演方法則與“剪切—粘貼”（Cut－And－Paste，CAP）方法［8］相似，搜索不同震源深度的最佳雙力偶解，然后對比不同深度反演誤差，得到最佳震源深度和雙力偶解．遠震和近震震源區(qū)地殼地震波速度結構模型都采用初動反演所用模型．

前期研究表明，由于蘆山地震為逆沖型地震，其遠震SH 波部分信噪比較P波部分低［12］，因此未參與反演．遠震P 波經去除傾斜，去除儀器響應等預處理后，根據IASP91計算的理論到時進行截取，窗口長度為60s，濾波頻段選為0．05～0．1 Hz．四川盆地和青藏高原東緣速度模型差異較大［13－16］，所以近震波形中主要采用青藏高原東緣臺站，且利用傳播效應較為簡單的徑向和豎直向長周期S 波和面波記錄，為了減少速度模型橫向變換的影響，主要選取了震中西部和南部較遠距離臺站．原始數據經預處理后，根據震中位置，將水平記錄旋轉徑向和切向，并根據S波理論到時進行截取，窗口長度為60 s，濾波至0．05～0．1Hz．遠震和近震數據間權重根據陳偉文等［7］提出的方法進行選取．

我們計算了2～20km 間10個深度下的各臺站格林函數，破裂時間設為10s，首先在各深度對斷層走向、傾角、滑動角以10°間隔進行搜索，得到的最佳矩心深度為12km，最佳雙力偶機制解為：節(jié)面1：211°/47°/90°；節(jié)面2：30°/43°/89°．根據中國地震臺網快報目錄（data．earthquake．cn）給出的短期余震分布，我們認為節(jié)面1為可能的發(fā)震斷層．在此基礎上，將搜索間隔縮小為2．5°，得到解為：節(jié)面1：216°/47°/93°；節(jié) 面2：31°/43°/86°，矩 震 級 為6．7，其理論波形和實際記錄對比見圖3．由波形對比圖可見波形擬合較好，遠震P波部分互相關系數均不小于0．8，近震面波窗口的徑向和豎直向波形互相關系數在0．8左右．

圖3 蘆山地震CAPjoint反演實際波形與理論波形對比紅線為理論波形，黑線為實際波形，波形下方數字分別為互相關系數百分比（%）和對應時間偏移（s），臺站名上方數字為方位角（°）和震中距（km）．Fig．3 CAPjoint modeling for Lushan earthquakeRed line shows theoretical seismogram while black line shows data．Numbers below waveform are cross－correlation（%）in percent（%）and time shift（s）．Numbers above station name are azimuth（°）and epicenteral distance（km）

4 討論與結論

本文根據P 波初動和波形反演得到的斷層面參數為212°/44°/92°和212°/47°/93°，基本符合震源區(qū)構造應力場［17－19］，但是斷層面三個參數中，走向和滑動角與其他機構結果相近，但傾角較大．對于逆沖斷層而言，傾角越大，上盤向上錯動所需能量越大，以發(fā)生在同一斷裂帶的汶川地震為例，GlobalCMT（www．globalcmt．org）給出的斷層傾角為35°，而蘆山地震傾角明顯較大．為了研究斷層面傾角的不確定性，需要分離斷層面走向、滑動角、深度等模型參數影響．采用波形反演得到斷層面的走向（216°）和滑動角（93°），首先計算了12km 深度下不同傾角的波形反演誤差，結果見圖4a，隨著傾角遠離最佳解（47°），誤差迅速增加，表明這一傾角是可靠的，而固定斷層面參數得到的不同深度反演誤差也表明這一深度結果較為可靠（圖4b）．

波形反演過程中地震波速度模型對反演結果也有一定影響，雖然本文未選取四川盆地內部臺站，但是為了考察不同模型對斷層傾角的影響，選取了三個較為典型的速度模型進行分析，分別是震中處由CrustModel2．0 模型基礎上的修改模型（表1 中CM2）［12］、較為極端的盆地模型［19－20］（表1中SC）和低速體模型［21］（表1中LVZ）．在已有機制解的斷層走向（200°～230°）和滑動角（80°～110°）附近搜索最佳傾角，結果見圖5．三個速度模型下，大部分機制解的傾角較高（＞40°），但是低速層模型由于低速層引起較強反射，增加了波形高頻成分，導致誤差分布跳躍較大，且震源深度出現了10km 和16km兩個局部最優(yōu)解．

采用同樣的模型，利用200km 內臺站的體波和300km 內臺站面波，利用CAP 方法對蘆山地震震后4個小時內的較強余震機制解進行了初步研究，表2給出了雙力偶解中沿龍門山斷裂帶走向的節(jié)面參數、震源深度和震級．大部分余震機制解與主震相近，其中3 號地震較為特殊，但另一個節(jié)面解（36°/82°/73°）也為高角度逆沖，具體發(fā)震斷層面仍有待進一步分析．同時這一區(qū)域的地質考察也表明當地存在較高傾角斷裂［22］．

表2 蘆山地震短期較強余震震源機制解Table 2 Focal mechanisms of short term aftershocks of Lushan earthquake

如此高角度的逆斷層錯動，只有在極大的差應力和長期應變能積累后才能發(fā)生［23］．龍門山斷裂帶位于青藏高原東南緣，由于青藏高原隆升引起的地殼物質東流，在此遇到四川盆地阻擋，因此應變能長期積累，定量計算模擬表明，整個龍門山斷裂帶上地殼底部都是應力高增長率區(qū)［24］．在汶川地震破裂過程中，主要破裂集中在龍門斷裂帶北段，而南部基本沒有發(fā)生錯動，余震活動也較小，因此其上積累的應變能并沒有得到釋放．汶川地震對周圍斷層的影響較大，根據完全彈性和黏彈性情況下計算得到的庫侖應力變化，蘆山地震震源區(qū)的庫侖應力增加量約為0．01MPa［25－27］．近期在龍門山斷裂帶南段進行的原位地應力測量結果也表明這一區(qū)域的最大水平主應力已達斷層活動應力臨界下限值，斷裂活動進入臨界狀態(tài)［16］，蘆山地震正是發(fā)震斷層面上剪切應力超過閾值引起的斷層錯動．

本文綜合利用近遠震地震數據，通過P 波初動反演和波形反演方法，確定了蘆山地震的震源機制解，并結合短期余震的機制解，認為發(fā)震斷層為一高角度逆沖斷層．由于在反演過程中采用的點源近似，分層均勻地震波速度模型等與真實情況存在一定差異．地震機制解雖然有助于探討發(fā)震斷層，但具體發(fā)震斷層，仍需要通過震后地質考察、高精度主震和余震定位等方法得到．如此高角度逆沖地震的力學機制則需要通過多種方法得到震源區(qū)構造應力場、介質模型、斷層形態(tài)和摩擦系數等參數進行定量分析得到．

致 謝 中國地震局地球物理研究所國家數字測震臺網數據備份中心為本研究提供地震波形數據，GSN 臺網數據由IRIS DMC 提供，本文圖件采用Generic Mapping Tools繪制，在此一并表示感謝．

（References）

［1］ 張勇，許力生，陳運泰．蘆山4．20地震破裂過程及其致災特征初步分析．地球物理學報，2013，56（4）：1408－1411，doi：10．6038/cjg20130435．

Zhang Y，Xu L S，Chen Y T．Rupture process of the Lushan 4．20earthquake and preliminary analysis on the disastercausing mechanism．ChineseJ．Geophys．（in Chinese），2013，56（4）：1408－1411，doi：10．6038/cjg20130435．

［2］ 王衛(wèi)民，郝金來，姚振興．2013年4月20日四川蘆山地震震源破裂過程反演初步結果．地球物理學報，2013，56（4）：1412－1417，doi：10．6038/cjg20130436．

Wang W M，Hao J L，Yao Z X．Preliminary result for rupture process of Apr．20，2013，Lushan Earthquake，Sichuan，China．ChineseJ．Geophys．（in Chinese），2013，56（4）：1412－1417，doi：10．6038/cjg20130436．

［3］ 鄧起東，張培震，冉勇康等．中國活動構造基本特征．中國科學（D輯），2002，32（12）：1020－1030．

Dong Q D，Zhang P Z，Ran Y K，et al．Basic feature of China active structures．Sci．ChinaSer．D－EarthSci．（in Chinese），2002，32（12）：1020－1030．

［4］ 許忠淮，汪素云，黃雨蕊等．由大量地震資料推斷的我國大陸構造應力場．地球物理學報，1989，32（6）：636－647．

Xu Z H，Wang S Y，Huang Y R，et al．The tectonic stress field of Chinese continent deduced from a great number of earthquakes．ChineseJ．Geophys．（in Chinese），1989，32（6）：636－647．

［5］ 鄭秀芬，歐陽飚，張東寧等．“國家數字測震臺網數據備份中心”技術系統(tǒng)建設及其對汶川大地震研究的數據支撐．地球物理學報，2009，52（5）：1412－1417．

Zheng X F，Ouyyang B，Zhang D N，et al．Technical systerm construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake．ChineseJ．Geophysi．（in Chinese），2009，52（5）：1412－1417．

［6］ 鄭勇，馬宏生，呂堅等．汶川地震強余震（Ms≥5．6）的震源機制解及其與發(fā)震構造的關系．中國科學（D 輯），2009，39（4）：413－426．

Zheng Y，Ma H W，Lu J，et al．Source mechanism of strong aftershock for seismotectonic．Sci．ChinaSer．D－EarthSci．（in Chinese），2009，39（4）：413－426．

［7］ Hardebeck J L，Shearer P M．A new method for determining first－motion focal mechanisms．BulletinoftheSeismological SocietyofAmerica，2002，92（6）：2264－2276．

［8］ Zhu L，Helmberger D V．Advancement in source estimation techniques using broadband regional seismograms．Bulletin oftheSeismologicalSocietyofAmerica，1996，86（5）：1634－1641．

［9］ 陳偉文，倪四道，汪貞杰等．2010年高雄地震震源參數的近遠震波形聯(lián)合反演．地球物理學報，2012，55（7）：2319－2328．

Chen W W，Ni S D，Wang Z J，et al．Joint inversion with both local and teleseismic waveforms for source parameters of the 2010 Kaohsiung earthquake．ChineseJ．Geophys．（in Chinese），2012，55（7）：2319－2328．

［10］ Kikuchi M，Kanamori H．Inversion of complex body waves．BulletinoftheSeismologicalSocietyofAmerica，1982，72（2）：491－506．

［11］ Zhu L，Rivera L A．A note on the dynamic and static displacements from a point source in multilayered media．GeophysicalJournalInternational，2002，148（3）：619－627．

［12］ 曾祥方，韓立波，石耀霖．四川蘆山Ms7．0地震震源機制解初步研究．科技導報．2013，31（14）：5－8．

Zeng X F，Han L B，Shi Y L．Focal mechanism of LushanMs7．0Earthquake，April 20，2013．ScienceandTechnology Servies．2013，31（14）：5－8．

［13］ 劉啟元，李昱，陳九輝等．汶川Ms8．0地震：地殼上地幔S波速度結構的初步研究．地球物理學報，2009．52（2），309－319．

Liu Q Y，Li Y，Chen J H，et al． WenchuanMS8．0 earthquake：preliminary study of the S－wave velocity structure of the crust and upper mantle．ChineseJ．Geophys．（in Chinese），2009，52（2）：309－319．

［14］ 吳建平，黃媛，張?zhí)熘械龋氪∕S8．0級地震余震分布及周邊區(qū)域P 波三維速度結構研究．地球物理學報，2009，52（2）：320－328．

Wu J P，Huang Y，Zhang T Z，et al．Aftershock distribution of theMs8．0 Wenchuan earthquake and three dimensional Pwave velocity structure in and around source region．Chinese J．Geophys．（in Chinese），2009，52（2）：320－328．

［15］ 雷建設，趙大鵬，蘇金蓉等．龍門山斷裂帶地殼精細結構與汶川地震發(fā)震機理．地球物理學報，2009，52（2）：339－345．

Lei J S，Zhao D P，Su J R，et al．Fine seismic structure under the Longmengshan fault zone and the mechanism of large Wenchuan earthquake．ChineseJ．Geophys．（in Chinese），2009，52（2）：339－345．

［16］ Li Z W，Ni S D，Hao T Y，et al．Uppermost mantle structure of the eastern margin of the Tibetan plateau from interstation Pn traveltime difference tomography．Earthand PlanetaryScienceLetters，2012．335－336：195－205．

［17］ 鐘繼茂，程萬正．由多個地震震源機制解求川滇地區(qū)平均應力場方向．地震學報，2006，28（4）：337－346．

Zhong J M，Cheng W Z．Determination of directions of the mean stress field in Sichuan－Yunnan region from a number of focal mechanism solutions．ActaSeismologicaSinica（in Chinese），2006，28（4）：337－346．

［18］ 秦向輝，陳群策，譚成軒等．龍門山斷裂帶西南段現今地應力狀態(tài)與地震危險性分析．巖石力學與工程學報，2013，32（增1）：2870－2876．

Qin X H，Chen Q C，Tan C X，et al．Analysis of current geostress state and seismic risk in southwest segment of Longmenshan fracture belt．ChineseJournalofRockMechanics andEngineering（in Chinese），2013，32（S1）：2870－2876．

［19］ 萬柯松，倪四道，趙里．龍門山周邊地區(qū)地震波速度結構的有限頻率反演．國際地震動態(tài)，2010，6：16－16．

Wan K S，Ni S D，Zhao L．Finite frequency waveform inversion for velocity structure in Longmengshan area．Recent DevelopmentsinWorldSeismology（in Chinese），2010，6：16－16

［20］ 謝軍，倪四道，曾祥方．四川盆地中部淺層地殼一維剪切波速度結構初步研究．四川地震，2012，143（2）：20－24．

Xie J，Ni S D，Zeng X F．1Dshear wave velocity structure of the shallow upper crust in central Sichuan Basin．Earthquake ResearchinSichuan（in Chinese），2012，143（2）：20－24．

［21］ 王椿鏞，吳建平，樓海等．川西藏東地區(qū)的地殼P 波速度結構．中國科學（D 輯），2003，33（增刊）：181－189．

Wang C Y，Wu J P，Lou H，et al．P－wave velocity structure in western Sichuan and eastern Tibetan region．Sci．China Ser．D－EarthSci．，2003，46（Suppl．）：254－265．

［22］ 楊曉平，蔣溥，宋方敏等．龍門山斷裂帶南段錯斷晚更新世以來地層的證據．地震地質，1999，21（4）：341－345．

Yang X P，Jiang P，Song F M，et al．The evidence of the south Longmenshan fault zones cutting late Quaternary stratum．SeismologyandGeology（in Chinese），1999，21（4）：341－345．

［23］ 周永勝，何昌榮．汶川地震區(qū)的流變結構與發(fā)震高角度逆斷層滑動的力學條件．地球物理學報，2009，52（2）：474－484．

Zhou Y C，He C R．The rhelogical structure of crust and mechanics of high－angle reverse fault slip for WenchuanMs8．0 earthquake．ChineseJ．Geophys．（in Chinese），2009，52（2）：474－484．

［24］ 柳暢，朱伯靖，石耀霖．粘彈性數值模擬龍門山斷裂帶應力積累及大震復發(fā)周期．地質學報，2012，86（1）：157－169．

Liu C，Zhu B J，Shi Y L．Stress accumulation of the Longmenshan Fault and recurrence interval of Wenchuan Earthquake based on viscoelasticity simulation．ActaGeologica Sinica（in Chinese），2012，86（1）：157－169．

［25］ 單斌，熊熊，鄭勇等．2008年5月12日Mw7．9汶川地震導致的周邊斷層應力變化．中國科學（D 輯），2009，39（5）：527－545．

Shan B，Xiong X，Zheng Y，et al．Stress changes on major faults caused byMw7．9 Wenchuan earthquake，May 12，2008．，Sci．ChinaSer．D－EarthSci．，2009，52（5）：593－601．

［26］ 邵志剛，周龍泉，蔣長勝等．2008年汶川MS8．0地震對周邊斷層地震活動的影響．地球物理學報，2010，53（8）：1784－1795．

Shao Z G，Zhou L Q，Jiang C S，et al．The impact of WenchuanMs8．0earthquake on the seismicity of surrounding faults．ChineseJ．Geophysi．（in Chinese），2010，53（8）：1784－1795．

［27］ 石耀霖，曹建玲．庫侖應力計算及應用過程中若干問題的討論—以汶川地震為例．地球物理學報，2010，53（1）：102－110．

Shi Y L，Cao J L．Some aspects in static stress change calculation－case study on Wenchuan earthquake．ChineseJ．Geophys．（in Chinese），2010，53（1）：102－110．