GNSS实时卫星钟差估计在地震监测中的应用.pdf
《GNSS实时卫星钟差估计在地震监测中的应用.pdf》由会员分享,可在线阅读,更多相关《GNSS实时卫星钟差估计在地震监测中的应用.pdf(9页珍藏版)》请在咨信网上搜索。
1、第1 0卷 第3期2 0 2 3年5月导航定位与授时N a v i g a t i o nP o s i t i o n i n ga n dT i m i n gV o l.1 0N o.3M a y2 0 2 3d o i:1 0.1 9 3 0 6/j.c n k i.2 0 9 5-8 1 1 0.2 0 2 3.0 3.0 1 2G N S S实时卫星钟差估计在地震监测中的应用王浩浩,郝 明,庄文泉(中国地震局第二监测中心,西安7 1 0 0 5 4)摘 要:为快速、有效地获取地震发生阶段震源周边地区站点的动态位移,为地震预警系统提供高可靠性的地表形变信息,利用全球导航卫星系统(g
2、l o b a ln a v i g a t i o ns a t e l l i t es y s t e m,GN S S)高频观测数据,基于非差估计法对多模GN S S卫星钟差进行实时估计及性能分析,并将其应用于精密单点定位(p r e c i s ep o i n tp o s i t i o n i n g,P P P)实时计算2 0 2 1年漾濞Mw6.4地震和玛多Mw7.4地震的地面动态形变。结果表明,GN S S四系统实时估计卫星钟差的标准差(s t a n d a r dd e v i a t i o n,S T D)均值为0.1 4 2n s,其多系统组合P P P动态解的
3、平均标准差在水平方向达到0.5c m,高程方向达到1.0c m,计算得到的地震动态位移波形相对G P S单系统更为稳定,而且能够获得较为准确的同震形变。关键词:全球导航卫星系统(GN S S);卫星钟差;实时估计;精密单点定位;地震监测中图分类号:P 2 2 8 文献标志码:A 文章编号:2 0 9 5-8 1 1 0(2 0 2 3)0 3-0 1 0 8-0 9A p p l i c a t i o no fG N S Sr e a l-t i m e s a t e l l i t e s c l o c ko f f s e te s t i m a t i o n i ns e i
4、s m i cm o n i t o r i n gWAN G H a o h a o,HAO M i n g,Z HUAN G W e n q u a n(T h eS e c o n dM o n i t o r i n ga n dA p p l i c a t i o nC e n t e r,C E A,X i a n7 1 0 0 5 4,C h i n a)A b s t r a c t:I no r d e r t oq u i c k l ya n de f f e c t i v e l yo b t a i n t h ed y n a m i cd i s p l a
5、c e m e n t o f t h e s t a t i o n s a r o u n dt h ee p i c e n t e rd u r i n gt h ee a r t h q u a k e,a n dp r o v i d eh i g h l yr e l i a b l es u r f a c ed e f o r m a t i o ni n f o r m a t i o nf o r t h ee a r t h q u a k ee a r l yw a r n i n gs y s t e m,m u l t i-GN S Sr e a l-t i m
6、es a t e l l i t e sc l o c ko f f s e t i se s t i m a t e da n dt h ep e r f o r m a n c e i s e v a l u a t e db a s e do n t h eu n d i f f e r e n c e dm o d e l u s i n gg l o b a l n a v i g a t i o ns a t e l l i t es y s t e m(GN S S)h i g h-r a t eo b s e r v a t i o nd a t a,w h i c ha r
7、e a p p l i e d t op r e c i s ep o i n t p o s i t i o n i n g(P P P)t oc a l c u l a t e t h es u r f a c ed y n a m i cd i s p l a c e m e n to f t h e2 0 2 1Y a n g b iMw6.4e a r t h q u a k ea n dM a d o iMw7.4e a r t h q u a k e i nr e a l-t i m e.R e s u l t s s h o w s t h a t t h e a v e r
8、a g e ds t a n d a r dd e v i a t i o n s(S T D)o f t h e r e-a l-t i m ee s t i m a t e ds a t e l l i t ec l o c ko f f s e t so ft h ef o u rGN S Ss y s t e m sa r e0.1 4 2n s,t h ep o s i t i o n i n gp e r f o r m a n c eo fm u l t i-GN S SP P Pk i n e m a t i cs o l u t i o ns h o w st h a ta
9、b o u t0.5c mc a nb ea c h i e v e di nh o r i z o n t a l d i r e c t i o na n dg e n e r a l l y1.0c mi nv e r t i c a ld i r e c t i o n.M e a n w h i l e,t h ec a l c u l a t e ds e i s m i cd y n a m i cd i s p l a c e m e n tw a v e f o r m so f t h em u l t-GN S Ss y s t e ma r em o r e s t a
10、 b l e t h a n t h o s eo fG P Ss i n-g l es y s t e m,a n dc a no b t a i nm o r ea c c u r a t ec o s e i s m i cd e f o r m a t i o n.K e yw o r d s:G l o b a ln a v i g a t i o ns a t e l l i t es y s t e m(GN S S);S a t e l l i t ec l o c ko f f s e t;R e a l-t i m ee s t i m a-t i o n;P r e c
11、i s ep o i n tp o s i t i o n i n g;E a r t h q u a k em o n i t o r i n g 收稿日期:2 0 2 2-1 0-1 1;修订日期:2 0 2 3-0 2-0 1基金项目:科技部第二次青藏高原综合科学考察研究子专题(2 0 1 9 Q Z KK 0 9 0 1)作者简介:王浩浩(1 9 9 5-),男,硕士,助理工程师,主要从事G N S S精密数据处理方面的研究。通信作者:郝明(1 9 8 2-),男,博士,研究员,主要从事地壳形变监测与动力学机理方面的研究。第3期G N S S实时卫星钟差估计在地震监测中的应用0 引言近
12、年来,随着全球导航卫星系统(g l o b a ln a v i-g a t i o ns a t e l l i t es y s t e m,GN S S)高频定位技术的不断发展,其在地壳运动监测领域发挥着越来越重要的作用,为地震监测提供了一种有效的新型技术途径1-2。其中,利用差分相对定位技术可以计算得到监测站相对于某一固定参考站毫米级的动态位移,但监测区域内有时因地质条件等客观因素难以布设观测环境良好、稳定的基准站,而且坐标精确已知的基准站在强震发生阶段可能会发生移动,导致动态解算的定位精度显著降低3-4。而采用国际GN S S服务(i n t e r n a t i o n a lG
13、N S Ss e r v i c e,I G S)组织发布的精密卫星轨道和钟差等产品的精密单点定位(p r e c i s ep o i n tp o s i t i o n i n g,P P P)技术,具有不依赖于某一特定参考基准站、实时性强等优势,仅利用单台GN S S接收机即可获得国际地球参考框架下高精度的“绝对位置”,具备准确捕捉地震位移波形的能力,更 适 合 长 距 离、大 范 围 的 地 壳 形 变 监测3-4。文献5 实验结果表明,卫星钟差的采样率越高,利用P P P动态解捕捉远场站点形变信息的优势越 明 显。文 献 6 利 用 武 汉 大 学 自 主 研 发 的P AN D
14、A软件对高频GN S S观测信号进行P P P后处理,能够很好地获取2 0 1 6-1 1-1 3新西兰Mw7.8地震产生的动态位移特征。地震瞬时地表动态位移的实时高可靠性监测对地震预警系统而言至关重要,能够为震中以及震级的快速确定等研究工作提供关键信息3。然而,实时P P P高精度动态解的实现取决于精密卫星轨道和钟差的实时估计性能。目前超快速实时预报轨道产品已经能够满足其时效性和可靠性等要求,但卫星钟差的精确预报极易受到自身时频特性以及复杂太空环境的影响,高可靠的厘米级GN S S实时卫星钟差则需利用地面站观测数据进行实时估计得到7。文献8 实时估计的多模GN S S卫星钟差与武汉大学最终精
15、密钟差互差优于0.2n s。文献9-1 0 利用均方根信息滤波实现GN S S四系统实时卫星钟差的快速估计,且与最终精密钟差产品符合性较好。文献1 1 提出了一种基于序贯最小二乘的在线质量控制方法,GN S S四系统实时估计卫星钟差的标准差均值优于0.1n s。鉴于此,本文以长安大学北斗分析中心为平台,基于多模全球导航卫星系统实验跟踪网(m u l t i-GN S Se x p e r i m e n t,MG E X)监测站1H z的GN S S观测数 据,采 用 模 型 严 密、精 度 较 高 的 非 差 估 计法1 2,实现了多模GN S S卫星钟差的实时估计和性能评估,并将其用于高频
16、动态P P P实时获取2 0 2 1年漾濞Mw6.4地震和玛多Mw7.4地震发生时段的地表形变波形,具体分析震时点位的运动变化情况。1 数据处理1.1 钟差估计本文采用无电离层组合观测值进行GN S S卫星钟差的实时估计,非差伪距和载波相位的观测方程为vl c=tr-ts+br-bs+Nsr+msrTr+l cvp c=tr-ts+dr-ds+msrTr+p c(1)式中,vl c和vp c分别为相位和伪距的观测值残差;r和s分别表示接收机和卫星;tr和ts分别表示接收机和卫星的钟误差;br和bs分别表示接收机端和卫星端的相位硬件延迟偏差;dr和ds分别表示接收机端和卫星端的伪距码延迟偏差;N
17、sr表示无电离层组合观测值的相位模糊度;Tr和msr分别为对流层天顶湿延迟以及与之对应的投影函数;l c和p c分别为载波相位和伪距观测量的综合噪声。由于公式中的偏差参数具有很强甚至完全的线性相关性,将其作为未知数进行估计,不仅会引入大量待估参数,还会减弱卫星钟差的估计精度。其中,卫星端的码偏差参数具有较高的时间稳定性,可以在观测方程中将其合并到卫星的钟误差参数中,而且相位延迟偏差可以被相应的相位模糊度参数吸收。进而可得观测模型vl c=tr-ts+Nsr+msrTr+l cvp c=tr-ts+msrTr+p c(2)式中,tr=tr+dr为接收机钟差参数,ts=ts+ds为卫星钟差参数,N
18、sr=Nsr+br-bs-(dr-ds)为模糊度参数。对于全球定位系统(g l o b a l p o s i t i o n i n gs y s t e m,G P S)、格洛纳斯卫星导航系统(g l o b a ln a v i g a t i o ns a t e l l i t es y s t e m,G L ONA S S)、伽利略卫星导航系统(G a l i l e os a t e l l i t en a v i g a t i o ns y s t e m,G a l i l e o)以及我 国 的 北 斗 卫 星 导 航 系 统(B e i D o un a v i g
19、 a t i o ns a t e l l i t es y s t e m,B D S)而言,不同卫星导航系统所采用的时空基准以及信号体制不一致,导致卫星信901 导航定位与授时2 0 2 3年5月号在GN S S接 收 机 内 部 产 生 系 统 间 偏 差(i n t e rs y s t e mb i a s,I S B)以及G L ONA S S所特有的频率间偏差(i n t e rf r e q u e n c yb i a s,I F B)1 1。若以G P S时间系统为基准,则顾及I S B/I F B参数的GN S S卫星钟差实时估计模型为vGl c=tr-ts,G+Ns,G
20、r+ms,GrTr+Gl cvRkl c=tr-ts,Rk+I F BRk-G+Ns,Rkr+ms,RkrTr+Rkl cvCl c=tr-ts,C+I S BC-G+Ns,Cr+ms,CrTr+Cl cvEl c=tr-ts,E+I S BE-G+Ns,Er+ms,ErTr+El cvGp c=tr-ts,G+ms,GrTr+Gp cvRkp c=tr-ts,Rk+I F BRk-G+ms,RkrTr+Rkp cvCp c=tr-ts,C+I S BC-G+ms,CrTr+Cp cvEp c=tr-ts,E+I S BE-G+ms,ErTr+Ep c(3)式中,G、R、C、E分别表示G P
21、S、G L ONA S S、B D S和G a l i l e o系统;G L ONA S S各颗卫星的载波相位频率可以通过频率因子k计算得到;I S BC-G=dCr-dGr;I S BE-G=dEr-dGr;I F BRk-G=dRkr-dGr。为进一步提高GN S S实时卫星钟差的估计效率,通常将G L ONA S S系统的I F B参数忽略不计,以避免对其进行参数估计而 极大增加实 时钟差解 算 的 计 算负担1 3。1.2 处理策略无电离层组合模型是GN S S精密卫星钟差估计、精密单点定位等数据处理常用的函数模型,表1总结了GN S S卫 星 钟 差 实 时 估 计 以 及 实 时
22、 动 态P P P基于该模型的数据处理策略。对于GN S S卫星钟差的实时估计,则是在观测方程中将卫星钟误差作为白噪声进行估计以免受到钟跳的影响,将卫星轨道和 测 站 坐 标 作 为 已 知 值 进 行 改 正。对 于GN S S实时动态P P P,则是将测站的位置坐标作为待估参数进行估计,将精密卫星轨道和卫星钟差作为已知值进行改正。表1 数据处理策略T a b.1 D a t ap r o c e s s i n gs t r a t e g y项目参数处理方法观测模型观测量非差无电离层伪距和相位观测值频点G P S:L 1/L 2;G L ONA S S:G 1/G 2;B D S:B 1
23、/B 2;G a l i l e o:E 1/E 5 a截止高度角7 观测量加权高度角()定权(p),p=1,3 0 2 s i n(),3 0误差改正天线相位中心和变化I G S_1 4.a t x差分码偏差(d i f f e r e n t i a l c o d eb i a s,D C B)产品改正潮汐、相位缠绕、相对论效应模型改正参数估计估计方法序贯最小二乘卫星轨道固定,超快速轨道产品卫星钟差钟差估计:估计,白噪声P P P:固定,精密卫星钟差测站位置钟差估计:固定,I G S周解P P P:估计,白噪声接收机钟差估计,白噪声对流层延迟S a a s t a m o i n e n
24、模型和GMF映射函数1 4I S B常数估计,零均值条件约束相位模糊度常数估计,连续弧段为一个常数2 多模G N S S实时卫星钟差估计分析在MG E X中选取6 0个均匀分布在全球的连续运行跟踪站。从I G S下载2 0 2 1年5月1 9日2 0 2 1年5月2 2日连续4d所选测站1s采样间隔的GN S S观测数据用于卫星钟差的实时估计。在实时011 第3期G N S S实时卫星钟差估计在地震监测中的应用估计多模GN S S卫星钟差的过程中,测站始终处于静止状态且位置已知,可以将测站坐标固定到I G S周解。现阶段,超快速实时预报轨道产品与最终精密轨道产品的精度量级相当,卫星轨道可以固定
25、为武汉大学提供的包含G/R/E/C四系统轨道的6h超快速解,以减少待估参数的个数。同时,选择某一测站的接收机钟差作为参考基准钟进行约束1 5,基于双频非差消电离层组合的载波相位和伪距观测值,采用序贯最小二乘平差和验后残差质量控制算法,最终估计得到历元间隔为5s的多模GN S S实时卫星钟差。2.1 估计时间分析当利用6 0个测站构成的地面跟踪站网进行多模GN S S卫星钟差实时估计时,每个历元大约需要处理30 0 0个观测值并对将近20 0 0维的矩阵进行求逆。每个历元需要处理的观测值数目和法方程维数如图1所示。为加快高维矩阵的运算速度,提高钟差估计算法的计算效率,本文采用L A P A C
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- GNSS 实时 卫星 估计 地震 监测 中的 应用
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。