1 研究背景

地震发生后，人们往往关心的是后续地震趋势，有无更大余震发生。要对地震做出相应的趋势预测，主要是基于对发震地区地质构造特征、历史地震活动性和地震序列特征进行统计分析。蒋海昆等（2006）对中国大陆地区M≥5.0地震序列特征已进行系统的研究，为震后地震趋势的分析积累了大量地震序列震例资料。赵小艳等（2021）对漾濞M_S 6.4地震序列特征及强余震判定进行研究，通过分析历史地震序列类型、h值等判定地震序列为前震—主震—余震型，并利用云南地区6级地震最大余震发震时间统计和震级差分布特征判定最大余震已经发生，为地震序列后续发展趋势判断提供有力的理论依据。因此，大地震发生后，对其序列类型的判定和余震特征的分析具有较强的研究意义。

2022年1月2日云南省丽江市宁蒗县M_S 5.5地震发生后，出现人员受伤现象，产生了一定社会影响。及时进行余震趋势研判，对于抗震减灾具有重要意义。

2 研究内容及方法

历史地震序列类比统计分析是判定地震序列类型及余震活动趋势的重要手段之一。为分析本次地震序列特征，分别对滇西北和整个云南地区5.0—5.9级地震序列类型进行统计。基于历史地震活动类比、地震目录、序列活动特征等，根据序列类型早期判定依据，对地震序列类型及其余震活动特征进行分析。

3 研究结果 3.1 序列类型判定

将有历史记录以来，滇西北和云南地区5.0—5.9级地震序列类型进行统计，结果见图 1。由图 1可见，滇西北地区M_S 5.0—5.9地震主余型序列占比约70%，震群型（含双震）占比约28%，孤立型序列占比约2%，主余型和震群型地震序列均明显高于整个云南地区。根据序列中最大地震释放能量E_max与全序列地震释放能量E_total之比R_E和最大地震与次大地震的震级差ΔМ = M_max - M_次来进行序列的早期类型判断：2022年1月2日至1月13日宁蒗地震序列，最大地震M_S 5.5地震占整个序列能量的90.46%，90%≤R_E＜99.99%，且本次宁蒗地震序列中最大地震与次大地震（M_S 4.6地震）之间震级差ΔМ = 0.9，0.6≤ΔМ≤2.4，均符合主余型地震系列的判别标准。

统计分析2022年宁蒗M_S 5.5地震震中半径100 km范围内历史地震序列特征，结果显示：1985年至2022年9月发生8组（10次）M≥5.0地震，其中7组地震为主余型序列，仅1组为震群型，主余型与震群型（含双震）分别占比87.5%、12.5%，相较于滇西北地区，此次宁蒗地震序列周边地区主余型序列比例更高。综合分析认为，2022年1月2日宁蒗M_S 5.5地震序列属于主余型的概率较大。

3.2 余震特征分析

一般，序列最大余震震级M₁与主震震级M₀基本成正比。根据巴特定律，序列主震震级M₀与最大余震震级M₁的震级差ΔМ = M₀ - M₁为一个常数，即ΔМ = 1.2，但ΔМ有一定区域差异。统计分析宁蒗M_S 5.5地震震中半径100 km范围内历史地震序列主震与最大余震变化范围，其0.8＜ΔМ＜2.0。结合序列M—T图和日频度图（图 2）来看，主震发生后，余震序列活动频度较弱，仅主震当天发生1次3.9级余震，4月16日M_S 4.6强余震发生后余震数量显著增多，且后续仍有4级强余震发生。

使用Zmap程序计算得到主震发生后24小时内最小完整性震级M_c为1.2。基于计算的完整性震级M_c，分别计算宁蒗M_S 5.5地震后整个序列和4月16日M_S 4.6余震后序列的h、b值（图 3）。1985年以来宁蒗地震震中100 km范围内M_S≥5.0地震的b值基本维持在0.5—1.15之间；h值基本维持在0.9—1.2之间。序列参数类比显示，本次地震b值为0.82、h值为0.42[图 3(a)]，均低于该区域平均水平。4月16日M_S 4.6地震发生后，序列参数进一步升高，h值为1.47，b值为1.04[图 3(b)]。值得注意的是，h值出现上翘现象，表明余震序列存在起伏活动。

综合分析认为，本次宁蒗M_S 5.5地震后最大余震可能已经发生，但根据序列活动特征和参数分析，仍需注意后续4级左右地震发生的可能。

4 结束语

地震序列类型判定是震后趋势预测的重要内容和基础性工作，对后续地震活动的分析具有较好的指导意义。本次宁蒗M_S 5.5地震发生在滇西北地区，1985年以来该区域内发生的地震多以主余型为主。结合历史地震序列类比及参数分析，认为本次宁蒗地震为主余型地震，且最大地震可能已经发生。但整个序列h值和b值拟合效果较差，低于该区域平均水平，且序列余震活动起伏明显，仍需注意后续4级地震发生的可能。