1、第 卷 第期 年月大 地 测 量 与 地 球 动 力 学 ，收稿日期：项目来源：河北省地震科技星火计划（）；河北省地震局雄安新区流动重力加密观测项目（）。第一作者简介：杨雅慧，工程师，主要从事流动重力观测技术研究，：。：文章编号：（）河北地区 年重力场动态变化特征杨雅慧，刘洪良张展伟张庆扬昝书林河北红山巨厚沉积与地震灾害国家野外科学观测研究站，河北省邢台市隆尧县山口镇，河北省地震局保定地震监测中心站，河北省保定市天威中路 号，摘要：基于河北及邻区 年流动重力测网 期观测资料，采用绝对控制下平差算法，得到河北地区重力场动态变化图像，并采用多项式拟合法计算重力场平均年变率。结果显示，重力场动态变化
2、呈现明显的分区特征，西北部区域变化平稳，范围为 ；东南部区域 年存在微弱的正值累积并持续向北东方向扩散，至 年在天津邯郸一线形成 向的椭圆形显著异常区，异常区北至燕山，西至太行山，正负差异最高达 ，平均年变率高值明显。结合区域内部分水位数据及地表沉降数据推测，东南部重力场高值区主要为地表沉降的结果，地下水影响相对被覆盖。关键词：河北地区；多项式拟合；重力平均年变率；地下水流失；地表沉降中图分类号：文献标识码：河北测区流动重力测网自 年开始观测，最早由邢台地区个闭合环组成，之后陆续建立华北火车测网、文霸汽车测网、邯安测网，年火车测网停测，文霸、邯安、邢台测网组合成一个横跨河北及周边省市的大型重力
3、测网并与全国测网连接，构成大华北流动重力测网的一部分，并于 年开始正式测量。近年来，通过本测区观测结果发现，河北东南部太行山以东区域重力场正异常累积明显。为进一步了解本测区重力场长期动态变化机制及影响因素，为地震前兆信息的辨别及中长期地震预报提供理论依据，本文基于河北测区 年积累的 期相对重力复测资料，并搜集周边相邻测区的相对重力复测数据，同时以区域内白家疃、太原、泰安、蓟县等绝对点同期重力点值作为控制，得到高精度的重力平差结果，获得重力场动态变化图像；然后采用多项式拟合方法计算重力场的平均年变率，结合已有数据及研究成果，对本测区重力场变化机制进行分析。观测数据与处理河北流动重力测网在 年每年
4、月、月各进行一期观测工作，每期观测用时约为。使用台 型相对重力仪进行观测，测量精度优于 ，读数分辨率为 。整理河北及周边测区重力资料及绝对点信息，使用 处理软件进行绝对重力点控制下的平差处理，得到河北及部分周边地区的平差结果。参与计算的各期测点总数为 个，测段数为 条，平差后平均点值精度优于 ，联测精度优于 （表）。表重力数据处理精度 年份期次测点个数绝对点相对点总数测段数条精度 平均点值 联测 参与计算的河北及周边地区重力资料的整体测网分布如图所示。其他测区使用仪器多为 型相对重力仪及 型相对重力仪。重力场动态特征重力场的动态变化可以反映重力场在研究时大 地 测 量 与 地 球 动 力 学
5、年月图 年河北及周边地区重力测网分布 段内的持续演变情况，由上述重力平差结果可得到重力场累积动态变化图（图）。可以看出，年重力场累积变化较小，大部分区域呈现低负值，仅在石家庄西部、邯郸北部、天津东北部出现小范围的低正值，整体重力在 之间变化。在 年的重力累积图像中，零值线呈 向分布，其中，蓟县北部出现轻微负异常累积，最高达 ，并在之后的 年间持续累积，至 年底负异常消失，之后变化平稳，未出现异常高值。另一方面，位于河北东南部、石家庄以东区域的天津邯郸一线自 年起出现正异常并持续累积，至 年高正值异常明显，并逐渐向北东方向延伸，年底形成走向近 的椭圆形高正值异常区。该高值区西部以太行山为界，北部
6、至燕山一带，东部与山东地区重力累积高值区相连接；零值线沿太行山山前断裂带展布，高梯度带 明 显，累 积 重 力 场 正 负 差 异 最 高 达 ，区域内重力场异常变化显著。年，在异常区高梯度附近相继发生临城 、永清 、丰南 、古冶 地震，重力异常仍在持续。其他地区重力场累积变化平缓，无明显异常。图重力场累积动态变化 平均年变率变化平均年变率可以剔除重力场的短期、局部影响因素，从而反映出由于地下活动或地下物质密度变化引起的重力场长时间变化趋势。为保证数据结果的可靠性，在进行重力场年变率相关计算之前，首先对参与计算的观测数据进行筛选，删除其中测量期次少于期的测点数据，避免部分测点由于数据量少造成误
7、差。假设重力观测数据（，）包含长期变化、周期性变化以及噪声信息，则（，）可写为：第 卷第期杨雅慧等：河北地区 年重力场动态变化特征（，）（）（）（）式中，为重力场变化速率，为时间差，、分别表示年周期、半年周期项，和为振幅。由于本文使用流动观测重力数据，观测周期长，周期变化不明显，故在实际计算中仅保留长期趋势项进行拟合。最终共有 个重力测点参与计算，其中河北地区有 个。部分高值区重力测点拟合情况如图所示，年重力平均年变率情况如图所示。图重力测点拟合结果 图（）（）依次按测点地理位置由北至南排序，平均年变率范围为 。图（）中测点拟合相关系数最小，其余个测点均在 以上，拟合效果较好。重力值均呈现整体
8、上升的累积态势。从图可以看出，河北北部、西部区域重力平均年变率较小，多在 区间内变图重力平均年变率 化，居多，重力场多为季节性等短期变化，长期趋势变化无明显异常；而东南部天津邯郸一 线 平 均 年 变 率 高 值 显 著，最 大 值 可 达 ，重力持续正向累积。高值区域整体走向近 向，与累积重力变化图中高值区域形态及位置大致相符，但整体范围相对缩小。区域内年变率拟合相关系数 居多，拟合结果较好，这些区域一直处于较规律的变动状态。为排除测量期间测点发生变动对重力平差结果的影响，对河北东南部重力平均年变率高值区域内重力测点逐一进行点位变动情况核实。经统计，共核实重力测点 个，在 年测量期间点位发生
9、变动的测点为个。基于其发生变动前的部分观测数据进行拟合，然后与全部数据拟合结果对比，发现两者重力平均年变率差值在 之间。根据已有研究，河北地区重力变化主要受地表沉降及地下水流失的影响。地表沉降会使区域重力值上升，地下水流失则会使区域重力值下降。年河北东南部的天津邯郸一线重力累积变化及平均年变率均呈现明显高值，表明地表沉降的影响覆盖了地下水流失的影响，过度使用地下水的情况可能有所改善，与以往研究结果相符。为进一步研究地下水对重力场的影响，由北至南选取北京南部的河间观测井（井）、石家庄观测井（井）以及石家庄南部永年观测井（井）的水位数据（具体井位置见图）与观测井最近的重力测点（测点、）数据进行对比
10、，结果如图所示。从图（）中看出，在 年口观测井的水位均呈不同程度的下降趋势，造成重力值的下降；而图（）中重力点值呈现上下波动，大 地 测 量 与 地 球 动 力 学 年月与水位的持续下降无明显对应。故短期来看，地下水流失对重力的影响相对被覆盖。在长期变化趋势方面，除测点呈现轻微下降趋势外，测点、均为整体上升趋势。由于河北地区地势多平坦，对于观测井附近一定区域，可以按照无限平面近似公式计算其水位变化对重力场的影响：（）（）式中，为给水密度，为水位变化。根据已有研究，观测井的取 ，观测井和的均取 。若不考虑其他影响因素，将重力平均年变率全部归于地下水流失及地表沉降的影响，则个观测井的计算结果如表所
11、示。图水位和重力测点变化 表水位变化及地表沉降影响 选取地点重力平均年变率（）水位变化地表沉降（）（）（）（）井区域 井区域 井区域 由表可见，除井水位变化对重力平均年变率的影响稍大外，其他口井对重力平均年变率的影响均相对较小；地表沉降对重力平均年变率的影响明显高于地下水流失。根据已有地表沉降速率研究结果，天津邯郸一带重力场变化高值区多位 于 沉 降 中强 发 育 地 区，沉 降 速 率 在 以上，其中强发育地区达到 以上，仅 地 表 沉 降 则 至 少 可 以 引 起 重 力 场 的变化，强发育地区至少可以引起重力场 的变化。而重力平均年变率计算结果中，该部分区域重力变化多在 ，两者具有较好
12、的一致性。结语）为提高数据处理的精度，除本测区观测资料外，另外添加周边测区相对重力观测资料及绝对重力值进行整体平差控制，使最终点值精度控制在 以内，数据结果更加精准；在计算年变率时，首先对测点进行筛选，使结果能更好地反映重力场的真实变化。）基于长时间的观测结果，可以分析得出重力场长期的背景趋势变化情况。在重力累积动态变化图中，河北东南部天津邯郸一线重力场正异常持续累积，高值区逐渐扩大，最终在 年底形成走向近 的椭圆形高值异常区；北至燕山附近，西至太行山，正负差异最高可达 。在重力平均年变率结果中，正高值区域分布与重力累积动态变化结果基本一致，且该区域拟合相关系数较高。说明除去季节性变化等短期、
13、局部影响因素外，重力变化呈现长期规律上升态势。）综合上述结论及河北地区高值区域内选取的部分水位数据、相关地表沉降研究结果得出，相对于地表沉降因素，重力场受地下水流失的影响较小，综合表现为地表沉降引起的正异常累积，沉降发育严重区域影响至少可达 。参考文献张娜，龚燕民，苏树朋，等 年河北省流动重力资料整理地震地磁观测与研究，（）：（，第 卷第期杨雅慧等：河北地区 年重力场动态变化特征 ，（）：）刘洪良，王青华，张展伟，等河北唐山古冶 级地震前的重力 变 化 华 南 地 震，（）：（，（）：）曹金国，王来鹏，翟广卿，等 重力仪及应用北京：解放军出版社，（，：，）邢乐林，李辉，玄松柏，等 和地面重力测量监测到的中国大陆长期重力变化地球物理学报，（）：（，（）：），（）：，（）：李文鹏，王龙凤，杨会峰，等华北平原地下水超采状况与治理对策建议中国水利，（）：（，（）：）郭海朋，李文鹏，王丽亚，等华北平原地下水位驱动下的地面沉降现状与研究展望水文地质工程地质，（）：（，（）：）贾民育，游泽霖，万素凡，等地下水活动对精密重力测量的影响及排除方法地壳形变与地震，（）：（，（）：）贾金生，刘昌明华北平原地下水动态及其对不同开采量响应的计算：以河北省栾城县为例地理学报，（）：（，：，（）：），y ，：，；，：；：，；，：，：