孟令媛，史海霞，周龍泉

0 引言

喜馬拉雅造山帶包含喜馬拉雅弧、東構(gòu)造結(jié)和西構(gòu)造結(jié)3個基本組成部分，是大陸碰撞后印度板塊繼續(xù)向北移動，并向西藏高原下俯沖產(chǎn)生的構(gòu)造變形系統(tǒng)向南凸出的弧形，東、西構(gòu)造結(jié)的變形過程在時間、空間上的差異，是西藏東部下地殼向東、東南流動的部分原因 (滕吉文等，2006;馬曉靜，高祥林，2011)。依據(jù)構(gòu)造地質(zhì)學(xué)劃分喜馬拉雅造山帶區(qū)域，東構(gòu)造結(jié)為北東走向的馬蹄形同軸弧，主要是由主邊界逆斷層?xùn)|端和緬甸弧構(gòu)成，還包含幾條大致平行的斷層弧，最顯著的構(gòu)造特征是環(huán)繞東構(gòu)造結(jié)的大尺度順時針旋轉(zhuǎn)運(yùn)動，以及緬甸弧東邊的寬闊右旋剪切變形帶，其中實(shí)皆 (Sagaing)斷裂是緬甸弧及東構(gòu)造結(jié)的重要組成部分，也是整個喜馬拉雅造山帶的東邊界 (尹安，2006)。實(shí)皆斷裂是緬甸中部近南北向高角度右旋走滑斷裂，長約1 200 km，GPS資料研究表明，該斷裂右旋走滑速率為15～18 mm/a(宋健等，2011)，是東構(gòu)造結(jié)強(qiáng)震最為活躍的斷裂帶，且靠近我國滇西南地區(qū)，自1900年以來發(fā)生過多次M≥7地震，最近一次為2012年11月11日的緬甸MS7.0地震，其震中位置為發(fā)生在實(shí)皆斷裂帶上的1946年MS7.8地震和1956年MS7.0地震中間的空段內(nèi)，震源深度約20 km，震中位置距離我國云南省邊境約200 km①http://earthquake.usgs.gov/earthauakes/eqarchives/epic/result.php．。

近年來，國際地震界廣泛關(guān)注強(qiáng)震產(chǎn)生的庫侖破裂應(yīng)力變化及其對周圍斷層的影響，主要涉及后續(xù)地震的發(fā)震位置和時間 (Harris，1998;Stein，1999)。即某一斷層發(fā)生地震時，不僅在該發(fā)震斷層上產(chǎn)生應(yīng)力降及應(yīng)變能的釋放，同時應(yīng)力還會被傳遞至其周圍的斷層上，進(jìn)而使得這些斷層發(fā)生應(yīng)力的改變和調(diào)整 (Harris et al.，1995;Toda et al.，1998)。萬永革等 (2002)針對中國大陸發(fā)生的4次大震震后應(yīng)力變化情況指出，后續(xù)破裂事件均受到前面發(fā)生的地震事件的影響，即前面的地震事件對后續(xù)地震事件均存在某種程度的觸發(fā)作用。斷層之間的相互作用不僅可以通過觸發(fā)來體現(xiàn)，而且當(dāng)某一斷層發(fā)生地震時，對其周圍斷層應(yīng)力的影響或者有增加或者有減少，即加速或者延緩地震的發(fā)生時間，產(chǎn)生觸震和緩震效應(yīng) (King et al.，1994;張秋文等，2003)，如1857年Ft.Tejon大地震產(chǎn)生的緩震效應(yīng)對南加州后續(xù)地震的抑制作用十分明顯 (Harris，Simpson，1996)。此外，針對2008年汶川地震后對周邊斷層地震活動性影響的研究中指出，汶川地震的發(fā)生使得鮮水河斷裂道孚—康定段及則木河斷裂的庫侖破裂應(yīng)力增加，金沙江斷裂的庫侖破裂應(yīng)力減小 (萬永革等，2009;邵志剛等，2010)。

本文針對2000年以來喜馬拉雅東構(gòu)造結(jié)發(fā)生的3次7級左右地震后周圍應(yīng)力場的變化進(jìn)行分析，計(jì)算其對我國滇西南地區(qū)斷層靜態(tài)庫侖破裂應(yīng)力的具體影響，討論滇西南斷層間的相互影響及加卸載效應(yīng)，并估計(jì)該區(qū)域潛在的地震危險性。

1 計(jì)算方法

地震的發(fā)生主要由地下巖石的錯動造成，地震發(fā)生前后較短時間內(nèi)實(shí)現(xiàn)了彈性能的快速釋放，假設(shè)地球介質(zhì)為半無限空間且空間內(nèi)各向同性，若已知地震位錯的幾何參數(shù)和錯動大小，即可求出半無限空間內(nèi)部的位移場和應(yīng)變場?；趲靵龆桑瑪鄬用嫔蠋靵銎屏褢?yīng)力變化的表達(dá)式為

式中:Δτ為剪切應(yīng)力變化;Δσn為正應(yīng)力變化;f為斷層面介質(zhì)摩擦系數(shù);Δp為孔隙壓力變化，實(shí)際上Δp調(diào)節(jié)著斷層面上的有效正應(yīng)力，Δp減少摩擦系數(shù)的效應(yīng)可以用f'=f(1－B)表示，其中B為Skemptons系數(shù)，取值范圍為0～1，當(dāng)巖石應(yīng)力的改變遠(yuǎn)遠(yuǎn)快于巖石中Δp的擴(kuò)散，則Δp可以通過B值在式 (1)中得到體現(xiàn)，且轉(zhuǎn)換為

其中，f'包含了孔隙流體和斷層面上介質(zhì)的特性，一般取值為0.2～0.8(King et al.，1994)。

基于式 (1)和 (2)，采用Toda等 (1998)開發(fā)的Coulomb3.3軟件對庫侖破裂應(yīng)力變化進(jìn)行計(jì)算。計(jì)算中地殼的剪切模量取3.3×104MPa，泊松比取值 0.25，f'取值 0.4，該值為 King等(1994)通過不同f'取值相應(yīng)庫侖破裂應(yīng)力圖像變化情況給出的經(jīng)驗(yàn)取值。通常引起應(yīng)力變化的地震和破裂面稱為“源地震”和“源斷層”，區(qū)域中的其它先存在的斷層為“接收斷層”，Coulomb3.3軟件的計(jì)算過程中，接收斷層的ΔCFS的計(jì)算結(jié)果受到源斷層的走向、傾角和滑動角，源地震的同震位錯，以及接收斷層的幾何特征和摩擦系數(shù)的綜合影響 (King et al.，1994)。

活動斷層庫侖破裂應(yīng)力的變化對其地震的活動性存在著一定的影響，簡單的對應(yīng)關(guān)系即為應(yīng)力增加，則地震活動性增強(qiáng)，最后會慢慢回歸到背景的水平，而應(yīng)力減小則可能使周圍的地震活動性的進(jìn)程受到延緩。斷層的演化過程一般可分為3個階段:穩(wěn)定滑動階段、加速階段、瞬態(tài)失穩(wěn)階段。第一階段即為斷層的閉鎖階段，通常占到整個演化階段的約90%的時間進(jìn)程，而同震過程則對應(yīng)于第三階段，即斷層的動態(tài)破裂，通常在幾十秒到上百秒內(nèi)完成。彈性位錯 (Elastic Dislocation)理論下的觸震效應(yīng)和緩震效應(yīng)的時間變化量Δt可表示為

式中，τ'為剪應(yīng)力加載速率，式 (3)被廣泛應(yīng)用于對地震加卸載區(qū)成因的討論和解釋 (Harris，1998;Toda et al.，1998;Parsons et al.，2000)。剪應(yīng)力加載速率τ'無法直接測得，需要基于理論計(jì)算求得，通常根據(jù)剪應(yīng)力加載速率τ'和靜態(tài)應(yīng)力降Δτs及地震回復(fù)周期tr之間的近似關(guān)系 (Dieterich，1994)求得:

式中，Δτs為主震靜態(tài)應(yīng)力降，可以表示為(Knopoff，1958)

其中:μ為剪切模量，且μ=3×1010Pa;W為發(fā)震斷層寬度;D為地震的平均滑動位移。

在計(jì)算過程中，選取喜馬拉雅東構(gòu)造結(jié)靠近滇西南地區(qū)2000年以來發(fā)生的2000年6月7日緬甸MS6.9地震、2011年3月24日緬甸MS7.2地震和2012年11月11日緬甸MS7.0地震作為源地震(圖1)。把滇西南地區(qū)靠近邊境的大盈江斷裂、龍陵—瑞麗斷裂、南汀河西支斷裂、南汀河?xùn)|支斷裂、孟連—瀾滄斷裂、打洛—景洪斷裂、怒江斷裂、木嘎斷裂和龍陵—瀾滄斷裂9條斷層設(shè)定為接收斷層，并將9條斷層各自的走向、傾角和滑動角作為相應(yīng)接收斷層的幾何輸入?yún)?shù) (圖1)。根據(jù)以上方法和程序，計(jì)算這3次地震引起的周圍應(yīng)力的變化情況，并針對滇西南地區(qū)大盈江斷裂等9個斷層所受到的ΔCFS進(jìn)行計(jì)算，從而進(jìn)一步分析喜馬拉雅東構(gòu)造結(jié)附近強(qiáng)震對滇西南地區(qū)斷層的加卸載影響。

2 地震參數(shù)選取及周圍斷層參數(shù)的設(shè)定

2000年6月7日緬甸MS6.9地震、2011年3月24日緬甸MS7.2地震和2012年11月11日緬甸MS7.0地震3次地震的發(fā)震斷層產(chǎn)狀以及位錯參數(shù)見表1。為計(jì)算滇西南地區(qū)庫侖破裂應(yīng)力的變化，需要明確計(jì)算區(qū)域內(nèi)斷層的幾何參數(shù)及運(yùn)動學(xué)特征，筆者主要參照鄧起東等 (2002)依據(jù)地質(zhì)和地震信息提供的活動斷層相關(guān)數(shù)據(jù)來構(gòu)建斷層模型參數(shù) (表2)，這些斷層即為應(yīng)用Coulomb3.3軟件的計(jì)算過程的接收斷層。圖1中給出的F1～F9斷層主要考慮滇西南地區(qū)靠近緬甸3次7級左右強(qiáng)震的位置及影響范圍，實(shí)皆斷裂為近NS向的走滑型，因此，F(xiàn)1～F6均為NE向，其中F6打洛—景洪斷裂展布相對較長，且靠近2011年緬甸MS7.2地震，F(xiàn)7為距離2000年緬甸MS6.9、2012年緬甸MS7.0兩次地震較近，且與實(shí)皆斷裂近乎平行的怒江斷裂，F(xiàn)8、F9則為靠近2011年地震的兩條NW向斷裂，其中F9為龍陵—瀾滄斷裂的南段部分，本文考慮的計(jì)算長度約占總長的1/4。

表1 喜馬拉雅東構(gòu)造結(jié)附近3次7.0級左右地震的發(fā)震斷層幾何參數(shù)及位錯參數(shù)Tab.1 Geometric and dislocation parameters on the seismogenic fault of three M7.0 earthquakes around the eastern Himalayan tectonic syntax is

表2 滇西南地區(qū)接收斷層幾何參數(shù)Tab.2 Geometric parameters of the receiver faults in Southwest Yunnan area

3 計(jì)算結(jié)果

根據(jù)已設(shè)定的地震斷層模型和活動斷層性質(zhì)及滑動特性，計(jì)算投影到斷層面和滑動方向上的庫侖破裂應(yīng)力變化 (圖2)。本研究所選取的3次7級左右地震對滇西南地區(qū)活動斷層的影響主要以應(yīng)力卸載為主，其中大盈江斷裂 (F1)、龍陵—瑞麗斷裂 (F2)、南汀河西支斷裂 (F3)、南汀河?xùn)|支斷裂 (F4)、孟連—瀾滄斷裂 (F5)和打洛—景洪斷裂 (F6)呈NE向展布，與滇西南地區(qū)西側(cè)靠近實(shí)皆斷裂的整體構(gòu)造相對一致，庫侖破裂應(yīng)力均為負(fù)值，受應(yīng)力卸載的影響，量值范圍在－0.06～－0.01 bar之間，這6條NE向斷裂的ΔCFS由北 至 南 依 次 為 － 0.01 bar、 － 0.04 bar、－0.02 bar、 －0.02 bar、 －0.03 bar和 0.06 bar;怒江斷裂 (F7)、木嘎斷裂 (F8)和龍陵—瀾滄斷裂 (F9)為NW向展布，其中怒江斷裂 (F7)為應(yīng)力卸載影響，ΔCFS為 －0.02 bar，木嘎斷裂(F8)和龍陵—瀾滄斷裂 (F9)為應(yīng)力加載影響，ΔCFS分別為0.002 bar和0.004 bar。

參照圖1，2000年緬甸MS6.9地震發(fā)生之后至2011年緬甸MS7.2地震之間，滇西南地區(qū)M≥5.0地震主要集中分布在大盈江斷裂 (F1)和龍陵—瑞麗斷裂 (F2)附近，相對平靜的斷層周圍為南汀河西支斷裂 (F3)、南汀河?xùn)|支斷裂 (F4)和怒江斷裂 (F7);緬甸2011年MS7.2地震發(fā)生之后截至2012年12月31日，滇西南地區(qū)M≥5.0地震仍分布在F1和F2附近區(qū)域。對比圖2，可將滇西南靠近緬甸一側(cè)的9條活動斷層大體分為3個區(qū)域:第一個區(qū)域 (A)包含 NE向的 F1、F2和NNE向的F7;第二個區(qū)域 (B)包含2個NE向的F3和F4;第三個區(qū)域 (C)則由NE向的F5、F6和NW向的F8、F9組成。A區(qū)和B區(qū)均為庫侖破裂應(yīng)力卸載區(qū)域，A區(qū)持續(xù)有M≥5.0地震活動，但沒有達(dá)到M≥6.0，B區(qū)自2000年緬甸MS6.9地震之后持續(xù)M≥5.0地震平靜;C區(qū)為庫侖破裂應(yīng)力加載區(qū)域，且自2011年緬甸MS7.2地震之后至2012年底，C區(qū)的地震活動同樣處于持續(xù)的M≥5.0地震平靜之中。

針對本研究涉及的9條活動斷層，定義其分別作為特征斷層歷史距今發(fā)生M≥5.0地震的時間間隔為視復(fù)發(fā)周期，即tr*，最近一次M≥5.0地震的主震的靜態(tài)應(yīng)力降為Δτrs，則基于式 (4)則有τ'=Δτrs/tr*。利用圖2和式 (3)，可以近似給出本文所涉及的9條活動斷層的觸震和緩震效應(yīng)推算結(jié)果 (表3)及示意圖 (圖3)。本研究給出的結(jié)果是假設(shè)庫侖破裂應(yīng)力的擾動發(fā)生在斷層演化的第二階段，即加速階段，則表3中利用式 (3)給出的Δt計(jì)算結(jié)果均為假設(shè)所分析的活動斷層已處在加速階段的前提下，受到ΔCFS擾動的觸震和緩震影響。

表3 滇西南地區(qū)接收斷層觸震/緩震計(jì)算結(jié)果Tab.3 Results of the trigging and delaying on the receiver faults in Southwest Yunnan area

實(shí)際上，關(guān)于視復(fù)發(fā)周期的定義，一方面，特征斷層上歷史上距今最近一次的大地震和前一次大地震的間隔可以定義為視復(fù)發(fā)周期;另一方面，如果距今僅發(fā)生過一次較大地震，由于本研究給出的結(jié)果是假設(shè)庫侖破裂應(yīng)力的擾動發(fā)生在斷層演化的第二階段，則可以認(rèn)為這一次地震距離現(xiàn)在的時間間隔為視復(fù)發(fā)周期。

圖3中C區(qū)則為觸震效應(yīng)的影響區(qū)，對比“M7專項(xiàng)”針對南北帶南段—中—長期大地震危險性的研究成果，指出騰沖—瑞麗地區(qū)為M≥7.0地震中—長期危險性值得注意地區(qū) (M7專項(xiàng)工作組，2012;聞學(xué)澤等，2012)。對比圖3中A區(qū)和騰沖—瑞麗注意地區(qū)不包含怒江斷裂 (F7)，主要以大盈江斷裂 (F1)和龍陵—瑞麗斷裂 (F3)為主，B區(qū)和C區(qū)孟連—瀾滄斷裂 (F5)的附近區(qū)域?yàn)榈卣鹂諈^(qū)部分，C區(qū)龍陵—瀾滄斷裂附近區(qū)域?yàn)閺?qiáng)震破裂區(qū)。

4 討論與結(jié)論

本研究主要涉及喜馬拉雅東構(gòu)造結(jié)2000年以來的3次7級左右地震對我國滇西南地區(qū)空間應(yīng)力變化及活動斷層靜態(tài)庫侖應(yīng)力變化的影響，這3次7級左右地震震中位置均在我國西南邊境境外區(qū)域。同時，由于東構(gòu)造結(jié)范圍較廣，且實(shí)皆斷裂展布較長，本研究計(jì)算區(qū)域大于1×106km2，受到計(jì)算能力的限制，在針對滇西南地區(qū)活動斷層計(jì)算庫侖破裂應(yīng)力變化時，沒有對單一活動斷層進(jìn)行分段化處理，而是將每一個活動斷層作為整體計(jì)算ΔCFS取值。3次7級左右地震的震源深度范圍大致在20 km(表1)，可以認(rèn)為該區(qū)域地殼的脆性破裂向下不超過25 km。在庫侖破裂應(yīng)力變化計(jì)算過程中，將斷層面深度范圍取值地表至地下20 km處，根據(jù)Coulomb3.3軟件的算法原理，現(xiàn)有的計(jì)算結(jié)果為地表以下10 km深度庫侖破裂應(yīng)力的變化情況。結(jié)果顯示，選取的3次7級左右地震對滇西南地區(qū)活動斷層的影響主要以應(yīng)力卸載為主，其中大盈江斷裂、龍陵—瑞麗斷裂、南汀河西支斷裂、南汀河?xùn)|支斷裂、孟連—瀾滄斷裂和打洛—景洪斷裂呈現(xiàn)NE向展布，與滇西南地區(qū)西側(cè)靠近實(shí)皆斷裂的整體構(gòu)造相對一致，庫侖破裂應(yīng)力均為負(fù)值，受到應(yīng)力卸載的影響，量值范圍在－0.06～－0.01 bar;NW向的怒江斷裂為應(yīng)力卸載影響，ΔCFS為－0.02 bar，而木嘎斷裂和龍陵—瀾滄斷裂為應(yīng)力加載影響，ΔCFS分別為+0.002 bar和 +0.004 bar。

實(shí)際上，在本文所選取的9條接收斷層境內(nèi)的附近區(qū)域2000年以來也有6級左右的地震發(fā)生，這些地震同樣會對這些接收斷層產(chǎn)生一定的庫倫破裂應(yīng)力影響，本文僅選取了境外的3次7級左右地震作為源地震進(jìn)行研究，主要是考慮了這3個震例分布由北至南，基本展布滇西南地區(qū)的境外空間。基于現(xiàn)有的計(jì)算結(jié)果，得出初步結(jié)論為:2000年的緬甸MS6.9地震和2012年的緬甸MS7.0地震距離滇西南地區(qū)接收斷層距離較遠(yuǎn)，發(fā)震斷層空間展布以近南北向?yàn)橹?(圖1)，二者的影響較小且主要以應(yīng)力卸載為主;而2011年的緬甸MS7.2地震，距滇西南地區(qū)，近且發(fā)震斷層展布為NE向，其對周圍接收斷層的加卸載影響則均相對略大。這一初步結(jié)論，對這一區(qū)域未來發(fā)生的7級左右強(qiáng)震對滇西南地區(qū)斷層的庫侖破裂應(yīng)力影響有一定的參照意義。

通常情況下觸震和緩震的作用是相對的，并且M≥5.0地震的成核過程也要相對復(fù)雜，孕育時間也不是絕對的。3次7級左右地震對大盈江斷裂、龍陵—瑞麗斷裂、南汀河西支斷裂、南汀河?xùn)|支斷裂、孟連—瀾滄斷裂、打洛—景洪斷裂、怒江斷裂的緩震效應(yīng)影響在－0.69～－0.21 a，對龍陵—瀾滄斷裂的觸震影響為+0.35 a(表3)。表3中給出的9條活動斷層的ΔCFS的絕對數(shù)值也較小，主要是由于所選擇的這3個7級左右地震距離9條活動斷層相對較遠(yuǎn)，但庫侖破裂應(yīng)力的影響是存在的，由式 (3)推算出來的觸震和緩震效應(yīng)的Δt也是在假設(shè)擾動發(fā)生在斷層演化的加速階段中。因此，表3中給出的Δt的計(jì)算結(jié)果對于設(shè)定區(qū)域短臨期 (1 a左右尺度)地震活動性的研究具備一定的參考價值。

本文在撰寫過程中，劉桂萍研究員、劉杰研究員、蔣海昆研究員和張永仙研究員給予了有益的指導(dǎo)與幫助，作者在此一并表示感謝。

鄧起東，于貴華，葉文華.1992.地震地表破裂參數(shù)與震級關(guān)系的研究［A］//國家地震局地質(zhì)研究所.活動斷裂研究(2).北京:地震出版社，247－264.

鄧起東，張培震，冉永康，等.2002.中國活動構(gòu)造基本特征［J］．中國科學(xué)(D輯)，32(12):1020－1030.

馬曉靜，高祥林.2011.橫跨喜馬拉雅造山帶的構(gòu)造運(yùn)動轉(zhuǎn)換與變形分配［J］．地球物理學(xué)報，54(6):1528－1535.

邵志剛，周龍泉，蔣長勝，等.2010.2008年汶川MS8.0地震對周邊斷層地震活動的影響［J］．地球物理學(xué)報，53(8):1784－1795.

宋健，唐方頭，鄧志輝，等.2011.喜馬拉雅東構(gòu)造結(jié)周邊地區(qū)主要斷裂現(xiàn)今運(yùn)動特征與數(shù)值模擬研究［J］．地球物理學(xué)報，54(6):1536－1548.

滕吉文，王謙身，王光杰，等.2006.喜馬拉雅“東構(gòu)造結(jié)”地區(qū)的特異重力場與深部構(gòu)造［J］．地球物理學(xué)報，49(4):1045－1052.

萬永革，沈正康，盛書中，等.2009.2008年汶川大地震對周圍斷層的影響［J］．地震學(xué)報，31(2):128－139.

萬永革，吳忠良，周公威，等.2002.幾次復(fù)雜地震中不同破裂事件之間的“應(yīng)力觸發(fā)”問題［J］．地震學(xué)報，22(6):568－576.

尹安.2006.喜馬拉雅造山帶新生代構(gòu)造演化:沿走向變化的構(gòu)造幾何形變、剝露歷史和前陸沉積的約束［J］．地學(xué)前緣，13(5):416－515.

張秋文，張培震，王乘，等.2003.鮮水河斷裂帶斷層間相互作用的觸震與緩震效應(yīng)［J］．地震學(xué)報，25(2):143 －153.

Dieterich J.H..1994.A Constitutive Law for Rate of Earthquake Production and its Application to Earthquake Clustering［J］．J.Geophys.Res.，99(B2):2601 －2618.

Gkarlaouni C.，Papadimitriou E.E.，Karakostas V.G.，et al..2008. 斷層相互作用對中國川滇地區(qū)地震危險評估的意義［J］．地震學(xué)報，30(2):176－197.

Harris R.A.，Simpson R.W.，Reasenberg P.A..1995.Influence of Static Stress Changes on Earthquake Locations in Southern California［J］．Nature，375(6528):221 －224.

Harris R.A.，Simpson R.W..1996.In the Shadow of 1957—the Effect of the Great Ft Tejon Earthquake on Subsequent Earthquake in Southern California［J］．Geophys Res Lett，23(3):229 － 232.

Harris R.A..1998.Introduction to Special Section:Stress Triggers，Stress Shadows，and Implications for Seismic Hazard［J］．J Geophys Res，103(B10):24347－24358.

King G.C.P.，Stein R.S.，Lin J..1994.Static Stress Changes and the Triggering of Earthquakes［J］．Bull.Seismol.Soc.Amer.，84(3):935－953.

Knopoff L..1958.Energy Release in Earthquakes［J］．Geophys.J.，1(1):44－52.

M7專項(xiàng)工作組.2012.中國大陸大地震中—長期危險性研究［M］．北京:地震出版社，234－248.

Parsons T.S.，Toda R.S.，Stein A.，et al..2000.Heightened Odds of Large Earthquakes Near Istanbul:an Interaction-based Probability Calculation［J］．Science，288(5466)，661 －665.

Stein R.S..1999.The Role of Stress Transfer in Earthquake Occurrence 1999［J］．Nature，402(6762):605 － 609.

Toda S.，Stein R.S.，Reasenberg P.A.，et al..1998.Stress Transferred by the 1995 M(w)=6.9 Kobe，Japan，Shock-Effect on Aftershocks and Future Earthquake Pobabilities［J］．Journal of geophysical research，103(B10)，24543 －524565.