梁姍姍　鄒立曄　劉艷瓊

摘要利用中國(guó)地震臺(tái)網(wǎng)記錄的寬頻帶波形資料，采用近震全波形反演方法得到2022年10月1 日—2023年2月28日發(fā)生在中國(guó)大陸地區(qū)的 M≥4.0共46次地震的震源機(jī)制解以及2023年1月30日新疆沙雅 M6.1地震的矩心位置和矩心偏移時(shí)間。結(jié)果顯示逆斷型15次，走滑型15次，正斷型14次，未知型2次。

關(guān)鍵詞震源機(jī)制；震源參數(shù)

中圖分類號(hào)： P315.3+3文獻(xiàn)標(biāo)識(shí)碼： A文章編號(hào)：2096-7780（2023）04-0185-07

doi：10.19987/j.dzkxjz.2023-019

Determination of focal mechanism solutions of the earthquakes with M≥4.0occurred in the mainland of China during October 2022 to February 2023

Liang Shanshan，Zou Liye，Liu Yanqiong

（China Earthquake Networks Center， Beijing 100045， China）

AbstractIn this paper，the regional full waveform inversion using the broadband waveforms recorded by China Seismic Network were conducted，and the focal mechanism solutions of the 46 earthquakes with M≥4.0 occurred in the mainland of China from October 2022 to February 2023 were obtained. The centroid location and time shift of the M6.1 Shaya，Xinjiangearthquakeweredeterminatedbyagridsearchinthehorizontalplane. Thetypesofthesefocal mechanism solutions show 15 reverse faulting，15 strike-slip faulting，14 normal faulting and 2 odd earthquakes.

Keywordsfocal mechanism; source parameters

引言

據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定，北京時(shí)間2022年10月1 日0時(shí)—2023年2月28日24時(shí)，中國(guó)大陸地區(qū)共發(fā)生 M≥4.0地震47次（表1）。其中，6.0級(jí)以上地震1次，為2023年1月30日新疆阿克蘇地區(qū)沙雅縣 M6.1地震。本文對(duì)這些地震進(jìn)行震源機(jī)制反演，進(jìn)一步豐富中國(guó)測(cè)震臺(tái)網(wǎng)地震編目產(chǎn)出內(nèi)容，以期為地球科學(xué)研究人員提供更為詳實(shí)的基礎(chǔ)成果資料[1-5]。

1 資料與方法

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

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

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

2 結(jié)果分析

經(jīng)過(guò)對(duì)這些地震事件波形的預(yù)處理和震源機(jī)制反演，本文得到了46次 M≥4.0地震事件的震源機(jī)制解，具體震源參數(shù)如表2所示。需說(shuō)明的是，由于2022年11月22日新疆皮山M4.2地震臺(tái)站分布不夠理想，且波形信噪比不高，很難得到該次地震的穩(wěn)定最佳震源機(jī)制結(jié)果。我們利用震源機(jī)制三角形分類方法將本文所得到的震源機(jī)制解進(jìn)行分類，其中，逆斷型15次，走滑型15次，正斷型14次，未知型2次（圖2）。

為減少初始震中和有限斷層長(zhǎng)度對(duì)大地震震源機(jī)制反演結(jié)果的影響，我們采用平面網(wǎng)格搜索法來(lái)進(jìn)一步確定2023年1月30日沙雅 M6.1地震的最佳矩心位置和矩心偏移時(shí)間。綜合固定震中位置確定的沙雅 M6.1地震最佳矩心深度、節(jié)面走向以及烈度分布特征（www.cenc.ac.cn），初步判斷走向161°的節(jié)面為沙雅地震可能的發(fā)震斷層面。我們以初始定位震中和反演矩心深度為參考點(diǎn)，對(duì)走向161°的水平面震源區(qū)進(jìn)行網(wǎng)格化，平面長(zhǎng)度、寬度分別為40 km 和20 km，網(wǎng)格步長(zhǎng)5 km×5 km，共劃分45個(gè)網(wǎng)格點(diǎn)。圖3為水平面內(nèi)網(wǎng)格點(diǎn)分布及搜索的最佳震源機(jī)制解。由圖可知，13號(hào)網(wǎng)格點(diǎn)對(duì)應(yīng)波形擬合系數(shù)最大，對(duì)應(yīng)的最佳震源機(jī)制解參數(shù)：節(jié)面Ⅰ、節(jié)面Ⅱ?qū)?yīng)的走向、傾角和滑動(dòng)角分別為255°、66°、11°和160°、80°、156°;最佳矩心位置為（40.0246°N，82.3454°E）；矩心偏移時(shí)間1.78 s；方差減小量 VR=0.66。

致謝

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

參考文獻(xiàn)

［1］梁姍姍，鄒立曄，劉艷瓊，等.2021年12月中國(guó)大陸地區(qū) M≥4.0地震震源機(jī)制解測(cè)定[J].地震科學(xué)進(jìn)展，2022，52（1）：40-44

Liang S S，Zou L Y，Liu Y Q， et al. Determination the focal mechanism solutions of earthquakes with M≥4.0 occurred in the mainland of China in December 2021[J]. Progress in Earthquake Sciences，2022，52（1）：40-44

［2］梁姍姍，鄒立曄，劉艷瓊，等.2022年1月中國(guó)大陸地區(qū) MS ≥4.0地震震源機(jī)制解測(cè)定[J].地震科學(xué)進(jìn)展，2022，52（2）：89-94

Liang S S，Zou L Y，Liu Y Q，et al. Determination of the focal mechanism solutions of the earthquakes with MS ≥4.0 occurred in the mainland of China in January 2022[J]. Progress in Earthquake Sciences，2022，52（2）：89-94

［3］梁姍姍，鄒立曄，劉艷瓊，等.2022年2—3月中國(guó)大陸地區(qū) MS ≥4.0地震震源機(jī)制解測(cè)定[J].地震科學(xué)進(jìn)展，2022，52（4）：183-189

Liang S S，Zou L Y，Liu Y Q，et al. Determination of the focal mechanism solutions of the earthquakes with MS ≥4.0 occurred in the mainland of China during February to March 2022[J]. Progress in Earthquake Sciences，2022，52（4）：183-189

［4］梁姍姍，鄒立曄，劉艷瓊，等.2022年4—6月中國(guó)大陸地區(qū) MS ≥4.0地震震源機(jī)制解測(cè)定[J].地震科學(xué)進(jìn)展，2022，52（7）：338-343

Liang S S，Zou L Y，Liu Y Q，et al. Determination of focal mechanism solutions of the earthquakes with MS ≥4.0 occurred in the mainland of China during April to June 2022[J]. Progress in Earthquake Sciences，2022，52（7）：338-343

［5］梁姍姍，鄒立曄，劉艷瓊，等.2022年7—9月中國(guó)大陸地區(qū) M≥4.0地震震源機(jī)制解測(cè)定[J].地震科學(xué)進(jìn)展，2022，52（11）：547-552

Liang S S，Zou L Y，Liu Y Q，et al. Determination of focal mechanism solutions of the earthquakes with M≥4.0 occurred in the mainland of China during July to September 2022[J]. Progress in Earthquake Sciences，2022，52（11）：547-552

［6］Sokos E N，Zahradnik J. ISOLA a Fortran code and a Matlab GUI to perform multiple-point source inversion of seismic data[J]. Computers and Geosciences，2008，34（8）：967-977

［7］SokosE， ZahradnikJ. Evaluatingcentroid-moment-tensoruncertaintyinthenewversionofISOLAsoftware[J]. Seismological Research Letters，2013，84（4）：656-665

［8］Bassin C GL，Laske G，Masters G. The current limits of resolution for surface wave tomography in North America[J]. Eos Transactions American Geophysical Union，2000，81（48）：F897

［9］Kennett B L N，Kerry N J. Seismic waves in a stratified half space[J]. Geophysical Journal International，1979，57（3）：557-583

［10］Bouchon M. A simple method to calculate Green s functions for elastic layered media[J]. Bulletin of theSeismologicalSociety of America，1981，71（4）：959-971