2023 年8—10月中国大陆地区M≥4.0 地震震源机制解测定＊

梁姗姗 邹立晔 刘艳琼 任 枭

（中国地震台网中心，北京 100045）

0 引言

据中国地震台网测定，北京时间2023 年8 月1 日0 时—2023 年10 月31 日24 时，中国大陆地区共发生M≥4.0 地震26 次（表1）。本文对这些地震进行震源机制反演，进一步丰富中国测震台网地震编目产出内容，以期为地球科学研究人员提供更为详实的基础成果资料[1-7]。

表1 2023 年8—10 月中国大陆地区M≥4.0 地震基本参数Table 1 Events list of the earthquakes with M≥4.0 in the mainland of China

1 资料与方法

本文使用中国地震台网固定台站记录的宽频地震波形资料，选取台站方位角覆盖好、波形质量高、距震中500 km 以内的台站参与震源机制反演计算。图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]，对中国大陆地区M≥4.0 地震进行震源机制反演，并确定其最佳震源机制解和矩心深度。反演时，采用中国地震台网中心速报正式结果作为参考参数，以震中为起始点，浅源地震和中源地震矩心深度的搜索步长分别为1 km 和10 km，搜索范围分别为1～30 km 和10～200 km，地震矩心时间偏移的搜索范围为发震时刻前后2.25 s，时间步长为0.1 s。文中采用偏矩张量反演模式，基于单点源模型在时间和空间范围内搜索震源机制最优解，并将震源机制反演得到的波形互相关系数作为震源深度的函数，反演不同深度的震源机制解，以最大波形拟合系数对应的矩心深度和震源机制解为最佳结果。

本文震源机制反演使用的一维速度模型参考了Crust2.0 全球速度模型[10]，应用离散波数法[11-12]计算得到格林函数，采样频率为1 Hz。

2 结果测定

经过对这些地震事件波形的预处理和震源机制反演，本文得到了26 次M≥4.0 地震事件的震源机制解，具体震源参数如表2 所示。随后，我们利用震源机制三角形分类方法将本文所得到的震源机制解进行分类，其中，逆断型9 次，走滑型14 次，正断型3 次（图2）。

图2 2023 年8—10 月中国大陆地区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 月中国大陆地区M≥4.0 地震的震源机制结果Table 2 Focal mechanisms of the earthquakes with M≥4.0 in the mainland of China

3 2023 年山东平原MS5.5 地震仪器烈度计算

随着国家地震烈度速报与预警工程的实施和试运行，加速度计和烈度计数据已逐步得到应用。2023 年8 月6 日2 时33 分，山东德州市平原县发生MS5.5 地震，我们使用震中周边地区部分基本站和一般站记录的加速度波形资料，在进行基线校正、记录转换和数字滤波后，计算峰值加速度和峰值速度，最后根据仪器烈度计算的国家标准《中国地震烈度表》（GB/T 17742—2020）得到单个测点的仪器烈度值。图3 为山东平原MS5.5 地震仪器烈度分布，可以看出我们所得到的结果与现场考察结果[13]吻合度较好，二者在一定程度上的对应，说明仪器烈度可以较为客观地反映震区灾害情况。

致谢

感谢国家海洋环境预报中心徐志国正高级工程师和应急管理部国家自然灾害防治研究院张广伟副研究员在震源机制解结果验证方面和文章成文过程中给予的帮助。