福建水口及仙游地震序列的譜振幅相關(guān)系數(shù)變化特征

袁麗文，李強(qiáng)，陳智勇

（福建省地震局，福州350003）

福建水口及仙游地震序列的譜振幅相關(guān)系數(shù)變化特征

袁麗文，李強(qiáng)，陳智勇

（福建省地震局，福州350003）

利用福建省數(shù)字地震臺(tái)網(wǎng)記錄的波形數(shù)據(jù)，計(jì)算了水口和仙游地震序列的譜振幅相關(guān)系數(shù)，結(jié)果顯示水口序列譜振幅相關(guān)系數(shù)在0.90-0.98之間變化，相關(guān)系數(shù)高，序列地震震源機(jī)制解相似程度高；仙游序列譜振幅相關(guān)系數(shù)在0.72-0.81之間變化，相關(guān)系數(shù)略低，震源機(jī)制解整體相似程度不高。結(jié)合福建地區(qū)同期地震活動(dòng)背景水平分析認(rèn)為，較高的譜振幅相關(guān)系數(shù)與福建地區(qū)較高的背景應(yīng)力積累水平存在一定相關(guān)性。水口和仙游序列譜振幅相關(guān)系數(shù)時(shí)序變化顯示，序列最大地震發(fā)生前的譜振幅相關(guān)系數(shù)存在上升變化，顯示較大地震發(fā)生前震區(qū)應(yīng)力水平有增強(qiáng)跡象。

地震序列；譜振幅相關(guān)系數(shù)；震源機(jī)制一致性；應(yīng)力水平

0 引言

基于多種孕震物理模式的研究發(fā)現(xiàn)，在大地震孕育過(guò)程中，震前在震源區(qū)及鄰近地區(qū)應(yīng)力水平有增強(qiáng)過(guò)程。在區(qū)域應(yīng)力水平增強(qiáng)的背景下，震源區(qū)介質(zhì)的微裂隙呈現(xiàn)某一優(yōu)勢(shì)方向上的定向排列，介質(zhì)各向異性增強(qiáng)，強(qiáng)震前的前震序列震源機(jī)制解應(yīng)趨于一致。我國(guó)研究學(xué)者對(duì)國(guó)內(nèi)外多次強(qiáng)震前小震的震源機(jī)制解進(jìn)行了研究，顯示出主震前中小地震震源機(jī)制解趨于一致的現(xiàn)象明顯[1-5]。因此，利用小震震源機(jī)制解一致性判定小震群后續(xù)地震活動(dòng)趨勢(shì)，，進(jìn)而判定小震群序列類型是可行的。然而，受觀測(cè)臺(tái)網(wǎng)布局和臺(tái)站密度等的觀測(cè)條件的限制，有時(shí)直接求解小震群中大量中小地震的震源機(jī)制解存在一定困難，其求解精度也難以保證。因此，Lund等[6]提出了微震體波譜振幅相關(guān)分析方法，將求解大量中小地震震源機(jī)制解是否相似的問(wèn)題轉(zhuǎn)化為求解其零頻譜振幅是否相關(guān)的。

國(guó)內(nèi)許多研究人員利用譜振幅相關(guān)系數(shù)方法開(kāi)展了一系列地震序列的研究工作，獲得了譜振幅相關(guān)系數(shù)與序列地震震源機(jī)制一致性、區(qū)域應(yīng)力水平和地震活動(dòng)趨勢(shì)判定等的研究成果。崔子健[7]等對(duì)云南西北部5組小震群序列開(kāi)展了震例的回溯性檢驗(yàn)工作，通過(guò)譜振幅相關(guān)系數(shù)的研究顯示，前兆震群的譜振幅相關(guān)系數(shù)較高，一般性的非前兆震群活動(dòng)其譜振幅相關(guān)系數(shù)較低，提出可利用譜振幅相關(guān)系數(shù)的大小判定小震群的前兆性質(zhì)，從而進(jìn)行區(qū)域地震危險(xiǎn)性預(yù)測(cè)研究。黃浩等[8]對(duì)滇西地區(qū)6組地震序列（含4種不同序列類型）進(jìn)行譜振幅相關(guān)系數(shù)研究，分析認(rèn)為其較高的譜振幅相關(guān)系數(shù)與區(qū)域較高的應(yīng)力積累水平有關(guān)，利于后續(xù)中強(qiáng)地震的發(fā)生，但不能把譜振幅相關(guān)系數(shù)較高作為判定前震序列的充分條件。王瑩等[9]對(duì)陜西及周邊地區(qū)4個(gè)中小地震叢集序列進(jìn)行了研究，指出譜振幅相關(guān)系數(shù)較高的地震序列所在區(qū)域，其后續(xù)強(qiáng)震危險(xiǎn)性較大。除了上述利用強(qiáng)震發(fā)生前的小震序列譜振幅相關(guān)系數(shù)進(jìn)行中強(qiáng)地震活動(dòng)趨勢(shì)預(yù)測(cè)之外，也有研究通過(guò)對(duì)主震后余震序列譜振幅相關(guān)系數(shù)的時(shí)序分析，指出余震活動(dòng)時(shí)段的譜振幅相關(guān)系數(shù)上升或持續(xù)高值可能與強(qiáng)余震活動(dòng)或震源區(qū)周邊中強(qiáng)地震活動(dòng)有關(guān)[10-11]，若余震序列的譜振幅相關(guān)系數(shù)下降，則后續(xù)余震活動(dòng)發(fā)生更大地震可能性不大[12]。張致偉等[13]將譜振幅相關(guān)系數(shù)和地震視應(yīng)力研究結(jié)合起來(lái)，對(duì)四川自貢和長(zhǎng)寧兩個(gè)注水誘發(fā)的地震序列進(jìn)行了研究，指出譜振幅相關(guān)系數(shù)的高值不一定對(duì)應(yīng)區(qū)域處于高應(yīng)力水平，也可能是區(qū)域應(yīng)力場(chǎng)受到流體作用影響。此外，亦有研究人員利用振幅相關(guān)系數(shù)的聚類分析法研究序列地震震源機(jī)制變化特征，由震源機(jī)制一致性增強(qiáng)或減弱的變化來(lái)判定區(qū)域未來(lái)發(fā)生中強(qiáng)地震的危險(xiǎn)性[14-16]。

本文利用譜振幅相關(guān)系數(shù)方法，對(duì)近年來(lái)福建省最顯著的水口和仙游地震序列活動(dòng)進(jìn)行分析研究，類比已有震例的研究結(jié)果，分析兩個(gè)序列譜振幅相關(guān)系數(shù)變化與震源機(jī)制一致性、序列顯著地震活動(dòng)和區(qū)域應(yīng)力水平的相關(guān)性。

1 原理及方法

地震記錄位移譜包含了震源效應(yīng)、地震與臺(tái)站之間的傳播路徑效應(yīng)、臺(tái)站場(chǎng)地效應(yīng)、臺(tái)站附近的地面運(yùn)動(dòng)噪聲、地震儀的儀器響應(yīng)和臺(tái)站附近地表自由表面效應(yīng)等信息，可將其表達(dá)為：

其中，Uij（f）表示地震i在臺(tái)站j記錄到的位移譜，f為頻率；i為地震；j為臺(tái)站；Si（f）為地震i包含輻射圖型因子的震源譜；Pij（f）為地震i與臺(tái)站j之間的傳播路徑效應(yīng)；Nj（f）為臺(tái)站j附近的地面運(yùn)動(dòng)噪聲；L＇j（f）為臺(tái)站j處的局地場(chǎng)地效應(yīng)；Ij（f）為臺(tái)站j的儀器響應(yīng)；Surj為臺(tái)站j附近地表的自由表面效應(yīng)。

Si（f）為地震i包含輻射圖型因子的震源譜，可表達(dá)為：

式（2）中，S＇i（f）為地震i的震源譜；φij為地震i的震源輻射圖型因子。震源輻射圖型因子包含了震源處地震的震源機(jī)制信息（斷層的走向、傾角、滑動(dòng)角）和震源位置信息（震源到臺(tái)站的方位角、離源角）。

根據(jù)Brune模型[17]，地震i的震源譜可表達(dá)為：

式（3）中，Ω＇0i為地震i的零頻極限值；fc為拐角頻率。

將（3）式帶入（2）式中有：

式（4）中，Ω＇0iφij表示臺(tái)站j記錄到的地震i的包含輻射圖型因子的震源譜零頻極限值，即譜振幅Ω0ij。

地震i與臺(tái)站j之間的傳播路徑效應(yīng)可表示為：

式（5）可見(jiàn)，地震與臺(tái)站間的傳播路徑效應(yīng)受到幾何擴(kuò)散因子Gij、震源距Rij、介質(zhì)品質(zhì)因子Q（f）和地震波波速v的影響。當(dāng)兩次地震被同一個(gè)臺(tái)站記錄到，且兩次地震之間的距離足夠接近時(shí)，其傳播路徑效應(yīng)幾乎相等。

當(dāng)?shù)卣饃和地震y的距離足夠接近時(shí)，地震射線到同一臺(tái)站的傳播路徑幾乎一樣，即Pxj（f）與Pyj（f）幾乎相等，根據(jù)公式（4）可見(jiàn)，臺(tái)站j記錄到的地震x和地震y的地面運(yùn)動(dòng)位移譜Uxj（f）、Uyj（f）與Ω＇0x、Ω＇0y和φxj、φyj有關(guān)。當(dāng)震源機(jī)制相同的情況下，即φxj=φyj時(shí)，地面運(yùn)動(dòng)位移譜只與Ω＇0x、Ω＇0y的相對(duì)大小有關(guān)；在震源機(jī)制不同的時(shí)，則還與表征震源機(jī)制信息的φxj、φyj有關(guān)。因此，通過(guò)地震記錄反演譜振幅時(shí)，若震源機(jī)制相同，反演得到的譜振幅Ω0xj和Ω0yj的相對(duì)大小只與Ω＇0x和Ω＇0y的相對(duì)大小有關(guān)；否則，還與震源機(jī)制有關(guān)。因此，可由臺(tái)站記錄反演地震x與y的譜振幅，通過(guò)計(jì)算相關(guān)系數(shù)rxy來(lái)描述兩次地震震源機(jī)制是否相似的問(wèn)題。

地震x與y的譜振幅相關(guān)系數(shù)rxy可表示為[6]：

式（9）中：xjl和yjl分別表示臺(tái)站j記錄的地震x和y的l分量波形資料反演所得到的譜振幅Ω0xj和Ω0yj的對(duì)數(shù)值；l為地震波形的分量，包含P波垂向、徑向和S波垂向、徑向和切向共5個(gè)分量，對(duì)地震波形的5分量分別計(jì)算譜振幅，而后進(jìn)行互相關(guān)計(jì)算，可以充分利用地震波形所包含的震源機(jī)制信息，提高研究結(jié)果信度。

對(duì)一組地震序列按地震發(fā)生時(shí)間順序進(jìn)行排列，序列第m（m≥2）次地震與前面m-1次地震編為一組，計(jì)算組內(nèi)兩兩地震之間的譜振幅相關(guān)系數(shù)rxy，將所有兩兩地震間的相關(guān)系數(shù)求算數(shù)平均值，則得到第m個(gè)地震時(shí)的譜振幅相關(guān)系數(shù)。以滑動(dòng)步長(zhǎng)為1進(jìn)行滑動(dòng)，依次完成每個(gè)地震時(shí)刻譜振幅相關(guān)系數(shù)的計(jì)算，則可得到隨時(shí)間變化的譜振幅相關(guān)系數(shù)。

根據(jù)上述原理分析，對(duì)水口和仙游地震序列采取如下計(jì)算步驟：

（1）收集水口和仙游序列地震波形資料。

（2）將南北、東西和垂直向的三分量地震記錄旋轉(zhuǎn)到徑向、切向和垂直向上，采用延遲時(shí)間窗方法[18-19]計(jì)算地震位移譜。圖1為本文使用的一個(gè)地震的旋轉(zhuǎn)后的波形及其位移譜圖像。

（3）分別利用水口和仙游周邊地區(qū)背景地震資料（圖2），采用多臺(tái)聯(lián)合反演方法[20]反演水口和仙游地區(qū)的Q值，得到水口地區(qū)QP=34.8f0.906，QS=83.5f0.829；仙游地區(qū)QP=47.7f0.843，QS= 210.5f0.631。

（4）采用Moya等[21]引入的遺傳算法反演各臺(tái)站的場(chǎng)地效應(yīng)。

（5）獲得經(jīng)過(guò)傳播路徑、場(chǎng)地響應(yīng)、噪聲和儀器響應(yīng)校正后的譜振幅，并計(jì)算譜振幅相關(guān)系數(shù)。

2 計(jì)算結(jié)果

2.1 2007至2008年水口地震序列

福建水口地區(qū)自1993年7月開(kāi)始小震活動(dòng)，之后地震活動(dòng)呈斷續(xù)起伏狀態(tài)，但震級(jí)水平不高，但自2007年10月開(kāi)始，水口地區(qū)地震活動(dòng)頻次和強(qiáng)度水平顯著上升，截止2008年12月共發(fā)生ML2.5級(jí)以上地震45次，其中ML4.0級(jí)以上地震3次，最大地震為2008年3月6日福建水口ML4.8級(jí)地震（26.37°N，118.67°E）。本研究利用南平臺(tái)（NAP）、福州臺(tái)（FUZ）、閩清臺(tái)（MIQ）、古田臺(tái)（GUT）、尤溪臺(tái)（YOX）、寧德臺(tái)（NID）、德化湯頭臺(tái)（DHTT）等7個(gè)地震臺(tái)站記錄到的33次ML2.5級(jí)以上地震（圖3）的波形資料計(jì)算譜振幅相關(guān)系數(shù)。水口序列地震的譜振幅相關(guān)系數(shù)在0.90～0.98之間變化（圖4），相關(guān)系數(shù)數(shù)值高，表明序列小震震源機(jī)制解一致性好。

圖1 永定臺(tái)（YOD）記錄到的一次地震的波形（a）及計(jì)算所得的相應(yīng)觀測(cè)位移譜（b）Fig.1 Waveform of one earthquake recorded by YOD station（a），and the calculated spectrum of the observed displacement（b）

圖2 地震與臺(tái)站射線分布圖Fig.2 Distribution of earthquakes and stations

2.2 2012至2014年仙游地震序列

福建仙游地區(qū)自2010年10月開(kāi)始小震群活動(dòng)，隨后地震活動(dòng)斷續(xù)起伏，截止2014年12月序列仍有ML2.0級(jí)以上地震活動(dòng)。序列共發(fā)生ML4.0級(jí)以上地震5次，其中最大地震為2013年9月3日福建仙游ML5.0（MS4.8）級(jí)地震（25.63°N， 118.75°E）。本研究利用安溪長(zhǎng)坑臺(tái)（AXCK）、安溪參內(nèi)臺(tái)（AXCN）、德化唐寨臺(tái)（DHTZ）、福清東張臺(tái)（FQDZ）、福州城門(mén)臺(tái)（FZCM）、福州瑯岐臺(tái)（FZLQ）、福州日溪臺(tái)（FZRX）、惠安凈峰臺(tái)（HAJF）、閩侯竹岐臺(tái)（MHZQ）、閩清東橋（MQDQ）、莆田埭頭（PTDT）、莆田荔城臺(tái)（PTLC）、莆田湄洲臺(tái)（PTMZ）、莆田南日臺(tái)（PTNR）、泉港前黃臺(tái)（QGQH）、泉州臺(tái)（QZH）、仙游西苑臺(tái)（XYXY）、永春石鼓臺(tái)（YCSG）、永泰葛嶺臺(tái)（YTGL）、永泰蓋洋臺(tái)（YTGY）、尤溪坂面臺(tái)（YXBM）等21個(gè)地震臺(tái)記錄到的52次ML2.5級(jí)以上地震（圖5）的波形資料計(jì)算譜振幅相關(guān)系數(shù)。仙游序列地震譜振幅相關(guān)系數(shù)在0.72～0.81之間變化（圖6），相關(guān)系數(shù)略低，說(shuō)明序列小震震源機(jī)制一致性不高。

圖3 2007—2008年水口地震序列及周邊臺(tái)站分布Fig.3 Earthquakes and stations around Shuikou area between 2007 and 2008

圖4 水口地震序列譜振幅相關(guān)系數(shù)和M-t圖Fig.4 The correlation coefficient of spectral amplitude and M-t plot of Shuikou earthquake sequence

3 討論與結(jié)論

3.1 譜振幅相關(guān)系數(shù)與序列地震參數(shù)

譜振幅相關(guān)系數(shù)從震源物理角度提供了對(duì)震源機(jī)制解一致性、區(qū)域應(yīng)力場(chǎng)狀態(tài)的描述，為判定地震序列活動(dòng)趨勢(shì)提供了重要信息。在我國(guó)的地震分析預(yù)報(bào)實(shí)踐中，許多學(xué)者利用地震目錄，通過(guò)多大量震例的分析研究，總結(jié)出了利用表征震群能量釋放均勻度U值、震群歸一化信息熵K值、序列地震時(shí)間叢集程度ρ值、序列頻次衰減快慢程度h值和序列大小地震比例關(guān)系b值等對(duì)序列活動(dòng)特征及趨勢(shì)進(jìn)行判定，并提出了一套關(guān)于前兆震群判定的指標(biāo)體系[22-24]。當(dāng)其中任一參數(shù)符合前兆震群閾值范圍，則判定震群為前兆震群，通常多指標(biāo)判定效果優(yōu)于單一指標(biāo)。鑒于水口和仙游地震序列持續(xù)活動(dòng)時(shí)間較長(zhǎng)，序列地震活動(dòng)活躍-平靜交替現(xiàn)象明顯，根據(jù)序列參數(shù)定義，對(duì)兩個(gè)序列進(jìn)行分時(shí)段序列參數(shù)的計(jì)算。U值、K值、ρ值、b值和h值計(jì)算結(jié)果見(jiàn)表1。

圖5 2012—2014年仙游地震序列及周邊臺(tái)站分布Fig.5 Earthquakes and stations around Xianyou area between 2012 and 2014

圖6 仙游地震序列譜振幅相關(guān)系數(shù)和M-t圖Fig.6The correlation coefficient of spectral amplitude and M-t plot of Xianyou earthquake sequence

表1 水口、仙游地震序列測(cè)震學(xué)參數(shù)統(tǒng)計(jì)表Table 1 The seismometry parameters of Shuikou and Xianyou earthquake sequence

從表1中看出，編號(hào)1的水口序列h值，ρ值符合前震序列判斷閾值，該序列具有前兆性質(zhì)，且該時(shí)段內(nèi)譜振幅相關(guān)系數(shù)高，尤其在2007年10月為研究時(shí)段的最高值，顯示這一時(shí)期區(qū)域應(yīng)力場(chǎng)增強(qiáng)，小震震源機(jī)制解一致性程度增強(qiáng)，隨后發(fā)生2008年3月6日水口ML4.8級(jí)地震（序列最大地震），該地震后水口序列持續(xù)活動(dòng)；編號(hào)2地震序列參數(shù)計(jì)算結(jié)果僅ρ值還存在一定的前兆性質(zhì)，同期譜振幅相關(guān)系數(shù)雖然數(shù)值高，但整體呈下降趨勢(shì)（圖4），隨后水口序列雖然仍持續(xù)了長(zhǎng)時(shí)間的活動(dòng)，但整體地震活動(dòng)水平呈衰減狀態(tài)。編號(hào)3的仙游序列除K值不符合前震序列判斷閾值外，其余參數(shù)均表現(xiàn)出明顯的前兆性質(zhì)，該時(shí)段發(fā)生的最大地震為2012年4月15日仙游ML4.1級(jí)地震，譜振幅相關(guān)系數(shù)為研究時(shí)段內(nèi)最高，且ML4.1級(jí)地震前相關(guān)系數(shù)逐漸增大，顯示小震震源機(jī)制解趨于一致；編號(hào)4與編號(hào)5的仙游序列U值和ρ值符合前震序列判斷閾值，顯示序列具有一定的前兆性質(zhì)，但對(duì)比編號(hào)3的前兆性質(zhì)不明顯，同期譜振幅相關(guān)系數(shù)處于研究時(shí)段較低水平（圖6）。綜上所述，水口和仙游序列分時(shí)段的序列參數(shù)計(jì)算均顯示出一定的前兆性質(zhì)，兩個(gè)序列最大地震發(fā)生前的序列參數(shù)其前兆性質(zhì)更為明顯，相應(yīng)的譜振幅相關(guān)系數(shù)也相對(duì)較高。

3.2 譜振幅相關(guān)系數(shù)與震源機(jī)制解

利用P初動(dòng)求解了研究時(shí)段內(nèi)水口序列ML3.0級(jí)以上地震震源機(jī)制，顯示具有較好的一致性，與譜振幅相關(guān)系數(shù)在高值波動(dòng)（0.90～0.98）所顯示出的震源機(jī)制解一致性好是吻合的；利用Hash方法求解了研究時(shí)段內(nèi)仙游地區(qū)150次ML2.0級(jí)以上震源機(jī)制，其震源機(jī)制解一致性較好，但其譜振幅相關(guān)系數(shù)在0.72～0.81之間變化，數(shù)值略低，顯示出震源機(jī)制解一致性程度不高，分析認(rèn)為這種差異性的存在可能與應(yīng)力水平的高低有關(guān)。

3.3 譜振幅相關(guān)系數(shù)與應(yīng)力水平

從譜振幅相關(guān)系數(shù)的時(shí)序變化來(lái)看，在兩個(gè)序列最大地震發(fā)生前，譜振幅相關(guān)系數(shù)均存在高值上升而后下降的變化，顯示出序列顯著地震前局地的區(qū)域應(yīng)力場(chǎng)有增強(qiáng)的現(xiàn)象。從譜振幅相關(guān)系數(shù)數(shù)值來(lái)看，水口序列高于仙游序列。本文水口序列譜振幅相關(guān)系數(shù)研究時(shí)段（2007—2008年）前后，福建省正處于ML4.0級(jí)地震集中活躍時(shí)期，2007—2009年福建省陸域及近海發(fā)生10次ML4.0級(jí)以上地震（含水口地區(qū)3次），是我省1971年以來(lái)陸域ML4級(jí)地震最活躍時(shí)期，其中9次發(fā)生在以水口為中心半徑200 km范圍內(nèi)，水口序列的譜振幅相關(guān)系數(shù)的高值與福建地區(qū)整體處于高應(yīng)力背景水平有關(guān)。仙游序列譜振幅相關(guān)系數(shù)數(shù)值略低，仙游序列ML4.0級(jí)地震活動(dòng)前福建省ML4.0級(jí)以上地震持續(xù)平靜超過(guò)3年，且ML3級(jí)地震活動(dòng)水平偏低，仙游序列地震活動(dòng)平靜后，福建地區(qū)地震活動(dòng)也逐漸平靜，截止2015年底超過(guò)2年福建地區(qū)未發(fā)生ML4.0級(jí)以上地震。由此可見(jiàn)，仙游序列發(fā)生前后福建地區(qū)中等地震活動(dòng)水平低，福建地區(qū)處于較低的應(yīng)力背景水平，與其譜振幅相關(guān)系數(shù)偏低一致。仙游序列中等地震活動(dòng)主要受制于小區(qū)域的應(yīng)力調(diào)整活動(dòng)，不存在福建整體區(qū)域的高應(yīng)力背景。

致謝：中國(guó)地震局預(yù)測(cè)研究所崔子健博士為本研究提供了計(jì)算程序，并在程序使用等方面提供了悉心的幫助與指導(dǎo)；福建數(shù)字地震臺(tái)網(wǎng)為本研究提供了高質(zhì)量的觀測(cè)數(shù)據(jù)。作者在此一并表示感謝。

[1]陳颙.用震源機(jī)制一致性作為描述地震活動(dòng)性的新參數(shù)[J].地球物理學(xué)報(bào)，1978，21（2）：142-459.

[2]刁桂苓，于新昌.唐山地震前后京、津、唐、張地區(qū)的綜合斷層面解[J].西北地震學(xué)報(bào)，1980，2（3）：39-47.

[3]刁桂苓，于利民，李欽祖.強(qiáng)震前后震源區(qū)應(yīng)力場(chǎng)變化一例[J].地震學(xué)報(bào)，1994，16（1）：64-69.

[4]王俊國(guó)，刁桂苓.千島島弧大震前哈佛大學(xué)矩心矩張量（CMT）解一致性的預(yù)測(cè)意義[J].地震報(bào)，2005，27（2）：178-183.

[5]刁桂苓，徐錫偉，陳于高，等.汶川MW7.9和集集MW7.6地震前應(yīng)力場(chǎng)轉(zhuǎn)換現(xiàn)象及其可能的前兆意義[J].地球物理學(xué)報(bào)，2011，54（1）：128-136.

[6]Lund B，Boevarsson R.Correlation of microearthquake body-wave spectral amplitudes[J].Bull.Seism.Soc.Am.，2002，92（6）：2 419-2 433.

[7]崔子健，李志雄，陳章立，等.判別小震群序列類型的新方法研究——譜振幅相關(guān)分析法[J].地球物理學(xué)報(bào)，2012，55（5）：1 718-1 724.

[8]黃浩，付虹.2008年以來(lái)滇西地區(qū)地震序列的譜振幅相關(guān)系數(shù)變化特征[J].地震學(xué)報(bào)，2014，36（4）：631-639.

[9]王瑩，趙韜，劉春.陜西及周鄰地區(qū)幾個(gè)地震叢序列的震源譜振幅相關(guān)研究[J].高原地震，2015，27（1）：1-7.

[10]黃浩，付虹，王培玲，等.云南魯?shù)?.5級(jí)地震序列的譜振幅相關(guān)系數(shù)變化特征[J].地震研究，2014，37（4）：532-536.

[11]崔子健，李志雄，陳章立，等.云南景谷MS6.6、云南滄源MS5.5地震譜振幅相關(guān)系數(shù)特征分析[J].地震研究，2015，38（4）：535-540.

[12]崔子健，李志雄，陳章立，等.中國(guó)東北2次中強(qiáng)地震譜振幅相關(guān)系數(shù)研究[J].中國(guó)地震，2105，31（2）：218-226.

[13]張致偉，喬慧珍，吳朋，等.注水誘發(fā)地震的譜振幅相關(guān)系數(shù)及視應(yīng)力研究[J].地震研究，2015，38（1）：42-50.

[14]朱航，劉杰，陳天長(zhǎng).采用體波譜振幅相關(guān)系數(shù)方法研究地震序列的震源機(jī)制變化過(guò)程[J].地震，2006，26（2）：1-11.

[15]李霞，董翔，劉希強(qiáng)，等.采用體波譜振幅相關(guān)系數(shù)方法研究地震集中區(qū)震源機(jī)制變化過(guò)程[J].西北地震學(xué)報(bào)，2010，32（3）：231-236.

[16]劉澤民，劉東旺，夏仕安，等.體波譜振幅相關(guān)系數(shù)法在肥東地震序列震后判定中的應(yīng)用[J].地震地磁觀測(cè)與研究，2010，31（5）：13-16.

[17]Brune J N.Tectonic stress and the spectra of seismic shearwavesfromearthquake[J].JGeophyRes，1970，75（26）：4 997-5 009.

[18]Chael E P.Spectral scaling of earthquakes in the Miramichi region of New Brunswick[J].Bull.Seism.Soc.Am.，1987，77（2）：347-365.

[19]黃玉龍，鄭斯華，劉杰，等.廣東地區(qū)地震動(dòng)衰減和場(chǎng)地響應(yīng)的研究[J].地球物理報(bào)，2003，46（1）：54-61.

[20]Atkinson G M，Mereu R F.The shape of ground motion attenuationcurvesinsoutheasternCanada[J].Bull Seismol Soc Am，1992，82（5）：2 014-2 031.

[21]Moya A，Aguirre J，Irikura K.Inversion of source parameters and site effects from strong ground motion records using geneticalgorithms[J].Bull.Seism.Soc.Am.，2000，90（4）：977-992.

[22]國(guó)家地震局預(yù)測(cè)預(yù)防司.測(cè)震學(xué)分析預(yù)報(bào)方法[M].北京：地震出版社，1997.

[23]梅世蓉，馮德益，張國(guó)民，等.中國(guó)地震預(yù)報(bào)概論[M].北京：地震出版社，1993.

[24]吳忠良.關(guān)于b值應(yīng)用于地震趨勢(shì)預(yù)測(cè)的討論[J].地震學(xué)報(bào)，2001，23（5）：548-551.

Variation Characteristics of the Correlation Coefficient of Spectral Amplitude of Shuikou and Xianyou Earthquake Sequences in Fujian

YUAN Liwen，LI Qiang，CHEN Zhiyong
（Earthquake Administration of Fujian Province，F(xiàn)uzhou 350003，China）

By using digital seismic waves recorded by Fujian Digital Seismic Network，the paper calculated the correlation coefficient of spectral amplitude of Shuikou and Xianyou earthquake sequences.The results shows that the correlation coefficient of spectral amplitude of Shuikou varies between 0.90-0.98，the high value of correlation coefficient indicates that the focal mechanism of the Shuikou small earthquakes are more consistent; the correlation coefficient of spectral amplitude of Xianyou varies between 0.72-0.81，the lower value of correlation coefficient indicates that the focal mechanism of the Shuikou small earthquakes are less consistent. Considering the seismicity background in Fujian area，the high correlation coefficient of spectral amplitude is related to the high background stress level in Fujian.The rising changes of correlation coefficient of spectralamplitude before maximum earthquake may be caused by the enhance of regional stress level.

Earthquake sequence；Correlation coefficient of spectral amplitude；Consistency of focal mechanism；Stress level

P315.31

A

1001-8662（2016）03-0035-08

10.13512/j.hndz.2016.03.005

2016-07-20

福建省地震局青年科技基金（Y201302）資助

袁麗文（1982-），女，工程師，碩士，主要從事地震活動(dòng)性及數(shù)字地震學(xué)研究工作.

E-mail：lwysz@126.com.

袁麗文，李強(qiáng)，陳智勇.福建水口及仙游地震序列的譜振幅相關(guān)系數(shù)變化特征[J].華南地震，2016，36（3）：35-42.YUAN Liwen，LI Qiang，CHEN Zhiyong.Variation Characteristics of the Correlation Coefficient of Spectral Amplitude of Shuikou and Xianyou Earthquake Sequences in Fujian[J].South china journal of seismology，2016，36（3）：35-42.]