劉博研，史保平，雷建設(shè)

1中國地震局地殼應(yīng)力研究所（地殼動力學(xué)重點實驗室），北京 100085

2中國科學(xué)院大學(xué)地球科學(xué)學(xué)院，北京 100049

1 引言

前人對于動態(tài)和靜態(tài)應(yīng)力變化的研究表明，同震過程導(dǎo)致了斷層間的相互作用和地震的觸發(fā).靜態(tài)庫侖應(yīng)力變化常常被用于解釋區(qū)域地震活動性速率的變化（Reasenberg and Simpson，1992；Simpson and Reasenberg，1994；Toda et al.，1998；Harris and Simpson，1998）、主震斷層外余震的發(fā)生（Rybicki，1973；Das and Scholz，1981；Stein et al.，1992，1994；Dieterich，1972，1994；Dieterich and Kilgore，1996；Gross and Kisslinger，1997；Gross and Bürgmann，1998）和即將失穩(wěn)的斷層地震發(fā)生概率的變化（Toda et al.，1998）；而動態(tài)模型被用來研究由動態(tài)應(yīng)力變化引起的斷層間相互作用和地震的觸發(fā)（Harris and Day，1993；Gomberg et al.，1997，1998）.1992年美國Landers地震為我們提供了一些小范圍內(nèi)斷層相互作用和遠距離地震活動性觸發(fā)的實例.King等（1994）計算得到 Landers地震對后續(xù)M6.5Big Bear地震產(chǎn)生的庫侖破裂應(yīng)力增量為0.2MPa左右，認為是Landers地震加速了Big Bear地震的發(fā)生.Spudich等（1995）還討論了關(guān)于Big Bear余震是由1992年Landers地震觸發(fā)的問題，并指出不論用動態(tài)觸發(fā)還是靜態(tài)觸發(fā)，都很難解釋這3.5個小時的時間延遲.Gomberg等（1998）解釋了摩擦失穩(wěn)模型中近距離（靜態(tài)和動態(tài)）觸發(fā)和遠距離（動態(tài)）觸發(fā)的主要特征，并討論了靜態(tài)和動態(tài)加載是怎樣改變了斷層的失穩(wěn)時間的.而Harris（1998）詳細總結(jié)了應(yīng)力觸發(fā)、應(yīng)力影區(qū)和其對地震危險性的影響.關(guān)于應(yīng)力觸發(fā)的相關(guān)研究受到了國內(nèi)外學(xué)者的重視.特別是2013年蘆山地震發(fā)生后，很多學(xué)者認為，此次地震與2008年汶川地震有密切關(guān)系（陳運泰等，2013；單斌等，2013；杜方等，2013；雷興林等，2013）.Parsons等（2008）初步探討了汶川地震產(chǎn)生的應(yīng)力變化和對周圍斷層的影響.劉博研等（2013）定量計算了汶川地震后，雅安地區(qū)庫侖應(yīng)力變化導(dǎo)致的發(fā)震概率的增加.

蘆山地震發(fā)生在汶川地震五年后，是汶川地震后發(fā)生在龍門山斷裂帶上最強的一次地震.該地震發(fā)生在龍門山斷裂帶的南端，在汶川地震發(fā)生時這一地段沒有發(fā)生破裂滑動.與此相似的是，2008年3月21日6時33分02秒，在新疆于田、策勒與西藏交界的昆侖山阿什庫勒火山南側(cè)發(fā)生過一次MS7.3地震，該地震發(fā)生在距印度板塊和歐亞板塊碰撞帶數(shù)百公里的青藏高原北緣，位于東北向的阿爾金斷裂和東西向的普魯斷裂之間.在六年后，2014年2月12日17時19分，新疆和田地區(qū)于田縣又發(fā)生了一次MS7.3級地震，震中位置為36.1°N，82.5°E，震源深度12km.該次強地震位于阿爾金斷裂帶西端部馬尾狀分支斷裂上，為帶NE向左旋走滑的正斷層地震（徐錫偉和于貴華，2014）.兩次地震的震中相距約120km.青藏高原處在7至8級地震活動的叢集期，從1997年至今發(fā)生在中國大陸的強地震基本都分布在巴顏喀拉地塊周緣邊界上，但由于大地震的發(fā)生需要長時間的能量積累，不論是2001年11月14日昆侖山地震（MS8.1）、2008年3月21日于田地震（MS7.3）、2008年5月12日汶川地震（MS8.0），還是2010年4月21日玉樹地震（MS7.1），盡管同在巴顏喀拉地塊邊界發(fā)生，但彼此距離都在600km以上.然而，2013年在蘆山發(fā)生的MS7.0地震和2014年MS7.3于田地震，都發(fā)生在5～6年前曾經(jīng)發(fā)生過大地震的周邊，這必將引起我們的重視.如果說蘆山地震第一次提醒人們大地震過后周圍斷層的地震危險性并不僅僅由余震造成，那么2014年的于田地震則再一次為人們敲響了警鐘.

2008年于田地震的發(fā)生必定帶來了周邊斷層的庫侖應(yīng)力變化.究竟2008年于田地震與2014年于田地震有什么關(guān)系，是否是2008年的地震觸發(fā)了2014年的于田地震？2014年于田地震后是否還會有類似的強地震在周圍地區(qū)發(fā)生？這才是研究大震后區(qū)域庫侖應(yīng)力變化的最主要目的和最核心的問題.

本文利用滑移速率和狀態(tài)相依賴的摩擦定律（Dieterich，1994），結(jié)合2008年于田地震前后的地震活動性水平，深入闡述了2014年于田地震的靜態(tài)應(yīng)力觸發(fā)過程，從地震發(fā)生概率的變化上理解2014年于田地震的發(fā)震根源，并計算周邊斷層經(jīng)歷兩次于田地震后的發(fā)震概率變化，從而為區(qū)域地震活動性以及危險性評價提供重要的參考.

2 滑移速率和狀態(tài)相依賴的摩擦定律

滑移速率和狀態(tài)相依賴的摩擦本構(gòu)關(guān)系為我們提供了一個對斷層屬性進行復(fù)雜的定量實驗觀測的基本框架（Gomberg et al.，1998；Dieterich，1979，1981；Ruina，1983）.最簡單的表述滑移速率和狀態(tài)相依賴的斷層應(yīng)力是由Ruina（1983）基于Dieterich（1994）的理論提出的.如果將其推廣為多態(tài)變量（Rice and Gu，1983），可以寫作：

其中，τ和σN分別為剪切應(yīng)力和正應(yīng)力，δ·是滑移速率，θi是狀態(tài)變量；參數(shù)μ0、A和Bi為通過實驗得到的系數(shù)；具有星號的項為標準化常量（Kanamori and Brodsky，2004）.假設(shè)斷層的摩擦過程由一維彈簧－滑塊模型描述（圖1），剪切摩擦應(yīng)力τ＝－k（δδ0），其中k為彈簧的有效彈性系數(shù)，δ和δ0分別為滑動位移和初始滑動位移.若剪應(yīng)力變化率τ·為常數(shù)且不等于0，其加載過程由τ（t）＝τ0＋描述，可得到滑動速率為（Kanamori and Brodsky，2004）：

圖1 彈簧－滑塊模型示意圖，接觸面上的剪切應(yīng)力為：τ＝－k（δ－δ0）Fig.1 Spring－slide block model.The shear stress on surface is：τ＝－k（δ－δ0）

進一步，如果假設(shè)第n次子事件的失穩(wěn)時間tf＝nΔt，其中Δt是n－1次事件到n次事件的時間間隔，那么可以得到第n次子事件的滑動速率：

在初始時刻t0時由于主震使得周圍斷層突然有一個Δτ的應(yīng)力加載，那么基于彈簧－滑塊模型（圖1），有kδ0＝τ0（t＜t0）和kδ0＝τ0＋Δτ（t＞t0）.為保證在t＝t0時滑動位移和滑動速率是連續(xù)的，有δ（t＝t0－ε）＝δ（t＝t0＋ε）和（t＝t0＋ε）＝（t＝即滑動速率增加了倍.此時失穩(wěn)時間變?yōu)椋?/p>

將公式（4）代入公式（5）可得地震（事件）次數(shù)：

此時，n和tf（n）都是離散變量.假設(shè)很多類型的地震是通過背景場地震活動性的擾動引起，而這個擾動又來自之前的地震造成的應(yīng)力場狀態(tài)的改變.將瞬態(tài)地震活動性速率R定義為R＝dn／dtf（n），并把n和tf（n）都看作連續(xù)函數(shù)，當參考剪切應(yīng)力變化率τ·

r與實際剪切應(yīng)力變化率τ·相等時，有其中，r是在參考場或背景場下區(qū)域內(nèi)剪切應(yīng)力變化率τ·

r下的常數(shù)穩(wěn)態(tài)背景地震活動性；ta為余震的持續(xù)時間，即余震地震活動性恢復(fù)到背景場的特征時間，與剪應(yīng)力的加載速率相關(guān)，當t＜ta時，公式（7）符合大森余震衰減定律，而t＞ta時，地震活動性恢復(fù)到背景場r的水平.需要說明的是，公式（7）中，R和r都為震級M的函數(shù).

由于地震活動性與區(qū)域應(yīng)力場有關(guān)，那么區(qū)域應(yīng)力場的改變必將導(dǎo)致相應(yīng)的地震危險性概率的變化.假定地震序列在t時間內(nèi)，對震級M以上的發(fā)震概率服從泊松模型（Dieterich and Kilgore，1996）：

那么，利用公式（8）的概率分布，結(jié)合公式（7），可以計算在Δτ的應(yīng)力擾動下造成大于M級地震的發(fā)震概率.

然而，類似于2013年蘆山地震和2014年于田地震這種情況，即在大震后不久又發(fā)生大地震時，就不能單一地用公式（8）計算，因為第一次地震之后由于應(yīng)力的調(diào)整，地震活動性會發(fā)生顯著變化，所以，對于區(qū)域內(nèi)多次地震的情況，應(yīng)考慮發(fā)震概率的變化.

3 發(fā)震概率計算

新疆于田2008年3月21日7.3級地震位于西昆侖地震帶康西瓦斷裂、阿爾金地震帶和東昆侖地震帶可可西里斷裂交匯區(qū)域，地質(zhì)構(gòu)造發(fā)育而復(fù)雜.震源斷錯類型以走滑型為主（高國英等，2005；王盛澤和高國英，1992；朱令人，2002）.2008年于田7.3級地震震源區(qū)及其附近歷史中強震主壓應(yīng)力軸方位以近NS向為主，震源區(qū)斷錯性質(zhì)以走滑為主.昆侖山地震帶和阿爾金地震帶具備發(fā)生7級以上地震的構(gòu)造條件，1900年以來先后發(fā)生了1924年民豐7.3級雙震，1973年瑪尼7.3級地震，1996年喀喇昆侖山7.1級地震，1997年瑪尼7.5級地震，以及2001年昆侖山口西8.1級地震等5組7級以上地震（王瓊等，2009）.2008年于田7.3級地震后，至2013年11月，已有M＞6.0的地震3個，M＞5.0的地震23個，M＞4.0的地震193個，M＞3.0的地震383個（中國地震臺網(wǎng)中心）.圖2顯示了2008年于田7.3級地震前（1970年1月1日至2008年3月20日）后（2008年3月21日至2013年11月30日）的地震分布以及M＞3.0地震的年發(fā)生率.震后于田周邊地區(qū)的應(yīng)力狀態(tài)發(fā)生了明顯的變化（王瓊等，2009；萬永革等，2010），從而導(dǎo)致了地震活動性分布的顯著不同，并影響了未發(fā)震斷層的發(fā)震概率.

3.1 余震持續(xù)時間估算

圖2 2008年于田地震前（a）后（b）的M＞3.0地震分布以及地震年發(fā)生率；網(wǎng)格大小為0.02°×0.02°，右側(cè)色標顯示的年發(fā)生率為網(wǎng)格內(nèi)M＞3.0的年發(fā)生率.紅色圓圈代表M＞5.0的地震，黑色圓圈代表3.0＜M≤5.0的地震，灰白色實線代表斷層的位置Fig.2 Seismicity rate before（a）and after（b）Yutian earthquake（M＞3.0），represented by red（M＞5.0）and black circles（3.0＜M≤5.0）；the seismicity on colorbar on the right is in the grid of 0.02°×0.02°.The gray line is the fault.

Toda等（1998）在對1995年日本阪神地震余震的研究中，采用對數(shù)坐標線性擬合單位時間內(nèi)余震個數(shù)隨時間的變化關(guān)系，并以擬合線回到背景地震活動性所用時間作為余震持續(xù)時間.我們統(tǒng)計得到2008年于田地震后單位時間內(nèi)（月）余震數(shù)目，并進行對數(shù)線性擬合，當余震發(fā)生率回歸到背景場水平時則認為余震序列不再衰減，而回歸所需時間則認為是余震序列持續(xù)時間.考慮到早期背景場地震數(shù)據(jù)中小震數(shù)目不足以及中小震數(shù)據(jù)嚴重缺失，我們選取從1987年開始的M＞3.0和M＞4.0地震數(shù)據(jù)進行分析.從1987年1月1日到2008年3月20日，共發(fā)生M＞3.0地震311次，M＞4.0地震195次，背景場M＞3.0地震月發(fā)生率為1.22個，M＞4.0地震月發(fā)生率為0.76個.

圖3中圓圈是中國地震臺網(wǎng)中心的地震目錄統(tǒng)計得到的2008年于田地震后的地震月發(fā)生率，實線是用對數(shù)線性擬合得到的地震月發(fā)生率，虛線為背景場地震月發(fā)生率，實線與虛線的交點橫坐標即為地震月發(fā)生率回歸到背景場所需要的時間，即余震序列持續(xù)時間.

圖3 2008年于田地震后單位時間內(nèi)余震數(shù)目隨時間衰減關(guān)系（a）M＞3.0，（b）M＞4.0.圓圈是實際地震數(shù)據(jù)統(tǒng)計得到的地震月發(fā)生率，實線是用對數(shù)線性擬合得到的地震發(fā)生率，虛線為背景場地震月發(fā)生率.Fig.3 Estimation of the aftershock duration ta，defined as the time until the seismicity rate returns to the background rate before 2008Yutian earthquakeThe attenuation of earthquake rate by month for（a）M＞3.0and（b）M＞4.0.The background rate of seismicity is shown as the dashed line.The circles represent the observed rate of earthquakes during the aftershock sequence.

從圖3中可見，主震后地震發(fā)生率逐月衰減，擬合曲線與衰減趨勢具有一定的符合度，基本反映出地震發(fā)生率的衰減趨勢.M＞3.0地震月發(fā)生率回歸到背景場水平所需時間為15.2個月，M＞4.0地震發(fā)生率回歸到背景場水平所需時間為14.9個月，綜合兩者結(jié)果，用對數(shù)線性擬合得到的余震序列持續(xù)時間約為15個月左右.

3.2 康西瓦斷裂東段的發(fā)震概率計算

阿爾金斷裂帶的滑動速率，由于它與青藏高原形成演化模型之間的密切關(guān)系而備受矚目.王峰等（2004）通過對阿爾金斷裂帶西段的三個典型斷錯地貌點進行研究，得到了阿爾金斷裂帶西段的滑動速率為11.4±2.5mm·a－1.2008年于田地震在地表破裂帶測量到的最大左旋走滑位移Dm＝1.8m（徐錫偉等，2011），震源深度12km（王凡等，2011），從而可以推斷此處發(fā)生7級以上地震的復(fù)發(fā)周期tr約為129～202年.假設(shè)斷層半寬度a＝10km，那么應(yīng)力降Δσ＝2μDm／（3a）＝3.6MPa.假設(shè)區(qū)域內(nèi)發(fā)生7級以上的平均復(fù)發(fā)周期為160年，那么M＞7.0的地震在該地區(qū)的年發(fā)生率根據(jù)背景場地震活動性速率公式r＞M＝N／Δt（其中N為發(fā)生大于震級M的地震個數(shù)，Δt為時間間隔）可以求出：r＞7.0＝0.0056個／年，設(shè)A＝0.01，σN＝10MPa（Dieterich and Kilgore，1996），余震持續(xù)時間ta＝15個月，在不同應(yīng)力擾動的情況下，根據(jù)公式（6）和（7）計算10年內(nèi)發(fā)生大于設(shè)定震級M7.0的概率如圖4所示.從圖4中可見，大應(yīng)力擾動可導(dǎo)致地震發(fā)生率的快速上升，因此受到大應(yīng)力擾動的地區(qū)地震發(fā)生的時間會大幅提前，黑色實線為應(yīng)力擾動造成的發(fā)震概率上限.根據(jù)萬永革等（2010）對2008年于田地震主斷層周圍庫侖破裂應(yīng)力變化的研究結(jié)果，康西瓦斷裂東段庫侖破裂應(yīng)力增加了約0.8MPa，而2014年于田地震的發(fā)震斷層恰為阿爾金斷裂帶西端部馬尾狀分支斷裂康西瓦斷裂東段上.假設(shè)區(qū)域內(nèi)各條斷層7.0級以上地震的復(fù)發(fā)周期相同，那么各條斷層從某一時刻開始的發(fā)震概率變化也相同.根據(jù)圖5中黑色和紅色曲線，受到2008年于田地震的影響，康西瓦斷裂東段在2014年發(fā)生7.0級以上的概率從3.3%增加到5.83%，康西瓦斷裂東段6.5級以上的發(fā)震概率從23.9%增加到38.6%，6.0級以上的發(fā)震概率從42.0%增加到62.3%.

圖4 在不同應(yīng)力擾動下，10年內(nèi)發(fā)生大于設(shè)定震級M7.0的概率變化Fig.4 The probability of earthquake occurrence of M＞7.0 of different shear stress disturbance in 10years

如果2014年康西瓦斷裂未發(fā)生地震，那么10年后（到2024年）康西瓦斷裂東段發(fā)生6.0級以上、6.5級以上和7.0級以上地震的概率分別應(yīng)為84.8%、61.0%和11%，如圖5藍色曲線所示；由于康西瓦斷裂東段在2014年發(fā)生失穩(wěn)，所以斷層又重新進入了下一個孕震階段，6.0級以上、6.5級以上和7.0級以上地震的發(fā)震概率迅速減小，而后隨時間推延又逐漸增長，到2024年，康西瓦斷裂東段的發(fā)震概率增加到了63.9%、39.9%和6.09%.同樣，20年后（到2034年）康西瓦斷裂東段的發(fā)震概率也分別從93.9%、75.3%和15.8%變化到85.5%、61.9%和11.2%.從圖5中不難看出，無論是2008年于田地震導(dǎo)致的康西瓦斷裂東段的發(fā)震概率增加還是2014年康西瓦斷裂東段的失穩(wěn)導(dǎo)致的未來幾十年的發(fā)震概率的減小，6.0級以上和6.5級以上的發(fā)震概率的變化都比7.0級以上的變化顯著.而在未來20年以后，6.5級以上的地震發(fā)震概率在2008年和2014年震前震后的差別比6.0級以上的地震更加顯著.7.0級以上地震的發(fā)震概率在2014年震后以近似線性方式增長.

3.3 周邊斷層發(fā)震概率計算

圖5 康西瓦斷裂東段50年內(nèi)發(fā)生大于不同設(shè)定震級的概率變化橫坐標是從2008年于田地震開始計算，所以三條淺藍色豎線分別代表2014年、2024年和2034年.黑色曲線代表2008年于田地震庫侖應(yīng)力變化后的發(fā)震概率，紅色曲線代表無應(yīng)力擾動下的發(fā)震概率，藍色曲線代表2014年發(fā)震后的發(fā)震概率.Fig.5 The probability of earthquake occurrence of Kengxiwar fault in 50years The X－coordinate is the time from 2008Yutian earthquake.The light blue vertical lines point out the year of 2014，2024and 2034.The black curves are the probability of earthquake after the stress perturbations in 2008；the red curves are the probability of earthquake without any perturbations；the blue curves are the probability of earthquake after 2014.

事實上，2008年于田地震之后不僅康西瓦斷裂東段的庫侖應(yīng)力增加，其周邊的斷層也有庫侖應(yīng)力的明顯變化（萬永革等，2010），庫侖應(yīng)力在貢嘎錯斷裂的西南段（GGS）和東北段（GGN）、康西瓦斷裂西段（KXW）、普魯斷裂西段（PLW）、龍木錯—邦達錯西段（LBW）、阿爾金斷裂（ATF）和西昆侖山前逆沖斷裂（WKS）的增加分別達6.9×104Pa、4.4×103Pa、2.8×103Pa、2.3×104Pa、5.7×103Pa、1.6×103Pa和1.1×103Pa，而在貢嘎錯斷裂中段（GGC）、康西瓦斷裂中段（KXC）、龍木錯—邦達錯東段（LBE）、普魯斷裂中段（PLC）和喀喇昆侖北段（KKN）斷層上的庫侖破裂應(yīng)力減小分別達3.9×105Pa、5.8×105Pa、1.5×104Pa、1.5×104Pa和1.1×103Pa，斷層位置如圖6所示（萬永革等，2010）.對于這些斷層而言，除了受2008年于田地震的庫侖應(yīng)力影響外，在2014年于田地震發(fā)生后，還會受到庫侖應(yīng)力變化的再次影響.中國地震局地震預(yù)測研究所給出了2014年2月12日新疆于田地震震源和庫侖應(yīng)力觸發(fā)初步研究結(jié)果：貢嘎錯斷裂中段（GGC）和西南段（GGS）受到了小于0.6MPa的庫侖應(yīng)力增加，貢嘎錯斷裂北段（GGN）和阿爾金斷裂（ATF）受到了約0.01MPa的庫侖應(yīng)力增加，康西瓦斷裂中段（KXC）受到了小于0.6MPa的庫侖應(yīng)力減少.上述五條斷層在兩次于田地震后的庫侖應(yīng)力變化見表1.

表1 五條斷層庫侖應(yīng)力變化Table 1 The Coulomb stress change of the five faults

圖6 新疆于田地區(qū)主要斷裂分布黑線為該區(qū)域的主要活動斷層，斷層名以縮寫表示.Fig.6 Regional fault map in Yutian area Black lines are the major active faults with the name described by abbreviations.

由于貢嘎錯斷裂北段（GGN）和阿爾金斷裂（ATF）距離較遠，受到這兩次于田地震的庫侖應(yīng)力變化不大，在此不予討論，下面對貢嘎錯斷裂中段（GGC）、南段（GGS）和康西瓦斷裂中段（KXC）受到兩次應(yīng)力擾動后的發(fā)震概率變化進行分析.

貢嘎錯斷裂中段（GGC）在2008年于田地震后，庫侖應(yīng)力減小了0.39MPa，使得發(fā)震概率增幅大幅度減緩，發(fā)生7.0級以上地震的概率為0.14%，發(fā)生6.5級以上地震的概率下降到1.14%，發(fā)生6.0級以上地震的概率下降到2.28%（見圖7）；直到2014年于田再次地震后，該斷層受到小于0.6MPa的庫侖應(yīng)力增加，發(fā)震概率增幅加速，未來50年的發(fā)震概率由表2給出.值得注意的是，表2和圖7中給出的是按照0.6MPa庫侖應(yīng)力增加計算的值，所以應(yīng)該為該斷層未來50年的發(fā)震概率上限.受到兩次于田地震的影響，貢嘎錯斷裂中段（GGC）7.0級以上震級的發(fā)震概率超越95%需要530年；6.5級以上震級的發(fā)震概率超越95%需要62年，而6.0級以上震級的發(fā)震概率超越95%只需要約29年，見圖7.

表2 2014年于田地震后貢嘎錯斷裂中段（GGC）發(fā)生中強地震的概率Table 2 The earthquake probability of the central segment of Gonggacuo fault（GGC）in 50years after 2014Yutian earthquake

2008年于田地震后，在貢嘎錯斷裂西南段（GGS）庫侖應(yīng)力增加了0.069MPa，使得發(fā)震概率增幅加快，發(fā)生7.0級以上地震的概率增加到4.60%，發(fā)生6.5級和6.0級以上地震的概率分別增加到31.78%和53.47%（圖8）；2014年于田再次地震后，該斷層又受到小于0.6MPa的庫侖應(yīng)力增加，發(fā)震概率增幅突然加速，未來50年的發(fā)震概率由表3給出.同樣，表3和圖8中給出的是按照0.6MPa庫侖應(yīng)力增加計算的值，所以也應(yīng)該為該斷層未來50年的發(fā)震概率上限.受到兩次于田地震的影響，若今后周邊無庫侖應(yīng)力擾動，貢嘎錯斷裂西南段（GGS）7.0級、6.5級和6.0級地震的發(fā)震概率超越95%分別需要494年、26年和0.8年，見圖8.

表3 2014年于田地震后貢嘎錯斷裂西南段（GGS）發(fā)生中強地震的概率Table 3 The earthquake probability of the southwest segment of Gonggacuo fault（GGS）in 50years after 2014Yutian earthquake

圖7 貢嘎錯斷裂中段（GGC）在2008年于田地震后50年內(nèi)發(fā)生大于不同設(shè)定震級的概率變化黑色曲線代表經(jīng)歷了兩次于田地震庫侖應(yīng)力變化的發(fā)震概率，紅色曲線代表無應(yīng)力擾動下的發(fā)震概率.黑色豎線指示貢嘎錯斷裂中段（GGC）2014年于田地震發(fā)震時，不同設(shè)定震級的發(fā)震概率變化.Fig.7 The probability of earthquake occurrence of the central segment of Gonggacuo fault（GGC）in 50years after 2008Yutian earthquake The black curves are the probabilities of earthquake with the stress perturbations in 2008and 2014，while the red curves are the probability of earthquake without any perturbations.The black vertical line points out the earthquake probability change when the M7.3Yutian earthquake happened on Feb 12th，2014.

2008年于田地震后，在康西瓦斷裂中段（KXC）庫侖應(yīng)力降低了0.58MPa，使得發(fā)震概率增幅幾乎停滯，發(fā)生7.0級以上地震的概率降低到0.02%，發(fā)生6.5級以上地震的概率下降到0.18%，發(fā)生6.0級以上地震的概率下降到0.35%（見圖9）；2014年于田再次地震后，該斷層又受到小于0.6MPa的庫侖應(yīng)力下降，發(fā)震概率增幅繼續(xù)減緩，未來50年的發(fā)震概率由表4給出.表4和圖9中給出的是按照0.6MPa庫侖應(yīng)力降低計算的值，所以也應(yīng)該為該斷層未來50年的發(fā)震概率下限.受到兩次于田地震的影響，若之后周邊無庫侖應(yīng)力擾動，康西瓦斷裂中段（KXC）相當不活躍，7.0級、6.5級和6.0級地震的發(fā)震概率超越95%分別需要556年、87年和54年，見圖9.

表4 2014年于田地震后康西瓦斷裂中段（KXC）發(fā)生中強地震的概率Table 4 The earthquake probability of the central segment of Kengxiwar fault（KXC）in 50years after 2014Yutian earthquake

圖8 同圖7，但為貢嘎錯斷裂西南段（GGS）在2008年于田地震后50年內(nèi)發(fā)生大于不同設(shè)定震級的概率變化Fig.8 Same as Fig.7，but for the probability of earthquake occurrence of the southwest segment of Gonggacuo fault（GGS）in 50years after 2008Yutian earthquake

圖9 同圖7，但為康西瓦斷裂中段（KXC）在2008年于田地震后50年內(nèi)發(fā)生大于不同設(shè)定震級的概率變化Fig.9 Same as Fig.7，but for the probability of earthquake occurrence of the central segment of Kengxiwar fault（KXC）in 50years after 2008Yutian earthquake

4 討論與結(jié)論

蘆山M7.0強地震的發(fā)生重新引起了對汶川地震的熱議，而2014年MS7.3于田地震再一次為我們敲響了警鐘.大震后周圍應(yīng)力場變化導(dǎo)致的區(qū)域地震活動性的重新估計成為了震后庫侖應(yīng)力計算的重要方面.本文利用滑移速率和狀態(tài)相依賴的摩擦定律（Dieterich，1994），結(jié)合2008年于田地震前后的地震活動性水平，從地震發(fā)生概率的變化上解釋了2014年于田地震的發(fā)震根源.康西瓦斷裂東段在2014年發(fā)生7.0級以上的概率從3.3%增加到5.83%.康西瓦斷裂東段受到2008年于田地震主斷層庫侖破裂應(yīng)力變化的影響，6.5級以上震級的發(fā)震概率從23.9%增加到38.6%，6.0級以上震級的發(fā)震概率從42.0%增加到62.3%.而后，又對2014年康西瓦斷裂進行未發(fā)震假設(shè)，對10年后（到2024年）和20年后（到2034年）康西瓦斷裂東段發(fā)生6.0級以上、6.5級以上和7.0級以上地震的概率與其在2014年失穩(wěn)后的概率進行比較，發(fā)現(xiàn)無論是2008年于田地震導(dǎo)致的康西瓦斷裂東段的發(fā)震概率增加還是2014年康西瓦斷裂東段的失穩(wěn)導(dǎo)致的未來幾十年的發(fā)震概率的減小，6.0級以上和6.5級以上的發(fā)震概率的變化都比7.0級以上的變化顯著.而在未來20年以后，6.5級以上的地震發(fā)震概率在2008年和2014年震前震后的差別比6.0級以上的地震的更加顯著.

此外，本文對周邊受到兩次于田地震影響較大的三條斷層也進行了發(fā)震概率的計算.經(jīng)歷了兩次于田地震的區(qū)域庫侖應(yīng)力狀態(tài)變化，貢嘎錯斷裂中段（GGC）、貢嘎錯斷裂西南段（GGS）和康西瓦斷裂中段（KXC）分別經(jīng)歷了發(fā)震概率先降后升、先升后升和先降后降兩個階段，充分顯示了庫侖應(yīng)力的細微變化造成的周邊斷層的危險性的變化.這三條斷裂發(fā)生7.0級以上地震的發(fā)震概率超越95%均需要500年左右；貢嘎錯斷裂西南段（GGS）發(fā)生中強地震的可能性較大，6.5級和6.0級地震的發(fā)震概率超越95%分別需要26年和0.8年；而康西瓦斷裂中段（KXC）活躍度較低，6.5級和6.0級地震的發(fā)震概率超越95%分別需要87年和54年.

然而，在計算過程中，存在很多不確定性.首先，參數(shù)的選取可能帶來一定誤差，比如A取0.01，σN取10MPa（Dieterich and Kilgore，1996）等，與實際不完全一致；其次，庫侖應(yīng)力狀態(tài)的計算過程依賴模型的選取，本文僅選用萬永革等（2010）和中國地震局地震預(yù)測研究所的結(jié)果進行計算也會帶來一定的誤差；第三，地震目錄特別是小震目錄的缺失使得對背景場地震活動性的估計不夠準確，從而影響余震持續(xù)時間的估計；此外，在本文對庫侖應(yīng)力的計算過程中未考慮數(shù)十年的較大時間尺度內(nèi)介質(zhì)黏彈效應(yīng)的影響，這對于庫侖應(yīng)力的估計也會有一定的影響.所以，應(yīng)用條件概率模型會更具可靠性.但是由于缺乏發(fā)震斷層的滑移歷史、發(fā)震的準周期以及上次地震至今的時間，條件概率模型在目前的研究中，尤其對我國西南地區(qū)來講，仍然難以應(yīng)用.

Chen Y T，Yang Z X，Zhang Y，et al.2013.From 2008Wenchuan earthquake to 2013Lushan earthquake.ScientiaSinicaTerrae（in Chinese），43（6）：1064－1072.

Das S，Scholz C H.1981.Theory of time－dependent rupture in the Earth.J.Geophys.Res.，86（B7）：6039－6051.

Dieterich J H.1972. Time－dependent friction as a possible mechanism for aftershocks.J.Geophys.Res.，77（20）：3771－3781.

Dieterich J H.1979.Modeling of rock friction：1.Experimental results and constitutive equations.J.Geophys.Res.，84（B5）：2161－2168.

Dieterich J H.1981.Constitutive properties of faults with simulated gouge.／／ Carter N L，F(xiàn)riedman M，Logan J M eds.Mechanical Behavior of Crustal Rocks：The Handin Volume.Geophysical Monograph Series，Volume 24.Washington D.C.：AGU，103－120.

Dieterich J H.1994.A constitutive law for rate of earthquake production and its application to earthquake clustering.J.Geophys.Res.，99（B2）：2601－2618.

Dieterich J H，Kilgore B D.1996.Implications of fault constitutive properties for earthquake prediction.Proc.Natl.Acad.Sci.U.S.A.，93（9）：3787－3794.

Du F，Long F，Ruan X，et al.2013.TheM7.0Lushan earthquake and the relationship with theM8.0Wenchuan earthquake in Sichuan，China.ChineseJ.Geophys.（in Chinese），56（5）：1772－1783，doi：10.6038／cjg20130535.

Gao G Y，Wen H P，Nie X H.2005.Analysis of focal－mechanism solution of moderately strong earthquakes in Xinjiang from 1991to 2002.Earthquake（in Chinese），25（1）：81－87.

Gomberg J，Blampied M L，Beeler N M.1997.Transient triggering of near and distant earthquakes.Bull.Seismol.Soc.Am.，87（2）：294－309.

Gomberg J，Beeler N M，Blampied M L，et al.1998.Earthquake triggering by transient and static deformations.J.Geophys.Res.，103（B10）：24411－24426.

Gross S，Kisslinger C.1997.Estimating tectonic stress rate and state with Landers aftershocks.J.Geophys.Res.，102（B4）：7603－7612.

Gross S，Bürgmann R.1998.Rate and state of background stress estimated from the aftershocks of the 1989Loma Prieta，California，earthquake.J.Geophys.Res.，103（B3）：4915－4927.

Harris R A，Day S M.1993.Dynamics of fault interaction：Parallel strike－slip fault.J.Geophys.Res.，98（B3）：4461－4472.

Harris R A.1998.Introduction to special section：Stress triggers，stress shadows，and implications for seismic hazard.J.Geophys.Res.，103（B10）：24347－24358.

Harris R A，Simpson R W.1998.Suppression of large earthquakes by stress shadows：a comparison of Coulomb and rate－and－state failure.J.Geophys.Res.，103（B10）：24439－24451.

Kanamori H，Brodsky E E.2004.The physics of earthquakes.Rep.Prog.Phys.，67（8）：1429－1496.

King G C P，Stein R S，Lin J.1994.Static stress changes and the triggering of earthquakes.Bull.Seismol.Soc.Am.，84（3）：935－953.

Lei X L，Ma S L，Su J R，et al.2013.Inelastic triggering of the 2013Mw6.6Lushan earthquake by the 2008Mw7.9Wenchuan earthquake.SeismologyandGeology（in Chinese），35（2）：411－422.

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.

Parsons T，Ji C，Kirby E.2008.Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin.Nature，454（7203）：509－510.

Reasenberg P A，Simpson R W.1992.Response of regional seismicity to the static stress change produced by the Loma Prieta earthquake.Science，255（5052）：1687－1690.

Rice J R，Gu J C.1983.Earthquake aftereffects and triggered seismic phenomena.PureAppl.Geophys.，121（2）：187－219.Ruina A L.1983.Slip instability and state variable friction laws.J.Geophys.Res.，88（B12）：10359－10370.

Rybicki K.1973.Analysis of aftershocks on the basis of dislocation theory.Phys.EarthPlanet.Inter.，7（4）：409－422.

Shan B，Xiong X，Zheng Y，et al.2013.Stress changes on major faults caused by 2013Lushan earthquake，and its relationship with 2008Wenchuan earthquake.ScienceChina：Earth Sciences（in Chinese），43（6）：1002－1009.

Simpson R W，Reasenberg P A.1994.Earthquake－induced static stress changes on central California faults.U.S.Geol.Surv.Prof.Pap.，1550－F：F55－F89.

Spudich P，Steck L K，Hellweg M，et al.1995.Transient stresses at Parkfield，California，produced by theM7.4Landers earthquake of June 28，1992：Observations from UPSAR dense seismograph array.J.Geophys.Res.，100（B1）：675－690.

Stein R S，King G C P，Lin J.1992.Change in failure stress on the southern San Andreas fault system caused by the 1992 magnitude＝7.4Landers earthquake.Science，258（5086）：1328－1332.

Stein R S，King G C P，Lin J.1994.Stress triggering of the 1994M＝6.7Northridge，California，earthquake by its predecessors.Science，265（5177）：1432－1435.

Toda S，Stein R S，Reasenberg P A，et al.1998.Stress transferred by the 1995MW＝6.9Kobe，Japan，shock：Effect on aftershocks and future earthquakes probabilities.J.Geophys.Res.，103（B10）：24543－24565.

Wan Y G，Shen Z K，Sheng S Z，et al.2010.The mechanical effects of the 2008MS7.3Yutian，Xinjiang earthquake on the neighboring faults and its tectonic origin of normal faulting mechanism.ChineseJ.Geophys.（in Chinese），53（2）：280－289，doi：10.3969／j.issn.0001－5733.2010.02.006.

Wang F，Xu X W，Zheng R Z，et al.2004.Late quaternary sliprate on the Altun fault west to the Qarqan river.Seismology andGeology（in Chinese），26（2）：200－208.

Wang F，Chen W T，Wang M，et al.2011.Coseismic and postseismic deformation ofMS7.3Yutian earthquake derived from GPS data.ChineseJ.Geophys.（in Chinese），54（9）：2250－2255，doi：10.3969／j.issn.0001－5733.2011.09.007.

Wang Q，Nie X H，Wen H P.2009.Preliminary research on the characteristics of regional seismic activity and stress state before the 2008MS7.3earthquake in Yutian，Xinjiang of China.ActaSeismologicaSinica（in Chinese），31（3）：235－244.

Wang S Z，Gao G Y.1992.The regional characteristics of tectonic stress field in Xinjiang and its neighboring areas.Acta SeismologicaSinica（in Chinese），14（S）：612－620.

Xu X W，Tan X B，Wu G D，et al.2011.Surface rupture features of the 2008YutianMS7.3earthquake and its tectonic nature.SeismologyandGeology（in Chinese），33（2）：462－471.

Xu X W，Yu G H.2014.Earthquakes structural map ofM7.3in Yutian County，Hotan Prefecture，Xinjiang Uygur Autonomous Region（latitude 36.1°，longitude 82.5°）（in Chinese）.http：／／www.eq－igl.ac.cn［2014－02－16］.

Zhu L R.2002.Xinjiang Annals（Volumes 11，Earthquake）（in Chinese）.Urumqi：Xinjiang People′s Publishing House.

附中文參考文獻

陳運泰，楊智嫻，張勇等.2013.從汶川地震到蘆山地震.中國科學(xué)：地球科學(xué)，43（6）：1064－1072.

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

高國英，溫和平，聶曉紅.2005.1991—2002年新疆中強震震源機制解分析.地震，25（1）：81－87.

雷興林，馬勝利，蘇金蓉等.2013.汶川地震后中下地殼及上地幔的粘彈性效應(yīng)引起的應(yīng)力變化與蘆山地震的發(fā)生機制.地震地質(zhì)，35（2）：411－422.

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

單斌，熊熊，鄭勇等.2013.2013年蘆山地震導(dǎo)致的周邊斷層應(yīng)力變化及其與2008年汶川地震的關(guān)系.中國科學(xué)：地球科學(xué)，43（6）：1002－1009.

萬永革，沈正康，盛書中等.2010.2008年新疆于田7.3級地震對周圍斷層的影響及其正斷層機制的區(qū)域構(gòu)造解釋.地球物理學(xué)報，53（2）：280－289，doi：10.3969／j.issn.0001－5733.2010.02.006.

王峰，徐錫偉，鄭榮章等.2004.阿爾金斷裂帶西段車爾臣河以西晚第四紀以來的滑動速率研究.地震地質(zhì)，26（2）：200－208.

王凡，陳為濤，王敏等.2011.利用GPS觀測資料分析2008年于田MS7.3地震的同震位移及震后形變.地球物理學(xué)報，54（9）：2250－2255，doi：10.3969／j.issn.0001－5733.2011.09.007.

王瓊，聶曉紅，溫和平.2009.2008年新疆于田MS7.3地震前區(qū)域地震活動和應(yīng)力狀態(tài)特征初步研究.地震學(xué)報，31（3）：235－244.

王盛澤，高國英.1992.新疆及其鄰近地區(qū)現(xiàn)代構(gòu)造應(yīng)力場的區(qū)域特征.地震學(xué)報，14（增刊）：612－620.

徐錫偉，譚錫斌，吳國棟等.2011.2008年于田MS7.3地震地表破裂帶特征及其構(gòu)造屬性討論.地震地質(zhì)，33（2）：462－471.

徐錫偉，于貴華.2014.新疆維吾爾自治區(qū)和田地區(qū)于田縣（北緯36.1°，東經(jīng)82.5°）7.3級地震發(fā)震構(gòu)造圖.http：／／www.eqigl.ac.cn［2014－02－16］.

朱令人.2002.新疆通志（十一卷，地震志）.烏魯木齊：新疆人民出版社.