張勇，陳運(yùn)泰，，許力生，魏星，金明培，張森

1 北京大學(xué)地球與空間科學(xué)學(xué)院，北京 100871

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

3 中國地震臺網(wǎng)中心，北京 100045

4 云南省地震局，昆明 650041

1 引言

2014年8月3日發(fā)生在云南魯?shù)榈腗W6.1（MS6.5）地震是一次中等強(qiáng)度的地震事件，但余震分布卻相當(dāng)復(fù)雜.根據(jù)最新的余震精定位結(jié)果，魯?shù)榈卣鸬挠嗾鸱植荚贜NW—SSE（以下簡稱近南北向）和ENE—WSW（以下簡稱近東西向）兩個方向上都有一定尺度的延伸（王未來等，2014；張廣偉等，2014），表明此次地震的發(fā)震構(gòu)造和破裂過程可能具有一定的復(fù)雜性.這一現(xiàn)象引起了地震學(xué)界的廣泛關(guān)注（徐錫偉等，2014；張勇等，2014；劉成利等，2014；張振國等，2014；王未來等，2014；張廣偉等，2014）.如何判別這次地震的發(fā)震斷層，成為了相關(guān)研究工作的一大熱點(diǎn).目前，烈度分布和破裂過程反演等研究都更傾向于近南北向的小河—包谷垴斷裂是此次地震的發(fā)震構(gòu)造（徐錫偉等，2014；張勇等，2014；劉成利等，2014）.而視震源時間函數(shù)分析表明，魯?shù)榈卣鹂赡苁枪曹棓鄬酉群笃屏训慕Y(jié)果（許力生等，2014）.另外，張廣偉等（2014）根據(jù)余震分布和震源機(jī)制解結(jié)果也提出可能存在共軛破裂的情況.為更定量地研究兩條斷層的破裂屬性，有必要在破裂過程反演中，進(jìn)一步考慮共軛斷層破裂的可能性.

實(shí)際上，主震破裂和余震呈現(xiàn)出復(fù)雜分布的地震并不鮮見.比如2011年盈江地震（Lei et al.，2012）和2012年蘇門答臘北部海域MW8.6地震.后者的破裂過程中，由于多達(dá)三條共軛斷層先后發(fā)生了破裂（Meng et al.，2012；H.Zhang et al.，2012；Yue et al.，2012），余震分布呈現(xiàn)出比較散亂的特征.然而，主震破裂與余震分布之間也不總是嚴(yán)格地相互對應(yīng).發(fā)生破裂的區(qū)域一般伴隨著余震分布，但有余震分布的區(qū)域不一定有破裂發(fā)生.一個典型的例子是2010年青海玉樹地震.該地震是一次以單側(cè)破裂為主的地震事件，破裂主要朝東南方向傳播（Zhang et al.，2010），但在主要破裂區(qū)域的另一端，即震中西北約40km處，也存在著比較集中的余震分布（Wang et al.，2013）.對于2014年魯?shù)榈卣?，余震在近東西向和近南北向都呈現(xiàn)出明顯的條帶狀分布，是兩條斷層都發(fā)生了破裂，還是只有一條斷層發(fā)生破裂，是震源研究需要關(guān)注和解決的問題.

本文搜集了距震中250km范圍內(nèi)強(qiáng)震與近震寬頻帶的全波形數(shù)據(jù)，以及遠(yuǎn)震寬頻帶體波數(shù)據(jù)，分別基于單一斷層和兩條交叉的共軛斷層模型，對魯?shù)榈卣鸬钠屏堰^程進(jìn)行了反演.通過比較與分析反演結(jié)果與反演得到的波形殘差，確定魯?shù)榈卣鹗窃诒蔽飨蛑鲏簯?yīng)力方向與北東向主張應(yīng)力方向的統(tǒng)一應(yīng)力場下發(fā)生的兩條共軛的斷層先后破裂的一次復(fù)雜地震事件.最后，我們根據(jù)得到的震源破裂模型，探討了共軛斷層破裂與形成復(fù)雜余震分布的原因.

2 數(shù)據(jù)

本文選用了13個強(qiáng)震臺站（圖1青色三角形）和較近距離上11個未限幅的寬頻帶臺站（圖1黃色三角形）的三分量全波形記錄，作為參與反演的近震數(shù)據(jù).在分別對強(qiáng)震數(shù)據(jù)的加速度記錄積分兩次與對寬頻帶數(shù)據(jù)的速度記錄積分一次后，得到相應(yīng)的位移記錄.采用張勇等（2014）使用的遠(yuǎn)震數(shù)據(jù)，經(jīng)過進(jìn)一步篩選，去掉了信噪比較低的2個P波臺站和1個SH波臺站，得到了26條垂直向P波波形和31條SH波波形（圖2），作為參與反演的遠(yuǎn)震數(shù)據(jù).與近震數(shù)據(jù)一樣，將遠(yuǎn)震數(shù)據(jù)積分得到位移記錄.為去掉近震數(shù)據(jù)存在的低頻基線漂移，并且考慮到寬頻帶儀器的有效頻段響應(yīng)范圍，我們對近震和遠(yuǎn)震數(shù)據(jù)都進(jìn)行了0.02～0.5Hz的帶通濾波，保證在去除近震數(shù)據(jù)基線漂移的同時，得到的數(shù)據(jù)信息位于寬頻帶儀器有效頻段響應(yīng)范圍內(nèi).

本文采用 Wang（1999）的方法計算格林函數(shù).魯?shù)榈卣鹫鹬兴诘拇ǖ岬貐^(qū)是地震多發(fā)區(qū)域，其地下速度結(jié)構(gòu)得到了廣泛的關(guān)注和比較系統(tǒng)的研究（熊紹柏等，1993；王椿鏞等，2002；Wang et al.，2003）.在反演中，我們選用王未來等（2014）給出的區(qū)域速度結(jié)構(gòu)，作為近震數(shù)據(jù)反演的介質(zhì)模型；遠(yuǎn)震數(shù)據(jù)反演采用了AK135全球速度結(jié)構(gòu)模型（Kennett et al.，1995）.與觀測記錄一樣，對格林函數(shù)位移波形記錄也進(jìn)行了相同的帶通濾波處理.

3 反演

3.1 反演方法與參數(shù)

根據(jù)王未來等（2014）對魯?shù)榈卣鹦蛄械亩ㄎ唤Y(jié)果，以及張廣偉等（2014）反演區(qū)域地震波資料得到的震源機(jī)制解，我們給定了兩條共軛斷層的位置、走向、傾角等參數(shù).其中，近南北向斷層為：走向162°／傾角70°／滑動角－14°，近東西向斷層為：走向257°／傾角77°／滑動角－159°；主壓應(yīng)力軸（P軸）方位角121°／俯角24°；主張應(yīng)力軸（T軸）方位角28°／俯角5°；中間主應(yīng)力軸（B軸）方位角288°／俯角66°.主壓應(yīng)力軸方向與魯?shù)橹苓厖^(qū)域的水平運(yùn)動方向一致，但斷層面走向與昭通斷裂的走向有所差別（聞學(xué)澤等，2013），表明魯?shù)榈卣鸬倪@一特殊的震源機(jī)制可能是在區(qū)域應(yīng)力背景下，因昭通斷裂的走向在魯?shù)楦浇l(fā)生局部改變而形成的.另一方面，由這兩個斷層面參數(shù)確定的斷層位置與余震分布（王未來等，2014；張廣偉等，2014）相當(dāng)一致，表明震源機(jī)制結(jié)果是可靠的.與作者之前的工作一樣（張勇等，2014），兩個斷層面的尺度都給定得足夠大，即長42km、寬20km，并將其劃分為21×10＝210個2km×2km大小的子斷層.初始破裂點(diǎn)均位于沿走向方向第11個、沿傾向方向向下第6個子斷層處，在兩個斷層面上的深度分別約為10km和11km.這與中國地震臺網(wǎng)中心得到的結(jié)果（12km）相仿，但略淺于精定位得到的15km和13km（王未來等，2014；張廣偉等，2014）.

本文采用Y.Zhang等（2012）的線性反演方法進(jìn)行反演.在這一方法中，無需預(yù)先假定子斷層震源時間函數(shù)的形式.相反地，子斷層震源時間函數(shù)是作為待定參數(shù)由波形反演自動給出.因此，可根據(jù)反演結(jié)果探討斷層兩盤相對錯動過程的時間復(fù)雜性.同時，作為約束條件，這一方法需要給定破裂速度和破裂持續(xù)時間上限.鑒于此次地震規(guī)模較小，震源尺度不大，我們給定破裂速度上限為3km·s－1，即不考慮超剪切破裂的情況，并限定最大破裂持續(xù)時間不超過10s.在這兩個參數(shù)相對寬松的約束下，斷層面上距震源（初始破裂點(diǎn)）dkm處破裂速度大于d／（d／3＋10）km·s－1小于3km·s－1的滑動信息都可以被順利提取.同時，這一相對比較寬松的破裂速度范圍可在一定程度上減弱震源位置（包括震源深度）的不確定性對破裂分布求解的干擾.

在破裂過程反演中，本文考慮了三種可能的斷層破裂模式：①破裂只發(fā)生在近南北向斷層上；②破裂只發(fā)生在近東西向斷層上；③破裂可同時在兩條共軛斷層上自由傳播.根據(jù)這三種斷層破裂模式進(jìn)行反演，得到不同的破裂模型，結(jié)合反演得到的波形擬合情況，可以推測真實(shí)的發(fā)震構(gòu)造和破裂過程是否具備多斷層特征.

圖2 震中與遠(yuǎn)震資料臺站位置（a）震中與P波臺站位置分布；（b）震中與SH波臺站位置分布.Fig.2 The epicenter and teleseismic stations（a）Epicenter and seismic stations of P waves.（b）Epicenter and seismic stations of SH waves.

3.2 近震數(shù)據(jù)反演

首先采用近震數(shù)據(jù)，分別針對以上三種情形進(jìn)行反演.反演結(jié)果如圖3所示.首先，從反演得到的滑動分布和震源時間函數(shù)看，共軛斷層破裂模型反演結(jié)果都更接近于近南北向斷層結(jié)果.表明在完全自由的條件下，破裂更傾向于發(fā)生在近南北向斷層上，意味著近南北向斷層應(yīng)當(dāng)是此次地震的主要發(fā)震斷層.其次，采用近南北向和近東西向斷層進(jìn)行反演得到的波形殘差分別為0.60和0.64，表明如果只考慮單一斷層破裂的情形，則近南北向斷層更有可能是真實(shí)的發(fā)震斷層，這也與現(xiàn)有的研究結(jié)果相符（徐錫偉等，2014；許力生等，2014；張勇等，2014；劉成利等，2014）.相比之下，共軛斷層模型波形殘差為0.57，意味著除近南北向斷層外，近東西向斷層的破裂對波形擬合也存在一定貢獻(xiàn).在共軛斷層模型反演結(jié)果中，近南北向斷層釋放了62%的地震矩，其對波形擬合的貢獻(xiàn)約60%；近東西向斷層釋放了38%的地震矩，其對波形擬合的貢獻(xiàn)約為40%.可見，盡管近東西向斷層上發(fā)生的破裂規(guī)模不占優(yōu)勢，但單位地震矩對資料擬合的貢獻(xiàn)程度還略大于近南北向斷層，表明它在地震破裂過程中可能也發(fā)生了破裂.

需要說明的是，以上結(jié)果（圖3）與作者之前采用遠(yuǎn)震資料反演得到的結(jié)果（張勇等，2014）在主要特征上一致：近南北向斷層反演結(jié)果中，破裂主要朝震中以南傳播，而近東西向斷層面的反演結(jié)果則顯示，震中以東的滑動分布更占優(yōu)勢.不同的是，在遠(yuǎn)震資料反演結(jié)果中，由于地震資料的信噪比較差等因素，只能根據(jù)得到的滑動分布與余震分布的一致性來推斷真實(shí)發(fā)震斷層；但在本文開展的近震資料反演中，由于數(shù)據(jù)相對更高的空間分辨能力，我們可以通過資料殘差來確定發(fā)震斷層.

3.3 聯(lián)合反演

圖3 近震數(shù)據(jù)反演得到的滑動量分布與震源時間函數(shù)從左至右分別為：近南北向斷層反演結(jié)果，近東西向斷層反演結(jié)果，以及兩條共軛斷層反演結(jié)果.每個反演結(jié)果中，從上到下依次為斷層面滑動量分布在地表的投影，震源時間函數(shù)（插圖），以及斷層面上的靜態(tài)滑動量分布.地圖中的紅線表示斷層（聞學(xué)澤等，2013）.Fig.3 Inversion results of local seismic wavesFrom left to right are the inversion results of NNW－SSE fault，ENE－WSW fault，and two conjugated fault.From top to bottom for each panel of inversion results are the surface projection of fault slip distribution，source time function （STF，see subset），and fault slip distribution.Red lines on the upper maps denote the faults（Wen et al.，2013）.

在近震數(shù)據(jù)反演基礎(chǔ)上，我們增加了遠(yuǎn)震地震數(shù)據(jù)，分別對單一斷層和共軛斷層模型進(jìn)行了聯(lián)合反演.由于近震與遠(yuǎn)震數(shù)據(jù)的聯(lián)合使用提高了對斷層的離源角覆蓋范圍，理論上可以增強(qiáng)反演在深度方向上的空間分辨能力（Zhang et al.，2014）.另外，一些強(qiáng)震數(shù)據(jù)位于土層上，波形記錄受臺站下方局部效應(yīng)的影響，其幅度有可能被放大.在格林函數(shù)計算中，這種非常局部的臺站效應(yīng)很難予以考慮.因此，近震數(shù)據(jù)和遠(yuǎn)震數(shù)據(jù)的聯(lián)合，還可更好地約束矩震級.在聯(lián)合反演中，兩種數(shù)據(jù)被賦予了相同的權(quán)重.反演結(jié)果如圖4所示.在主要破裂特征上，聯(lián)合反演結(jié)果與近震數(shù)據(jù)單獨(dú)反演結(jié)果基本相似.只是在矩震級和最大滑動量上，聯(lián)合反演結(jié)果稍小于近震數(shù)據(jù)反演結(jié)果.上文的近震資料反演與分析表明，近東西向斷層也可能發(fā)生了破裂，因此我們傾向于基于共軛斷層模型的聯(lián)合反演結(jié)果，為本文的最終破裂模型.根據(jù)該模型，地震的矩震級約為MW6.1.近南北向斷層上最大滑動量約為0.5m，發(fā)生在震中以南約8km處；近東西向斷層上最大滑動量約為0.3m，位于震中附近.如果只考慮滑動量大于0.2m的區(qū)域（橙色—紅色區(qū)域），則近南北向斷層上的破裂尺度約為12km，近東西向斷層僅約為3 km，前者遠(yuǎn)大于后者.

圖5顯示了采用共軛斷層模型進(jìn)行聯(lián)合反演得到的滑動速率時空分布圖像.由圖可見，近東西向斷層上的破裂開始較早，表現(xiàn)為從深部到淺部的擴(kuò)展過程，其破裂在6s左右結(jié)束.近南北向斷層上的破裂開始稍晚一些，但規(guī)模更大，主要破裂存在從深部向淺部和南南東方向傳播的過程.綜合兩條斷層上的破裂特征，在破裂開始后0～2s，滑動主要發(fā)生在近東西向斷層上；2s之后，近南北向斷層開始破裂.由此推測，在破裂過程的早期，近南北向斷層的破裂可能是受近東西向斷層上的破裂觸發(fā)所致.從破裂開始后2s直至破裂結(jié)束，兩條斷層上的滑動速率變化分布比較獨(dú)立，沒有表現(xiàn)出明顯的相干性.其中在破裂開始后4～7s，破裂抵達(dá)淺地表處，并在兩條斷層上形成了各自的最大滑動區(qū)域.破裂開始8s之后，主要滑動區(qū)域已經(jīng)形成，雖然從震源時間函數(shù)看，地震矩率的大小還比較可觀，但斷層面上已經(jīng)沒有比較集中的滑動分布，意味著破裂已經(jīng)進(jìn)入愈合階段.

圖6比較了以上破裂模型的觀測與合成波形.總體上，近震和遠(yuǎn)震記錄的主要波形特征都得到了較好地解釋.其中，一些臺站處的近震記錄受臺站效應(yīng)的影響，表現(xiàn)出一定的高頻振蕩特征.由于不是有效的震源破裂信號，且局部場地效應(yīng)很難在地球模型中予以考慮，這部分波形的擬合相對較差.在遠(yuǎn)震波形擬合方面，聯(lián)合反演與遠(yuǎn)震記錄單獨(dú)反演（張勇等，2014）基本相當(dāng).這意味著在需要兼顧近震資料擬合的情況下，遠(yuǎn)震記錄的資料解釋并未明顯變差，表明我們從近震和遠(yuǎn)震記錄中有效提取了二者所共有的震源破裂信息.

圖4 近震和遠(yuǎn)震數(shù)據(jù)聯(lián)合反演結(jié)果（其他說明參見圖3）Fig.4 Same as Fig.3，but for joint inversion of local and teleseismic data

圖5 基于共軛斷層模型的聯(lián)合反演結(jié)果（a）斷層面上滑動量累積分布在地面投影的時空變化；（b）斷層面上滑動速率的時空分布.Fig.5 Results of joint inversion based on the two conjugated faults（a）Spatial and temporal variations of surface projection of cumulative fault slips；（b）Spatial and temporal distribution of slip－rate on the faults.

4 討論與結(jié)論

本文考慮了單一斷層破裂和共軛斷層同時破裂的幾種可能的情形，分別采用近震資料，以及近震與遠(yuǎn)震資料對2014年魯?shù)榈卣鸬钠屏堰^程進(jìn)行了反演.近震資料單獨(dú)反演表明，魯?shù)榈卣鸬膬蓷l共軛斷層——近南北向斷層和近東西向斷層——都發(fā)生了破裂，但主要地震矩釋放過程發(fā)生在近南北向的斷層上，且破裂主要朝南南東方向擴(kuò)展.近、遠(yuǎn)震數(shù)據(jù)反演顯示，破裂最初可能發(fā)生在近東西向斷層上，并在破裂開始后2～3s觸發(fā)了與其共軛的近南北向斷層的破裂.此后，兩條斷層上的破裂相對獨(dú)立地擴(kuò)展.近東西向斷層的破裂主要向淺地表處延伸，而近南北向斷層上的破裂則朝地表和南南東方向傳播，最終形成了復(fù)雜的滑動分布.本文得到的這一破裂特征與區(qū)域距離上獨(dú)立得到的視震源時間函數(shù)分析結(jié)果可以很好地相互印證（許力生等，2014）.

需要注意的是，在共軛斷層模型的近震資料反演和聯(lián)合反演中，破裂在近南北向斷層上的傳播距離與余震分布比較吻合；但在近東西向斷層上的西段，破裂延伸范圍明顯小于余震分布的尺度，不能解釋該方向上余震密集分布的現(xiàn)象.根據(jù)本文反演結(jié)果，我們認(rèn)為近南北向斷層上位于震中以南的破裂可能是形成這一異常余震活動的主要原因.由于該段破裂以左旋走滑為主，且是整個地震中占主導(dǎo)的破裂，導(dǎo)致在破裂區(qū)域以東形成向北的推擠作用（圖7），增大了近東西向斷層?xùn)|段上下盤間的法向壓力與摩擦阻力，造成其進(jìn)一步閉鎖，因而該方向沒有密集余震發(fā)生；而在近南北向斷層以西，情況正好相反，主震破裂形成了向南的拉伸作用，導(dǎo)致近東西向斷層西段在一定程度上解鎖.由于地震前該斷層上閉鎖程度（如斷層面粗糙度）的非均勻性，以及積累的應(yīng)變能的有限性，這種解鎖作用沒有導(dǎo)致大規(guī)模破裂的發(fā)生，而是形成了密集的小規(guī)模破裂事件（余震）.這是近東西向斷層西段沒有破裂發(fā)生但卻存在密集的余震分布的主要原因.

圖6 基于共軛斷層模型的聯(lián)合反演得到的觀測地震波形（黑線）與合成地震波形（紅線）的比較（a）近震數(shù)據(jù)擬合；（b）遠(yuǎn)震數(shù)據(jù)擬合.Fig.6 Comparisons between observed（black lines）and synthetic（red lines）waves of joint inversion based on the conjugate－fault model（a）Local seismic data；（b）Teleseismic data.

另外，精定位結(jié)果也顯示，余震分布盡管呈 “倒L”型，但在不同方向上的發(fā)生時間存在差別（王未來等，2014；張廣偉等，2014）.早期余震主要沿近南北向斷層分布，表明是在主震破裂區(qū)域附近介質(zhì)應(yīng)力集中情況下發(fā)生的后續(xù)小規(guī)模破裂事件；震中以西（尤其是最西端）余震活動開始于震后1.8h，且在震后1.8h至震后1天的活動性都較弱，可能是主震破裂在近東西向斷層上的解鎖作用所致（圖7）.可見這種共軛斷層的余震觸發(fā)模式可以與兩個方向上余震發(fā)生的先后次序很好地對應(yīng).實(shí)際上，位于主震破裂區(qū)域之外的余震發(fā)生較晚的現(xiàn)象也可以從其他地震中找到例證，比如2010年青海玉樹MW6.9地震，震中東南且位于主震破裂區(qū)域周邊的余震發(fā)生較早，而位于主震破裂區(qū)域之外（震中西北約40km處）的余震則大都發(fā)生較晚（王未來等，2012）.

圖7 共軛斷層破裂與余震關(guān)系示意圖Fig.7 Sketch of unlocking process，which was caused by the left lateral strike－slips on the south part of the NNW－SSE fault，and impacts the west part of the ENE－WSW fault

綜上所述，結(jié)合本文確定的地震共軛斷層破裂模型，以及魯?shù)榈卣鹩嗾鹁ㄎ唤Y(jié)果（王未來等，2014；張廣偉等，2014），我們確定魯?shù)榈卣鹗窃诒蔽鳌蠔|向主壓應(yīng)力方向與北東—南西向主張應(yīng)力方向的統(tǒng)一應(yīng)力場下發(fā)生的兩條共軛的斷層先后破裂的一次復(fù)雜地震事件，該地震破裂過程和余震活動可以概括為以下5個階段：

① 破裂開始后1～2s，開始階段，近東西向斷層上發(fā)生破裂，且其破裂一直集中在震中附近；

② 破裂開始后2～8s，主要破裂階段，近南北向斷層和近東西向斷層同時破裂，其中近南北向斷層上的破裂朝南傳播，并釋放了大部分地震矩；

③ 破裂開始后8～17s，破裂愈合階段，斷層面上破裂分布比較零散；

④ 破裂開始17s到震后1.8h，受主震破裂的影響，余震活動主要沿近南北向斷層分布；

⑤ 地震發(fā)生1.8h后，由于主震破裂導(dǎo)致的解鎖作用，余震逐漸開始在震中以西的近東西向斷層上分布.

致謝 本文所用遠(yuǎn)震寬頻帶地震數(shù)據(jù)取自IRIS數(shù)據(jù)中心；強(qiáng)震數(shù)據(jù)來自西南強(qiáng)震動臺網(wǎng)中心；近震寬頻帶數(shù)據(jù)來自云南省地震局.審稿專家對本文提供了有益的修改意見；中國科學(xué)院測量與地球物理研究所刁法啟博士對本文的數(shù)據(jù)處理工作提供了有益的意見.作者謹(jǐn)向以上單位和個人表示由衷的感謝.

Deng Q D，Zhang P Z，Ran Y K，et al.2003.Basic characteristics of active tectonics of China.Sci.ChinaSer.D.，46（4）：356－372.

Kennett B L N，Engdahl E R，Buland R.1995.Constraints on seismic velocities in the earth from traveltimes.Geophys.J.Intl.，122（1）：108－124.

Lei J S，Zhang G W，Xie F R，et al.2012.Relocation of the 10 March 2011Yingjiang，China，earthquake sequence and its tectonic implications.Earthq.Sci.，25（1）：103－110.

Liu C L，Zheng Y，Xiong X，et al.2014.Rupture process ofMS6.5 Ludian earthquake constrained by regional broadband seismograms.ChineseJ.Geophys.（in Chinese），57（9）：3028－3037，doi：10.6038／cjg20140927.

Meng L，Ampuero J P，Stock J，et al.2012.Earthquake in a Maze：Compressional rupture branching during the 2012MW8.6 Sumatra earthquake.Science，337（6095）：724－726.

Wang C Y，Mooney W D，Wang X L，et al.2002.Study on 3－D velocity structure of crust and upper mantle in Sichuan－Yunnan region，China.ActaSeismol.Sinica（in Chinese），24（1）：1－16.

Wang C Y，Chan W W，Mooney W D.2003.Three－dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications.J.Geophys.Res.，108（B9）：2442.

Wang R J.1999.A simple orthonormalization method for stable and efficient computation of Green′s functions.Bull.Seismol.Soc.Am.，89（3）：733－741.

Wang W L，Wu J P，F(xiàn)ang L H，et al.2013.Relocation of the Yushu MS7.1earthquake and its aftershocks in 2010from HypoDD.ScienceChina：EarthSciences，56（2）：182－191，doi：10.1007／s11430－012－4450－z.

Wang W L，Wu J P，F(xiàn)ang L H，et al.2014.Double difference location of the LudianMS6.5earthquake sequences in Yunnan province in 2014.ChineseJ.Geophys.（in Chinese），57（9）：3042－3051，doi：10.6038／cjg20140929.

Wen X Z，Du F，Yi G X，et al.2013.Earthquake potential of the Zhaotong and Lianfeng fault zones of the eastern Sichuan－Yunnan border region.ChineseJ.Geophys.（in Chinese），56（10）：3361－3372，doi：10.6038／cjg20131012.

Xiong S B，Zheng Y，Yin Z X，et al.1993.The 2－D structure and it′s tectonic implications of the crust in the Lijiang－Panzhihua－Zhehai region.ChineseJ.Geophys.（in Chinese），36（4）：434－444.

Xu L S，Zhang X，Yan C，et al.2014.Analysis of the Love waves for the source complexity of the LudianMS6.5earthquake.ChineseJ.Geophys.（in Chinese），57（9）：3006－3017，doi：10.6038／cjg20140925.

Xu X W，Jiang G Y，Yu G H，et al.2014.Discussion on seismogenic fault of the LudianMS6.5earthquake and its tectonic attribution.ChineseJ.Geophys.（in Chinese），57（9）：3060－3068，doi：10.6038／cjg20140931.

Yue H，Lay T，Koper K D.2012.Enéchelon and orthogonal fault ruptures of the 11April 2012great intraplate earthquakes.Nature，490（7419）：245－250，doi：10.1038／nature11492.

Zhang G W，Lei J S，Liang S S，et al.2014.Relocations and focal mechanism solutions of the 3August 2014Ludian，YunnanMS6.5earthquake sequence.ChineseJ.Geophys.（in Chinese），57（9）：3018－3027，doi：10.6038／cjg20140926.

Zhang H，Chen J，Ge Z.2012.Multi－fault rupture and successive triggering during the 2012Mw8.6Sumatra offshore earthquake.GeophysicalResearchLetters，39，L22305，doi：10.1029／2012GL053805.

Zhang Y，Xu L S，Chen Y T.2010.Source process of the 2010 Yushu，Qinghai，earthquake.ScienceChina：EarthSciences，53（9）：1249－1251.

Zhang Y，F(xiàn)eng W P，Chen Y T，et al.2012.The 2009L′AquilaMW6.3earthquake：a new technique to locate the hypocentre in the joint inversion of earthquake rupture process.Geophys.J.Int.，191（3）：1417－1426，doi：10.1111／j.1365－246X.2012.05694.x.

Zhang Y，Wang R J，Chen Y T，et al.2014.Kinematic rupture model and hypocenter relocation of the 2013Mw6.6Lushan earthquake constrained by strong－motion and teleseismic data.Seismol.Res.Lett.，85（1）：15－22.

Zhang Y，Xu L S，Chen Y T，et al.2014.Rupture process of the 3 August 2014Ludian，Yunnan，MW6.1（MS6.5）earthquake.ChineseJ.Geophys.（in Chinese），57（9）：3052－3059，doi：10.6038／cjg20140930.

Zhang Z G，Sun Y C，Xu J K，et al.2014.Preliminary simulation of strong ground motion for Ludian，Yunnan earthquake of 3 August 2014，and hazard implication.ChineseJ.Geophys.（in Chinese），57（9）：3038－3041，doi：10.6038／cjg20140928.

附中文參考文獻(xiàn)

劉成利，鄭勇，熊熊等.2014.利用區(qū)域?qū)掝l帶數(shù)據(jù)反演魯?shù)镸S6.5級地震震源破裂過程.地球物理學(xué)報，57（9）：3028－3037，doi：10.6038／cjg20140927.

王椿鏞，Mooney W D，王溪莉等.2002.川滇地區(qū)地殼上地幔三維速度結(jié)構(gòu)研究.地震學(xué)報，24（1）：1－16.

王未來，吳建平，房立華等.2014.2014年云南魯?shù)镸S6.5地震序列的雙差定位.地球物理學(xué)報，57（9）：3042－3051，doi：10.6038／cjg20140929.

聞學(xué)澤，杜方，易桂喜等.2013.川滇交界東段昭通、蓮峰斷裂帶的地震危險背景.地球物理學(xué)報，56（10）：3361－3372，doi：10.6038／cjg20131012.

熊紹柏，鄭曄，尹周勛等.1993.麗江—攀枝花—者海地帶二維地殼結(jié)構(gòu)及其構(gòu)造意義.地球物理學(xué)報，36（4）：434－444.

許力生，張旭，嚴(yán)川等.2014.基于勒夫波的魯?shù)镸S6.5地震震源復(fù)雜性分析.地球物理學(xué)報，57（9）：3006－3017，doi：10.6038／cjg20140925.

徐錫偉，江國焰，于貴華等.2014.魯?shù)?.5級地震發(fā)震斷層判定及其構(gòu)造屬性討論.地球物理學(xué)報，57（9）：3060－3068，doi：10.6038／cjg20140931.

張廣偉，雷建設(shè)，梁姍姍等.2014.2014年8月3日云南魯?shù)镸S6.5級地震序列重定位與震源機(jī)制研究.地球物理學(xué)報，57（9）：3018－3027，doi：10.6038／cjg20140926.

張勇，許力生，陳運(yùn)泰等.2014.2014年8月3日云南魯?shù)镸W6.1（MS6.5）地震破裂過程.地球物理學(xué)報，57（9）：3052－3059，doi：10.6038／cjg20140930.

張振國，孫耀充，徐建寬等.2014.2014年8月3日云南魯?shù)榈卣饛?qiáng)地面運(yùn)動初步模擬及烈度預(yù)測.地球物理學(xué)報，57（9）：3038－3041，doi：10.6038／cjg20140928.