基金项目:国家自然科学基金项目(41604055,41676057);广东省科技计划项目(2019B020208014 ,2018B020207011);南方海洋科学与工程广东省实验室环南海地质过程与灾害项目(99147-42080024)联合资助
作者简介:邓志辉(1962- ),男,博士,研究员,主要研究方向为地震异常信息挖掘与预测方法研究。E-mail:deng6789@163.com.
(1.广东省地震局,广州 510070;2.南方海洋科学与工程广东省实验室,广东 珠海 519000)
(1.Guangdong Earthquake Agency, Guangzhou 510070,China;2.Southern Marine Science and Engineering Guangdong Laboratory, Zhuhai 519000, China)
Huazhou-Beiliu earthquake;Bobai-Maoming fault; Seismic activity;Flowing gravity;Underground fluid; Borehole strain
DOI: 10.13512/j.hndz.2020.02.01
2019年10月12日,在广东省化州市与广西壮族自治区北流市交界地区发生5.2级地震,这是华南地区进入20世纪以来第Ⅴ活跃幕之后发生的第3次5级以上地震,也是粤桂交界地区第Ⅴ活跃幕以来第2次5级以上地震。此地震发生北西向巴马—博白—茂名断裂带的陆川—高州段上,属6级潜在震源区,基本烈度为Ⅵ度区。震中区在历史上属于少震弱震区,但从2014年以来,地震活动出现间歇性增强,多项地球物理观测数据也时有异常表现。广东省地震局在2018年度地震会商时把此区域确定为2019年度地震危险重点注意地区。10月12日5.2级地震的发生验证了年度地震危险重点注意地区的判定,为以后的地震预测工作提供了可参考的分析依据。
On October 12th,2019, an MS5.2 earthquake occurred at the junction of Huazhou City,Guangdong Province and Beiliu City,Guangxi Zhuang Autonomous Region,which is the third earthquake with MS≥5 after the South China region entered the second active period and the fifth active scene,and is also the second earthquake with MS≥5 since the fifth active scene in the border area between Guangdong and Guangxi.This earthquake occurred on the Luchuan-Gaozhou section of the Bama-Bobai-Maoming fault zone in the northwest direction, which is a potential source of magnitude 6 with a basic intensity of Ⅵ.The epicentral area is a weak earthquake areahistorically,but since 2014,seismic activity has increased intermittently,and many geophysical observations also have anomalous performance. The Guangdong earthquake agency identified the area as a keyearthquake hazard area in 2019in the annual earthquake conference(2018).The occurrence of the MS5.2 earthquake has verified the judgment of the key earthquake hazard areaof the year and provided a reference for the future earthquake prediction.
经中国地震台网测定,2019年10月12日22时55分在广东省化州市与广西壮族自治区北流市交界地区(北纬22.18°,东经110.51°)发生5.2级地震,震源深度10.0km ,震中烈度Ⅵ度,等震线椭圆长轴北西走向,共造成广东和广西5个乡镇受灾。这是华南地区进入20世纪以来第Ⅴ活跃幕之后发生的第3次5级以上地震,也是粤桂交界地区第Ⅴ活跃幕以来第2个5级以上地震。在第五代《中国地震动参数区划图》中属于茂名6级潜在震源区。震中附近的地震区划基本烈度为Ⅵ度。
此次地震震中现场调查等震线为椭圆型分布,总体走向为北西向,震中烈度为Ⅵ度,与第Ⅴ代地震动参数区划图基本地震烈度一致,共造成26间房屋破坏,59人受灾,无人伤亡 [1]。震源机制解为第一个节面走向103.12°,倾角77.50°,滑动角-179.92°;第二个节面走向13.10°,倾角89.92°,滑动角-12.50°;P轴走向327.42°,倾伏角8.86°;T轴的走向58.80°,倾伏角8.75°;B轴的走向192.74°,倾伏角77.50°[2]。与区域应力场方向总体一致。
截止到2020年3月31日,震中50km范围内共发生1级以上余震96次,其中1.0~1.9级81次,2.0~2.9级13次,3.0~3.9级2次,最大余震为2019年10月12日的3.2级地震。
此次地震发生在广东化州与广西北流交界地区,地震距离茂名市约68 km。本区地质构造比较复杂,区域构造总体格局为一组北西向的活动断层和一组北东向的活动断层互相切割,形成网格状构造(图1),其中北西向断层主要有发震构造巴马—博白—茂名断裂(F9)、西南侧那丽—合浦—沈塘断裂(F10)、东侧丰头河断裂(F8);北东向断层裂主要有中部廉江—信宜断裂带(F3)、西北侧合浦—北流断裂带(F2)和东南侧四会—吴川断裂带(F4)[3]。
此次地震主要发震构造为巴马—博白—茂名断裂,该断裂带东南始于广东省茂名一带,往西北经高州、广西博白、横县、昆仑关、大化、巴马,而后进入贵州省境内,总体走向310°~330°,全长达800多公里。倾向以北东为主,倾角40°~85°。此断裂带最早形成于华力西期构造旋回,印支期强烈活动,表现出左旋剪切—挤压性质。断裂切割寒武系至第三系,在新生代以来和第四纪时期具强烈的活动[4]。1958年灵山级地震就发生在该断裂带上,1936年灵山级地震也与此断裂带的关系密切。该断裂带由多条断裂组成,断续分布,自西北而东南包括:高田—横县断裂、镇龙断裂、单塘—石南断裂、平山—永安断裂、六万山断裂、乌石—石鼓断裂及陆川—高州断裂等。陆川—高州断裂是此次化州—北流5.2级的发震构造。
陆川—高州断裂西北起于陆川往南东至宝圩以南,经高州而达茂名,地表断续延伸百多公里,切割震旦系、古生界、加里东混合花岗岩及白垩系,在茂名地区控制第三纪断陷盆地的发育。断裂出露数米宽的糜棱岩及破碎角砾岩带。在高州南侧可见宽达200 m的硅化岩带。新构造时期以来,断裂在茂名地区使下第三系逆冲于上第三系之上。在陆川一带,卫片分析断裂具有左旋走滑运动。该断裂曾发生过两次5级地震和一些有感地震[5]。
有地震历史记录以来,震中周边100 km范围内共发生5级以上地震7次[6](图1,表1),其中5.0~5.9级地震4次,6.0~6.9级地震3次。最大地震为1605年7月19日广东廉江附近的6.5级地震,距离约为67 km。距离此次地震最近的5级以上地震为1778年10月29日广西玉林东南级地震,约为36 km。1970年以来震中周边未记录到4级以上地震。除1605年连续发生多次6级以上地震外,总体上表现为间发型分布(图2、图3),间隔时间81~112年,显示出时间间隔增加的趋势,最近两次时间间隔为129年(1890—2019),超过历史最大值[7-8]。
表1 震中附近历史地震记录(MS≥5.0)
Table 1 Historical earthquake records near the epicenter (MS≥5.0)M
图2 历史地震M-t图(MS≥5.0)
Fig.2 M-t map of historical earthquakes (MS≥5.0)
从1850年至2018年,区域中强地震应变能的累积释放曲线显示该区积累了发生M≥5.0级地震所需的应变能(图4)。1970—2018年ML≥2.0级地震的LgN-M曲线显示,该区可能发生ML5.2(M4.8)级地震(图4)。实际地震5.2级略大于事前的估计。
在2019年度广东省地震活动危险性研究报告中,在广东地区确定了4个2019年度重点注意地区,粤桂交界中南段是其中之 一[9]。2019年10月12日广西北流—广东化州5.2级地震就发生在这个重点注意区内。这一地区被确定为年度重点注意地区在历史上是第一次,这一地区从2014年开始就陆续出现了多项地震活动和地球物理观测数据异常,主要分布在震中附近100 km范围内。
从1970年1月1日至2019年10月11日,在震中区附近50 km范围内,共发生ML1.0级以上地震473次,其中1.0~1.9级361次,2.0~2.9级100次,3.0~3.9级11次,4.0~4.9级1次,最大为2013年12月6日4.2级地震(图5、图6、图7)。
从ML1.0级以上地震M-T图和N-T图(图6、图7)来看,该区地震活动从2007年才开始活跃,主要发生在信宜与高州交界地区的高州水库,而从2014年之后,地震主要活跃区域往西迁移,频度进一步上升。2016—2017年该区地震活动相对减弱,但2018年下半年地震活动再次增强,频度显著增加(图5、图8、图9、图 10)。
图5 震中附近50km范围内ML≥1.0级地震分布图(1970年至2019年10月11日)
Fig.5 Distribution map ofearthquakes with ML≥1.0 within 50 km of the epicenter (from 1970 to October 11th,2019)
图6 震中附近50 km 范围内1970年以来ML≥1.0级地震M-T图
Fig.6 M-T map of earthquakes with ML≥1.0 within 50 km of the epicenter since 1970
图7 震中附近50 km 范围内1970年以来ML≥1.0级地震频度分布图
Fig.7 Frequency distribution of earthquakes with ML≥1.0within 50 km of the epicenter since 1970
图8 震中附近50 km 范围内2015年以来ML≥1.0级地震M-T图
Fig.8 M-T map of earthquakes with ML≥1.0 within 50 km of the epicenter since 2015
图9 震中附近50 km 范围内2015年以来ML≥1.0级地震频度分布图
Fig.9 Frequency distribution of earthquakes with ML≥1.0 within 50 km of the epicenter since 2015
图 10 区域ML≥1.0级地震活动频度变化分布图
Fig.10 Distribution of seismic activity frequency changeof earthquakes with ML≥1.0 in the region
b值是Gutenberg - Richter(1944)地震震级-频度经验关系式
log10 N=a-bM
中的拟合系数。式中,N为震级M ≥ MC( MC为最小完整性震级)地震的累积频度,a代表某一区域某个时期内地震活动率,b代表统计时段内不同大小地震之间的比例。
b值本质上反映的是不同震级地震之间的比例关系。研究表明,b值主要取决于震源区应力状况和介质性质,与应力呈反比关系[10-11],b值大,表示区域平均应力低, b值小,则区域平均应力高,因此, b值可以作为衡量区域应力水平的一个指标[12],用于地震危险区判定。此外, b值还受到构造条件、温度、流体、受力方式等因素的控制[13]。
利用粤西及邻区1970年以来ML1.0级以上地震目录,网格大小为0.005,扫描半径为2km,进行b值空间扫描。扫描结果显示,低b值区域主要集中在雷州半岛东侧、北部湾地区和广西沿海。广东化州至广西北流交界及其附近区域在5.2级地震前近3年b值呈下降变化(图 11),显示应力持续积累。
研究表明,多个地震参数异常的同时出现对中等—中强震具一定的前兆意义。在2016年7月31日苍梧5.4级地震和2019年10月12日广西北流—广东化州5.2级地震前粤桂交界地区都出现了显著的多参数异常(图 12)。
在2016年苍梧5.4级地震前,10个地震学参数除GL值、η值和b值外,其余参数均出现不同程度的异常。2016年苍梧5.4级地震后,只有C值和D值异常在短时间内恢复,其它参数仍维持苍梧地震前的异常状态,C值和D值在2018年下半年开始又出现异常,地震活动增强的地震频度N值和地震活动度S值不仅由苍梧震前开始持续异常,而且从2019年之后还逐步升高,说明粤桂琼交界及其邻海地震活动并没有因为苍梧5.4级地震而减弱,2019年10月12日发生广西北流—广东化州5.2级地震,多参数异常可能是此次地震前后的异常反应。
作为统计物理与地震学的交叉,PI方法的出发点是:地震活动可以作为受稳恒速度连续驱动的“阈值系统”;系统中的能量耗散在长期平均的意义上接近稳定,但常有随机的波尔兹曼涨落;相关性是判断系统是否具有“合作行为”的重要判据;PI方法通过对地震活动的“涨落”的分析,在中长期时间尺度上估计未来地震的发生概率[14]。
选取广东全省ML2.0级以上地震目录,计算PI空间分布,第一时段为1996年9月30日至2010年9月30日,第二时段为2016年9月30日至2018年9月30日,以第一时段PI值作为背景,第二时段空间PI值减去背景值,发现茂名地区存在高值异常点(图 13)[17-18]。
广西北流—广东化州5.2级地震前的2018年在化州—信宜—阳江地区流动重力测值表现为显著的重力正异常,年累积变化达33~42μgal,重力等值线异常分布表现为较为明显的四象限分布[15-16],且存在密集的重力等值线,反映了重力变化存在高梯度带(图 14)。区内的阳江、程村、儒洞等测点的时序曲线自2018年以来表现为由下降转折上升(图 15)。
图 14 粤西地区流动重力场年变化异常分布图(2018-03~2018-08)
Fig.14 Anomalous distribution of the annual variation of the mobile gravity field in western Guangdong (2018.3-2018.8)
广东信宜地震台1号井流量自1986年1月开始观测,观测资料连续、平稳。该井流量观测每年5至10月为流量高值,12月至次年3月常为低值,具有较清晰的年变动态,但也受降雨影响[17]。
图 16是信宜1号井月均值流量图,2016年以来流量呈上升趋势,2018年12月至2019年3月流量在降雨量减少的情况下没有下降,形态与往年不同(图 16),出现破年变现象。该井自观测以来,震中距250 km范围内发生过M4级以上地震共4次(图 16),分别为2001年阳江4.2级、2004年阳江4.9级、2016年广西苍梧5.4级和本次2019年广西北流—广东化州5.2级地震。本研究以流量5个月为窗长,1个月为步长,滑动计算,得到结果显示:除2004年阳江4.9级地震外,其余3次地震前半年左右流量都出现了破年变现象。
通过现场调查、气象因素分析、震例分析、相关数据对比分析,未找到明显的干扰源,认为2019年信宜1号井流量破年变异常,可能是由区域构造变形引起的,反映了区域应力的变化,与区域应力变化有关,有可能是10月12日广西北流—广东化州5.2级地震的前兆异常。
广东信宜地震台钻孔应变从2014年12月开始观测,具有良好的固体潮汐变化,2018年8月16日至19日期间出现了固体潮各分量不同幅度的“上升—转折—下降—恢复”畸变异常(图 17),其中北南分量和北西分量持续时间较长,幅度较大,北南分量变化幅度达0.15%,北西分量变化幅度为0.07%。东西分量和北东分量持续时间小于1 d且幅度较小。
经现场调查,可排除观测系统、远震响应、气象、爆破等因素的影响,同一观测井的水位测项也出现同步变化。分析认为2018年信宜地震台钻孔应变异常是5.2级地震前的异常表 现[18]。
综上所述,本研究通过对广西北流—广东化州5.2级地震特征的初步分析研究,得到如下几点认识:
(1)2019年10月12日广东省化州市与广西壮族自治区北流市交界地区发生5.2级地震,这是华南地区已经进入20世纪以来第Ⅴ活跃幕又一重要事件。
(2)此地震发生北西向巴马—博白—茂名断裂带上,此断裂带最早形成于华力西期构造旋回,印支期强烈活动,表现出左旋剪切—挤压性质。断裂切割寒武系至第三系,在新生代以来和第四纪时期具强烈的活动。1958年灵山53/4级地震发生在该断裂带上,1936年灵山63/4级地震也与此断裂带的关系密切。
(3)震中区在历史上属于少震弱震区,但震中周边100 km范围内历史共记录到5级以上地震7次,最大为1605年7月19日广东廉江附近的6.5级地震,距离67 km。区域上历史地震表现为间发型分布,间隔时间81~112年,并显示出时间间隔增加的趋势,最近两次时间间隔为129年(1890—2019),超过历史最大值。
(4)从2014年以来,地震活动出现间歇性增强现象,并且多项地球物理观测数据也时有异常表现。2019年10月12日5.2级地震前先后出现了中小地震活动频度增强、b值异常、地震多参数异常、地震活动PI异常、流动重力异常、温泉流量异常和钻孔应变异常等现象。
广东省地震局在2018年度地震会商时把此地震及其附近区域确定为2019年度地震危险重点注意地区。10月12日5.2级地震的发生验证了年度会商的判定,为以后的地震预测工作提供了可参考的分析依据。
本研究过程中,也发现了一些问题,例如:①由于震中附近地区地震活动背景较弱,又未发现晚更新世以后的活动断裂,地震强度估计偏低(4~5级);②流动重力观测异常偏低(33~42μgal),没有达到通常异常判定的50μgal;③震中附近地震观测台站密度较低,综合地震观测台只有信宜1个台站,而且在地震发生前,台站附近正在进行大规模的房地产开发,对信宜地震台的观测资料产生了很大的影响;等等。这些问题将在今后的工作中不断解决和改进。