2023 年8—10月中國大陸地區(qū)M≥4.0 地震震源機制解測定＊

梁姍姍 鄒立曄 劉艷瓊 任 梟

（中國地震臺網(wǎng)中心，北京 100045）

0 引言

據(jù)中國地震臺網(wǎng)測定，北京時間2023 年8 月1 日0 時—2023 年10 月31 日24 時，中國大陸地區(qū)共發(fā)生M≥4.0 地震26 次（表1）。本文對這些地震進行震源機制反演，進一步豐富中國測震臺網(wǎng)地震編目產(chǎn)出內(nèi)容，以期為地球科學(xué)研究人員提供更為詳實的基礎(chǔ)成果資料[1-7]。

表1 2023 年8—10 月中國大陸地區(qū)M≥4.0 地震基本參數(shù)Table 1 Events list of the earthquakes with M≥4.0 in the mainland of China

1 資料與方法

本文使用中國地震臺網(wǎng)固定臺站記錄的寬頻地震波形資料，選取臺站方位角覆蓋好、波形質(zhì)量高、距震中500 km 以內(nèi)的臺站參與震源機制反演計算。圖1 給出了所使用臺站的分布。

圖1 26 次M≥4.0 以上地震震中位置和所用臺站分布Fig.1 Map of the epicenters of the 26 earthquakes with M≥4.0 and the used stations

本文采用ISOLA 近震全波形方法[8-9]，對中國大陸地區(qū)M≥4.0 地震進行震源機制反演，并確定其最佳震源機制解和矩心深度。反演時，采用中國地震臺網(wǎng)中心速報正式結(jié)果作為參考參數(shù)，以震中為起始點，淺源地震和中源地震矩心深度的搜索步長分別為1 km 和10 km，搜索范圍分別為1～30 km 和10～200 km，地震矩心時間偏移的搜索范圍為發(fā)震時刻前后2.25 s，時間步長為0.1 s。文中采用偏矩張量反演模式，基于單點源模型在時間和空間范圍內(nèi)搜索震源機制最優(yōu)解，并將震源機制反演得到的波形互相關(guān)系數(shù)作為震源深度的函數(shù)，反演不同深度的震源機制解，以最大波形擬合系數(shù)對應(yīng)的矩心深度和震源機制解為最佳結(jié)果。

本文震源機制反演使用的一維速度模型參考了Crust2.0 全球速度模型[10]，應(yīng)用離散波數(shù)法[11-12]計算得到格林函數(shù)，采樣頻率為1 Hz。

2 結(jié)果測定

經(jīng)過對這些地震事件波形的預(yù)處理和震源機制反演，本文得到了26 次M≥4.0 地震事件的震源機制解，具體震源參數(shù)如表2 所示。隨后，我們利用震源機制三角形分類方法將本文所得到的震源機制解進行分類，其中，逆斷型9 次，走滑型14 次，正斷型3 次（圖2）。

圖2 2023 年8—10 月中國大陸地區(qū)M≥4.0 地震類型圖解Fig.2 The ternary diagram of the focal mechanisms in the mainland of China during August to October 2023

表2 2023 年8—10 月中國大陸地區(qū)M≥4.0 地震的震源機制結(jié)果Table 2 Focal mechanisms of the earthquakes with M≥4.0 in the mainland of China

3 2023 年山東平原MS5.5 地震儀器烈度計算

隨著國家地震烈度速報與預(yù)警工程的實施和試運行，加速度計和烈度計數(shù)據(jù)已逐步得到應(yīng)用。2023 年8 月6 日2 時33 分，山東德州市平原縣發(fā)生MS5.5 地震，我們使用震中周邊地區(qū)部分基本站和一般站記錄的加速度波形資料，在進行基線校正、記錄轉(zhuǎn)換和數(shù)字濾波后，計算峰值加速度和峰值速度，最后根據(jù)儀器烈度計算的國家標(biāo)準(zhǔn)《中國地震烈度表》（GB/T 17742—2020）得到單個測點的儀器烈度值。圖3 為山東平原MS5.5 地震儀器烈度分布，可以看出我們所得到的結(jié)果與現(xiàn)場考察結(jié)果[13]吻合度較好，二者在一定程度上的對應(yīng)，說明儀器烈度可以較為客觀地反映震區(qū)災(zāi)害情況。

致謝

感謝國家海洋環(huán)境預(yù)報中心徐志國正高級工程師和應(yīng)急管理部國家自然災(zāi)害防治研究院張廣偉副研究員在震源機制解結(jié)果驗證方面和文章成文過程中給予的幫助。