大渡河瀑布溝水電站水庫地震活動(dòng)特征分析

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

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

0 引 言

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

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

1.1 瀑布溝水電工程

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

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

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

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

水庫地震遙測(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ù)采集器，通過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ū)見圖1。

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

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

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

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

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

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

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

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

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

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

表2 2009～2017年瀑布溝水庫蓄水后地震事件統(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日水庫庫尾以西的石棉ML4.0地震。根據(jù)地震強(qiáng)度分析，水庫地震多屬微震，震級(jí)小，少數(shù)中等地震，地震強(qiáng)度呈現(xiàn)隨機(jī)變化特征，沒有明顯的周期性規(guī)律。

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

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

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

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

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

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

根據(jù)水庫地震活動(dòng)分布特征，對(duì)于庫尾鮮水河斷裂帶、安寧河斷裂帶、大涼山斷裂帶交匯部位的地震活動(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)系)，通過理論地震個(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]，通過對(duì)實(shí)際震級(jí)-頻度分布曲線與理論震級(jí)-頻度分布曲線的比較，發(fā)現(xiàn)當(dāng)震級(jí)不完整時(shí)，小震擬合效果存在著較大的差異。而GFT法正是通過對(duì)這些差異的分析，來判別地震目錄的完整性。首先假設(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)坪，白日照著，有一圈沒事可做的閑人，皆為一件小小熱鬧粘合在那里。

通過計(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)，沒有明顯的加速上升或減速變緩的趨勢(shì)，表明水庫蓄水后對(duì)地震數(shù)量沒有明顯的影響，地震與庫水位變化沒有直接相關(guān)的關(guān)系。

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

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

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

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

瀑布溝水庫附近地震活動(dòng)頻次呈均勻變化，水庫水位周期性升降變化與地震活動(dòng)頻次無明顯相關(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)為，水庫蓄水后主要為區(qū)域應(yīng)力場(chǎng)擾動(dòng)造成的中小型構(gòu)造地震，并非水庫蓄水引起的誘發(fā)地震。

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

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

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

3 結(jié) 語

雖然水庫蓄水后沒有明顯的誘發(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ù)研究，為水庫地震監(jiān)控和工程安全穩(wěn)定運(yùn)行提供技術(shù)保障。