重慶三江口水庫誘發地震分析研究①

吳建超， 張麗芬， 汪旭濤， 蔡永建

(1.中國地震局地震研究所(地震大地測量重點實驗室)，湖北 武漢 430071； 2.武漢地震工程研究院，湖北 武漢 430071)

重慶三江口水庫誘發地震分析研究①

吳建超1，2， 張麗芬1， 汪旭濤1，2， 蔡永建2

(1.中國地震局地震研究所(地震大地測量重點實驗室)，湖北 武漢 430071； 2.武漢地震工程研究院，湖北 武漢 430071)

通過野外地震地質調查并參考前人研究成果，分析三江口水庫區的地質構造背景、地震活動性及水文地質條件等資料，對該水庫誘發地震的可能性進行分析。構造類比法分析表明：蓄水后發生構造型水庫誘發地震的可能性較小，但有可能發生巖溶塌陷型水庫誘發地震。概率預測法分析表明：庫首段(新灘子—獅貍彎)發震概率較小，僅為0.02；庫中段(獅貍彎—牛鼻子)巖溶不太發育，誘震可能性較小，不發震的概率為0.96；庫尾段(牛鼻子—峽馬口)有可能誘發微震，發震(M<3.0)概率為0.10。

三江口水庫； 構造類比法； 概率預測法； 水庫誘發地震； 發震概率

0 引言

三江口水庫壩址位于重慶市彭水縣普子河中下游段的新灘子附近，水庫正常蓄水位306 m，庫容0.68×108m3。大壩為混凝土重力壩，壩高71.5 m，壩頂長218.4 m，裝機容量18.9 MW，屬III型中等峽谷型水庫。

水庫蓄水導致庫區地震活動增多，這已被國內外大量研究所證實[1-4]。水庫誘發地震會對工程設施和居民生產生活產生不同程度的破壞和影響。三江口水庫于2008年開始設計論證，水庫蓄水后能否誘發地震、地震對庫區環境的影響等與水庫所處的地震構造環境有密切關系。本文在分析三江口水庫區地質構造背景、地震活動性及水文地質條件等資料的基礎上，對庫區誘發地震的可能性進行分析計算，為進一步開展水庫誘發地震的監測研究和防震減災工作奠定基礎。

1 庫區地質背景

1.1 地質構造

三江口水庫在大地構造上位于揚子準地臺上揚子臺褶帶黔江拱褶斷束內，其北西側以巫山—金佛山基底斷裂為界，與四川臺坳川東褶皺束相接[5]。區域構造以寬緩背斜、向斜或箱狀復式背斜及伴生斷裂為特征。水庫壩址位于普子復向斜東南翼，庫區及外圍NE-NNE向褶皺和斷裂呈有規律的分布，構成了庫區內的基本構造格局。

庫區出露的斷裂主要為東南緣的郁山斷裂，其次為西緣的馬武斷裂(圖1、表1)。其中，郁山斷裂發育在郁山背斜西北翼下奧陶統南津關組中，地貌上多呈現洼地、斷層崖、埡口。斷層帶由斷層角礫巖、碎裂巖組成，膠結較緊密，斷層面不平整，略粗糙，未見有擦痕。斷裂活動具有垂直分量，為正斷層。斷裂南北兩段的測年數據(TL法和ESR法)集中在30～50×104a，表明該斷裂最新活動時代為中更新世中期，晚更新世以來沒有活動的跡象。馬武斷裂為前第四紀斷裂。

表1 三江口水庫區主要斷裂概況

1.2 地層巖性

水庫區出露的地層比較簡單，主要為奧陶系、二疊系和三疊系地層。第四系零星發育在沖溝溝口、支流與干流匯合處、高陡岸坡腳和河流階地上。奧陶系主要分布于普子復向斜西北翼和郁江背斜的核部。巖性以灰黑色薄至中層灰巖、泥質灰巖、瘤狀灰巖為主，局部含燧石結核，并夾有少量褐色厚層結晶灰巖及灰綠、黃綠色頁巖。二疊系地層分布于普子復向斜的兩翼，巖性以深灰色生物灰巖和泥灰巖為主，頂部為含燧石條帶灰巖等。三疊系地層構成了普子復向斜的核部，巖性主要為泥質灰巖、白云質灰巖和灰巖等。

1.3 地形地貌

三江口水庫區以中低山和小型山間盆地地貌為主，山脈連綿，溝谷發育。庫區兩岸邊坡高陡，地勢北東高，南西低，呈階梯狀降低。普子河蜿蜒曲折，“V”型和“U”型谷相間分布。由于受地勢的影響，河流由北東向南西依勢而流，地貌類型多受巖性的控制。在碎屑巖分布區以侵蝕、剝蝕地貌為主，多呈脊狀低中山地貌，溪溝縱橫交錯，相對高差500～1 000 m，平均自然坡度在40°～60°。而在碳酸鹽巖分布地區，地形切割較強烈，多以溶蝕地貌為主，如在二疊系和三疊系等灰巖分布地區，地貌形態多以峰林谷地、峰林洼地和峰叢槽谷為主，在垂直方向上可見有明顯的臺階或夷平面。

2 庫區地震活動背景

三江口水庫區歷史和現今地震活動水平較低，近幾十年來庫區20 km范圍內有26次小震活動(圖1)。其中，M1.0～1.9地震21次，M2.0～2.9地震4次，M3.0以上地震僅有1次，發生于1979年12月31日，距三江口壩址約2 km。這些小震、微震活動大部分分布于郁山斷裂的上盤，呈現出NE-SW向的條帶狀展布特征，與本地區主要地質構造方向基本一致。

3 水庫誘發地震危險性分析

3.1 水文地質條件分析

水庫誘發地震是復雜的科學問題，目前還處在探索過程中。根據前人的研究成果[6]，水庫誘發地震的類型可分為構造型、巖溶塌陷型和應力局部調整型三類，各類水庫誘發地震又各具不同的發震條件。構造型水庫誘發地震多與斷裂和庫水的聯合作用相關；巖溶型水庫誘發地震多發生在大面積碳酸鹽巖出露、現代巖溶作用強烈、具備向深部滲流條件的地區；應力局部調整型水庫誘發地震的發震條件為河谷岸坡相對高差大、巖體微裂隙較發育、有利于庫水滲透的水文地質結構面。從已知地震經驗類比總結出來的發震條件是分析判斷水庫是否會誘發地震的基礎[7-8]。

庫區巖性條件對庫水向地下滲透起著重要作用。在我國已發生的27例水庫地震中，發生在碳酸鹽巖地層中的有18例，占66%；發生在花崗巖類巖石中的有5例，約占18%[9]。巖性條件是產生水庫誘發地震的直接誘因，同時也是討論水庫誘發地震的重要依據。

三江口水庫區地下水類型以裂隙溶洞水為主，孔隙潛水和基巖裂隙水次之。由于溶洞、暗河、泉水、巖溶管道系統發育，裂隙溶洞水常形成地下水匯集和排泄通道，多以層狀出現。由于水庫區的中段和庫尾段巖石性質和巖石類型略有不同，故巖溶發育程度不一。如圖2所示，牛鼻子-峽馬口一帶地層為三疊系灰巖，透水性較強，溶洞、暗河和泉水等均較發育，是地表水與地下水動力聯系的重要通道，也是庫水向深部滲透的良好條件，有利于產生誘發地震。而在牛鼻子-獅貍彎庫段，主要為二疊系地層，巖溶發育略弱，多以峰叢、槽谷、溶溝等為主，地表水與地下水聯系不如前者，因此產生誘發地震的機率較小。

圖2 三江口水庫區碳酸鹽巖透水層與碎屑巖隔水層分布圖Fig.2 Distribution of the carbonate permeable layer and clastic rock impermeable layer in Sanjiangkou reservoir area

綜上所述，三江口庫區不存在活動斷裂，幾乎不存在發生構造型水庫地震的可能。但蓄水后誘發巖溶塌陷型水庫地震的可能性較大，主要集中在二疊系和三疊系碳酸鹽巖大面積分布、巖溶管道發育的庫段。

3.2 構造類比分析

該方法根據對已發震水庫資料的分析，總結歸納出誘發地震的主要條件。這些條件在擬建水庫中越齊備越典型，該水庫誘發地震的可能性就越大。

在與水庫規模、地質條件等相近的四川銅街子、云南魯布革、湖北前進、湖南黃石等已發震水庫進行對比后，認為三江口水庫蓄水后發生構造型水庫誘發地震的可能性較小。其理由是：庫區沒有較大的活動斷裂通過，僅在庫區東南緣和西北緣分布有兩條非活動斷裂，即郁山斷裂和馬武斷裂。庫區位于普子復向斜內，沒有展布大型的地球物理異常帶。本文認為水庫蓄水后有可能發生巖溶塌陷型水庫地震。其理由是：庫區內多為中厚層灰巖出露，巖溶較發育，庫區內的褶皺軸與水庫區展布基本相似，容易形成地下水管流系統。蓄水后庫水的荷載作用和滲透作用容易引起溶洞內巖體塌陷和裂隙擴展，進而誘發巖溶塌陷型水庫地震。

3.3 概率預測分析

該方法是Packer和Beacher提出來的，是在概率意義上對一個新建水庫的誘發地震危險性進行預測，從現有水庫與誘發地震有密切關系的因素進行統計分析[11-12]。

水庫誘發地震受多種因素的影響，一般把水庫誘發地震過程看成一個多因子函數，將影響因子看作隨機因子，真實全面地獲得水庫誘發地震的組合環境資料，通過概率統計等數學模型來計算和處理，就可以得到條件和結果之間的函數關系[13-14]。

水庫誘發地震的影響因素很多，不同的研究者所選取的誘震因子集也往往不同。考慮到三江口庫區的具體條件，依據蘇錦星等[14]提出的誘震因子，選用庫水深度(D)、庫區巖性條件(G)、區域應力狀態(S)、斷層活動性(F)、巖溶發育程度(K)、斷層透水深度(H)、與庫水溝通關系(C)和地震活動背景(E)等八個誘震因子進行震級概率預測計算 (表2)。八因子方案將預測水庫地震震級分為“強、中強、弱、微、無震”5種狀態。具體為：M4代表強震(M≥6.0)，M3代表中強震(6.0>M≥4.5)，M2代表弱震(4.5>M≥3.0)，M1代表微震(M<3.0)，M0代表無震(不發震)。

表2 水庫誘發地震影響因素及其狀態[13]

根據貝葉斯條件概率理論，預測水庫誘發地震的統計模型可表達為：

(1)

式中：P(Mi/D，G，S，F，K，C，H，E)為不同震級地震類別的先驗概率。根據國內外251座大型水庫資料統計的發震概率，發生強震、中強震、弱震、微震和不發震的先驗概率分別為：P(M4)=0.02；P(M3)=0.04；P(M2)=0.05；P(M1)=0.07；P(M0)=0.82。P(D，G，S，F，K，C，H，E/Mi)為各影響因素組合條件下不同震級的條件概率，即P(D，G，S，F，K，C，H，E/Mi)=P(D/Mi)·P(G/Mi)·P(S/Mi)·P(F/Mi)·P(K/Mi)·P(H/Mi)·P(C/Mi)·P(E/Mi)。表3為八項誘震因子在不同震級下的條件概率統計結果。

表3 誘震因子不同狀態發震先驗概率統計表[13]

根據三江口庫區各庫段誘發地震的地質構造條件和誘發地震環境，本文將庫區自下游至上游劃分為三段：庫首段(新灘子-獅貍彎)、庫中段(獅貍彎—牛鼻子)和庫尾段(牛鼻子—峽馬口)。利用MATLAB軟件編程，計算得出三個庫段分屬五個發震震級的概率，取其中概率值最大或較大者對應的地震震級分級區間，即為預測的可能發震強度，見表4。

由計算結果可知，庫首段(新灘子-獅貍彎)為水庫誘發地震發震概率較小的庫段，發震概率僅為0.02，說明庫壩區較穩定。庫中段(獅貍彎-牛鼻子)雖然地表出露二疊系灰巖，但巖溶不太發育，不發震的概率為0.96，為誘震可能性較小的地段。庫尾段(牛鼻子-峽馬口)雖不發震的概率約0.83，但誘震M1(M<3.0)的發震概率為0.10，說明庫尾段有可能誘發微震。

4 討論與結論

本文采用構造類比和概率預測兩種方法對三江口水庫誘發地震問題進行分析研究，得出的預測結果基本一致。

庫首段(新灘子-獅貍彎)的河谷地貌形態以中低山、丘陵的曲流寬谷與狹谷相間排布為特征，岸坡大多較緩，組成河谷的地層為隔水層，滲透條件不暢，且無較大斷裂通過，缺乏地震活動。確定性評價認為不易發生水庫誘發地震，計算不發震的概率為0.98。

庫中段(獅貍彎-牛鼻子)的基本條件與庫首區相同，但庫區巖性為弱透水層，局部地段巖溶發育，有利于庫水滲透。確定性評價認為發生水庫誘發地震的可能性較小，計算不發震的概率為0.96。

庫尾段(牛鼻子—峽馬口)最顯著的特征是：河谷呈狹窄的V型，局部為V型障谷、隘谷，岸坡陡峻，幾近直立，三疊系灰巖中巖溶發育，溶洞、暗河地下水量較大，易于滲透。確定性評價認為。該段發生巖溶塌陷型水庫地震的可能性較大。概率計算得到的誘發M1(M<3.0)的發震概率為0.10。綜合評價庫尾段(牛鼻子—峽馬口)水庫誘發地震的可能性較大，最大震級小于3.0級。

表4 誘震因子狀態及統計模型預測結果

水庫地震誘發地震危險性評估方法有很多，如類比法、成因模式法、概率預測法、模糊數學法、灰色系統模型法、灰色聚類方法等，它們是預測水庫誘發地震危險性的主要方法[15-16]。今后應進一步采用多種誘震預測評價方法計算成果和各種誘震因素分析，判斷三江口水庫在未來關鍵時段內的最大誘震震級。

References)

[1] Gupta H K.A Review of Recent Studies of Triggered Earthquakes by Artificial Water Reservoirs with Special Emphasis on Earthquakes in Koyna, India[J].Earth-Science Reviews，2002(3):279-310.

[2] 賀為民，秦建增，劉明軍,等.小浪底水庫誘發地震預測研究[J].西北地震學報，2001，23(2)：164-168.HE Wei-min，QIN Jian-zeng，LIU Ming-jun，et al.Forecast on Induced Earthquake for Xiaolangdi Reservoir[J].Northwestern Seismological Journal，2001，23(2)：164-168.(in Chinese)

[3] 馬文濤，藺永，范京立，等.水庫誘發地震的震例比較與分析[J].地震地質，2013，35(12)：914-921.MA Wen-tao，LIN Yong，FAN Jing-li，et al.Comparison and Analysis on the Basic Features of Reservoir-induced Seismicity[J].Seismology and Geology，2013，35(12)：914-921.(in Chinese)

[4] 陳蜀俊.三峽庫首區蓄水前后構造應力場數值模擬及地震危險性研究[D].武漢：中國地質大學，2005：23-37.CHEN Shu-jun.Numerical Modeling of Tectonic Stress Field Before and After Reservoir Impounding in the Head Area of the Three Gorges and Studies on the Seismic Risk[D].Wuhan：China University of Geoscience，2005：23-37.(in Chinese)

[5] 潘桂棠，肖慶輝，陸松年，等.中國大地構造單元劃分[J].中國地質，2009，36(1)：6-13.PAN Gui-tang，XIAO Qing-hui，LU Song-nian，et al.Subdivision of Tectonic Units in China[J].Geology in China，2009，36(1)：6-13.(in Chinese)

[6] 王秋良，姚運生，夏金梧，等.統計預測模型對三峽庫區誘發地震預測的應用[J].地震地質，2009，31(2)：287-294.WANG Qiu-liang，YAO Yun-sheng，XIA Jin-wu，et al.Application of Statistical Forecasting Model to the Prediction of the Three Gorges Reservoir Induced Seismicity[J].Seismology and Geology，2009，31(2)：287-294.(in Chinese)

[7] 賀為民，楊杰.中波水電站水庫誘發地震預測研究[J].自然災害學報，2013，22(5)：105-108.HE Wei-min，YANG Jie.Research on Forecast of Earthquake Induced by Reservoir of Zhongbo Hydropower Station[J].Journal of Natural Disasters，2013，22(5)：105-108.(in Chinese)

[8] 蘇錦星，汪雍熙.水庫誘發地震可能性的定量評價[J].工程地質學報，2000，8(增刊)：70-73.SU Jin-xing，WANG Yong-xi.Quantitative Evaluation of Reservoir Induced Earthquake Possibility[J].Journal of Engineering Geology，2000，8(Supp)：70-73.(in Chinese)

[9] 楊清源，胡毓良，陳獻程，等.國內外水庫誘發地震目錄[J].地震地質，1996，18(4)：453-461.YANG Qing-yuan，HU Yu-liang，CHEN Xian-cheng，et al.Reservoir Induced Earthquake Catalogue at Home and Abroad.[J].Seismology and Geology，1996，18(4)：453-461.(in Chinese)

[10] Beacher G B.Reservoir Induced Seismicity，Ⅱ：Probability of Risk at Deep And/or Very Large Reservoirs[J].G S A，1979(7)：3-9.

[11] Packer D R.Reservoir Induced Seismicity，Ⅱ：Summary and Implications for Occurrence of Risk[J].G S A，1979(7)：12-18.

[12] 吳建超，張麗芬，蔡永建，等.重慶羊頭鋪水庫誘發地震分析研究[J].地震研究，2012，35(3)：338-340.WU Jian-chao，ZHANG Li-fen，CAI Yong-jian，et al.Research on Earthquake Induced by Yangtoupu Reservoir in Chongqing[J].Journal of Seismological Research，2012，35(3)：338-340.(in Chinese)

[13] 蘇錦星，夏其發，汪雍熙，等.大藤峽水庫誘發地震初步預測[J].中國地質災害防治學報，1993，4(4)：48-50.SU Jin-xing，XIA Qi-fa，WANG Yong-xi，et al.Preliminary Prediction on Datengxia Reservoir Induced Seismicity[J].The Chinese Journal of Geological Hazard and Control，1993，4(4)：48-50.(in Chinese)

[14] 劉先鋒，李碧雄，鄧建輝，等.水庫誘發地震預測方法研討[J].地震工程學報，2014，36(1)：171-176.LIU Xian-feng，LI Bi-xiong，DENG Jian-hui，et al.Discussion on the Prediction Methods for Reservoir-induced Seismicity[J].China Earthquake Engineering Journal，2014，36(1)：171-176.(in Chinese)

[15] 廖建忠.誘震八因子概率統計模型法在水庫誘發地震危險性概率評價中的應用[J].西北水電，2013(2)：11-14.LIAO Jian-zhong.Eight Factors Induce Earthquake Probability Statistical Model Method in the Application of Reservoir Induced Earthquake Risk Probability Assessment[J].Northwest Hydropower，2013(2)：11-14.(in Chinese)

Reservoir-induced Earthquakes for the Sanjiangkou Reservoir in Chongqing

WU Jian-chao1，2， ZHANG Li-fen1， WANG Xu-tao1，2， CAI Yong-jian2

(1.KeyLaboratoryofEarthquakeGeodesy，InstituteofSeismology，CEA，Wuhan，Hubei430071，China；2.InstituteofEarthquakeEngineeringofWuhan，Wuhan，Hubei430071，China)

Based on field investigations and previous research results，this study analyzes the tectonic background，seismic activity，and hydrogeology condition of the Sanjiangkou reservoir area.The Yushan fault is the main fault located in the southeast margin of the Sanjiangkou reservoir area.The other major fault is termed the Mawu fault，which is exposed in the west margin of the Sanjiangkou reservoir area.The active age of the Yushan fault is mid-Pleistocene and the Mawu fault is pre-Quaternary.The level of historic and current seismicity is low in the reservoir area.Employing structural analogy and probabilistic statistic methods，the probability of the induced earthquake in the Sanjiangkou reservoir area is discussed.These results are used to compare the Sanjiangkou reservoir with other reservoirs that have experienced earthquakes.The structural analogy method analysis showed that the possibility of tectonic-type reservoir-induced earthquakes is small following impoundment.However，there is the possibility of a karst subsidence-type reservoir-induced earthquake.This event could occur because of the existence of a large amount of thick layers of limestone in the reservoir area that could foster karst development.The loading and infiltration effects of water may easily cause rock collapse and crack propagation in the karst cave.This would induce a karst subsidence-type reservoir-induced earthquake.Considering the conditions of the Sanjiangkou reservoir area eight induced earthquake factors were selected to conduct the probabilistic prediction of the earthquake magnitude.These factors were reservoir water depth，lithology，regional state of stress，fault activity，degree of karst development，fault permeable depth，communication with the reservoir water，and seismic activity background.Based on the geological conditions and induced earthquake environment of the Sanjiangkou reservoir area，the reservoir area was divided into three segments: the head (Xintanzi-Shiliwan)，middle (Shiliwan-Niubizi)，and end (Niubizi-Xiamakou) regions.Using the MATLAB software programming，the probability of five earthquake magnitude was calculated for the three segments.The larger earthquake magnitude grading range was selected that corresponded to the relative probability value.The grading range was considered as the possible earthquake magnitude.The probabilistic statistic method analysis showed that the reservoir’s head region (Xintanzi-Shiliwan) had a small possibility of an induced earthquake with a possibility of earthquake occurrence of 0.02.The reservoir’s middle region (Shiliwan-Niubizi)，where the karst development is unlikely，had a possibility of induced earthquake of 0.96.The end region (Niubizi-Xiamakou) may induce microseismic activity with the possibility of inducedML(M< 3.0) of 0.10.

the Sanjiangkou reservoir； structural analogy method； probabilistic statistic method； reservoir-induced earthquake； possibility of earthquake occurrence

2014-06-19

中國長江三峽集團公司科研專項基金(SXSN/3354)

吳建超(1985-)，男，工程師，碩士，主要從事地震工程和水庫誘發地震研究.E-mail: jianchaowu85@gmail.com

P315

A

1000-0844(2015)02-0512-06

10.3969/j.issn.1000-0844.2015.02.0512