一种提高对流层天顶延迟估计精度的方法技术

本发明专利技术提供了一种提高对流层天顶延迟估计精度的方法，属于导航定位技术领域。其技术方案为：一种提高对流层天顶延迟估计精度的方法，包括针对GNSS实测信号计算结果作为对流层天顶延迟的参考值，记为ZTD

【技术实现步骤摘要】
一种提高对流层天顶延迟估计精度的方法


[0001]本专利技术涉及导航定位
，尤其涉及一种提高对流层天顶延迟估计精度的方法。

技术介绍

[0002]对流层天顶延迟是影响GNSS高精度定位的关键因素之一，忽视其2m多的误差会直接影响定位结果。此外，甚长基线干涉测量、星载多普勒定轨定位系统以及合成孔径雷达干涉等技术也受到了对流层延迟的影响，目前削弱对流层延迟误差应用最广泛的方法就是建立高精度的对流层天顶延迟模型。
[0003]通常根据模型是否需要实时气象参数，对流层天顶延迟模型分为基于实测气象参数的对流层天顶延迟模型和无需实测气象参数的对流层天顶延迟经验模型。对于基于实测气象参数的对流层天顶延迟模型，温度，压力和水汽压等都是使用过程中不可或缺的重要参数。
[0004]然而实际工作中，很多时候都无法实时获取气象参数，进而导致基于实测的模型无法广泛应用；当前已经建立无需实测气象参数的大量区域经验模型，且经过验证发现其精度并不低于基于实测气象参数的模型。但经验模型的精度也并非不具有提升空间，如何基于经验模型提高对流层天顶延迟估计精度，已经是实现快速定位、提高定位效率的关键所在。

技术实现思路

[0005]针对上述现有技术中的问题，本专利技术的目的在于提供一种提高对流层天顶延迟估计精度的方法，在GPT2w模型的基础上进行改进，建立了高精度的单站点对流层天顶延迟经验模型，通过联合GNSS实测信号解算的毫米级对流层天顶延迟数据，确保经验改进模型的精度和可靠性。
[0006]本专利技术是通过如下技术方案实现的：一种提高对流层天顶延迟估计精度的方法，包括如下步骤：
[0007]S1、针对GNSS实测信号进行对流层天顶延迟解算，将计算结果作为对流层天顶延迟的参考值，记为ZTD
GNSS
；
[0008]S2、GPT2w模型的对流层天顶延迟计算，GPT2w模型提供1
°
和5
°
两种分辨率的网格化参数信息文件，具体如下：
[0009]S21、根据GPT2w模型提供的网格化参数信息文件计算站点位置气象参数；
[0010]S22、结合Saastamoinen模型和Askne and Nordius模型并代入所述站点位置气象参数，分别计算天顶干延迟ZHD和天顶湿延迟ZWD；
[0011]S23、计算GPT2w模型对流层天顶延迟ZTD
GPT2w
，ZTD
GPT2w
＝ZHD+ZWD；
[0012]S3、基于模型误差补偿的思想，根据已有的站点数据，计算对流层天顶延迟残差dZTD：
[0013]dZTD＝ZTD
GNSS
‑
ZTD
GPT2w
；
[0014]S4、分析对流层天顶延迟残差序列中的周期信号后，构建考虑年、半年和季节周期项的残差周期模型dZTD
model
，利用最小二乘法对模型中的未知参数进行估计；
[0015]S5、构建对流层天顶延迟经验改进模型ZTD
model
：ZTD
model
＝ZTD
GPT2w
+dZTD
model
。
[0016]进一步，所述dZTD
model
具体为：
[0017][0018]其中：a0为对流层天顶延迟残差dZTD的平均值；a1～a4和分别代表年、半年以及120天和90天的振幅系数及初始相位。
[0019]进一步，所述S1具体为：通过高精度解算软件进行对流层延迟的解算，对于GNSS观测的实际接收信号，以仿实时/事后PPP处理，即使用IGS记录的轨道和时钟产品，估算得到的对流层天顶延迟作为对流层延迟的参考值ZTD
GNSS
，具体处理策略可根据实际观测情况。
[0020]进一步，所述S21具体为：
[0021]S211、首先确定站点位置周围的四个格网点，根据所述网格化参数信息文件获取格网点高度的气象参数；
[0022]S212、通过计算将所述格网点高度的气象参数垂直修正至站点高度；
[0023]S213，利用双线性插值法在水平方向将四个所述网点高度的气象参数插值到站点位置；
[0024]S214，经过计算后得到所述站点位置的气象参数，包括P、e、T
m
、λ分别表示气压、水汽压、加权平均温度与水汽递减因子。
[0025]进一步，所述S212的具体计算公式为：
[0026]T
V
＝T0*(1+0.6077Q)
[0027]e0＝Q*P0/(0.622+0.378*Q)
[0028][0029]T＝T0+dT*(h
‑
h0)
[0030]e＝e0(P*100/P0)
λ+1
[0031]其中，P0、e0、T0分别表示格网点高度的气压、水汽压和温度，Q为比湿度；P,T和e分别表示站点高度的气压、温度和水汽压；dT表示温度递减率；g
m
表示重力，取值为9.80665m/s2；T
v
表示虚温；dM和R
g
分别表示干空气的摩尔质量，值分别为28.965
×
10
‑3kg/mol和8.3143J/K/mol。
[0032]进一步，基于Saastamoinen模型计算所述天顶干延迟ZHD，公式如下：
[0033][0034]其中：H分别表示站点纬度和大地高；P表示地表气压。
[0035]进一步，根据Askne and Nordius模型计算所述天顶湿延迟ZWD，公式如下：
[0036][0037]其中：k'2、k3、R
d
均为模型公式中的常数。
[0038]一种高精度对流层天顶延迟测量系统，包括基于权利要求1～7任意一个所述的提高对流层天顶延迟估计精度的方法所构建的对流层天顶延迟经验改进模型ZTD
model
。
[0039]本专利技术的有益效果为：本专利技术构建了高精度的单站点对流层天顶延迟经验模型，通过联合GNSS实测信号解算的毫米级对流层天顶延迟数据，确保经验改进模型的精度和可靠性；相比GPT2w模型有效提高了计算精度，对提高GNSS定位精度尤其是高程方向的定位精度具有重要的实际工程应用价值；可以提高定位效率，有效降低定位过程中估计对流层延迟带来的收敛时间。
附图说明
[0040]图1为本专利技术方法步骤流程图。
[0041]图2为本专利技术几种方法的ZTD估计值对比。
[0042]图3为本专利技术中残差周期模型拟合效果。
具体实施方式
[0043]为能清楚说明本方案的技术特点，下面通过具体实施方式，对本方案进行阐述。
[0044]实施例一，参见图1，本专利技术是通过如下技术方案实现的：一种提高对流层天顶延迟估计精度的方法，包括如下步骤：
...

【技术保护点】

【技术特征摘要】
1.一种提高对流层天顶延迟估计精度的方法，其特征在于，包括如下步骤：S1、针对GNSS实测信号进行对流层天顶延迟解算，将计算结果作为对流层天顶延迟的参考值，记为ZTD
GNSS
；S2、GPT2w模型的对流层天顶延迟计算，具体如下：S21、根据GPT2w模型提供的网格化参数信息文件计算站点位置气象参数；S22、结合Saastamoinen模型和Askne and Nordius模型，并代入所述站点位置气象参数，分别计算天顶干延迟ZHD和天顶湿延迟ZWD；S23、计算GPT2w模型对流层天顶延迟ZTD
GPT2w
，ZTD
GPT2w
＝ZHD+ZWD；S3、根据已有的站点数据，计算对流层天顶延迟残差dZTD：dZTD＝ZTD
GNSS
‑
ZTD
GPT2w
；S4、分析对流层天顶延迟残差序列中的周期信号后，构建考虑年、半年和季节周期项的残差周期模型dZTD
model
；S5、构建对流层天顶延迟经验改进模型ZTD
model
：ZTD
model
＝ZTD
GPT2w
+dZTD
model
。2.根据权利要求1所述的提高对流层天顶延迟估计精度的方法，其特征在于，所述dZTD
model
具体为：其中：a0为对流层天顶延迟残差dZTD的平均值；a1～a4和分别代表年、半年以及120天和90天的振幅系数及初始相位。3.根据权利要求1所述的提高对流层天顶延迟估计精度的方法，其特征在于，所述S1具体为：通过高精度解算软件进行对流层延迟的解算，对于GNSS观测的实际接收信号，以仿实时/事后PPP处理，估算得到的对流层天顶延迟作为对流层延迟的参考值ZTD
GNSS
。4.根据权利要求1所述的提高对流层天顶延迟估计精度的方法，其特征在于，所述S21具体为：S21...

【专利技术属性】
技术研发人员：齐静静，刘佳林，李家宁，刘福兴，徐兴雨，侯增鹏，张乐，史小东，孙晨，王婷婷，徐涛，张凤，
申请(专利权)人：中国石油化工股份有限公司胜利油田分公司技术检测中心胜利油田检测评价研究有限公司，
类型：发明
国别省市：