关于相对重力测量数据处理的探讨.pdf
《关于相对重力测量数据处理的探讨.pdf》由会员分享,可在线阅读,更多相关《关于相对重力测量数据处理的探讨.pdf(6页珍藏版)》请在咨信网上搜索。
1、第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关
2、于相对重力测量数据处理的探讨易天阳1 郑 兵1 李亚伦1 温军军1 易松泉1 张 稔11 四川省地震局地壳形变观测中心,四川省雅安市上坝路1 3 9号,6 2 5 0 0 0摘 要:按照 地震重力测量规范 相关理论和细则,对目前所采用的“重力测量电子记簿系统”数据处理软件存在、实际工作中急需解决的几个问题,如仪器高改正、段差计算、环闭合差计算等进行讨论。关键词:格值表;气压改正;高度改正;潮汐改正;零漂改正中图分类号:P 3 1 2 文献标识码:A 相对重力测量是指测站间的重力差测量,可采用动力法和静力法,现在普遍采用基于静力法的弹簧重力仪测定重力差值。重力测量包括野外数据采集和内业平差计算两
3、个方面,中国地震局各野外测量队外业数据采集所使用的“重力测量电子记簿系统”由中国地震局地震研究所研制,可用于地面重力仪各种观测数据的记录,包括数据采集、限差校核、实时预处理;而外业电子记簿系统导出的交换文件为内业平差计算所需的数据格式,为最终的数据处理作准备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
4、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
5、手簿工程文件到平差交换文件根据“重力测量电子记簿系统”平差交换文件(表2)输出的数据内容,结合参与重力计算的相关观测数据并与表1比较发现,最终会参与到归化改正的关键数据中缺失仪器高、经纬度及海拔高度。虽然经纬度和海拔高度可根据点号联合点之记文件计算,但仪器高不可或缺,而与归化改正无关的环境温度值却被要求记录。在一个观测点进行多次(多周期)观测仪器安置时,仪器顶面到观测台的顶面高度不一定完全相同,而且多年后可能由于观测台倾斜或顶面不平整需要添加托盘支架才能水平安置仪器,造成观测周期间同一测点的仪器安置高度不一定相同。另外,同一小组、同一观测周期使用多台仪器进行观察,仪器间的高差也可能不一样。为了
6、进行零漂率计算和不同周期间数据的比较分析,将每一次观测值改算到观测台顶面非常必要,因此量测仪器高度是必不可少的。表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
7、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)的计算,即是在“改后段差”中将潮汐改正、气压改正、高度
8、改正、零漂改正这几项均纳入了一次项系数的计算。这种做法是不合理的,因为仪器的一次项系数仅与观测仪器有关,而与潮汐、气压、高度、零漂改正无关。例如已知如下数据: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 e
9、d 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
10、 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
11、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为计算
12、时刻太阳地心天顶距;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
13、 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
14、 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 仪器读数转重力值依据仪器制
15、造厂家给定的仪器格值表(单位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=
16、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
17、(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时刻按式(
18、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 计算实例按照规范要求,经过固体潮改正后进行漂移率的计算,但往
19、返测各测点的重力值不应该是仅进行了固体潮和零漂改正之后的重力值。因为流动周期观测任务中,同一测点在不同时间进行观测,其气压值和仪器高度不一定完全相同,各测点的重力值还应归算到测点同一高度和经过气压改正之后的重力值,这样往返测量结果的较差才有意义。表5是对一个测段且往返测在同一天完成的简单重力测量线路进行零漂计算,计算的段差结果与表3一致。虽然这种做法不会对测点间的段差大小产生较大影响,但会明显影响测点的重力值计算(表6,单位mG a l)。如果测线涉及多个测段在多天完成,或同一天同一测点二次设站,或线路中较常规计划线路有支线或新增测点,则线路拟合零漂率的计算较为复杂。表6 段差(往返同天)计算
20、比较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
21、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个测点初测重力值之差获得后续测点的归一重力值
22、,在当日起点时间基础上加上测点当日时间与当天第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气压改正预处理值改
23、正后值改后段差#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 1
24、20 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.
25、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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 关于 相对 重力 测量 数据处理 探讨
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。