熊 維 譚 凱 劉 剛 王 偉

1)中國地震局地震研究所(地震大地測量重點(diǎn)實驗室)，武漢 430071

2)中國地震局地殼應(yīng)力研究所武漢科技創(chuàng)新基地，武漢430071

地震是應(yīng)力積累到一定程度突然釋放的過程，而應(yīng)力并不會完全消失，而是通過傳輸轉(zhuǎn)移到其他地區(qū)，導(dǎo)致應(yīng)力聚集進(jìn)而誘發(fā)后續(xù)地震。這即是地震的應(yīng)力觸發(fā)理論［1－2］。

國內(nèi)外學(xué)者對于地震應(yīng)力觸發(fā)方面的研究已經(jīng)有許多。Stein［3－4］在對圣安地列斯斷層的歷史地震資料和1992年發(fā)生在加州的蘭德斯地震的研究中認(rèn)為，地震之間的相互作用對地震危險性會產(chǎn)生很大影響，并認(rèn)為蘭德斯地震升高了圣安地列斯斷層南部的地震發(fā)生概率，降低了其北部(洛杉磯附近)的發(fā)震概率［4］。Parsons［5－6］計算了1999年伊茲米特地震對伊斯坦布爾附近斷層的作用，認(rèn)為該地震增加了伊斯坦布爾未來30 a 內(nèi)發(fā)生強(qiáng)震的概率。在汶川地震應(yīng)力觸發(fā)的研究方面，Parsons［7］對汶川地震產(chǎn)生的應(yīng)力變化及其對周邊斷層的影響作了初步計算。Toda［8］通過對汶川地震釋放應(yīng)力的計算認(rèn)為，在未來10 a，鮮水河斷裂、東昆侖斷裂和岷江斷裂的地震發(fā)生概率升高了1 倍。但是因為地質(zhì)背景信息和震源模型的不足，Toda 的結(jié)果存在一定的不確定性［1］。萬永革［9］、單斌［10］和解朝娣等［1］也計算了汶川地震導(dǎo)致的周邊斷層應(yīng)力變化。這些研究大部分集中于庫侖應(yīng)力的變化及分布，并未定量分析斷層上的發(fā)震概率變化。

汶川地震發(fā)生后，我們利用近場506 個大地測量點(diǎn)位及InSAR 資料反演了此次地震的滑動分布及斷層的幾何形態(tài)，得到了精細(xì)的同震破裂模型［11］。本文以此破裂模型為基礎(chǔ)，計算周邊斷層的庫侖應(yīng)力變化，得到地震危險性變化明顯的區(qū)域，并結(jié)合斷層發(fā)震概率的計算，定量地給出了研究區(qū)斷層在未來特定時段的發(fā)震概率。

1 原理與方法

1.1 背景地震發(fā)生概率和地震活動率的計算

在不考慮應(yīng)力擾動的條件下，Ward［12］提出GPS速度場計算地震矩率進(jìn)而計算地震發(fā)生概率的方法。根據(jù)此方法，本文利用汶川地震發(fā)生前在川滇地區(qū)的GPS 觀測數(shù)據(jù)試算了部分活動斷裂在特定時段的地震發(fā)生概率(表1)，其結(jié)果與張效亮等［13］計算結(jié)果近似。表1 中另外加入了參考文獻(xiàn)中用滑動速率方法計算的地震復(fù)發(fā)間隔作為計算結(jié)果的驗證。

地震發(fā)生率轉(zhuǎn)換為標(biāo)準(zhǔn)泊松概率的公式為［8］:

將用泊松概率P表示為:

若此時P 為未發(fā)生應(yīng)力擾動時的斷層發(fā)震概率，則r 即為背景地震活動率。

1.2 庫侖應(yīng)力的計算

使巖石破裂失穩(wěn)的庫侖應(yīng)力為:

式中，Δτ 為斷層上的剪切應(yīng)力變化(沿斷層滑動方向為正)，Δσn為正應(yīng)力變化(使斷層解鎖為正)。ΔP 為斷層區(qū)的孔隙壓力變化(壓縮為正)，μ 為摩擦系數(shù)(范圍為0 ～1)。破裂在Δσf為正時被驅(qū)使，為負(fù)時被抑制;增加剪切應(yīng)力和使斷層解鎖這兩種情況均可驅(qū)使破裂［3］。ΔP 與Δσn作用的趨勢相反，經(jīng)常在式(3)中引入“有效”摩擦系數(shù)μ'而合并成一項［14］，即

1.3 庫侖應(yīng)力變化對地震發(fā)生概率影響的計算

Dieterich［15］推導(dǎo)了應(yīng)力擾動后的地震活動率關(guān)于時間t 的函數(shù)R(t):

其中，A是斷層結(jié)構(gòu)參數(shù)，σ是全部正應(yīng)力，實際應(yīng)用中通常取Aσ 為0.1 ～0.5 bar，在此參考Toda等［5］的研究取為0.4 bar。r是背景地震活動率，Δσf是庫侖應(yīng)力變化，τa是應(yīng)力擾動的持續(xù)時間。在此參考解朝娣等［1，15］的研究，取值為290.5 a。要計算R(t)，首先要得到背景地震活動率和庫侖應(yīng)力變化。前者由式(2)得到，后者由式(4)得到。

令發(fā)生應(yīng)力擾動時的斷層發(fā)震概率為P'，則由式(1)、(2)、(5)有:

1.4 庫侖應(yīng)力計算中參數(shù)的選取

庫侖應(yīng)力的計算主要需要確定有效摩擦系數(shù)、接收斷層傾角、滑動角以及研究深度等參數(shù)。

由于幾乎不可能通過實際觀測來測定斷層的有效摩擦系數(shù)，所以目前主要是依據(jù)經(jīng)驗選?。?6］。有效摩擦系數(shù)的選取主要參考單斌等［10］的研究，而接收斷層的傾角和滑動角主要參考鄧起東等［17］。研究區(qū)有效摩擦系數(shù)、接收斷層傾角和滑動角的選取見表2。另外，本文計算的庫侖應(yīng)力變化面深度為10 km。

表1 川滇地區(qū)部分活動斷裂在汶川地震前的發(fā)震概率Tab.1 The earthquake probability of main active faults in Chuandian region before Wenchuan earthquake

圖1 庫侖應(yīng)力分布Fig.1 Distribution of Coulomb stress in the faults around

2 計算結(jié)果與分析

2.1 周邊斷層上的應(yīng)力變化分析

采用Coulomb3.3 軟件對汶川地震周邊斷層的庫侖應(yīng)力進(jìn)行計算，得到的結(jié)果(圖1、表3)與單斌等［10］、萬永革等［9］的結(jié)果接近，其差異主要來源于震源模型的選擇。單斌等［10］主要采用的是Ji 和Hayes 以及Nishimura 和Yagi 提供的震源破裂模型，這些模型主要源于遠(yuǎn)場地震波的數(shù)據(jù)反演［10］，并沒有加入足夠的野外考察結(jié)果。相比之下，加入了更多大地測量觀測數(shù)據(jù)作為約束條件的模型對于近場形變擬合效果更優(yōu)，更加適用于計算周邊斷層的庫侖應(yīng)力變化。

單斌等［10］認(rèn)為，鮮水河斷裂道孚-康定段大部分應(yīng)力增量在0.1 bar 以上，東昆侖斷裂、青川斷裂等處應(yīng)力增加，虎牙斷裂、撫邊河斷裂等處應(yīng)力降低。萬永革等［9］認(rèn)為，鮮水河斷裂庫侖應(yīng)力增加達(dá)0.01 ～0.6 bar，東昆侖斷裂、龍門山斷裂地震斷層兩端應(yīng)力均明顯增加。我們同樣發(fā)現(xiàn)，鮮水河斷裂道孚-康定段有明顯的應(yīng)力增加現(xiàn)象，其中最大應(yīng)力增加值達(dá)到0.67 bar。雖然該應(yīng)力變化峰值有可能是引入粗差導(dǎo)致的計算錯誤，但是次最大值也達(dá)到了0.29 bar，這說明鮮水河斷裂道孚-康定段受到汶川地震釋放的應(yīng)力影響很大。鮮水河斷裂活動頻繁、地震多發(fā)，自公元1747年以來共發(fā)生破壞性地震10 余次。受到汶川地震影響后，該地區(qū)地震發(fā)生概率大幅提升。Toda 等［8］認(rèn)為，鮮水河斷裂從上一次地震至今，已經(jīng)積累了1.1 ～2.3 m 的滑動不足量?？梢哉J(rèn)為，受汶川地震應(yīng)力影響，該斷層更加接近于破裂。

表2 研究區(qū)主要活動斷裂計算參數(shù)選取Tab.2 The parameters of main active faults in the area

表3 汶川地震導(dǎo)致的周邊斷層靜態(tài)庫侖應(yīng)力變化Tab.3 Coulomb stress change in the faults around caused by Wenchuan Earthquake

其他應(yīng)力增加的斷層有東昆侖斷裂東段、嘉黎察隅斷裂、岷江斷裂南段、馬爾康斷裂南段、成都德陽斷裂東北段等。東昆侖斷裂曾發(fā)生過1997年Mw7.6 瑪尼地震和2001年Mw7.8 昆侖山西口地震等多次強(qiáng)震，而瑪曲-南坪段則是著名的地震空區(qū)，有長時間的應(yīng)力積累，有可能在未來發(fā)生強(qiáng)烈地震。岷江斷裂南端是地震易觸發(fā)區(qū)域，曾于1933年發(fā)生過M7.5地震［10］。而且岷江斷裂斷層參數(shù)的不確定程度較高，在未來需要結(jié)合更多的資料確定斷層參數(shù)，以獲得更加可靠的結(jié)果。

應(yīng)力積累明顯降低的區(qū)域有鮮水河斷裂康定-瀘定段、岷江斷裂北段、成都德陽斷裂西南段和虎牙斷裂。這些區(qū)域未來的地震發(fā)生概率將明顯降低。

其他斷層上的應(yīng)力變化不明顯，很多區(qū)域只有10－3bar 量級的變化，對于斷層發(fā)震概率的影響可以忽略不計。

2.2 周邊斷層發(fā)震概率的估計

鮮水河斷裂道孚-康定段是地震易觸發(fā)區(qū)域，且受汶川地震釋放應(yīng)力的影響很大，該區(qū)域在震前和震后特定時段的地震發(fā)生概率如表4、圖2 所示。

表4 鮮水河斷裂道孚-康定段的地震發(fā)生概率Tab.4 The earthquake probability in Daofu-Kangding segment of Xianshuihe fault

顯而易見，在汶川地震發(fā)震后，該區(qū)域的地震發(fā)生概率比震前有明顯升高。如未來10 a年內(nèi)7級地震的發(fā)生概率由震前的8.76% 升高到震后的15%，升高了近1 倍。

圖2 鮮水河斷裂道孚-康定段的地震發(fā)生概率Fig.2 Histogram showing the earthquake probability in Daofu-Kangding segment of Xianshuihe fault

考慮汶川地震的影響，在未來100 a 內(nèi)，該區(qū)域發(fā)生6.5級地震的概率達(dá)96%，發(fā)生7級地震的可能性也接近75%。該結(jié)果說明，鮮水河斷裂道孚-康定段已成為極易觸發(fā)地震的區(qū)域，此處發(fā)生8級強(qiáng)震的概率約8%。

另外，由式(3)可得到應(yīng)力擾動后地震發(fā)生概率隨時間變化的趨勢，取研究時限為震后1 000 a，結(jié)果如圖3。該結(jié)果與解朝娣［1］接近。

圖3 鮮水河斷裂道孚-康定段的地震發(fā)生概率隨時間的變化Fig.3 Changes of earthquake probability with time in Daofu-Kangding segment of Xianshuihe

3 結(jié) 論

在計算震源周邊地區(qū)最優(yōu)取向破裂面上的應(yīng)力變化中發(fā)現(xiàn)，其應(yīng)力空間分布與余震的空間分布區(qū)域大致吻合。這種計算方法有其局限性:首先，最優(yōu)取向破裂面的破裂方向并不一定與實際存在的斷層破裂方向一致，有時甚至相反，這種情況下相應(yīng)區(qū)域的庫侖應(yīng)力變化的計算結(jié)果就失去了意義;第二，主震是發(fā)生在斷層上的，而震源周邊存在很多不包含活動斷層的區(qū)域，在這些區(qū)域上的計算或許對余震的預(yù)報有一定的參考價值，但是對主震發(fā)生概率的研究意義并不大。相較之下，對特定接收斷層上的計算，由于考慮了具體斷層的走向、傾角、滑動角、有效摩擦系數(shù)等參數(shù)，其結(jié)果更加切合特定斷層的情況，因此更加可靠。

受汶川地震影響，庫侖應(yīng)力增加最明顯的區(qū)域是鮮水河斷裂道孚-康定段。Toda 等［8］認(rèn)為，鮮水河斷裂在上一次地震至今，已經(jīng)積累了1.1 ～2.3 m的滑動不足量。受汶川地震應(yīng)力影響，該斷層更加接近于破裂。

庫侖應(yīng)力結(jié)果受到多種因素的影響。等效摩擦系數(shù)的變化會改變正應(yīng)力在庫侖應(yīng)力中的權(quán)重，對最終結(jié)果產(chǎn)生線性影響;接收斷層參數(shù)(走向、傾角和滑動角)的變化會改變地震應(yīng)力張量在斷層破裂面上的投影結(jié)果，對庫侖應(yīng)力計算產(chǎn)生非線性影響。因此，必須結(jié)合各方面的資料，縮小各類參數(shù)的不確定范圍，減小參數(shù)不確定性對結(jié)果的影響。

參考Ward 等［12］的研究，并結(jié)合張效亮等［13］的工作，總結(jié)了由GPS 速度場計算研究區(qū)活動斷層地震發(fā)生概率的方法，計算所得地震發(fā)生概率可以作為研究區(qū)地震背景發(fā)生率的計算依據(jù)。

1 解朝娣.Ms8.0 汶川地震產(chǎn)生的應(yīng)力變化空間分布及其對地震活動性的影響［J］.中國科學(xué):地球科學(xué)，2010，40(6):655－698.(Xie Chaodi.Pattern of stress change and its effect on seismicity rate caused by Ms8.0 Wenchuan earthquake［J］.Sci China Earth Sci，2010，40(6):688－698)

2 Stein R S.Earthquake conversations［J］.Sci Am，2003，288:72－79.

3 Stein R S.The role of stress transfer in earthquake occurrence［J］.Nature，1999，402:605－609.

4 Stein R S，King G C P，Lin J.Change in failure stress on the southern San Andreas fault system caused by the 1992 Magnitude=7.4 Landers earthquake［J］.Science，1992，258:1 328－1 332.

5 Barka A，Toda S，Stein R S，et al.Influence of the 17 August 1999 Izmit earthquake on seismic hazards in Istanbul［C］.The 1999 Izmit and Duzce Earthquakes:Preliminary Results，Istanbul，2000.

6 Toda P S，Stein R S，Barka A，et al.Heightened odds of large earthquakes near Istanbul:An interaction-based probability calculation［J］.Science，2000，288:661－665.

7 Parsons T，Ji C，Kirby E.Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin［J］.Nature，2008，doi:1011038.

8 Toda S，Lin J，Meghraoui M，et al.1 2 May2 0 0 8 M= 7.9 Wenchuan，China，earthquake calculated to increase failure stress and seismicity rate on three major fault systems［J］.Geophys Res Lett，2008，35:L17305.

9 萬永革，沈正康，盛書中，等.2008年汶川大地震對周圍斷層的影響［J］.地震學(xué)報，2009，31:128－139.(Wan Yongge，Shen Zhengkang，Sheng Shuzhong，et al.The influence of 2008 Wenchuan earthquake on surrounding faults［J］.Acta Seismologica Sinica，2009，31(2):128－139)

10 單斌，熊熊，鄭勇，等.2008年5月12日Mw7.9 汶川地震導(dǎo)致的周邊斷層應(yīng)力變化［J］.中國科學(xué)D 輯:地球科學(xué)，2009，39:537－545.(Shan B，Xiong X，Zheng Y，et al.Stress changes on major faults caused by Mw7.9 Wenchuan earthquake，May12，2008［J］.Sci China Ser D-Earth Sci，2009，52(5):593－601)

11 Wang Qi，Qiao Xuejun，Lin J.Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmenshan［J］.Nature，2011.DOI:10.1038/NGEO1210.

12 Steven N Ward.A multidisciplinary approach to seismic hazard in southern California［J］.Bull Seismol Soc Am，1994，84:1 293－1 309.

13 張效亮，謝富仁.利用GPS 數(shù)據(jù)和實時概率模型評估川滇南部中長期地震危險性［J］.大地測量與地球動力學(xué)，2009，29(3):38－41.(Zhang Xiaoliang，Xie Furen.Middle and long-term seismic risk evaluation of southern Sichuan-Yunnan region based on GPS data and real-time probabilistic model［J］.J of Geodesy and Geodynamic，2009，29(3):38－41)

14 King G C P，Stein R S，Lin J.Static stress changes and the triggering of earthquakes［J］.Bull Seismol Soc Am，1994，84:935－953.

15 Dieterich J.A constitutive law for rate of earthquake production and its application to earthquake clustering［J］.J Geophys Res，1994，99:2 601－2 618.

16 Parsons T，Stein R S，Simpson R W，et al.Stress sensitivity of fault seismicity:a comparison between limited-offset oblique and major strike-slip faults［J］.J Geophys 2002 Res，1999，104:20 183－20 202.

17 鄧起東，張培震，冉勇康，等.中國活動構(gòu)造基本特征［J］.中國科學(xué)D 輯，2002，32(12):1 020－1 030.(Deng Qidong，Zhang Peizhen，Ran Yongkang，et al.China activity structure basic characteristics［J］.Science in China(Series D)，2002，32(12):1 020－1 030)