虎雄林，解朝娣，劉麗芳

(1.云南省地震局，云南昆明650224;2.云南大學(xué)地球物理系，云南昆明650091)

0 引言

2003年7月21日、10月16日在云南省大姚縣先后發(fā)生了6.2、6.1級(jí)地震，震中位于滇中構(gòu)造塊體內(nèi)部。該區(qū)歷史上地震活動(dòng)較弱，1993年前，僅記錄到M≥5.0地震5次。然而自1993年以來，該區(qū)地震活動(dòng)明顯增強(qiáng)，中強(qiáng)地震頻發(fā)，相繼發(fā)生了1993年2月1日大姚5.3級(jí)地震、1993年8月14日姚安5.6級(jí)地震、2000年1月15日姚安5.9級(jí)、6.5級(jí)地震和2003年大姚6.2級(jí)、6.1級(jí)地震。云南地區(qū)6級(jí)以上地震主要沿著一些主要斷裂帶分布，地震的發(fā)生都能找到對(duì)應(yīng)的發(fā)震構(gòu)造，但上述這6次地震沒有發(fā)生在明顯的斷裂構(gòu)造附近(蘇有錦，2004)。為了探討該地區(qū)強(qiáng)震發(fā)生的動(dòng)力背景及雙震發(fā)生的機(jī)制，張建國等(2009)，華衛(wèi)等(2006)，周龍泉等(2007)從地震構(gòu)造、地殼速度結(jié)構(gòu)及震源特征等方面對(duì)這6次地震進(jìn)行了研究。

謝富仁等(1994，2011)研究表明，地殼中的應(yīng)力狀態(tài)與地震活動(dòng)關(guān)系密切，活動(dòng)斷裂帶的構(gòu)造應(yīng)力方向、類型變化的部位是易發(fā)生強(qiáng)震的地點(diǎn);均勻區(qū)域應(yīng)力場背景上的局部應(yīng)力變化是強(qiáng)震活動(dòng)相對(duì)集中的地區(qū)。構(gòu)造應(yīng)力場是強(qiáng)震孕育、發(fā)生的動(dòng)力學(xué)來源，而斷層之間相互作用引起的應(yīng)力變化也能觸發(fā)強(qiáng)震。本文旨在探討2003年大姚兩次強(qiáng)震的構(gòu)造應(yīng)力處于什么狀態(tài)?是否與區(qū)域構(gòu)造應(yīng)力一致?兩次地震的構(gòu)造應(yīng)力是否發(fā)生了改變?2003年7月21日大姚6.2級(jí)地震引起的應(yīng)力變化是否觸發(fā)了10月16日大姚6.1級(jí)地震。

構(gòu)造應(yīng)力的基礎(chǔ)資料包括震源機(jī)制解、斷層滑動(dòng)反演、鉆孔崩落、水壓致裂和應(yīng)力解除資料等。其中斷層滑動(dòng)反演、鉆孔崩落、水壓致裂和應(yīng)力解除資料反映的是地殼淺部的應(yīng)力狀態(tài)，而震源機(jī)制解資料主要反映了中、上地殼的應(yīng)力狀態(tài)，地震的發(fā)生一般在十至數(shù)十公里深的地殼內(nèi)部，震源機(jī)制解資料能更好的反映構(gòu)造應(yīng)力環(huán)境。通常對(duì)震源機(jī)制解進(jìn)行統(tǒng)計(jì)分析得到P、T和N軸的優(yōu)勢方位，這個(gè)優(yōu)勢方位代表該地區(qū)的構(gòu)造應(yīng)力場方向。但P、B、T軸并不等同于構(gòu)造應(yīng)力場的3個(gè)主應(yīng)力σ1、σ2和σ3軸，統(tǒng)計(jì)分析得到P、T和N軸的優(yōu)勢方位只能對(duì)構(gòu)造應(yīng)力場進(jìn)行大致的方向性描述，很難準(zhǔn)確描述某一地區(qū)的構(gòu)造應(yīng)力場特征。

筆者利用2003年7月21日、10月16日在大姚6.2、6.1級(jí)地震余震序列的震源機(jī)制解資料，采用滑動(dòng)方向擬合法(許忠淮，1985)反演出兩次強(qiáng)震的構(gòu)造應(yīng)力，對(duì)其進(jìn)行對(duì)比分析，定量給出兩次地震的構(gòu)造應(yīng)力變化量。

地震是應(yīng)力積累到一定程度，地殼巖層突然破裂、錯(cuò)動(dòng)的一種自然現(xiàn)象(徐紀(jì)人等，2008)。強(qiáng)震的同震位錯(cuò)會(huì)引起附近區(qū)域的應(yīng)力變化，其通過持久性地改變施加在斷層上的應(yīng)力來影響后續(xù)地震，應(yīng)力變化將影響斷裂的破裂失穩(wěn)條件，當(dāng)應(yīng)力變化為正值時(shí)，加速了周圍斷裂或臨近區(qū)域的應(yīng)力累積，使下次地震提前發(fā)生，這種作用即為增震作用，也稱為地震觸發(fā)效應(yīng);相反，應(yīng)力變化為負(fù)值時(shí)，周圍斷裂或臨近區(qū)域的應(yīng)力累積進(jìn)程會(huì)被延緩滯后(萬永革等，2000，2002;劉桂萍，傅征祥，2001;張國民，李麗，1997;Harris，1998;吳小平等，2007)。華衛(wèi)等(2006)對(duì)2003年大姚6.2、6.1級(jí)地震庫侖破裂靜應(yīng)力變化研究表明，7月21日6.2級(jí)地震對(duì)10月16日6.1級(jí)地震有明顯的觸發(fā)作用，但其在計(jì)算地震庫侖破裂靜應(yīng)力變化時(shí)沒有考慮構(gòu)造應(yīng)力的影響。

筆者以反演得到的2003年大姚6.2、6.1級(jí)地震構(gòu)造應(yīng)力為基礎(chǔ)，分別計(jì)算兩次強(qiáng)震后產(chǎn)生的靜態(tài)庫侖破裂應(yīng)力變化，探討2003年7月21日大姚6.2級(jí)地震引起的應(yīng)力變化是否觸發(fā)了10月16日大姚6.1級(jí)地震。

1 震源區(qū)構(gòu)造應(yīng)力場反演

1.1 計(jì)算方法

(1)斷層面上剪應(yīng)力的方向

設(shè)控制斷層滑動(dòng)的應(yīng)力張量為T，則作用于斷層面上的力(圖1a)

式中，n為斷層面單位法向矢量。f'可以分解成沿?cái)鄬用娴募魬?yīng)力和垂直斷層面的正應(yīng)力

則斷層面上的剪應(yīng)力

應(yīng)力張量T可以用含有歸一化偏應(yīng)力張量D的形式來表示(Angelier，1979;謝富仁等，1994;許忠淮，戈澍謨，1984)

式中，I為單位矩陣;t1(＞0)，t2為任意數(shù)。由于只有偏應(yīng)力張量D才產(chǎn)生剪應(yīng)力

所以τ與τ'的方向相同，如圖1b所示。

圖1 斷層面上應(yīng)力示意圖(n為斷層面單位法向矢量，f'為斷層面上應(yīng)力，τ'為斷層面上剪應(yīng)力，σn'為斷層面上正應(yīng)力，S為斷層面滑動(dòng)方向，α為斷層面上剪應(yīng)力與滑動(dòng)方向間的夾角;f，σn和τ分別是偏應(yīng)力張量作用于斷層面上的應(yīng)力、正應(yīng)力和剪應(yīng)力)Fig.1 Schematic diagram of stress on the fault plane(n is unit normal vector，f'is stress，τ'is shear stress，σn'is the normal stress，S is the sliding direction，α is angle between shear stress and slip direction on the fault plane;f，σnand τ is the stress，normal stress and shear stress which are deviator stress tensor acting on the fault plane)

滑動(dòng)方向擬合法提出3個(gè)基本假定:(1)斷層沿?cái)嗝嫔霞魬?yīng)力的方向滑動(dòng);(2)一組斷層所分布的區(qū)域內(nèi)存在著平均應(yīng)力場;(3)斷層滑動(dòng)相互獨(dú)立。在此基本假定條件下，滑動(dòng)方向擬合法就是尋找一個(gè)適當(dāng)?shù)钠珣?yīng)力張量D，使得作用于一組斷層面上的剪應(yīng)力τi(i=1，2，…，m;m為斷層數(shù))與斷層滑動(dòng)單位矢量Si之間的夾角α最小(謝富仁等，1994;許忠淮，戈澍謨，1984)(圖1a)。由于t1、t2的任意性，只能由一組斷層滑動(dòng)單位矢量確定應(yīng)力張量T的6個(gè)獨(dú)立分量中的4個(gè)，即3個(gè)主應(yīng)力的方向和1個(gè)表示主應(yīng)力相對(duì)大小的應(yīng)力形因子

(2)斷層面上剪應(yīng)力的相對(duì)大小

若要使斷層真正產(chǎn)生滑動(dòng)，剪應(yīng)力的大小應(yīng)符合摩擦滑動(dòng)準(zhǔn)則

其中，μ是摩擦系數(shù)，C0為聚合強(qiáng)度。對(duì)于斷層C0可以忽略不計(jì)(Brace，Kohlstedt，1980;Zoback，1992)，式(7)可改寫成

由于只有偏應(yīng)力張量D才產(chǎn)生沿?cái)鄬用娴募魬?yīng)力，球應(yīng)力張量部分只對(duì)正應(yīng)力有影響，由式(3)、(4)和(5)推出(圖1b)

也就是說，當(dāng)偏應(yīng)力張量作用在斷層面上剪應(yīng)力與正應(yīng)力的大小之比滿足

時(shí)，斷層才能產(chǎn)生滑動(dòng)。

1.2 基礎(chǔ)資料

據(jù)云南數(shù)字地震臺(tái)網(wǎng)測定，2003年7月21日23時(shí)16分30秒，云南大姚發(fā)生了M6.2地震，震中位置(25.95°N，101.23°E)，震源深度為6 km。之后，10月16日20時(shí)28分03秒，在該區(qū)再次發(fā)生M6.1地震，震中位置(25.92°N，101.30°E)，震源深度為5 km，位于第1次地震東南向約8 km處(圖2)。震后，云南省地震局在地震現(xiàn)場布設(shè)了臨時(shí)地震臺(tái)網(wǎng)，自震后第二天進(jìn)行了為期一個(gè)月的近源觀測，監(jiān)測余震的發(fā)展趨勢，獲得了豐富的數(shù)字化余震記錄。讀取臺(tái)網(wǎng)接收到的波形資料中的P波初動(dòng)和振幅比，反演得到了兩次主震和49個(gè)余震的震源機(jī)制解(表1、表2)。其中大姚6.2級(jí)地震有37個(gè)余震震源機(jī)制解，大姚6.1級(jí)地震有12個(gè)余震震源機(jī)制解。

圖2 大姚地區(qū)地震構(gòu)造圖Fig.2 Seismotectonic map of Dayao area

表1 大姚M6.2、M6.1地震震源機(jī)制解Tab.1 Focal mechanism solution of Dayao M6.2，M6.1 earthquakes

1.3 區(qū)域構(gòu)造應(yīng)力場

由大姚6.2級(jí)地震的37個(gè)余震、6.1級(jí)地震的12個(gè)余震的震源機(jī)制解資料，利用Max Wyss和Stefan Wiemer的震源機(jī)制應(yīng)力反演軟件ZMAP程序，分別反演了大姚6.2、6.1級(jí)地震時(shí)，在大姚附近地區(qū)的區(qū)域構(gòu)造應(yīng)力場(表3)。從圖3a可以看出，大姚6.2級(jí)地震所在地區(qū)的區(qū)域最大主壓應(yīng)力以水平分量為主(傾角為28°)，方位角呈南西向(方位角為 240°)，代表應(yīng)力值相對(duì)大小的R值為0.4。從圖3b可以看出，大姚6.1級(jí)地震所在地區(qū)的區(qū)域最大主壓應(yīng)力以水平分量為主(傾角為6°)，方位角呈南西向(方位角為250°)，代表應(yīng)力值相對(duì)大小的R值為0.6。

?

?

表3 大姚6.2、6.1級(jí)地震構(gòu)造應(yīng)力場反演結(jié)果Tab.3 Inversion results of tectonic stress field of Dayao MS6.2，MS6.1 earthquakes

圖3 反演得到的大姚6.2級(jí)(a)、6.1級(jí)(b)地震區(qū)域構(gòu)造應(yīng)力場Fig.3 Regional tectonic stress field of Dayao MS6.2(a)，MS6.1(b)earthquakes obtained in inversion

崔效鋒和謝富仁(1999)通過利用震源機(jī)制解反演我國西南及鄰區(qū)構(gòu)造應(yīng)力，對(duì)應(yīng)力場進(jìn)行了分區(qū)，大姚兩次強(qiáng)震所在的區(qū)域位于應(yīng)力場分區(qū)的D區(qū)(表4)。對(duì)比反演得到的大姚兩次地震構(gòu)造應(yīng)力與應(yīng)力場分區(qū)D區(qū)構(gòu)造應(yīng)力，可以發(fā)現(xiàn)兩次地震構(gòu)造應(yīng)力最大主應(yīng)力軸S1方位分別相差21°和11°，仰角分別相差5°和17°;R值分別相差0.01和19。雖然最大主應(yīng)力軸S1方位分別相差21°和11°，但其方位均為北北東向，反演得到大姚兩次地震構(gòu)造應(yīng)力的最大主應(yīng)力軸S1方位與區(qū)域構(gòu)造應(yīng)力基本一致。

表4 滇西南地區(qū)D區(qū)構(gòu)造應(yīng)力場(崔效鋒，謝富仁，1999)Tab.4 Tectonic stress field of D-area in Southwestern Yunnan(Cui，Xie，1999)

3 最優(yōu)取向斷層面投影的庫侖破裂應(yīng)力變化計(jì)算

3.1 原理和方法

根據(jù)庫侖破裂準(zhǔn)則，使巖石破裂失穩(wěn)的庫侖破裂應(yīng)力為(King et al，1994)

式中，σβ為斷層面上的正應(yīng)力，τβ為切應(yīng)力，p為孔隙流體壓力，μ為斷層面介質(zhì)的摩擦系數(shù)。

如果斷層面與主應(yīng)力σ1軸的夾角為β(圖4)，則正應(yīng)力和切應(yīng)力可以表示為

其中，σ1是最大主應(yīng)力;σ3是最小主應(yīng)力。式(11)為

對(duì)式(14)求β的偏微分，可以發(fā)現(xiàn)β滿足cot2β= －時(shí)，庫侖破裂應(yīng)力取最大值。

圖4 最優(yōu)取向斷層面上應(yīng)力的坐標(biāo)系示意圖Fig.4 Sketch map of coordinate system on the fault plane stress in optimal orientation

如式(11)所示，孔隙流體壓力p調(diào)制著斷層面上的有效正應(yīng)力，當(dāng)巖石應(yīng)力的改變遠(yuǎn)遠(yuǎn)快于巖石中的流體壓力擴(kuò)散，則流體壓力p可以通過Skemptons系數(shù)B在式(11)中得到反映，并取視摩擦系數(shù)μ'=μ(1－b)，則式(11)變?yōu)?/p>

μ'包括了孔隙流體和斷層面上的介質(zhì)特性，一般μ'取0.2～0.8。由于難以對(duì)地殼應(yīng)力值直接進(jìn)行精確測量，因此利用庫侖破裂應(yīng)力變化可估計(jì)應(yīng)力的相對(duì)變化，式(15)可改寫為

式中，Δσβ為斷層面上的正應(yīng)力變化，Δτβ為切應(yīng)力變化。

庫侖破裂應(yīng)力變化定義于具體的斷層面，所謂最優(yōu)取向斷層面是指計(jì)算庫侖破裂應(yīng)力變化時(shí)，某一產(chǎn)狀接收斷層面上計(jì)算得到的庫侖應(yīng)力變化大于同一地點(diǎn)其他產(chǎn)狀的接收斷層面上得到的應(yīng)力變化幅值，這類產(chǎn)狀的接收斷層即是最優(yōu)取向斷層。以最優(yōu)取向斷層面投影計(jì)算得到的應(yīng)力變化空間分布，可以解釋余震分布情況，也可以預(yù)測未來后續(xù)地震活動(dòng)的空間分布。

要計(jì)算最優(yōu)取向斷層面上的庫侖破裂應(yīng)力變化，首先用Okada(1992)給出的靜態(tài)位移和應(yīng)變的解析表達(dá)式，計(jì)算震源斷層在接收斷層面處產(chǎn)生的應(yīng)變分量，由胡克定律計(jì)算出應(yīng)力分量。然后將已計(jì)算出的和區(qū)域構(gòu)造應(yīng)力疊加得到總應(yīng)力:

此時(shí)最優(yōu)取向斷層面的方向ψ0通過(θ±β)得以確定。

雖然最優(yōu)取向斷層面的方向性通過總應(yīng)力確定，但是該斷層面上的正應(yīng)力和切應(yīng)力變化值卻只取決于震源斷層產(chǎn)生的應(yīng)力變化，因此在最優(yōu)取向斷層面上的庫侖破裂應(yīng)力變化表示為

求得 σ33和 τ13的相對(duì)變化 Δσ33與 Δτ13，進(jìn)而最優(yōu)取向投影方向上的庫侖破裂應(yīng)力變化(King et al，1994)

3.2 計(jì)算步驟

筆者使用的大姚地區(qū)的地殼模型見表5，該模型參考了國家地震局82測深工程有關(guān)資料和闞榮舉和林中洋(1986)，王椿鏞等(2002)，吳建平等(2001)，胡家富等(2003)，周龍泉等(2007)有關(guān)云南地殼結(jié)構(gòu)的研究，并在上部增設(shè)了一個(gè)彈性半空間，表5中S波速度VS由經(jīng)驗(yàn)關(guān)系式VP=1.73VS確定，密度由周真恒等(2001)對(duì)云南地殼和上地幔的巖石學(xué)結(jié)構(gòu)研究給出。

研究區(qū)域以大姚6.2級(jí)地震震中為中心的(100×100)km2的區(qū)域，如圖5a所示，圖中中間粗線表示地震斷層在地面的投影。

利用P波初動(dòng)和振幅比對(duì)大姚M6.2、M6.1級(jí)地震序列求解得到的震源機(jī)制解(表1)，結(jié)合兩次地震的等震線圖(圖2)，得到兩次地震的斷層參數(shù)，如表6所示，斷層長度據(jù)華衛(wèi)等(2006)運(yùn)用地震標(biāo)定律對(duì)兩次地震研究得到的結(jié)果，分別取15.1 km和13.8 km。

表5 大姚地區(qū)的地殼模型Tab.5 The crustal model in Dayao area

表6 大姚主震斷層參數(shù)Tab.6 Fault parameters of Dayao MS6.2、MS6.1 earthquakes

筆者利用Robinson和McGinty(2000)根據(jù)Okada(1992)方法編制的GNStress程序計(jì)算地震產(chǎn)生的庫侖破裂靜應(yīng)力變化。取μ'=0.75，計(jì)算結(jié)果見圖5，圖中藍(lán)色區(qū)域代表庫侖破裂應(yīng)力變化為負(fù)的地區(qū)，紅色區(qū)域代表庫侖破裂應(yīng)力變化為正的地區(qū)，空心圓圈代表地震。

3.3 庫侖破裂靜應(yīng)力變化

根據(jù)表5中地殼模型和表6中大姚6.2級(jí)地震斷層參數(shù)，以反演得到的大姚6.2級(jí)地震震源區(qū)構(gòu)造應(yīng)力場為基礎(chǔ)建立最佳斷層投影面。由于大姚6.2級(jí)地震震源深度在5 km。因此，我們計(jì)算此次地震在深度為5 km處的地震庫侖破裂靜應(yīng)力變化(圖5a)。由圖可見，6.2級(jí)地震位于庫侖破裂靜應(yīng)力變化為正值的地區(qū)，庫侖破裂靜應(yīng)力變化大約在0.16～0.179 MPa之間，靜態(tài)庫侖破裂應(yīng)力變化觸發(fā)閾值一般為0.1 MPa，這表明6.1級(jí)地震可能是由6.2級(jí)地震觸發(fā)的結(jié)果。

根據(jù)地殼模型(表5)和大姚6.1級(jí)地震斷層參數(shù)(表6)，以反演得到的大姚6.2、6.1級(jí)地震震源區(qū)構(gòu)造應(yīng)力場為基礎(chǔ)建立最佳斷層投影面。由于10月16日6.1級(jí)地震序列震源深度在2～10 km之間，平均深度6 km。因此，我們計(jì)算的是6.2、6.1級(jí)地震在深度為6 km處的地震庫侖破裂靜應(yīng)力變化(圖5b)。由圖可見，庫侖破裂靜應(yīng)力變化為正值的地區(qū)最大值為0.098 1 MPa，小于靜態(tài)庫侖破裂應(yīng)力變化觸發(fā)閾值0.1 MPa。這表明6.1級(jí)地震發(fā)生后，由6.2、6.1級(jí)地震共同產(chǎn)生的靜態(tài)庫侖破裂應(yīng)力為正值的地區(qū)雖然面積較大，但其最大值為0.098 1，小于靜態(tài)庫侖破裂應(yīng)力變化觸發(fā)閾值，其對(duì)后續(xù)地震的發(fā)生作用較弱。

圖5 地震庫侖破裂靜應(yīng)力變化計(jì)算結(jié)果(a)大姚MS6.2地震;(b)大姚MS6.1地震Fig.5 Calculation results of variation of seismic coulomb static stress(a)Dayao MS6.2 earthquake;(b)Dayao MS6.1 earthquake

4 結(jié)論與討論

震源機(jī)制解滑動(dòng)方向擬合法結(jié)果不僅具有實(shí)際物理意義，并且是定量化的，避免了利用震源機(jī)制解的P、T軸分析帶來的系統(tǒng)偏差。另外，由于震源機(jī)制解真實(shí)反映了深部地殼現(xiàn)今應(yīng)力狀況和斷層構(gòu)造的運(yùn)動(dòng)特性，因而其結(jié)果可能比其它資料得到的結(jié)果可靠性更高。

由大姚M6.2、M6.1地震余震震源機(jī)制解，反演得到的大姚兩次地震震源區(qū)的區(qū)域構(gòu)造應(yīng)力場結(jié)果顯示:大姚M6.2地震震源區(qū)最大主壓應(yīng)力S1方位角為240°，仰角為28°，相對(duì)大小的R值為0.4;大姚6.1級(jí)地震震源區(qū)最大主壓應(yīng)力S1方位角為250°，仰角為6°，相對(duì)大小的R值為0.6。大姚6.1與6.2級(jí)地震S1方位角和仰角相比，大姚M6.1地震 S1方位角相差10°，仰角相差22°，R值相差0.2。這說明兩次地震發(fā)生后構(gòu)造應(yīng)力場發(fā)生了微小改變，震源區(qū)構(gòu)造應(yīng)力場最大主壓應(yīng)力S1作用仍為北東東向水平作用。

此次大姚6.2、6.1級(jí)雙震，震中位置十分靠近，余震的空間分布呈不均勻單側(cè)展布，其余震序列空間分布優(yōu)勢方向北西西—南東東的展布方向與震源機(jī)制節(jié)面B、地震的等震線長軸走向基本一致，表明兩次地震受控于同一發(fā)震斷層。大姚6.1級(jí)地震余震序列位于6.2級(jí)地震時(shí)斷層破裂不充分的東南端，表明大姚6.1級(jí)地震是6.2級(jí)地震在同一斷層向東南端延伸的另一次破裂。地震庫侖破裂靜應(yīng)力變化的計(jì)算結(jié)果表明，7月21日大姚6.2級(jí)地震對(duì)10月16日6.1級(jí)地震有明顯的觸發(fā)作用，這種觸發(fā)作用可能與兩次地震受控于同一斷裂相關(guān)。

崔效鋒，謝富仁.1999.利用震源機(jī)制解對(duì)中國西南及鄰區(qū)進(jìn)行應(yīng)力分區(qū)的初步研究［J］.地震學(xué)報(bào)，21(5):513－522.

胡家富，蘇有錦，朱雄關(guān)，等.2003.云南的地殼S波速度與泊松比結(jié)構(gòu)及其意義［J］.中國科學(xué)(D輯:地球科學(xué))，33(8):714－712.

華衛(wèi)，劉杰，鄭斯華，等.2006.2003年云南大姚6.2、6.1級(jí)地震序列特征分析及地震觸發(fā)研究［J］.中國地震，22(1):10－23.

闞榮舉，林中洋.1986.云南地殼上地幔構(gòu)造的初步研究［J］.中國地震，2(4):50～61.

劉桂萍，傅征祥.2001.海原大地震對(duì)古浪大地震的靜應(yīng)力觸發(fā)研究［J］.地球物理學(xué)報(bào)，44(z1):107－115.

蘇有錦.2004.2003年7月21日、10月16日云南大姚6.2級(jí)和6.1級(jí)地震預(yù)測預(yù)報(bào)回顧與討論［J］.國際地震動(dòng)態(tài)，(1):18－21.

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

萬永革，吳忠良，周公威，等.2002.地震應(yīng)力觸發(fā)研究［J］.地震學(xué)報(bào)，24(5):533－551.

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

吳建平，明躍紅，王春槦，等.2001.云南數(shù)字地震臺(tái)下方的S波速度結(jié)構(gòu)研究［J］.地球物理學(xué)報(bào)，44(2):228 －237.

吳小平，虎雄林，Michel Bouchon，等.2007.云南瀾滄—耿馬MS7.6地震的完全庫侖破裂應(yīng)力變化與后續(xù)地震的動(dòng)態(tài)、靜態(tài)應(yīng)力觸發(fā)［J］.中國科學(xué)(D 輯:地球科學(xué))，37(6):746－752.

謝富仁，劉光勛，梁海慶.1994.滇西北及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場［J］.地震地質(zhì)，16(4):329－33.

謝富仁，張紅艷，崔效峰，等.2001.中國大陸現(xiàn)代構(gòu)造應(yīng)力場與強(qiáng)震活動(dòng)［J］.國際地震動(dòng)態(tài)，(1):4－12.

徐紀(jì)人，趙志新，石川有三.2008.中國大陸地殼應(yīng)力場與構(gòu)造運(yùn)動(dòng)區(qū)域特征研究［J］.地球物理學(xué)報(bào)，51(3):770－781.

許忠淮，戈澍謨.1984.用滑動(dòng)方向擬合法反演富蘊(yùn)地震斷裂帶應(yīng)力場［J］.地震學(xué)報(bào)，6(4):395－404.

許忠淮.1985.用滑動(dòng)方向擬合法反演唐山余震區(qū)的平均應(yīng)力場［J］.地震學(xué)報(bào)，7(4):349－361.

張國民，李麗.1997.強(qiáng)震成組孕育、成組發(fā)生過程中相互間影響的研究［J］.地震，17(3):221－231.

張建國，劉麗芳，李西，等.2009.姚安、大姚中強(qiáng)震區(qū)地震構(gòu)造初析［J］.地震地質(zhì)，31(3):536 －540.

周龍泉，劉杰，張曉東，等.2007.2003年大姚6.2和6.1級(jí)地震前三維波速結(jié)構(gòu)的演化［J］.地震學(xué)報(bào)，29(1):20－30.

周真恒，向方英，楊昆杰.2001.云南地區(qū)地殼和上地幔的巖石學(xué)結(jié)構(gòu)［J］.地震地質(zhì)，21(1):69－78.

Angelier J.1979.Determination of the mean principal direction of stresses for a given fault population［J］.Tectonophysics，56(3 －4):17 － 26.

Brace W F，Kohlstedt D.1980.Limits on lithospheric stress imposed by laboratory experiments［J］.JGR，85(B11):6 248 － 6 252.

Harris R A.1998.Introduction to special section:stress triggers，stress shadows，and implications for seismic hazard［J］.JGR，103(10):24 347－24 358.

King G，Stein R S，Lin J.1994.Static stress changes and the triggering of earthquake［J］.BSSA，84(3):935 －953.

Okada Y.1992.Internal deformation due to shear and tensile faults in a half-space［J］.BSSA，82(2):1018 －1040.

Robinson R，McGinty P J.2000.The enigma of the Arthur's Pass，New-Zealand，earthquake 2:the aftershock distribution and its relation to regional and induced stress fields［J］.JGR，105(B7):16 139 ～16 150.

Zoback M l.1992.First-and second-order patterns of stress in the lithosphere:the world stress map project［J］.JGR，97(B8):11 703 －11 728.