2013-12-16巴东MS5.1地震前后波速比异常特征

陈丽娟 吴绍宇 董 蕾 魏红梅

1 重庆市地震局，重庆市红黄路339号，401147

当地下介质受到应力作用出现微破裂、扩容、塑性硬化及相变等变化时，通过地壳介质的地震P波速度变化较大，S波速变化很小，从而导致波速比出现变化，这是利用波速比研究介质物性的重要依据[1]。波速比的变化能在一定程度上反映地下介质的变化，也能间接反映构造应力的作用，所以常被用于地震研究。已有震例表明，典型的震前波速比异常呈现下降-回升-发震的变化[2-6]，某些地震前波速比呈现下降-发震的变化[7-9]，也有震例中主震前波速比出现下降-回升-下降-发震的复杂变化[10]。

2013-12-16湖北巴东(110.45°E,31.1°N)发生MS5.1地震，是三峡水库附近蓄水以来发生的最大地震。前人针对巴东地震的地震类型[11]、发震构造与成因机制[12]等开展了许多研究。早期对三峡库区的波速比研究表明，库区ML2.0～ML3.0地震发生前波速比出现震前降低-震后回升的变化[13]，但目前尚未有针对2013年巴东MS5.1地震前后波速比变化特征的研究。本文分别利用单震多台及单台多震和达法计算波速比，探讨波速比在巴东地震前后的时空演化特征。

1 计算方法

假设某个地震被n个台站记录到，通过观测报告提取P波走时tPi和S波走时tSi，计算纵横波到时差Δti=tSi-tPi，对这组tPi和Δti利用最小二乘法拟合直线求得斜率，加上1即可得到波速比。图1中白色圆为2008～2013年地震，紫色圆为2014～2021年地震，黑色三角为台站，红色线框内为单震多台法选用的地震。整理图1中2008～2021年的地震观测报告，提取出初至P波、S波的走时数据，参考文献[9]的方法计算波速比VP/VS及相关系数R、计算误差γ等。

图1 研究区域及地震分布Fig.1 Earthquakes distribution and the study area

2 计算结果

2.1 单台多震和达法波速比

利用单台多震和达法计算研究区内各个台站单台波速比，选取台站周边半径170 km内ML≥1.5地震作为研究资料，选取计算误差γ≥0.05、相关系数R≥0.95、70～90个地震的走时数据为一组资料计算得到1个波速比，在图2中以黑色空心圆表示；为提高结果的稳定性和可靠性，以10个波速比求平均、1个波速比滑动得到滑动平均值，在图2中以蓝色线段表示。结果表明，有8个台站的单台波速比出现异常变化，这8个台站分布如图3中红色三角所示。据已有研究，巴东MS5.1地震为走滑兼逆冲地震，震源机制解如图3中红色震源球所示。前人得到的震源机制及余震序列展布表明，发震断裂走向为NEE，震源区受到NEE近EW向的水平挤压作用[12]。而从图3可知，波速比异常台站分布长轴呈近EW向展布，与发震断裂的走向吻合，单台波速比异常台站分布反映了孕震区的大概范围。

图2 8个异常台站波速比随时间的变化Fig.2 The temporal variations of wave velocity ratio of the 8 abnormal stations

图3 8个异常台站分布Fig.3 Distribution of the 8 abnormal stations

由图2可见， 8个台站波速比在巴东MS5.1地震前变化形态较为相似，但略有差异。波速比随时间呈现明显的下降-回升变化，从2009～2011年主要为下降变化，2012～2013年主要为回升变化。在2017年、2018年2次秭归4级以上地震前波速比也有短时间的小幅下降或者下降-回升过程，但异常幅度和持续时间均不如2013年巴东MS5.1地震前。JIP、SHL、YCH、JIZ、YDU等台站单台波速比在2013年巴东MS5.1地震后短时间内仍持续升高，这可能是因为主震破裂后，库水大面积涌入，造成P波速度恢复，从而导致波速比增大。表1统计了上述8个台站单台波速比异常变化形态。可以看出，各个单台波速比的变化特征基本一致，均表现为下降-回升-发震，但下降及回升的持续时间有差异，这种差异与震中距的大小没有明显的关系。

表1 单台波速比异常特征

2.2 单震多台和达法波速比

利用单震多台和达法计算图1中红色线框区域内ML≥1.0地震的波速比，计算时选取震中距100 km范围内的台站记录到的纵横波走时，计算误差γ≥0.05，以600 d作为时间窗长计算出波速比平均值，并以30 d作为时间步长滑动，得到研究区内多台波速比随时间的变化特征，如图4所示，图中虚线为1倍均方差。可以看出，多台波速比的平均值约为1.69，波速比在2009～2013年整体较低，2009年约为1.68，之后开始下降，2012年降至最低水平(约为1.65)，此后2 a波速比一直回升，直至2013-12发生巴东MS5.1地震，波速比升至最高值(约为1.72)。地震后波速比略有下降但仍维持高值，2014-03秭归分别发生MS4.3、MS4.7地震。2015年波速比开始下降，至2016年降至平均值，之后小幅回升，2017-06该区秭归、巴东分别发生MS4.3、MS4.1地震。之后波速比再次出现下降-回升变化， 2018-10秭归发生MS4.5、MS4.1地震。此后波速比平稳变化，截至目前研究区内再无显著地震发生。综上，2013年巴东MS5.1地震前研究区内多台波速比呈现明显的2 a下降-2 a回升变化过程，与典型的波速比异常形态较为一致，异常持续时间约为4 a。2017年与2018年2次MS4.0以上地震发生前波速比也经历下降-回升变化过程，但下降幅度及持续时间均小于2013年巴东MS5.1地震前。

图4 多台波速比随时间的变化Fig.4 The temporal variations of wave velocity ratio obtained by the multi-station method

根据多台波速比随时间的变化，以0.3°×0.3°为空间窗口、0.1°×0.1°为空间窗长滑动，分3个时段给出图1中红色线框内波速比的空间分布，如图5所示。

图5 多台波速比空间演化特征Fig.5 The spatial evolution characteristic of wave velocity ratio obtained by the multi-station method

从图5看出，2008～2009年巴东周边波速比有部分高值，2010～2012年巴东周边波速比在大范围内呈现明显的低值，2014～2015年波速比升至高值，且高于2008～2009年。波速比低值集中分布在巴东MS5.1地震周边区域，反映震中周边地区介质在构造应力作用下产生微裂隙，P波速度下降，从而导致波速比出现大面积低值异常；震后波速比出现大面积高值是因为库水渗透、P波速度恢复所致。

3 讨 论

波速比在强震前的异常变化一般可分为2种形态，分别是下降-回升-发震和下降-发震，巴东MS5.1地震为前者。波速比在强震前的下降-回升异常现象可通过扩容-流体模式[14-15]解释，这种模式中流体起到关键作用，震中附近构造应力的积累造成微裂隙增加，引发扩容，P波速度下降导致波速比下降，而后期流体渗入裂隙中使P波速度恢复，从而导致波速比升高。

巴东MS5.1地震发生在长江三峡库区，该区岩溶发育，有利于库水持续向地下渗透，地震发生在三峡水库175 m高水位缓慢下降阶段，震中距离长江干流4.3 km。随着构造应力的积累，深部岩石扩容，波速比下降，随着蓄水容量的增加，水向下渗透速度加快，岩石孔隙液压上升，波速比开始明显回升，所以2013年巴东MS5.1地震前波速比才会出现下降-回升的变化形态，异常持续时间约为4 a，变化幅度约为4%。1962年新丰江水库MS6.1地震前波速比也呈现出类似的下降-回升变化[4]，但异常持续时间仅约1 a，变化幅度约为11%。

统计几次强震前波速比的异常变化形态，如表2所示，可以看出，波速比的异常持续时间、异常幅度与震级大小似乎并不存在相关关系。可能是因为这些地震处于不同的构造区域，区域地下介质物性不同，受到的构造应力作用方式不同，所以反映到波速比上的异常幅度及持续时间也不尽相同。

表2 波速比异常总结

研究区内在2017年、2018年分别发生数次MS4.0以上地震，这些地震前波速异常现象不再显著。新丰江水库6.1级地震后曾发生过数次MS4.0～5.0地震，这些地震前波速比异常也不显著[4]，与本文结果相似。这种现象可能是因为岩石中的裂隙水已经处于饱和或近饱和状态，即使构造应力再次积累，形成新裂隙，库水也能快速注入，所以波速比的异常就不再显著，说明水库地区波速比受库水的影响较大。

4 结 语

本文采用2008～2021年巴东周边地区地震的纵横直达波走时数据，分别采用单台多震、单震多台和达法计算波速比，分析2013年巴东MS5.1地震前波速比的时空演化特征。多台波速比的空间演化特征表明，震前波速比低值主要分布在巴东MS5.1地震震中周边地区，地震后波速比大范围回升至较高值；时间演化特征表明，巴东MS5.1地震前波速比出现2 a下降-2 a回升的变化形态。2013年巴东MS5.1地震前震中东西两侧8个台站的单台波速比出现不同程度的下降-回升变化。2种方法计算得到的波速比结果随时间的变化形态较为相似，可以相互印证。