庫(kù)水載荷對(duì)水庫(kù)觸發(fā)地震的機(jī)制研究＊

王秋月 朱守彪

(中國(guó)地震局地殼應(yīng)力研究所,北京 100085)

庫(kù)水載荷對(duì)水庫(kù)觸發(fā)地震的機(jī)制研究＊

王秋月 朱守彪

(中國(guó)地震局地殼應(yīng)力研究所,北京 100085)

利用有限元方法,研究庫(kù)水載荷產(chǎn)生的應(yīng)力場(chǎng)及其在斷層面上引起的靜態(tài)庫(kù)侖破裂應(yīng)力變化。通過計(jì)算發(fā)現(xiàn),庫(kù)水載荷的加載位置對(duì)地震的觸發(fā)作用有重要的影響：當(dāng)庫(kù)水加載于斷層正上方或斷層上盤時(shí),對(duì)正斷層型的地震有促進(jìn)作用,對(duì)逆斷層型的地震有抑制作用。但是,當(dāng)庫(kù)水載荷作用于斷層的下盤時(shí),對(duì)于傾角較大的正斷層與逆斷層地震,庫(kù)水都起抑制作用;而對(duì)傾角很小的斷層起促進(jìn)作用,但不會(huì)產(chǎn)生很強(qiáng)的水庫(kù)地震。

水庫(kù)地震;有限元;庫(kù)水載荷;靜態(tài)庫(kù)侖應(yīng)力變化;斷層

1 引言

水庫(kù)在造福人類的同時(shí),也改變了周邊環(huán)境,甚至誘發(fā)地震。自從 Carder[1]在 1945年指出人工水庫(kù)蓄水可能誘發(fā)地震以來,水庫(kù)觸發(fā)地震已經(jīng)在全球范圍內(nèi)被發(fā)現(xiàn)[2]。1962年中國(guó)的新豐江,1963年贊比亞和津巴布韋邊界處的 Kariba,1966年希臘的Kremasta,1967年印度的 Koyna都發(fā)生過由于水庫(kù)運(yùn)行誘發(fā)的震級(jí)超過 6級(jí)的中強(qiáng)地震。其中,Koyna大壩建成后,從 1962年開始蓄水,地震就開始不斷發(fā)生,并且頻度大,強(qiáng)度高,1967年的 Koyna地震達(dá)到 6.3級(jí),是截至目前為止最強(qiáng)烈也是危害最嚴(yán)重的水庫(kù)誘發(fā)地震,致使約 200人喪生,1 500人受傷,數(shù)以千計(jì)的人無家可歸[2]。

此后,水庫(kù)誘發(fā)地震得到普遍關(guān)注,國(guó)內(nèi)外許多學(xué)者對(duì)水庫(kù)誘發(fā)地震的特征及機(jī)制進(jìn)行了研究[1-10],通過大量震例分析認(rèn)為水庫(kù)誘發(fā)地震的主要原因是：庫(kù)水載荷作用,孔隙壓變化,水對(duì)巖體及斷層的腐蝕、軟化、潤(rùn)滑作用等。實(shí)際上,水庫(kù)地震發(fā)生的機(jī)制非常復(fù)雜,不同的水庫(kù),由于其庫(kù)盆的構(gòu)造應(yīng)力環(huán)境、水文地質(zhì)條件、巖石及斷層性質(zhì)等的差異,影響地震發(fā)生的主要因素不同,并且隨著時(shí)間不斷變化。如在蓄水初期,庫(kù)水載荷可能起主要作用;但隨著水位變化、水在巖石裂隙中的不斷滲流,孔隙壓的改變對(duì)誘發(fā)地震可能起著決定性作用。

由于水庫(kù)誘發(fā)地震的機(jī)制十分復(fù)雜,為突出主要問題,本文只討論庫(kù)水載荷的觸發(fā)作用。

不少學(xué)者利用有限元方法,計(jì)算了庫(kù)水載荷所產(chǎn)生的應(yīng)力場(chǎng)分布及其對(duì)地震的觸發(fā)[11-14],他們的研究成果對(duì)認(rèn)識(shí)水庫(kù)觸發(fā)地震的機(jī)理起到一定的積極作用。但由于庫(kù)區(qū)及周邊地區(qū)初始應(yīng)力狀態(tài)不為所知,利用計(jì)算的應(yīng)力來分析斷層的運(yùn)動(dòng)狀態(tài),進(jìn)而推測(cè)地震風(fēng)險(xiǎn)往往遇到很多困難。

數(shù)十年來,利用斷層面上的庫(kù)侖破裂應(yīng)力變化(ΔCFS)來分析地震危險(xiǎn)性已成為國(guó)際地震學(xué)研究的熱點(diǎn)[15-19],并成功地運(yùn)用于許多地震活動(dòng)性的研究之中。本文將嘗試?yán)糜邢迒卧治龇椒?計(jì)算庫(kù)水的應(yīng)力場(chǎng)及其在斷層面上產(chǎn)生的庫(kù)侖破裂應(yīng)力變化,由此來探索庫(kù)水載荷觸發(fā)地震的規(guī)律。

2 庫(kù)水產(chǎn)生的應(yīng)力場(chǎng)及位移場(chǎng)

2.1 有限元模型及邊界條件

考慮一個(gè)狹長(zhǎng)型水庫(kù),假定水庫(kù)水面寬為 20 km,斷層走向與庫(kù)體方向一致,現(xiàn)取庫(kù)體中的一橫截面來進(jìn)行研究。為盡量減小邊界效應(yīng)對(duì)計(jì)算結(jié)果的影響,文中建立的有限元模型比實(shí)際的尺寸要大得多。圖 1為有限元模型的幾何形狀、網(wǎng)格剖分及邊界條件分布圖。圖中模型的長(zhǎng)取 200 km,深為100 km,水載荷范圍為 20 km(圖中上方黑色箭頭所示),水頭取 50 m;模型中共有三角形單元 2 808個(gè),節(jié)點(diǎn) 1 476個(gè);介質(zhì)的楊氏模量取 4.0×1010Pa,泊松比為 0.25。

模型的邊界條件選取如下：水載荷作用在上表面20 km的范圍內(nèi)(載荷大小為50 m水深產(chǎn)生的壓力),上表面其他地方自由;底面的垂向位移為零,水平方向自由;兩側(cè)邊界的水平方向位移為零,垂向自由。計(jì)算中采用平面應(yīng)變的力學(xué)模型。

圖1 有限元模型及網(wǎng)格劃分Fig.1 Finite elementmodel and its grid separating

2.2 載荷應(yīng)力場(chǎng)及位移場(chǎng)

圖 2 庫(kù)水載荷產(chǎn)生的應(yīng)力場(chǎng)等值線分布圖Fig.2 Contours of stress field caused by the load of reservoirwater

圖 2是庫(kù)水載荷產(chǎn)生的應(yīng)力場(chǎng)分量的等值線分布圖。圖 2(a)是水平向應(yīng)力 (σx)的等值線分布,由圖 2(a)可見,水平向應(yīng)力主要集中在庫(kù)水作用的范圍,且其深度影響范圍很小 (主要在深度小于15 km的范圍內(nèi))。圖 2(b)為垂直方向應(yīng)力 (σy)等值線分布圖,圖 2(b)顯示,垂向應(yīng)力影響的深度比水平向大得多。圖 2(c)為剪應(yīng)力 (τxy)等值線分布圖,由圖 2(c)可見,剪應(yīng)力基本呈對(duì)稱分布。

圖 3為庫(kù)水載荷產(chǎn)生的位移場(chǎng)等值線分布圖。圖 3顯示,庫(kù)水產(chǎn)生的位移場(chǎng)為對(duì)稱分布,在庫(kù)水作用中心最大,最大值約0.5 m。但在水平向及垂向迅速衰減,當(dāng)超過庫(kù)水作用距離約 60 km后,無論水平向還是垂向其位移都趨向于零。

圖 3 庫(kù)水載荷產(chǎn)生的位移場(chǎng)等值線分布圖Fig.3 Contours of displacements caused by the load of reservoirwater

3 斷層面上庫(kù)侖破裂應(yīng)力變化

3.1 庫(kù)侖破裂應(yīng)力變化

通過有限元數(shù)值計(jì)算得到庫(kù)水載荷產(chǎn)生的應(yīng)力場(chǎng)后,地下任意一點(diǎn)的應(yīng)力狀態(tài)變化也就唯一確定。由此就可以計(jì)算具有任一傾角斷層面上的正應(yīng)力與剪應(yīng)力變化,進(jìn)而計(jì)算此斷層面上的庫(kù)侖破裂應(yīng)力變化。設(shè)斷層面的傾角為θ,則其上變化的正應(yīng)力與剪應(yīng)力分別為[20]:

其中,σx、σy、τxy分別為 x、y向及剪應(yīng)力分量。根據(jù)庫(kù)侖定律,區(qū)域內(nèi)某一斷層面上沿滑動(dòng)方向的庫(kù)侖破裂應(yīng)力變化 (ΔCFS)可以寫成[15,21,22]:

其中,Δp是孔隙壓變化。通過一些假定與簡(jiǎn)化處理,式 (2)可以進(jìn)一步簡(jiǎn)化為[21-24]:

上式中μ’為視摩擦系數(shù),本文計(jì)算中其值取為0.4[21-23]。

3.2 斷層面上的庫(kù)侖應(yīng)力變化

若斷層面上的庫(kù)侖破裂應(yīng)力變化為正時(shí),庫(kù)水將促進(jìn)斷層滑動(dòng),容易產(chǎn)生地震;庫(kù)侖破裂應(yīng)力變化為負(fù)值時(shí),庫(kù)水會(huì)抑制地震的發(fā)生。為了研究庫(kù)侖應(yīng)力變化與加載位置、斷層類型、斷層位置的關(guān)系,文中讓斷層類型、斷層傾角發(fā)生變化,分 6種情況進(jìn)行研究。這 6種情況分別是：1)水庫(kù)位于正斷層正上方;2)水庫(kù)位于逆斷層正上方;3)水庫(kù)位于正斷層上盤;4)水庫(kù)位于逆斷層上盤;5)水庫(kù)位于正斷層下盤;6)水庫(kù)位于逆斷層下盤。下面將計(jì)算在每種情況下,當(dāng)斷層傾角分別取 53°、45°、37°、25°、14°時(shí),分別研究斷層面上的庫(kù)侖破裂應(yīng)力變化分布情況(圖 4～6中,從左向右,斷層傾角的變化均相同,下文同)。

圖 4 水庫(kù)位于斷層正上方時(shí),斷層面上的庫(kù)侖破裂應(yīng)力變化分布Fig.4 Static Coulomb stress changes on the specified receiver faultswhen the reservoir locates right over the faults

圖 4表示水庫(kù)位于斷層正上方時(shí),斷層面上的庫(kù)侖破裂應(yīng)力變化分布。圖 4中方框尺度為長(zhǎng)度60 km,深度為 50 km(即圖 1中方框表示的地區(qū));圖 4上方 5幅圖表示的是第 1)種情況,即斷層為正斷層,其中圖 4(a)～ (e)分別對(duì)應(yīng)斷層傾角為 53°、 45°、37°、25°、14°的情況;下方 5幅圖表示第 2)種情況,即斷層為逆斷層。

圖 4清楚地顯示,對(duì)于正斷層,庫(kù)水載荷產(chǎn)生的庫(kù)侖破裂應(yīng)力變化基本上都是正值,即庫(kù)水對(duì)正斷層型的地震有觸發(fā)作用,但當(dāng)斷層深度很淺、傾角很小時(shí),觸發(fā)作用不明顯。對(duì)于逆沖型斷層,圖 4結(jié)果顯示,庫(kù)水對(duì)地震的發(fā)生全是抑制作用。若某一地區(qū),庫(kù)區(qū)本底地震機(jī)制為逆沖型,在水庫(kù)蓄水后,庫(kù)水作用會(huì)使該地區(qū)的地震的活動(dòng)性降低。

圖 5顯示的是水庫(kù)位于正斷層上盤與水庫(kù)位于逆斷層上盤的情況。由圖 5(a)～(e)可知,當(dāng)庫(kù)水作用于正斷層上盤時(shí),位于庫(kù)岸的正斷層面上的庫(kù)侖應(yīng)力變化大于零,即庫(kù)水有促使地震發(fā)生的作用,但隨著傾角的減小、距離庫(kù)岸距離的變遠(yuǎn),觸發(fā)能力越來越弱。圖 5(f)～(i)是庫(kù)水作用于逆斷層上盤的情況,由圖 5(f)～(i)可見,對(duì)于深度小于～15 km的逆斷層,其斷層面的庫(kù)侖應(yīng)力變化小于零,即庫(kù)水對(duì)地震有抑制作用,但當(dāng)斷層深度大于 15 km后,庫(kù)水對(duì)庫(kù)岸的逆斷層型地震有促進(jìn)作用。同時(shí)由圖 5(f)～(i)還可以看出,對(duì)于庫(kù)岸的逆斷層,不管傾角如何,只要深度較淺,庫(kù)水對(duì)其都有抑制作用。由此可見,當(dāng)庫(kù)水作用于斷層上盤時(shí),庫(kù)岸很難發(fā)生震源深度較小的逆斷層型地震。

圖 5 水庫(kù)位于斷層上盤時(shí),斷層面上的庫(kù)侖破裂應(yīng)力變化Fig.5 Static Coulomb stress changes on faultswhen the reservoir locates on the hangingwall

圖 6 庫(kù)水位于斷層下盤時(shí),斷層面上庫(kù)侖破裂應(yīng)力的變化Fig.6 Static Coulomb stress changes for specified receiver faultwhen the reservoir locates on the footwall

圖 6顯示的是庫(kù)水作用于斷層下盤的情況。6 (a)～(e)表示正斷層的情況,6(f)～(i)表示逆斷層情況。由圖 6(a)～(e)可見,當(dāng)斷層傾角較大時(shí),庫(kù)水載荷無論對(duì)正斷層還是逆斷層有抑制作用。但隨著斷層傾角的減小,庫(kù)水載荷對(duì)地震作用的性質(zhì)會(huì)發(fā)生變化,由抑制向促進(jìn)轉(zhuǎn)變;當(dāng)傾角小到一定程度后,會(huì)促進(jìn)地震的發(fā)生,但這時(shí)發(fā)生的地震,其震源深度會(huì)很小,因此地震的震級(jí)也會(huì)很小。由此可知,當(dāng)庫(kù)水作用在斷層下盤時(shí),不可能觸發(fā)震級(jí)很高的地震。

4 討論與結(jié)論

利用有限元方法,計(jì)算了庫(kù)水載荷產(chǎn)生的應(yīng)力場(chǎng)、位移場(chǎng)及其對(duì)水庫(kù)下方及庫(kù)岸斷層中地震的觸發(fā)作用。從計(jì)算的結(jié)果看,庫(kù)水產(chǎn)生的庫(kù)侖應(yīng)力變化最大值可達(dá) 0.06 MPa,這大大超過了典型靜態(tài)應(yīng)力觸發(fā)的閾值 0.01 MPa[21-24]。由此可知,庫(kù)水對(duì)地震的觸發(fā)是有效的,這就是我們能夠觀測(cè)到大量的快速響應(yīng)性水庫(kù)觸發(fā)地震的主要原因(如前蘇聯(lián)塔吉克斯坦的 Nurek水庫(kù)地震)[25]。文中的水頭取為 50m,實(shí)際上在蓄水初期,庫(kù)水在斷層面產(chǎn)生的庫(kù)侖應(yīng)力變化有可能超過斷層破裂極限,形成地震[2,26];隨著水庫(kù)蓄水水位的不斷升高,庫(kù)水產(chǎn)生的庫(kù)侖破裂應(yīng)力變化不斷增大,可能會(huì)觸發(fā)更多、更大的地震。但也有不少水庫(kù)蓄水水位很高,但沒有觸發(fā)地震的情況,這可能是當(dāng)?shù)貥?gòu)造應(yīng)力場(chǎng)還遠(yuǎn)離斷層失穩(wěn)極限。因此庫(kù)水能否觸發(fā)地震一方面與庫(kù)水產(chǎn)生的庫(kù)侖應(yīng)力變化有關(guān),同時(shí)還取決于當(dāng)?shù)氐臉?gòu)造應(yīng)力狀態(tài)。

從計(jì)算結(jié)果看,庫(kù)水載荷不僅對(duì)水庫(kù)底部的斷層有觸發(fā)作用,對(duì)庫(kù)岸斷層也有影響,所以我們看到很多水庫(kù)地震不僅發(fā)生在庫(kù)底,而且也發(fā)生在庫(kù)岸[4-10]。

從圖 4～6可知,一般情況下,庫(kù)水對(duì)正斷層型的地震有觸發(fā)作用,因此正斷層型的水庫(kù)地震較多。如：1963年的卡里巴地震,1966年的科列馬斯塔地震,1975年的奧洛維爾地震等都為正斷層型地震[28];王墩等[29]結(jié)合我國(guó)長(zhǎng)江三峽地區(qū)發(fā)生水庫(kù)地震的具體情況,選取長(zhǎng)江三峽數(shù)字遙測(cè)臺(tái)記錄的水庫(kù) 2000—2006年蓄水后的波形數(shù)據(jù),得到水庫(kù)重點(diǎn)監(jiān)視區(qū)內(nèi)庫(kù)區(qū)震源機(jī)制解和庫(kù)岸小震綜合節(jié)面解。計(jì)算結(jié)果顯示：庫(kù)區(qū) 25次水庫(kù)地震震源機(jī)制解除兩次情況為無解外,斷層性質(zhì)多為正斷層,其次為走滑型斷層,僅有 3次為逆斷層情況,P軸仰角較大,P軸方位、T軸方位、水平最大主壓應(yīng)力方向均雜亂無章;巴東庫(kù)段沿江 5 km范圍內(nèi) 1≤ML≤2.4小震綜合節(jié)面解顯示庫(kù)岸地震機(jī)制 80%為正斷層,其余為其他情況。這表明文中的計(jì)算結(jié)果與實(shí)際情況較為一致,對(duì)實(shí)際地震的分析預(yù)測(cè)有著現(xiàn)實(shí)意義。

庫(kù)水對(duì)逆斷層型地震起抑制作用。但對(duì)庫(kù)岸逆斷層型的地震,當(dāng)深度達(dá)到一定范圍時(shí)也有觸發(fā)作用。

當(dāng)庫(kù)水在斷層下盤加載時(shí),無論對(duì)正斷層還是逆斷層都起到抑制地震發(fā)生的作用。但是當(dāng)斷層傾角很小時(shí),庫(kù)水載荷對(duì)地震發(fā)生也有促進(jìn)作用。通常這種情況下產(chǎn)生的水庫(kù)地震,震級(jí)不會(huì)太大。因此,我們?cè)诮ㄔ焖畮?kù)時(shí),若斷層無法回避,盡量讓庫(kù)水加載于斷層下盤,也有利于避免庫(kù)水觸發(fā)地震。這對(duì)于大型水庫(kù)建設(shè)有重要的參考價(jià)值,希望引起有關(guān)方面的重視。

此外,文中沒有討論庫(kù)水載荷對(duì)走滑型斷層的影響,從力學(xué)分析來看,庫(kù)水載荷不會(huì)對(duì)走滑型地震產(chǎn)生影響[30]。因此,我們看到的走滑型水庫(kù)地震不是庫(kù)水載荷觸發(fā)的,而是由于庫(kù)水的滲流、軟化、孔隙壓的變化等誘發(fā)產(chǎn)生的。作者正在進(jìn)行這方面的深入研究。

水庫(kù)水位、斷層的性質(zhì)及其位置與庫(kù)水載荷觸發(fā)地震的關(guān)系非常復(fù)雜。能否根據(jù)這一現(xiàn)象,選擇合適的地區(qū),建立一座中型水庫(kù),作為水庫(kù)地震研究的實(shí)驗(yàn)室。通過人為調(diào)節(jié)水位(升高或降低及水位變化的速率),人為地控制地震的發(fā)生 (強(qiáng)度、數(shù)量、地點(diǎn))來深入研究水庫(kù)觸發(fā)地震及地震預(yù)測(cè)問題,為成功實(shí)現(xiàn)地震預(yù)報(bào)尋找突破口。

綜合文中的結(jié)果,我們得出以下幾點(diǎn)初步認(rèn)識(shí)：

1)庫(kù)水載荷可以改變水庫(kù)底部及庫(kù)岸的應(yīng)力場(chǎng),并導(dǎo)致斷層面上的庫(kù)侖破裂應(yīng)力變化,從而觸發(fā)地震。

2)一般來說,庫(kù)水載荷對(duì)正斷層型地震有觸發(fā)作用,對(duì)逆斷層型地震起抑制作用;但對(duì)庫(kù)岸深度較深的逆斷層也有觸發(fā)作用。

3)當(dāng)庫(kù)水載荷作用于斷層的下盤時(shí),無論是正斷層還是逆斷層地震,都起抑制作用。

1 Carder D S.Seismic investigations in Boulder dam area, 1940-1944,and the influence of reservoir loading on earthquake activity[J].Bull. Seis. Soc.Am.,1945, (35)：175-192.

2 Gupta and Harsh K.A review of recent studies of triggered earthquakes by artificialwater reservoirswith special emphasis on earthquakes in Koyna,India[J].Earth-Science Reviews,2002,(58)：279-310.

3 Grouph D I and Gouph W I.Load-induced earthquakes at Lake Kariba[J].Geophys.J.,1970,21：79-101.

4 王妙月,等.新豐江水庫(kù)地震的震源機(jī)制及其成因的初步探討[J].中國(guó)科學(xué),1976,1：85-97.

5 Roeloffs E A.Fault stability changes induced beneath a reservoirwith cyclic variations in water level[J].J.Geophys. Res.,1988,93(B3)：2 107-2 124.

6 Reyners M.Reservoir induced seismicity at Lake Pukaki, New Zealand[J].Geophys.J.,1988,93：127-135.

7 高士鈞,等.長(zhǎng)江三峽地區(qū)地殼應(yīng)力場(chǎng)與地震[M].北京：地震出版社,1992.

8 Hu Y L,et al.Induced seis micity atWujiangdu Reservoir, China：a case induced in KarstArea[J].Pure Appl.Geophys.,1996,147：409-418.

9 王清云,張秋文,李峰.長(zhǎng)江三峽工程庫(kù)首區(qū)誘發(fā)地震危險(xiǎn)性研究[J].大地測(cè)量與地球動(dòng)力學(xué),2003,(2)：101 -104.

10 陳蜀俊,姚運(yùn)生,曾佐勛.三峽水庫(kù)蓄水對(duì)庫(kù)區(qū)孕震環(huán)境及潛在震源影響研究[J].大地測(cè)量與地球動(dòng)力學(xué), 2005,(3)：116-120.

11 劉素梅,徐禮華.丹江口水庫(kù)區(qū)構(gòu)造應(yīng)力場(chǎng)的數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報(bào),2004,23(23)：4 017-4 021.

12 王紹晉,秦嘉政,龍曉帆.漫灣水庫(kù)蓄水前后庫(kù)區(qū)地震活動(dòng)性與構(gòu)造應(yīng)力場(chǎng)分析[J].地震研究,2005,28(1)：53-57.

13 齊文浩,薄景山,李貞德.江口水庫(kù)誘發(fā)地震的有限元分析[J].地震工程與工程振動(dòng),2006,26(6)：7-10.

14 李錦輝.庫(kù)區(qū)構(gòu)造應(yīng)力場(chǎng)對(duì)水庫(kù)誘發(fā)地震的影響[J].武漢大學(xué)學(xué)報(bào) (工學(xué)版),2007,40(6)：51-55.

15 Harris R A.Introduction to special session：Stress triggers, stress shadows,and implications for seismic hazard[J].J. Geophys.Res.,1998,103：24 347-24 358.

16 石耀霖.關(guān)于應(yīng)力觸發(fā)和應(yīng)力影概念在地震預(yù)報(bào)中應(yīng)用的一些思考[J].地震,2001,21(3)：1-7.

17 BelardinelliM E,BizzarriA and CoccoM.Earthquake triggering by static and dynamic stress changes[J].J.Geophys. Res., 2003, 108 (B3), doi：10. 1029/ 2002JB001779.

18 Steacy S,Gomberg J and CoccoM.Introduction to special section：Stress transfer,earthquake triggering,and time-dependent seismic hazard[J].J.Geophys.Res.,2005,110, doi：10.1029/2005JB003692.

19 Parsons T,Chen Ji and Kirby E.Stress changes from the 2008Wenchuan earthquake and increased hazard in the Sichuan basin[J].Nature, 2008, doi：10.1038/nature07177.

20 徐芝綸.彈性力學(xué)[M].北京：高等教育出版社,2008.

21 Stein R S,Barka A A and Dieterich J H.Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering[J].Geophys.J. Int.,l997,128：594-604.

22 Reasenberg P A and Simpson R W.Response of regional seis micity to the static stress change produced by the Loma Prieta earthquake[J].Science,1992,255：168721 690.

23 Stein R S,King G C P and Lin J.Change in failure stress on the southern San Andreas fault system caused by the 1992 magnitude=7.4 Landers earthquake[J].Science, 258：1 328-1 332.

24 King G C P,Stein R S andLin J.Static stress changes and the triggering of earthquakes[J].Bull.Seis m.Soc.A-mer.,1994,84：935-953.

25 Simpson D W and Negmatullaev S.Induced seis micity at Nurek reservoir,Tadjikistan,USSR[J].Bull. Seismol. Soc.Am.,1981,71：1 561-1 586.

26 Si mpson D W.Seismicity changes associated with reservoir loading[J].Eng.Geol.,1976,10：2-4.

28 張倬元.工程地質(zhì)分析原理 [M].北京：地質(zhì)出版社, 2005.

29 王墩,等.三峽水庫(kù)重點(diǎn)監(jiān)視區(qū)蓄水前后震源機(jī)制解研究[J].大地測(cè)量與地球動(dòng)力學(xué),2007,(5)：103-107.

30 Scholz C H.The mechanics of earthquakes and faulting [M].Cambridge University Press,1990.

STUDY ON TRIGGERING M ECHANISM OF RESERVO IR EARTHQUAKE BY WATER IM POUNDM ENT

WangQiuyue and Zhu Shoubiao
(Institute of Crustal Dynam ics,CEA,B eijing 100085)

Using finite elementmethod.The stresse field and the static Coulomb failure stress changeswas calculated.The calculated results show that the location of the water impoundmentplays an important role in triggering reservoir earthquake.When the reservoirwater is loaded right over the fault and on the hangingwall of the fault,it will promote the earthquakes for normal faults and inhibit the seis micity for thrust faults.However,the earthquakes will all be inhibited both for normal and thrust faultswith high dip angleswhen the load of reservoirwater is located on the footwall of the faults.No strong reservoir-induced earthquakeswill occur despite the load of waterwill promote the earthquake for the faultswith small dip angleswhen the load is put on the footwall of the faults.

reservoir-induced earthquake;finite elementmodel;load of reservoirwater;static Coulomb failure stress change;fault

1671-5942(2010)Supp.(Ⅰ)-0017-06

2010-02-23

國(guó)家自然科學(xué)基金(40974020,40774024);北京市自然科學(xué)基金(8082024);地震動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室基金(LED2008B02)

王秋月,女,1982年生,碩士,主要從事地下流體及地球動(dòng)力學(xué)研究.E-mail：wangqiuyue09@sina.com

P315.72+8

A