水庫區(qū)地震精確定位研究

莊園旭　楊正剛　杜興忠　江曉濤　胡興堯

摘 要：利用北盤江地震臺網(wǎng)2012年2月至2012年12月記錄到的1513次地震的震相目錄數(shù)據(jù)，采用雙差定位法進(jìn)行精確定位，定位后到時均方根殘差減小至0.036s，獲得了北盤江流域1136次地震事件的精確定位結(jié)果。精確定位結(jié)果顯示，地震分為兩個地震叢集，呈現(xiàn)出條帶狀分布，震源深度主要集中在4-12km，震源深度大于5km的地震占83%，說明該區(qū)域地震以構(gòu)造地震為主。

關(guān)鍵詞：雙差定位；地震定位；精確定位

引言

北盤江屬珠江水系西江左岸支流。發(fā)源于云南省宣威馬雄山西北麓，流經(jīng)滇東、黔西南，于貴州省望謨縣蔗香雙江口與南盤江匯合后稱紅水河。北盤江中下游已建成光照水電站于2007年12月30日下閘蓄水、董箐水電站于2009年08月20日下閘蓄水。區(qū)內(nèi)主要斷裂有西部北西的埡都-紫云深斷裂、威寧-水城斷裂和北東向的石門坎斷裂，東部北東向的松桃-獨(dú)山斷裂，西南部的北東向的開遠(yuǎn)-平塘隱伏深斷裂以及橫穿貴州中部呈東西向的黔中深斷裂。

水庫蓄水可能伴隨地震產(chǎn)生，然而人們對水庫地震形成機(jī)制尚不清楚[1]，地震定位是地震研究中的基礎(chǔ)，因此精確定位對研究水庫地震十分重要。Waldhauser[2]于2000年應(yīng)用雙差地震定位算法對加州北海沃德斷層上的地震進(jìn)行了研究，該方法可以降低地震定位過程中對速度模型的依賴，從而提高定位精度。自該方法提出以來，被廣泛應(yīng)用于地震精確定位研究中。陳翰林等[3]利用該方法研究了龍灘水庫地震活動特征，李海鷗等[4]精確定位了三峽水庫泄灘西地區(qū)地震，吳朋等[5]對瀑布溝水庫區(qū)域地震進(jìn)行重新定位，均獲得了良好的結(jié)果。

1 理論方法

2 數(shù)據(jù)和一維速度模型

北盤江地震臺網(wǎng)包含15個固定臺站，臺站沿北盤江流域分布，如圖1所示，主要負(fù)責(zé)對光照水電站、董箐水電站水庫地震的監(jiān)測及簡評。地震臺網(wǎng)地震計(jì)均為DS-4A型短周期地震計(jì)，地震數(shù)據(jù)采集器為TDE-32CI型數(shù)據(jù)采集器。原始數(shù)據(jù)為臺站記錄到的所有地震的震相目錄，因此，需要在重定位前進(jìn)行數(shù)據(jù)篩選。我們選取出只在貴州范圍內(nèi)并且初始定位給出了震源深度的地震，由于目錄數(shù)據(jù)中有較多地震震相（PG、SG）數(shù)只有6條，為了盡可能多地利用地震，我們對震相數(shù)的限制為大于等于6條。在此基礎(chǔ)上，我們挑選出2012年2月至2012年12月的1654次發(fā)生在北盤江地震臺網(wǎng)監(jiān)測區(qū)域的地震進(jìn)行重新定位，范圍為104.4-106.2°E，25-26.6°N。

雙差定位過程中需要一維速度模型，我們在多次嘗試后，采用了鄭圻森[6]等的鄂湘黔地塊一維速度模型，如表1所示。

3 結(jié)果與分析

雙差預(yù)處理程序結(jié)果顯示，有17314個P波震相被選取，17348個S波震相被選取，并有1653個地震成功配對，得到6272個強(qiáng)連接地震對，平均連接距離為2.35km。雙差定位過程中，由于S波讀取精度較低，我們設(shè)定P波權(quán)重為1，S波權(quán)重為0.5，采用共軛梯度法迭代10次進(jìn)行求解，得到穩(wěn)定解。定位前到時平均均方根殘差為0.14s，定位后降低至0.036s。結(jié)果顯示研究區(qū)域內(nèi)有1513次地震參與重新定位，得到了1136次地震的精定位結(jié)果。參與重新定位的地震數(shù)目比成功配對的少，是由于有部分地震雖然形成地震對，但是為弱連接，在定位過程中會選取但是不參與計(jì)算。震源深度在0～25km范圍內(nèi)。重定位前后震中分布如圖2所示，重定位結(jié)果顯示，地震分布呈現(xiàn)出兩個地震叢集分布，均主要分布在河流兩側(cè)區(qū)域，地震叢集中心分別為（25.53°N，105.74°E，7.26km）和（26.03°N，105.16°E，7.78km）。

重定位前，震中分布零散，重定位后震源更集中，地震叢集更明顯。對定位前后NW方向分別做剖面，得到結(jié)果如圖3所示。定位前后震級分布如表2所示?？梢钥闯?，定位過程中有337個地震事件被剔除，每一個震級段都有事件被剔除，說明剔除條件與震級無關(guān)。我們對比圖3可以看出，重定位后邊緣的地震事件被剔除。結(jié)合圖1的臺站分布，我們認(rèn)為造成較多地震事件被剔除有兩個原因，一是由于原始數(shù)據(jù)中有較多6條震相記錄的數(shù)據(jù)，數(shù)據(jù)質(zhì)量較差，因此在定位過程中可能因?yàn)檎`差過大而被剔除；二是由于邊緣地震事件未被臺站包圍，重定位過程中誤差較大而被剔除。

從圖3可以看出，定位前地震分布為分層網(wǎng)格狀，是由于初始定位所使用的算法精度有限造成的。而使用雙差定位法進(jìn)行重定位后，定位結(jié)果得到較大改善，地震叢集分布特征更加明顯，展現(xiàn)出條帶狀分布特征。定位后震源深度在0-25km范圍內(nèi)，主要集中在4-12km范圍，83%的地震震源深度大于5km，并顯示出條帶狀分布。重定位后震源深度統(tǒng)計(jì)直方圖如圖4所示，橫坐標(biāo)為深度，縱坐標(biāo)為地震數(shù)目，地震次數(shù)隨深度變化類似正態(tài)分布，在7-8km處地震最多。

4 結(jié)束語

（1）本文利用雙差定位法對北盤江流域2012年2月至2012年12月發(fā)生的1654次地震進(jìn)行重定位，得到了1135次地震精確定位結(jié)果，精確定位后平均均方根殘差從0.14s降低至0.036s，定位結(jié)果得到較大改善。

（2）定位結(jié)果顯示，地震深度分布在0-25km范圍內(nèi)，主要分布在4-12km范圍。震源深度<5km的地震較少，大部分在5km以下，我們理解為該區(qū)域以構(gòu)造地震為主。

致謝：黎莎、趙立波、楊宜海、何富君、吳施楷為本文的寫作提供了建議和幫助咨詢，在此向他們表示衷心感謝。

參考文獻(xiàn)

[1]李碧雄，田明武，莫思特，等.水庫誘發(fā)地震研究進(jìn)展與思考[J].地震工程學(xué)報(bào)，2014，36（2）：380-386.

[2]Waldhauser F， Ellsworth W L. A double-difference earthquake location algorithm： method and application to the Northern Hayward Fault， California， from double-difference earthquake location. J. Geophys. Res.，2002，107（3）：1-15.

[3]陳翰林，趙翠萍，修濟(jì)剛，等.龍灘水庫地震精定位及活動特征研究[J].地球物理學(xué)報(bào)，2009，52（8）：20

35-2043.

[4]李海鷗，馬文濤，徐錫偉，等.三峽水庫泄灘西地區(qū)地震雙差定位及成因分析[J].地震地質(zhì)，2010，32（4）：56

4-569.

[5]吳朋，蘇金蓉，陳天長，等.瀑布溝水庫區(qū)域地震的重新定位[J].四川地震，2014（1）：1-4.

[6]鄭圻森，朱介壽，曹家敏，等.華南地區(qū)巖石圈地殼速度結(jié)構(gòu)數(shù)據(jù)處理[J].物探化探計(jì)算技術(shù)，2004，26（2）：97-100.