流动重力观测是地震监测预报的重要技术方法之一。早在1964年日本新潟、美国阿拉斯加等地震前后即观测到了重力变化。我国自1966年邢台大地震后开始开展地震流动重力监测工作，并观测到了1975年海城7.3级、1976年唐山7.8级地震前后重力变化异常^[1-3]。近年来，利用流动重力观测资料，中国地震局重力学科专家对2008年汶川8.0级、2013年四川芦山7.0级、2014年云南鲁甸6.5级、2016年青海门源6.4级、2016年新疆呼图壁6.2级、2017年四川九寨沟7.0级等强震均进行了较好的年度、年中预测^[4-6]。

近年来在华南块体中强地震集中活跃的背景下，连续发生了2013年9月4日福建仙游M_S4.8、2016年7月31日广西苍梧M_S5.4、2017年7月15日广西南丹M_S4.0、2018年3月20日广东阳西M_L4.2、2019年3月5日广东雷州M_L4.1、2019年8月20日海南三亚M_S4.2、 2019年 10月 12日广西北流M_S5.2、 2019年11月25日广西靖西M_S5.2与2021年8月4日广西德保M_S4.8等中强地震。华南地区粤、闽、桂、琼流动重力观测网络在每年的上、下半年各开展1期观测^[7]，通过对华南地区流动重力测网多期重力数据精细处理与分析，捕捉到广西南丹M_S4.0、广东阳西M_L4.2、广东雷州M_L4.1、广西北流M_S5.2、广西靖西M_S5.2和广西德保M_S4.8地震前震中附近的重力异常变化，并进行了较好的中期预测^[8-11，13]。目前许多学者先后提出过质量迁移、膨胀扩容、闭锁剪力等“模式”用于解释震前重力场变化^{[2-3，14-15]}，但受地球内部的“不可入性”、中强地震的“非频发性”、地震物理过程的复杂性等困难的制约^[16]，地震预报还处于经验预报阶段。因此，从实践中不断总结经验，再应用于指导地震监测预报实践，是提高流动重力地震分析预报能力的重要方法。随着典型震例资料的积累，目前对5级以上地震孕育相关重力变化过程、范围、量级等有了更多的认识^[6，17-20]，但由于观测技术和应用需求等方面的原因，对4、5级左右地震重力变化过程、范围、量级除有一定总结外^[18，20-21]，总体关注较少，相关研究不多。因此本文通过搜集相关文献资料，并结合笔者近年在流动重力地震分析预报实践工作，以华南地区近年来发生的多次4、5级左右震例中总结提炼地震孕育与重力变化异常时空演化特征之间的关系，分析其在流动重力地震分析预报的“时、空、强”定量依据，为其它中等地震构造区重力场异常指标分析提供参考。

目前华南地区流动重力观测网络主要由广西、广东、福建、海南4个省份构成。2016年7月31日发生广西苍梧M_S5.4地震，凸显出广西—广东交界区域流动重力观测点位空间分布不均、密度不足等问题，所得到重力信息较为片面，难以捕捉地震孕育发生过程中出现的完整前兆信息，这直接制约着地震分析预测能力^[7]。因此2016年9月起粤、闽、桂、琼开始优化交界区域重力网络布局，拼接孤立区域测网，进一步强化邻省公共测段观测，优化后的重力测网基本覆盖了粤、闽、桂、琼主要活动断裂带，旨在监测区域物质迁移及构造变形，为中等地震预测研究提供区域性构造运动前兆异常信息（图1），截止到2020年9月华南地区共开展9期流动重力观测。

流动重力观测资料处理过程：①采用LGADJ软件进行平差计算，考虑到测网内绝对重力点点位稳定，外界干扰影响小，且有高精度绝对重力观测值，因此采用经典平差计算方法；②平差计算时，先对多期重力观测资料计算结果进行整体分析，初步了解各台仪器观测精度后，合理确定各台仪器的先验方差，再重新平差计算，以得到最佳解算结果；③对平差结果中极个别周围环境改变或不重合的测点进行删除；④对重力观测数据进行拟合推估，以便突出显示构造因素的重力效应。各期观测资料及整体平差计算后重力点值精度统计情况见表1，9期重力观测资料点值平均精度在8.0×10^-8µGal左右。

图1 华南地区流动重力观测网络Fig.1 The mobile gravity monitoring network in South China

表1 华南地区流动重力测量情况及精度统计表Table 1 Statistical table of mobile gravity monitoring and accuracy in South China

一般情况下，未来地震的“时、空、强”与重力变化异常的“时、空、强”呈正相关，即重力变化异常的量级、范围越大，则对应震级越大，震级越大则震前重力场演化持续的时间可能越长，重力变化异常形态与未来发震地点密切相关。近年来随着典型震例资料的积累，许多地震工作者对5级以上地震孕育相关重力变化过程、范围、量级等有了更多的认识^[6，18-20]。贾民育等^[18]采用时变距S来描述重力变化异常范围，定义时变距S为震中至重力正异常区中心的距离，是重力特征异常区的半径；胡敏章等^[20]对时变距S做两点补充说明：当重力测网稀疏时，S取震中至重力正变化最显著点距离；重力异常变化出现四象限图像时，时变距S为震中至2个正变化中心距离的平均值。重力变化异常量级G定义为地震、断裂构造运动相关重力变化差异：祝意青等^[6]、胡敏章等^[20]将重力变化异常图像出现梯度带和四象限特征时，重力变化异常量级取正负变化的最大差值；仅有局部重力变化时，重力变化异常量级取最大变化值；由于中强震前一般有较长时间的重力场演化过程，而流动重力观测时间间隔较长，无法确定哪一期观测到的重力变化异常最具代表性，因此当某地震有多期变化值时，取累计变化较为显著时期平均值作为G值。通过提炼5级以上地震震级与重力变化异常范围和量级之间的统计关系^[6，20]，得到5、6、7、8级地震震级判定的重力变化异常量级（G）参考指标分别为50、70、90、120µGal，重力变化异常范围（两倍时变距S）参考指标分别为140、220、350、660 km。本文依据重力变化时变距与变化量级的定义，统计了华南地区4、5级左右震例对应的重力变化时变距S、量级G与震级M_S（M_L）间的关系见表2。可以看出震级为4级时，重力变化异常范围约为100 km，重力变化异常量级40µGal；震级为5级时，重力变化异常范围约为150 km，重力变化异常量级50µGal。

图6 雷琼地区震前重力场动态变化图像Fig.6 Image of dynamic variation of gravity field in Leiqiong area before earthquake

表2 华南4、5级左右地震及其重力变化时变距S与变化量级G Table 2 M4 and M5 earthquakes in South China and the time-varying distance S and magnitude G of their gravity variation

本文以华南地区近年来发生的多次4、5级左右地震为例，总结提炼地震孕育与重力变化异常时空演化特征之间的关系，分析流动重力地震分析预报的“时、空、强”定量依据，得到如下结论：

（1）地震从孕育到发生是一个从应力积累到能量释放的过程。在地震孕育过程中，随着震源区应力的不断积累，使得地壳内部产生变形和发生物质迁移，致使地壳内部形变和密度发生变化，从而使得该处地表变形和重力值发生变化，主要表现为地壳形变场空间分布的有序性变化和局部形变异常区及伴生的地壳形变高梯度带出现。4、5级左右地震活动本身积累的能量较小，伴随的重力场变化范围和量级小，在重力变化成因分析时更易受局部因素影响，导致震例验证机会较小。本文通过收集华南地区震例分析表明4、5级左右地震震前出现较为显著的局部重力场变化；通过分析地震前震中附近的重力异常变化特征，归纳得出地震易发生在与构造活动有关联的重力变化四象限中心地带或正、负异常区过渡的高梯度带附近，这与以往5级以上地震震例研究基本一致^[6，19-20]。

（2）以华南地区4、5级左右地震为例研究了4、5级左右地震及其重力变化时变距与变化量级关系。分析认为震级为4级时，重力变化异常范围约为100 km，重力变化异常量级约40µGal；震级为5级时，重力变化异常范围约为150 km，重力变化异常量级约50µGal。研究表明计算的5级左右地震时变距、重力变化异常量级与胡敏章等^[20]、祝意青等^[6]的总结基本一致；4级左右地震时变距、重力变化异常量级略大于文献^[20]总结结果（震级为4级时，重力变化异常范围约为70 km，重力变化异常量级约30µGal）。由于本文是利用华南地区震例样本量对其时变距与变化量级进行初步定量分析，为更好系统总结，在今后研究工作中纳入中国大陆地区更多4、5级左右震例，进一步梳理地震及其重力变化时变距与变化量级关系，建立4、5级左右地震在流动重力地震分析预报空间、强度方面的定量依据，为其它中等地震构造区重力场异常指标分析提供参考。

（3）目前流动重力观测多为1～2期/年的复测周期，未能对震前重力场演化过程进行较详细跟踪，因此尚不能得出与时间的统计关系，在时间上流动重力以中期地震预测为主。