2023年甘肃积石山MS6.2地震对周围断层的应力影响及对余震的触发作用

摘要：2023年12月18日，甘肃省临夏州积石山县发生MS6.2地震，震中位于拉脊山断裂带。为深度剖析本次地震对周围断层和后续地震的静态库仑破裂应力影响，文章基于地震破裂模型参数和弹性半空间模型，计算积石山MS6.2地震在周围主要断层上不同深度产生的静态库仑破裂应力变化，并评估这些应力变化对未来地震危险性的影响。结果显示：本次地震使青海南山—循化南山断裂、拉脊山北缘断裂和拉脊山南缘断裂部分区域的库仑破裂应力显著增加，青海南山—循化南山断裂东段的库仑破裂应力增加远超出静态应力触发阈值，并达到了0.022 MPa，表现出较高地震危险性；其他断层也呈现出不同程度的库仑应力变化，拉脊山北缘断裂和拉脊山南缘断裂部分区域的应力卸载量达到了万帕，青海南山—循化南山断裂除东段外的其他分段和庄浪河断裂的应力卸载量达到了百帕。对不同震源深度下断层库仑破裂应力变化进行对比分析，发现深度变化对此次地震库仑破裂应力变化模式影响很大。文章还计算了本次地震在较大余震断层面上的库仑破裂应力，结果表明：甘肃积石山MS6.2地震在后续三次MS≥4地震的断层面上产生的库仑破裂应力分别为0.024 MPa、0.033 MPa和0.034 MPa，均超过触发阈值（0.01 MPa），表明MS6.2地震对这三次地震具有明显的触发作用。文章可为该地区未来地震可能性与危险性的评估提供一定的参考。

关键词：甘肃积石山MS6.2地震; 库仑破裂应力; 地震危险性; 拉脊山断裂带

中图分类号： P315.727文献标志码：A文章编号： 1000-0844（2024）04-0942-08

DOI：10.20000/j.1000-0844.20240101001

Stress change on surrounding faults and triggering of aftershocks

by the 2023 Jishishan MS6.2 earthquake in Gansu， ChinaZHOU Mingyue WAN Yongge GUAN Zhaoxuan JIN Zhitong

（1. Institute of Disaster Prevention，Sanhe 065201， Hebei， China;

2.Hebei Key Laboratory of Earthquake Dynamics， Sanhe 065201， Hebei， China）Abstract：

On December 18， 2023， an MS6.2 earthquake occurred in Jishishan County， Linxia Prefecture， Gansu Province， China. The earthquake is located on the Lajishan fault. A comprehensive analysis of the impact of this earthquake on the static Coulomb failure stress （CFS） of surrounding faults and aftershocks was performed. The changes in static CFS at different depths along the surrounding faults due to the Jishishan MS6.2 earthquake were calculated based on the parameters of earthquake rupture and elastic half-space models， and the future potential seismic hazard in the study region was then estimated based on the CFS changes. Results indicate that the event substantially increased the CFS in some areas of the Qinghai Nanshan-Xunhua Nanshan fault and the northern and southern margin faults of Lajishan. Particularly， the eastern segment of the Qinghai Nanshan-Xunhua Nanshan fault experienced a CFS increase far exceeding the static stress triggering threshold， reaching 0.022 MPa， indicating a potentially high seismic hazard. Other faults also displayed various degrees of CFS changes， with stress unloading amounts reaching 10 000 Pa in some areas of the northern and southern margin faults of Lajishan and 100 Pa on the Zhuanglanghe fault and other segments of the Qinghai Nanshan-Xunhua Nanshan fault， with the exception of the eastern segment. Furthermore， a comparative analysis of CFS changes on faults at different depths revealed the considerable influence of depth variation on the CFS pattern of the earthquake. However， the increase in CFS along the Qinghai Nanshan-Xunhua Nanshan fault is crucial. This study provides a reference for evaluating the possibility and risk of future earthquakes in the study region. Moreover， the CFS changes due to the mainshock on the fault planes of some large aftershocks were calculated. The CFS on the fault planes of three subsequent MS≥4 earthquakes were estimated as 0.024， 0.033， and 0.034 MPa， and all exceeded the triggering threshold， displaying the triggering effect of the MS6.2 earthquake on the three events.

Keywords：Jishishan， Gansu MS6.2 earthquake; Coulomb failure stress; seismic hazard; Lajishan fault zone0引言

据中国地震台网中心测定，北京时间2023年12月18日23时59分，甘肃省临夏州积石山县发生MS6.2地震，震中位置为35.70°N ，102.79°E，震源深度10 km。震后发生了多次余震，属于典型的“主震-余震”型地震，截至2023年12月30日，共记录到3次MS≥4的地震。

该地震序列位于青藏高原东北缘的柴达木—祁连地块内部的拉脊山断裂带上。柴达木—祁连地块由阿尔金断裂带、祁连山—海原断裂带和东昆仑断裂带三条巨型左旋走滑断裂所围限［1］，内部孕育着多条断裂带，其中拉脊山断裂带是位于该地块内部的一条由拉脊山北缘断裂和拉脊山南缘断裂组成的大型断裂带［2］。李智敏等［3］研究表明该区域在史前存在着非常强烈的地震活动，这些强烈的古地震事件导致了喇家文化的毁灭。震源周围发育有青海南山—循化南山断裂、日月山断裂和西秦岭北缘断裂等断裂（图1）。因此，研究本次地震对周围断层的影响、及其对后续地震的触发关系，对评估未来可能的发震趋势、地震危险性，以及该地区的地震动力学研究具有重要意义。

前人研究表明，通过计算地震在周围断层上产生的库仑破裂应力变化可以评估地震对周围断层的影响。万永革等［4］计算出2008年汶川MW7.9地震造成的周围断层上的静态库仑破裂应力变化，发现多个断裂带均有明显的库仑破裂应力增加。盛书中等［5］计算了2015年尼泊尔MS8.1地震对中国大陆地区的静态应力触发情况，发现尼泊尔地震在中国大陆产生的应力变化量值很小。靳志同等［6］计算出2017年九寨沟MS7.0地震在周围主要活动断层上产生的库仑应力变化，发现九寨沟地震导致虎牙断裂中段的库仑破裂应力超过了应力触发阈值。冯淦等［7］分析2021年玛多MS7.4地震在震中附近主要断层上产生的同震库仑应力变化，探讨了玛多—甘德断裂和东昆仑断裂的后续地震危险性。许鑫等［8］通过剖析2021年漾濞MS6.4地震对周围断层的库仑破裂应力影响，发现漾濞地震对周围断裂带影响较小，均未达到应力触发阈值。这些研究通过计算地震在周围断层上产生的库仑破裂应力，为后续地震预测研究提供了数据基础。

黑色线条为断层；图中左下角给出了研究区域的位置，其中红色线条为块体边界，红色方框为研究区域的位置为深度分析积石山MS6.2地震对周围断层的影响，本研究将此次地震的破裂模型参数和Okada［9］给出的弹性半空间模型位错在周围产生应变的解析表达式作为数据资料和方法，首先，计算此次地震在不同深度下对周围主要断层产生的静态库仑破裂应力变化，并依据计算结果和相关地质资料对周围主要断层的地震危险性进行分析；然后，计算此次地震在后续较大余震断层面上产生的库仑破裂应力，以判断主震和较大余震的触发关系。

1库仑破裂应力计算方法

地震是由地球内部构造运动导致地壳中的应力积累超过岩石所能承受的上限而突然发生错动和破裂的自然现象，这种错动和破裂释放大量应变能导致周围地壳突然变形。地球介质可以一阶近似为半无限空间内均匀各向同性的完全弹性体，因此，通过地震破裂模型参数可计算出弹性体内部产生的应变场。本文利用北京大学张勇研究组提供的积石山MS6.2地震的断层破裂模型（https：//pku-geophysics-source.group），通过Okada给出的地震位错在弹性半空间模型中产生的应变场分量的解析表达式，计算积石山MS6.2地震产生的库仑破裂应力。断层面上库仑破裂应力变化可以由剪应力、正应力和视摩擦系数［10］表示：

ΔCFS=Δτ+μ′Δσ_n （1）

式中：ΔCFS为在断层上产生的库仑破裂应力变化［11-12］；Δτ为断层面滑动方向上的剪切应力变化;Δσn为断层面上的正应力变化;μ′为摩擦系数，一般取0.2～0.8［13-14］，本文仿照前人研究［15-16］，取μ′=0.4。ΔCFS=0时，表示地震没有对断层产生影响;0.01 MPagt;ΔCFSgt;0时，表明断层上受到的应力增加，但并未到达对后续地震具有明显触发作用的阈值（0.01 MPa）;ΔCFS≥0.01 MPa时，则对断层的破裂具有明显的促进作用;ΔCFSlt;0时，则抑制断层的破裂。

2主要断裂上的同震库仑应力变化

本研究首先从文献［3-4，17-19］中收集了此次地震震源区周围断裂的断层滑动参数（表1），利用张勇提供的破裂模型，计算甘肃积石山MS6.2地震在震源深度处周围主要断层面上产生的库仑破裂应力变化（表1、图2）。由于中国地震台网中心给出的震源深度为10 km，因此本研究选择10 km作为震源深度。

由图2所示，本次地震主要造成青海南山—循化南山断裂、拉脊山北缘断裂和拉脊山南缘断裂部分分段的库仑破裂应力增加。其中，青海南山—循化南山断裂东段是本次地震库仑应力增加最显著的区域，达到了0.022 MPa，超过静态应力触发阈值（0.01 MPa）;另外，拉脊山北缘断裂最大库仑破裂应力增加量为0.006 MPa，略低于静态应力触发阈值;拉脊山南缘断裂最大库仑破裂应力增加量为0.004 MPa，低于静态应力触发阈值。

此次地震也导致拉脊山北缘断裂、拉脊山南缘断裂部分区域的库仑破裂应力明显降低。其中，拉脊山北缘断裂的库仑破裂应力最大卸载量达0.028 MPa，拉脊山南缘断裂的库仑破裂应力最大卸载量达0.033 MPa。

此次地震对金城关断裂、庄浪河断裂、西秦岭北缘断裂、临潭—宕昌断裂、马衔山断裂和日月山断裂的影响较小。上述断裂库仑破裂应力的最大卸载量和最大增加量均为百帕，远低于触发的阈值。

在地震研究中，精确的震源深度很难确定。对积石山MS6.2地震，不同地震监测机构提供了不同的震源深度数据，中国地震台网中心给出的数据是10 km，全球矩心矩张量（GCMT）和德国地球科学研究中心（GFZ）提供的震源深度为19 km，而法属波利尼西亚探测与地球物理实验室（CPPT）给出的震源深度为16 km。为了探究深度对库仑破裂应力的影响，本文计算了在16 km和19 km深度下积石山MS6.2地震震源区主要断裂上的同震库仑应力变化情况。

研究结果显示（图3），改变震源深度对库仑破裂应力变化的总体模式产生了显著影响，与倾滑型地震产生的库仑破裂应力存在较大变化的现象一致。拉脊山北缘断裂和拉脊山南缘断裂的库仑破裂应力变化在不同深度下有很大区别，呈现出不同的变化模式：16 km下拉脊山北缘断裂南段最大库仑破裂应力达到了0.019 MPa，但在相同分段处，10，19 km下库仑破裂应力均为卸载状态；19 km下拉脊山南缘断裂南段库仑破裂应力变化模式和10，16 km下的结果不同，10和16 km下是产生库仑破裂应力卸载的分段，在19 km下表现为库仑破裂应力增加，增加量为0.004 MPa。但在不同深度下，青海南山—循化南山断裂东段的库仑破裂应力仍然明显增加，16、19 km下该断裂东段的库仑破裂应力最大增加量分别为0.014 MPa和0.012 MPa（图3），仍然超过了静态应力触发阈值。

此次地震的震中位置位于拉脊山北缘断裂和拉脊山南缘断裂之间，推测其发震断裂为拉脊山北缘断裂或拉脊山南缘断裂。当地震发生在发震断层上时断裂的应力会得到释放，但是，在16 km下拉脊山北缘断裂距离震中最近分段的库仑破裂应力增加量超过了静态应力触发阈值，表现为应力未得到释放的形式，在19 km下拉脊山南缘断裂也处于库仑破裂应力增加的区域。因此，本研究考虑了这些不同的数据来源，选择中国地震台网中心提供的震源深度（10 km），来计算本次地震对周围主要断层的库仑破裂应力变化是合理且可靠的。

为验证其他深度下青海南山—循化南山断裂东段的库仑破裂应力变化也超过了应力触发阈值，还计算了0～45 km范围下该断层的库仑破裂应力变化（表2）。如表2所列，0～33 km深度下的库仑破裂应力变化均超过应力触发阈值，但在33 km以上深度未超过应力触发阈值;在6～12 km深度下的库仑破裂应力变化远超应力触发阈值，这表明在6～12 km深度下青海南山—循化南山断裂东段的地震危险性最高。

根据库仑破裂应力计算结果可以看出，库仑破裂应力变化较大的区域位于青海南山—循化南山断裂东段，此处的库仑破裂应力得到了很大的积累，在不同深度下库仑破裂应力增加量均超过了静态应力触发阈值。青海南山—循化南山断裂近年来并未记录到明显的强震事件，但在本次地震后其库仑破裂应力变化较为显著，特别是部分区域的库仑破裂应力远远超过了静态应力触发阈值，这表明该断裂处于潜在的应力积累高值区域，其地震危险性值得关注。这种程度的库仑破裂应力增加可能会成为未来地震发生的重要因素。

拉脊山北缘和南缘断裂有部分分段区域的库仑应力增加量达到了千帕量级，低于静态应力触发阈值，但阈值只是前人根据经验给出的临界触发值。Ziv等［20］通过计算地震受到的累积静态应力变化，发现当库仑应力变化小于0.01 MPa时，也可能会导致区域累积应力变化的增加，进而影响地区的地震活动性。因此，该结果对该地区的危险性分析具有一定的参考意义。

本文在使用一级近似来估算地震对周围断层活动影响时，采用的是比较简单的弹性半空间位错模型。在实际情况中，库仑破裂应力变化还会受到地球内部的三维不均匀结构和地球的黏弹性松弛效应的影响[21-23]，但这种影响在震后数年时间内是可以忽略不计的［4］，因此，本研究只采用弹性半空间错位模型计算本次地震对断层的影响是合理的。地震发生后所产生的大量余震可能会对周围断层产生一定的影响［24］，但是本次地震产生的余震震级较小，最大余震震级为MS4.1。因此，余震所产生的应力变化并不会改变本研究结果的空间模式，本研究未对其进行考虑，而是在一级近似下计算了积石山MS6.2地震对周围主要断层的影响。

3对后续地震的触发作用

通过计算地震在后续地震断层面上的库仑破裂应力，可以揭示地震之间的触发关系［25］。为了讨论积石山MS6.2地震对其后三次较大震级地震触发作用，搜集了中国地震台网中心［26］和郭祥云（私人通讯）提供的震后三次MS≥4地震的震源机制解。由于使用的方法和数据不同，这些震源机制具有一定的离散度，为地震动力学分析或其他应用带来抉择的困难。因此，本研究整理了中国地震台网中心和郭祥云给出的震源机制解，根据万永革［27］ 提出的方法，求出与所有测定的震源机制的差别平方和最小的一个解，将其作为中心震源机制解，并将中心震源机制解中与断裂带走向大体一致的节面作为接收断层面（表3）。利用张勇提供的破裂模型，计算出18日MS6.2地震在19日0时36分MS4.0地震断层面的滑动方向上产生的剪切应力变化为0.022 MPa，正应力变化为-0.065 MPa，库仑破裂应力变化为0.024 MPa;计算出18日MS6.2地震在19日0时59分MS4.1地震断层面的滑动方向上产生的剪切应力变化为0.240 MPa，正应力变化为-0.006 MPa，库仑破裂应力变化为0.033 MPa;计算出18日MS6.2地震在21日4时2分MS4.1地震断层面的滑动方向上产生的剪切应力变化为0.228 MPa，正应力变化为-0.056 MPa，库仑破裂应力变化为0.034 MPa。

主震产生的应力张量在这三次地震的破裂面的滑动方向上产生的库仑破裂应力变化在主震破裂面附近分别为0.024 MPa、0.033 MPa和0.034 MPa，均超过了0.01 MPa的应力触发阈值，这表明积石山MS6.2地震发生后，主震产生的破裂促进了震后三次MS≥4地震的发生。

4结论

本研究基于甘肃省积石山县MS6.2地震的震源模型参数和弹性半空间模型，研究了甘肃积石山MS6.2地震对周围断层的影响，以及对后续地震的触发作用。

通过计算积石山MS6.2地震产生的库仑破裂应力，研究发现青海南山—循化南山断裂东段是本次地震最显著的影响区域，其库仑破裂应力增加量超过了静态应力触发阈值，其潜在的地震危险性值得高度关注。拉脊山北缘断裂和拉脊山南缘断裂等部分区域的库仑破裂应力也发生了显著变化，其最大库仑破裂应力增加量均为千帕量级，虽然库仑破裂应力略低于触发阈值，但考虑到库仑应力变化的累积效应，这样的结果仍然会对该地区的地震活动性产生影响。

针对震源深度的不确定性，本研究进行了不同深度下库仑破裂应力计算，结果显示青海南山—循化南山断裂东段的库仑破裂应力增加依旧显著，震源深度为0～33 km时，库仑破裂应力均超过了静态应力触发阈值。对于拉脊山北缘、南缘断裂，库仑破裂应力变化在不同深度下呈现出相异的模式。

最后，本文还讨论了积石山MS6.2地震对后续地震的触发作用，利用后续三次较大震级地震的中心震源机制解，计算了积石山MS6.2地震在这些地震断层面上产生的库仑破裂应力变化。结果显示，在这三次地震破裂面的滑动方向上，库仑破裂应力变化均超过了0.01 MP的阈值，说明积石山MS6.2地震对这些地震具有明显的促进作用。

致谢：感谢北京大学张勇研究组提供了甘肃积石山MS6.2地震的破裂模型。本文图形采用GMT软件［28］和MATLAB软件绘制，特此致谢。

参考文献（References）

［1］葛伟鹏，王敏，沈正康，等.柴达木—祁连山地块内部震间上地壳块体运动特征与变形模式研究［J］.地球物理学报，2013，56（9）：2994-3010.GE Weipeng，WANG Min，SHEN Zhengkang，et al.Intersiesmic kinematics and defromation patterns on the upper crust of Qaidam—Qilianshan Block［J］.Chinese Journal of Geophysics，2013，56（9）：2994-3010.

［2］袁道阳，张培震，雷中生，等.青海拉脊山断裂带新活动特征的初步研究［J］.中国地震，2005，21（1）：93-102.YUAN Daoyang，ZHANG Peizhen，LEI Zhongsheng，et al.A preliminary study on the new activity features of the Lajishan Mountain fault zone in Qinghai Province［J］.Earthquake Research In China，2005，21（1）：93-102.

［3］李智敏，李延京，田勤俭，等.拉脊山断裂古地震与喇家遗址灾变事件关系研究［J］.地震研究，2014，37（增刊1）：109-115.LI Zhimin，LI Yanjing，TIAN Qinjian，et al.Study on the relationship between the ancient earthquake of Lajishan fault and the catastrophic event of Lajia site［J］.Journal of Seismological Research，2014，37（Suppl01）：109-115.

［4］万永革，沈正康，盛书中，等.2008年汶川大地震对周围断层的影响［J］.地震学报，2009，31（2）：128-139.WAN Yongge，SHEN Zhengkang，SHENG Shuzhong，et al.The influence of 2008 Wenchuan earthquake on surrounding faults［J］.Acta Seismologica Sinica，2009，31（2）：128-139.

［5］盛书中，万永革，蒋长胜，等.2015年尼泊尔MS8.1强震对中国大陆静态应力触发影响的初探［J］.地球物理学报，2015，58（5）：1834-1842.SHENG Shuzhong，WAN Yongge，JIANG Changsheng，et al.Preliminary study on the static stress triggering effects on China Mainland with the 2015 Nepal MS8.1 earthquake［J］.Chinese Journal of Geophysics，2015，58（5）：1834-1842.

［6］靳志同，万永革，刘兆才，等.2017年九寨沟MS7.0地震对周围地区的静态应力影响［J］.地球物理学报，2019，62（4）：1282-1299.JIN Zhitong，WAN Yongge，LIU Zhaocai，et al.The static stress triggering influences of the 2017 MS7.0 Jiuzhaigou earthquake on neighboring areas［J］.Chinese Journal of Geophysics，2019，62（4）：1282-1299.

［7］冯淦，万永革，许鑫，等.2021年青海玛多MS7.4地震对周围地区的应力影响［J］.地球物理学报，2021，64（12）：4562-4571.FENG Gan，WAN Yongge，XU Xin，et al.Static stress influence of the 2021 MS7.4 Madoi，Qinghai earthquake on neighboring areas［J］.Chinese Journal of Geophysics，2021，64（12）：4562-4571.

［8］许鑫，万永革，李振月，等.2021年云南漾濞MS6.4地震序列静态应力触发研究［J］.地震工程学报，2022，44（4）：970-979.XU Xin，WAN Yongge，LI Zhenyue，et al.Static stress triggering of the 2021 Yangbi MS6.4 earthquake sequence in Yunnan Province［J］.China Earthquake Engineering Journal，2022，44（4）：970-979.

［9］OKADA Y.Internal deformation due to shear and tensile faults in a half-space［J］.The Bulletin of the Seismological Society of America，1992，82（2）：1018-1040.

［10］AMJAD A，SHAHID R M，SHAKIR U，et al.Evolution of Coulomb failure stress in Sulaiman Lobe，Pakistan，and its implications for seismic hazard assessment［J］.Earth and Planetary Physics，2023，7（2）：282-294.

［11］HARRIS R A.Introduction to special section：stress triggers，stress shadows，and implications for seismic hazard［J］.Journal of Geophysical Research：Solid Earth，1998，103（B10）：24347-24358.

［12］FREED A M.Earthquake triggering by static，dynamic，and postseismic stress transfer［J］.Annual Review of Earth and Planetary Sciences，2005，33：335-367.

［13］STEIN R S.The role of stress transfer in earthquake occurrence［J］.Nature，1999，402：605-609.

［14］COTTON F，COUTANT O.Dynamic stress variations due to shear faults in a plane-layered medium［J］.Geophysical Journal International，1997，128（3）：676-688.

［15］KING G C P，STEIN R S，LIN U J.Static stress changes and the triggering of earthquakes［J］.Bulletin of the Seismological Society of America，1994，84（3）：935-953.

［16］万永革，沈正康，盛书中，等.2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释［J］.地球物理学报，2010，53（2）：280-289.WAN Yongge，SHEN Zhengkang，SHENG Shuzhong，et al.The mechanical effects of the 2008 MS7.3 Yutian，Xinjiang earthquake on the neighboring faults and its tectonic origin of normal faulting mechanism［J］.Chinese Journal of Geophysics，2010，53（2）：280-289.

［17］曾宪伟，冯建刚，王晓山.利用小震分布和区域应力场确定临潭—宕昌断裂断层面参数［J］.地震工程学报，2013，35（4）：743-750.ZENG Xianwei，FENG Jiangang，WANG Xiaoshan.Determination of fault plane parameter of lintan-tanchang fault using distribution of small earthquakes and regional stress field［J］.China Earthquake Engineering Journal，2013，35（4）：743-750.

［18］侯珍清，王多杰，董志平.兰州金城关断裂的新活动性与地震［J］.西北地震学报，1985，7（2）：95-98.

HOU Zhenqing，WANG Duojie，DONG Zhiping.Aactivities and earthquakes of Lanzhou Jinchenguan fault since Holocene epoch［J］.Northwestern Seismological Journal，1985，7（2）：95-98.

［19］宋方敏，袁道阳，陈桂华，等.甘肃马衔山北缘断裂西北段几何结构及其新活动［J］.地震地质，2006，28（4）：547-560.SONG Fangmin，YUAN Daoyang，CHEN Guihua，et al.Geometric structures and recent acticvity along the northwest segment of north marginal fault of Maxianshan Mountains，Gansu Province［J］.Seismology and Geology，2006，28（4）：547-560.

［20］ZIV A，RUBIN A M.Static stress transfer and earthquake triggering：No lower threshold in sight？［J］.Journal of Geophysical Research：Solid Earth，2000，105（B6）：13631-13642.

［21］万永革，沈正康，曾跃华，等.青藏高原东北部的库仑应力积累演化对大地震发生的影响［J］.地震学报，2007，29（2）：115-129.WAN Yongge，SHEN Zhengkang，ZENG Yuehua，et al.Evolution of cumulative coulomb failure stress in northeastern Qinghai—Xizang （Tibetan） Plateau and its effect on large earthquake occurrence［J］.Acta Seismologica Sinica，2007，29（2）：115-129.

［22］万永革，沈正康，曾跃华，等.唐山地震序列应力触发的粘弹性力学模型研究［J］.地震学报，2008，30（6）：581-593.WAN Yongge，SHEN Zhengkang，ZENG Yuehua，et al.Study on visco-elastic stress triggering model of the 1976 Tangshan earthquake sequence［J］.Acta Seismologica Sinica，2008，30（6）：581-593.

［23］DENG J，GURNIS M，KANAMORI H，et al.Viscoelastic flow in the lower crust after the 1992 Landers，California，earthquake［J］.Science，1998，282（5394）：1689-1692.

［24］万永革，沈正康，兰从欣.兰德斯地震断层面及其附近余震产生的位移场研究［J］.地震学报，2005，27（2）：139-146.WAN Yongge，SHEN Zhengkang，LAN Congxin.Study on displacement field generated by aftershocks in landers seismic fault plane and its adjacent areas［J］.Acta Seismologica Sinica，2005，27（2）：139-146.

［25］关兆萱，万永革，黄少华，等.2022年墨西哥米却肯州7.6级地震序列的触发关系及发震动力学探讨［J］.中国地震，2023，39（3）：584-595.GUAN Zhaoxuan，WAN Yongge，HUANG Shaohua，et al.Study on the triggering mechanism and earthquake dynamics of the 2022 Michoacan MW7.6 earthquake sequence in Mexico［J］.Earthquake Research in China，2023，39（3）：584-595.

［26］杨志高，徐泰然，梁建宏等.中国地震台网中心大震震源机制CMT产品［EB/OL］.[2024-01-01].https：//https：//data.earthquake.cn/datashare/report.shtml？PAGEID=earthquake_dzzyjzamp;cmtype=review.YANG Zhigao，XU Tairan，LIANG Jianhong et al.China Earthquake Networks Center （CENC） moment tensor （CMT）product［EB/OL］.[2024-01-01].https：//https：//data.earthquake.cn/datashare/report.shtml？PAGEID=earthquake_dzzyjzamp;cmtype=review.

［27］万永革.同一地震多个震源机制中心解的确定［J］.地球物理学报，2019，62（12）：4718-4728.WAN Yongge.Determination of center of several focal mechanisms of the same earthquake［J］.Chinese Journal of Geophysics，2019，62（12）：4718-4728.

［28］WESSEL P，LUIS J F，UIEDA L，et al.The generic mapping tools version 6［J］.Geochemistry，Geophysics，Geosystems，2019，20（11）：5556-5564.

（本文编辑：贾源源）