米琦， 申文豪,2*， 史保平

1 中國科學(xué)院大學(xué)地球科學(xué)學(xué)院， 北京 100049 2 中國地震局地殼應(yīng)力研究所， 北京 100085

?

基于經(jīng)驗?zāi)Ｐ秃臀锢砟Ｐ脱芯?013MS7.0 蘆山地震余震序列

米琦1， 申文豪1,2*， 史保平1

1 中國科學(xué)院大學(xué)地球科學(xué)學(xué)院， 北京 100049 2 中國地震局地殼應(yīng)力研究所， 北京 100085

基于中國地震臺網(wǎng)中心2013MS7.0蘆山地震余震數(shù)據(jù)我們首先確定了余震空間分布范圍并根據(jù)G-R關(guān)系計算了主震后半小時內(nèi)的完備震級Mc=3.5，并且得到了ML≥3.5和ML≥3.0的地震在2001年至蘆山地震前的背景場地震發(fā)生率.通過Omori-Ustu經(jīng)驗定律和兩種Dieterich模型對蘆山地震余震發(fā)生率的擬合，我們發(fā)現(xiàn)階梯型Dieterich模型只能模擬p=1的情況，從而造成了模擬曲線與觀測數(shù)據(jù)的差別；前人研究表明震后滑移同樣是產(chǎn)生余震的原因，如果假設(shè)余震序列由主震靜態(tài)剪應(yīng)力Δτ和震后滑移共同作用所產(chǎn)生，我們數(shù)值模擬得到的對數(shù)型Dieterich模型能夠較好地推斷余震發(fā)生率R隨時間t增加而衰減的趨勢，能夠從物理機制上解釋MS7.0蘆山地震余震序列衰減指數(shù)大于1這一現(xiàn)象.通過對數(shù)型Dieterich模型的擬合并結(jié)合Andrews的方法，我們還得到MS7.0蘆山地震Aσ約為0.155 MPa，ta約為8.4年，這一值與前人研究結(jié)果十分接近.

余震衰減； 蘆山地震； Dieterich模型； 震后滑移

1 引言

余震是指主震發(fā)生后，圍繞主震斷層一定尺度區(qū)域范圍內(nèi)發(fā)生的一系列地震，主要表現(xiàn)為在該區(qū)域范圍內(nèi)地震活動性的顯著增加，且隨著時間增加而衰減.余震震級通常比主震要小，其影響范圍在空間上超過了主震斷層的破裂尺度.余震序列的持續(xù)時間通常為數(shù)年乃至數(shù)十年不等，在此期間區(qū)域地震活動性明顯高于背景場水平，并隨時間逐漸衰減.Omori(1894)最先發(fā)現(xiàn)余震發(fā)生率R正比于1/t，其中t為時間，自此這一用于擬合余震發(fā)生率隨時間衰減的經(jīng)驗?zāi)Ｐ图癠tsu (1961) 提出的改進(jìn)模型(Modified Omori Law, MOL)

(1)

已經(jīng)被成功用于全球范圍內(nèi)多次余震序列的擬合.地震學(xué)首先是一門觀測學(xué)科，在觀測現(xiàn)象未獲得明確物理解釋之前，像Omori-Ustu定律這種地震統(tǒng)計學(xué)觀測定律，或者說經(jīng)驗?zāi)Ｐ蛯τ谡J(rèn)知地震活動性、了解地震活動規(guī)律發(fā)揮著舉足輕重的作用(Ogata,1988).在公式(1)中，R(t)表示余震發(fā)生率，即單位時間里余震發(fā)生個數(shù)；公式右端t為主震后時間，K，c和p為待確定的系數(shù).在這三個參數(shù)中，p值的大小反映了余震發(fā)生率隨時間衰減的快慢，稱之為衰減指數(shù)，其取值一般在0.8～1.5之間，并且該參數(shù)與起算震級的選取沒有必然聯(lián)系(Utsu et al., 1995).也就是說，對于同一個主震引發(fā)的余震序列，不同起算震級得到的p值應(yīng)該是相同的.通過分析美國加利福尼亞地區(qū)地震數(shù)據(jù)，Ouillon和Sornette(2005)發(fā)現(xiàn)p值隨主震震級M的增加而增加.近期Hainzl和Marson(2008)利用非均勻的斷層面應(yīng)力變化模型和Dieterich 速率和狀態(tài)模型給出了對于上述現(xiàn)象可信的物理解釋.K值一般被視為余震序列產(chǎn)出量的大小(productivity)，與主震震級M存在正相關(guān)關(guān)系：K～10αM，其中α為大于零的常數(shù)(Hainzl and Marsan, 2008).在三個參數(shù)中，c值一直備受爭議，通常認(rèn)為c值表示主震后高水平余震發(fā)生率持續(xù)的時間(飽和期)，其值在0.01～1.0 day之間.Enescu等(2007)和Narteau等(2002)的研究認(rèn)為主震后初始階段出現(xiàn)的高水平余震發(fā)生率持續(xù)期 (即c值) 是余震序列衰減的真實特征，而Kagan(2004)則認(rèn)為該現(xiàn)象并非是真實的衰減特征，而是由于主震后較短時間內(nèi)地震目錄(尤其是中小地震目錄)缺失所導(dǎo)致的飽和現(xiàn)象.

隨著震源物理研究的深入，余震觸發(fā)的物理學(xué)機制在近些年來得到廣泛關(guān)注并取得一系列新的成果，許多物理模型開始被用于解釋余震發(fā)生率衰減現(xiàn)象.由于在統(tǒng)計地震學(xué)中的成功應(yīng)用，Omori-Ustu定律已經(jīng)成為這些物理模型的“試金石”，因此幾乎所有物理模型都試圖得到同Omori-Ustu定律相似樣式的衰減形式，并賦予K，c，p以明確的物理內(nèi)涵.在這些物理模型中，本文采用的是速率和狀態(tài)依賴性摩擦本構(gòu)關(guān)系并結(jié)合一維彈簧-滑塊模型所得到的理論解，通常稱為Dieterich 模型，選用該模型的原因是近些年來該模型在地震學(xué)中廣泛而成功的應(yīng)用.

汶川地震發(fā)生后，中國地震局及四川省地震局投入了大量的人力物力用于龍門山斷裂帶的地震活動性監(jiān)測，相較于汶川地震之前該區(qū)域的監(jiān)測能力有了明顯提高.2013年4月20日蘆山地震發(fā)生后，得益于監(jiān)測能力的改善，中國地震臺網(wǎng)中心為我們提供了一套相對完備的余震數(shù)據(jù).在本文中，借助這一余震數(shù)據(jù)體，我們將系統(tǒng)分析該余震序列的時空分布特征.首先，確定余震區(qū)域的范圍和主震后半小時(～0.02 day)內(nèi)截止震級的大小，以確保數(shù)據(jù)體的完備和有效.其次，通過Omori-Ustu定律和兩種Dieterich模型對MS7.0蘆山地震余震序列衰減的擬合，我們將分別討論經(jīng)驗?zāi)Ｐ秃碗A梯型Dieterich模型存在的優(yōu)缺點，進(jìn)而發(fā)現(xiàn)結(jié)合同震應(yīng)力變化和震后剪應(yīng)力變化率衰減的Dieterich模型可以很好地解釋和模擬MS7.0蘆山地震余震序列衰減p值大于1這一現(xiàn)象.另一方面，我們還利用這一模型討論了蘆山地震的余震持續(xù)時間及Dieterich摩擦定律中參數(shù)的確定.

2 2013 MS7.0 蘆山地震余震序列時空分布特征

準(zhǔn)確地描述余震序列的時空特征對于我們正確認(rèn)知斷層內(nèi)部演化狀態(tài)，了解構(gòu)造環(huán)境與地震活動性的關(guān)系都起著重要作用.2013MS7.0蘆山地震的發(fā)震斷層一度被認(rèn)為是龍門山斷裂帶南端的大川—雙石斷裂，而李傳友等(2013)和徐錫偉等(2013)的研究結(jié)果則認(rèn)為蘆山地震的發(fā)震斷層為一條現(xiàn)今尚未出露地表、其上斷點仍埋藏在地下9 km以下地殼中的一條盲逆斷層，該發(fā)震構(gòu)造是蘆山之下、大川—雙石斷裂和新開店斷裂之間的龍門山前緣滑脫帶，此滑脫帶在該段的運動導(dǎo)致了這次地震的發(fā)生，并可能帶動了它上面的大川—雙石和新開店等斷裂的活動.選取29.8°N—30.5°N，102.5°E—103.3°E作為研究區(qū)域，該區(qū)域從2001年到汶川地震發(fā)生之前7年間ML≥3.5 地震有4次，而在汶川地震發(fā)生后到蘆山地震發(fā)生之前5年間ML≥3.5 地震有3次，基本維持在背景場水平，也就是說汶川地震的發(fā)生并未造成該區(qū)域內(nèi)地震活動性的增加，蘆山地震后該區(qū)域內(nèi)的地震數(shù)目激增是由蘆山主震觸發(fā)的.圖1給出了蘆山地震后到2014年5月25日ML≥3.5 余震分布，可以看到余震呈密集條帶狀分布，其東北至西南方向走向展布約40 km，西北至東南向展布約25 km.圖1右下角的子圖給出的是MS7.0蘆山地震和MS8.0汶川地震的余震(ML≥4.0，2008-05-12—2014-05-25)分布圖，可以看到兩次地震的余震密集區(qū)相距約50 km (杜方等，2013).

此外，由于在主震發(fā)生后的數(shù)小時內(nèi)中低震級余震的數(shù)目缺失，我們必須求出在主震后短時間內(nèi)最小完備震級(minimum magnitude of complete recording)，Mc.過高的Mc有可能導(dǎo)致一些可用數(shù)據(jù)的丟失，而如果低估了Mc則會得到錯誤的地震活動性參數(shù)值，從而導(dǎo)致分析中出現(xiàn)偏差(Mignan and Woessner, 2012).Mc通?？捎晒诺潜?里克特(G-R)定律擬合觀測到的震級-頻度分布(FMD)來估計，并且取兩者相符部分的起始震級作為Mc的估值(Zúniga and Wyss, 1995).表1給出了主震后半小時內(nèi)(～0.02 day)余震序列目錄(中國地震臺網(wǎng)中心，2014).

表1 MS 7.0蘆山地震后半小時內(nèi)(～0.02 day)主震及余震序列目錄

利用表1的余震序列數(shù)據(jù)，我們估算了在主震后半小時內(nèi)(～0.02 day)余震序列的完備震級Mc的取值為3.5，并得到了G-R關(guān)系中a和b的取值

圖1 MS7.0蘆山地震震中位置及其余震序列分布. 統(tǒng)計時間為2013年4月20日到2014年5月25日, 震級ML≥3.5藍(lán)色五角星為主震震中位置，紅色圓點為余震位置.Fig.1 Map of the location of MS 7.0 Lushan earthquake and its aftershocks distribution from the mainshock to May 25, 2014 for ML≥3.5 lue star indicated epicenter location of the mainshock, and red dots indicated aftershock locations.

圖2 (a) 2001年至主震發(fā)生前研究區(qū)域震級-頻度分布圖; (b) MS7.0蘆山地震序列主震后半小時內(nèi)震級-頻度分布圖Fig.2 (a) Magnitude-frequency distribution in the study area from 2001 to the mainshock; (b) Magnitude -frequency distribution of aftershocks for MS 7.0 Lushan earthquake in the first half hour after the main shock

(圖 2b).同樣，根據(jù)2001年至蘆山地震前的背景場地震數(shù)據(jù)我們得到了這段時間的完備震級為1.4 (圖 2a)，明顯小于主震后半小時內(nèi)余震序列的完備震級Mc= 3.5，說明在背景場地震數(shù)據(jù)中ML≥3.5的地震目錄是完整的，根據(jù)G-R關(guān)系我們估算得到ML≥3.5的地震在2001年至蘆山地震前的背景場發(fā)生率為0.0016/day.

3 利用經(jīng)驗?zāi)Ｐ蛿M合蘆山地震余震序列

選取研究區(qū)域內(nèi)2013年4月20日主震發(fā)生后至2014年5月25日ML≥3.5的余震數(shù)據(jù)，并采用Narteau等(2002)給出的計算方法，得到了該段時間余震發(fā)生率隨時間的變化(圖3).從圖3中可以看到實測余震發(fā)生率隨時間快速衰減，但是在震后約1.5天和震后20天左右出現(xiàn)了兩次異常突增，前一次地震活動性的突增主要是由于在2013年4月21日17時5分左右發(fā)生的一次ML=5.4的余震引起的，此次地震也是蘆山地震迄今為止最大的余震；第二次則是由于在2013年5月11日一天內(nèi)發(fā)生了三次4.0級以上的余震，引起了該時間段內(nèi)地震活動性的突增.通過最小二乘法的擬合，我們得到公式(1)與觀測數(shù)據(jù)擬合最優(yōu)時K，c，p的取值分別為28.12，1.21和0.05 (圖3).從圖3中可以看到利用Omori-Ustu經(jīng)驗公式可以很好地擬合實測數(shù)據(jù)，但是其參數(shù)c的取值在學(xué)術(shù)界存在較大爭議.Vidale等(2003)發(fā)現(xiàn)通過對在主震后數(shù)分鐘的波形數(shù)據(jù)進(jìn)行濾波得到的余震事件數(shù)目要數(shù)倍于地震目錄記錄到的余震事件，因此c值是震后地震目錄缺失產(chǎn)生的.Kagan(2004)通過計算地震矩釋放率發(fā)現(xiàn)c的取值接近0，并且應(yīng)該是一個負(fù)值.這些發(fā)現(xiàn)都更新了研究者對Omori-Ustu定律的認(rèn)識，而正如前文中提到的那樣，Omori-Ustu定律最大的問題在于其雖然對物理現(xiàn)象描述很好但是終究缺乏明確的物理機制解釋.

4 利用物理模型擬合蘆山地震余震序列

近些年來隨著地震科學(xué)的發(fā)展，越來越多的研究者致力于從震源物理機制上解釋主震后余震活動性的時空衰減規(guī)律，這些物理機制包括震后蠕變，流體擴散，速率和狀態(tài)依賴性摩擦定律，應(yīng)力腐蝕，損傷力學(xué)和亞臨界裂紋擴展機制等.在本研究中針對MS7.0蘆山地震我們使用的是速率和狀態(tài)依賴性摩擦定律即Dieterich模型,該模型為我們提供了一個對斷層屬性進(jìn)行復(fù)雜定量試驗觀測的基本框架(劉博研等，2013).如果給定不同的選取參數(shù)，Dieterich既能模擬斷層摩擦滑動的穩(wěn)定狀態(tài)又能模擬非穩(wěn)定狀態(tài).在失穩(wěn)狀態(tài)下，Dieterich模型能夠提供不同應(yīng)力加載形式與地震活動性之間的對應(yīng)關(guān)系，而這種關(guān)系可以用于預(yù)測不同應(yīng)力變化所觸發(fā)的余震發(fā)生率.這些應(yīng)力變化的來源是多方面的，包括主震造成的靜態(tài)應(yīng)力或動態(tài)應(yīng)力，震后滑移，巖漿入侵或火山噴發(fā)造成的瞬態(tài)形變，慢地震或者潮汐.

Dieterich(1994)引入一個應(yīng)力演化參數(shù)γ用來描述斷層上加載的剪應(yīng)力隨時間變化的過程(假設(shè)正應(yīng)力不變)，

(2)

圖3 MS7.0蘆山地震主震后ML≥3.5余震發(fā)生率隨時間變化示意圖Fig.3 Observations and modeling of seismicity rate during Lushan MS7.0 earthquake aftershocks sequence for earthquake ML≥3.5 Circles indicated observed data, the green solid line indicated the fitting curve of Omori-Ustu law, the red and purple solid line indicated fitting curve of Eqs. (5) and (7) respectively, and the blue solid line indicated the curve fitting by the logarithmic stressing Dieterich model.

(3)

4.1 階梯型Dietrich模型擬合

假設(shè)余震序列僅由主震造成的靜態(tài)剪應(yīng)力Δτ產(chǎn)生，則主震后的應(yīng)力演化參數(shù)γ初值可以寫為：

(4)

(5)

(6)

將公式(6)重新整理可以得到：

(7)

針對MS7.0蘆山地震余震序列我們分別用公式(5)和公式(7)來擬合，在這兩個公式中需要調(diào)整兩個參數(shù)值，一個是ta，另一個是Δτ/Aσ的值.通過網(wǎng)格搜索法，我們得到擬合優(yōu)度最高時ta=2.3×104day，約63年；Δτ/Aσ=13.5(圖3).從圖3中可以看到在余震序列的前一部分公式(5)和公式(7)與觀測數(shù)據(jù)符合較好，而在后一部分明顯與觀測數(shù)據(jù)產(chǎn)生分離.由Omori-Ustu定律的擬合我們知道MS7.0蘆山地震余震發(fā)生率的衰減指數(shù)p=1.21，而公式(5)和公式(7)卻只能模擬p=1的情況，因此造成了模擬曲線與觀測數(shù)據(jù)的分離，事實上這一點正是階梯型Dieterich模型在應(yīng)用中所存在的最大問題.大量的地震觀測數(shù)據(jù)表明衰減指數(shù)p值介于0.8～1.5(Utsu et al., 1995)，階梯型Dieterich模型顯然不能解釋p≠1的情況.為解決這一問題Hainzl和Marson(2008)，Helmstetter和Shaw(2006)等研究者引入了運動學(xué)有限斷層模型(如k-2模型)，將主震后的靜態(tài)應(yīng)力變化視為不均勻的(有正有負(fù))而非單一值，這樣就可以得到p<1的階梯型Dieterich模型，但是這種模型依然不能解釋p>1的情況，因此在該模型下得到的余震持續(xù)時間往往會被高估，進(jìn)而導(dǎo)致余震地震危險性分析的誤差.

4.2 對數(shù)型Dietrich模型擬合

除了余震現(xiàn)象之外，在主震后還會發(fā)生明顯的震后滑移，而大量研究則表明余震的衰減與震后位移的衰減遵從相同的時間衰減規(guī)律(Perfettini and Avouac, 2007)，這一現(xiàn)象使得許多研究者開始關(guān)注震后滑移與余震之間的關(guān)系，并且越來越多的證據(jù)表明震后滑移能夠觸發(fā)余震(Perfettini and Avouac,2007; Hsu et al., 2006; Bourouis and Bernard, 2007).Dieterich(1994)，Helmstetter和Shaw(2009)以及Hainzl和Marsan(2008)等研究給出了震后滑移引起的剪應(yīng)力變化率解析形式：

(8)

(9)

(10)

(11)

圖4 基于對數(shù)型Dieterich模型得到的地震發(fā)生率示意圖.Aσ=0.1 MPa，Δτ=0.92 MPa，t*=1.0 day，分別取-1.5，-0.5，0，0.5和1.5.m= 0時表示僅由同震靜態(tài)應(yīng)力誘發(fā)余震Fig.4 Seismicity rate triggered by a continuous stress change given by Aσ=0.1 MPa，Δτ=0.92 MPa，t*=1.0 r =10-6 MPa/day and for different amplitudes of the postseismic stress change, from m= -1.5 (bottom) to m=1.5 (top). The case m=0 represents seismicity rate triggered only by the coseismic stress step

圖6 MS7.0蘆山地震主震后ML≥3.0余震發(fā)生率隨時間變化示意圖Fig.6 Observations and modeling of seismicity rate during Lushan Ms7.0 earthquake aftershocks sequence for earthquake ML≥3.0 Circles indicated observed data, the green solid line indicated the fitting curve of Omori-Ustu law, the red and purple solid line indicated fitting curve of Eqs. (5) and (7) respectively, and the blue solid line indicated the curve fitting by the logarithmic stressing Dieterich model.

除了ML≥3.5的地震之外，我們還給出了ML≥3.0余震序列實測數(shù)據(jù)與經(jīng)驗?zāi)Ｐ图拔锢砟Ｐ偷膶Ρ?圖6)，根據(jù)該區(qū)域G-R關(guān)系我們計算得到ML≥3.0余震序列背景場地震發(fā)生率為0.0037/day.在使用階梯型Dieterich模型和對數(shù)型Dieterich模型進(jìn)行擬合時，所用到的參數(shù)除背景場地震發(fā)生率r外其他參數(shù)與ML≥3.5余震序列擬合中采用參數(shù)完全一致.另外，用Omori-Ustu經(jīng)驗定律對ML≥3.0余震序列擬合時K值調(diào)整為58.12.對比圖3和圖6，可以看到兩種截止震級得到的余震衰減指數(shù)一致，均為p=1.21，并未隨截止震級改變而改變，不同的是ML≥3.0余震序列在主震后初期由于數(shù)據(jù)不完整，發(fā)生率明顯低于模型預(yù)測結(jié)果.從圖3和圖6中可以看到對數(shù)型Dieterich模型預(yù)測的余震發(fā)生率衰減趨勢與觀測數(shù)據(jù)符合較好，并且衰減指數(shù)與經(jīng)驗?zāi)Ｐ鸵恢?，因此我們可以推測MS7.0蘆山地震余震序列是由主震靜態(tài)應(yīng)力和震后滑移共同作用產(chǎn)生的.

5 結(jié)果與討論

為分離各個參數(shù)對模擬曲線的影響，我們設(shè)定在每個方案中只變化一個參數(shù)的量，其他參數(shù)采用表2給出的標(biāo)準(zhǔn)值.每個方案所調(diào)整的參數(shù)名和值如表3所示，每一方案的最終結(jié)果由圖7表示.

表2 對數(shù)型Dieterich模型模擬余震衰減使用的標(biāo)準(zhǔn)參數(shù)

表3 參數(shù)劃分方案

無論是階梯型Dieterich模型還是對數(shù)型Dieterich模型，又或是其他基于速率和狀態(tài)依賴性摩擦定律的物理模型，斷層摩擦本構(gòu)參數(shù)A和正應(yīng)力σ的取值都起著非常關(guān)鍵的作用.現(xiàn)有研究表明，A值主要與溫度和斷層巖性有關(guān)，而正應(yīng)力σ則與深度、區(qū)域應(yīng)力場、斷層方位和孔隙壓力等因素有關(guān)系.在實際應(yīng)用中Aσ一般被視為常數(shù)，其中A的參考值0.005～0.012是在實驗室得到的，而在天然地震應(yīng)用中如何取值一直沒有明確限制.在許多基于速率和狀態(tài)依賴性摩擦定律的研究中，Aσ的估算值為0.01～0.1 MPa(Hainzl et al., 2010)，Perfettini和Avouac(2004,2007)針對1992年LandersMW7.3地震和1999年Chi-ChiMW7.6地震分別得到Aσ約為0.47～0.53 MPa和0.34～1.5 MPa.Cochran等(2004)在對潮汐誘發(fā)地震的研究中給出Aσ的區(qū)間為0.048～0.11 MPa.在本文中針對MS7.0蘆山地震，基于階梯型Dieterich模型，得到Δτ/(Aσ)=13.5，根據(jù)上述已經(jīng)求出Δτ約為2.0 MPa，所以Aσ=0.148 MPa，基于對數(shù)型Dieterich模型，求得Aσ=0.155 MPa.如果取A約為0.005～0.012，可以得到斷層面上的正應(yīng)力σ約為12.3～31.0 MPa.

圖7 不同參數(shù)對對數(shù)型Dieterich模型余震衰減特征的影響Fig.7 Comparison of the decay for afterslip triggering Dieterich Model based on different parameters

Dieterich(1994)發(fā)現(xiàn)淺源地震的余震持續(xù)時間ta在10.2年左右，Parson(2002)的研究發(fā)現(xiàn)全球范圍內(nèi)主震誘發(fā)的余震持續(xù)時間ta約為7～11年.在圖3和圖6中，對數(shù)型Dieterich模型和Omori-Ustu定律預(yù)測的ta值為8.4年，這一值與上述結(jié)果十分接近，而階梯型Dieterich模型的結(jié)果明顯過大.另外從圖中還可以看到，在衰減到背景場地震發(fā)生率r后，對數(shù)型Dieterich模型預(yù)測R會繼續(xù)衰減一定時間，然后又逐漸增加恢復(fù)到背景場水平，也就是說在主震發(fā)生8.4年后研究區(qū)域有可能出現(xiàn)地震發(fā)生率低于背景場水平的現(xiàn)象，形成地震空區(qū)，而階梯型Dieterich模型則不能有效解釋地震空區(qū)現(xiàn)象.

盡管Dieterich 模型在實際應(yīng)用中取得巨大的成功，但是由于其本身來源于室內(nèi)的巖石物理學(xué)實驗，在該模型中一些重要本構(gòu)參數(shù)的取值(如A)是否能直接應(yīng)用于天然地震仍存在較大爭議，另外這些本構(gòu)參數(shù)是否為常數(shù)又或者與構(gòu)造環(huán)境有關(guān)系也存在不確定性.在運用該模型進(jìn)行模擬時我們還忽略了流體壓力、動態(tài)應(yīng)力變化等因素的影響，這也有可能給結(jié)果帶來不確定性.

6 結(jié)論

通過經(jīng)驗性的Omori-Ustu定律和兩種Dieterich模型對蘆山地震余震發(fā)生率的擬合我們發(fā)現(xiàn)階梯型Dieterich模型只能模擬p=1的情況，因此造成了模擬曲線與觀測數(shù)據(jù)的分離；而如果假設(shè)余震序列由主震靜態(tài)剪應(yīng)力Δτ和震后滑移共同作用所產(chǎn)生，這樣得到的對數(shù)型Dieterich模型能夠較好預(yù)測蘆山余震發(fā)生率衰減趨勢，與觀測數(shù)據(jù)符合較好，能夠從物理機制上解釋蘆山地震余震序列衰減指數(shù)大于1這一現(xiàn)象.除此之外對數(shù)型Dieterich模型還可以合理解釋地震空區(qū)現(xiàn)象，有效地彌補了階梯型Dieterich模型的缺陷.綜上所述，通過經(jīng)驗?zāi)Ｐ秃臀锢砟Ｐ偷膶Ρ?，我們有理由相信主震靜態(tài)應(yīng)力和震后滑移共同作用產(chǎn)生了蘆山地震余震序列，我們得到的模型參數(shù)結(jié)果對于研究蘆山地區(qū)地震危險性以及分析該區(qū)域未來地震活動性具有重要參考意義.

Bourouis S, Bernard P. 2007. Evidence for coupled seismic and aseismic fault slip during water injection in the geothermal site of Soultz (France), and implications for seismogenic transients.GeophysicalJournalInternational, 169(2): 723-732, doi: 10.1111/j.1365-246X.2006.03325.x.

Cochran E S, Vidale J E, Tanaka S. 2004. Earth tides can trigger shallow thrust fault earthquakes.Science, 306(5699): 1164-1166.

Dieterich J H. 1994. A constitutive law for rate of earthquake production and its application to earthquake clustering.JournalofGeophysicalResearch, 99(B2): 2601-2618.

Du F, Long F, Ruan X, et al. 2013. The M7.0 Lushan earthquake and the relationship with theM8.0 Wenchuan earthquake in Sichuan, China.ChineseJournalofGeophysics(in Chinese), 2013, 56(5): 1772-1783, doi: 10.6038/cjg20130535.

Enescu B, Mori J, Miyazawa M. 2007. Quantifying early aftershock activity of the 2004 mid-Niigata Prefecture earthquake (Mw6.6).JournalofGeophysicalResearch:SolidEarth, 112(B4), doi: 10.1029/2006JB004629.

Hainzl S, Steacy S, Marson D. 2010. Seismicity models based on Coulomb stress calculations.CommunityOnlineResourceforStatisticalSeismicityAnalysis, doi: 10.5078/corssa-32035809. Available at http:∥www.corssa.org.Hainzl S, Marsan D. 2008. Dependence of the Omori-Utsu law parameters on main shock magnitude: Observations and modeling.JournalofGeophysicalResearch:SolidEarth, 113: B10309, doi: 10.1029/2007JB005492.

Helmstetter A, Shaw B E. 2006. Relation between stress heterogeneity and aftershock rate in the rate-and-state model.JournalofGeophysicalResearch, 111: B07304, doi: 10.1029/2005JB004077.

Helmstetter A, Shaw B E. 2009. Afterslip and aftershocks in the rate-and-state friction law.JournalofGeophysicalResearch, 114: B01308, doi: 10.1029/2007JB005077.

Hsu Y J, Simons M, Avouac J P, et al. 2006. Frictional afterslip following the 2005 Nias-Simeulue earthquake, Sumatra.Science, 312(5782): 1921-1926, doi: 10.1126/science.1126960.Kagan Y Y. 2004. Short-term properties of earthquake catalogs and models of earthquake source.BulletinoftheSeismologicalSocietyofAmerican, 94(4): 1207-1228.

Li C Y, Xu X W, Gan W J, et al. 2013. Seismogenic structures associated with the 20 April 2013Ms7.0 Lushan earthquake, Sichuan province.SeismologyandGeology(in Chinese), 35(3): 671-683.

Liu B Y, Shi B P, Lei J S. 2013. Effect of Wenchuan earthquake on probabilities of earthquake occurrence of Lushan and surrounding faults.ActaSeismologicaSinica(in Chinese), 35(5): 642-651. Mignan A, Woessner J. 2012. Estimating the magnitude of completeness for earthquake catalogs. Community Online Resource for Statistical Seismicity Analysis, doi:10.5078/corssa-00180805. Available at http:∥www.corssa.org.

Narteau C, Shebalin P, Holschneider M. 2002. Temporal limits of the power law aftershock decay rate.JournalofGeophysicalResearch:SolidEarth, 107(B12): ESE 12-1-ESE 12-14, doi: 10.1029/2002JB001868.

Ogata Y. 1988. Statistical models for earthquake occurrences and residual analysis for point processes.JournaloftheAmericanStatisticalAssociation, 83(401): 9-27.

Omori F. 1894. On after-shocks of earthquakes.JournaloftheCollegeofScience,ImperialUniversityofTokyo, 7: 111-200.

Ouillon G, Sornette D. 2005. Magnitude-dependent Omori law: empirical study and theory.JournalofGeophysicalResearch:SolidEarth, 110: B04306, doi: 10.1029/2004JB003311.

Parsons T. 2002. Global Omori law decay of triggered earthquakes: Large aftershocks outside the classical aftershock zone.JournalofGeophysicalResearch, 107(B9): ESE 9-1-ESE 9-20, doi: 10.1029/2001JB000646.

Perfettini H, Avouac J P. 2004. Postseismic relaxation driven by brittle creep: A possible mechanism to reconcile geodetic measurements and the decay rate of aftershocks, application to the Chi-Chi earthquake, Taiwan.JournalofGeophysicalResearch, 109: B02304, doi: 10.1029/2003JB002488.

Perfettini H, Avouac J P. 2007. Modeling afterslip and aftershocks following the 1992 Landers earthquake.JournalofGeophysicalResearch, 112: B07409, doi: 10.1029/2006JB0041399. Ripperger J, Mai P M. 2004. Fast computation of static stress changes on 2D faults from final slip distributions.GeophysicalResearchLetter, 31(18): L18610, doi: 10.1029/2004GL020594. Shen W H, Qiu Z, Shi B P. 2014. Strong ground motion simulation for the 2013 LushanMw6.6 earthquake, Sichuan, China, based on the inverted and synthetic slip models.EarthquakeScience, 27(4): 377-389, doi: 10.1007/s11589-013-0057-5.

Utsu T. 1961. A statistical study on the occurrence of aftershocks.GeophysicalMagazine, 30: 521-605.

Utsu T, Ogata Y, Ritsuko S, et al. 1995. The centenary of the Omori formula for a decay law of aftershock activity.JournalofPhysicsoftheEarth, 43(1): 1-33.

Vidale J E, Cochran E S, Kanamori H, et al. 2003. After the lightning and before the thunder; non-Omori behavior of early aftershocks? Eos Trans AGU Fall Meet Suppl Abstract S31A-08, 84(46).

Wang W M, Hao J L, Yao Z X. 2013. Preliminary result for rupture process of Apr. 20, 2013, Lushan earthquake, Sichuan, China.ChineseJournalofGeophysics(in Chinese), 56(4): 1412-1417, doi: 10.6038/cjg20130436. Xu X W, Chen G H, Yu G H, et al. 2013. Seismogenic structure of Lushan earthquake and its relationship with Wenchuan earthquake.EarthScienceFrontiers(in Chinese), 20(3): 11-20.

Zhang P Z, Xu X W, Wen X Z, et al. 2008. 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.ChineseJournalofGeophysics(in Chinese), 51(4): 1066-1073, doi: 10.3321/j.issn:0001-5733.2008.04.015.

Zhang Y, Xu L S, Chen Y T. 2013. Rupture process of the Lushan

4.20 earthquake and preliminary analysis on the disaster-causing mechanism.ChineseJournalofGeophysics(in Chinese), 56(4): 1408-1411, doi: 10.6038/cjg20130435.Zúniga F R, Wyss M. 1995. Inadvertent changes in magnitude reported in earthquake catalogs: Their evaluation through b-value estimates.BulletinoftheSeismologicalSocietyofAmerican, 85(6): 1858-1866.

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

杜方, 龍鋒, 阮祥等. 2013. 四川蘆山7.0級地震及其與汶川8.0地震的關(guān)系. 地球物理學(xué)報, 56(5): 1772-1783, doi: 10.6038/cjg20130535.

李傳友, 徐錫偉, 甘衛(wèi)軍等. 2013. 四川省蘆山Ms7.0地震發(fā)震構(gòu)造分析. 地震地質(zhì), 35(3): 671-683.

劉博研, 史保平, 雷建設(shè). 2013. 汶川地震對蘆山地震及周邊斷層發(fā)震概率的影響. 地震學(xué)報, 35(5): 642-651.

王衛(wèi)民, 郝金來, 姚振興. 2013. 2013年4月20日四川蘆山地震震源破裂過程反演初步結(jié)果. 地球物理學(xué)報, 56(4): 1412-1417, doi: 10.6038/cjg20130436.

徐錫偉, 陳桂華, 于貴華等. 2013. 蘆山地震發(fā)震構(gòu)造及其與汶川地震關(guān)系討論. 地學(xué)前緣, 20(3): 11-20.

張培震, 徐錫偉, 聞學(xué)澤等. 2008. 2008年汶川8.0級地震發(fā)震斷裂的滑動速率、復(fù)發(fā)周期和構(gòu)造成因. 地球物理學(xué)報, 51(4): 1066-1073, doi: 10.3321/j.issn:0001-5733.2008.04.015.

張勇, 許力生, 陳運泰. 2013. 蘆山4. 20地震破裂過程及其致災(zāi)特征初步分析. 地球物理學(xué)報, 56(4): 1408-1411, doi: 10.6038/cjg20130435.

中國地震臺網(wǎng)中心. 2014. http:∥www.csndmc.ac.cn/newweb/catalog_direct_link.htm.

(本文編輯 何燕)

Aftershock decay of the 2013 LushanMS7.0 earthquake derived from the empirical and physical models

MI Qi1， SHEN Wen-Hao1，2*， SHI Bao-Ping1

1CollegeofEarthScience,UniversityofChineseAcademyofSciences,Beijing100049,China2InstituteofCrustalDynamics,ChinaEarthquakeAdministration,Beijing100085,China

Accurately predicting aftershocks decay after the main shock over time is conducive to learn the aftershocks triggering physical mechanism. On the other hand, the predicted duration of aftershocks is also an important aspect of seismic hazard analysis. In this study, by fitting the seismicity decay with both the empirical and physical models, we try to study Lushan earthquake triggering mechanism, and figure out how coseismic stress, afterslip and other parameters affect aftershocks decay.First, based on the seismicity data from the China Earthquake Network Center, the magnitude of completenessMcfor the early part was estimated using a maximum likelihood procedure. We then compute the observed seismicity decay and fit it with both the empirical and physical models. By using a grid search method, parameters of predication models could be determinated.The magnitude of completenessMcfor the early part (～0.02 day) was estimated at 3.5 using a maximum likelihood procedure. We found that, with the stress step formed Dieterich model, apvalue larger than 1 could not be explained, which was in contradiction with empirical observations ofp>1 aftershock decays. We have shown that by modeling a coseismic step and a postseismic stress change with a constant background stress rate, the Omori exponent larger than one for theMS7.0 Lushan Earthquake aftershock sequence could be well explained and fitted. With Andrews′ method, we calculated the distributed stress changes on the fault plane, and we inferredAσ～0.155 MPa, also, the logarithmic stressing Dieterich model predicts a more reasonable aftershocks durationta～8.4 years, which is close to former researchers′ result.By showing the merits and limitations of the modified Omori law and the step function formed Dieterich model, we find that the stress step formed Dieterich model cannot explain an Omori law decay withp>1, thus, this physical model may overestimate the aftershock duration. We show that modeling by a coseismic step, a postseismic stress change with a constant background stress rate, the Omori exponent larger than one forMS7.0 Lushan earthquake could be explained.

Aftershock decay; Lushan earthquake; Dieterich model; Afterslip

10.6038/cjg20150608.

中國科學(xué)院創(chuàng)新團(tuán)隊項目(KZZD-EW-TZ-19)資助.

米琦，男，博士研究生，主要從事地應(yīng)力場數(shù)值模擬研究. E-mail: miqi3007@163.com

*通訊作者 申文豪，男，博士研究生. E-mail: shenxiaoxi200212@163.com

10.6038/cjg20150608

P315

2014-06-27，2014-10-29收修定稿

米琦，申文豪，史保平. 2015. 基于經(jīng)驗?zāi)Ｐ秃臀锢砟Ｐ脱芯?013MS7.0 蘆山地震余震序列.地球物理學(xué)報，58(6):1919-1930,

Mi Q, Shen W H, Shi B P. 2015. Aftershock decay of the 2013 LushanMS7.0 earthquake derived from the empirical and physical models.ChineseJ.Geophys. (in Chinese),58(6):1919-1930,doi:10.6038/cjg20150608.