2022年10月—2023年2月中國大陸地區 M≥4.0地震震源機制解測定

梁姍姍　鄒立曄　劉艷瓊

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

關鍵詞震源機制；震源參數

中圖分類號： P315.3+3文獻標識碼： A文章編號：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

引言

據中國地震臺網測定，北京時間2022年10月1 日0時—2023年2月28日24時，中國大陸地區共發生 M≥4.0地震47次（表1）。其中，6.0級以上地震1次，為2023年1月30日新疆阿克蘇地區沙雅縣 M6.1地震。本文對這些地震進行震源機制反演，進一步豐富中國測震臺網地震編目產出內容，以期為地球科學研究人員提供更為詳實的基礎成果資料[1-5]。

1 資料與方法

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

本文采用 ISOLA 近震全波形方法[6-7]，對中國大陸地區 M≥4.0地震進行震源機制反演，并確定其最佳震源機制解和矩心深度。反演時，采用中國地震臺網中心速報正式結果作為參考參數，以震中為起始點，淺源地震和中深源地震矩心深度的搜索步長分別為1 km 和10 km，搜索范圍分別為1—30 km 和 500—600 km，M4.0—5.9和 M6.0以上地震矩心時間偏移的搜索范圍為發震時刻前后2.25 s 和5.0 s，時間步長為0.1 s。文中采用偏矩張量反演模式，基于單點源模型在時間和空間范圍內搜索震源機制最優解，并將震源機制反演得到的波形互相關系數作為震源深度的函數，反演不同深度的震源機制解，以最大波形擬合系數對應的矩心深度和震源機制解為最佳結果。

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

2 結果分析

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

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

致謝

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

參考文獻

［1］梁姍姍，鄒立曄，劉艷瓊，等.2021年12月中國大陸地區 M≥4.0地震震源機制解測定[J].地震科學進展，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月中國大陸地區 MS ≥4.0地震震源機制解測定[J].地震科學進展，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月中國大陸地區 MS ≥4.0地震震源機制解測定[J].地震科學進展，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月中國大陸地區 MS ≥4.0地震震源機制解測定[J].地震科學進展，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月中國大陸地區 M≥4.0地震震源機制解測定[J].地震科學進展，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