周 斌， 李 蕾， 孫 峰， 史水平

(1.廣西壯族自治區(qū)地震局，廣西 南寧 530022；2.中國(guó)石油大學(xué)（華東），山東 青島 266580）

0 引言

水庫(kù)誘發(fā)地震（Reservoir-induced seismicity，簡(jiǎn)稱RIS）是由于水庫(kù)蓄水或排水過程引發(fā)的在一定時(shí)間內(nèi)庫(kù)區(qū)及其周邊不大區(qū)域范圍內(nèi)所發(fā)生的地震活動(dòng)[1]。迄今為止，全球已報(bào)道的RIS確切震例約有130例，其中，大于6.0級(jí)的有4例，5.0～5.9級(jí)的有15例，其余占85%以上的均小于5.0級(jí)。雖然大部分RIS的震級(jí)不高，但由于其震中位置一般鄰近重要的水利工程設(shè)施，且震源淺、震中烈度高，往往具有很大的破壞性，可造成大壩及附近建筑物的破壞和人員傷亡，甚至?xí)鸹?、坍塌等?yán)重的次生災(zāi)害危及下游安全。從全球已發(fā)生的RIS震例來看，雖然發(fā)震水庫(kù)占大中型水庫(kù)總數(shù)的比例較少，但在全球分布比較廣泛，構(gòu)造地震活動(dòng)比較頻繁的地區(qū)和相對(duì)平靜的少震、弱震地區(qū)都有發(fā)生。為了進(jìn)一步研究庫(kù)水加卸載及滲透過程中，RIS的動(dòng)態(tài)響應(yīng)的力學(xué)機(jī)制，對(duì)已觀測(cè)到的確切震例的時(shí)空活動(dòng)特征進(jìn)行多方面的了解和認(rèn)識(shí)是十分必要的。此外，RIS在許多方面具有與天然構(gòu)造地震不同的特征，了解RIS的空間分布規(guī)律、誘發(fā)地震活動(dòng)與水庫(kù)蓄、排水過程中水位變化之間的響應(yīng)關(guān)系，特別是RIS時(shí)空分布與庫(kù)區(qū)地質(zhì)構(gòu)造環(huán)境關(guān)系，是RIS機(jī)理與預(yù)測(cè)研究最基本也是最重要的條件。本文廣泛收集了國(guó)內(nèi)外已確定的RIS震例資料，利用統(tǒng)計(jì)學(xué)的方法，開展了此方面的研究。

1 RIS活動(dòng)的空間分布特征

1.1 RIS區(qū)域分布特征及其與地震活動(dòng)背景的關(guān)系

在RIS研究早期，一些學(xué)者就試圖從誘發(fā)地震與天然地震活動(dòng)性的關(guān)系、誘發(fā)地震的新構(gòu)造環(huán)境、 地 殼應(yīng)變 速 率等 方 面闡釋 RIS的 成 因[2、3]。Castle等[4]對(duì)38個(gè)RIS震例進(jìn)行了統(tǒng)計(jì)研究，在33個(gè)精確定位的震例中，發(fā)生在無震區(qū)的有10個(gè)，低震區(qū)11個(gè)，低－中等震區(qū)2個(gè)，中等震區(qū)4個(gè)，中等－強(qiáng)震區(qū)4個(gè)，強(qiáng)震區(qū)2個(gè)。筆者對(duì)全球99例MS≥3.0級(jí)的RIS震例進(jìn)行了統(tǒng)計(jì)，有80個(gè)震例可以確定其所在區(qū)域的天然地震活動(dòng)水平，其中，發(fā)生在無震區(qū)的有24個(gè)，弱震區(qū)23個(gè)，中等震區(qū)8個(gè)，強(qiáng)震區(qū)24個(gè)，距強(qiáng)震活動(dòng)區(qū)12 km的1個(gè)。

從全球130例水庫(kù)誘發(fā)地震的震中分布情況可以看出，RIS的震中位置并不是都分布于地震活動(dòng)構(gòu)造帶上。印度半島是世界上水庫(kù)誘發(fā)地震較多的地區(qū)之一，著名的柯依納（Koyna）MS6.5級(jí)破壞性水庫(kù)誘發(fā)地震就發(fā)生在這一地區(qū)。在區(qū)域地震構(gòu)造環(huán)境上，印度半島屬前寒武紀(jì)地盾，是世界上最穩(wěn)定的地質(zhì)構(gòu)造區(qū)之一。根據(jù)印度歷史記載和現(xiàn)代地震記錄，印度的強(qiáng)震活動(dòng)主要分布在北部和西部地區(qū)，而在1967年12月10日柯依納水庫(kù)誘發(fā)地震以前，印度半島極少有地震活動(dòng)[5]。

北美的水庫(kù)誘發(fā)地震，絕大多數(shù)是發(fā)生于區(qū)域地震活動(dòng)水平低的地區(qū)。例如喬卡茜 （Jocasse）、蒙蒂塞洛（Monticello）等發(fā)震水庫(kù)集中分布的南卡羅來納州，歷史地震活動(dòng)性處于較低的水平。加拿大的發(fā)震水庫(kù)馬尼克-3（Manicouagan-3）和LG-3，蓄水前庫(kù)區(qū)均無有感地震記錄[6]。

一些著名的水庫(kù)誘發(fā)地震，如位于贊比亞－津巴布韋邊界的卡里巴水庫(kù)（Kariba）、埃及阿斯旺水庫(kù)（Aswan）、塔吉克斯坦的努列克水庫(kù)（Nurek）等，蓄水后誘發(fā)地震的震級(jí)和烈度均超過當(dāng)?shù)貧v史最高水平。

目前，我國(guó)已經(jīng)發(fā)生的RIS震例約有46個(gè)。在這些震例中，發(fā)生于強(qiáng)震構(gòu)造區(qū)的有10個(gè)，包括：新疆強(qiáng)震構(gòu)造區(qū)的克孜爾水庫(kù)和恰甫其海水庫(kù)，青藏高原強(qiáng)震構(gòu)造區(qū)的紫坪鋪水庫(kù)、漫灣水庫(kù)、大橋水庫(kù)、盛家峽水庫(kù)、二灘水庫(kù)、茄子山水庫(kù)和小灣水庫(kù)，臺(tái)灣強(qiáng)震構(gòu)造區(qū)的曾文水庫(kù)；發(fā)生于中等強(qiáng)度地震構(gòu)造區(qū)的有7個(gè)，分別是龍羊峽水庫(kù)、李家峽水庫(kù)、新店水庫(kù)、銅子街水庫(kù)、佛子嶺水庫(kù)、云鵬水庫(kù)和參窩水庫(kù)；其余的29個(gè)震例基本上分布在華南弱震區(qū)。

以上統(tǒng)計(jì)結(jié)果似乎表明，水庫(kù)誘發(fā)地震與區(qū)域地震活動(dòng)性并無密切的相關(guān)性，相對(duì)于天然地震活動(dòng)背景而言，水庫(kù)誘發(fā)地震的發(fā)震比率在缺震區(qū)和弱震區(qū)反而占有一定的優(yōu)勢(shì)。

1.2 RIS震中分布規(guī)律

統(tǒng)計(jì)分析表明，水庫(kù)誘發(fā)地震的震中大多分布于水庫(kù)周緣10 km范圍之內(nèi)，并局限在特定的庫(kù)段叢集分布。例如，我國(guó)的新豐江、大化巖灘、三峽、水口、烏溪江以及國(guó)外的柯依納、阿斯旺等水庫(kù)的震中分布。

構(gòu)造地震的震中遷移是人們所熟知的一種自然地質(zhì)現(xiàn)象，特別是局部地段局部時(shí)間內(nèi)的震中遷移已被廣泛研究。RIS震例統(tǒng)計(jì)表明，當(dāng)某些水庫(kù)區(qū)具備一定的地質(zhì)構(gòu)造與水文地質(zhì)條件時(shí)，誘發(fā)地震活動(dòng)亦會(huì)隨著水庫(kù)蓄、放水過程中的水位變化而發(fā)生遷移。概括起來，RIS遷移的情況大致有4種：一種是在水庫(kù)蓄水初期，隨著水位的快速上升，庫(kù)區(qū)外圍的地震活動(dòng)遷移至庫(kù)壩附近，在已掌握的RIS震例中，此類情況比較普遍；另一種發(fā)生在喀斯特地貌發(fā)育的庫(kù)區(qū)或庫(kù)段，如湖北省丹江口水庫(kù)，該水庫(kù)在蓄水初期，地震主要集中分布在庫(kù)體覆蓋區(qū)周圍，在1971年l1月至1972年4月3日，當(dāng)水位上升達(dá)到150 m后，隨著水庫(kù)淹沒區(qū)范圍的擴(kuò)大，震中逐漸向遠(yuǎn)處的林茂山和宋灣兩個(gè)灰?guī)r峽谷區(qū)遷移，形成了兩個(gè)地震密集區(qū)[7]，這種情況誘發(fā)地震的震源深度一般較淺，且多為巖溶坍塌型。第三種情況跟與水庫(kù)有水力聯(lián)系的斷裂構(gòu)造有關(guān)，往往表現(xiàn)為水庫(kù)經(jīng)歷了長(zhǎng)時(shí)間的蓄水以后，誘發(fā)地震活動(dòng)沿?cái)嗔训淖呦蛳蜻h(yuǎn)離庫(kù)岸的方向遷移（離開庫(kù)緣達(dá)10 km甚至更多），此類誘發(fā)地震一般表現(xiàn)為高震級(jí)、大范圍和有較深的震源深度；第四種情況的震中遷移一般發(fā)生在處于擠壓構(gòu)造應(yīng)力環(huán)境的庫(kù)區(qū)，且受控于水庫(kù)水位的動(dòng)態(tài)變化過程，如塔吉克斯坦的努列克水庫(kù)，該水庫(kù)位于南天山東南部的塔吉克凹陷，由于受?。瓉喆箨懪鲎驳挠绊懀F(xiàn)今構(gòu)造應(yīng)力場(chǎng)以N－NW向擠壓為主。該水庫(kù)蓄水后庫(kù)區(qū)地震活動(dòng)明顯增強(qiáng)，特別是在1972年底（水位上升到105 m）和1976年底（水位上升到215 m）兩個(gè)重要的蓄水階段，隨著水庫(kù)水位的快速抬升，地震震中由水庫(kù)的西南側(cè)沿水庫(kù)軸線向東北側(cè)遷移[8]。我國(guó)四川省紫坪鋪水庫(kù)蓄水后誘發(fā)小震活動(dòng)的遷移現(xiàn)象也屬于此類[9、10]。

RIS的震中遷移特征與水庫(kù)加卸載過程中附加水頭壓力的擴(kuò)散密切相關(guān)，1985年，Talwani等在假定庫(kù)水壓力傳導(dǎo)至震源即發(fā)生地震的基礎(chǔ)上，通過計(jì)算誘震庫(kù)區(qū)介質(zhì)水力傳導(dǎo)率的方法試圖解釋震中遷移的現(xiàn)象[11]，不同研究者以這種方法計(jì)算了20多例誘震庫(kù)區(qū)介質(zhì)的孔隙壓力擴(kuò)散率。周斌等基于地下應(yīng)力場(chǎng)與滲流場(chǎng)耦合理論，以紫坪鋪水庫(kù)為例，探討誘發(fā)地震活動(dòng)時(shí)空演化與庫(kù)水加卸載及滲透過程的關(guān)系，闡釋RIS震中遷移的力學(xué)機(jī)制[9]。

1.3 RIS震源深度分布與遷移

與天然構(gòu)造地震相比，水庫(kù)誘發(fā)地震具有震源深度淺的特征。筆者對(duì)掌握震源深度數(shù)據(jù)的43個(gè)水庫(kù)震例資料做統(tǒng)計(jì)分析，其中，震源深度在0～4.9 km的占震例總數(shù)的46.5%，震源深度在5.0～9.9 km的占震例總數(shù)的34.9%（圖1）。由此可見，水庫(kù)誘發(fā)地震主要發(fā)生在地下10 km以上的脆性上地殼中。

圖1 水庫(kù)誘發(fā)地震震源深度統(tǒng)計(jì)直方與累計(jì)百分比圖Fig.1 Focal depth statistical histogram and cumulative percentage figure of reservoir-induced seismicity

但對(duì)于不同的水庫(kù)，其誘發(fā)地震活動(dòng)的深度范圍亦存在著較大的差異。如，我國(guó)湖北省三峽、貴州省烏江渡、巴西的帕萊布納等水庫(kù)誘發(fā)的溶巖或塌陷型小震群活動(dòng)，其震源深度基本上在2 km以內(nèi)，最淺的僅有幾百米深。我國(guó)廣東省新豐江、廣西壯族自治區(qū)龍灘、巖灘、四川省紫坪鋪以及印度柯依納、希臘克里馬斯塔、贊比亞卡里巴、泰國(guó)斯林那加林、美國(guó)奧洛維爾等絕大多數(shù)水庫(kù)誘發(fā)地震活動(dòng)的深度優(yōu)勢(shì)分布在5～10 km之間。到目前為止，全球記錄到最深誘發(fā)地震的實(shí)例為埃及的阿斯旺水庫(kù)，該水庫(kù)誘發(fā)地震的震源深度定位結(jié)果表明[12]，有深源地震（優(yōu)勢(shì)分布在15～25 km）和淺源地震（優(yōu)勢(shì)分布在 5～10 km）兩類，與世界上其它地區(qū)水庫(kù)誘發(fā)地震的震源深度相比，具有其獨(dú)特的方面。

水庫(kù)誘發(fā)地震的另一個(gè)重要的特征是，震源深度隨水庫(kù)蓄水的時(shí)間歷程有向深部遷移的現(xiàn)象。如廣東省新豐江水庫(kù)，在1961年7月20日至1962年3月19日期間，共記錄到56次MS≥3.0級(jí)地震，其震源深度優(yōu)勢(shì)分布在5.5 km左右，且3 km左右的極淺源地震占有比較大的比例；主震發(fā)生以后，在1963年3月19日至1987年10月3日期間，共記錄到201次MS≥3.0級(jí)地震，此時(shí)震源深度優(yōu)勢(shì)分布在6.5 km左右，3 km左右的極淺源地震減少，而震源深度在8～10 km左右的地震明顯增多，顯示了誘發(fā)地震活動(dòng)向深部的遷移[13]。

再如湖南省東江水庫(kù)，在蓄水初期，誘發(fā)地震的震源深度一般為2～2.5 km，而后期發(fā)生的地震一般都在5 km以上，少數(shù)地震已接近10 km。這充分說明庫(kù)水逐漸沿著庫(kù)區(qū)的斷裂構(gòu)造帶和溶洞向庫(kù)底深部滲透，使地下巖石發(fā)生物理和化學(xué)變化，從而誘發(fā)地震活動(dòng)[14]。應(yīng)該說，震源遷移的物理力學(xué)機(jī)制及判定方法研究對(duì)于水庫(kù)地震趨勢(shì)預(yù)測(cè)具有重要的意義，如果通過震源的遷移特征可以判斷誘發(fā)地震的震源分布叢集于一定深度，并且隨水位波動(dòng)不再向深部遷移，那么就有理由認(rèn)為后續(xù)發(fā)生中強(qiáng)以上誘發(fā)地震的概率減小，這方面蒙蒂塞洛水庫(kù)就是一個(gè)很好的例子。

水庫(kù)誘發(fā)地震震源分布的第三個(gè)特征是其空間分布受與水庫(kù)具有水力聯(lián)系的斷裂控制，此方面易立新[13]已有詳盡的論述，在此不贅述。

2 水庫(kù)蓄水過程與誘發(fā)地震響應(yīng)的關(guān)系

2.1 RIS時(shí)間響應(yīng)的基本特征與類型

RIS的直接起因與水庫(kù)蓄水過程有關(guān)。全球已報(bào)道的水庫(kù)誘發(fā)地震活動(dòng)觀測(cè)震例表明，誘發(fā)地震活動(dòng)在不同地區(qū)或不同水庫(kù)的時(shí)間響應(yīng)是有差別的。例如，有的誘發(fā)地震活動(dòng)發(fā)生在水庫(kù)蓄水后的數(shù)年才響應(yīng)，有的誘發(fā)地震活動(dòng)可能在蓄水過程中數(shù)天或數(shù)周內(nèi)就會(huì)立即發(fā)生，有的水庫(kù)誘發(fā)地震活動(dòng)與水位的周期性變化相關(guān)，兼有前兩者的特點(diǎn)，使得庫(kù)區(qū)的誘發(fā)地震處于長(zhǎng)期持續(xù)不斷的活動(dòng)過程中。

在全球99個(gè)MS≥3.0級(jí)的RIS震例中，有85例可以確定水庫(kù)初始蓄水時(shí)間與初次發(fā)生響應(yīng)性地震活動(dòng)的間隔時(shí)間。通過對(duì)這85個(gè)震例進(jìn)行統(tǒng)計(jì)分析表明：在初次蓄水1個(gè)月以內(nèi)就發(fā)生響應(yīng)性地震活動(dòng)的有22例，2～6個(gè)月發(fā)生的有19例，7～12個(gè)月發(fā)生的有11例，13～36個(gè)月發(fā)生的有16例，其余17例的發(fā)生時(shí)間間隔都在3年以上。初次發(fā)震時(shí)間間隔在1年以內(nèi)的占震例總數(shù)的61.18%，這說明在水庫(kù)蓄水后1年以內(nèi)，RIS活動(dòng)的初次響應(yīng)具有優(yōu)勢(shì)分布的特征（圖2）。有90個(gè)震例可以確定RIS主震與水庫(kù)初始蓄水的間隔時(shí)間，其中，間隔時(shí)間為1個(gè)月的有4例，2～6個(gè)月的有7例，7～12個(gè)月的有10例，13～36個(gè)月的有 26例，3年以上的有43例（圖3）。由發(fā)震頻次累計(jì)百分比圖可以看出，在水庫(kù)初始蓄水后1年以內(nèi)發(fā)生主震活動(dòng)的僅占震例總數(shù)的23.33%，而絕大部分誘發(fā)主震活動(dòng)都發(fā)生在水庫(kù)初始蓄水1年以后。

圖2 RIS初次響應(yīng)時(shí)間統(tǒng)計(jì)直方與累計(jì)百分比圖Fig.2 Statistical histogram and cumulative percentage figure of initial response time

圖3 RIS主震與水庫(kù)初始蓄水時(shí)間間隔統(tǒng)計(jì)直方與累計(jì)百分比圖Fig.3 Statistical histogram and cumulative percentage figure of the intervals between the largest reservoir-induced earthquake occurrance and the initial impounding reservoir

Simpson等根據(jù)水庫(kù)蓄水后誘發(fā)地震發(fā)生的時(shí)間響應(yīng)，把誘發(fā)地震分為快速響應(yīng)型和滯后發(fā)生型[15]。快速響應(yīng)型表現(xiàn)為，水庫(kù)初始蓄水，地震活動(dòng)立即增加，或水位迅速變化后地震活動(dòng)急劇變化，這類誘發(fā)地震活動(dòng)具有震級(jí)低、震源淺、多以震群發(fā)生的特點(diǎn)；而滯后發(fā)生型的主要地震活動(dòng)發(fā)生于水庫(kù)蓄水歷程的較晚階段，具有震級(jí)高、震源深、多以前震－主震－余震型發(fā)生的特點(diǎn)。從水庫(kù)誘發(fā)地震時(shí)間活動(dòng)現(xiàn)象分析，可將誘發(fā)地震活動(dòng)進(jìn)一步分為快速響應(yīng)型、滯后響應(yīng)型和混合響應(yīng)型3大類。圖4給出了這3種類型的典型例子：廣西龍灘水庫(kù)的水位變化與誘發(fā)地震活動(dòng)表現(xiàn)出與初始蓄水過程強(qiáng)烈的相關(guān)性，屬快速響應(yīng)型；阿斯旺水庫(kù)1964年開始蓄水，直到1981年底才誘發(fā)了大于3級(jí)的地震活動(dòng)，屬滯后響應(yīng)型[16]；美國(guó)米德湖水庫(kù)從1935年開始蓄水，之后一年多就不斷發(fā)生水庫(kù)誘發(fā)地震，至1939年5月發(fā)生5級(jí)的最大誘發(fā)地震，屬混合響應(yīng)型[17]。這三種類型僅是某些水庫(kù)誘發(fā)地震活動(dòng)表現(xiàn)出的特點(diǎn)，要嚴(yán)格區(qū)分它們的活動(dòng)規(guī)律仍然是困難的，總體來說，混合響應(yīng)型應(yīng)是水庫(kù)誘發(fā)地震活動(dòng)表現(xiàn)形式的主體，其它兩類誘發(fā)地震活動(dòng)僅僅是在一定的時(shí)間段內(nèi)反應(yīng)的階段性現(xiàn)象。

2.2 RIS響應(yīng)與水庫(kù)水位變化之間的相關(guān)性

2.2.1 RIS響應(yīng)與水庫(kù)水位變化的正相關(guān)性

圖4 典型水庫(kù)誘發(fā)地震活動(dòng)的時(shí)間響應(yīng)特征Fig.4 The time response characteristics of typical reservoir-induced seismicities

RIS對(duì)水庫(kù)水位變化響應(yīng)的一種形式是庫(kù)區(qū)地震活動(dòng)與水庫(kù)水位變化之間存在著正相關(guān)關(guān)系，即水位上升伴隨著誘發(fā)地震活動(dòng)的增強(qiáng)，而水位下降則對(duì)應(yīng)于地震活動(dòng)的減弱。例如：新豐江水庫(kù)蓄水后不久，庫(kù)區(qū)地震活動(dòng)便開始頻繁發(fā)生，且隨著水位上升，地震活動(dòng)逐漸增強(qiáng)。在1962年3月19日水庫(kù)水位達(dá)110.5 m時(shí)，發(fā)生了MS6.1級(jí)主震活動(dòng)，主震后至1964年期間，地震活動(dòng)頻度隨著水庫(kù)水位的下降而快速降低；巖灘水庫(kù)蓄水后庫(kù)區(qū)地震活動(dòng)與水庫(kù)水位升降的關(guān)系也非常密切，較密集的震群活動(dòng)或較大震級(jí)地震的發(fā)生均與水位到達(dá)峰值同步[18]；丹江口水庫(kù)蓄水初期水位上升比較緩慢，庫(kù)區(qū)沒有震級(jí)大于ML1.2級(jí)的地震發(fā)生，隨著水位逐漸上升，地震頻度和強(qiáng)度不斷升高[19]；南水水庫(kù)在蓄水初期的第一、二個(gè)高水位階段，地震活動(dòng)與庫(kù)水位的對(duì)應(yīng)關(guān)系較好，較高的水位或較大的水位變化，伴隨著較高水平的地震活動(dòng)[20]。在其他國(guó)家，如印度的柯依納水庫(kù)、希臘克里馬斯塔水庫(kù)、贊比亞－津巴布韋交界的卡里巴水庫(kù)、加拿大的馬尼克-3水庫(kù)和塔吉克斯坦的努列克水庫(kù)等，庫(kù)區(qū)地震活動(dòng)都表現(xiàn)出了與水庫(kù)水位變化的正相關(guān)對(duì)應(yīng)關(guān)系。

2.2.2 RIS響應(yīng)與水庫(kù)水位變化的負(fù)相關(guān)性

一些庫(kù)區(qū)的誘發(fā)地震響應(yīng)與水庫(kù)水位變化之間存在著負(fù)相關(guān)性，表現(xiàn)為在高水位下降后的低水位時(shí)段反而出現(xiàn)了誘發(fā)地震活動(dòng)增強(qiáng)的現(xiàn)象。珊溪水庫(kù)就是一個(gè)比較典型的例子，該水庫(kù)位于浙江省南部溫州市飛云江干流的上游河段，在大地構(gòu)造上處在浙東南褶皺帶內(nèi)的溫州－臨海拗陷帶的南部，為壓剪性的局部構(gòu)造應(yīng)力環(huán)境。庫(kù)區(qū)斷層構(gòu)造大多沿NW、NE兩組走向共軛發(fā)育，多為陡傾角的逆斷層、逆走滑斷層，一般長(zhǎng)達(dá)10 km以上，深達(dá)5 km以上，切穿基底變質(zhì)巖系（圖5a）。自從2000年5月水庫(kù)下閘蓄水以來，庫(kù)水位共有三個(gè)時(shí)段超過了132 m高水位，分別是2002年 5～9月、 2004年 9～12月和 2005年 5～11月。水庫(kù)誘發(fā)地震活動(dòng)主要發(fā)生在高水位過后，由高水位下降到低水位過程中，且地震活動(dòng)與高水位出現(xiàn)時(shí)間有一至數(shù)月的滯后時(shí)間，只有極少一些地震發(fā)生在水位上升時(shí)期（圖5b）。

呈負(fù)相關(guān)關(guān)系的震例還有青海的盛家峽水庫(kù)和臺(tái)灣的曾文水庫(kù)。盛家峽水庫(kù)自1980年汛期蓄水以來，當(dāng)1981年3～4月放水排空時(shí)，隨即誘發(fā)了地震，以后每年秋季蓄水滿庫(kù)時(shí)無震，而到翌年春天灌溉放水排空時(shí)，隨之誘發(fā)微震活動(dòng)，1984年3月7日發(fā)生了最大為MS3.6級(jí)的誘發(fā)地震。而臺(tái)灣的曾文水庫(kù)在滿庫(kù)時(shí)地震活動(dòng)減弱，庫(kù)水位下降后則地震活動(dòng)相對(duì)增強(qiáng)，但總的趨勢(shì)是蓄水后比蓄水前減弱。在1973年曾文水庫(kù)蓄水前，庫(kù)區(qū)發(fā)生過5～6級(jí)的中強(qiáng)地震，蓄水初期還頻繁發(fā)生過的3～4級(jí)地震，但當(dāng)水庫(kù)達(dá)到高水位后，震級(jí)頻次明顯降低。在1974～1976年水庫(kù)處于高水位期間沒有MS3級(jí)以上地震活動(dòng)，微震也極少；當(dāng)1977～1978年庫(kù)水位低于40 m后，震級(jí)隨之增大（MS3.7級(jí)），微震的頻次也隨之增高。

圖5 珊溪水庫(kù)斷裂構(gòu)造、震中分布(a)及地震序列與水位變化的關(guān)系(b)Fig.5 The fracture structure and epicenter distribution of Shanxi reservoir(a)and the relationship between earthquake sequence and the change of water level(b)

2.2.3 較高震級(jí)RIS響應(yīng)對(duì)水庫(kù)水位變化速率的依賴性

RIS與水庫(kù)水位變化之間的響應(yīng)關(guān)系還表現(xiàn)在較高震級(jí)的誘發(fā)地震活動(dòng)對(duì)水位變化速率的依賴性方面。Gupta在對(duì)印度柯依納水庫(kù)MS≥4.0級(jí)地震與水位變化速率之間關(guān)系的研究后發(fā)現(xiàn)，水位周變化率超過12 m/周是柯依納庫(kù)區(qū)發(fā)生MS≥5.0級(jí)地震的一個(gè)必要條件[21]。Simpson等也報(bào)道了在努列克水庫(kù)類似的發(fā)現(xiàn)[22]。在我國(guó)也不乏此類震例，如福建省的水口水庫(kù)。該水庫(kù)的水位變化與庫(kù)區(qū)地震活動(dòng)之間的對(duì)應(yīng)關(guān)系大體分為兩種情況：一是在蓄水初期水位尚未達(dá)到正常運(yùn)行水位的時(shí)段（1994年3月前），水位急劇升降隨后發(fā)震，如1993年7月23日3.3級(jí)和10月12日3.6級(jí)地震，這些地震發(fā)震時(shí)間滯后水位變化時(shí)間較短；二是水位基本上達(dá)到60 m以上的正常水位后，在高水位下急劇升降時(shí)誘發(fā)較大地震，如1995年3月16日3.2級(jí)和1996年4月21日4.1級(jí)地震，發(fā)震時(shí)間滯后于水位的變化較前一種情況稍長(zhǎng)。彭美鳳等對(duì)水口水庫(kù)l5次3級(jí)以上地震進(jìn)行統(tǒng)計(jì)，發(fā)現(xiàn)水位平均變化率達(dá)0.75～1.08 m/d，且連續(xù)變化幅度在6.5～10.5 m范圍內(nèi)發(fā)震。說明該水庫(kù)水的“沖擊波”式誘震作用的“靈敏性”要比柯依納水庫(kù)高[23]。因此，通過適當(dāng)控制水庫(kù)水位變化的速率，可能會(huì)避免未來局部破壞性誘發(fā)地震的發(fā)生。

2.2.4 RIS響應(yīng)與水庫(kù)水位變化之間關(guān)系的復(fù)雜性

震例研究表明，無論是正相關(guān)還是負(fù)相關(guān)，僅僅是特定水庫(kù)在特定的時(shí)間段內(nèi)反應(yīng)的階段性現(xiàn)象，對(duì)長(zhǎng)時(shí)間的蓄水歷程而言，RIS響應(yīng)與水庫(kù)水位變化的關(guān)系是十分復(fù)雜多樣的。

一些震例在水庫(kù)蓄水初期為正相關(guān)，隨著蓄水歷程的增長(zhǎng)，會(huì)逐漸轉(zhuǎn)變?yōu)樨?fù)相關(guān)關(guān)系。如新豐江水庫(kù)蓄水初期，水位變化與地震活動(dòng)之間存在著較好的正相關(guān)性，其中，MS6.1級(jí)主震和兩次MS5級(jí)以上強(qiáng)余震都發(fā)生在高水位階段。而1965年后這一規(guī)律已不明顯，ML≥3.0級(jí)以上地震多數(shù)發(fā)生在水位下降的過程中或由下降轉(zhuǎn)為平穩(wěn)波動(dòng)的過程中。這反映了蓄水?dāng)?shù)年后庫(kù)區(qū)積累的能量由于地震活動(dòng)的不斷釋放而達(dá)到新的平衡，水庫(kù)卸荷階段地震活動(dòng)增強(qiáng)，則可能是由于過?？紫秹毫υ斐伤矔r(shí)破裂引起的。

更多的震例表現(xiàn)為水庫(kù)蓄水初期正相關(guān)，隨著蓄水歷程的增長(zhǎng)，地震活動(dòng)呈逐步衰減的趨勢(shì)。如廣東省南水水庫(kù)，該水庫(kù)在蓄水初期（1969～1970年底）的第一、二個(gè)高水位階段，地震活動(dòng)與庫(kù)水位的對(duì)應(yīng)關(guān)系較好，較高的水位或較大的水位變化，往往伴隨較高水平的地震活動(dòng)。但是，自1976年第三個(gè)高水位階段之后，此對(duì)應(yīng)關(guān)系已基本消失。無論是高水位抑或較大的水位變化，均未能改變地震活動(dòng)逐步衰減的總趨勢(shì)。

還有一些震例，在水庫(kù)蓄水初期誘發(fā)地震活動(dòng)與水庫(kù)水位變化之間具有較好的對(duì)應(yīng)關(guān)系，但隨著水庫(kù)蓄水時(shí)間的延長(zhǎng)，庫(kù)區(qū)地震活動(dòng)似乎表現(xiàn)出了受控于區(qū)域應(yīng)力場(chǎng)擾動(dòng)的趨勢(shì)。例如水口水庫(kù)，上文已闡述，在該水庫(kù)蓄水初期，誘發(fā)地震活動(dòng)與水庫(kù)水位變化的對(duì)應(yīng)關(guān)系比較明顯；但在1999年以后，庫(kù)區(qū)地震活動(dòng)的起伏與周邊地震活動(dòng)的起伏表現(xiàn)出了一定的一致性，似乎說明此時(shí)段庫(kù)區(qū)地震活動(dòng)主要受控于區(qū)域應(yīng)力場(chǎng)的擾動(dòng)。

3 RIS時(shí)空分布與庫(kù)區(qū)地質(zhì)環(huán)境的關(guān)系

3.1 庫(kù)區(qū)地形地貌及地質(zhì)構(gòu)造對(duì)RIS空間分布的控制

統(tǒng)計(jì)分析表明，RIS多發(fā)生于深水峽谷區(qū)，而在寬谷地帶很少有誘發(fā)地震分布。這可能是因?yàn)閸{谷區(qū)岸坡是蓄水過程中地下水位變動(dòng)最劇烈的地帶，卸荷裂隙發(fā)育，有利于地應(yīng)力的集中，從而表現(xiàn)出了地形地貌對(duì)淺部誘發(fā)地震活動(dòng)的控制作用。

震例統(tǒng)計(jì)資料也表明：絕大部分構(gòu)造型中、強(qiáng)水庫(kù)地震的發(fā)生均與活動(dòng)斷裂或由于蓄水而活化了的斷裂活動(dòng)有關(guān)。因此，庫(kù)區(qū)存在與有水庫(kù)有水力聯(lián)系的斷裂構(gòu)造，特別是第四紀(jì)以來有過活動(dòng)的斷裂，是發(fā)生中、強(qiáng)水庫(kù)誘發(fā)地震的重要前提條件。RIS空間分布受與水庫(kù)有水力聯(lián)系斷裂構(gòu)造控制的典型例子是埃及的阿斯旺水庫(kù)，該水庫(kù)誘發(fā)地震活動(dòng)主要集中在兩個(gè)較小的區(qū)域內(nèi)；一個(gè)是馬拉瓦山附近，這一區(qū)域是水庫(kù)誘發(fā)地震最活躍、最集中的地區(qū)，MS5.3級(jí)主震就發(fā)生在這一地區(qū)；另一個(gè)區(qū)域是瓦迪卡拉沙附近。地質(zhì)證據(jù)表明，阿斯旺水庫(kù)誘發(fā)地震活動(dòng)的密集分布帶完好地沿卡拉沙斷層的三個(gè)斷層段展布，該斷層位于卡拉沙灣底，與水庫(kù)具有密切的水力聯(lián)系[12]。

研究還表明，RIS空間分布受局部褶皺構(gòu)造的控制。如塔吉克斯坦的努列克水庫(kù)，該水庫(kù)區(qū)中新生代地層由海相沉積層、大陸地臺(tái)沉積層和造山陸相磨拉石建造組成，巖鹽構(gòu)造發(fā)育，并且無震區(qū)和發(fā)震區(qū)的地層發(fā)育情況基本相同。但在局部構(gòu)造環(huán)境上，無震區(qū)為向斜構(gòu)造區(qū)，地下水很難向深部滲透；而發(fā)震區(qū)位于背斜構(gòu)造區(qū)，地下水可以向深部滲透引起孔隙壓升高從而誘發(fā)了地震活動(dòng)。

3.2 RIS空間分布對(duì)庫(kù)區(qū)巖性條件的依賴

RIS與庫(kù)區(qū)的巖性條件關(guān)系密切，早在水庫(kù)誘發(fā)地震研究的初期，Rothé等[24、25]就提出了預(yù)先存在的裂隙化、非均質(zhì)、脆性巖體是水庫(kù)誘發(fā)地震先決條件的觀點(diǎn)。對(duì)全球105個(gè)誘震最大震級(jí)在1.0級(jí)以上水庫(kù)庫(kù)壩區(qū)的巖性資料進(jìn)行了歸納和統(tǒng)計(jì)，結(jié)果顯示：在所有誘震最大震級(jí)ML≥1.0級(jí)的水庫(kù)中，庫(kù)壩區(qū)巖性為石灰?guī)r、花崗巖和片麻巖的占了76%之多，加上玄武巖、白云巖，接近占到90%；而發(fā)震庫(kù)區(qū)巖性為砂巖、泥巖、千枚巖和大理巖的僅僅占到10.5%。其中，變質(zhì)巖巖性以片麻巖為主，還有少量的千枚巖；沉積巖巖性幾乎全是石灰?guī)r和白云巖，只有非常少的砂巖和泥巖；巖漿巖巖性以花崗巖占絕大多數(shù)，還有部分玄武巖（圖 6）。

之所以會(huì)出現(xiàn)這種現(xiàn)象，主要是由巖體力學(xué)性質(zhì)和滲透特性兩方面所決定的。首先就巖體力學(xué)性質(zhì)而言，由石灰?guī)r、花崗巖、片麻巖和玄武巖等組成的脆性巖體，較延性巖體更容易聚集變形能量，原位地應(yīng)力高于延性巖體，在外力作用下易引起脆性斷裂。其次就巖體滲透特性而言，脆性巖體往往構(gòu)造節(jié)理、裂隙發(fā)育，滲透性能較好且表現(xiàn)出非均一性和各向異性，特別是碳酸鹽體巖溶的發(fā)育為水的滲透循環(huán)創(chuàng)造了更為有利的條件，這樣在水庫(kù)附加水頭壓力的擴(kuò)散作用下，容易引起孔隙壓力在特定方向和特定部位的增高，從而誘發(fā)地震活動(dòng)。因此，RIS對(duì)庫(kù)區(qū)巖性條件的依賴也同時(shí)決定了誘發(fā)地震活動(dòng)空間分布的規(guī)律?？v觀水庫(kù)誘發(fā)地震空間分布特征，震中往往集中在庫(kù)體重力場(chǎng)和地表水體向深部滲流影響范圍內(nèi)的脆性巖體分布帶，表現(xiàn)出了RIS對(duì)庫(kù)區(qū)巖性條件強(qiáng)烈的依賴性。

圖6 水庫(kù)庫(kù)壩區(qū)巖性條件統(tǒng)計(jì)圖Fig.6 The lithology statistical chart of reservoir dam area

3.3 庫(kù)區(qū)水文地質(zhì)條件對(duì)RIS時(shí)間響應(yīng)的控制

水庫(kù)蓄水后，庫(kù)體重力荷載作用和水庫(kù)附加水頭壓力向深部的擴(kuò)散作用是誘發(fā)庫(kù)底巖體或斷層破裂失穩(wěn)的兩個(gè)重要的外在因素，RIS的時(shí)間響應(yīng)和演化明顯受控于這兩個(gè)作用相互制約與促進(jìn)的動(dòng)態(tài)耦合過程[9、10]，而地表水體向深部滲透的程度和范圍主要取決于庫(kù)區(qū)水文地質(zhì)條件，從而決定了水文地質(zhì)條件對(duì)RIS時(shí)間序列的控制作用。通過對(duì)已發(fā)震水庫(kù)庫(kù)區(qū)水文地質(zhì)條件的總結(jié)不難發(fā)現(xiàn)：當(dāng)庫(kù)區(qū)巖體的含水性和導(dǎo)水性差，特別是當(dāng)庫(kù)區(qū)巖體為裂隙不發(fā)育的脆性花崗巖、玄武巖體時(shí)，一般在蓄水以后不久就會(huì)發(fā)生快速響應(yīng)型誘發(fā)地震活動(dòng)。例如，蒙蒂塞洛水庫(kù)分布于片巖、片麻巖并侵入有花崗巖、花崗閃長(zhǎng)巖的地區(qū)。該區(qū)地層普遍含水性差，侵入體裂隙不發(fā)育，為隔水巖體。水庫(kù)蓄水后誘發(fā)地震活動(dòng)的快速響應(yīng)可能反映了庫(kù)水荷載效應(yīng)起到了主導(dǎo)作用。此外，馬尼克-3、卡里巴、克里馬斯塔等水庫(kù)也屬于這一情況。

當(dāng)庫(kù)區(qū)巖體的含水性和導(dǎo)水性好，特別是發(fā)育有與水庫(kù)有密切水力聯(lián)系的導(dǎo)水性斷裂構(gòu)造時(shí)，一般會(huì)發(fā)生混合響應(yīng)型和滯后響應(yīng)型的誘發(fā)地震活動(dòng)。如美國(guó)的米德湖水庫(kù)，該水庫(kù)蓄水后，地震主要發(fā)生在水庫(kù)西南部的Boulder盆地。從水文地質(zhì)條件分析，Boulder盆地第三紀(jì)火成巖系和寒武紀(jì)巖漿巖系中發(fā)育有沉積巖系，力學(xué)性質(zhì)上是兩層脆性巖體夾一層塑性巖體，含水性上是兩層隔水層夾一層含水層，并有導(dǎo)水性斷層發(fā)育。從地震序列分析，該水庫(kù)誘發(fā)的地震活動(dòng)屬于混合型響應(yīng)類型，蓄水初期的地震可能與庫(kù)體荷載作用有關(guān)，而后期發(fā)生的地震可能是庫(kù)水?dāng)U散至沉積巖系與寒武系巖漿巖系接觸面上引起的孔隙壓力升高造成的，這種解釋與該區(qū)的震源深度也基本吻合[13]。再如我國(guó)的新豐江水庫(kù)，幾個(gè)震中密集區(qū)都位于滲透性好的原生節(jié)理和構(gòu)造節(jié)理極端發(fā)育的花崗巖體內(nèi)，而在滲透性差的紅色碎屑巖建造、復(fù)理石建造中，地震極少或沒有，誘發(fā)地震活動(dòng)時(shí)間序列表現(xiàn)出初期快速響應(yīng)及后期滯后響應(yīng)的混合響應(yīng)的特點(diǎn)。

事實(shí)上，庫(kù)區(qū)局部地質(zhì)構(gòu)造、巖性條件和水文地質(zhì)條件對(duì)RIS時(shí)空分布的控制作用并不是截然分開的，它們之間存在著相互依存的、辯證統(tǒng)一的關(guān)系。易立新等[26]將斷裂構(gòu)造、巖性條件和水文地質(zhì)條件概括為水文地質(zhì)結(jié)構(gòu)，既包含了斷層導(dǎo)水性能的好壞，又涵蓋了巖性條件的差異，同時(shí)指出了斷層與圍巖不同組合形式對(duì)RIS的影響，這一歸納就比較恰當(dāng)。

4 結(jié)論

本文在廣泛收集國(guó)內(nèi)外水庫(kù)震例資料的基礎(chǔ)上，利用統(tǒng)計(jì)學(xué)的方法開展了RIS空間分布與時(shí)間響應(yīng)特征及其與庫(kù)區(qū)地質(zhì)構(gòu)造環(huán)境關(guān)系的研究，得到以下幾點(diǎn)主要認(rèn)識(shí)：

（1）RIS在空間分布上主要集中在水庫(kù)周緣10 km范圍之內(nèi)，并局限于特定的庫(kù)段叢集分布。當(dāng)庫(kù)區(qū)具備一定的地質(zhì)構(gòu)造與水文地質(zhì)條件時(shí)，誘發(fā)地震活動(dòng)會(huì)隨著水庫(kù)蓄、放水過程中的水位變化而發(fā)生遷移。RIS遷移現(xiàn)象可概括為4種情況：①是在水庫(kù)蓄水初期，庫(kù)區(qū)外圍的地震活動(dòng)向庫(kù)壩區(qū)遷移；②是隨著水庫(kù)淹沒區(qū)范圍的擴(kuò)大，地震向遠(yuǎn)離庫(kù)壩的灰?guī)r峽谷區(qū)遷移；③是地震沿與水庫(kù)有水力聯(lián)系的斷裂構(gòu)造向遠(yuǎn)離庫(kù)岸的方向遷移；④是在處于擠壓構(gòu)造應(yīng)力環(huán)境的庫(kù)區(qū)，受水庫(kù)水位動(dòng)態(tài)變化的影響，地震活動(dòng)發(fā)生動(dòng)態(tài)遷移。此外，與天然構(gòu)造地震相比，水庫(kù)誘發(fā)地震具有震源深度淺的特征，一般分布在地下10 km以上的范圍，但對(duì)于不同的水庫(kù)，其誘發(fā)地震震源深度亦存在著較大的差異。

（2）震例統(tǒng)計(jì)分析表明，水庫(kù)蓄水后1年以內(nèi)，RIS活動(dòng)的初次響應(yīng)具有優(yōu)勢(shì)分布的特征，而絕大部分誘發(fā)主震活動(dòng)都發(fā)生在水庫(kù)初始蓄水1年以后。誘發(fā)地震活動(dòng)在不同地區(qū)或不同水庫(kù)的時(shí)間響應(yīng)存在著差異，從RIS時(shí)間響應(yīng)現(xiàn)象分析，可分為快速響應(yīng)、滯后響應(yīng)和混合響應(yīng)3種類型?？傮w來說，混合響應(yīng)型應(yīng)是水庫(kù)誘發(fā)地震活動(dòng)表現(xiàn)形式的主體，其它兩類誘發(fā)地震活動(dòng)僅僅是在一定的時(shí)間段內(nèi)反應(yīng)的階段性現(xiàn)象。

（3）RIS時(shí)間響應(yīng)與水庫(kù)水位變化之間關(guān)系密切，一些庫(kù)區(qū)地震活動(dòng)與水庫(kù)水位變化之間存在著正相關(guān)關(guān)系，而一些庫(kù)區(qū)則存在著負(fù)相關(guān)性。但無論是正相關(guān)還是負(fù)相關(guān)，僅僅是特定水庫(kù)在特定的時(shí)間段內(nèi)反應(yīng)的階段性現(xiàn)象，對(duì)長(zhǎng)時(shí)間的蓄水歷程而言，RIS響應(yīng)與水庫(kù)水位變化的關(guān)系是可以發(fā)生變化的。震例研究同時(shí)表明，較高震級(jí)的誘發(fā)地震活動(dòng)依賴于水位變化的速率。

（4）RIS時(shí)空分布受發(fā)震庫(kù)區(qū)三個(gè)方面的先存內(nèi)因的控制：①是地形地貌與局部地質(zhì)構(gòu)造條件。在地形地貌上，一般深水峽谷區(qū)易發(fā)生卸荷型淺層誘發(fā)地震活動(dòng)；在地質(zhì)構(gòu)造條件上，絕大部分構(gòu)造型中、強(qiáng)水庫(kù)地震的發(fā)生均與活動(dòng)斷裂或由于蓄水而活化了的斷裂活動(dòng)有關(guān)，同時(shí)，局部褶皺構(gòu)造也表現(xiàn)出了對(duì)RIS空間分布的控制作用。②是庫(kù)區(qū)巖體的巖性組成，統(tǒng)計(jì)結(jié)果顯示，誘發(fā)地震往往發(fā)生在脆性巖體分布帶，表現(xiàn)出了RIS對(duì)庫(kù)區(qū)巖性條件的依賴性，這可能主要是由巖體力學(xué)性質(zhì)和滲透特性兩方面所決定的。③是庫(kù)區(qū)水文地質(zhì)條件，由于地表水體向深部滲透的程度和范圍主要取決于庫(kù)區(qū)水文地質(zhì)條件，從而決定了水文地質(zhì)條件對(duì)RIS時(shí)間序列的控制作用。事實(shí)上，庫(kù)區(qū)局部地質(zhì)構(gòu)造、巖性條件和水文地質(zhì)條件對(duì)RIS時(shí)空分布的控制作用并不是截然分開的，它們之間存在著相互依存的、辯證統(tǒng)一的關(guān)系。

[1]秦嘉政，劉麗芳，錢曉東.水庫(kù)誘發(fā)地震活動(dòng)特征及其預(yù)測(cè)方法研究[J].地震研究，2009，32(2)：105-113.

[2]Nikolaev N I.Tectonic conditions favourable for causing earthquakes occurring in connection with reservoir filling[J].Eng.Geol.， 1974， 8： 171.

[3]Gupta H K,Rastogi B K,Narain H.Common features of the reservoir associated seismic activities[J].Bull.Seismol.Soc.Am.， 1972， 62： 481～492.

[4]Castle R O,Clark MM,Grantz A，et al.Tectonic state,its significanceand characterization in theassessmentof seismic effects associated with reservoir impounding[J].Eng.Geol.， 1980， 15(1-2)： 53～99.

[5]Guha S K,Patil D N.Large Water-Reservoir-related Induced Seismicity[J].Gerland Beitrage zur Geophysik，1990， 99： 265.

[6]Pei D Z.Characteristics of small earthquakes and stress field:a study of reservoir induced sSeismicity in quebec[J].Soil Dynamics&Earthquake Eng， 1992， 11(4)： 193～202.

[7]王清云，高士鈞.丹江口水庫(kù)誘發(fā)地震研究[J].大地測(cè)量與地球動(dòng)力學(xué)，2003，23(1)：103-106

[8]Simpson D W,Negmatullaev S K.Induced seismicity at Nurek reservoir[J].Bull.Seismol.Soc.Am.， 1981， 71：1561.

[9]周斌，薛世峰，鄧志輝，等.水庫(kù)誘發(fā)地震時(shí)空演化與庫(kù)水加卸載及滲透過程的關(guān)系-以紫坪鋪水庫(kù)為例[J].地球物理學(xué)報(bào)，2010，53(11)：2651～2670.

[10]周斌.水庫(kù)誘發(fā)地震時(shí)空演化特征及其動(dòng)態(tài)響應(yīng)機(jī)制研究-以紫坪鋪水庫(kù)為例 [J].國(guó)際地震動(dòng)態(tài)，2010，6：41-44.

[11]Talwani P,Acree S.Pore pressure diffusion and the mechanism of reservoir～I(xiàn)nduced seismicity[J].Pure and Applied Geophysics.， 1985， 122(6)： 947～965.

[12]Mohamed A,Megume M.Earthquake activity in the Aswan region,Egypt[J].Pure and Applied Geophysics，1995， 145(1)： 69-86.

[13]易立新.水庫(kù)誘發(fā)地震理論及三峽水庫(kù)誘發(fā)地震預(yù)測(cè)研究[D].北京：中國(guó)地震局地質(zhì)研究所，2001.

[14]曹華.東江水庫(kù)誘發(fā)地震探討[J].大壩與安全，2006， 1： 21-22.

[15]Simpson D W,Leith W S,Scholz C H.Two types of reservoir-induced seismicity[J].Bull.Seismol.of America， 1988， 78(6)： 2025～2040

[16]Mekkawi M,Grasso J R,Schnegg P A.A long-lasting relaxation of seismicity atAswan Reservoir,Egypt，1982-2001[J].Bull.Seismo.Soc.Am.,2004， 94(2)：479～492.

[17]Carder D S.Seismic investigations in the Boulder Dam area,1940-1944,and the influence of reservoir loading on earthquake activity [J].Bull.Seismol.Soc.Am.，1945， 35： 175～192

[18]光耀華.大化巖灘梯級(jí)水庫(kù)誘發(fā)地震特征[J].水力發(fā)電學(xué)報(bào)， 1996， 4： 45-53.

[19]劉素梅.丹江口水庫(kù)續(xù)建工程誘發(fā)地震研究 [D].武漢：武漢理工大學(xué)，2005.

[20]肖安予.南水水庫(kù)地震及其發(fā)展趨勢(shì)[J].華南地震，1990， 10(2)： 68～77.

[21]Gupta H K.Induced seismicity hazard mitigation through water level manipulation at koyna,India:a suggestion[J].Bull.Seismol.Soc.Am.,1983,73(2):679～682.

[22]Simpson D W,Leith W S.Induced seismicity at Toktogul Reservoir,Soviet entral Asian.U.S[J].Geo.Surv.,1988,No.14-08-0001-G1188,32.

[23]彭美鳳，林世敏，林松鍵.水口水庫(kù)地震及其活動(dòng)特征[J].華南地震， 1997， 17(2)： 83～89.

[24]Rothé J P.Seismic artificiels (man-made earthquakes)[J].Tectonophysics,1970,9:215～238.

[25]Carder D S.Reservoir loading and local earthquakes[J].Geol.Soc.Am.Bull.,1970,81(8):51～61

[26]易立新，王廣才，李榴芬.水文地質(zhì)結(jié)構(gòu)與水庫(kù)誘發(fā)地震[J].水文地質(zhì)工程地質(zhì)， 2004， 2： 29～32.