收稿日期：2022-10-12基金项目：国家自然科学基金（U1901602）；广东省科技计划项目（2019B020208014）；南方海洋科学与工程广东省实验室（珠海）创新团队项目（99147-42080024）联合资助。作者简介：郜怀龙（1988-），男，硕士，工程师，研究方向：测绘工程精密测量与变形监测方向。E-mail：798061889@qq.com通信作者：李翀（1993-），男，汉族，硕士，助理工程师，研究方向：地球物理学重力测量方向。E-mail：452692262@qq.com
Analysis on Intensity Characteristics of Cross-fault Deformation Anomaly in Eastern Guangdong
GAO Huailong，LI Chong，LIAO Guijin，SHAN Heqiu
Guangdong Earthquake Agency ， Guangzhou 510070， China
Eastern Guangdong;Cross-fault deformation;Intensity of deformation anomaly;Precursory anomaly
DOI: 10.13512/j.hndz.2023.02.05
备注
摘要
全文
图/表
参考文献
对粤东地区“榕江断裂”和“博美—潮安—三饶断裂”的跨断层形变数据进行整理。在形变趋势累计的基础上，进一步对形变速率处理方法改进，解算出跨断层形变异常强度值，并做归一化处理。将形变异常强度指标引入至流动跨断层测量资料的处理分析中，应用于跨断层形变曲线的异常辅助识别，并结合粤东地区震例开展对比研究。分析发现粤东地区地震多发在异常强度高或者由高降低的时段内，对粤东地区的地震前兆识别有一定指示意义。
The cross-fault deformation data of "Rongjiang fault" and "Bomei-Chao'an-Sanrao fault" in eastern Guangdong were sorted out. On the basis of the cumulative deformation trend， the deformation rate processing method is further improved， and the abnormal intensity value of cross-fault deformation is calculated and normalized. The abnormal intensity index of deformation is introduced into the processing and analysis of flow cross-fault measurement data，which is applied to the abnormal auxiliary identification of cross-fault deformation curve， and the comparative study is carried out with the earthquake cases in eastern Guangdong. It is found that earthquakes in eastern Guangdong occur frequently in the period of high abnormal intensity or decrease from high， which has certain indicative significance for earthquake precursor identification in eastern Guangdong.
引言
跨断层水准测量作为地震监测预报重要手段之一，通过对布设在断裂两侧的控制点进行周期性监测，可以直监测到断层的相对运动情况，粤东地区是我国东南沿海区域新构造运动最强地区之一，孕育有多条深大断裂带，其中布设有跨断层监测的断裂带包括：榕江大断裂和博美—潮安—三饶断裂。对跨断层监测数据的处理，先后有很多位学者进行了研究，将断层视为刚性块体研究断层的三维运动；构建基于混合回归理论的动态数据处理模型；利用断层位移速率的空间分布特征拟合逼近震源位置[1-3]；使用自适应控制理论分析形变数据，构建慢时变前兆系统结构[4]，以上的各种研究对地震异常前兆的提供了理论依据，对地震的监测预报提供了一定的参考意义。本文在该领域学者研究的基础上，利用跨断层形变异常强度指标，研究提取前兆异常信息，并结合粤东地区实际震例进一步分析，验证了使用跨断层形变异常强度指标辅助识别前兆异常的可行性。
-
1 研究方法与原理
通过绘制跨断层形变异常时空动态演化曲线，以辅助识别地震前兆异常信息。此方法以断层总体的平均变化速率绝对量为基值，根据历年的跨断层水准监测资料计算相对运动速率，以相对运动速率在基值中的占比反映形变异常强度，根据异常强度值绘制出异常时空动态演化曲线[5-7]。
首先，断层的相对变化速率是用公式（1）求解，i为观测量的序列号，y为断层监测点的位移量，t为监测的间隔时间，v是断层位移量与监测间隔时间的偏导数，表示断层形变速率值。
该方法用于去整体周期趋势化[8-10]，可以减小环境因素的干扰，能够更真实的反映出断层的相对运动情况。式（1）中先解算相同观测月份（序列号为k）不同观测年份（序列号为j）速率v'j,k的最或然值，式（3）将每年每月份的速率值vj , k减去相应月份的速率最或然值
，得到总体性去周期值v'j,k，其中m为总的观测年份。此方法是借鉴了大地水准面的计算方法，与相邻观测年份同观测周期的数据做差，可以一定程度上消除环境等外界因素对数据的影响。
解算出总体性去周期值v'j,k后，该值与最或然值
比率值，即为异常强度值ki，该值是一个无量纲的量。
为统一不同测区间的断层形变异常强度值ki，将该值进一步做归一化处理，得到归一化后的形变异常强度k'i，k'i的取值由0到1表示异常强度逐渐递增，由式（6）可知当k'i=0时，形变异常强度为零，表示相邻两期断层形变速率相同，断层的活动状态趋于稳定；k'i=1时，形变异常强度为最大，表示相邻两期断层形变速率强烈，断层的活动状态趋于增强，预示着前兆异常的生成，有孕育地震形成的可能。
但需要注意的是，形变速率的变化不仅与震中的距离、深度、震级有关，也会受到震前出现的闭锁、微破裂、预滑移和测区降水、日照、气候等的影响，因此需要通过收集断层区域的震例，统计震中位置和震级等信息，绘制出监测区域的震级时序图，与形变异常强度对比分析其相关性，进一步探究形变异常强度诱发地震的规律。
按照以上的研究思路，收集整理粤东地区的跨断层水准测量资料，结合断层运动速率值和形变异常强度值，绘制粤东区域的断层形变异常强度时空演变曲线，结合监测区域中的震级时序曲线图，研究其相关性。粤东地区涉及到的断层监测相关信息见表1。
-
2 跨断层监测场地情况
本次选用了粤东地区的榕江大断裂和博美—潮安—三饶断裂两处场地为研究对象，其中博美—潮安—三饶断裂贯穿饶平、潮安、普宁陆丰等地，在广东境内陆地部分长达210 km，总体走向北东40°～60°，倾角70°～80°，断裂性质为压扭性，穿行于上三叠—下侏罗统砂页岩及燕山三、
四期花岗岩之中，形成于燕山运动第三幕，断面沿走向及倾向均成舒缓波状，表现为多期的岩浆侵入，该断裂在历史发生过破坏性地震，此断裂表现为活动断层；榕江大断裂北起丰顺县北斗，沿榕江走向北西320°，倾角70°～80°，长约100 km，宽约200 m，由多条断裂组成线性构造清晰，线性构造清晰属正断层，该断裂上震例多发，历史上有发有过中大型震例，如1995年11月丰顺4.7级地震、1895年揭阳6级地震，表示该断裂活动频繁且有较大活动性。
表1 粤东地区跨断层监测信息表Table 1 Cross-fault monitoring information table in eastern Guangdong
利用布设在该两处断裂上的流动水准测线开展实验分析，其中广太场地位于粤东地区的普宁市广太镇，跨博美—潮安—三饶构造带广太段（图1~2），该场地共布设4个测段，北西向测线跨断层，北东向测线与断层平行，测线端点均埋设主、副点土层基本水准标石；玉湖场地位于揭阳市揭东县玉湖镇，跨越榕江断裂玉湖段（图3~4），该场地共布设5个测段，北东向测线垂直于断层，北西向测线与断层斜交20°左右，测线端点均埋设有主、副点土层基本水准标石（表2）。两出测量场地始建于1976年4月，于1976年9月开始第一期水准测量，每年复测6期，每年的单月份为复测时间， 1979年5月观测完第六期后便停测，直至1995年7月恢复观测，至2021年12月为止，已观测165期，为地震前兆异常的分析和地震监测预报研究提供了丰富的宝贵数据。
图1 揭阳市广太跨断层水准测量场地Fig.1 Cross-fault leveling site in Guangtai Town，Jieyang
图2 广太跨断层水准测线示意图Fig.2 Schematic diagram of cross-fault leveling line in Guangtai
图3 揭阳市玉湖断层水准测量场地Fig.3 Fault leveling site in Yuhu Town，Jieyang
-
3 形变异常强度与震例分析
3.1 震例选取在此次相关震例选择时，主要选用距离跨断层监测场地周边200 km范围内的地震（图5），以确保断层形变异常强度为前兆异常的可靠性[11]。由于地震诱发的复杂性及远距离强震的影响，对地震前兆异常的识别带来很大困难，文中在分析跨断层形变曲线和形变累计曲线的基础上，研究形变异常强度和异常强度累计率辅助识别地震异常前兆信息[12-13]，并进一步探讨其映震规律和相关性。
图4 玉湖跨断层水准测线示意图Fig.4 Schematic diagram of cross-fault leveling line in Yuhu
图5 广太、玉湖测区地震时序图Fig.5 Time sequence diagram of earthquakes in Guangtai and Yuhu survey area
3.2 博美—潮安—三饶断裂异常强度分析博美—潮安—三饶断裂共布设水准监测点5个，可用于分析的跨断层测段共有GTBF-GTN、GTBF-GTNF、 GTBF-GTX、 GTB-GTN、 GTB-GTNF、GTB-GTX六个测段，通过对近12年监测数据的对比分析和参照历年的地震趋势会商资料，选取稳定性能良好的测段GTB-GTN形变异常强度为研究对象。整理汇总该场地自2009年以来的监测数据，该测段的相邻期变量，与总体累计变量见表2。
表2 广太场地GTB-GTN测段数据表Table 2 Data table of GTB-GTN survey section in Guangtai site
图6中GTB-GTN测段形变曲线为同比临期形变量差值，数据2009年至2018年围绕均值波动， 2018年至2021年围绕均值缓慢向上波动，但在2010年、 2013年、 2021年出现异常波动，其中2010年3月达到异常波动低值-1.72 mm，2013年5月达到异常波动低值-1.52 mm，2021年5月达到异常波动高值1.68 mm；该测段的形变异常强度曲线在2010年3月达到异常高值1.00，2011年至2013年初断层活动较弱， 2013年5月达到异常高值0.88，2013年末至2018年初断层活动又转换为相对较弱， 2018年3月和11月达到异常高值0.76， 2021年5月达到异常高值0.98，表现为2018年后断层运动逐步活跃，形变异常强度所反映的断层活动特征与形变异常波动数据相吻合，说明GTB-GTN测段形变异常强度可用于辅助识别形变异常波动特征（图6）。
2013年5月GTB-GTN测段形变曲线震荡低值达到-1.52 mm，形变异常强度达到局部高值0.88，异常强度累计率达到局部高值0.3后逐步恢复稳定，反映广太场地在2013年5月份区域应力处于局部活跃高峰值，同年12月份在广东丰顺发生2.8级地震，发生在异常高值的后7个月。
2015年9月广东陆丰发生3.8级地震，广太场地GTB-GTN测段形变曲线和形变异常强度曲线未出现明显的异常变化，但按照形变异常强度计算方法，将计算基质替换为形变异常累积值，解算出形变异常强度累计率自2014年中旬出现上升趋势，2015年9月份达到最高值0.41后转折下降趋势，表明在2015年前后断层活跃度高，且该区域同时段发生过多次小规模震群，与揭阳地区地震时序图相吻合。
2020年1月广东丰顺发生3.7级地震，广太场地GTB-GTN测段形变曲线和形变异常强度曲线呈继承性规律波动，未出现明显波动异常，但该测段形变强度累计率在2019年5月打破原有波动规律，呈震荡上升波动趋势，此次前兆异常出现在地震前7个月。
3.3 榕江大断裂异常强度分析榕江大断裂共布设水准监测点6个，可用于分析的跨断层测段共有 YHBB-YHNB、 YHBB-YHNc、 YHBc-YHNB、 YHBc-YHNc、 YHDF-YHNB、YHDF-YHNc、YHD-YHNB、YHD-YHNc八个测段，通过对近12年监测数据的对比分析和参照历年的地震趋势会商资料，选取数据完整和稳定的测段YHNB-YHBc、YHNB-YHD形变异常强度为研究对象。整理汇总该场地自2009年以来的监测数据，该测段的相邻期变量，与总体累计变量见表3、表4。
YHNB-YHBc测段形变数据整体呈现规律性的波动变化，在-3~4 mm周期性震荡，并于每年的第1期达到震荡低值，第4期达到震荡高值，表现为该断层夏季垂直型变量增大，冬季垂直形变量减小的特征；形变异常累计率在2009年至2019年在小范围内震荡，其中2010年、2019年达到局部的异常高值分别为0.61、0.56，2020年开始加速震荡上升达到0.98的异常高值，表明该测段断层在2010年和2019年附近活动频繁，且2020年开始活动频繁加剧；该测段的均值累计线分别在2011年、2015年、2018年出现低值转折，2013年和2016年出现高值转折，总体在-0.8 mm附近波动，但2020年到达波动低值时，并未出现转折向上震荡，并进一步维持向下震荡趋势，打破了原有的震荡趋势和震荡低值（图7）。
图6 广太场地GTB-GTN测段形变数据图Fig.6 Deformation data map of GTB-GTN survey section in Guangtai site
表3 玉湖场地YHNB-YHBc测段数据表Table 3 Data table of YHNB-YHBc survey section in Yuhu site
表4 玉湖场地YHNB-YHD测段数据表Table 4 Data table of YHNB-YHD survey section in Yuhu site
YHNB-YHBc测段形变曲线在2013年附近为显示无特殊变化，但均值累计曲线走势在2013年出现转折点，异常强度累计率在2013年9月出现局部高值达到0.55，局部异常高值出现在丰顺地震的前3个月（图7）。
2015年9月广东陆丰发生3.8级地震，玉湖场地YHNB-YHBc测段形变曲线未发生明显异常变化，但形变异常强度累计率出现局部高值，在同年份该测段均值累计出现趋势转折变化。该震例反映出在测段形变异常变化无明显规律时，可通过形变异常强度和形变异常强度累计率来辅助识别和提取前兆异常（图7）。
玉湖场地YHNB-YHBc测段形变曲线呈现出夏高东低的波动规律，在2020年丰顺地震前后波形未显示出明显的异常变化，但形变异常强度累计率在2018年下旬和2019年1月出现了两次震荡高值，表明该断层运动开始增强，并于12个月后出现3.7级地震，同时该测段均值累计值在2019年9月打破原波动趋势（图7）。
图7 玉湖场地YHNB-YHBc测段形变数据图Fig.7 Deformation data map of YHNB-YHBc survey section in Yuhu site
YHNB-YHD测段形变数据围绕0.00 mm值附近波动，但没有出现明显的季节性波动规律，表明该测段受季节性因素影响较小，但测段在2011年1月出现局部的震荡低值达到-0.50 mm，2013年5月、2014年5月出现震荡高值分别达到0.71 mm、0.91 mm，在2020年达到震荡低值达到-2.74 mm；该测段形变异常强度累计率在2009年至2012年呈震荡下降趋势，2013年2019年在0.45 mm附近波动，2019年至2020年由震荡低值突发转变为局部震荡上升，整体测段上分别在2011年1月、2013年、2014年、2015年、2020年出现局的震荡高值，表明在这几个时间段的断层活动速率增强，属断层运动的活跃期，形变异常强度累计率反映的运动特征部分与测段的形变数据吻合；测段的均值累计在2012年，2020年达到波动高值分别为-2.14 mm、-2.07 mm，并且在2011年打破了原曲线走向趋势（图8）。
2020年1月广东丰顺发生3.7级地震，同年1月YHNB-YHD测段均值累计出现高值达-2.07 mm，然后转折向下波动，前兆异常出现在地震异常的前1个月时间，玉湖场地YHNB-YHBc测段和YHNB-YHD测段异常前兆呈现出一定的相关性，也再次表明该断层监测场地对地震前兆异常反应的可靠性。
图8 玉湖场地YHNB-YHD测段形变数据图Fig.8 Deformation data map of YHNB-YHD survey section in Yuhu site
-
4 结论
本文利用断层形变异常强度和断层形变异常强度累计率方法对粤东地区“榕江断裂”和“博美—潮安—三饶断裂”的跨断层形变数据进行研究，并结合断裂附近区域的地震序列分析得出如下结论：
（1）“榕江断裂”和“博美—潮安—三饶断裂”之间的跨断层形变数据、断层形变异常强度、断层形变异常强度累积率和形变均值累积量呈现出一定的相关性，均会在震前7～18个月期间出现不同程度的前兆异常，通常前兆异常表现为形变异常强度波动幅度的增大，形变异常强度和形变异常强度累计率的震荡高值和局部高值，形变均值累积继承性趋势的变化和震荡高低值的趋势转折变化。
（2）榕江大断裂的YHNB-YHBc测段形变数据受季节性影响较大，整体呈现夏季高冬低的波动变化，因此该测段前兆异常的判别和提取有一定困难，通过结合对形变异常强度和形变异常强度累计率判读，并利用YHNB-YHD测段数据加以佐证，可以一定程度上排除季节变化对形变数据的影响，进一步提高前兆异常的判别和提取的准确性。
（3）粤东地区断层形变异常强度和断层形变异常强度累计率能够很好地反应监测断层的能量聚集程度和断层活跃强度，在断层形变无明显异常时，用于辅助识别提取前兆异常具有一定的可靠性。
- [1]葛计划，宁斌，孙军，等.断层形变异常强度在跨断层水准测量中的应用[J].防灾科技学院学报，2011，13（4）：19-24.
- [2]李杰，唐廷梅，荆强，等.跨断层形变测量异常特征分析[J].地震，2010，30（2）：100-110.
- [3]张超，王世进，王宁，等.垂直形变趋势累计率在川滇地区的应用[J].震灾防御技术，2017，12（2）：329-337.
- [4]贾鹏，张希，路珍，等.利用最小二乘配置识别川滇地区跨断层形变异常的研究[J].大地测量与地球动力学，2016，36 （11）：972-976.
- [5]吴少华，程朋根，付博.改进的最优非负变权组合模型及其应用[J].大地测量与地球动力学，2016，36（11）：968-971
- [6]李福利.宁夏跨断层水准测量资料震前异常指标分析[J].防灾减灾学报，2013，29（3）：35-38.
- [7]陆明勇，于海生，李志雄.华北地区跨断层流动形变监测场地优化初探[J].中国地震，2018，34（4）：837-844.
- [8]王喜龙，杨梦尧，孔祥瑞，等.辽宁朝阳—北票断裂跨断层水准与断层土壤气地球化学特征对比分析[J].大地测量与地球动力学，2021，41（6）：600-604.
- [9]王永江，刘永强，李子涛.辽宁地区跨断层场地水准测量与地震监测预测[J].防灾减灾学报，2011，27（3）：50-55.
- [10]周海涛，薄万举.利用跨断层形变资料求解断层运动及断裂带应变方法研究[J].大地测量与地球动力学，2013， 33（6）：31-35.
- [11]郜怀龙，王发胜，禹小军，等.韶关台VS型垂直摆映震效能分析[J].华南地震，2020，40（3）：84-91.
- [12]曹志磊，周琼，葛计划.郯庐断裂带安徽段及“霍山窗”断层活动特征与地震关联性研究[J].大地测量与地球动力学，2019，39（7）：681-685.
- [13]洪敏，张勇，王伶俐，等.跨断层扰动协调比异常提取方法与效能评估及其在云南地区的应用[J].地震研究， 2023，46（1）：99-106.
