广西北流-广东化州交界5.2级地震前RTL异常特征

郑永通，李文杰，刘愫昀，王绍然

（1.福建省地震局龙岩地震台,福建 龙岩 364000；2.漳平市防震减灾中心,福建 漳平 364400）

0 引言

2019年10月12日22 时55 分在广西玉林市北流市（22.18°N，110.51°E）发生5.2 级地震，震源深度10 km。该地震位于巴马-博白断裂带东南端，归属东南沿海地震带西段。震中区北面由印支花岗岩基体组成的大容山、望军山，最高峰海拔分别为1 275 m 及1 147 m；南面的天堂山最高峰海拔为1 274 m，地貌上皆属中低山与丘陵区。该区域发生的地震均为中、小地震，未发生过6 级以上地震，最大的一次地震为1860年发生在玉林—陆川之间的5.5 级地震。

一般情况下，中强地震前在围绕未来主地震断层的大范围内会出现地震平静或增强现象。Sobolev等据此提出一种检测地震活动水平相对变化的区域-时间-长度算法（简称RTL 算法，下同）[1-2]，对于给定点在t时刻的RTL 值可定义为震中距、时间以及破裂尺度函数分别除以各自标准差后的乘积；蒋海昆等[3]对该方法进行了改进，并利用改进后的RTL 算法对华北5 级以上地震进行回顾性预测研究，结果表明，改进后的RTL 算法所得到的异常大多具有短期或中短期特征；张双凤等利用RTL 时间扫描方法对1980—2005年发生在冀鲁豫交界地区ML5.0 以上地震进行了回顾性检验，发现震前都达到了VRTL异常指标，说明RTL 算法对冀鲁豫交界地区的5 级以上地震活动具有一定的预报能力[4]。

本文利用蒋海昆等[3]改进后的RTL 算法对2019年10月12日发生在广西北流-广东化州交界5.2 级地震进行回顾性检验。使用该方法对震中附近进行空间扫描分析，以判定该方法在本次地震前震中及附近是否存在异常及有何特征。

1 改进后的RTL 方法简介

对任意一个时空坐标点(x,y,z,t)，假定t之前发生的所有地震对位置(x,y,z,t)处的介质变形均产生影响，并且当一个已发生地震的时空坐标距(x,y,z,t)越近，其影响的权系数越大。经蒋海昆等人改进后定义的震中距函数R(x,y,z,t)、时间函数T(x,y,z,t)、破裂尺度函数L(x,y,z,t)及VRTL（x,y,z,t）函数分别为：

式中：ri为第i个地震到(x,y,z)的距离；li和ti为第i个地震的破裂尺度和发生时间；Rtr、Ttr和Ltr分别为R、T和L的动态期望；n为满足式(2)所表达条件的地震数。

当所发生地震在空间上距(x,y,z)足够遥远(ri→∞+)、或在时间上距t时刻足够久远(t-ti→∞+)、或其破裂尺度足够小(li→0)，则其对(x,y,z)的贡献均趋于0。在实际计算中引入约束条件：

式中：Mi为第i个地震的震级；Mmin为震级下限；r0为特征距离；t0为特征时间跨度；di为第i个地震的震源深度；d0为截止深度。

第i次地震的破裂尺度li由式（3）获得：

式中：Mi为地震震级。

由式(1)可见，R、T、L及VRTL均为无量纲函数，VRTL以0 为期望，其数值表征了t时刻地震活动相对于动态背景变化率的偏离，偏离程度以其标准差σ=σR·σT·σL为单位进行表征。VRTL下降意味着(x,y,z)附近地震活动相对于背景变化率的降低，在一定程度上表征了地震活动的平静；VRTL上升则表征(x,y,z)附近地震活动水平相对于背景变化率的增高，意味着地震活动增强。

蒋海昆等[3]统计表明，当VRTL曲线明显偏离0 线、偏离幅度大于3.5σ（所有震例平均偏离程度的下限）以上，且VRTL极值大于150 时，可初步判定出现VRTL异常。震前VRTL变化型态可分为两类：I 类具有相对完整的变化形态，大多呈现从0 开始的上升－转折－下降或下降－转折－上升形态；II 类变化形态不完整，一般从0 开始急剧上升或下降，无明显转折，主震发生在VRTL极值附近。在I 类RTL 异常中，VRTL曲线从极值转折进入1σ的时间开始进行主震时间估计，在6 个月内发震的震例数约占85%，而1年内发震的比例达到96%，即多数地震是在极值转折到1σ后的一年内发生主震。

2 资料处理

2.1 广西北流-广东化州交界5.2 级地震震中RTL异常特征

在本文中采用中国地震台网测定的地震目录，根据LgN-M图，可以确定广西北流-广东化州（经度范围：109.0°～112.0°E；纬度范围：20.5°～23.5°N）的最小完整性震级Mmin取ML2.0（图1）。公式（1）、（2）式中统一取t0= 12 个月、r0= 60 km；d0不加以限制，x、y取本次地震震中位置，资料时段开始取震前6年、研究时段截止于主震前1 天（即2015年10月11日至2019年10月10日），其中4年为可预测时段（初期的资料时段Tmax= 2t0= 24 个月为不预测时段），时间t以10 天为步长进行滑动，计算后根据数据出现的具体情况再进行时段调整。

图1 广西北流-广东化州最小完整性震级的确定（2015—2019年）

根据设定的条件计算出广西北流-广东化州交界5.2 级地震前震中的Mmin值（图2），其σ值、3σ值和RTL 最大值如表1。

表1 震中RTL 的σ 值和RTL 最大值

图2 震中RTL 曲线图

在图2 中，VRTL曲线一开始在0 值附近，从2016年11月下旬开始下降波动，最低值达到-0.19，到2018年7月恢复到0 值附近，2018年9月下旬开始上升，2019年4月上、中旬上升到最大，之后开始下降，下降到0 值附近（0.23）发生地震，是一个比较完整的VRTL曲线，即曲线在0 值附近波动→上升（或下降）→最大（或最小）→下降（或上升）→发生地震。

由图2 也可以看出，震前VRTL曲线明显偏离0 值线，偏离幅度大于3.5σ（即3.5σ≤1），在2019年4月上、中旬VRTL曲线达到峰值，RTL 最大值为2 022（大于150），根据蒋海昆等的统计结果可判定此次地震震中位置震前RTL 存在异常。VRTL曲线最后降到0.23 时发生地震，0.23 小于1σ（0.257），VRTL曲线具有相对完整的变化形态，VRTL曲线为I 类。

2.2 广西北流-广东化州交界5.2 级地震震中附近RTL 异常特征

在进行地震预测预报时很难确定未来的震中位置。即便震前出现异常，也是在一个区域范围，无法固定在一个确定的点。因此，通常会采用空间扫描的方法分析某一个参数在某一个区域内是否存在异常，若存在异常，则判定这个区域内或附近有可能发生中强地震或大地震。根据以上观点，利用RTL 方法对广西北流-广东化州交界5.2 级地震的震中附近进行分析，看该区域是否存在震前异常。

在式（1）、（2）中，x、y取震中附近区域的经、纬度。即经度范围：109.0°～112.0°E；纬度范围：20.5°～23.5°N，以0.5°为步长，共49 个点进行RTL 分析。分别计算各点的RTL 值，并根据数据具体情况进行时段调整，尽量使RTL 以最大限度异常成立为目标进行时段调整，计算结束时间尽可能截止到2019年10月10日。

根据以上设定条件计算了广西北流-广东化州交界5.2 级地震震中附近各点VRTL值（图3），并计算各点曲线的σ值、3σ值和RTL 最大值，判定各点曲线的类型，根据异常是否成立的判据判定每点的RTL 异常是否成立(表2)。

图3 震中附近RTL 曲线图

表2 震中附近RTL 各点的σ 值和RTL 最大值

续表2

在表2 中，粗体字是3σ值≤1、RTL 最大值≥150，即符合RTL 异常成立条件的点。从表2 中可知，共有13 个点符合条件，这13 个点VRTL曲线类型都属于I 类。再将这13 个异常点的类型进行细化，即：异常快结束时发生主震的为IA 类，异常结束后经过一段时间再发生主震的为IB 类（表3），将符合IA 类和IB 类的点绘制到地图上（图4）。

图4 震中附近RTL 异常成立点分布图

由图4 看出，这13 个点中有8 个点分布于震中及附近南北两侧，异常类型均为IA 类。由图3 也可以看出，这8 个点的曲线都是先达到峰值，之后开始下降，在即将接近0 值线时发生地震。

在图4 中，另有5 个点分布于广东省南部海外和广东省阳江市，这5 个点的异常类型都是IB 类，这5 个点是否是广西北流-广东化州交界5.2 级地震前的异常，需进一步分析。

通过分析表3，发现这5 个点的异常结束时间均在2018年3月20日以后，因此认为这5 个点的异常是受广东省阳江市2018年3月20日M3.7 地震影响。蒋海昆等的统计表明，RTL 异常结束后1年内发生地震的比例为96%，而这5 个点异常结束到发生广西北流-广东化州交界5.2 级地震的时间达到17.8～18.5 个月，异常结束到发震已超过1年时间，因此判定这5 个点的异常已经失效。如果将截止时间定为2018年5月31日，向前推6年，起始时间定为2014年10月19日，计算这5 个点的RTL 值，则这5 个点的VRTL曲线如图5，5 个点VRTL曲线的σ值、3σ值和RTL 最大值见表4。在表4 中，除了（21.0°N，112.0°E）这个点异常成立外，其余4 个点的异常不成立。加上这5 个点距离广西北流-广东化州交界5.2 级地震震中较远，所以这5 个点的异常不能作为广西北流-广东化州交界5.2 级地震前的异常。

表3 RTL 异常成立点的异常信息

表4 广东省海外和阳江市5 个点RTL 的σ 值、RTL 最大值

图5 广东海外和阳江市5 个点截止到2018年5月31日RTL 曲线图

3 结论

从利用RTL 算法对广西北流-广东化州交界5.2 级地震震中及其附近区域的计算结果分析显示，震前震中及其南北两侧共有8 个点存在RTL 异常。这8 个点的曲线均表现为初始上升，当曲线到达峰值后转折下降，在曲线即将接近0 值线时发生广西北流-广东化州交界5.2 级地震，此次震中就位于RTL 异常区内。因此认为，RTL 方法可以作为本次地震前的一项异常指标。