论文部分内容阅读
全球导航卫星系统(Global Navigation Satellite System, GNSS)是二十世纪人类最伟大的科技成果之一,且已经在大地测量、地球动力学以及大气研究等多个领域得到广泛应用。近年来,随着我国北斗导航定位系统(BDS)二代的全面运行、欧洲伽利略系统(GALILEO)的不断发展以及美国GPS系统的现代化,GNSS已经从单系统发展为多系统并存、双频服务发展为至少三个频率服务的局面,这给目前的GNSS数据处理方法带来了极大的挑战。一方面,目前所有公开的GNSS产品均是基于无电离层组合观测值或差分观测值的数据处理方法获得,然而这种方法己经无法完全满足当前多系统多频率环境下的高精度GNSS数据处理要求了。首先为了充分利用各个系统之间的优势,得到更稳定可靠的结果,并为其他用户提供兼容统一的改正数信息,必须将多系统多频率的数据在观测值层面进行统一处理;其次为了尽可能保留原始观测值中的信息,必须在原始观测值而不是差分或组合观测值上建立观测方程。针对此,本文详细研究了一种基于原始观测值的数据处理方法,任意GNSS系统任意频率的数据均能统一采用这种方法进行数据处理。另一方面,计算效率低一直是GNSS领域中未能很好解决的问题。即使仅针对GPS系统,采用组合或差分观测值进行数据处理的方法,几乎目前所有著名GNSS软件都难以或无法同时对含有数百甚至上千个测站的GNSS网进行整网数据处理。然而相对当前组合观测值或差分观测值数据处理方法,原始观测值数据处理方法耗时更为严重,即使针对含有约100个测站的GNSS网,也难以进行日常常规处理。虽然先分子网解算,再合并子网结果的策略可在一定程度缓解这种矛盾,但由于子网问的公共站多次被使用,影响了整网解的严密性。为此,本文提出了一整套能加速GNSS数据处理的策略,不仅可使基于无电离层组合观测值的数据方法可以在短时间内处理含有数百甚至上千测站的GNSS网,还能有效加速原始观测值数据处理。本文的主要工作和贡献为:(1)归纳总结了当前GNSS数据处理方法面临的挑战;系统地研究了基于原始观测值的多系统多频率GNSS数据处理模型;给出了分离钟差参数、系统偏差/频率偏差参数、DCB参数等多类参数的方法;讨论了电离层参数、DCB参数,模糊度之间的相关性;并探讨了在不施加电离层时空约束条件下的模糊度固定方法。(2)利用原始观测值数据处理方法初步探讨了GPS和BDS频率之间的兼容性。实验发现:GPS L5的残差约为L1/L2残差的3倍,可达9~10mm,而且连续数天内其变化量级和变化趋势基本一致,这意味着GPS L5可能含有某种与L1/L2不兼容的系统偏差。而BDS三个频率的残差均在2-3mm,没有发现类似的情况。实验还发现L5码伪距上的DCB估值在数天甚至一个月内都十分稳定,月平均值的标准差优于3cm。(3)评估了在使用两个频率且不施加电离层约束的条件下,基于原始观测值的GNSS数据处理方法与基于无电离层组合观测值的数据处理方法得到的产品的异同。结果发现两种方法计算得到的轨道、钟差、对流层以及系统间偏差均无明显差异,这表明当使用两个频率且不施加电离层时空约束时,基于原始观测值的数据方法与无电离层组合观测值的数据方法基本等价。(4)提出了利用PPP模糊度技术在统一基准下生成Carrier-range的新方法,并在此基础上设计了一整套高效率GNSS数据处理流程。本文从PPP定位的基本观测方程出发,证明了利用单站模糊度技术在统一基准下生成Carrier-range的可行性,并设计了整套高效率GNSS数据流程。主要包括六个步骤:采用全球均匀分布的适合的网计算精密轨道和卫星钟差;利用定轨过程中得到的浮点模糊度计算卫星端的UPD;固定卫星轨道和钟差,逐一测站进行PPP定位,并利用生成的UPD进行单站模糊度固定;然后将载波相位观测值转换成Carrier-range,并将其写入新的RINEX中;重复上两步直至所有测站均生成了新的RINEX文件;利用含有Carrier-range的新RINEX文件进行最终的网解。最终网解不再需要估计或者仅需要估计少量模糊度,其计算效率将大幅度提高。最终解算可采取两种模式,一种仅采用Carrier-range进行计算,此时UPD参数将完全被卫星钟差吸收,不再需要估计,卫星钟差也由于吸收了UPD,成了“整数钟”。另一种方法是同时采用伪距和Carrier-range,此时为了使伪距和Carrier-range钟差定义一致,在Carrier-range上需要估计一组UPD参数,此时钟差绝对基准由伪距决定,与IGS定义保持一致。(5)统计分析了UPD参数的稳定性,提出了UPD参数估计方法。通过对实验数据的分析,发现非经历地影的BLOCK IIA GPS卫星的窄巷UPD在24h内变化平稳,其一天平均值的标准差大都不大于0.05周,而经历地影的BLOCK IIA卫星在出地影后,窄巷UPD可能发生0.3~0.4周的跳变。若去掉这种跳变,其天均值的标准差降为0.03~0.04左右。因此本文建议在求解最终估计UPD的时候,至少需要对经历地影的BLOCK IIA卫星在出地影后额外估计一个UPD参数。(6)发现并指出Carrier-range可改善数据的连续性,进而可提高结果质量。通过比较分析和实例验证,发现相对于传统方法,Carrier-range策略具有更好的数据连续性,其等价于将同一测站同一卫星的所有模糊度均连起来了,从而拥有更好的数据连续性。实验发现采用同样的数据,采用Carrier-range策略得到的轨道的天与天之间的重复轨道RMS值相对传统方法提高约3~4mm。(7)提出了加快无电离层组合观测值数据处理策略,解决了大规模GNSS网数据处理效率低的难题。实验发现:当解算含有460个测站的网时,新策略耗时仅14分钟,而传统策略耗时约82分钟,而且新策略所需的计算时间随着测站数的增长近似呈线性增长,而传统策略耗时随着测站数增长近似呈指数增长。(8)设计了事后和实时高频整数钟差的高效率计算流程,有效提高了事后和实时高频钟差或“整数钟差”计算的效率。实验证明,当采用约250个测站事后计算一天内30s钟差时,耗时仅约32分钟,而当采用约200个测站实时计算钟差时,在所有模糊度固定了的情况下,一个历元耗时不到1s,且直接利用Carrier-range计算的得到的钟差包含了UPD,可直接用于单站模糊度固定。(9)提出了基于原始观测值的数据高效处理方法,有效提高了计算效率,使得采用该方法进行数据日常常规处理成为可能。实验显示,当处理含约100个测站的GNSS网时,在不施加电离层时空约束的情况下,新策略可使单次参数估计时间从97分钟缩短为31分钟,在施加一简单的电离层时空约束后,新策略可使单次参数估计时间从387分钟减少为206分钟。这表明新策略后可有效加速原始观测值的数据处理。