2022年门源MS 6.9地震序列特征及强余震判定

余 娜 张晓清 胡维云

（中国西宁 810000 青海省地震局）

0 引言

2022 年1 月8 日1 时45 分青海海北州门源县发生MS6.9 地震（37.77°N，101.26°E），震源深度10 km，震后发生了大量余震并于1 月12 日发生MS5.2 强余震。根据蒋海昆等（2015）地震序列类型定义，门源MS6.9 地震序列为主—余型地震。国内外在余震预测方面已开展了一些研究工作，全球 “地震可预测性合作研究”（CSEP）计划是在全球不同地区开展可检验的地震短期余震预测模型（Jordan，2006），且认为“传染型余震序列”模型（Ogata，1988，1989，2001）是目前最接近真实情况的模型之一，基于该模型参数以及“瘦化算法”（Lewisetal，1979；Ogata，1981）可进行余震概率预测。蒋长胜等（2015）基于ETAS 模型参数进行“瘦化算法”对2014 年云南鲁甸6.5 级地震余震进行“向前预测”并取得了较好效果；郭安宁等（1999）采用柯里奥利力效应研究了兰州地区强余震的预测；刘正荣（1979，1984，1986）对地震频度衰减进行分析并给出了最大余震震级；毕金孟等（2017）、余娜等（2020）采用“瘦化算法”分别对晋冀蒙交界地区和青海地区进行余震短期发生概率预测并给出了预测效能评估。本文拟对门源MS6.9 地震序列进行时空特征分析，并基于历史地震类比、序列自身的活动特征给出地震序列类型及强余震活动特征。为了对比分析的一致性，对于ML≥4.0 地震采用MS震级，ML＜4.0 地震采用ML震级。

1 门源MS 6.9 地震序列基本情况

1.1 概况

门源MS6.9 地震是祁连地震带中东段自2016 年1 月21 日门源MS6.4 地震后发生的一次强震，是继2021 年5 月22 日玛多MS7.4 地震后青海地区近半年内发生的第2 个强震。祁连地区监测能力相对较强，门源MS6.9 地震震中50 km 范围内仅有1 个测震台，100 km范围内有5 个测震台（图1）。据青海地震台提供的青海省地震监测能力分布情况和已有研究结果（冯建刚等，2012；韩立波等，2012；余娜等，2020），震区ML≥1.5 地震目录基本完整。据青海地震台网定位结果，门源MS6.9 地震序列从2022 年1 月8 日1 时45 分MS6.9 主震开始，截至2 月8 日，共发生ML≥4.0 地震14 次（表1），最大余震震级为MS5.2，与主震震级差为1.7。

图1 门源MS 6.9 地震序列震中及周边地震台站分布Fig.1 Distribution of seismic stations and the epicenter of MS 6.9 earthquake of Menyuan

表1 门源MS 6.9 地震序列（ML ≥4.0）Table 1 Catalogue of Menyuan MS 6.9 earthquake sequence（ML ≥4.0）

1.2 最小完备震级Mc

最小完备震级Mc指一定空间、时间内100%被检测到的地震事件的最小震级（Rydelek etal，1989；Taylor et al，1990；Wiemer etal，2000）。为确保地震序列的完整性，采用“震级—序号”法（蒋长胜等，2014）讨论门源MS6.9 地震序列的最小完备震级，根据地震发生的先后顺序进行排序，将地震分布密度较大区域的连线作为最小完备震级Mc的时序变化。采用“震级—序号”法的优势是可以避免强震发生后短期内较小余震的“淹没”、余震区甚至更大范围内地震监测能力显著降低的现象（Akaike，1974）。图2 是门源MS6.9 地震序列目录的完整性分析结果，最小完备震级Mc选取为ML1.8，该结果与冯建刚等（2012）、韩立波等（2012）提供的监测能力分布情况相符。

图2 门源MS 6.9 地震序列目录的完整性分析红色五角星为MS 6.9 地震；黑色圆点为用自然边界法获得的地震序列；红色圆点为背景地震Fig.2 Catalogue completeness analysis for the Menyuan MS 6.9 earthquake sequence

2 地震序列时空特征

2.1 空间展布

根据精定位结果，门源MS6.9 地震序列整体上呈西段NWW 走向、东段SE 走向的带状分布。序列西段近NWW 向带状分布的长轴约19 km，短轴约7 km；序列东段呈SE 向带状分布的长轴约15 km，短轴约8 km（图3）。由图3 可见，MS6.9 主震震中位于序列西段的东南部，而随后的2 个5.0 级以上余震都先后发生在主震震中的两侧，MS5.1 余震发生在序列西段中东部，MS5.2 余震发生在序列东段、呈SE 走向展布的最东南端，距主震震中较远。MS6.9 主震震中西侧余震密集，展布较宽，东侧余震较稀疏。1 月10 日前余震震中主要集中在西段和东段的中西部，1 月12 日MS5.2 地震发生后，该区域地震活动较活跃。

图3 门源MS 6.9 地震序列精定位后的震中分布Fig.3 Distribution of precision positioning of the Menyuan earthquake sequence

2.2 时间分布特征

根据上述序列完整性分析可知，门源MS6.9 地震序列的最小完备震级是ML1.8，因此，采用青海地震台网ML≥1.8 地震序列目录来分析序列的时空特征。图4 为由ETAS 模型给出的门源MS6.9 地震序列ML≥1.8 地震的条件强度曲线和M—t图。由图4（a）可知，门源MS6.9 地震序列的余震活动衰减较缓慢，在主震后4 天内序列衰减较平稳，其后由4 级以上余震引起的条件强度曲线的变化较明显。主震后24 min 后发生了第1 个MS5.1 余震，4 天15h 后发生了MS5.2 的最大震级余震，绝大多数ML≥4.0 余震发生在主震当天和震后第4 天，其后以ML≤4.0 小震活动为主，小震活动频度维持在17—42次［图4（b）］。利用ETAS 模型对门源MS6.9 地震序列进行最大似然估算，得出α=1.882 6、p=0.960 8，相较于2016 年门源MS6.4 地震序列参数α=1.402 5、p=1.303 5，发现此次地震激发高阶余震的能力强于2016 年门源MS6.4 地震，序列衰减速率较慢；与蒋海昆等（2006）给出的中国大陆7 级以上地震序列的平均序列参数α=0.509 相比较，此次门源MS6.9 地震激发高阶余震的能力较弱。

图4 ETAS 模型给出的门源MS 6.9 地震序列ML ≥1.8 地震的条件强度曲线和M-t 图（a）地震序号ML ≥1.8 地震的条件强度曲线；（b）M—t 图Fig.4 Temporal variation of the conditional intensity from fitting the ETAS model to the Menyuan MS 6.9 earthquake sequence with cutoff magnitude ML1.8 and M-t plot

3 地震序列类型判定及强余震早期预测

3.1 序列类型判定

强震后地震序列类型判定是震后趋势估计的一项重要内容和基础性工作。根据蒋海昆等（2015）提出的序列类型判定方法，基于历史地震类比、序列自身的活动特征及地震序列参数等对门源MS6.9 地震进行分析。

3.1.1 历史地震类比分析。门源MS6.9 地震发生在祁连山断裂带中东段，1970 年以来在距此次地震震中100 km 范围内5 级以上地震共6 次（表2）。门源MS6.9 地震序列震中附近区域6 组历史地震中，有5 组为主—余（或前—主—余）型序列，1 组为双震型。1 组双震型地震序列2 次主震发震时间仅间隔7 min，截至2022 年2 月8 日，门源MS6.9 地震已发生1 个月，与上述双震型差异很大。类比该构造的历史地震序列，此次地震与1986 年门源MS6.5 地震序列的发展相似，后者在ML4.0 地震平静19 天后有1 次起伏活动，即发生了MS5.2 最大震级余震，之后该地震序列正常衰减。根据历史地震序列类比和最大余震震级差，判定门源MS6.9 地震为主—余型地震序列。

表2 距震中100 km 范围内M ≥5.0 地震序列类型Table 2 Sequence types of earthquakes with M ≥5.0 within 100 km of the epicenter

3.1.2 基于序列演化特征的序列类型判定。地震序列类型判定更多地是从序列本身进行定性或半定量的判定（蒋海昆等，2007），主要方法有h值等具有一定物理意义并且在一定程度上表征序列特征的参数法和等待时间法等。

（1）h值。h值主要基于频度衰减的大森公式进行计算，根据衰减系数h值大于1、等于1 或小于1 来判定地震序列是前震序列还是余震序列（刘正荣等，1979，1984；蒋海昆等，2007）。选取ML1.8 为门源MS6.9 地震序列的最小完备震级来计算h值，计算结果显示，震后2 天h值为1.8，震后3 天为1.75，震后4 天为1.6，震后5—7 天为1.4，其后h值稳定在1.34 左右（图5）。此次地震序列的衰减系数特征与蒋海昆等（2020）给出的中国大陆75 次地震序列特征一致，即强震后前几天h值波动较大，一般4—8 天后趋于稳定。因此判断门源MS6.9 地震序列为主—余型。

图5 门源MS 6.9 地震序列h 值Fig.5 h value of the of Menyuan MS 6.9 earthquake sequence

（2）等待时间法。将连续发生的2 个强余震之间的时间间隔作为后一个强余震的“等待时间”，即第i个强余震的“等待时间”Δti=ti-ti-1，i为序列中第i个强余震的序号，ti为第i个强余震的发震时刻（由主震起算）。利用最小二乘法进行拟合，发现主—余型地震序列随着t的增大，较强余震的Δti也随着增大；对于主—余型地震或双主震型地震在前一个主震发生到后一个主震发生前的一段时间里，lgΔt-lgt的线性关系不明显（谷继成等，1979）。根据门源MS6.9 地震后ML≥4.0 强余震的发震时间及其时间间隔，得到强余震的位置，其在双对数坐标中大体呈线性分布（图6），表明门源MS6.9 地震后ML≥4.0 地震时间间隔正常，序列为主—余型。

图6 门源MS 6.9 地震Δt—t 曲线Fig.6 Δt-t curve of Menyuan MS 6.9 earthquake

3.2 强余震预测

由表2还可见，ΔM波动范围较大，为0.3—2.3。50%的最大震级余震发震在主震后1 天内，33%的为1—7 天，没有超过30 天的震例，且最大震级余震发生在主震当天的概率最大。2022年1月8日门源MS6.9地震后至2月8日，序列发生的最大震级余震为MS5.2，其发生在主震后4 天 15 h，震级差ΔM为1.7。根据历史地震类比和序列自身活动特征分析认为，门源MS6.9 地震的最大震级余震已经发生，即为1 月12 日18 时20 分的MS5.2 地震。

4 结论与讨论

（1）选取截止震级ML1.8 以上的地震进行精定位，发现门源MS6.9 地震的余震序列震中整体上呈西段NWW、东段SE 的带状分布，且主震震中西侧余震密集，展布较宽，东侧余震较稀疏。

（2）采用ETAS 模型计算门源MS6.9 地震序列参数，发现此次地震序列衰减较缓慢，在主震后4 天内序列衰减较平稳，其后由4 级以上余震引起的条件强度曲线的变化较明显。主震后24 min 后发生了第1 个MS5.1 余震，4 天15 h 后发生了MS5.2 的最大震级余震，绝大多数ML≥4.0 余震发生在主震当天和震后第4 天。

（3）基于构造、历史地震类比及序列自身活动特征对门源MS6.9 地震进行分析。结果显示，祁连山地震带中东段5 级以上地震序列类型多以主—余型地为主，只有1 次双震型；门源MS6.9 地震后h值≥1 且相对稳定，表明主震当天小震频度较高；等待时间法结果显示，门源MS6.9 地震后ML≥4.0 余震的发震时间和时间间隔大体呈线性分布，表明门源MS6.9 地震后ML≥4.0 地震时间间隔正常，为主—余型。根据祁连地震带5 级以上地震最大余震发生时间统计和震级差特征，分析认为门源MS6.9 地震最大震级余震已经发生，即为1 月12 日18 时20 分的MS5.2 地震。

门源MS6.9 地震的余震序列震中形成了西段NWW、东段SE 的带状分布，主震及余震序列的震源机制解以走滑型为主，也存在部分其他类型，反映了局部复杂的构造形态，后续还需进一步开展震源破裂特征和该区域地震危险性研究。