关于相对重力测量数据处理的探讨.pdf

第4 3卷 第9期2 0 2 3年9月大 地 测 量 与 地 球 动 力 学J o u r n a l o fG e o d e s ya n dG e o d y n a m i c sV o l.4 3N o.9S e p t.,2 0 2 3收稿日期:2 0 2 2-1 0-3 0第一作者简介:易天阳,高级工程师,主要从事地震监测研究,E-m a i l:s u n s k y-y i 1 6 3.c o m。D O I:1 0.1 4 0 7 5/j.j g g.2 0 2 3.0 9.0 1 7文章编号:1 6 7 1-5 9 4 2(2 0 2 3)0 9-0 9 7 4-0 6关于相对重力测量数据处理的探讨易天阳1 郑 兵1 李亚伦1 温军军1 易松泉1 张 稔11 四川省地震局地壳形变观测中心,四川省雅安市上坝路1 3 9号,6 2 5 0 0 0摘 要:按照 地震重力测量规范 相关理论和细则,对目前所采用的“重力测量电子记簿系统”数据处理软件存在、实际工作中急需解决的几个问题,如仪器高改正、段差计算、环闭合差计算等进行讨论。关键词:格值表;气压改正;高度改正;潮汐改正;零漂改正中图分类号:P 3 1 2 文献标识码:A 相对重力测量是指测站间的重力差测量,可采用动力法和静力法,现在普遍采用基于静力法的弹簧重力仪测定重力差值。重力测量包括野外数据采集和内业平差计算两个方面,中国地震局各野外测量队外业数据采集所使用的“重力测量电子记簿系统”由中国地震局地震研究所研制,可用于地面重力仪各种观测数据的记录,包括数据采集、限差校核、实时预处理;而外业电子记簿系统导出的交换文件为内业平差计算所需的数据格式,为最终的数据处理作准备1。对重力网进行计算之前,必须对重力仪读数进行归化和改正。对与固体潮、极移和大气压力变化有关的重力瞬时变化,可借助模型计算。就数据采集而言,“重力测量电子记簿系统”所记录的关键信息(表1)非常完善,但是在数据处理上存在值得讨论的地方。表1 数据记录信息T a b.1 D a t ar e c o r d i n f o r m a t i o n仪器号点号点名时间读数值/mG a l气压/h P a仪器高/mm海拔/m经度纬度振幅G 1 1 4 95 1 0 5 2 5 0 3#0 0 12 0 2 1-0 6-0 71 2:0 122 9 1.4 1 708 4 82 1 015 6 1.60G 1 1 4 95 1 0 5 2 4 0 2#0 0 22 0 2 1-0 6-0 71 2:4 523 1 3.4 9 408 6 42 1 313 3 6.00G 1 1 4 95 1 0 5 2 5 0 3#0 0 32 0 2 1-0 6-0 71 4:1 322 9 1.4 6 908 4 62 1 115 6 1.601 手簿工程文件到平差交换文件根据“重力测量电子记簿系统”平差交换文件(表2)输出的数据内容,结合参与重力计算的相关观测数据并与表1比较发现,最终会参与到归化改正的关键数据中缺失仪器高、经纬度及海拔高度。虽然经纬度和海拔高度可根据点号联合点之记文件计算,但仪器高不可或缺,而与归化改正无关的环境温度值却被要求记录。在一个观测点进行多次(多周期)观测仪器安置时,仪器顶面到观测台的顶面高度不一定完全相同,而且多年后可能由于观测台倾斜或顶面不平整需要添加托盘支架才能水平安置仪器,造成观测周期间同一测点的仪器安置高度不一定相同。另外,同一小组、同一观测周期使用多台仪器进行观察,仪器间的高差也可能不一样。为了进行零漂率计算和不同周期间数据的比较分析,将每一次观测值改算到观测台顶面非常必要,因此量测仪器高度是必不可少的。表2 平差交换文件T a b.2 D o c u m e n t so fa d j u s t m e n t e x c h a n g e仪器号 测区名时间点号读数值/mG a l气压/h P a温度/G 1 4 9 G C YM2 0 2 1-0 6-0 71 2:0 05 1 0 5 2 5 0 322 9 1.4 8 4 82 1G 1 4 9 G C YM2 0 2 1-0 6-0 71 2:4 45 1 0 5 2 4 0 223 1 3.5 8 6 42 2G 1 4 9 G C YM2 0 2 1-0 6-0 71 4:1 35 1 0 5 2 5 0 322 9 1.5 8 4 62 02 数据的实时预处理流动重力野外数据采集的仪器读数并不是重力值,需要利用仪器厂家给定的格值表进行转化。在格值转换(读数值到重力值的计算)过程中,按照 地震重力测量规范2的要求,格值转换后加入一次项系数的计算再进行潮汐改正,零漂改正在潮汐改正之后。而“重力测量电子记簿系统”中测线段 第4 3卷第9期易天阳等:关于相对重力测量数据处理的探讨差在“改前段差”到“改后段差”这一过程加入一次项系数(表3,单位m G a l)的计算,即是在“改后段差”中将潮汐改正、气压改正、高度改正、零漂改正这几项均纳入了一次项系数的计算。这种做法是不合理的,因为仪器的一次项系数仅与观测仪器有关,而与潮汐、气压、高度、零漂改正无关。例如已知如下数据:1)仪器号:G 1 1 4 9;零漂值:0.0 0 00 0;零漂率:-0.0 0 69 1 7;2)一次项系数:1.0 0 00 3 8;3)格值表(表4,单位m G a l),可结合表1观测数据完成相关计算。表3 观测手簿完成(往返同天)的重力联测计算T a b.3 C a l c u l a t i o no fg r a v i t y j o i n t s u r v e yc o m p l e t e d i no n ed a y点名日期读数/潮汐改正/高度改正/零漂改正/改前段差/编号时间mG a l气压改正预处理值改正后值改后段差#0 0 12 0 2 1-0 6-0 722 9 1.4 1 70/0.1 2 96/0.0 6 48/5 1 0 5 2 5 0 31 2:0 123 4 2.2 6 920.0 0 2623 4 2.4 6 6223 4 2.4 6 62#0 0 22 0 2 1-0 6-0 723 1 3.4 9 40/0.1 1 60/0.0 6 57/0.0 0 51/2 2.5 6 70/5 1 0 5 2 4 0 21 2:4 523 6 4.8 4 600.0 0 0423 6 5.0 2 8023 6 5.0 3 302 2.5 6 79#0 0 12 0 2 1-0 6-0 722 9 1.4 6 90/0.0 6 15/0.0 6 51/0.0 1 52/5 1 0 5 2 5 0 31 4:1 323 4 2.3 2 240.0 0 2023 4 2.4 5 1023 4 2.4 6 60 注:后期纳入一次项系数,改前段差一次项系数=改后段差,表中2 2.5 6 701.0 0 00 3 8=2 2.5 6 79。表4 格值T a b.4 G r i dv a l u e t a b l e计数器读数值重力值内插因子22 0 0.0 022 4 8.7 8 81.0 2 25 823 0 0.0 023 5 1.0 4 61.0 2 26 82.1 潮汐改正计算重力固体潮汐理论值时,应使用测点经、纬度值并结合当前日期、时间,具体可采用直接法或分波叠加法。采用直接法较为简单,从未作谐波展开的引潮位公式出发,计算重力固体潮理论值。计算公式为3:g=-1 6 5.7F()CR 3c o s2Z-13 -1.3 7F2()CR 4(5 c o s2Z-3)-7 6.0 8F()CsRs 3c o s2Zs-13 (1)式中,F()=0.9 9 83 2 7+0.0 0 16 7 6c o s 2,为计算点的地理纬度;Z为计算时刻月亮地心天顶距;Zs为计算时刻太阳地心天顶距;CR为月地平均距离和瞬时距离之比;CsRs为日地平均距离和瞬时距离之比。为计算g,必须使用5个有关的天文参数:月亮平黄经S,太阳平黄经H,月亮近地点经度p,月亮升交点经度N,近日点经度Ps,这5个参数分别有各自的涉及儒略世纪数T的计算公式,具体见文献3。表1中5 1 0 5 2 5 0 3测点1 2:0 1时刻的潮汐改正值依据式(1)计算为0.1 3 08mG a l(表5,单位mG a l),与表3潮汐改正数0.1 2 96mG a l略有差异,可能是由于“重力测量电子记簿系统”所用算法不是直接算法而造成的。表5 重力观测(往返同天)数据改正计算T a b.5 C o r r e c t i o nc a l c u l a t i o no fg r a v i t yo b s e r v a t i o nd a t ac o m p l e t e d i no n ed a y点号时间读数重力值潮汐改正气压改正仪器高改正归一时顺归一重力dtdR零漂改正改正后重力5 1 0 5 2 5 0 32 0 2 1-0 6-0 71 2:0 122 9 1.4 1 7 23 4 2.3 5 82 0.1 3 08 0.0 0 26 0.0 6 48 1 2.0 1 67 23 4 2.5 5 64023 4 2.5 5 645 1 0 5 2 4 0 22 0 2 1-0 6-0 71 2:4 523 1 3.4 9 4 23 6 4.9 3 59 0.1 1 72 0.0 0 04 0.0 6 57 1 2.7 5 00 23 6 5.1 1 920.0 0 51 23 6 5.1 2 435 1 0 5 2 5 0 32 0 2 1-0 6-0 71 4:1 322 9 1.4 6 9 23 4 2.4 1 14 0.0 6 27 0.0 0 20 0.0 6 51 1 4.2 1 67 23 4 2.5 4 12-0.0 1 522.20.0 1 52 23 4 2.5 5 64K=-0.0 0 69 0 92.2 仪器读数转重力值依据仪器制造厂家给定的仪器格值表(单位1 0G a l)将仪器读数转换为相应的重力值。重力值换算按下式进行3:g=F1+(R-R1)F2F-T D(2)式中,g为相应仪器读数的重力值,F1为仪器读数每1 0 0个分划间隔的整数所对应的重力值,R为仪器读数,R1为仪器读数以1 0 0分划为间隔的最大整数,F2为仪器格值间隔因子,F为仪器格值一次项系数,T D为潮汐改正值。表1中5 1 0 5 2 5 0 3测点1 2:0 1时刻读数转换重力值按式(2)的计算为:(22 4 8.7 8 8+(22 9 1.4 1 70-22 0 0.0 0)1.0 2 25 8)1.0 0 00 3 8+0.1 3 08=579大 地 测 量 与 地 球 动 力 学2 0 2 3年9月23 4 2.4 8 90(mG a l)而表3中数据转换重力值计算为:(22 4 8.7 8 8+(22 9 1.4 1 70-22 0 0.0 0)1.0 2 25 8)=23 4 2.2 6 92(mG a l)2.3 气压改正大气压力变化的直接和间接效应会引起重力变化。考虑到测点上实际大气压p(气压计测量的精度为1h P a)和正常大气压pn之间的差异,应采用下式进行归一化3:gp=3(p-pn)(h P a)n m s-2(3)式中,pn为根据测站高度H计算的标准大气压值:pn=10 1 3.2 51-0.0 0 65H(m)2 8 8.1 5 5.2 5 59h P a(4)表1中5 1 0 5 2 5 0 3测点1 2:0 1时刻按式(4)进行气压改正:gp=3(8 4 8-10 1 3.2 51-0.0 0 6515 6 1.6 02 8 8.1 5 5.2 5 59 n m s-2=2 6.4 3n m s-2=0.0 0 26mG a l表3与表5的气压改算结果一致。2.4 高度改正正常重力和高程的关系通常使用垂向导数描述,展开到f阶并归化到=4 5 为3:h 0=-30 8 6n m s-2/m=-3.0 8 6mG a l/m(5)表1中5 1 0 5 2 5 0 3测点1 2:0 1时刻按式(5)进行高度改正:-30 8 6n m s-22 1 0mm10 0 0mm=-6 4 8n m s-2=-0.0 6 48mG a l表4与表5的高度改算结果一致。2.5 零漂改正要完成零漂改正,首先需要完成重力仪漂移率的计算,而重力漂移率的计算方法可采用全网最或然漂移率,亦可求取测线拟合漂移率。目前较为常用的是对测线进行漂移率计算,具体公式为3:K=dRidtidtidti=(g i-gi)(t i-ti)(t i-ti)2(6)式中,K为漂移率,dRi为第i点经固体潮改正的往返观测重力值之差,dti为第i点往返观测时间之差。3 计算实例按照规范要求,经过固体潮改正后进行漂移率的计算,但往返测各测点的重力值不应该是仅进行了固体潮和零漂改正之后的重力值。因为流动周期观测任务中,同一测点在不同时间进行观测,其气压值和仪器高度不一定完全相同,各测点的重力值还应归算到测点同一高度和经过气压改正之后的重力值,这样往返测量结果的较差才有意义。表5是对一个测段且往返测在同一天完成的简单重力测量线路进行零漂计算,计算的段差结果与表3一致。虽然这种做法不会对测点间的段差大小产生较大影响,但会明显影响测点的重力值计算(表6,单位mG a l)。如果测线涉及多个测段在多天完成,或同一天同一测点二次设站,或线路中较常规计划线路有支线或新增测点,则线路拟合零漂率的计算较为复杂。表6 段差(往返同天)计算比较T a b.6 C o m p a r i s o no f s e g m e n td i f f e r e n c ec a l c u l a t i o n s c o m p l e t e d i no n ed a y点名F=1.0 0 00 3 8先期纳入一次项系数后期纳入一次项系数改后重力段差改前重力改前段差改后重力改后段差重力值较差段差较差5 1 0 5 2 5 0 323 4 2.5 5 6423 4 2.4 6 6223 4 2.5 5 520.0 0 125 1 0 5 2 4 0 223 6 5.1 2 432 2.5 6 7923 6 5.0 3 322 2.5 6 7023 6 5.1 2 312 2.5 6 790.0 0 120 以往返测不同天完成的测线为例,表7(单位mG a l)中5 1 0 2 6 4 0 0测点结束于1 1日下午1 6:0 1,1 2日上午8:2 3接着该点出发,2次读数肯定不同,经过格值转换的重力值也不一定完全相等,即便经过潮汐改正、气压改正、高度改正,重力值也不一定完全相等。但是该点这2d的重力值理论上应该相等,因此将第2天的开始时间归算到前1天的结束时间,将第2天的开始重力值归算到前一天的结束重力值,在当日起点归算重力值的基础上加上当日测线上测点初测重力值与当天第1个测点初测重力值之差获得后续测点的归一重力值,在当日起点时间基础上加上测点当日时间与当天第1个测点起测时间之差获得归一时间,具体见表8。表8为同一测点Ni的2次(同天或不同天观测)设站情况。一般而言,2次设站读数值经过679 第4 3卷第9期易天阳等:关于相对重力测量数据处理的探讨表7 观测手簿完成(往返不同天)的重力联测计算T a b.7 C a l c u l a t i o no fg r a v i t y j o i n t s u r v e yc o m p l e t e d i nd i f f e r e n td a y s点名日期读数/潮汐改正/高度改正/零漂改正/改前段差/编号时间mG a l气压改正预处理值改正后值改后段差#0 1 02 0 2 1-0 6-1 120 2 1.3 6 00/0.1 5 17/0.0 6 57/5 1 0 3 4 0 0 11 2:4 320 6 6.1 3 750.0 0 2420 6 6.3 5 7320 6 6.3 5 73#0 1 12 0 2 1-0 6-1 120 1 5.5 2 00/0.1 6 28/0.0 6 48/0.0 0 19/-5.9 6 60/5 1 0 2 6 4 0 01 3:3 220 6 0.1 6 660.0 0 2420 6 0.3 9 6720 6 0.3 9 86-5.9 6 62#0 1 22 0 2 1-0 6-1 120 8 1.4 7 90/0.1 5 00/0.0 6 51/0.0 0 46/6 7.4 2 97/5 1 0 2 6 3 0 01 4:4 121 2 7.6 0 310.0 0 2721 2 7.8 2 0921 2 7.8 2 556 7.4 3 22#0 1 12 0 2 1-0 6-1 120 1 5.5 7 00/0.1 0 02/0.0 6 54/0.0 0 76/5 1 0 2 6 4 0 01 6:0 120 6 0.2 1 780.0 0 2120 6 0.3 8 5520 6 0.3 9 31#0 1 12 0 2 1-0 6-1 220 1 5.7 1 40/-0.0 6 28/0.0 6 63/0.0 0 76/5 1 0 2 6 4 0 00 8:2 320 6 0.3 6 500.0 0 4520 6 0.3 7 3020 6 0.3 9 31#0 1 02 0 2 1-0 6-1 220 2 1.5 2 80/-0.0 3 49/0.0 6 63/0.0 0 96/5 1 0 3 4 0 0 10 9:1 320 6 6.3 0 920.0 0 3620 6 6.3 4 4320 6 6.3 6 63潮汐改正后的重力值不一定相等,而理论上同一个测点的重力值应该相等,因此为方便计算漂移率,需要将2次设站时间和重力值统一到上一次设站,即=,G i=G i,而后续测点的时间(注意将时分转换为时小数计算):=+-,重力值:Gi+1=G i+G3+(7)+(8)+(9)-G2-(4)-(5)-(6)。依据 表8所 描 述 的 计 算 方 法,表9(单 位mG a l)中5 1 0 3 4 0 0 1测点返测相关计算如下:1)5 1 0 3 4 0 0 1测点返测观测的归一时顺时间为:1 6.0 1 67+9+1 3/6 0-8-2 3/6 0=1 6.8 5 00(h)。2)5 1 0 3 4 0 0 1测点返测观测的归一重力为:20 6 0.4 6 45+20 6 6.3 8 77-0.0 3 43+0.0 0 36+0.0 6 63-20 6 0.4 4 33+0.0 6 23-0.0 0 45-0.0 6 63=20 6 6.4 3 60(mG a l)。依据式(6),表9中的漂移率K值计算为:K=(2.4 8 34(-0.0 1 11)+4.1 3 33(-0.0 0 05)/(2.4 8 342+4.1 3 332)=-0.0 0 12 7(m G a l/)h。将K值再乘以各个观测点与第1个观测点时间差,得到各观测点的零漂改正,最后得到各观测点的改正后重力值。先期纳入一次项系数与后期纳入一次项系数的测点重力值、测点间段差的结果有没有差别是不同方法应该比较的问题。将表9与表7的计算结果进行对比,结果见表1 0(单位mG a l)。通过表1 0可以看出,不同计算方法的段差较差基本保表8 重力观测(往返不同天)数据改正计算表模板T a b.8 T e m p l a t e f o rc o r r e c t i o nc a l c u l a t i o no fg r a v i t yo b s e r v a t i o nd a t ac o m p l e t e d i nd i f f e r e n td a y s点号年月日时分读数重力值潮汐改正 气压改正 仪器高改正 归一时顺 归一重力dtdR零漂改正改正后重力N iG1(1)(2)(3)G i(1 0)N iG2(4)(5)(6)G i(1 1)Ni+1G3(7)(8)(9)Gi+1(1 2)K=表9 重力观测(往返不同天)数据改正计算T a b.9 C o r r e c t i o nc a l c u l a t i o nT a b l e.o fg r a v i t yo b s e r v a t i o nd a t ac o m p l e t e d i nd i f f e r e n td a y s点号时间读数重力值潮汐改正气压改正仪器高改正归一时顺归一重力dtdR零漂改正改正后重力5 1 0 3 4 0 0 12 0 2 1-0 6-1 11 2:4 320 2 1.3 6 0 20 6 6.2 1 60 0.1 5 24 0.0 0 24 0.0 6 57 1 2.7 1 6720 6 6.4 3 650.0 0 00 20 6 6.4 3 655 1 0 2 6 4 0 02 0 2 1-0 6-1 11 3:3 220 1 5.5 2 0 20 6 0.2 4 49 0.1 6 34 0.0 0 24 0.0 6 48 1 3.5 3 3320 6 0.4 7 560.0 0 10 20 6 0.4 7 665 1 0 2 6 3 0 02 0 2 1-0 6-1 11 4:4 120 8 1.4 7 9 21 2 7.6 8 40 0.1 5 05 0.0 0 27 0.0 6 51 1 4.6 8 3321 2 7.9 0 230.0 0 25 21 2 7.9 0 485 1 0 2 6 4 0 02 0 2 1-0 6-1 11 6:0 120 1 5.5 7 0 20 6 0.2 9 61 0.1 0 09 0.0 0 21 0.0 6 54 1 6.0 1 6720 6 0.4 6 45 2.4 8 34-0.0 1 11 0.0 0 42 20 6 0.4 6 875 1 0 2 6 4 0 02 0 2 1-0 6-1 20 8:2 320 1 5.7 1 4 20 6 0.4 4 33-0.0 6 23 0.0 0 45 0.0 6 63 1 6.0 1 6720 6 0.4 6 450.0 0 42 20 6 0.4 6 875 1 0 3 4 0 0 12 0 2 1-0 6-1 20 9:1 320 2 1.5 2 8 20 6 6.3 8 77-0.0 3 43 0.0 0 36 0.0 6 63 1 6.8 5 0020 6 6.4 3 60 4.1 3 33-0.0 0 05 0.0 0 52 20 6 6.4 4 12K=-0.0 0 12 7779大 地 测 量 与 地 球 动 力 学2 0 2 3年9月表1 0 段差(往返不同天)计算结果比较T a b.1 0 C o m p a r i s o no f s e g m e n td i f f e r e n c ec a l c u l a t i o n s c o m p l e t e d i nd i f f e r e n td a y s点名F=1.0 0 00 3 8先期纳入一次项系数后期纳入一次项系数改后重力段差改前重力改前段差改后重力改后段差重力值较差段差较差5 1 0 3 4 0 0 120 6 6.4 3 6520 6 6.3 5 7320 6 6.4 3 580.0 0 075 1 0 2 6 4 0 020 6 0.4 7 66-5.9 5 9920 6 0.3 9 86-5.9 5 8720 6 0.4 7 69-5.9 5 89-0.0 0 03-0.0 0 105 1 0 2 6 3 0 021 2 7.9 0 486 7.4 2 8221 2 7.8 2 556 7.4 2 6921 2 7.9 0 646 7.4 2 95-0.0 0 16-0.0 0 135 1 0 2 6 4 0 020 6 0.4 6 87-6 7.4 3 6120 6 0.3 9 31-6 7.4 3 2420 6 0.4 7 14-6 7.4 3 50-0.0 0 28-0.0 0 125 1 0 3 4 0 0 120 6 6.4 4 125.9 7 2520 6 6.3 6 635.9 7 3220 6 6.4 4 485.9 7 34-0.0 0 34-0.0 0 08持在同一水平,但测点重力值较差随着作业的先后在逐渐增大,表明“重力测量电子记簿系统”所采用的后期纳入一次项系数计算会影响测点间的段差值及测点的重力值。究其原因,是测前段差包含了潮汐改正、气压改正和仪器高改正,也是“重力测量电子记簿系统”将潮汐改正、气压改正和仪器高改正加入一次项系数计算,而仪器一次项系数仅与仪器有关,将潮汐改正、气压改正和仪器高改正纳入一次项系数计算是不合适的。4 闭合环综合评价 相对重力联测资料质量评比评分细则4在环闭合差的评价方面有相关得分要求,而现有“重力测量电子记簿系统”尚不能完成闭合环的搜寻和环闭合差统计,需要增加上述功能并在野外数据成环测段进行环闭合差均值的计算。环闭合差均值的计算公式4为:A6=Nn=1|Bn|Cn/5/N(7)式中,N为闭合环个数,Bn为第n个闭合环的闭合差,Cn为第n个闭合环测段数(或节点数)。假定有不同测段数的4个闭合环,每个环闭合差及相关计算结果见表1 1(单位mG a l)。表1 1 环闭合差均值计算比较T a b.1 1 C o m p a r i s o no fm e a no f l o o pc l o s u r e e r r o r环号N测段数Cn闭合差Bn Cn/5|Bn|Cn5/NBn/Cn|Bn|/Cn5/N11 40.0 0 331.6 7 33 20.0 0 13 8 00.0 0 02 3 60.0 0 04 9 321 40.0 5 321.6 7 33 20.0 2 22 5 50.0 0 38 0 00.0 0 79 4 831 1-0.0 3 081.4 8 32 40.0 1 14 2 1-0.0 0 28 0 00.0 0 51 9 143 80.1 2 002.7 5 68 10.0 8 27 0 40.0 0 31 5 80.0 1 08 8 2A60.1 1 78 0 00.0 2 45 0 0 由表1 1可以看出,本文最后得出的环闭合差均值为0.1 1 78mG a l,其值远大于文献4 要求的最大值0.0 3 6mG a l,然而通过Bn/Cn的比值发现,每一个环内测段的平均闭合差相差并不大,其根本原因在于环内的测段数多,造成乘积的结果偏差大。假设某个测区只有1个如表1 1中的3号闭合环,1 1个测段的环闭合差0.0 3 08mG a l,按式(7)计算出的环闭合差均值为0.0 4 57mG a l;另外假设3号环只有6个测段,其环闭合差也是0.0 3 08mG a l,计算出的环闭合差均值为0.0 3 37mG a l。通过1 1个测段与6个测段的闭合环比较可以看出,在环闭合差相同的情况下,测段数多的闭合环精度较测段数少的闭合环精度差,这显然是不合适的。由此可见,文献4 中的算法会造成结果值偏大,尤其对于测段数较多的闭合环来说其评价结果是非常不利的。因此,|Bn|后的乘号更改为除号更合适,这样表1 1中4个环的环闭合差均值为0.0 2 45mG a l。5 结 语1)“重力测量电子记簿系统”中工程文件到平差交换文件缺少仪器高数据,造成野外量测仪器高似乎是多余的,然而仪器架设在每一测点或同一测点往返测的高度并不相同,为将重力值归算到观测墩台面,仪器高改正不可缺少,因此工程文件到平差交换文件需要完善仪器高的输出。2)一次项系数仅与仪器有关,而与气压改正、仪器高改正、潮汐改正、零漂改正等无关,因此一次项系数参与计算应该在格值转换后完成,而目前所采用的“重力测量电子记簿系统”后期纳入一次项系数计算会影响测点间的段差值和测点的重力值成果。3)简单线路的各项改正计算较为明确和简单,但涉及测点间歇或线路上增设测点后漂移率的计算较为复杂。879 第4 3卷第9期易天阳等:关于相对重力测量数据处理的探讨4)对文献4 中成环的环平均闭合差计算公式存疑,建议将公式中|Bn|后的乘号修正为除号。参考文献1 中国地震局地震研究所.C G-5重力仪的地震重力测量应用手册Z.2 0 0 7(I n s t i t u t eo fS e i s m o l o g y,C E A.A p p l i c a-t i o nM a n u a l o fS e i s m i cG r a v i t yM e a s u r e m e n to fC G-5G r a-v i m e t e rZ.2 0 0 7)2 国家地震局.地震重力测量规范S.北京:地震出版社,1 9 7 9(C h i n aE a r t h q u a k eA d m i n i s t r a t i o n.S p e c i f i c a t i o nf o rS e i s-m i cG r a v i t yS u r v e yS.B e i j i n g:S e i s m o l o g i c a lP r e s s,1 9 7 9)3 T o r g e W.重 力 测 量 学 M.北 京:地 震 出 版 社,1 9 9 3(T o r g eW.G r a v i m e t r yM.B e i j i n g:S e i s m o l o g i c a lP r e s s,1 9 9 3)4 中国地震局地震研究所.相对重力联测资料质量评比评分细则 Z.2 0 2 2(I n s t i t u t eo fS e i s m o l o g y,C E A.D e t a i l e dR u l e s f o rQ u a l i t yE v a l u a t i o na n dG r a d i n go fR e l a t i v eG r a v i-t yJ o i n tS u r v e yD a t aZ.2 0 2 2)D i s c u s s i o no nD a t aP r o c e s s i n go fR e l a t i v eG r a v i t yS u r v e yY IT i a n y a n g1 ZHENGB i n g1 L IY a l u n1 WENJ u n j u n1 Y IS o n g q u a n1 ZHANGR e n11 C r u s t a lD e f o r m a t i o nM o n i t o r i n gC e n t e ro fS i c h u a nE a r t h q u a k eA g e n c y,1 3 9S h a n g b aR o a d,Y aa n6 2 5 0 0 0,C h i n aA b s t r a c t:I na c c o r d a n c ew i t h t h e r e l e v a n t t h e o r i e s a n d s p e c i f i c a t i o n s o f t h eS t a n d a r df o rE a r t h q u a k eG r a v i t yS u r v e y,w e d i s c u s s s e v e r a l p r o b l e m s u r g e n t l yn e e d t ob e s o l v e d i np r a c t i c a lw o r k,a n d e x i s t i n g i n t h e c u r r e n t-l yu s e dd a t ap r o c e s s i n gs o f t w a r eo f t h eg r a v i t ys u r v e ye l e c t r o n i c r e c o r d i n gs y s t e m,s u c ha s i n s t r u m e n th e i g h tc o r r e c t i o n,s e c t i o nd i f f e r e n c e c a l c u l a t i o n,r i n gc l o s u r ed i f f e r e n c e c a l c u l a t i o n,e t c.K e yw o r d s:g r i dv a l u e t a b l e;b a r o m e t r i c c o r r e c t i o n;h e i g h t c o r r e c t i o n;t i d a l c o r r e c t i o n;z e r od r i f t c o r r e c t i o nA b o u t t h e f i r s t a u t h o r:Y IT i a n y a n g,s e n i o re n g i n e e r,m a j o r s i ne a r t h q u a k em o n i t o r i n g,E-m a i l:s u n s k y-y i 1 6 3.c o m.(上接第9 4 9页)G r a v i t yF i e l dC h a n g e s i nH a i n a nI s l a n db e f o r ea n da f t e rS a n y aMS4.2E a r t h q u a k ea n dI t sW a v e l e tM u l t i s c a l eD e c o m p o s i t i o nL IS h e n g1 HU ANGZ h a n g r o n g1 HUJ i u c h a n g1 X I EX i a o j i n g1 J I AW e i11 H a i n a nE a r t h q u a k eA g e n c y,1 3M e i y u a nR o a d,H a i k o u5 7 0 2 0 3,C h i n aA b s t r a c t:U s i n gt h er e p e a t e dm e a s u r e m e n td a t a f r o m H a i n a n i s l a n dm o b i l eg r a v i t ys u r v e yn e t w o r ko f1 3p e r i o d s f r o m2 0 1 6 t o2 0 2 2,w e a n a l y z e t h ed y n a m i c c h a n g e so f d i f f e r e n t i a l g r a v i t y f i e l da n dc u m u l a-t i v eg r a v i t y f i e l d i n t h e s t u d ya r e a,a n de m p h a t i c a l l ya n a l y z e t h e s p a t i a l a n d t e m p o r a l c h a n g e so f g r a v i-t yf i e l d i n t h e s o u t h w e s t o fH a i n a n i s l a n db e f o r e a n da f t e r t h eS a n y