2021 年5 月21 日漾濞4 次MS≥5.0 地震前后b 值变化＊

曾宪伟 赵小艳 李蒙亚 李文君唐 浩 袁 媛

1） 中国银川 750001 宁夏回族自治区地震局

2） 中国昆明 650224 云南省地震局

引言

2021 年5 月21 日21 时21 分至22 时31 分，云南大理州漾濞县接连发生4 次MS≥5.0 地震，分别是MS5.6，MS6.4，MS5.0 和MS5.2 地震（以下简称漾濞地震），其中MS6.4 地震打破了2014 年10 月7 日云南景谷MS6.6 地震后云南地区长达6 年多的M6.0 以上地震平静。统计发现，此次漾濞地震震源区为少震区，且地震活动强度不大。震中西北30 km 左右的洱源地区1970 年以来偶有M5.0 地震活动，最大震级为MS5.5。

漾濞MS5.6 地震前，2021 年5 月18 日18 时至5 月21 日21 时，漾濞地区接连发生了13 次ML3.0—4.9 地震和上百次ML3.0 以下地震，其中ML4.0—4.9 地震3 次，最大地震为5 月19 日20 时5 分的ML4.5 地震。根据震群的定义（国家地震局科技监测司，1990），此次地震序列已然构成震群（以下简称漾濞震群）。同时，漾濞震群尚在持续活动时发生了漾濞4 次MS≥5.0 地震。与以往漾濞地区地震活动特征对比可知，此次漾濞震群较以往震级明显偏大、频度明显偏高，1970 年以来该地区尚无类似强度和频度的小震震群活动，表明震源区及附近区域地壳介质应力水平较以往偏高。从统计分析的角度看，根据漾濞地区的历史地震活动水平，很难判定漾濞震群后区域地震危险性如何，也很难预测漾濞震群后漾濞将发生4 次MS≥5.0 地震。本文拟通过分析漾濞4 次MS≥5.0 地震前后的b值变化，重新认识漾濞震群与漾濞地震的关系。

古登堡—里克特关系式（Gutenberg，Richter，1944）中的b值具有较明确的物理含义。多数岩石破裂实验表明b值代表了介质内部应力水平的高低，并随应力增加而下降（Scholz，1968；张智等，1987；曾正文等，1995；刘力强等，2001；Amitrano，2003；李小军等，2010），即b值与应力水平呈反比（Wyss，1973；Urbancicet al，1992；Schorlemmer，2004）。许多震例研究表明强震前b值出现下降变化（马鸿庆，1978；Wiemer，Wyss，1997；Wysset al，2004；王辉等，2012；易桂喜等，2014；邵延秀等，2015；史海霞等，2018；张帆等，2018；韩佳东等，2019；曾宪伟等，2020，2021），因此，b值作为监视破坏性地震孕育过程的一种手段，可以反映一个地区承受平均应力和接近岩石破裂强度极限的程度（李全林等，1976）。这一认识与大多数实验及观测结果相符。

本文拟利用漾濞地震震源区及其附近地区2015 年以来记录的小震资料，通过分析2021 年漾濞4 次MS≥5.0 地震前后b值的时空变化特征，研究漾濞震群的前震意义，探讨震源区的应力变化过程。

1 资料选取与最小完整性震级分析

1.1 资料选取

本文以2021 年5 月21 日漾濞地震震源区及其附近区域为研究区（图1），选取2015 年1 月1 日至2021 年6 月4 日的地震资料，开展震前与震后b值空间变化的分析研究。该时段地震的震相报告通过全国统一编目系统（中国地震台网中心，2020）获得，其中2015 年1 月至2021 年3 月震相报告为全国台网正式观测报告，2021 年4 月以来为云南台网快报观测报告。

图1 2015 年1 月1 日至2021 年6 月4 日研究区震中分布Fig. 1 Distribution of epicenters from January 1，2015 to June 4，2021 in the studied area

2015 年至2021 年漾濞地震前，研究区主要以小震活动为主（图1，2），仅于2016 年5 月18 日和2017 年3 月27 日发生过云龙MS5.0 地震和漾濞MS5.2 地震。从历史地震活动情况看，漾濞地区鲜有M5.0 以上地震发生，此次漾濞地震的强度和频次均打破了以往对该地区中强地震活动的认识。在保证研究区完备震级的条件下（最小完整性震级Mc＝ML1.5，具体分析见2.2 节），选取截止震级Mcut-off＝ML1.5 绘制地震密度分布图（图3）。图3 为2.3 节b值空间扫描计算时同步绘制的图像，研究区网格划分为0.01°×0.01°，每个节点的地震密度为搜索区域内每平方千米的地震次数N。从研究区地震分布密度看，M5.0 以上地震震中区为小震密度最大的地区（图3）。

图2 研究区ML≥0.0 地震M-t 图（a）和震级频次图（b）Fig. 2 M-t （a） and magnitude-frequency （b） diagrams of the ML≥0.0 earthquakes from January 1，2015 to June 4，2021 in the studied area

图3 研究区ML≥1.5 地震密度分布图Fig. 3 Distribution of seismic density with ML≥1.5 for the studied area

1.2 最小完整性震级分析

最小完整性震级Mc的大小与测震台站的分布密度相关，故Mc往往存在区域差异。研究区Mc分布结果显示（图4），研究区中北部云龙以东地区Mc为ML1.0 左右，中南部漾濞地震附近区域Mc为ML1.8 左右，其它地区Mc基本介于ML1.2—1.4，故研究区Mc基本介于ML1.0—1.8。

图4 研究区最小完整性震级Mc 的空间分布Fig. 4 Spatial distribution of minimum complete magnitude Mc

1.3 扫描参数设置

进行b值空间扫描时，参数设置一般采用固定搜索半径或固定计算样本量的方式。前者适合于地震密度较高且分布较均匀的区域，若地震稀疏则样本量不足将导致无法计算b值；后者根据地震分布密度调整每个节点的搜索半径，适合地震分布不均匀的地区。本文采用固定计算样本量的设置方式，将研究区划分为0.01°×0.01°的网格，每个节点的地震次数固定为100，并满足大于最小完整性震级Mc的地震次数至少为20，分析b值空间扫描的所有节点搜索半径占比（图5），图5 显示80%以上的节点搜索半径在5—17 km。因此，大部分地区未出现因搜索半径过大而降低b值空间分辨率的情况，说明b值扫描参数的设置是合适的。

图5 b 值空间扫描时不同扫描半径所占节点数的比例Fig. 5 Proportion of nodes with different search during b value space scanning

2 结果分析

2.1 重新定位

为了得到更可靠的b值空间分析结果，本文对研究区的地震进行重新定位，地壳速度结构取自王帅军等（2015）的深地震探测研究结果。首先利用Hyp2000 定位方法（Klein，2007）对研究区内的地震进行绝对定位，再利用双差定位法（Waldhauser，Ellsworth，2000）进行相对定位。由于双差定位后往往会造成部分地震丢失，为了满足b值计算的需要，本文采用将双差定位地震与Hyp2000 定位地震进行合并的方法（曾宪伟等，2021）以保证地震目录的完整性。重新定位的地震共计6 097 次，其中ML0.0—0.9 地震1 193 次，ML1.0—1.9 地震3 746 次，ML2.0—2.9 地震958 次，ML3.0—3.9 地震164 次，ML4.0—4.9 地震29 次，ML5.0—5.9 地震6 次，ML6.0 地震1 次。重定位误差结果显示，双差定位水平向和垂直向误差均小于0.1 km （相对震群矩心的相对误差）的地震占98%以上，Hyp2000 定位水平向误差均值为1.8 km，垂直向误差均值为2.3 km。两种方法重新定位后的地震震源更加精确，基础数据更加可靠，b值计算的可靠性得到保证。

2.2 背景b 值计算

研究表明，低b值阈值大小与震源深度（Mori，Abercrombie，1997）以及震源机制类型（Schorlemmeret al，2005）有关。因此，不同构造区往往低b值阈值的大小不同（王辉等，2012；易桂喜等，2014；邵延秀等，2015；张帆等，2018；韩佳东等，2019）。为了确定研究区低b值的阈值大小，首先需要计算研究区背景b值的大小。本文选取2015 年以来收集到的研究区的全部地震资料，利用极大似然法（Utsu，1966；Woessner，Wiemer，2005）计算区域平均b值及其标准差，以此作为研究区的背景b值。

首先需要分析b值计算时截止震级Mcut-off对b值的影响。图6 显示，当Mcut-off≤ML1.4时，b值标准差随Mcut-off的增大而减小，b值计算结果浮动范围较大；当ML1.5≤Mcut-off≤ML1.9 时，b值标准差最小，且b值计算结果稳定；当Mcut-off≥ML2.0 时，b值标准差随Mcut-off的增大而增大，b值计算结果存在起伏。因此，计算b值时取ML1.5≤Mcut-off≤ML1.9 是合适的。震级-频度曲线（图7）的拟合结果显示，最小完整性震级Mc＝ML1.5，与b值计算结果稳定时的截止震级吻合，此时研究区b值及其标准差分别为0.65 和0.04，b值大小可表示为0.65±0.04。本文将0.65 作为低b值异常的阈值。

图6 b 值及标准差随截止震级的变化Fig. 6 b-value and its standard deviation changes with cut-off magnitude

图7 研究区地震震级-频度关系曲线（2015 年1 月1 日—2021 年6 月4 日）Fig. 7 Magnitude-frequency curve of earthquakes from January 1，2015 to June 4，2021 in the studied area

2.3 b 值空间变化

2021 年5 月18 日至5 月21 日21 时，漾濞地区发生震群活动，其中ML3.0—3.9 地震10 次，ML4.0—4.9 地震3 次，ML3.0 以下地震上百次。本文以此次震群活动时间为依据，对比分析2015 年1 月1 日至2021 年5 月17 日（A时段）、2015 年1 月1 日至2021 年5 月21 日21 时（B时段）和2015 年1 月1 日至2021 年6 月4 日（C时段）三个时段的b值空间变化，即分别对应漾濞震群前、漾濞震群后漾濞地震前和全时段的b值空间分布。

进行b值空间扫描时，按照1.3 节的参数设置，利用极大似然法计算每个节点的b值，分别绘制以上三个时段的b值平面分布图（图8），以及B时段相对A时段及C时段相对B时段的b值变化分布图（图9），即漾濞震群前后和漾濞地震前后的b值变化。结果显示：① 低b值区基本处于地震分布较密集的地区（图3、图8 和图9），这些区域搜索半径偏小，b值空间分辨率较高；②A时段漾濞震群前，漾濞地震震源区b值偏高（图8a），漾濞地震西北出现两处低b值异常区，分别发生过2016 年5 月18 日云龙MS5.0 地震和2017 年3 月27 日漾濞MS5.1 地震，文后将对这两次地震与低b值异常的关系作具体分析；③B时段漾濞震群后漾濞地震前，漾濞地震震源区出现明显的低b值异常（图8b），b值在3 天内快速下降（图9a），可能预示着局部地壳介质强度接近临界状态，区域孕震过程出现临震信号，这一认识与岩石破裂试验结果相符（Scholz，1968；张智等，1987；曾正文等，1995；刘力强等，2001；Amitrano，2003；李小军等，2010）；④C时段漾濞地震后，震源区低b值异常减弱（图8c），相比B时段震源区b值明显回升（图9b），这与区域应力释放、地震危险性降低相吻合。同时震源区小范围低b值异常依旧存在（图8c），显示余震可能持续，但震级不大。另外，还需关注震源区周围出现零星小范围b值下降现象（图9b），应与局部应力调整有关，b值下降区可能出现一些小震活动。

图8 研究区不同时段b 值空间分布Fig. 8 Spatial distribution of b-value in different periods in the studied area

图9 研究区b 值空间变化图（数字编号代表的地震与图8 相同，下同）Fig. 9 Spatial distribution of b-value variation （The numbers present the same earthquakes as in Fig. 8）

漾濞地震前后b值变化特征再次表明，将b值平面分布图（低b值异常区预示高应力积累区）和变化图（b值快速下降区预示强震危险）结合起来，可作为一种区域地震危险性判定的可行方法（曾宪伟等，2020）。

3 关于b 值剖面分布的讨论

选取B时段（2015-01-01—2021-05-21 21:00）的地震资料，沿图1 中地震展布方向绘制NW 向震源深度剖面（图10a），并将深度剖面划分为1 km×1 km 的网格，b值扫描参数设置同2.3 节，计算得到b值剖面分布图（图10b）。结果显示：①MS5.0 以上地震基本发生在低b值异常区的边缘或内部，尤其是漾濞4 次MS≥5.0 地震（图10b 中③④⑤⑥）均沿低b值异常区的边缘分布；② 低b值异常区深度基本在15 km 以浅，且异常区分布对中强地震的震源位置具有指示意义；③ 漾濞震群发生后，b值剖面西部出现20 km×20 km 的低b值异常区。漾濞MS6.4 地震震源机制为右旋走滑机制（中国地震局地球物理研究所，2021），可采用华北地区走滑型地震的震级与震源破裂长度的关系式MS＝1.86lgL＋3.821 以及震级与震源破裂面积的关系式MS＝0.954lgA＋4.134 估算震级大小（龙锋等，2006），将长度L＝20 km 和面积A＝20 km×20 km 分别代入以上两式，估算震级分别为MS6.2 和MS6.6，与实际发生震级MS6.4 基本吻合；④b值剖面中部存在一个长约35 km 的低b值异常区，该区域曾发生2017 年3 月27 日漾濞MS5.2 地震，但低b值异常尚不能完全与漾濞MS5.2 地震相对应。原因有两个方面：一是MS5.2 地震前震中附近出现低b值异常区（图11b），而震后直至2021 年6 月低b值异常仍在持续，且异常区较大，显示局部应力水平一直偏高（图7，10）；二是将低b值异常区长度L＝35 km 和面积A＝35 km×15 km （图10）分别代入以上震级与震源破裂长度和震级与震源破裂面积的关系式，估算震级均为MS6.7。因此，漾濞西北40 km 处未来存在发生强震的风险，且震源深度应在15 km 以浅；⑤b值剖面西部存在一个长约15 km 的低b值异常区，该区域曾发生2016 年5 月18 日云龙MS5.0 地震，同样低b值异常也不能完全对应云龙MS5.0 地震。分析认为：一是MS5.0 地震前震中附近未出现低b值异常（图11a），而低b值异常出现在震后（图7，10，11）；二是将低b值异常区长度L＝15 km 和面积A＝15 km×18 km（图10）分别代入上文两个关系式，估算震级分别为MS6.0 和MS6.4。因此，漾濞西北70 km 处未来存在发生6 级左右地震的风险，且震源深度应在20 km 以浅。

图10 沿断裂方向的地震深度剖面图（a）和b 值剖面分布图（b）（地震资料时段为2015-01-01—2021-05-21 21:00）Fig. 10 Depth profile （a） and b-value profile （b） along the fault（The seismic data is recorded from January 1，2015 to 21:00，May 21，2021）

图11 研究区②号地震前（a）与①号地震前（b）的b 值空间分布Fig. 11 Spatial distribution of b-value before earthquake ② （a）and earthquake ① （b） in the studied area

研究表明，低速体与高速介质的同时存在有利于应力集中而孕育地震（滕吉文，2010）。诸多震例分析也印证了地震往往发生于高速区与低速区的过渡带（Lees，Malin，1990；孙若昧，刘福田，1995；王椿镛等，2002；陈九辉等，2005；曾宪伟等，2014，2017）。贾佳（2020）利用双差层析成像法研究了洱源—漾濞地区0.25°×0.25°网格和0.15°×0.15°网格的三维P 波速度精细结构。前者沿断裂方向（NW 向）的P 波速度结构剖面结果显示，12—20 km 深度存在低速体，且低速体范围较广，在100 km 左右；后者沿断裂方向（NW 向）的P 波速度结构剖面结果显示，在14—20 km 深度存在两处显著的低速异常体，2017 年3 月27 日漾濞MS5.2 地震和2021 年漾濞地震震源位置与这两处低速体位置吻合。由此可见，研究区存在低速异常体为中强地震孕育提供了介质条件，意味着该区域20 km 以浅未来存在发生中强地震的危险。这一认识与前文基于b值剖面的地震危险性分析得到的结论是一致的。

4 结论

本文选取2015 年1 月1 日至2021 年6 月4 日漾濞地震震源区及其附近区域记录到的地震资料，利用极大似然法计算了研究区b值的背景大小为0.65±0.04，并将0.65 作为低b值异常的阈值。然后分析了漾濞震群前、漾濞震群后漾濞地震前和全时段三个不同时段的b值空间分布特征，以及漾濞震群前后和漾濞地震前后的b值空间变化特征，主要得到以下认识：

1） 漾濞震群发生后，震源区出现b值快速下降，可能预示着局部地壳介质强度接近临界状态。漾濞地震发生后，震源区b值明显回升，与区域应力释放、地震危险性降低相吻合。同时，震源区依旧存在小范围低b值异常，应与余震持续活动有关，但震级不大。震源区周围出现零星小范围b值下降现象，推测与局部应力调整有关。漾濞地震前后b值变化特征分析结果再次证明，综合分析b值平面分布图和变化图可作为一种区域地震危险性判定的有效方法。

2） 漾濞震群发生后，漾濞4 次MS≥5.0 地震均沿剖面低b值异常区的边缘分布，反映了低b值异常区的分布对中强地震的震源位置具有一定的指示意义。根据剖面低b值异常区的尺度以及震级与震源破裂长度和震级与震源破裂面积的经验关系式，推测漾濞地震震源区孕震震级为MS6.2 和MS6.6，与实际发生震级MS6.4 基本吻合，推测漾濞地震西北40 km 和70 km 处未来存在发生中强地震的风险，且震源深度应在15—20 km 以浅。

本文的b值计算和绘图程序来自zmap 程序包（Wiemer，2001），审稿专家提出了富有建设性的意见和建议，作者在此一并表示感谢。