王明振, 李 珊, 高 霖*

(1.重慶文理學(xué)院土木工程學(xué)院, 重慶 402160；2.重慶文理學(xué)院土木工程防災(zāi)減災(zāi)研究所, 重慶 402160)

長江三峽工程位于長江西陵峽中段，壩址在湖北宜昌三斗坪，是當(dāng)今世界上最大的綜合性水利樞紐工程，也是具有防洪、供水和航運等長江水資源開發(fā)利用及治理功能的長江水流控制性工程。

根據(jù)相關(guān)資料可知，長江三峽工程庫區(qū)(以下簡稱庫區(qū))及其附近存在仙女山斷裂、黔江斷裂、高橋斷裂、巫山-金佛山斷裂、九畹溪斷裂、都勻-豐都斷裂、大茅田-建始斷裂、長壽-遵義斷裂和華鎣山斷裂等斷裂帶。歷史上距長江三峽工程大壩壩址300 km 范圍內(nèi)共計發(fā)生4次6級左右地震，最大地震記錄為1856年黔江6.25級地震，震中烈度8度，地震形成7 600萬平方米堰塞湖。且開工前近40年內(nèi)，長江三峽工程及鄰區(qū)共記錄到2 000余次地震。綜合歷史震例數(shù)據(jù)表明，庫區(qū)是中國中等強度地震活動區(qū)[1]。

在長江三峽工程建設(shè)和運行的不同階段，庫區(qū)范圍內(nèi)地震發(fā)生頻次具有明顯變化。并且，觀測數(shù)據(jù)表明庫區(qū)水位高度變化亦會對地震發(fā)生頻次產(chǎn)生影響。圖1表示了長江三峽工程開建以來的歷年最高水位、各級地震發(fā)生年頻次變化情況[2]。

圖1 庫區(qū)歷年最高水位及地震發(fā)生年頻次Fig.1 The highest water level and earthquake frequency per year of reservoir region

由圖1可知，庫區(qū)地震主要表現(xiàn)為微震、極微震的活動形式，庫區(qū)范圍內(nèi)所發(fā)生最大震級地震為2013年12月16日湖北省恩施土家族苗族自治州巴東縣5.1級地震，震源深度5 km。雖然庫區(qū)強烈地震發(fā)生頻次不高，因地震災(zāi)害所產(chǎn)生的直接經(jīng)濟損失相對較小，但由于庫區(qū)存在大量條狀突出的山嘴、高聳孤立的山丘、非巖石和強風(fēng)化巖的陡坡、河岸和邊坡邊緣等不利地形，且上述不利地形在地震動作用下易產(chǎn)生地震動放大效應(yīng)，繼而易間接引起滑坡、崩塌、庫岸失穩(wěn)等地質(zhì)災(zāi)害。因此，有必要對長江三峽工程庫區(qū)范圍內(nèi)的地震活動性模型進行研究，以期對當(dāng)?shù)胤罏?zāi)減災(zāi)規(guī)劃、地震地質(zhì)災(zāi)害預(yù)防和預(yù)測提供參考和指導(dǎo)。

丁仁杰等[1]根據(jù)歷史震害資料對庫區(qū)地震災(zāi)害孕育條件、地震危險性、地震監(jiān)測預(yù)報以及震害預(yù)防等問題進行分析研究，提出了庫區(qū)蓄水對水庫地震、地質(zhì)災(zāi)害的潛在危險性等具體問題。王儒述[3-4]基于中外水庫誘發(fā)地震資料研究了水庫誘發(fā)地震的分布范圍、發(fā)震時間、發(fā)震趨勢、發(fā)震特點、誘發(fā)條件以及地震分類等問題，并根據(jù)長江三峽工程水庫的具體情況研究了庫區(qū)地震發(fā)生的危險性情況。戴苗等[5]針對長江三峽工程壩前水位變化情況與庫區(qū)地震活動性規(guī)律之間因果關(guān)系研究這一科學(xué)問題，結(jié)合壩前水位、壩前水位變化速率、地震活動性等監(jiān)測數(shù)據(jù)和地震現(xiàn)場調(diào)查資料，分析了庫區(qū)斷層滲透結(jié)構(gòu)與地震發(fā)生之間的聯(lián)系，總結(jié)分析出水庫地震發(fā)生頻次與壩前水位變化速率相關(guān)的結(jié)論，確定蓄水后庫區(qū)誘發(fā)地震多為礦山塌陷、滑坡型地震和巖溶型地震，構(gòu)造性地震占比不占主要地位。張燕等[6]通過分析三峽大壩壩址附近斷層活動性與庫區(qū)水位變化之間的關(guān)系，發(fā)現(xiàn)長江三峽水庫蓄水對周圍斷層活動具有較為明顯的影響，但這種影響主要是由水載荷產(chǎn)生的應(yīng)力場與斷層性狀耦合作用所致。楊陽等[7]利用遙感影像數(shù)據(jù)采用負荷格林函數(shù)積分對長江三峽工程庫區(qū)地面重力雖蓄水高度變化而變化情況進行了模擬計算，計算結(jié)果表明隨著水位上升大壩上游地面重力增大而下游呈減小趨勢。魏東等[8]從區(qū)域地質(zhì)背景、震中位置分布規(guī)律、庫區(qū)蓄水水位與地震發(fā)生頻次關(guān)系等角度入手，研究了三峽庫區(qū)巴東庫岸段地震震源機制解與震級-頻度關(guān)系系數(shù)。綜上可知，現(xiàn)有研究多集中在庫區(qū)水庫地震與蓄水水位的關(guān)系、庫首區(qū)域地震變遷規(guī)律、庫區(qū)地震的震源機制分析等方面，但對水庫范圍內(nèi)地震活動性規(guī)律等問題缺乏深入研究。

選取1996年1月—2017年12月期間發(fā)生在庫區(qū)范圍內(nèi)的地震為研究樣本，采用概率地震危險性評估理論，對庫區(qū)地震發(fā)生規(guī)律進行了分析研究。以期對庫區(qū)進行地震危險性分析和防震減災(zāi)工作提供指導(dǎo)和參考。

1 地震活動性模型

1968年Cornell提出了以概率來評價和表達一個地區(qū)或區(qū)域未來可能遭遇地震活動性大小的概率地震危險性方法(probabilistic seismic hazard analysis，PSHA方法)。PSHA方法自提出以來被廣泛應(yīng)用于地震小區(qū)劃、地震區(qū)劃圖以及工程地震易損性分析等研究工作[9]。

由于地震是隨機事件，即未來地震的發(fā)生時間、震中位置和震級大小等因素都具有不確定性，而PSHA方法建立地震危險性分析模型是須考慮地震發(fā)生時間、震中位置和震級大小的聯(lián)合概率分布，在建立聯(lián)合概率分布過程中，PSHA方法設(shè)定有三個基本假定：①潛在震源區(qū)內(nèi)地震發(fā)生時間服從泊松分布，且地震年平均發(fā)生率是常數(shù)；②潛在震源區(qū)內(nèi)地震發(fā)生的空間分布(震中位置)為均勻分布；③潛在震源區(qū)內(nèi)不同震級的發(fā)生頻度服從古登堡-里克特(Gutenberg-Richter，G-R)關(guān)系[9]。采用PSHA方法分析一個地區(qū)地震活動性時，需要用到的物理參數(shù)有：震級-頻度關(guān)系系數(shù)b、起算震級M0、震級上限Mu、地震年平均發(fā)生率ν等。由于本文所研究的地震活動區(qū)為長江三峽工程庫區(qū)，故在進行地震危險性分析過程中不考慮地震空間分布函數(shù)及其相關(guān)系數(shù)。

對于震級-頻度關(guān)系系數(shù)b的計算，根據(jù)G-R關(guān)系可按式(1)求得。

lgN=a-bM

(1)

式(1)中，M為震級；N為震級大于或等于M的地震次數(shù)；a和b為經(jīng)過統(tǒng)計分析得到的經(jīng)驗常數(shù)。其中b稱為震級-頻度關(guān)系系數(shù)，該系數(shù)與地震區(qū)帶內(nèi)的巖層應(yīng)力狀態(tài)和地殼巖石破裂強度有關(guān)。

起算震級M0是指潛在震源區(qū)內(nèi)可能發(fā)生的對工程場地產(chǎn)生影響的最小震級。根據(jù)研究目標(biāo)和潛在震源區(qū)的不同，M0應(yīng)根據(jù)實際情況確定。

震級上限Mu是指潛在震源區(qū)內(nèi)未來可能發(fā)生最大地震的震級。對于水庫區(qū)域內(nèi)最大震級Mu，本文的確定原則如下：取按式(2)計算確定出的數(shù)值[10]和區(qū)域范圍內(nèi)歷史最大地震震級向上取整所確定數(shù)值的較大值。

(2)

式(2)中，S為庫水面積，km2；H為最大庫深，m；V為庫容，106m3。

地震年平均發(fā)生率ν是指潛在震源區(qū)內(nèi)平均每年發(fā)生地震的次數(shù)，一般情況下為等于或大于起算震級M0的地震次數(shù)。

依據(jù)上述4個地震活動性參數(shù)，可計算得到潛在震源區(qū)震級分布概率、地震發(fā)生次數(shù)概率，計算公式如下：

(3)

式(3)中，β=bln10；ΔM為震級分檔間隔；Mj為第j個震級分檔。

(4)

式(4)中，n為一年內(nèi)地震發(fā)生次數(shù)。

2 庫區(qū)地震活動性概況

根據(jù)長江三峽工程誘發(fā)地震監(jiān)測系統(tǒng)連續(xù)監(jiān)測數(shù)據(jù)，《長江三峽工程生態(tài)與環(huán)境監(jiān)測公報1997—2018》[11]和《長江三峽工程運行實錄2003—2017》[12]分年度對庫區(qū)范圍內(nèi)不同震級分檔的地震發(fā)生次數(shù)進行了統(tǒng)計和分析，文獻[11-12]所公布庫區(qū)地震數(shù)據(jù)的震級分檔為0.0～0.9、1.0～1.9、2.0～2.9、3.0～3.9、4.0～4.9、5.0～5.9，1996年以來庫區(qū)內(nèi)各震級分檔的地震發(fā)生次數(shù)如圖2所示。

圖2 1996—2017年庫區(qū)不同級別地震年頻次圖Fig.2 Different levels earthquake frequency diagram of reservoir region in 1996—2017

由圖2可知，三峽工程庫區(qū)近期發(fā)生的地震主要以微震、極微震為主，最大地震為2013年12月16日13時4分發(fā)生在湖北省巴東縣的M5.1級地震。自2003年正式蓄水以來庫區(qū)地震發(fā)生頻次明顯升高，且2008年汛后開始175 m實驗性蓄水后地震活動年度總頻次達到峰值。

3 庫區(qū)地震震級-頻度關(guān)系模型

根據(jù)1996年1月1日—2017年12月31日地震震級與地震次數(shù)監(jiān)測數(shù)據(jù)樣本，按照式(1)可擬合得到震級-頻度關(guān)系系數(shù)b=0.795，系數(shù)a=4.304。擬合結(jié)果表明，擬合所得公式的計算結(jié)果與樣本數(shù)據(jù)的相關(guān)系數(shù)和均方差分別為0.999 999 999 6和0.097，結(jié)果表明所得參數(shù)具有較高可信度。

當(dāng)震級-頻度關(guān)系系數(shù)b=0.795時，則說明長江三峽工程庫區(qū)地震震級每增加一級，則地震發(fā)生次數(shù)約減少84%。

此外，從1996—2017年23年間的長江三峽庫區(qū)震級-頻度關(guān)系系數(shù)b=0.795來看，庫區(qū)內(nèi)地震的地震特征介于構(gòu)造地震與水庫誘發(fā)地震之間[13-15]，可定義為構(gòu)造“觸發(fā)型”水庫地震[8,16]。

4 庫區(qū)地震活動性規(guī)律初探

在建立長江三峽庫區(qū)地震活動性規(guī)律研究的過程中，需要確定起算震級、震級上限、震級分檔、地震年平均發(fā)生率。

起算震級M0取值為0級。

針對長江三峽工程庫區(qū)，S=1 084 km2、H=175 m、V=393 000×106 m3，則按式(2)進行計算可得Mu=3.47，而庫區(qū)近23年來發(fā)生最大地震震級為5.1級，按照前文所定原則最終確定Mu=6.0。

本文取震級分檔為0.2級。

根據(jù)圖2所示的地震統(tǒng)計數(shù)據(jù)可得庫區(qū)0、1.0、2.0、3.0、4.0、5.0級等不同震級分檔及以上地震年平均發(fā)生率ν0=681.5，ν1=146.8，ν2=41.6，ν3=3.7，ν4=0.5，ν5=0.087。

根據(jù)式(3)可得長江三峽工程庫區(qū)震級分布概率曲線，如圖3所示。

圖3 震級概率分布曲線Fig.3 Magnitude probability distribution curve

由圖3可知，長江三峽工程庫區(qū)M1級、M2級、M3級等震級分檔的地震發(fā)生概率分別為5.9%、0.95%、0.15%。

根據(jù)式(4)可得長江三峽工程庫區(qū)地震發(fā)生次數(shù)概率，如圖4所示。

圖4 地震年發(fā)生次數(shù)概率分布Fig.4 The distribution of earthquake occurrence probability per year

由圖4可知，對于長江三峽庫區(qū)不同震級1年發(fā)生頻次的概率估計結(jié)果如下：M5級及以上地震1年發(fā)生1次的概率為8.2%，實際平均每年發(fā)生0.087次；M4級及以上地震年發(fā)生1、2、3次的概率分別為30.3%、7.5%、1.3%，實際平均每年發(fā)生0.5次；M3級及以上地震1年發(fā)生3次的概率為20.9%，實際平均每年發(fā)生3.7次；M2級及以上地震1年發(fā)生41次的概率為6.2%，實際平均每年發(fā)生41.6次；M1級及以上地震1年發(fā)生146次的概率為3.3%，實際平均每年發(fā)生146.8次；M0級及以上地震1年發(fā)生681次的概率為1.5%，實際平均每年發(fā)生681.5次。

由以上結(jié)果可知，使用式(4)進行地震發(fā)生次數(shù)估計時，當(dāng)震級處于M5級及以上、M4級及以上和M3級及以上時，估計結(jié)果與實際統(tǒng)計結(jié)果較為吻合，但震級較小時估計結(jié)果與實際結(jié)果差別較大。因此，針對M2級及以下震級的地震發(fā)生頻率估計，不宜采用式(4)的泊松模型。

5 結(jié)論

針對長江三峽工程庫區(qū)地震發(fā)生規(guī)律，得到以下研究結(jié)論。

(1)總結(jié)了長江三峽工程自1996年開工以來至2017年12月31日之間庫區(qū)范圍內(nèi)地震發(fā)生情況，根據(jù)歷史地震記錄可知長江三峽工程庫區(qū)地震震級以微震與極微震為主，且?guī)靺^(qū)地震發(fā)生頻次與庫區(qū)蓄水水位的變化具有一定的正相關(guān)關(guān)系，因此未來應(yīng)注重監(jiān)測和研究庫區(qū)水位變化與庫區(qū)地震以及地質(zhì)災(zāi)害發(fā)生頻次變化之間的關(guān)系。

(2)采用古登堡-里克特關(guān)系模型估計出長江三峽工程庫區(qū)震級-頻度關(guān)系系數(shù)b=0.795，并根據(jù)參數(shù)b初步判斷長江三峽工程庫區(qū)地震以構(gòu)造“觸發(fā)型”水庫地震為主，構(gòu)造“觸發(fā)型”水庫地震的誘發(fā)與庫區(qū)消落區(qū)工程地質(zhì)條件和水位變化規(guī)律等因素相關(guān)。

(3)使用概率地震危險性方法，對長江三峽工程庫區(qū)的震級分布概率、地震發(fā)生次數(shù)概率等地震發(fā)生規(guī)律進行了研究，結(jié)果表明當(dāng)震級稍大時地震發(fā)生次數(shù)基本服從泊松分布，但極微震條件下地震發(fā)生次數(shù)不宜使用泊松分布進行概率估算。