李 俊，陳智偉，李金龍，鄒振軒，于俊誼，俞鐵宏

(1.浙江省地震局，浙江 杭州 310013；2.溫州地震臺，浙江 溫州 325015）

0 前言

隨著國家 “九五”、 “十五”和 “十一五”計劃期間，浙江省完成了幾十個地震臺站數(shù)字化建設(shè)，這為地震工作和地震學(xué)研究提供了良好的基礎(chǔ)和數(shù)據(jù)。地震數(shù)字化和地震臺站的加密建設(shè)對中小地震監(jiān)測能力有著很大的提高，也大大的豐富了地震事件的記錄。

巖石物理實驗研究和數(shù)值模擬實驗都表明，斷裂帶在地震前后經(jīng)歷了較大的應(yīng)力變化，會造成斷裂帶巖石中微裂隙的增加[1～4]。而地震學(xué)研究中對波速變化測量對地震預(yù)報預(yù)測工作也有有著重要的意義，而波速變化的測量是往往通過重復(fù)測量相同路徑上的地震波走時來實現(xiàn)的，只有保證震源的重復(fù)性，才能進(jìn)行有效的波速變化測量。除了利用人工震源之外，近年的研究表明重復(fù)地震也可為我們進(jìn)行波速變化測量提供了一個途徑，因此近年來國際國內(nèi)對重復(fù)地震的研究也越來越廣泛。

1 研究現(xiàn)狀

“重復(fù)地震”很早就已提出[5]，近年來開始引起越來越多的關(guān)注，被廣泛應(yīng)用于地下介質(zhì)物性變化測量[6]、地震預(yù)測[7]、斷層結(jié)構(gòu)和震源物理研究等[8、9]。 “重復(fù)地震”對于地震臺網(wǎng)定位能力的評估和改善具有重要的實際意義[10～13]，不僅可用來估計臺網(wǎng)的定位精度，而且在地震"復(fù)發(fā)"和滑動速率估計方面[14～16]，在檢測地下介質(zhì)性質(zhì)的時間變化方面[17～22]都顯示出值得關(guān)注的應(yīng)用潛力。尤其值得注意的是，近年來的一個重要發(fā)現(xiàn)是， “重復(fù)地震”在地震活動中看來并不是罕見的特例。例如，一個有代表性的工作是Schaff和Richards[23]給出的中國大陸及周邊地區(qū)1987～2000年的1301次地震組成的950對 “重復(fù)地震對”，這些地震占同期同地全部地震總數(shù)的10%，而國內(nèi)有人近年的研究結(jié)果則是更多。

重復(fù)地震在地震斷層附近多有發(fā)生，一般具有相似的震級和極為相似的波形和震源機(jī)制[24、25]。多數(shù)研究人員認(rèn)為，重復(fù)地震的發(fā)生同介質(zhì)穩(wěn)滑區(qū)域內(nèi)的凹凸體有關(guān)[26]。在美國的Parkfield斷層[27]、Calaveras斷層[28]、Hayward斷層[29]、 土耳其的North Aantolian斷層[30]以及日本的的俯沖帶[31、32]都發(fā)現(xiàn)有重復(fù)地震。中國唐山震區(qū)、云南小江斷裂帶以及臺灣東部地區(qū)也都發(fā)現(xiàn)大量重復(fù)地震[33～36]。Schaff and Richards[25]的研究表明，中國境內(nèi)發(fā)生的地震約有10%為重復(fù)地震，Schaff等對重復(fù)地震的定義即一對地震事件被至少一個臺站記錄到的波形的相關(guān)系數(shù)不小于0.8。本文與之相關(guān)研究中的重復(fù)地震均采用這一定義，并稱其為波形相似意義上的重復(fù)地震。重復(fù)地震被廣泛應(yīng)用于地下介質(zhì)物性變化測量,地震預(yù)測和臺網(wǎng)定位精度評估[37～39]。在斷層行為和震源物理研究中，重復(fù)地震也得到了廣泛的應(yīng)用。

2 地震數(shù)據(jù)

珊溪水庫于2000年5月下閘蓄水，2002年7月28日庫區(qū)發(fā)生3.5級 （為ML震級，下同）地震后，每年都有地震活動，到2012年止，庫區(qū)記錄到ML≥1.0級地震共1 300多次，其中2.0～2.9級243次，3.0～3.9級46次，4.0級以上13次，最大為2006年2月9日4.6級。根據(jù)浙江省數(shù)字地震臺網(wǎng)測定，地震全部發(fā)生在大壩上游的文成縣珊溪鎮(zhèn)與泰順縣包垟鄉(xiāng)交界處水庫庫區(qū)，震中集中分布在一個由北西走向和北東走向斷裂圍限的小菱形塊體（約 7×9 km）內(nèi)， 震源深度為 1～10 km(圖 1)[40～43]。

當(dāng)兩個地震事件在同一臺站顯示幾乎相同的地震波形圖時，便構(gòu)成重復(fù)地震，兩者不僅應(yīng)具有相同的震源位置，還要具有一樣的震源機(jī)制，本文認(rèn)為，如果在同一地點發(fā)生，震源深度相當(dāng)，波形又相似，則可以認(rèn)為震源機(jī)制相同或者相近，因此本文只從波形的相似度上來考慮地震的相似性和重復(fù)性。

3 地震對整理

本文對珊溪水庫的地震事件進(jìn)行截取、挑選、震相標(biāo)記和濾波處理后，再對事件波形進(jìn)行兩兩互相關(guān)處理。波形互相關(guān)方法[24]可用來檢測地震波形之間的相似性，Schaf等也因此將在P波之前5 s至 波之后40 s的時間窗內(nèi)具有0.8以上相關(guān)系數(shù)的兩地震事件定義為重復(fù)地震。

互相關(guān)計算公式：

上式中，f1(t)和f2(t)分別為重復(fù)地震對的兩個不同的波形時間序列，τ為互相關(guān)函數(shù)計算中的時間延遲，C（τ）為計算的互相關(guān)函數(shù)。

本文選擇了記錄歷史較長的HT(黃坦)臺和WEZ(溫州)臺的記錄進(jìn)行處理和重復(fù)地震挑選，最終得到這兩個臺站記錄的重復(fù)地震對。

由于WEZ臺距離震中區(qū)域200多公里，因此能清晰記錄的地震事件波形對震級要求相對高一些，對WEZ臺的地震數(shù)據(jù)進(jìn)行篩選，得到了55個ML3.0級左右的地震事件，對這些事件進(jìn)行處理和互相關(guān)分析后，得到了圖3的結(jié)果。選取相關(guān)系數(shù)為0.8以上的地震事件對，則55個地震中有40個，分為13組相似地震對，如果設(shè)相關(guān)系數(shù)閥值為0.9時，則有21個，分為8組相似地震對。

圖1 珊溪水庫地震震中分布圖Fig.1 Epicenter distribution of Shanxi reservoir earthquakes

圖2 兩地震事件波形互相關(guān)方法示意圖Fig.2 Diagram of waveform cross-correlation method for two earthquake events

圖3 WEZ(溫州)臺重復(fù)地震簇結(jié)果Fig.3 Clusters of repeated earthquakes recorded by WEZ station

黃坦臺建臺較早，在本研究中，我們挑選了黃坦臺數(shù)據(jù)質(zhì)量好，震級在ML2.0以上的223個地震事件，經(jīng)過Pg震相標(biāo)定，濾波等常規(guī)處理后，對223個地震事件進(jìn)行了兩兩互相關(guān)計算，然后根據(jù)聚類算法，挑選得到重復(fù)地震對。經(jīng)計算，其中設(shè)定相關(guān)系數(shù)閥值為0.85時，有湊對的相似地震144個，占總數(shù)223個的64.5%；對比相似地震對的原始地震波形發(fā)現(xiàn)，基本上相似地震對的垂直向P波初動方向也是一致的。圖4(a)則為相關(guān)系數(shù)閥值為0.9的重復(fù)地震對，共有106個，占47.5%，圖4(b)為其中典型的地震對波形比較；說明珊溪水庫地震群的相似度比較高，地震重復(fù)率比較高。

圖5為不同時期的地震對比較，其中圖5(a)為2010年一組相似地震，5(b)為2006年和2007年的相似地震，結(jié)合圖4，可以看出在不同時期發(fā)生的相似地震組，其波形不一樣，這說明了這些地震組的發(fā)震機(jī)制或者區(qū)域介質(zhì)有著較大的變化，可根據(jù)震源機(jī)制解和水庫庫區(qū)的水位等數(shù)據(jù)作進(jìn)一步的分析和推斷。

4 分析與討論

圖4 黃坦臺地震對結(jié)果（a)和典型地震對波形圖(b)Fig.4 Clusters of repeated earthquakes recorded by HT station（a)and seismic waveforms of typical repeated earthquakes(b)

圖5 HT臺相似地震對比較Fig.5 Comparison of similar earthquakes recorded by HT station

從本次珊溪水庫中小地震事件的重復(fù) (相似)地震的處理結(jié)果來看，發(fā)震時間相近的地震相似度較高，時間尺度相距較長則相似率較低。

重復(fù)地震的研究主要為了進(jìn)一步得到介質(zhì)的狀態(tài)變化，主要為波速的變換，但是由于記錄時間不夠長，地震數(shù)目也比較少，特別是2008年后發(fā)生的數(shù)量較少，因此暫時得不到較長時間段的介質(zhì)變化，有待更長周期的研究。另外，在研究時間范圍內(nèi)，臺站有改造搬遷，其中黃坦 (HT)臺曾經(jīng)在2007年移過位置，與原先臺站距離1 km左右，這對結(jié)果也會造成一定的影響。WEZ(溫州)臺由于臺站建立時間較長，隨著儀器老化信噪比降低，監(jiān)測能力有所下降，導(dǎo)致不能清晰記錄2009年后的2～3級地震波形。

本文的研究雖然結(jié)合了相似地震對的P波初動來進(jìn)一步判斷重復(fù)性，但是如果今后能結(jié)合事件震源機(jī)制解的研究，對判斷和研究重復(fù)地震將更有意義，結(jié)合尾波干涉的方法更可較精確的測定波速的變化，從而對地震預(yù)報預(yù)測工作提供更大的幫助。

致謝：感謝中國地震局地球物理研究所王偉濤同志的幫助與指導(dǎo)。

[1]Aki K.Asperities,barriers,characteristic earthquakes and strong motion prediction[J].Geophys.Res.,1984,89:5867-5872.

[2]J.R.Rice.Fault stress states,pore pressure distributions,and the weakness of the san andreas fault[A].In B.Evans and T.-F.Wong,editors,Fault Mechanics and Transport Properties of Rocks[C].Academic,San Diego,Calif.,1992,475-503

[3]James H.Dieterich.Time-dependent friction and the mechanics of stick-slip[J].Pure and Applied Geophysics,1978,116(4):790--806.

[4]C.H.Scholz.The Mechanics of Earthquakes and Faulting[M].Cambridge Univ.Press,New York,2002.

[5]Isacks B L， Sykes I R， Oliver J.Spatial and temporal clustering of deep and shallow earthquakes inthe Fiji-Tonga kermadec region[J].Bull Seism Soc Amer， 1967， 57(5)： 935-958.

[6]Poupinet G， Ellsworth W L， Frechet J.Monitoring velocity variations in the crust using earthquakedoublets：an application to the Calaveras Fault， California[J].J Geophys Res， 1984， 89(B7)： 5719-5731.

[7]Tthatcher W.The earthquake deformation cycle,recurrence and the time-predictable model[J].Geophys.Res.,2004,89(B7):5674-5680.

[8]R.M.Nadeau,W.Foxall,and T.V.McEvilly.Clustering and periodic recurrence of microearthquakes on the san andreas fault at parkfield,california[J].Science,1995,267(5197):503--507.

[9]Robert M.Nadeau and Lane R.Johnson.Seismological studies at parkfield vi:Moment release rates and estimates of source parameters for small repeating earthquakes[J].Bulletin of the Seismological Society of America,1998,88(3):790--814.

[10]Allan M.Rubin.and Gillard,D.2000.Aftershock asymmetry/rupture directivity among central San Andreas fault microearthquakes[J].Geophys.Res.,2000,105:19095-19109.

[11]蔣長勝、吳忠良，由 “重復(fù)地震”給出的中國地震臺網(wǎng)的定位精度估計 [J].中國地震，2005，,21(2)：147-154。

[12]David P.Schaff and Gregory C.Beroza.Coseismic and postseismic velocity changes measured by repeating earthquakes[J].Geophys.Res.,2004,109:B10302.

[13]David P.Schaff and Paul G.Richards.Repeating seismic events in china[J].Science,2004,303(5661):1176-1178.

[14]Li Yong-Gang and John E.Vidale.Healing of the shallow fault zone from 1994-1998 after the 1992 m7.5 landers,california,earthquake[J].Geophys.Res.Lett.,2001,28(5):2999-3002.

[15]Snieder,R.,and M.Vrijlandt,Constraining Relative Source Locations with Coda Wave Interferometry:Theory and Application to Earthquake Doublets in the Hayward Fault,California[J].Geophys.Res.,2005,110;B04301.

[16]Toshihiro Igarashi,Toru Matsuzawa,and Akira Hasegawa.Repeating earthquakes and interplate aseismic slip in the northeastern japan subduction zone[J].Geophys.Res.,2003,108:2249.

[17]Fenglin Niu,Paul G.Silver,Thomas M.Daley,Xin Cheng,and Ernest L.Majer.Preseismic velocity changes observed from active source monitoring at the parkfield safod drill site[J].Nature,2008,454(7201):204--208.

[18]林建民，王寶善，葛洪魁，等.重復(fù)地震及其在人工探測中的潛在應(yīng)用[J].中國地震，2006，22(1)：1-9.

[19]武敏捷，徐 平，林向東，等.利用重復(fù)地震觀測灤縣地區(qū)的地殼介質(zhì)變化 [J].華北地震科學(xué)，2011， 29(3)： 23-27.

[20]武敏捷，徐平，馬士振，等.利用重復(fù)地震研究灤縣地區(qū)的地殼介質(zhì)變化[J].地震地磁觀測與研究，2011， 32(5)： 20-24.

[21]李宇彤，蔣長勝.波形相關(guān)意義 “重復(fù)地震”研究綜述[J].中國地震，2011，27(4)：335-347.

[22]王偉濤，王寶善.重復(fù)地震研究綜述[J].西北地震學(xué)報，2011，33(3)：305-311，F(xiàn)0003.

[23]Schaff D P， Richards P G.Repeating seismic events in China[J].Science， 2004， 303： 1176-1178.

[24]G.Poupinet,W.L.Ellsworth,and J.Frechet.Monitoring velocity variations in the crust using earthquake doublets:An application to the calaveras fault,california[J].Geophys.Res.,1984,89:6719-5731.

[25]David P.Schaff and Gregory C.Beroza.Coseismic and postseismic velocity changes measured by repeating earthquakes[J].Geophys.Res.,2004,109:B10302.

[26]Tracy Johnson.Time-dependent friction of granite:implications for precursory slip on faults[J].Geophys.Res.,1981,86:6017-6028.

[27]W.L Ellsworth.Characteristic earthquakes and long-term earthquake forecasts,implications of central california seismicity[A].In F.Y.Cheng and M.S.Sheu,editors,Urban Disaster Mitigation:the Role of Science and Technology[C].Elsevier Science Ltd,1995，1-14.

[28]Dennise C.Templeton,Robert M.Nadeau,and Roland Burgmann.Behavior of repeating earthquake sequences in central california and the implications for subsurface fault creep[J].Bulletin of the Seismological Society of America,2008,98(1):52-65.

[29]Roland Burgmann,D.Schmidt,R.M.Nadeau,M.d'Alessio,E.Fielding,D.Manaker,T.V.McEvilly,and M.H.Murray.Earthquake potential along the northern hayward fault,california[J].Science,2000,289(5482):1178-1182.

[30]Peng Zhigang,John E.Vidale,Chris Marone,and Allan Rubin.Systematic variations in recurrence interval and moment of repeating aftershocks[J].Geophys.Res.Lett.,2005,32:L1531.

[31]Toshihiro Igarashi,Toru Matsuzawa,and Akira Hasegawa.Repeating earthquakes and interplate aseismic slip in the northeastern japan subduction zone[J].Geophys.Res.,2003,108:2249.

[32]Naoki Uchida,Toru Matsuzawa,Akira Hasegawa,and Toshihiro Igarashi.Recurrence intervals of characteristic m4.81 earthquakes off-kamaishi,ne japan--comparison with creep rate estimated from small repeating earthquake data[J].Earth and Planetary Science Letters， 2005， 233(1-2):155--165.

[33]李 樂，陳棋福.利用重復(fù)地震估算北京平原地區(qū)隱伏斷裂深部的活動速率 [J].地震地質(zhì)，2010，32(3)：508-519.

[34]王偉濤，王寶善.基于聚類分析的多尺度相似地震快速識別方法及其在汶川地震東北端余震序列分析中的應(yīng)用[J].地球物理學(xué)報，2012，55(6)：1952-1962.

[35]Li,L.,Q.-F.Chen,X.Cheng,and F.Niu,Spatial clustering and repeating of seismic events observed along the 1976 Tangshan fault,north China[J].Geophys.Res.Lett.,2007,34.

[36]李 樂 陳棋福.利用重復(fù)地震估算北京平原地區(qū)隱伏斷裂深部的活動速率 [J].地震地質(zhì)，2010，32(3)：508-519.

[37]李宇彤，吳忠良，蔣長勝，等.利用遼寧區(qū)域地震臺網(wǎng)記錄分析 “重復(fù)地震”[J].地震學(xué)報，2008，30(4)： 383-396.

[38]張致偉，程萬正，吳 朋，等.利用重復(fù)地震P波走時差觀測地殼介質(zhì)的微動態(tài) [J].中國地震，2010， 28(4)： 357-372.

[39]蔣長勝，吳忠良，李宇彤.首都圈地區(qū) “重復(fù)地震事件”及其在區(qū)域地震臺網(wǎng)定位精度評價中的應(yīng)用.[J].地球物理學(xué)報， 2008， 51(3)： 817-827

[40]鐘羽云，張帆，趙冬.珊溪水庫ML4.6震群精確定位與發(fā)震構(gòu)造研究 [J].地震研究，2011，34(2)：158-165.

[41]鐘羽云，張震峰，闞寶祥.溫州珊溪水庫地震重新定位與速度結(jié)構(gòu)聯(lián)合反演[J].中國地震，2010(3)：265-272.

[42]鐘羽云，周昕，張帆，等.2006年溫州珊溪水庫地震序列特征[J].華南地震，2007，27(1)：21-30.

[43]李金龍，鄒振軒，張震峰.珊溪水庫震群的地震定位[J].地震地磁觀測與研究，2007，28(1)：14-18.