内容简介
结合当前多频GNSS发展趋势,《多模多频GNSS精密单点定位理论与方法》主要对多频精密单点定位中的若干关键技术进行介绍,包括GNSS时空基准、PPP基本模型与误差处理,重点围绕PPP钟差估计、相位偏差估计、模糊度解算策略,以及区域参考站与低轨卫星增强的快速PPP等方面进行介绍。相关模型和算法优化创新,提高了PPP定位性能。
目录
目录
前言
第1章绪论1
1.1背景介绍1
1.2国内外现状2
1.2.1精密轨道和钟差产品2
1.2.2多模多频PPP函数模型4
1.2.3多模多频模糊度固定方法5
1.2.4多模多频PPP偏差估计7
1.2.5LEO增强的GNSS精密定位8
参考文献10
第2章GNSS时空基准17
2.1时间基准17
2.1.1常用时间系统17
2.1.2GNSS时间系统19
2.2空间基准20
2.2.1常用坐标系统21
2.2.2GNSS坐标系统与转换25
参考文献27
第3章精密单点定位模型与方法29
3.1GNSS精密单点定位函数模型29
3.1.1传统双频无电离层组合模型30
3.1.2统一多频无电离层组合模型32
3.1.3多个双频无电离层组合模型34
3.1.4基于原始观测值的非组合模型35
3.1.5多系统组合PPP模型37
3.2GNSS精密单点定位随机模型39
3.3参数估计策略40
3.3.1*小二乘估计40
3.3.2卡尔曼滤波估计41
3.3.3抗差卡尔曼滤波44
参考文献44
第4章精密单点定位系统误差处理46
4.1与卫星有关的误差46
4.1.1星历误差46
4.1.2相对论效应48
4.1.3卫星差分码偏差48
4.1.4频间钟差53
4.1.5BDS-2星端伪距偏差53
4.1.6GLONASS伪距IFB54
4.1.7卫星端天线改正54
4.1.8卫星姿态建模55
4.2与信号传播有关的误差63
4.2.1对流层延迟63
4.2.2电离层延迟64
4.2.3地球自转改正64
4.2.4多路径效应65
4.3与接收机有关的误差65
4.3.1接收机钟差65
4.3.2地球潮汐改正67
4.3.3接收机端天线改正68
参考文献68
第5章精密单点定位钟差估计方法72
5.1精密卫星钟差估计模型72
5.1.1基于非差模型的精密卫星钟差估计72
5.1.2基于星间单差模型的精密卫星钟差估计73
5.1.3基于组合差分模型的精密卫星钟差估计74
5.1.4卫星钟差估计实验分析75
5.2基于广播星历的实时轨道钟差综合补偿方法86
5.2.1单站综合误差计算方法86
5.2.2多站综合误差加权87
5.2.3综合误差提取与验证分析88
5.3多频钟差频间偏差(IFCB)估计102
5.3.1IFCB估计方法103
5.3.2IFCB建模方法104
5.3.3IFCB实验分析105
参考文献116
第6章多频精密单点定位小数偏差改正117
6.1FCB估计基本模型117
6.2双频无电离层模型FCB估计118
6.3双频/多频非组合模型FCB估计120
6.4不同FCB之间的转换关系122
6.4.1双频无电离层与非组合模型FCB转换123
6.4.2多频MW组合与非组合模型FCB转换125
6.5多频FCB实验分析126
6.5.1多频FCB单天稳定性128
6.5.2多频FCB验后残差分布133
6.5.3双频无电离层与非组合FCB一致性分析134
参考文献137
第7章多频精密单点定位模糊度固定138
7.1多频非组合PPP模糊度固定138
7.1.1整体固定策略139
7.1.2分步固定策略139
7.1.3两种固定策略的比较141
7.1.4固定解参数更新方法142
7.1.5多频非组合PPP固定解实验分析142
7.2多频弱电离层组合PPP模糊度固定153
7.2.1多频观测值组合153
7.2.2多频弱电离层(IR)PPP定位模型154
7.2.3基于多频IRPPP的宽巷解实验分析155
7.2.4基于多频IRPPP的窄巷解实验分析160
参考文献164
第8章区域参考站增强的精密单点定位166
8.1基准站大气信息提取方法166
8.1.1基准站大气延迟提取166
8.1.2基准站间大气延迟互检166
8.2用户站增强信息生成方法167
8.2.1线性组合模型168
8.2.2线性内插模型168
8.2.3基于距离的线性内插模型169
8.2.4低阶*面拟合模型169
8.3大气增强信息约束方法170
8.3.1常数约束法171
8.3.2时空约束法171
8.3.3逐步松弛约束法172
8.3.4误差函数法172
8.4区域参考站增强的PPP实验分析172
8.4.130s采样率数据实验分析172
8.4.21s采样率数据实验分析181参考文献186
第9章低轨卫星增强的精密单点定位188
9.1仿真数据生成方法188
9.2纯仿真LEO/GNSS数据实验分析189
9.2.1PPP浮点解结果191
9.2.2PPP固定解结果194
9.2.3定位精度对比196
9.3仿真LEO+实测GNSS数据实验分析198
9.4LEO与精密大气联合增强的PPP201
9.4.1LEO/GNSS增强信息精度评估201
9.4.2LEO与大气联合增强的PPP定位性能评估202
参考文献205
试读
第1章 绪论
1.1 背景介绍
全球导航卫星系统 (Global Navigation Satellite System,GNSS)凭借其全天候工作、全球覆盖率高、定位精度高以及实时高效性等诸多优势,成为人们获取精准时空信息的主要手段之一,被广泛应用于导航定位、变形监测、精密授时、精密定轨等诸多领域。作为国家的重大基础设施,也是国家定位导航授时 (Positioning Navigation and Timing,PNT)体系中极为重要的一环。目前,国际公认的四大 GNSS系统包括美国的 GPS、俄罗斯的 GLONASS、欧盟的 Galileo以及我国自主研发的北斗卫星导航系统 (BDS)[1-3]。近年来,随着 GPS的现代化、 GLONASS和 Galileo系统的不断完善以及我国 BDS-3系统的全球组网,新一代的导航卫星能够在三个乃至更多的频点播发导航信号,GNSS正式进入了多模多频兼容互操作的时代[4]。多系统多频观测信息为卫星定位精度、可靠性等性能的提升提供了新的技术条件,进一步激发了卫星定位新理论、新方法研究的热潮。
实时动态定位 (Real-Time Kinematic,RTK)和精密单点定位 (Precise Point Positioning,PPP)技术是 GNSS高精度定位的代表性技术[5-7]。其中,RTK是一种基于载波差分的相对定位技术,经历了由单站 RTK向网络 RTK技术的转变,通过密集的参考站对区域大气延迟进行高精度建模,可为用户提供实时厘米级的定位,该项技术较为成熟,已在国内各大省市实现工程化推广应用[8,9]。不过,由于其建设成本高以及数据处理中心的运算负荷限制,难以实现广域的高精度定位。相比之下,PPP技术通过精确考虑各种误差改正,基于单台接收机即可在全球范围内实现绝对高精度定位,具有定位灵活、便捷性高的优点。不过,其往往受限于外部的精密轨道和钟差产品,且模型待估参数较多且相关性强,通常需要较长的时间才能达到厘米级的精度,同时,由于存在小数周偏差 (Fractional Cycle Bias,FCB)的影响,常规的PPP数据处理仅能得到浮点解,定位的可靠性较低,即便在对 FCB改正的基础上,可靠的模糊度固定也需要相当的初始化时间。随着自动驾驶等新型位置服务行业的兴起,除了定位精度以外,用户更加注重定位的实时性、连续性和可靠性,尤其是在城市峡谷等复杂应用场景下,受限于GNSS信号自身的脆弱性,一旦发生周跳或信号失锁,PPP将面临着重收敛或重初始化的问题,这些难点也成为PPP规模化推广应用的瓶颈。
随着多模多频逐步成为GNSS领域新的发展趋势,相比于单系统,多模 GNSS意味着更多的可视卫星数与更优的几何分布,多频信息意味着更多特性优良的观测值组合,这些都可以加强 PPP的模型强度,一定程度上缩短其收敛和初始化时间,提高定位的连续性和可靠性,但距离实时应用还有一定差距。为了实现PPP模糊度的准实时固定,通常需要引入额外的增强信息,其中较常用的为区域参考站生成的高精度大气信息,通过附加约束或直接改正的方式,减弱模型待估参数之间的相关性,实现模糊度的瞬时固定。除此之外,随着 SpaceX Starlink、“鸿雁”等低轨(Low Earth Orbit,LEO)星座的相继提出与建设,传统的地基增强 GNSS高精度定位将面临新的挑战,相比于 GNSS中高轨道卫星,LEO卫星具有几何图形变化快的优势,理论上可进一步缩短 PPP的初始化时间,通过LEO增强GNSS实现高中低轨优势互补,形成星地联合增强的新一代 PNT体系,对增强GNSS的精度、完好性、连续性和可用性具有重要意义。
1.2国内外现状
PPP技术自1997年提出发展至今,经历了由单系统单/双频向多模多频的发展,致力于实现快速可靠的PPP,国内外相关研究机构和人员在多系统精密产品生成、多模多频解算模型、偏差估计与模糊度解算等方面均作了较为丰富的研究,同时随着LEO星座的建设发展,LEO增强的快速PPP也是目前的研究热点,为此,下面将对相关研究内容进行阐述与总结。
1.2.1精密轨道和钟差产品
实时精密轨道和钟差产品是PPP高精度定位的先决条件,其连续稳定性将极大地影响 PPP的定位性能。传统的 PPP通常采用国际 GNSS服务(International GNSS Service,IGS)各分析中心提供的昀终产品,其标称的轨道精度通常优于2.5cm,钟差精度优于 75ps,不过,其存在数十天的滞后性,无法应用于实时 PPP。为了解决精密产品的时延问题,推广 PPP的实时应用,IGS早在 2001年就成立了实时产品相关的工作小组 (Real-Time Working Group,RTWG),并于 2013年正式启动实时试验计划 (Real-Time Pilot Project,RTPP),在这之前,所谓的实时 PPP主要是基于 IGS超快速预报星历 (IGS Ultra-Rapid Products,IGU)进行,不过,由于 IGU星历的钟差产品精度较差,实际使用中,往往还需要配合区域的参考站进行自主钟差估计,因此并非严格意义的实时 PPP[10,11]。随着 RTPP的正式运行,基于广播星历状态空间表示 (State Space Representation,SSR)改正数的实时 PPP逐渐成为主流。王胜利等对 IGS实时数据流的精度进行评估,表明实时轨道精度优于 5cm,钟差精度优于 0.2ns,满足实时 PPP厘米级的定位需求[12]。随后,同济大学的刘志强等使用 IGS提供的 SSR改正数据实现了静态平面 2cm,高程 4cm的定位精度[13],武汉大学夏凤雨等采用多个挂载点的 SSR数据实现了仿动态 PPP单天解 E、N、U方向优于 20cm的定位精度[14]。此外,东南大学程良涛、张浩还将实时 GPS/GLONASS双系统 PPP算法进行了嵌入式平台的移植[15,16]。
随着四大全球导航系统的建设基本完成,GNSS进入多模多频时代,精密产品也从传统的单 GPS或 GPS/GLONASS双系统发展为同时兼容 GPS、GLONASS、 Galileo、BDS-2/BDS-3以及部分区域增强系统。武汉大学施闯等采用全国北斗地基增强网的数据,研制了广域精密定位增强服务系统,其发布的BDS-2实时轨道钟差改正数产品的等效用户距离精度优于 2.0cm[17]。此外,针对多 GNSS并存的局面,德国地学研究中心 (Deutsches Geoforschungszentrum Potsdam,GFZ)的李星星教授、葛茂荣教授等研发了 GPS/GLONASS/BDS-2/Galileo实时精密定位服务系统[18]。进一步,Li等通过固定非差模糊度,生成了可直接用于 PPP模糊度固定的多系统整数相位钟产品,其钟差精度优于 0.1ns,相比于传统的浮点解钟差产品,各系统相位钟的精度提高了 33.9%~63.7%[19]。类似地,Geng等通过在网解过程中对载波观测值的模糊度和相位偏差进行改正,提出了一种改进的相位钟差估计方法[20]。针对大桥变形监测,Tang等估计了 1Hz的高采样率 GPS钟差改正数[21]。在 BDS-3精密产品估计方面,Li等指出,与 BDS-2的动偏 /零偏不同,对 BDS-3 MEO卫星采用连续偏航的姿态控制模型可提高其轨道的重叠精度[22]。长安大学的相关团队也发布了 BDS-3卫星的 SSR改正数产品,其实时轨道和钟差精度可达 12cm和 0.5ns[23,24]。张勤等通过对 BDS卫星天线相位中心模型精化,进一步提高了 BDS-2/BDS-3联合定轨的精度[25]。此外,欧洲定轨中心 (Center for Orbit Determination in Europe,CODE)和欧洲空间局 (European Space Agency,ESA)等分析中心也陆续提供包括 BDS-3在内的多系统轨道钟差改正数,其中 BDS-2/BDS-3的钟差解算策略通常基于二者共有的 B1I/B3I观测值。
当前实时精密轨道的精度可达厘米级,精密钟差的精度可达亚纳秒级,基本能够满足实时 PPP的精度需求[26],不过其往往存在一定的滞后性和不确定性,连续性和稳定性较差,主要体现在服务系统临时故障或网络延迟导致的数据中断或较大数据时延。尽管卫星轨道改正数具有较高的时域相关性,但卫星钟差随时间变化较快,其时域相关性较弱,对于厘米级定位而言,使用不同步卫星钟差进行定位一般只能维持少数几秒钟。当精密星历中断数十秒至数分钟时,用户终端定位的实时性和可靠性将难以保障,为解决该问题,当实时改正数中断或时延较大时,通常采用对轨道和钟差进行预报的方法。 Hauschild等分析了 GNSS卫星钟差的短期稳定性,并初步评估了不同采样率钟差的内插误差及其对 PPP的影响[27]。 Senior等和 Heo等在对卫星钟差的周期特性进行分析的基础上,在预报模型中引入了相应的周期项消除其影响,提高了模型精度[28,29]。随后,Huang等进一步在模型中同时引入周期项和自适应权函数,取得了优于 0.55ns的实时预测精度[30]。 Hadas等在对实时数据流可用性和时延进行评估的过程中,明确指出在高精度定位中不应采用已过时的改正数,同时,他们采用了多项式拟合的方法对钟差进行预报,预报 8min的精度优于 10cm[31]。El-Mowafy等采用多项式和正弦项相结合的预报模型,通过多历元的实时数据流进行模型系数的拟合,在几分钟到数小时不等的预报区间内,取得 3D优于 10cm的精度[32]。在此基础上,Nie等通过改进预报模型系数的确定方法,仅采用 1个历元的数据即可进行钟差预测,扩展了模型的可用性[33]。上述分别对轨道和钟差进行预测的方法一定程度解决了数据时延与短时中断的问题,不过考虑到轨道误差在信号传播方向的投影与钟差相关性较强,因此可进一步将二者合并为一项综合误差进行提取与播发,并且在实际应用中也可减小数据播发量[34,35]。
1.2.2多模多频 PPP函数模型
在双频 PPP数据处理中,由于无电离层组合可通过双频观测值消除电离层一阶项延迟,同时计算量较小,是应用和研究昀多的函数模型。张小红等基于自主研发的 TriP软件,通过亚太区域多模 GNSS试验(Multi-GNSS EXperiment,MGEX)站点的实测 GPS/BDS-2数据对双频 PPP浮点解收敛性能和定位精度进行了评估,结果表明单 BDS-2通常需要 80~100min收敛才能达到厘米级的精度,不论静态还是仿动态,与单 BDS-2相比,双系统组合可大幅缩短收敛时间,减小定位结果的波动[36]。丁赫等和刘金健在对 BDS-2/GPS/GLONASS组合 PPP的研究中也得到类似的结论,同时指出三系统组合动态 PPP定位性能昀好,平均收敛时间为 37min[37,38]。基于武汉大学的广域轨道钟差改正数,韩啸等在嵌入式平台对 BDS-2/GPS/GLONASS三系统组合 PPP在静态、仿动态和车载动态等不同条件下的定位性能进行评估,均取得了水平和高程厘米级的定位精度[39,40]。Li等在 GPS/BDS-2/Galileo/GLONASS组合 PPP研究中指出,多系统融合可显著改善恶劣环境下 GNSS的定位性能,即使在高度截止角为 40°的情形下,依然能够取得较稳定的定位解,此外,多系统融合一般可将收敛时间由单系统的近 20min缩短至约 10min[18]。除了双频无电离层模型外,2010年,中国科学院测量与地球物理研究所张宝成博士通过对电离层和差分码偏差 (Differential Code Bias,DCB)进行重参化,提出了基于原始观测值的非组合 PPP模型,由此可得到高精度的电离层观测值,用于全球电离层建模[41,42]。张小红