黃 會(huì) 寶， 莊 園 旭， 江 德 軍， 宋 澄

(1.國(guó)能大渡河流域水電開(kāi)發(fā)有限公司，四川 成都 610041；2.四川省地震局，四川 成都 610000；3.四川大學(xué)水利水電學(xué)院，四川 成都 610065)

0 引 言

水庫(kù)誘發(fā)地震是由于水庫(kù)蓄水或水位變化而引發(fā)的一種獨(dú)特的地震類型，主要有水庫(kù)蓄水而引起的地震和在水庫(kù)鄰近地區(qū)的地震。水庫(kù)誘發(fā)地震問(wèn)題始于1940 年，特別是1962 年中國(guó)新豐江水庫(kù)6.1 級(jí)、1963 年贊比亞-津巴布韋卡里巴(Kariba)水庫(kù)6.25級(jí)、1966年希臘克里瑪斯塔(Kremasta)水庫(kù)6.3級(jí)和1967年印度柯依納(Koyna)水庫(kù)6.5級(jí)地震，引起了人們足夠的重視[1]。水庫(kù)誘發(fā)地震研究至今已有80年，不少地震研究學(xué)者對(duì)水庫(kù)誘發(fā)地震研究結(jié)果也有不同見(jiàn)解，對(duì)地震是否與水庫(kù)蓄水有關(guān)的數(shù)據(jù)統(tǒng)計(jì)也有差異，但對(duì)部分水庫(kù)蓄水后會(huì)誘發(fā)地震活動(dòng)的事實(shí)達(dá)成了共識(shí)[2]。在分析水庫(kù)誘發(fā)地震機(jī)理時(shí)，學(xué)者認(rèn)為庫(kù)水對(duì)環(huán)境的載荷作用、孔隙壓力擴(kuò)散作用以及水的潤(rùn)滑作用是庫(kù)區(qū)附近發(fā)生地震的主要因素，揭示了地震與水庫(kù)蓄水的內(nèi)在聯(lián)系，提出在實(shí)際水庫(kù)誘發(fā)地震監(jiān)測(cè)和判識(shí)中應(yīng)分析水庫(kù)地震活動(dòng)特征。

1 工程及地質(zhì)背景

1.1 瀑布溝水電工程

大渡河瀑布溝水電站位于四川省漢源和甘洛兩縣交界處，為大渡河下游控制性水電工程。電站以發(fā)電為主，兼有防洪、攔沙等綜合利用的一等大(I)型工程，工程主要由礫石土心墻堆石壩，引水發(fā)電建筑物、溢洪道、泄洪洞，放空洞和尼日河引水建筑物組成。大壩最大壩高186 m，壩頂長(zhǎng)540.5 m，壩頂寬度14 m，水庫(kù)回水長(zhǎng)度72 km，正常蓄水位850.00 m，死水位790.00 m，消落深度60 m，總庫(kù)容53.37億m3，為不完全年調(diào)節(jié)水庫(kù)，電站裝機(jī)容量3 600 MW，大壩設(shè)防烈度8度。

1.2 庫(kù)區(qū)地質(zhì)構(gòu)造

水庫(kù)地處青藏高原與四川盆地的過(guò)渡地帶，位于揚(yáng)子地臺(tái)西部古老的康滇地軸北段東側(cè)，與現(xiàn)代強(qiáng)烈活動(dòng)的川滇菱形地塊毗鄰。水庫(kù)周圍分布了多條大型活動(dòng)斷裂帶，包括水庫(kù)以北的龍門山斷裂帶(2008年5月12發(fā)生了汶川8.0級(jí)地震)，水庫(kù)西北的鮮水河斷裂帶(包括1955年康定7.5級(jí)和1973年?duì)t霍7.6級(jí)在內(nèi)的多次7級(jí)以上地震)以及水庫(kù)以南的安寧河-則木河斷裂和大涼山斷裂，古地震考察顯示安寧河-則木河斷裂承載了多次7級(jí)地震的平均位錯(cuò)量[3]，大涼山斷裂帶的古地震震級(jí)也在7級(jí)以上[4]。從近場(chǎng)看庫(kù)區(qū)附近主要發(fā)育NW、NE向斷裂構(gòu)造，主要在第四系沉積區(qū)上，金坪斷裂穿過(guò)水庫(kù)蓄水區(qū)域，附近區(qū)域分布有美羅斷裂、順河斷裂、鑼鍋坪斷裂、紅花斷裂、黑馬鄉(xiāng)斷裂、烏斯河斷裂等多條小規(guī)模斷裂，自有記錄以來(lái)未發(fā)生過(guò)中等以上地震。

1.3 水庫(kù)地震監(jiān)測(cè)設(shè)施

水庫(kù)地震遙測(cè)臺(tái)網(wǎng)分布在北緯29.0°～29.6°N，東經(jīng)102.3°～103.1°E區(qū)域內(nèi)，共設(shè)置13個(gè)子臺(tái)和1個(gè)臺(tái)網(wǎng)中心。目前采用短周期速度地震計(jì)和高分辨率的數(shù)據(jù)采集器，通過(guò)CDMA路由器傳輸至臺(tái)網(wǎng)中心。本臺(tái)網(wǎng)監(jiān)測(cè)震級(jí)下限為ML0.5，臺(tái)站分布及臺(tái)網(wǎng)監(jiān)測(cè)能力、重點(diǎn)監(jiān)視區(qū)見(jiàn)圖1。

(a)臺(tái)站分布和監(jiān)測(cè)能力 (b)地震重點(diǎn)監(jiān)視區(qū)示意圖圖1 瀑布溝水庫(kù)地震臺(tái)網(wǎng)監(jiān)測(cè)能力及重點(diǎn)監(jiān)測(cè)區(qū)示意圖

2 地震活動(dòng)特征

分析水庫(kù)蓄水區(qū)域約10 km范圍內(nèi)地震情況，資料來(lái)源于水庫(kù)地震遙測(cè)臺(tái)網(wǎng)地震信息和地震活動(dòng)事件，其中蓄水前是2008～2009年水庫(kù)蓄水期間；蓄水后是指2009～2017年底的地震信息。

2.1 蓄水前的地震活動(dòng)特征

水庫(kù)蓄水前共記錄到地震事件1077次，主要發(fā)生在庫(kù)區(qū)中段漢源縣城附近，地震事件統(tǒng)計(jì)見(jiàn)表1。庫(kù)區(qū)中段漢源縣城附近在蓄水前已經(jīng)存在密集小震活動(dòng)，主要集中在0.0

表1 2008～2009年水庫(kù)蓄水前地震事件統(tǒng)計(jì)表

圖2 瀑布溝水庫(kù)中段漢源縣附近地震發(fā)生時(shí)段特征

2.2 蓄水后的地震活動(dòng)特征

2.2.1 地震活動(dòng)強(qiáng)度特征

水庫(kù)2009年11月1日下閘蓄水以來(lái)，共記錄到地震事件3 290次，平均日發(fā)生地震1.1次，最高日發(fā)生地震13次，年均發(fā)生地震事件數(shù)量較蓄水前有明顯減少，蓄水后地震事件統(tǒng)計(jì)見(jiàn)表2。

表2 2009～2017年瀑布溝水庫(kù)蓄水后地震事件統(tǒng)計(jì)表

地震活動(dòng)主要集中在ML≤2.9級(jí)范圍，占全部地震總數(shù)的96.44%，ML≥3.5級(jí)地震分別在2011年初、2012年中和2015年共發(fā)生3次，4.0≤ML≤4.9級(jí)地震僅發(fā)生1次，最大地震為2015年4月15日水庫(kù)庫(kù)尾以西的石棉ML4.0地震。根據(jù)地震強(qiáng)度分析，水庫(kù)地震多屬微震，震級(jí)小，少數(shù)中等地震，地震強(qiáng)度呈現(xiàn)隨機(jī)變化特征，沒(méi)有明顯的周期性規(guī)律。

2.2.2 地震活動(dòng)空間特征

水庫(kù)蓄水后地震活動(dòng)頻次不高，主要分布以庫(kù)尾的鮮水河斷裂帶、安寧河斷裂帶、大涼山斷裂帶交匯部位以及庫(kù)區(qū)中段漢源縣城附近，對(duì)庫(kù)區(qū)地震分布密度做空間分布統(tǒng)計(jì)(見(jiàn)圖3)，兩個(gè)地震密度相對(duì)較高的叢集區(qū)域滿足完整性震級(jí)條件的地震每km2約3～10次，其中地震密度相對(duì)較高的石棉縣附近，沿安寧河斷裂帶呈NNW-SSE向展布，整體與斷裂帶走向一致；在漢源縣城附近地震分布集中且基本在縣城區(qū)域，疑似人工生產(chǎn)活動(dòng)的可能性較高，地震分布見(jiàn)圖3。

圖3 瀑布溝水庫(kù)庫(kù)區(qū)地震活動(dòng)及地震密度分布圖(五角星為ML≥4.0地震；三角形為遙測(cè)地震臺(tái)站)

對(duì)水庫(kù)地震SW-NE和NW-SE兩個(gè)垂直剖面(見(jiàn)圖4)進(jìn)行深度分析，地震深度基本小于30 km范圍內(nèi)。其中SW-NE向剖面將兩個(gè)叢集區(qū)地震分離開(kāi)來(lái)，明顯看出深度差異，水庫(kù)庫(kù)尾附近叢集地震活動(dòng)深度在0～30 km范圍內(nèi)，相對(duì)于其他區(qū)域明顯較深；庫(kù)中段漢源縣城叢集地震深度在小于15 km范圍內(nèi)，水庫(kù)以東零散地震活動(dòng)大部分地震深度在5～25 km范圍內(nèi)。

圖4 瀑布溝水庫(kù)庫(kù)區(qū)內(nèi)地震深度垂直剖面投影

2.2.3 地震活動(dòng)的時(shí)間特征

根據(jù)水庫(kù)地震活動(dòng)分布特征，對(duì)于庫(kù)尾鮮水河斷裂帶、安寧河斷裂帶、大涼山斷裂帶交匯部位的地震活動(dòng)處于臺(tái)網(wǎng)監(jiān)測(cè)邊緣，該區(qū)域理論監(jiān)測(cè)水平覆蓋能力為1.0級(jí)。為更客觀地研究地震活動(dòng)規(guī)律，利用前期學(xué)者研究的ZMAP程序[5]中“擬合優(yōu)度檢驗(yàn)法(Goodness of Fit Test，簡(jiǎn)寫為GFT法)”，按最大似然法原理求解b值(震級(jí)-頻度關(guān)系)，通過(guò)理論地震個(gè)數(shù)和實(shí)際地震個(gè)數(shù)進(jìn)行擬合，認(rèn)為擬合優(yōu)度值越大地震目錄完整性越高，并據(jù)此分析地震目錄的最小完整震級(jí)Mc，盡可能地使用最多的地震記錄進(jìn)行分析。

擬合優(yōu)度檢驗(yàn)法(GFT 法)是在地震的震級(jí)-頻度在較大范圍內(nèi)滿足古登堡-里克特(G-R)震級(jí)頻率關(guān)系分析方法(lgN=a-bM)的基礎(chǔ)上提出的[6]，通過(guò)對(duì)實(shí)際震級(jí)-頻度分布曲線與理論震級(jí)-頻度分布曲線的比較，發(fā)現(xiàn)當(dāng)震級(jí)不完整時(shí)，小震擬合效果存在著較大的差異。而GFT法正是通過(guò)對(duì)這些差異的分析，來(lái)判別地震目錄的完整性。首先假設(shè)地震目錄自Mi是完整的，采用G-R公式擬合出該震級(jí)Mi所對(duì)應(yīng)的ai和bi值，將理論值與實(shí)際值進(jìn)行對(duì)比得出各個(gè)震級(jí)所對(duì)應(yīng)的擬合優(yōu)度值R：

式中Bi是指震級(jí)大于等于Mi的實(shí)際觀測(cè)地震個(gè)數(shù)；Si是理論推測(cè)的地震個(gè)數(shù)，Si=10ai-biM

北京城什剎海前海南頭，煤灰土新墊就一片場(chǎng)坪，白日照著，有一圈沒(méi)事可做的閑人，皆為一件小小熱鬧粘合在那里。

通過(guò)計(jì)算，實(shí)際產(chǎn)出監(jiān)測(cè)目錄的最小完整性震級(jí)Mc為0.9級(jí)，b值為0.872(圖5)。以最小完整性震級(jí)Mc為0.9級(jí)作為震級(jí)統(tǒng)計(jì)下限，可以看到地震累計(jì)發(fā)生頻度呈現(xiàn)均勻變化(圖6)，沒(méi)有明顯的加速上升或減速變緩的趨勢(shì)，表明水庫(kù)蓄水后對(duì)地震數(shù)量沒(méi)有明顯的影響，地震與庫(kù)水位變化沒(méi)有直接相關(guān)的關(guān)系。

圖5 瀑布溝水庫(kù)地震目錄擬合關(guān)系曲線

圖6 瀑布溝水庫(kù)地震累計(jì)發(fā)生頻度及時(shí)間關(guān)系

2.2.4 地震震源機(jī)理特征

水庫(kù)誘發(fā)地震主要有正傾錯(cuò)動(dòng)和走滑錯(cuò)動(dòng)兩種類型，近年來(lái)學(xué)者對(duì)高烈度地區(qū)水庫(kù)誘發(fā)地震的研究顯示逆沖型地震也占據(jù)了一定比例，因此，水庫(kù)地震錯(cuò)動(dòng)類型主要取決于庫(kù)區(qū)附近已有斷裂或裂隙的構(gòu)造特征。

瀑布溝水庫(kù)附近地震活動(dòng)頻次呈均勻變化，水庫(kù)水位周期性升降變化與地震活動(dòng)頻次無(wú)明顯相關(guān)性，地震活動(dòng)集中區(qū)域主要為石棉縣城附近，該區(qū)域中小地震以走滑地震為主，正傾和逆沖型地震稀少，與該區(qū)域主要斷裂構(gòu)造特征一致。利用該區(qū)域中小地震反演的應(yīng)力張量顯示主壓應(yīng)力為進(jìn)水平的NWW-SEE向，與大區(qū)域應(yīng)力場(chǎng)一致。因此，綜合分析認(rèn)為，水庫(kù)蓄水后主要為區(qū)域應(yīng)力場(chǎng)擾動(dòng)造成的中小型構(gòu)造地震，并非水庫(kù)蓄水引起的誘發(fā)地震。

2.3 地震活動(dòng)總體情況

水庫(kù)蓄水且水位快速上升至死水位期間，地震頻次較后續(xù)時(shí)段偏高(見(jiàn)圖7)，出現(xiàn)地震能量釋放階段，但震級(jí)均小于ML1.5，地震均為微小地震，后續(xù)地震頻次有所減少，各震級(jí)分布較為均勻，無(wú)突出高震集中分布，總體上水庫(kù)周邊地震活動(dòng)水平很低。

圖7 瀑布溝水庫(kù)地震發(fā)生頻度與庫(kù)水位關(guān)系

3 結(jié) 語(yǔ)

雖然水庫(kù)蓄水后沒(méi)有明顯的誘發(fā)地震現(xiàn)象，但附近存在大型構(gòu)造斷裂帶存在發(fā)生中強(qiáng)地震威脅的可能，仍需加強(qiáng)地震監(jiān)測(cè)工作，密切關(guān)注地震活動(dòng)狀態(tài)，及時(shí)整理分析地震監(jiān)測(cè)成果，同時(shí)加大地震監(jiān)測(cè)預(yù)測(cè)技術(shù)研究，為水庫(kù)地震監(jiān)控和工程安全穩(wěn)定運(yùn)行提供技術(shù)保障。