段 煒,朱智嬌,劉 杰,莊儒新,李 圣
(1.云南省地震局,云南 昆明 650041;2.云南震科抗震技術工程公司,云南 昆明 650041;3.中國電建集團昆明勘測設計研究院有限公司,云南 昆明 650000)
自20 世紀60 年代印度柯依納、中國新豐江、津巴布韋卡里巴和希臘克里馬斯塔相繼發生破壞性水庫誘發地震以來,水庫地震監測有了極大改進[1]。對某一水庫的專項研究、水庫誘發地震震例與水庫誘發地震特征、成因機理研究隨之進行[2]。受資料積累影響,水庫誘發地震研究還處于初級階段,對其孕震過程及機理處于爭論階段[3]。Simpson 等[4]根據水庫地震和水庫水位變化之間的時間關系,將水庫蓄水后的地震分為快速響應型和滯后響應型兩類。馬文濤等[5]認為水庫誘發地震b值往往較大,多在1 左右,地震序列以1~3 級微震居多,但也有少部分水庫地震b值與天然地震相類似,因為b值的高低反映的是大小破裂展布在空間上的幾何特征。水庫庫區及周邊發生的地震活動屬于構造地震還是水庫誘發地震的判定仍是個普遍被關注和有待解決的科學問題[6-7]。就中國而言,由“西部大開發”推進的大規模水電站建設在地質構造復雜、構造穩定性較差的中國西部,水庫誘發地震已經成為可能影響電站建設和運行安全的不可忽視的重要問題[3,8]。
魯地拉水電站位于云南省大理州賓川縣和麗江市永勝縣交界處的金沙江中游河段上,是金沙江中游水電規劃“一庫八級”(龍盤、兩家人、梨園、阿海、金安橋、龍開口、魯地拉、觀音巖)的第七級電站。水電站壩高約140 m,庫容約17.18×108m3,正常蓄水位1223 m,回水至龍開口水電站,庫長約99 km,設計烈度為Ⅸ度。魯地拉水電站于2013 年6 月20 日開始蓄水。庫區數字地震監測流動臺網、固定臺網分別于2012 年1 月20 日、2013 年8 月1 日 投入運行。魯地拉水庫專用臺網記錄了水庫蓄水前后的地震活動,記錄資料表明蓄水后水庫部分庫段小震及微震頻度明顯增加,并于2014 年4月4 日在永勝縣發生M4.2 地震,距壩址約24 km,距庫岸約1 km。鑒于水電工程的重要性,有必要研究水庫蓄水前后地震活動特征。本文主要從魯地拉水電站所處區域地震構造與地震活動背景、庫區專用地震臺網監測到的蓄水前后地震活動及其與水庫水位關系等方面,對水庫蓄水前后地震活動特征進行分析,進而分析臨近庫岸地區M4.2 有感地震與水庫蓄水的關系,為研究水庫影響區(即水庫區及其外延10 km 范圍,下文統稱為“研究區”)地震活動特性及地震活動趨勢提供參考。
魯地拉水電站工程區大地構造位置處于揚子準地臺康滇地軸西緣,新構造運動以整體掀斜隆升與斷塊間相對運動為主,塊體內差異運動相對較弱。影響場地的主要發震構造為距壩址最近距離約為18 km 的程?!e川活動斷裂帶,歷史地震對場地最高影響烈度為Ⅶ度。程?!e川斷裂,全長200 km,總體走向近南北,地表顯微波狀彎曲。傾向南西,傾角70°~85°。第四紀以來,斷裂活動十分強烈,是一條全新世活動斷裂[9]。
魯地拉水庫大致可分為四段,壩址到金江橋庫首段為碎屑巖區,水庫蓄水后不致誘發地震;金江橋附近碳酸鹽巖庫段有誘發地震的可能性,誘發最高震級為5 級,對壩址不會造成破壞性影響;金江街到朵美以上為松散巖庫段,水庫蓄水后也不致誘發地震;庫尾玄武巖庫段,誘發地震的可能性很?、僦袊卣鹁值刭|研究所.金沙江魯地拉水電站工程場地地震安全性評價和水庫誘發地震評價報告.北京:中國地震局地質研究所,2004.。
根據《中國歷史強震目錄(公元前23 世紀—公元1911 年)》[10]、《中國近現代地震目錄(公元1912—1990 年MS≥4.7)》[11]以及云南省地震目錄,自公元886 年至2017 年12 月魯地拉水庫 監測區(99.8127°~101.8173°E,25.3017°~27.1003°N)內發生了108 次M≥4.7 地震,其中4.7≤M<5.0 地震23 次,5.0≤M<6.0 地震63 次,6.0≤M<7.0 地震19 次,7.0≤M<8.0 級地震3 次。1515 年6 月27 日永勝西北紅石崖的7級、1803 年2 月2 日賓川州城附近6級、1959 年3月30日賓川東北5.5級、1959年4月26日賓川北5級、1964 年2 月13 日 賓 川 拉 烏5.4 級、1964 年5 月3 日 賓川 東 北4.7 級、1992 年12 月18 日永勝果園、東龍潭5.4 級、1992 年12 月22日永勝河口、清水5.1 級、2001 年10 月27 日永勝濤源6.0 級、2006 年1 月28 日永勝賓川交界東風鄉4.7 級、2009 年11 月2 日賓川5.0 級地震,均發生在程?!e川斷裂帶附近,距壩址分別約57 km、55 km、25 km、12 km、37 km、35 km、26 km、28 km、22 km、0.23 km、25 km。距壩址最近地震是2006 年1 月28 日發生在永勝賓川交界的東風鄉4.7級地震,該地震距壩址僅約230 m。近期最大地震為2009 年11 月2 日發生在賓川的5 級地震,該地震距壩址約25 km。研究區內發生了5 次M≥4.7 地震,其中4.7≤M<5.0 地震2 次,5.0≤M<6.0 地震2 次,6.0≤M<7.0地震1次,無M≥7.0地震活動。最大地震為2001年10 月27 日永勝濤源的6.0 級地震,該地震距壩址約22 km,距庫岸約5 km,震中烈度為Ⅶ度。
根據云南和四川兩省區域地震臺網觀測數據,1965—2017年魯地拉研究區在49年期間記錄到儀測小地震2.0≤M≤4.6 地震184 次,其中2.0≤M<3.0 地震151 次,3.0≤M<4.0 地震29 次,4.0≤M≤4.6 地震4 次,地震主要分布在程?!e川斷裂帶附近,大部分地震距金沙江庫岸0~3 km。從時間上分析,研究區范圍內約每12 年會發生一次4.5級左右的有感地震。最近一次為2014年4月4日發生在庫中段的永勝M4.2地震。
魯地拉水電站于2013 年6 月20 日開始蓄水。因蓄水前水庫專用臺網記錄的時間較短,故選用了全國地震編目系統提供的地震目錄作為補充,時間為2008 年1 月—2012 年1 月,范圍為研究區ML≥0 地震142 次;魯地拉水庫專用臺網產出的地震目錄,時間為2012 年2 月至2017 年12 月,范圍為研究區ML≥-0.8 地震2584次。將以上目錄合并,根據蓄水時間進行分段統計,將蓄 水 前2008 年1 月1 日—2013 年6 月19 日研究區內的168 次地震,蓄水后2013 年6月20 日—2017 年12 月31 日研究區內的2558 次地震作為此次研究樣本。
2008 年1 月—2017 年12 月研究區內共發生地震2726次,其中ML≤0.9地震1858次,1.0≤ML<2.0 地 震777 次,2.0≤ML<3.0 地 震78 次,3.0≤ML<4.0 地震11 次,4.0≤ML<5.0 地震2 次,無5.0級以上地震。最大地震為2014 年4 月4 日發生在庫中段的永勝M4.2地震。
圖1 魯地拉水庫研究區蓄水前后地震分布Fig.1 Earthquake distribution before and after the water storage in Ludila reservoir of the study area
蓄 水 前2008 年1 月1 日—2013 年6 月19 日近五年半的時間里,研究區內共記錄到地震168 次,日均發生地震0.08 次,其中ML≤0.9 地震0 次,1.0≤ML<2.0 地震161 次,2.0≤ML<3.0 地震7次,無3.0級以上地震。最大地震為2013年6月10日發生在庫尾的鶴慶ML2.8地震,震源深度約10 km,距壩址約59 km,距庫岸約8 km。其它微震震源深度大部分為15 km 左右,表現為構造地震的特征。圖1 為研究區蓄水前后地震震中分布,由圖1(a)可知,蓄水前地震分布呈現為零散的自然分布,大部分地震分布在庫尾段。
蓄 水 后2013 年6 月20 日—2017 年12 月31日近四年半的時間里,研究區內共記錄到地震2558 次,日均發生地震2.34 次,與蓄水前相比地震頻次有較為明顯的增加,其中ML≤0.9 地震1858次,1.0≤ML<2.0地震616次,2.0≤ML<3.0地震71 次,3.0≤ML<4.0 地震11 次,4.0≤ML<5.0 地震2次,無5.0級以上地震。最大地震為2014年4 月4 日發生在庫中段的永勝M4.2 地震,震源深度約10 km,距壩址約24 km,距庫岸約1 km。蓄水后記錄了大量ML≤0.9 地震,說明研究區監測能力較蓄水前大大提高。從空間上看,由圖1(b)可知,蓄水后研究區內微小地震位置距庫岸約0~5 km,震源深度大部分在0~10 km,約占地震總數的75%以上,符合水庫誘發地震的“雙十”特征。庫壩區和庫中段地震活動明顯增多,尤其是庫中段地震呈現為叢集分布,且沿程?!e川斷裂展布,此庫段為碳酸鹽巖庫段,據魯地拉水電站工程地質條件評價報告分析,此庫段水庫蓄水后存在誘發地震的可能性,預測水庫誘發地震的震級為5.0級左右。
利用2008年1月1日—2017年12月31日魯地拉研究區地震目錄,采用G-R 關系計算了在此期間魯地拉研究區的b值,約為0.8431(圖2),相關系數為0.996,相關程度較高,最小完整性震級MC約為0.5 級。與李濤等12]基于分區空間掃描研究云南地區b值分布特征里計算得到的此區域的b值平均值為0.805 相比,b值有所增大,可能與水庫蓄水之后小震數量明顯增加有關。
圖2 研究區地震震級頻次關系Fig.2 The relationship between earthquake magnitude and frequency in the study area
圖3a、3c、3e 分別為研究區地震活動M-t圖、月頻度、日頻度與水位時序變化。由圖可見,蓄水前研究區內微小地震活動日均發生地震0.08 次,蓄水后研究區內微小地震活動較為頻繁,日均發生地震2.34 次,日頻度增加了近30倍。從水位曲線可知,從2013年6月19日開始蓄水,2013 年7 月蓄水至1210 m,維持此水位至2013 年11 月15 日,之后水位迅速下降,2014 年2 月6 日回落至蓄水前的水位,之后迅速上升,2014 年6 月3 日蓄水至1219 m,2014年8 月13 日蓄水至最高水位1223 m,2014 年6—8 月為汛期,之后一直維持在最高水位至2015 年6 月,之后進入汛期,有10 m 的落差,至2015 年8 月汛期結束后維持最高水位,2016年和2017 年的水位除了汛期有約10 m 的波動外,其余月份均維持在最高水位。從頻度圖可見,蓄水初期,2013 年6 月—2014 年3 月地震活動有所上升,但與蓄水前相差不大。2013 年11 月—2014 年6 月,水位先迅速下降至最低水位后又迅速上升至高水位,在迅速上升的中期地震活動頻次大幅增加,2014 年4 月4 日在庫中段發生M4.2 地震,最大日頻次為當日的101次。隨著水位的持續升高—降低—升高變化并維持在高水位,地震活動頻次和強度相對穩定,經過三年半的蓄水及高水位,研究區內微小構造已達到相對的穩定,地震活動強度維持在ML3.0~4.0 左右的水平,較蓄水前有增強的變化,在汛期仍然稍有波動。由月頻度圖亦可知,2014 年3 月開始月頻度開始增加,由蓄水前的不到3 次增加到70 次左右,最大月頻次為2014年4月的433次。由此可見,研究區地震活動與水庫水位變化有較為明顯的相關性。
圖3 魯地拉研究區蓄水前后地震M-t圖、月頻度、日頻度、震源深度分布與水位時序變化Fig.3 The relationship between M-t diagram,monthly frequency,daily frequency,focal depths distribution and the water level before and after water storage in Ludila reservoir of the study area
圖3(g)為研究區地震震源深度分布。由圖可見,魯地拉研究區蓄水后地震震源深度在5~10 km 的地震數量明顯增加,約占地震總數的53.64%,蓄水前僅為13.4%。蓄水前的地震震源深度主要在10~15 km,約占39.6%。由此可見,水庫蓄水對研究區10 km 深度范圍內的地震有顯著影響。圖4 為研究區地震震源深度沿剖面分布,由圖可知,庫尾段地震震源深度主要分布在10~20 km,庫中段微小地震震源深度主要分布在0~10 km,庫首段地震震源深度在5~15 km均有分布。
圖4 研究區地震震源深度沿剖面分布Fig.4 Distribution of focal depths along profile in study area
魯地拉水庫專用臺網庫區的地震監測能力可達到0.5級,但區域地震臺網的地震監測能力相對較弱。劉麗芳等[13]在云南地區地震目錄最小完整性震級研究中計算得到,云南“十五”數字化臺網建設階段(2007—2011 年),臺站數量由57 增加到63,云南省MC可達ML1.5 左右,如大理—麗江—鹽源區MC可控制在ML1.0~1.5。為了更準確地對比蓄水前后地震活動情況,對研究區內ML≥1.5 地震進行分析。2008 年1 月—2017 年12 月 研 究 區 共 發 生ML≥1.5 地 震330 次,其中1.5≤ML<2.0 地震239 次,2.0≤ML<3.0 地震78 次,3.0≤ML<4.0 地震11 次,4.0≤ML<5.0 地震2 次。蓄水前共發生ML≥1.5 地震98 次,日均發生地震0.05 次,其中1.5≤ML<2.0 地震91次,2.0≤ML<3.0 地震7 次。蓄水后共發生ML≥1.5地震232次,日均發生地震0.14次,其中1.5≤ML<2.0地震148次,2.0≤ML<3.0地震71次,3.0≤ML<4.0 地震11 次,4.0≤ML<5.0 地震2 次。蓄水后日頻次是蓄水前的2.8 倍。由圖1b 可見,蓄水后的地震主要分布在庫中段,程?!e川斷裂兩側。采用G-R 關系計算b值約為0.8052,與李濤等[12]計算得到的b值相同,屬正常構造地震活動。從圖3b、3d、3f 可見,蓄水后地震強度和頻度均有所上升,但與蓄水前差異不是特別明顯,并趨于穩定,但隨汛期有微弱的變化。從圖3h 可知,蓄水前10~15 km 深度范圍內的地震占36.08%,16~20 km 深度范圍內的地震占28.87%,蓄水后10~15 km深度范圍內的地震占56.9%,16~20 km 深度范圍內的地震占4.31%,蓄水前后地震均主要集中在10~20 km深度范圍內,由此可見,蓄水后研究區ML≥1.5地震主要反映為構造地震的特征。
雙差地震定位是一種相對定位方法,該方法利用信號的走時差反演震源位置,能夠有效地消除震源至臺站共同傳播路徑效應,對地殼速度模型的依賴性相對較小。它是由Waldhauser 和Ellsworth(2000)提出[14],目前該方法已經被國內外地震學家廣泛地應用到實際的地震定位中,并取得了很好的效果。
2014年3月28日10時—4月4日8時在云南省永勝縣持續記錄到微小地震,之后在2014 年4 月4 日9 時在庫中段的永勝縣發生M4.2 地震,距庫岸約1 km。2014 年3 月28 日—4 月30 日魯地拉水庫臺網記錄到446 次地震,其中ML≥0.5地震282 次。為了更準確地研究其震源位置,采用雙差定位方法對該段時間ML≥0.5 地震進行重定位,為了保證定位結果的可靠性,選擇了4個以上臺站定位的地震事件共266次。模型采用了吳建平等[15]計算得出的結果,如圖5 所示。圖6a、6b 為重新定位前后的震中分布,由圖可知,重新定位后的地震在空間分布上有更為明顯的成帶性,并沿著程?!e川斷裂的走向展布。圖6c、6d 為重新定位前后的震源深度分布,表1 為重新定位前后震源深度統計。由圖6c、6d 及表1 可知,重新定位前的地震主要分布在10~15 km,占比為59.77%。重新定位后的地震主要分布在5~10 km,占比為78.95%,其中絕大部分震源分布在7 km 附近,沒有隨時間發生較大的變化,10~17 km 范圍內也有零散分布。另外,重新定位后發現13~17 km 范圍內地震較之前有所增加。
表1 雙差定位前后震源深度統計(2014-03-28—2014-04-30)
圖5 雙差定位所用速度模型Fig.5 Velocity model used for double-difference location
圖6 利用雙差定位方法重新定位前后地震震中分布及震源深度分布Fig.6 Distribution of the seismic and focal depths before and after relocating by using double-difference location method
圖7 地震震源深度沿A1-A2、B1-B2剖面分布Fig.7 Distribution of focal depths along profile A1-A2 and B1-B2
由圖7 可知,地震震源沿A1-A2剖面在5~17 km 范圍內均有分布,斷層傾角約為70°,傾向南西,與安曉文等[9]結果基本一致。由地震震源沿B1-B2剖面的分布可見,斷裂東北段較淺,向西南段逐漸加深,表明距離庫岸近的位置淺源地震多,隨著距離庫岸位置增大,淺源地震逐漸減少,說明離庫岸越遠,水的滲透影響越小。
魯地拉水電站位于滇西北地震活動構造較為復雜的區域內,研究區內主要發震構造為距壩址最近距離約為18km的程?!e川活動斷裂帶,此斷裂帶歷史上最大地震為1515 年永勝西北紅石崖的級。近期最大地震為2009 年11 月2 日發生在賓川的M5.0 地震,該地震距壩址約25 km。據歷史地震統計,從時間上分析,研究區約每12 年會發生一次4.5 級左右的有感地震,最近四次分別發生在1975 年、1987 年、2001年和2014年。
通過蓄水前后對比分析可知,水庫蓄水對研究區內地震活動影響明顯,蓄水后日頻次是蓄水前的30 倍,特別是記錄了大量ML≤0.9 地震,其震源深度0~10 km 范圍占83.64%,說明研究區監測能力較蓄水前大大提高,亦說明ML≤0.9 地震可能為水庫誘發地震。蓄水前震源深度大部分為10~15 km,約占39.6%,蓄水后地震震源深度在5~10 km 范圍內的地震數量明顯增加,約占地震總數的53.64%,由此可見,水庫蓄水對研究區10 km 深度范圍內的地震有顯著影響。經過三年半的蓄水及高水位,研究區內微小構造已達到相對的穩定,地震活動強度維持在ML3.0~4.0級左右的水平,但在汛期仍然稍有波動。
對ML≥1.5 地震進行分析可知,蓄水后日頻次是蓄水前的2.8倍。蓄水后的地震主要分布在庫中段,沿程?!e川斷裂兩側展布。蓄水后地震強度和頻度均有所上升,但與蓄水前差異不是特別明顯,并趨于穩定,但隨汛期有微弱的變化。從震源深度分布圖可知,蓄水前后地震均主要集中在10~20 km 深度范圍內,由此可見,蓄水后研究區ML≥1.5 地震主要反映為構造地震的特征。
b值是震級-頻度線性關系中的斜率,即變化比例。研究表明b值隨時間與地區變化,以及與地震間的對應關系。每個地震活動區在一定時間內都有確定的b值,它隨時間變化。本文計算了研究區內b值為0.8431。
魯地拉水電站庫中段為碳酸鹽庫段,程?!e川斷裂在此庫段通過水庫,具有發生地震的應力條件,從構造、巖性、以及程?!e川斷裂帶的新活動性等綜合考慮,此庫段水庫蓄水后存在誘發地震的可能性。蓄水后,研究區內發生的地震呈叢集分布,主要分布于此庫段,并在水位急劇下降又急劇上升后,在庫中段發生M4.2 有感地震。但從歷史上看,研究區約每12 年會發生一次4.5 級左右的有感地震,此次有感地震發生時間雖與水庫蓄水相關,但也可能是水庫蓄水引起該庫段的溶巖作用和斷層活動共同觸發了構造地震。