引言

据中国地震台网中心（2022）测定，2022年1月8日1时45分27秒在青海省海北州门源县发生M_S6.9地震，其震中位于（101.26°E，37.77°N），震源深度为10 km。此次地震发生在青藏高原东北缘，处于欧亚板块与印度板块碰撞的远程效应区域。该区域的地质构造复杂，地表高程变化剧烈，区域地壳厚度由高原内部的65 km下降到边缘的55 km （Pan，Niu，2011；Deng et al，2015；Wang et al，2016）。区域断裂带多呈现NE−SE向的走滑断层和逆冲断层，以适应印度板块向高原的俯冲作用。

这次门源地震发生在祁连—海原断裂带上，该断裂是一个可以孕育M8强震的断裂带，在二十世纪发生过两次M8地震，分别是1920年海原大地震（郭增建，吴瑾冰，2000）和1927年古浪大地震（邹谨敞等，1994）。GPS观测显示海原断裂带的左旋速率为4—5 mm/a （Zheng et al，2017）。自二十世纪七十年代有地震仪器记录以来，本次门源M_S6.9地震附近50 km范围内共发生M≥5.0地震9次，其中M≥6.0地震3次，分别是1986年门源M_S6.4地震、2016年门源M_S6.4地震和此次地震。中国地震局（2022）组织的这次地震现场调查显示，Ⅷ度烈度区域为1143 km²，地震造成的地表破裂长约22 km，地表破裂性质以左旋走滑为主。

快速、可靠地获得强震的破裂过程对于震后灾害评估和救灾工作具有重要的指导意义，同时也是震源物理过程研究的重要内容。由于长周期远震波形数据反演破裂过程的时空分辨能力有限，只能给出M_W7.0以上地震的大致破裂模式，缺乏对破裂细节的约束，而高频近场强震观测波形数据可以更好地约束破裂过程的细节，因而能提供更多的地震破裂细节。为此，本文拟利用新布设的和改建的强震台网记录到的波形数据，基于迭代反褶积和叠加法（iterative deconvolution and stacking，缩写为IDS）（Zhang et al，2014）反演这次门源M_S6.9地震的破裂过程，并通过与其它机构的结果对比来展示本文破裂过程在反演的时效性和破裂过程细节方面的优势。

1.   方法与数据处理

本文采用IDS方法（Zhang et al，2014）基于强震数据进行破裂过程反演。IDS方法融合了经典的波形反演和反投影方法，基于近场观测数据，可以实现快速、稳定以及高度流程自动化的地震破裂过程反演。特别是利用近年来发展起来的强震观测台网和高频GPS观测台网的实时波形数据，地震破裂过程反演的时效性显著增强。该方法目前已经比较成熟，且在2014年南纳帕（South Napa）M_W6.1地震（Zhang et al，2015）和2019年九寨沟M_S7.0地震（郑绪君等，2017）的应用中都取得了较好的结果。

随着国家地震烈度速报与预警工程的实施，青海地区相继建成并完善了实时传输的强震台网，为本研究提供了数据基础。强震观测的主要目的是获取强震的近场观测记录，通常布设在活动断裂和人口稠密区域。这次门源地震震中150 km范围内有强震台34个，台站分布情况如图1所示。可见强震台主要分布在地震的南部，东北部稀缺。这样的台站分布可选用IDS方法来反演破裂过程。这是因为以往的震例显示IDS方法适用于地震台网方位角覆盖不完全的区域，例如2008—2017年期间日本东部海域M_W6.5以上地震的强震台站方位角覆盖不足180°，利用该方法也得到了稳定的破裂过程（Zheng et al，2020）。强震台网配备的地震仪器包括三分量TDA-33M型、JS-A2型和QA-2g型地震计，采样率为100 Hz，其对地面加速度的平坦响应频带范围分别是DC−200 Hz，DC−100 Hz和DC−170 Hz。

图  1  2022门源M_S6.9地震周边强震台站分布

紫色三角形表示最终参与反演的强震台站

Figure  1.  Distribution of strong motion stations near 2022 Menyuan M_S6.9 earthquake

Purple triangles are the strong motion stations used for inversion

下载: 全尺寸图片 幻灯片

门源M_S6.9地震的波形可通过中国地震台网中心布设的实时地震波形数据服务器获取，震源位置和发震时刻采用地震速报结果。由于使用大量近台数据，地震速报结果是可靠的。断层面参数选用中国地震台网中心的震源机制解（杨志高等，2022），该结果基于国内甚宽频带地震台和全球台网（global seismic network，缩写为 GSN）台站，震后20 min由 W 震相反演方法（Kanamori，Rivera，2008）获得。反演震源机制解的节面Ⅰ为走向102°、倾角69°、滑动角−2°，节面Ⅱ为走向193°、倾角89°、滑动角−159°，矩震级为M_W6.6，其中节面Ⅰ的走向与祁连—海原断裂带以及余震序列的长轴方位较为接近，因此选用节面Ⅰ参数作为震源破裂过程的断层面参数。破裂面尺度基于震级与破裂长宽之间的经验关系（Wells，Coppersmith，1994）计算得到，长度和宽度分别为34 km和14 km。子断层网格大小设为2 km×2 km，以此划分出沿走向和沿倾向的17×7个网格。格林函数使用Qseis （Wang，1999）程序，基于震源所处位置的crust1.0模型（Laske et al，2012）进行计算。地震数据筛选主要根据震中距、波形信噪比以及波形拟合度进行。具体的筛选步骤为：选取距震中150 km以内的台站；剔除存在明显漂移现象的波形数据，以避免积分后波形出现较大畸变；多次迭代过程中剔除拟合度较低（＜0.4）的台站。进行10次迭代运算后筛选出10个波形拟合度较高的台站（图1中紫色三角形所示），其观测波形与理论波形在振幅和形状上有较好的一致性（图2）。

图  2  2022年门源M_S6.9地震的理论波形与观测波形对比

台站名称标示在波形图的左上；黑线表示各个强震台的三分量观测波形，自上而下分别是EW，NS和UD分量的波形，红线表示基于破裂过程的理论地震图，观测地震图与理论地震图的拟合度标示在各分量波形左下

Figure  2.  Comparison of theoretical and observed waveforms of the 2022 Menyuan M_S6.9 earthquake

Station codes are marked on the upper left of waveforms. The black lines represent the three-component observation waveform observed by strong motion station. Each panel represents the east-west，north-south，and vertical components from top to bottom，respectively。The red lines represent the theoretical seismograms based on the inverted rupture process，and the cross correlation between the observed and the theoretical waveform is marked on the lower left of each waveform

下载: 全尺寸图片 幻灯片

2.   结果

图3给出断层面静态滑动量分布，可见：整体上地震的破裂以单侧为主，破裂主要自初始破裂点向ESE方位扩展，破裂中心位于震中偏ESE向约2—4 km；破裂中心的最大滑动量为3.6 m，深度约4 km；破裂整体长约20 km，向WNW向扩展约6 km，向SE向扩展约14 km。破裂过程纵向上展示出由深部向浅部扩展的特征，且浅部可见明显的滑动量。该静态滑动分布与基于大地测量方法的结果（Li et al，2022；Yang et al，2022）一致。根据震后现场考察结果（中国地震局，2022），地表可见22 km长的破裂带，烈度为Ⅸ度的区域面积约157 km²，烈度为Ⅷ度的区域为1143 km²。根据震级与Ⅷ度区域面积的经验关系（陈培善，刘家森，1975），当矩震级为M_W6.6时，Ⅷ度区域的面积为360 km²。这次地震引起的实际烈度Ⅷ度区域面积远大于经验值，可能是由于破裂由深部向浅部扩散并导致地表破裂带引起。震源时间函数主要包括一次持续8 s的主破裂和一次持续4 s的次级破裂。地震破裂过程持续时间约为14 s （图4），略大于M_W6.6地震的通常破裂持续时间11 s （Duputel et al，2012）。结合时空破裂过程结果（图5）可以看出：破裂速率在3—6 s达到峰值，主要集中在震中附近且接近地表；9 s之后滑动累积量趋于稳定，整个破裂中8 s前的主破裂事件占主导。

图  3  2022年门源M_S6.9地震断层面静态滑动分布（星形表示主震位置）

Figure  3.  Static slip distribution on the fault plane of the M_S6.9 Menyuan earthquake in 2022 （The red star indicates the epicenter of the main shock）

下载: 全尺寸图片 幻灯片

图  4  2022年门源M_S6.9地震震源时间函数

Figure  4.  Source time function of the M_S6.9 Menyuan earthquake in 2022

下载: 全尺寸图片 幻灯片

图  5  2022年门源M_S6.9地震时空破裂过程

（a） 滑动速率分布；（b） 累积滑动量分布

Figure  5.  The spatio-temporal rupture process of the 2022 Menyuan M_S6.9 earthquakee

（a） Distribution of slip rate；（b） Distribution of cumulative slip

下载: 全尺寸图片 幻灯片

相比过去触发式的强震台站存在较多到时不准、波形异常和仪器灵敏度更新不及时等问题，实时传输的波形数据弥补了触发台站的不足，可以快速检测到这些异常并实时进行更新和修正。通过强震波形反演破裂过程得到的矩震级为M_W6.6，这与利用长周期远震数据的W震相反演方法得到的结果一致，说明强震台仪器的灵敏度等参数配置妥当。

USGS （2022）基于长周期远震波形拟合反演的门源地震破裂过程如图6所示，其断层面模型的走向和倾向分别为289°和84°。对于高倾角走滑型地震而言，两个断层面比较接近，前者向南倾斜，后者向北倾斜。Zhang等（2015）对IDS方法的测试表明断层面参数在±15°范围内的变化不会显著影响主要破裂过程，我们使用USGS的断层面参数（289°和84°）计算后得到的破裂过程与本文参数（102°和69°）的结果几乎无差别，二者最大的不同之处在于破裂中心的位置。USGS结果显示破裂中心位于初始破裂偏西部，而本文结果显示中心位置在初始破裂以东。再者，远震波形拟合与强震波形拟合的破裂过程在空间分辨率上存在差别。远震波形拟合破裂过程的算法通常允许观测波形时间窗前后偏移以搜索最佳波形拟合度，但对破裂绝对位置的分辨能力弱；强震资料反演过程不允许波形窗偏移，有一定的破裂绝对位置定位能力（Zhang et al，2014）。

图  6  2022年门源M_S6.9地震静态滑动分布与震源时间函数（USGS，2022）

右图中红色虚线对应推定的地震破裂结束时间

Figure  6.  Static slip distribution and source time function of the 2022 M_S6.9 Menyuan earthquake （USGS，2022）

In the right panel the red dashed line represents the interpreted end of the event

下载: 全尺寸图片 幻灯片

时效性是强震数据反演破裂过程的优势之一。强震台网相较于GSN能够更早地接收到波形数据。此次反演使用距震中150 km范围内的台站，波形窗长度为50 s，震后2 min即可获得所需时间窗完整的波形数据。在准备好格林函数的情况下，以历史强震震源机制解作为断层面参数则可在震后3—5 min产出初步强震破裂过程。考虑到强震的震源机制通常与历史地震的震源机制有较好的一致性，因此初步破裂过程也具备一定的参考价值。对于M_W6.0左右的地震，通常在震后5—10 min内可以自动产出基于近震波形反演的震源机制解，半小时内可以产出基于远震波形（例如W震相方法）反演的震源机制解。以当前地震震源机制解为断层面参数输入时，可以在震后半小时持续更新破裂过程。

主震破裂的高滑动区域与低余震活动区域之间存在相关性（Beroza，1991；Das，Henry，2003；Chousianitis，Konca，2021），前者通常对应于较少的余震活动（Das，Henry，2003）。本次地震破裂过程的高滑动区域对应于较弱的余震活动（图7），这与Das和Henry （2003）报告的震例一致。解释这种现象的一个简单模型是，高滑动区域也是高强度区域，仅在主震期间滑动（Beroza，1991）并将应力完全释放（Chousianitis，Konca，2021）。

图  7  2022年门源M_S6.9地震序列（主震以白色星形表示，余震以黑色圆点表示）与2016年门源M_S6.4地震序列（主震以紫色星形表示，余震以紫色圆点表示）空间分布对比

Figure  7.  Spatial distribution comparison between the 2022 M_S6.9 Menyuan earthquake sequence （the white star stands for mainshock and the black dots indicate aftershocks） and the 2016 M_S6.4 Menyuan earthquake sequence （the purple star stands for mainshock and the purple dots indicate aftershocks）

下载: 全尺寸图片 幻灯片

2016年门源M_S6.4地震是距离本次地震最近的一次强震，两者的关系值得讨论。2022年门源地震为左旋走滑型事件，2016年地震为逆冲型事件，说明这两次地震的发震断层性质不同，而且发震断层并不完全相同：2016年地震是由冷龙岭断层北侧的断层引起（胡朝忠等，2016）。根据震级与破裂长度的经验关系，2016年门源M_S6.4地震的破裂长度为10 km左右。两次地震的位置相差近40 km，可以大体判断两次地震的破裂过程没有连通。截至北京时间2022年1月18日10时，已发生余震近千次且其空间分布展示出分段特征（图7）。震中附近余震空间分布的长轴方向与震源机制解节面吻合，而震中东南部余震的长轴方向与震源机制解节面走向存在一定角度。据现场考察此次地表破裂由两条破裂带组成，分别为NW向冷龙岭断裂西段和EW向托勒山断裂东段（Li et al，2022；Liu et al，2019），这两个断裂带位置与两次地震的余震分布带吻合。余震空间分布的分段性和历史震源机制解差异预示该区域存在强烈的构造转换，且本次地震破裂止于该处。赵凌强等（2022）指出冷龙岭断裂在南北两侧呈现重要物性差异带的特点，也从侧面指示了这一区域构造的转换性。由于这种构造差异性以及挤压应力的持续作用，该地区未来仍有发生强震的可能性。

3.   结论

基于实时强震连续波形数据和W震相反演的震源机制解，本文利用迭代反褶积和叠加法反演了2022年1月8日门源M_S6.9地震的破裂过程。破裂过程整体上表现为单侧为主的特征，主要破裂方向是ESE，破裂中心位于震中ESE向2—4 km处。破裂在纵向上自下而上达到地表，这与观测到的约22 km的地表破裂吻合。破裂持续时间约14 s，包含一次持续8 s的主破裂和一次持续4 s的次级破裂。破裂速率在3—6 s达到峰值，破裂集中在2—8 s。破裂过程展示的高滑动区域对应较少的余震活动。相较USGS （2022）基于远震波形拟合方法反演所得的破裂过程，本文结果在空间分辨率和时效性上更有优势。本次M_S6.9地震与2016年门源M_S6.4地震在空间上相距40 km，推测两次地震的破裂过程没有连通，且二者的发震断层不同。此次地震的余震空间分布的分段性和历史震源机制解差异预示发震区域存在强烈的构造转换，该地区未来仍然有发生强震的可能性。