黄雅 陈佳 金明培 李孝宾 茶文剑 王军 李四海
摘要：使用亚失稳台阵和固定台站组成的密集台网观测资料，对2021年漾濞MS6.4地震前滇西北地区完备震级以上地震的频度、断层总面积和b值等地震活动性参数进行时空扫描分析，发现研究区4级以上地震失稳前均有较为明显的增频、增总面积、低b值等共性特征，漾濞MS6.4地震前以小震活动为指针的断层协同化演化进程大致可分为3个阶段：2020年5—12月为应变释放点产生阶段，小震频次缓慢增加;2021年1月至5月17日为应变释放点扩展和增多阶段，小震频次逐步增多，3级以上地震释放点明显增强，且有向主震震中迁移现象，后期成带状扩展;2021年5月18—21日为应变释放连接段，中小地震活动快速增长并最终连接成片，甚至重叠，且4级以上地震集中释放。而亚失稳过程则应是从2021年1月至5月17日为准静态阶段，2021年5月18—20日为准动态阶段。最后，讨论了出现明显的亚失稳过程的原因和断层协同演化进程划分存在的问题和不足。
关键词：断层亚失稳阶段;协同化进程;地震活动性;漾濞MS6.4地震;亚失稳密集台阵;滇西北地区
0 引言
地震是一个地壳内应力累积并最终快速释放的过程，识别地震前断层所处的应力状态对于地震活动的研究极为关键。断层应力积累过程大致可以分为稳态、亚稳态、亚失稳态和失稳态，依次对应断层失稳过程中的线性、偏离线性、亚失稳和失稳4个阶段。断层的应力状态进入亚失稳的准动态阶段后，地震失稳将不可逆转地到来。因此，研究断层的亚失稳阶段特征，对于地震的短、临预测具有重要意义（王凯英等，2018;李世念等，2021;张华英等，2021）。马瑾（2016）在室内岩石实验基础上提出了断层亚失稳理论，并证明在室内岩石实验中可以准确识别亚失稳阶段。如何把室内简单岩样、清晰应力、可测参数的实验和研究成果，推广至构造应力多源、断层体系复杂、观测台站稀疏的天然野外环境中，如何识别断层亚失稳态的响应特征，给出具有可操作性的指标或判据，是地震工作者正在探索和实验的工作（宋春燕，2017）。
本文对漾濞MS6.4地震前滇西北地区（25°～27°N，99°～101°E）地震活动性特征进行时空扫描分析，计算研究区密集台阵观测编目结果的完备震级，分析漾濞MS6.4地震前断层面总面积、b值、ML≥2地震频度和ML≥3地震空间迁移等地震活动性特征，研究前震序列时空演化和震源机制解一致性等，以期为探索漾濞MS6.4地震前亚失稳过程的判定提出了地震活动参考依据。
1.1 台站分布情况
自2010年以来，维西—乔后—巍山断裂带附近发生了一系列中强地震，如2013年洱源MS5.5和MS5.0地震，2016年云龙MS5.0地震，2017年漾濞MS5.1地震。基于实验室的研究和识别结果，经过全国范围内的筛选和排查，马瑾院士研究团队最终选定滇西北地区包含维西—乔后—巍山断裂、红河断裂在内的显著应力释放区域作为断层亚失稳理论研究的天然野外观测实验区，并在此布设了由30个短周期台站组成的亚失稳野外观测研究台阵，结合该区已有的云南区域台网宽频带地震台站，在2021年漾濞MS6.4地震震中周围形成了包含42个台站、台间距约15 km的密集台阵。此密集台阵自2018年4月组网运行以来，记录了漾濞MS6.4地震序列孕育、发生和发展全过程，尤其是地震前3 a开始编目的大量微震（ML≥-1.4）数据，为本文的研究提供了详实、丰富的第一手地震活动性资料，为亚失稳态的地震活动性响应特征研究提供了不可多得的机会。
亚失稳密集台阵的观测站配备了英国GUALP公司生产的CMG-40T短周期地震计，其动态范围大于145 dB，频带宽度为2 s～100 Hz，配置有高精度时间和观测台点定位的GPS服务系统，该系统提供的最佳授时精度达1 μs，与地震计配套的数据采集器为REFTEK-130B。
1.2 地震目录完备性分析
最小完备震级是国际上普遍采用的评估地震台网检测能力的定量标准，即地震台网能以100%的概率检测到的地震的最小震级（Rydelek，Sacks，1989）。震级低于最小完备震级的地震可能会被台网漏检，从而导致地震目录不完备，进而可能会影响基于地震目录的相关研究结果的可信度（Huang et al，2001）。本文采用两种估算方法来综合估计研究区最小完备震级（MC），一是基于G-R关系（Gutenberg，Richter，1945）的震级-频度分布（FMD）方法，二是最大曲率法（冯建刚等，2012;侯金欣，王宝善，2017）。
本文选取2018年4月至2021年5月亚失稳密集台阵的地震目录，统计区域为滇西北地区（25°～27°N，99°～101°E）的最小完备震级（MC）。刘丽芳等（2012）计算的大理—丽江—盐源地区MC可控制在ML1.0～1.5，本文综合两种方法给出研究区MC为ML0.8（图2），可见亚失稳密集台阵和固定台组网明显提升了下限震级的检测能力。
2 漾濞MS6.4地震前亚失稳过程地震活动性响应特征
在实验室内观察到的断层失稳前的“协同化”进程，即应变释放区的扩展、加速扩展和连接的现象，能否在野外观测中重现，是检验实验过程和成果的重要手段。漾濞MS6.4地震发生在马瑾院士预先选定的地震重点危险区，本文尝试分析地震前滇西北地区中小震活动呈现的断层亚失稳过程响应。
选取2018年以来亚失稳密集台阵和固定台组网记录的研究区ML≥3.0地震，沿98°～102°E做时空迁移图（图3a），可以看出2018—2019年地震分布较均匀，2020年1—4月有短暂平静。2020年5—12月，99.5°E以东地区3、4级地震平静，该地段形成应变积累区，从2021年1月开始，100°E以西3、4级地震平静，该地段也形成应变积累区域，特别是2021年2月该区域3、4级地震加速活动。从图3b可以看到，2021年3—5月ML≥3.0地震有向主震震中收缩的现象，且呈现南北向条带状扩展，断层出现加速协同化过程，是应变主释放区。实验室结果显示失稳发生在应变积累区，漾濞MS6.4主震前的地震活动情况与其一致。
2.2 漾濞MS6.4地震前研究区地震活动情况
一个地区地震活动的频度和强度可以直观反映该区地震应变释放水平。本文对研究区ML≥2.0地震频度随时间演化特征进行分析，总结研究区漾濞MS6.4地震前4次ML≥4.0中强震（表1）失稳前共性特征。
图4给出了2018年4月至2021年漾濞MS6.4主震前研究区ML≥2.0地震频度图（统计窗长60 d，步长15 d）。从图中可以看出，滇西北地区ML≥4.0地震失稳前均出现ML2.0地震的显著增频现象，ML≥4.0地震前统计频次均超过均值线，ML≥4.0地震一般在上述频度超过均值线15～30 d内发生。可见ML≥2.0地震频次超过均值线可能是该地区开始以地震应变释放为主的一个重要判据。
2020年下半年至漾濞MS6.4地震前，ML≥2.0地震的频度已经多次超过均值线，至2021年1月前均未有失稳地震发生;从2021年1月开始，该区应变释放水平明显增强，随后发生了洱源ML4.1地震以及图3所示的多次ML≥3.0地震。
通常用地震频度和地震能量这两个参量作为地震带或地震区地震活动性高低的标志。由于小震的数目远远多于大震，所以地震频度实际上是由较小地震决定的，而地震的总能量却是由为数不多的几个较大地震所决定的，而小震的贡献往往被忽略了。为了克服这一矛盾，前人特定义了1个介于频度和能量之间的活动性参数：断层总面积A（t），它既反映了小地震的频度，也包含了地震的能量，从而能更全面地反映地震活动的强弱特征（宋德熹，时文勇，2005）。
断层总面积的计算如下：
A（t）=∑k-1k=k0N（k）Lk-k0 （L=4.5）（1）
式中：k0为所统计地震的能级下限;k为每个地震的能级。
本文对研究区的A（t）值进行了时间扫描，为便于同频度比较，窗长取60 d，滑动步长取15 d，同时为避免失稳地震对结果的影响，震级区间都取ML0.8～3.5。图5是断层总面积时间扫描曲线，本文取高于2 250作为异常起始指标，可得到4次震例中有3次出现了异常峰值超过2 250，占比75%，而永胜ML4.8地震前也有应变释放显示，但因偏离研究核心区，不如洱源—漾濞一带的特征显著。可见A（t）值具有比较好的预测效能。漾濞MS6.4地震前，从2021年1月开始A（t）曲线开始快速上升并高于指标值，虽然随后发生了洱源4.1级地震，但此次地震后曲线小幅回落又急剧上升，异常值高达4 536。结合地震频度和断层总面积的异常现象分析，漾濞MS6.4地震前，从2020年下半年失稳过程中应变释放协同化的标志开始呈现，即应变释放点开始产生（频度开始超过均值、A（t）值开始抬升），2021年1月底两参数均明显增大。
2.3 漾濞MS6.4地震前b值
地震是地壳应力累积到一定程度，能量快速释放的一个力学过程（尹祥础等，1987）。识别地震前断层所处的应力状态是研究地震发生过程的关键。而b值变化可以反应地下介质的应力状态变化（Scholtz，1968）。Scholtz（1968）和Mogi（1962）研究岩石破裂实验时提出，当介质的应力积累水平高时b值呈现低值，特别是岩石破裂前b值下降趋势显著，高应力闭锁断层是强震孕育的重要条件。一些室内试验也表明在岩石最终破裂前，b值会出现一个系统性的下降（Lei，Masuda，2004）。天然地震前同样观测到b值降低的现象（李全林，1979;Smith，1981;Nanjo et al，2012）。
最常用的b值计算方法有线性最小二乘拟合和最大似然估计，考虑到后者计算简便，不易受个别较大地震影响，且计算结果较为稳定，本文采用最大似然估计法求取b值，计算公式如下（Aki，1965）：
b=1lg（M-MC）（2）
σb=bN（3）
图6给出了滇西北地区ML≥4.0地震前b值扫描情况（统计时间窗长为60 d，步长为15 d，震级取ML0.8～3.5），漾濞MS6.4地震前的4次地震事件，有3次ML≥4.0地震發生前1个月内都出现低b值的现象，永胜ML4.8地震前出现时间大于3个月（表2），ML≥4.0地震失稳前均出现低b值现象，根据地震对应情况给出0.94的指标线。从2020年5月至震前，b值出现趋势性下降，特别从2021年1月开始，连续出现低b值情况，虽然其间发生了洱源ML4.1地震，但低b值仍然持续，直到发生漾濞MS6.4地震。说明2021年1月是漾濞MS6.4地震亚失稳过程的大致起始时段（对应图7中的O点），这和图3a中研究区ML3地震呈现协同化的时段相吻合，也与ML≥2.0增频和A（t）快速增加时段相吻合。
3 漾濞MS6.4地震亚失稳过程分析
由图7可知，从差应力时间上看：LM是应力积累为主阶段，地震活动呈现随机性，称作稳态;MO转为区域开始产生释放，地震开始活跃，称作亚稳态;OAB是应力释放为主阶段，也称为亚失稳阶段，其中OA段是准静态释放阶段，表现为应变释放扩展和增加，AB段是准动态释放阶段，表现为應变释放区的联结（马瑾等，2012;马瑾，郭彦双，2014;马瑾，2016）。
由图3～6可以看出漾濞MS6.4地震前存在断层亚失稳过程的地震活动性响应特征。断层变形的协同化过程一般包括应变释放区产生、释放区扩展和增加以及释放区的连接3个阶段。2020年5—12月是漾濞MS6.4地震前应变释放区产生的时段，体现在ML≥2.0 地震的60日频度多次超过19次;2021年1月起，进入应变释放区扩展和增加阶段，持续4个月。在偏离线性阶段开始时（图7中M时刻），孕震区域出现应变释放区和应变积累区;应变释放区的平稳扩展和增加与亚失稳初期阶段的准静态失稳有关;当孕震释放区足够多、已有释放区扩展足够大后，释放区间的相互作用就增强了，这时释放区加速扩展，断层进入亚失稳后期，即准动态阶段。
3.1 准动态阶段地震活动性特征
2021年5月18—21日漾濞MS6.4地震发生前，亚失稳密集台阵共记录到漾濞MS6.4地震前震1 906次，其中0级以下521次，0.0～0.9级1 009次，1.0～1.9级287次，2.0～2.9级69次，3.0～3.9级15次，4级以上3次，5级以上1次，如图8所示。
图9给出漾濞MS6.4地震两个典型前震序列，即5月18日MS4.2地震序列和5月19日MS4.4地震序列的时空展布、快速扩展和连接情况。从图9可以看出5月18日MS4.2地震序列（蓝色）高度集中在2 km×2 km的范围内，而5月19日MS4.4地震序列（红色）则朝着主震震中方向快速扩展，且展布范围无论在长度和宽度上都远超1个MS4.4的破裂尺度，呈现明确的加快扩展异常，且与前一释放区形成连接，甚至出现了重叠现象。这是亚失稳过程准动态阶段的典型空间特征。因此，5月18—19日，应该对应亚失稳过程的AB段，地震已经进入不可逆阶段。另一个不可逆的现象，是图9中5月19日MS4.4地震序列的破裂范围达8 km×4 km，展布面积是18日MS4.2的8倍。根据Wells 和Coppersmith（1994）建立的矩震级和破裂面积之间的经验关系式：
lg（RA）=a+b*MW（4）
式中：RA表示断层面面积;对走滑型地震a取-3.42±0.18、b取0.90±0.03;推算矩震级MW可达5.5左右，已经具备了前震的响应尺度特征。
3.2 准动态阶段震源机制一致性特征
地震是地下介质受到应力作用产生破裂的自然现象，震源机制解是研究构造应力场的重要基础资料，它反映了震源断层的力学特征，可以揭示地震破裂的力学机制，反映出地震等效释放应力场。陈鶞（1978）研究一些强震震例发现，强震前震源机制解趋于一致的现象时有发生，于是提出用前震震源机制的一致性来描述地震活动性异常。强震受到构造应力场的制约，一般情况下震源机制与应力场的方向吻合（李钦祖等，1982）。
本文选取2021年5月18—20日ML≥3.0的7次地震事件，采用CAP方法（Zhu，Helmberger，1996）反演出震源机制解的结果（表3），7次前震矩心深度均在4～7 km范围内，断层节面I的走向都是北西向，且走向最大差异均在10°以内，从滑动角（-163°～166°）可以判断这7次地震都是以走滑型为主。P轴显示了主压应力轴方位。图10中给出了7次前震的主压应力优势方向，结果表明这7次地震主压应力轴方向一致性较高，偏差在10°以内，并且仰角较小，接近水平方向。大震发生前孕震区内多个中小地震震源机制参数一致性反映了孕震区内变化应力场与背景应力场的关系（泽仁志玛等，2010）。这7次地震震源机制一致性较好，显示了构造应力场控制作用增强，是判定前震序列的参考判据之一，其后是否有更大地震还需要结合其它资料综合研判，但至少已经提供了1个临震的辅助判据。
4 结论与讨论
本文使用架设在滇西北地区的亚失稳密集台阵积累的观测资料，对漾濞MS6.4地震前亚失稳过程的地震活动性响应特征进行综合分析，得出以下结论：
（1）从时间扫描结果看，从2020年5月开始，研究区小震b值出现趋势性下降，至 2021年1月小于0.94，漾濞MS6.4地震前，连续4个月出现多次b值低于0.94的现象;同期地震活动频度开始逐步增强，经历多次起伏后，于2021年1月呈现频次增加、震级增大现象;断层总面积从2021年1月开始持续上升，并于2月出现异常高值。
（2）从空间扫描结果看，漾濞MS6.4地震前一年，研究区内ML≥3.0地震呈现朝震中收缩和迁移现象，随着失稳临近，2021年1月起，地震活动越来越密集，3—5月总体成南北向条状展布（与主震震源机制解P轴方向吻合）。2021年5月18—20日，两次4级前震序列在时间上加速释放，空间上快速扩展并最终连结甚至部分重叠，且5月19日MS4.4地震序列的尺度远超其应有的破裂尺度，具有了一定的前震预警功能。
（3）5月18—20日发生的7次M≥3.0前震震源机制解一致性较好，体现为走向、倾角、滑动角均在较小范围内波动，也为断层即将进入失稳阶段提供了重要的参考指标。
综上所述，本文依据马瑾院士团队在室内岩石实验得出的亚失稳阶段特征和断层协同化理论，以中小地震活动强弱为应变释放的主要标志，初步将2021年漾濞MS6.4地震的协同化进程分为3个阶段：2020年5—12月为MO段，应变释放点产生，小震频次缓慢增加;2021年1—5月17日为OA段，应变释放点扩展和增多，小震频次逐步增多，地震活动性增强，且3级以上主释放点明显增加，且向主震震中收拢，后期成南北向带状扩展;5月18—21日为AB段，中小地震活动快速增长并最终连接，甚至重叠，且出现4级以上地震集中和快速释放。其亚失稳过程则应是从2021年1月至5月17日为准静态阶段，而5月18至21日为准动态阶段。
漾濞MS6.4地震之所以呈现出较为明显的亚失稳演化特征，可能与其为走滑型震源机制为主有关，这与室内实验条件和马瑾等（2014）研究过的老虎山—毛毛山也属走滑机制相同，其余类型的地震能否识别到这一过程还需进一步观测和研究。
需要說明的是，本文研判仅仅依据地震活动性参数，加之时间较为仓促，因而还只是一个初步尝试和探索，今后还需综合同期部署的GNSS、基岩地温和其它地球物理前兆开展分析，给出更为可信的结果或佐证。
参考文献：
陈鶞.1978.用震源机制一致性作为描述地震活动性的新参数[J].地球物理学报，21（2）：142-159.
冯建刚，蒋长胜，韩立波，等.2012.甘肃测震台网监测能力及地震目录完整性分析[J].地震学报，34（5）：646-658+727.
侯金欣，王宝善.2017.2014年鲁甸MS6.5地震前后地震活动性[J].地球物理学报，60（4）：1446-1456.
李钦祖，靳雅敏，于新昌.1982.华北地区的震源机制与地壳应力场[J].地震学报，4（1）：55-61.
李全林.1979.地震频度—震级关系的时空扫描[M].北京：地震出版社.
李世念，马瑾，汲云涛，等.2021.亚失稳准动态及同震过程变形场时空演化特征——实验与分析[J].地震地质，43（1）：1-19.
刘丽芳，李志海，蒋长胜.2012.云南地区地震目录最小完整性震级研究[J].地震研究，35（4）：491-499.
马瑾，SHERMAN S I，郭彦双.2012.地震前亚失稳应力状态的识别——以5°拐折断层变形温度场演化的实验为例[J].中国科学：地球科学，42（5）：633-645.
马瑾，郭彦双.2014.失稳前断层加速协同化的实验室证据和地震实例[J].地震地质，36（3）：547-561.，
马瑾，马胜利，刘力强.1995.地震前异常的阶段性及其空间分布特征[J].地震地质，17（4）：363-371.
马瑾.2016.从“是否存在有助于预报的地震先兆”说起[J].科学通报，61（S1）：409-414.
任俊杰，张世民，侯治华，等.2007.滇西北通甸-巍山断裂中段的晚第四纪滑动速率[J].地震地质，29（4）：756-764.
宋春燕.2017.断裂亚失稳阶段及失稳部位特征的初步研究[D].北京：中国地震局地质研究所.
宋德熹，时文勇.2005.断层总面积理论在中强地震预测中的应用[J].河南理工大学学报（自然科学版），（3）：200-204.
王凯英，郭彦双，冯向东.2018.应力时空演化揭示出的汶川地震前亚失稳过程[J].地球物理学报，61（5）：1883-1890.
尹祥础，李世愚，李红，等.1987.从断裂力学观点探讨b值的物理实质[J].地震学报，4（4）：364-374+444.
宇津德治.1981.地震学[M].北京：地震出版社.
泽仁志玛，刁桂苓，李志雄，等.2020.大震前显示的地震震源机制趋于一致的变化[J].地震，30（1）：108-114.
张华英，杨建文，高琼.2021.2021年漾濞MS6.4地震前兆异常分析[J].华南地震，41（2）：176-182.
Seismic Characteristics in the Meta-instability Stageof the 2021 Yangbi，Yunnan MS6.4 Earthquake
HUANG Ya，CHEN Jia，JIN Mingpei，LI Xiaobin，CHA Wenjian，WANG Jun，LI Sihai
（Dali Center of China Earthquake Science Experimental Site，Dali 671000，Yunnan，China）
Abstract
Keywords：fault meta-instability stage;synergism process;seismicity;the Yangbi MS6.4 earthquake;meta-instability dense array;northwest Yunnan