基于爆炸地震波分析3次爆炸事件的主要特征*

郝春月，鄭 重

(中國地震局地球物理研究所，北京 100081)

2006年10月9日1時35分(世界時,UTC)，全球地震臺網(wǎng)(GSN)記錄到一次爆炸事件。根據(jù)美國地質(zhì)調(diào)查局所屬的國家地震信息中心提供的地震目錄(PDE)，此次爆炸事件發(fā)生的位置為北緯41.294°、東經(jīng)129.094°,震級為mb4.3， 深度為0 km[1]。該事件距離GSN的牡丹江地震臺(MDJ)約3.34°，方位角為186.4°。2009年5月25日0時54分(UTC)，GSN記錄到一次爆炸事件，位置為北緯41.303°、東經(jīng)129.037°，震級為mb4.7，深度為0 km[2]。該爆炸事件發(fā)生的位置距離牡丹江地震臺約3.34°，方位角為187.1°。 2013年2月12日2時57分(UTC)，GSN記錄到一次爆炸事件，位置為北緯41.308°、東經(jīng)129.076°，震級為mb5.1，深度為0 km[3]。該爆炸事件發(fā)生的位置距離牡丹江地震臺約3.33°，方位角為186.7°。這3次爆炸事件均位于朝鮮東北部咸鏡北道吉州郡附近的豐溪里試驗場。W.Y.Kim等[4]曾對2006年爆炸事件和爆炸事件所在區(qū)域的地震事件進行過判據(jù)方面的研究，范娜等[5]對這3次爆炸事件的震級進行過分析。由于這3次爆炸事件發(fā)生的位置較近，震級逐漸遞增，本文中主要對這3次爆炸事件的能量比、相關(guān)性和相對位置進行研究。本文中，對牡丹江地震臺3個分向(垂直向Z/北南向N/東西向E)記錄的3次爆炸事件進行了P波振幅比、體波和面波功率譜比、初動震相相關(guān)系數(shù)的計算與分析。為了更好地研究這3次爆炸事件的相對位置，還引進了GSN白家疃地震臺(BJT)記錄的波形數(shù)據(jù)。牡丹江地震臺位于中國黑龍江省東南部，白家疃地震臺位于北京西北部，同屬于GSN。為方便起見，用E06、E09和E13分別表示2006年、2009年和2013年發(fā)生的3次爆炸事件。

1 資 料

3次爆炸事件均被牡丹江和白家疃地震臺記錄到，牡丹江地震臺記錄的3次爆炸事件都很清晰,信噪比較高,見圖1(a)。白家疃地震臺由于距離較遠(震中距約1 000 km)，只清晰地記錄到了E09和E13。對于E06，經(jīng)過2～4 Hz濾波也能清晰地看到Sg震相，但Pg震相被淹沒在噪聲里,見圖1(b)。

圖1(a)顯示了120 s的波形數(shù)據(jù)。每幅圖上顯示的4個重要震相分別為Pn(首波，即Moho面繞射縱波)、Pg(直達縱波)、Sg(直達橫波)和瑞利波。3次爆炸事件的體波清晰，Z向體波最大振幅的比值如表1所示。由于震級的計算也是利用最大振幅，所以最大振幅的比值在一定程度上描述了3次爆炸事件爆發(fā)能量的關(guān)系。表1給出了最大體波振幅比和最大面波振幅比。根據(jù)表1，把最大體波振幅和最大面波振幅比進行平均，可得出E13最大振幅是E09最大振幅的2.3倍，而E13的最大振幅是E06最大振幅的10.1倍。

在白家疃地震臺記錄的3次爆炸事件中，面波不發(fā)育，E06的Pg震相不清晰，所以沒有進行能量比和相關(guān)系數(shù)的計算，引進白家疃地震臺數(shù)據(jù)的主要目的是對這3次爆炸事件的相對位置進行分析。

表1 牡丹江地震臺Z向記錄的3次爆炸事件的振幅比Table 1 Amplitude ratios between three explosion events recorded by Mudanjiang seismic station in vertical direction (Z)

圖1 2個地震臺垂直向(Z)記錄到的3次爆炸事件的波形Fig.1 Waveforms of three explosion events recorded by two seismic stations in vertical direction (Z)

2 功率譜比的計算

為了了解3次爆炸事件在頻率域的能量分布情況，對3次爆炸事件進行了功率譜的計算與分析。功率譜分析采用Welch平均周期圖法。Welch平均周期圖法是對直接法的改進，即把一長度為N的數(shù)據(jù)xN(n)分成L段(在分段時可允許每一段的數(shù)據(jù)有部分的重疊)，每一段的長度為M,分別求每一段的功率譜，然后加以平均。第i段的功率譜可由下式表示[6]：

(1)

(2)

在此，分別對2個時間窗進行了功率譜估計，一個為P波組的15 s時間窗，另一個為面波組的15 s時間窗。根據(jù)郝春月等[7]的研究結(jié)果可知，前2次事件體波主要能量集中在2～4 Hz，所以首先對3次事件的體波時間窗進行2～4 Hz頻段的帶通濾波，又由于瑞利波的主要頻段集中在1 Hz以下[8]，而后，對面波組時間窗進行了1 Hz以下的低通濾波。計算功率譜所用的15 s P波組和面波組數(shù)據(jù)均為600個采樣點(采樣率40 s-1)，利用Hanning窗，窗長為256點，重疊128點。3次事件的功率譜計算完畢后，可以得出3次事件中每對事件的功率譜點對點的比值。圖2給出了3次爆炸事件垂直向記錄中體波和面波的功率譜和功率譜比值。從圖2可以看出E06、E09、E13爆發(fā)能量是遞增的，譜比值見表2。

圖2 牡丹江地震臺記錄的3次爆炸事件Pg震相和瑞利波之間的功率譜與譜值比Fig.2 Power spectrum densities for the Pg and Rayleigh wave of the three explosion events and the ratios of them

計算功率譜是為了對3次爆炸事件釋放的能量進行對比，并且牡丹江地震臺的傳遞函數(shù)在3次事件之間沒有變化，所以沒有對波形進行去除儀器傳遞函數(shù)的計算。因為在計算3次事件的譜值比的過程中，會抵扣掉傳遞函數(shù)，所以省卻了這個步驟。省掉這個步驟的結(jié)果

表2 E13與E06、E09 Pg和瑞利波的功率譜比Table 2 Power spectrum ratios of E13 to E06 and E09

是圖2顯示的3次事件功率譜的縱軸單位不是(m/s)2,也不是表示能量衰減所用的dB，而是量綱為一。

表2中分頻段給出了3次爆炸事件的功率譜比值。2～4 Hz頻段(Pg震相的優(yōu)勢頻段)，E13與E06的 Pg震相功率譜比值P13/P06=13.9，而E13與E09的譜比值P13/P09=2.6。瑞利波的主要頻段為1 Hz以下，在此頻段，P13/P06= 13，P13/P09=2.7。功率譜代表了能量，把Pg與瑞利波震相的功率譜比平均，E13釋放的能量是E06的13.5倍，是E09的2.7倍。

3 3次爆炸事件的相關(guān)性分析

為了探索3次事件的相似性，對牡丹江地震臺記錄的3次事件的垂直向、北南向和東西向數(shù)據(jù)進行了相關(guān)性分析。為了對相似性進行定量地描述，對3次事件的相關(guān)系數(shù)進行了計算。首先對波形進行2～4 Hz的帶通濾波，然后，選用了3次爆炸頭2 s的數(shù)據(jù)進行了相關(guān)系數(shù)的計算。選擇相關(guān)系數(shù)計算的時間窗一般為1～2 s[9-10]，在此選擇2 s，使波形窗能包含更多的震相信息。

2個時間序列的互相關(guān)函數(shù)可以表示為[11]：

(3)

根據(jù)郝春月等[7]的計算結(jié)果，知道前2次事件的主要能量集中在2～4 Hz，所以對這3次事件進行了該頻段的濾波。再根據(jù)式(3)，計算得出牡丹江地震臺三分向記錄的3次爆炸事件波形相關(guān)系數(shù)的分析結(jié)果，由于篇幅限制，圖片只給出垂直向結(jié)果，見圖3。根據(jù)計算，得出E13與E06垂直向、北南向、東西向波形的最大相關(guān)系數(shù)分別為0.92、 0.88和0.92；E13與E09垂直向、北南向、東西向波形的最大相關(guān)系數(shù)為0.99、 0.99和0.98?？梢奅13與E09的相關(guān)系數(shù)較高，接近于1。E13與E06的相關(guān)系數(shù)較小，這也可能是由于E06的震級較小，波形信噪比較低造成的，因為低信噪比的波形計算會導致精度降低。3次事件的相關(guān)性分析表示，3次事件在它們的主要能量分布頻段高度相關(guān)。這個結(jié)果表明3次爆炸事件發(fā)生的位置相距很近，根據(jù)經(jīng)驗[12]，應不超過5 km。

圖3 E13與E06、E09垂直向記錄的相關(guān)系數(shù)計算結(jié)果Fig.3 Calculated correlation coefficients between E13 and E06，E09 recorded in vertical direction

4 事件的相對位置

3次事件具有特殊性，對它們進行精確定位具有重要意義。

為了精確定位地震位置，一般在地震活動區(qū)建立很多臺站，形成一個區(qū)域地震臺網(wǎng)。臺網(wǎng)中臺站的監(jiān)測能力能夠覆蓋該地區(qū)[13]。根據(jù)區(qū)域地震臺網(wǎng)管理規(guī)定，一個區(qū)域地震臺網(wǎng)的定位誤差小于5 km時，一般定義為Ⅰ類精度，也就是最好的定位結(jié)果。3次爆炸事件發(fā)生位置的周圍沒有區(qū)域地震臺網(wǎng)，而GSN的地震臺站分布稀疏，其中距離3次爆炸事件位置最近的是牡丹江地震臺，震中距為300多千米。在臺站稀疏并且沒有區(qū)域臺網(wǎng)的情況下很難對3次爆炸事件進行精確定位，也就是利用常規(guī)方法很難達到Ⅰ類精度的定位。鑒于確定絕對位置存在困難，本文中準備研究其相對位置。弄清其相對位置，對判斷該地區(qū)以后爆炸事件的發(fā)生位置、偵破該地區(qū)的地質(zhì)構(gòu)造都具有重要意義。

根據(jù)3次爆炸事件的相關(guān)性分析, 判斷3次爆炸事件組成的區(qū)域小于5 km, 而本文中要研究3次爆炸事件的相對位置, 就是要求精度比I類高。 由于E06震級較小，初動不清晰，初動讀數(shù)誤差將使定位精度大大降低，從而不能達到超越I類精度的要求，所以只針對E09和E13，計算其相對位置。

眾所周知，假設(shè)虛波速度不變，震中距越遠，則Sg與Pg的到時差tSg-tPg和Pg與Pn的到時差tPg-tPn就越大[14]。牡丹江地震臺記錄的波形出現(xiàn)了Pn震相(圖1(a))，Pg震相振幅最大。把牡丹江臺記錄的2次事件的Pn震相兩兩對齊，最大的Pg震相進行對比,見圖4(a)。白家疃地震臺距離事件位置1 000多千米,最明顯的震相是Pg和Sg(圖1(b))，把白家疃地震臺記錄的2次事件的Pg震相兩兩對齊，Sg震相的最大振幅進行對比，見圖4(b)。

圖4 牡丹江與白家疃地震臺記錄的E13、E09的相對位置分析Fig.4 Analysis of the relative location between E13 and E09 recorded by MDJ and BJT

圖5 臺站與事件的相對位置Fig.5 Distribution of stations and explosions

圖4(a)給出了牡丹江臺記錄的E13與 E09 Pg震相的波形對比?？梢钥闯觯珽13的Pg震相比E09的晚到，由于牡丹江地震臺位于爆炸事件位置的北偏東方向(圖5)，可以得出E13事件在E09事件的南側(cè)；圖4(b)給出了白家疃臺記錄的E13與E09間Sg震相的波形對比?？梢钥闯?，E13的Sg震相比E09的晚到，由于白家疃地震臺位于爆炸事件位置的西偏南方向(圖5)，可以得出E13在E09的東側(cè)。根據(jù)牡丹江、白家疃和2次爆炸事件的地理方位，可以判斷，E13位于E09的東南方向，見圖6。圖6只是示意圖，并不表示事件的真實位置。

本文結(jié)果與美國NEIC給出的PDE結(jié)果(圖6)不同，尤其是E13事件，牡丹江地震臺的波形很清晰，在Pn對齊后，E13的Pg震相比E09的晚到，即E13的震中距比E09的震中距大，又由于牡丹江地震臺位于爆炸事件的北偏東方向，所以E13位于E09的偏南方向。

圖6 3次爆炸事件的相對位置Fig.6 Relative location of three explosions

5 結(jié) 論

2013年爆炸事件(E13)爆發(fā)能量是2009年爆炸事件(E09)爆發(fā)能量的2倍多，是2006年爆炸事件(E06)爆炸能量的10多倍。

在2～4 Hz頻段內(nèi)，E13與E06和E09三分向波形的最大相關(guān)系數(shù)平均分別為0.90和0.99。分析3次爆炸事件的相關(guān)性可知，在2～4 Hz頻段，E13與E06和E09高度相關(guān)。這表明這3次事件發(fā)生位置相距很近,根據(jù)經(jīng)驗，總孔徑不超過5 km。

E13發(fā)生位置位于E09的東南方向,該結(jié)論與2013年2月5日鄭州晚報的報道[15]相符，該報道暗示2013年爆炸事件發(fā)生在豐溪里試驗場西側(cè)坑道(2009年爆炸事件發(fā)生地)東南向的南側(cè)坑道。

[1] USGS.Earthquakes: Search EQ archives[EB/OL].[2006-10-09].http:∥earthquake.usgs.gov/earthquakes/eqarchives/epic/ [2] USGS.Earthquakes: Search EQ archives[EB/OL].[2009-5-25].http:∥earthquake.usgs.gov/earthquakes/eqarchives/epic/ [3] USGS.Earthquakes: Search EQ archives[EB/OL].[2013-02-12].http:∥earthquake.usgs.gov/earthquakes/eqarchives/epic/ [4] Kim W Y, Richards P G.North Korean nuclear test: Seismic discrimination at low yield[J].Eos, Transactions American Geophysical Union， 2007,88(13/14):157-161.

[5] 范娜,趙連鋒,謝小碧,等.朝鮮核爆的Rayleigh波震級測量[J].地球物理學報,2013,56(3):906-915.Fan Na, Zhao Lian-feng, Xie Xiao-bi, et al.Measurement of Rayleigh-wave magnitudes for North Korean nuclear tests[J].Chinese Journal of Geophysics, 2013,56(3):906-915.

[6] 胡廣書.數(shù)字信號處理:理論、算法與實現(xiàn)[M].北京:清華大學出版社,2003:324-341.

[7] 郝春月,鄭重.2次爆炸事件的相關(guān)性與能量比研究[J].爆炸與沖擊,2010,30(5):535-540.Hao Chun-yue, Zheng Zhong.Relativity and energy ratio between two explosion events[J].Explosion and Shock Waves, 2010,30(5):535-540.

[8] 時振梁,張少泉,趙榮國,等.地震工作手冊[M].北京:地震出版社,1992:41-47.

[10] Mykkeltveit S, ?steb?l K, Doornbos D J, et al.Seismic array configuration optimization[J].Bulletin of the Seismological Society of America, 1983,73(1):173-186.

[11] Harjes H P.Design and siting of a new regional seismic array in Central Europe[J].Bulletin of the Seismological Society of America, 1990,80(6B):1801-1817.

[12] Ingate S F, Husebye E S, Christoffersson A.Regional arrays and optimum data processing schemes[J].Bulletin of the Seismological Society of America, 1985,75(4):1155-1177.

[13] 彼得·鮑曼.新地震觀測實踐手冊[M].中國地震局監(jiān)測預報司,譯.北京:地震出版社,2006:321-367.

[14] 傅淑芳,劉寶誠,李文藝.地震學教程[M].北京:地震出版社,1980:137-153.

[15] 鄭州晚報數(shù)字報:A12版:國際新聞[N/OL].[2013-02-05]http:∥zzwb.zynews.com/html/2013-02/05/content_440861.htm