石玉燕 張春鵬 顏 啟 苗慶杰

1 山東省地震局，濟(jì)南市文化東路20號， 250014

北京時間2022-01-15南太平洋島國湯加發(fā)生規(guī)模罕見的海底火山噴發(fā)，激發(fā)全球范圍的地震波、海嘯波、聲波和次聲波，被大量科學(xué)儀器所記錄，受到全球科學(xué)家的廣泛關(guān)注。火山噴發(fā)時釋放出巨大能量，其中部分能量以地震波形式傳播，即地震波輻射能量，簡稱地震輻射能或地震波能量。與火山噴發(fā)釋放的能量不同，地震波能量是一個可用地震學(xué)方法科學(xué)測量的物理量。本文利用山東地震臺網(wǎng)70套寬頻帶地震儀記錄的2022-01-15湯加火山噴發(fā)的原始地震波形數(shù)據(jù)，測算本次火山噴發(fā)事件的面波震級MS和能量震級Me，為地震臺網(wǎng)今后探尋火山噴發(fā)事件的參數(shù)速報工作提供一定參考。

1 數(shù)據(jù)與方法

1.1數(shù)據(jù)來源

湯加火山噴發(fā)后，截取山東臺網(wǎng)128個測震臺站記錄的該次事件地震波數(shù)據(jù)。對原始資料進(jìn)行去儀器響應(yīng)、濾波與仿真等預(yù)處理后，選用70個信噪比高、質(zhì)量好的臺站記錄資料進(jìn)行震級測定。在臺基校正時選用山東臺網(wǎng)2014～2021年記錄的15次區(qū)域內(nèi)天然地震和5次非天然地震事件(表1)。另外，在地震波記錄特征分析中還選取2022-01-27斐濟(jì)群島MS6.0地震的地震波記錄數(shù)據(jù)。用于計算的70個寬頻帶測震臺站分布如圖1所示，所有臺站的采樣率均為100 Hz。

圖1 研究使用的70個臺站分布Fig.1 Distribution of 70 seismic stations in study area

1.2 地震波記錄

本文選用的臺站到火山噴發(fā)口的震中距均在80°～90°范圍內(nèi)，相應(yīng)的地震波記錄震相發(fā)育良好并相互分開，可避免不同類型波相互交叉干擾[1-2]。圖2為SD/TCH臺三分向記錄的該震中距范圍內(nèi)兩次事件的地震波形圖，其中圖2(a)為2022-01-15湯加火山噴發(fā)的地震波形圖，圖2(b)為2022-01-27斐濟(jì)群島MS6.0地震的地震波形圖。圖中，三分向自上而下分別為垂直向UD、東西向EW和南北向NS；藍(lán)色P、S分別表示縱波和橫波到時；紅色LZ、LN和LE分別表示三分向最大振幅及最大振幅的周期值。

表1 事件參數(shù)

圖2 TCH臺三分向地震波記錄Fig.2 Three component waveforms record at TCH station

從圖2可以看出，火山噴發(fā)的垂直向P波初至比天然地震更加清晰與尖銳；火山噴發(fā)地震波垂直向比水平向最大振幅大、能量強(qiáng)，而天然地震水平向比垂直向最大振幅大、三分向地震波能量較均衡；雖然火山噴發(fā)比天然地震的面波最大振幅小、周期小，但面波波列的持續(xù)時間長、衰減慢。

為更清楚地了解火山噴發(fā)地震波頻率域特征，分別對以上兩條三分向地震波進(jìn)行時頻分析，得到火山噴發(fā)和天然地震的地震波時頻分布及頻譜曲線，結(jié)果如圖3所示。由圖可見，火山噴發(fā)比天然地震整個波列的頻率范圍大、主頻低。

圖3 火山噴發(fā)與天然地震的地震波時頻分析Fig.3 Seismic wave time-frequency analysis of volcanic eruption and earthquake

2 面波震級

2.1 測定方法

面波震級MS是用面波測定的一種震級標(biāo)度，采用中長周期地震儀記錄的地震波質(zhì)點運動最大值進(jìn)行測定[3-5]。我國在1966-01以后采用郭履燦等[6]提出的以北京白家疃地震臺為基準(zhǔn)的面波震級公式，并一直沿用至今。由于受我國地震臺網(wǎng)孔徑限制，在測定面波震級時均使用震中距為 130°以內(nèi)的臺站數(shù)據(jù)，計算公式為：

MS=lg(A/T)max+σ(Δ)，

(1°<Δ<130°)

(1)

對于較小的地震事件，遠(yuǎn)距離觀測資料信噪比較低，難以收集按方位角均勻分布的面波資料，從而無法有效消除輻射花樣和臺基的影響。利用區(qū)域事件可以對地震臺網(wǎng)進(jìn)行面波震級測量的量規(guī)處理，獲取臺基校正值，從而有效壓制網(wǎng)內(nèi)地震事件輻射花樣的影響。設(shè)地震臺網(wǎng)由L個臺站構(gòu)成，可利用N個量規(guī)事件，根據(jù)垂直分向面波資料，利用式(1)計算每個事件的面波震級。對每個臺站求取臺基校正值：

ΔMS(i)=

(2)

式中，m

2.2 湯加火山噴發(fā)面波震級測定

首先將湯加火山噴發(fā)事件原始波形數(shù)據(jù)進(jìn)行去儀器響應(yīng)與濾波預(yù)處理，使用區(qū)域地震臺網(wǎng)MSDP處理軟件人工拾取P、S波震相到時，然后將70個臺站的速度記錄仿真成位移記錄進(jìn)行最大面波振幅與周期標(biāo)注[7]，最后使用AK135全球速度模型進(jìn)行地震基本參數(shù)測算，得到本次火山噴發(fā)70個臺站校正前的面波震級值(表2)。將表1中20個區(qū)域內(nèi)地震事件速度記錄仿真成地動位移記錄，測量垂直向面波最大振幅，計算各臺站面波震級M臺，并與國家臺網(wǎng)中心編目震級M進(jìn)行比較，得到每個臺站的震級校正值M臺-M。將70個臺站的震級校正值繪制成面波震級校正值柱狀圖(圖4(a))。

表2為70個臺站校正前后測定的湯加火山噴發(fā)面波震級，由表可知，湯加火山噴發(fā)事件面波震級為MS5.674。依據(jù)表2結(jié)果繪制本次湯加火山噴發(fā)面波震級柱狀圖(圖4(b))。

3 地震波輻射能量與能量震級

3.1 測定方法

彈性波釋放的能量與地面運動速度的平方成正比，假設(shè)地震為點源，震源周圍為均勻介質(zhì)，則可利用遠(yuǎn)震 P 波的垂直向記錄測定地震能量[8]：

(3)

ES-MS的經(jīng)驗關(guān)系為：

lgES=1.5MS+4.4

(4)

能量震級Me定義為：

(5)

3.2 湯加火山噴發(fā)的地震波輻射能量估算

基于表2中70個寬頻帶測震臺站測定的面波震級MS，利用式(4)和式(5)分別得到湯加火山噴發(fā)每個臺站的地震波輻射能量與能量震級，然后將70個臺站測算的結(jié)果求均值得到湯加火山噴發(fā)的地震波輻射能量ES與能量震級Me。

圖5為70個寬頻帶臺站測算的本次湯加火山噴發(fā)的地震波輻射能量與能量震級圖。由圖可知，本次火山噴發(fā)的地震波輻射能量為8.155×1012J，能量震級為Me5.704。

4 與其他機(jī)構(gòu)結(jié)果對比

將本文測定結(jié)果與USGS和GFZ等國際機(jī)構(gòu)發(fā)布的結(jié)果(表3)進(jìn)行比較。

表2 70個臺站測算的湯加火山噴發(fā)震級參數(shù)

圖4 湯加火山噴發(fā)的面波震級柱狀圖Fig.4 Surface wave magnitude of Tonga volcanic eruption

圖5 湯加火山噴發(fā)的地震波輻射能量與能量震級Fig.5 Radiation energy and magnitude of seismic wave emitted by Tonga volcano

表3 湯加火山噴發(fā)事件國際地震機(jī)構(gòu)測定的地震基本參數(shù)

從表3可以看出，本文所得面波震級MS5.674比USGS結(jié)果小0.126，比GFZ結(jié)果小0.026。為進(jìn)一步查找本文測定震級偏小的原因，將表2中70個臺站測定的震級結(jié)果與USGS和GFZ結(jié)果進(jìn)行震級偏差計算，并統(tǒng)計不同震級偏差范圍的臺站數(shù)量，繪制震級偏差柱狀圖(圖6)。由圖可見，90%臺站測定的震級小于USGS結(jié)果。主要原因為山東地震臺網(wǎng)為區(qū)域臺網(wǎng)，臺網(wǎng)孔徑較小，臺站位置比較集中，而USGS使用全球臺網(wǎng)數(shù)據(jù)，臺站方位分布相對比較均衡。同時，地震多普勒效應(yīng)會導(dǎo)致某一方位的震動減弱[9]，這也可能是本文測定震級偏小的原因之一。

圖6 本文與國際機(jī)構(gòu)測定湯加火山噴發(fā)面波震級差值統(tǒng)計Fig.6 Statistics of surface wave magnitude difference measured by this paper and international seismic agency

綜上可知，本文使用山東地震臺網(wǎng)寬頻帶臺站測定湯加火山噴發(fā)的面波震級比較準(zhǔn)確可靠，可以在今后地震臺網(wǎng)速報火山噴發(fā)的地震參數(shù)工作中參考與使用。

5 結(jié) 語

本文是對火山噴發(fā)地震學(xué)震級測定的首次嘗試。通過對山東地震臺網(wǎng)70個寬頻帶測震臺站記錄的2022-01-15湯加火山噴發(fā)原始地震波形數(shù)據(jù)進(jìn)行濾波與仿真預(yù)處理，經(jīng)過面波震級的臺基校正，給出該次火山噴發(fā)的地震學(xué)震級與地震波輻射能量。主要得到以下結(jié)論：

1)本文計算的湯加火山噴發(fā)面波震級為MS5.674、能量震級為Me5.704，與GFZ測定震級基本一致，比USGS測定震級偏小0.126。臺網(wǎng)孔徑較小、臺站位置比較集中和地震多普勒效應(yīng)可能是本文測量震級偏小的原因。

2)本次湯加火山爆炸指數(shù)(VEI)為5級[10-11]，火山噴發(fā)猛烈，但本文測定的地震波輻射能量僅為8.155×1012J，可見該次火山噴發(fā)以地震波輻射方式釋放的能量較少。其主要原因為火山噴發(fā)時，絕大部分能量是以沖擊波的形式穿過大氣層迅速釋放到空中[12]，這種能量的快速釋放會抑制深部應(yīng)力的擠壓與拉伸。

3)本次火山噴發(fā)的地震波與類似震級、類似震中距的淺源地震相比，具有垂直向P波初動尖銳、面波衰減速度慢、持續(xù)時間長、頻率范圍大、主頻低等特點。

致謝：感謝中國地震臺網(wǎng)中心梁建宏研究員和山東省地震局崔鑫博士對本文研究給予的指導(dǎo)與幫助。