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