PPP在对流层延迟估计和固体潮测定中的应用研究.docVIP

PPP在对流层延迟估计和固体潮测定中的应用研究 announced the zenith of troposphere delay compare, precise point positioning estimate the reliability of the troposphere delay has quite, but also for further research of the tidal correction model have certain reference significance. 　　Keywords: precise point positioning; Without the ionosphere combination model; Troposphere delay; Tidal correction 　　中图分类号：P223文献标识码：A 文章编号： 　　 GPS气象学的重要内容是利用天顶对流层延迟来反演大气可降水量，常规的估计对流层延迟的方法是双差网解法，但该方法需要引入远距离测站来估计站点间的对流层绝对时延；[1，2]目前固体潮的测定主要是通过重力仪、倾斜仪和伸缩仪来观测，观测的点数和精度受到一定的限制。随着精密单点模型的逐步改进和定位精度的不断提高，特别是实时精密单点定位更有广阔的应用空间。精密单点定位技术在对流层延迟估计中有很多的优势：估计模型简单，站与站之间不相关，无需引入远距离测站即可估计绝对时延，处理大规模数据速度快等。[3]而且采用空间技术测定固体潮在精度和范围以及时间上(全天候连续观测)具有潜在的优势。 　　 本文基于PPP技术，编写了精密单点定位程序，通过大量的数据计算验证了该程序在估计测站天顶对流层延迟和计算固体潮对点位的影响方面的可靠性。 　　精密单点定位基本原理 　　 GPS非差精密单点定位一般是以单台双频GPS接收机采集数据，利用IGS提供的精密星历和卫星钟差，通过数据预处理和各种误差改正，采用非差观测模型进行的单点定位。其中数据预处理包括：粗差剔除、周跳的探测和修复、初始整周模糊度的确定、相位平滑伪距、精密星历和钟差的内插等。精密单点定位的传统消电离层组合模型，是采用双频GPS伪距和载波相位观测值的无电离层组合来构成观测模型，这种无电离层组合是减弱电离层影响最有名的公式，其观测模型的简化形式如下：[4] 　　（1.1） 　　 （1.2） 　　式中，为、无电离层伪距组合观测值；为、无电离层相位组合观测值（距离）；（）为载波频率；为卫星到接收机的几何距离；为光速，为接收机钟差，为对流层延迟改正参数；、分别为两种组合观测值的观测噪声和未被模型化的误差[5]；为无电离层折射组合相位观测值（观测）的模糊度。此外，观测模型还要做的误差改正有相对论效应、固体潮、海洋潮汐、极移、章动等，限于篇幅，不作详细论述。 　　精密单点定位参数估计方法 　　 将公式(1.1)和(1.2)分别在近似值处线性化，得到线性化的伪距和载波相位观测方程为： 　　 （2.1） 　　 （2.2） 　　 当观测到颗卫星，此时待估计参数为即： 　　 （2.3） 　　 待估参数： 为接收机位置坐标改正参数；为接收机钟差改正参数；为对流层延迟改正参数；为消电离层组合相位观测值的整周未知数相对应的距离改正参数。 　　 利用上述推导的线性化模型，再结合采用测码伪距和载波相位伪距观测值，就可以按照最小二乘法进行非差精密单点定位的平差计算。在解算过程中，静态情况下，坐标参数可以作为常数未知数处理；未发生周跳或修复周跳的情况下，整周模糊度当作常数处理，在发生周跳的情况，整周模糊度作为一个新的常数参数进行处理；因为接收机的钟不稳定，有明显的随机抖动。 　　3.精密单点定位在对流程延迟估计中的应用 　　 本文在采用精密单点定位技术估计对流层延迟参数时，采用消电离层影响的P码和载波相位组合观测值作为数据处理的观测值，并通过高度角定权的方式对观测值进行定权。在数据预处理得到“干净”的观测值的基础上，模糊度估计作为常数处理，发生周跳的地方则重置参数，采用第二节中的平差策略对参数进行估计。以下是计算上海佘山跟踪站观测数据，采用事后精密星历和精密钟差单天解估计天顶对流层延迟参数的结果（原始数据都是从IGS网站下载得到的）： 　　 使用上海佘山IGS跟踪站2011年5月4日30s间隔采样率的观测数据估计天顶对流层延迟随时间变化的情况： 　　 图1.1 2011年5月4日上海佘山跟踪站天顶对流层延迟变化情况 　　 使用上海佘山IGS跟踪站2011年7月6日30s间隔采样率的观测数据估计天顶对流层延迟随时间变化情况： 　　图1.2 2011年7月6日上海佘山跟踪站天顶对流层延迟变化情况