一种基于GRNN模型的区域NWP对流层延迟改正方法技术

本发明专利技术公开了一种基于GRNN模型的区域NWP对流层延迟改正方法，首先获取作业区域里选取所有连续运行参考站的连续一年的ZTD和该区域NWP提供的连续一年的再分析资料，利用积分法反演该区域上述连续运行参考站的ZTD，然后选取其中一部分站点为训练站点，其余部分站点为测试站点，利用GRNN拟合这些训练站点的NWP资料反演ZTD的残差，得到GRNN残差拟合模型，然后通过GRNN残差拟合模型对测试站点的NWP资料反演的ZTD做误差补偿，即得到测试站点精确的NWP_ZTD。本方法首次根据NWP_ZTD的残差随多种气象因素的变化规律，通过GRNN模型拟合NWP_ZTD的残差对NWP_ZTD做误差补偿，提高NWP反演对流层延迟精度。

A Regional NWP Tropospheric Delay Correction Method Based on GRNN Model

【技术实现步骤摘要】
一种基于GRNN模型的区域NWP对流层延迟改正方法
本专利技术涉及一种对流层延迟改正方法，尤其是一种基于广义回归神经网络(GeneralRegressionNeuralNetwork,简称GRNN)模型的区域数值天气预报(NumericalWeatherPrediction,简称NWP)模型对流层延迟改正方法，根据NWP_ZTD的残差随多种气象因素的变化规律，通过GRNN模型来拟合NWP_ZTD的残差，然后对NWP_ZTD做误差补偿，提高NWP反演对流层延迟精度，从而提高精密单点定位(PrecisePointPositioning,简称PPP)和长距离基线实时精密动态定位(RealTimeKinematic,简称RTK)的收敛速度和定位精度，属于全球卫星导航与定位

技术介绍
对流层延迟误差一直是影响高精度全球卫星导航系统(GlobalNavigationSatelliteSystem,简称GNSS)测量的主要因素，在无线电信号传播过程中，对流层从天顶方向到地平方向对无线电信号的延迟可达2m到20m，严重制约着GNSS精密单点定位和中长距离基线差分定位的模糊度收敛速度和定位精度。在GNSS数据处理的现有技术中，总结起来有三种方法可用来削弱对流层延迟：1.通过经验模型估计天顶对流层延迟(ZenithTroposphericDelay,简称ZTD)，但经验模型估计的ZTD的精度差，最大误差可达20cm；2.在数学模型中将ZTD作为参数与其他未知参数(如站点坐标和模糊度)共同求解。然而，对流层延迟与位置强相关，方程需要较强的卫星几何约束和较长时间的观测才能得到稳定解；3.利用数值天气预报模型的再分析或预报资料反演ZTD。现如今NWP预报资料反演ZTD的残差大概在±6cm上下浮动，平均残差和均方根偏差(RootMeanSquareDeviation，简称RMSD)分别在1cm和3cm左右，不能直接用于GNSS精密定位中。在GNSS精密定位解算中，通常将此ZTD作为初值，将ZTD的残差作为未知数，与模糊度、位置等参数共同求解。此时，NWP反演ZTD的精度将直接影响模糊度的收敛速度。要实现实时PPP或长距离RTK，需要提高NWP反演ZTD的精度，即改正NWP反演的ZTD的残差，针对这一问题，目前还未有切实可行的方法。
技术实现思路
本专利技术的目的在于，为提高卫星定位精度，对其中一项重要的误差源对流层延迟进行改正，提供了一种基于GRNN模型的区域NWP对流层延迟改正方法，用于解决现有的对流层延迟估计方法精度不足的技术问题。为了解决
技术介绍
所存在的问题，本专利技术提出了一种基于GRNN模型的区域NWP对流层延迟改正方法，包括以下步骤：第一步，参考数据获取。获取作业区域里所有连续运行参考站(ContinuousOperationalReferenceSystem，简称CORS)提供的连续一年的高精度ZTD产品，记做IGS_ZTD，并剔除粗差及数据缺失严重的站点，该数据可在各省市的CORS网中心网站获得；第二步，NWP数据获取。获取NWP模型提供的该区域与第一步同年的连续一年的再分析资料，其中包括气压、温度、相对湿度及位势高度，该数据可从各数值预报资料中心获得；第三步，NWP_ZTD计算。利用NWP模型的气象资料通过分层积分法反演该区域上述CORS站的ZTD，记做NWP_ZTD。采用分层积分法反演ZTD的计算公式为：上式中，ZTDgrid代表测站周围格网点(有4个)在测站所在高度上的ZTD值，H_IGS代表测站的高程，H_top代表NWP资料顶层高度，N代表大气折射率(无量纲)，Nn代表第n个积分区域范围内的大气折射率，△Hn代表第n个积分区域的高度。其中大气折射率N的表达公式如下：上式中，k1＝77.604K/mbar，k2＝64.79K/mbar，k3＝377600.0K2/mbar，P代表压强，e代表水汽压，单位均为mbar，Ntry代表干性气体的折射率，Nwet代表湿性气体的折射率。计算出测站周围4个格网点在测站所在高度上的ZTD后，通过双线性内插法内插出测站点的ZTD。第四步，残差计算。以IGS_ZTD作为参考，求出NWP_ZTD的残差；第五步，模型建立。选取作业区域内均匀分布的部分站点为训练站点，另一部分站点为测试站点，从NWP资料获取训练站点的温度和相对湿度作为模型输入，训练站点的NWP_ZTD的残差为模型输出，采用GRNN模型拟合残差，得到残差拟合模型；第六步，残差拟合。从NWP资料获取测试站点的温度和相对湿度输入上述GRNN残差拟合模型，拟合出测试站点的NWP_ZTD的残差；第七步，精确NWP_ZTD获得。测试站点的NWP_ZTD减去GRNN模型拟合出的测试站点的NWP_ZTD的残差，即可得到测试站点精确的NWP_ZTD。与现有技术相比，本专利技术的有益效果是：第一，首次根据测站点的温度和相对湿度通过GRNN模型来拟合NWP_ZTD的残差，然后对NWP_ZTD做误差补偿，提高NWP反演对流层延迟精度。第二，本专利技术提出的基于GRNN模型的区域NWP对流层延迟改正方法可以将网络RTK中的参考站作为训练数据进行建模，然后精确估计用户站与因信号遮挡或中断而导致高精度对流层产品缺失的参考站的对流层延迟，从而提高网络RTK的性能。第三，由于NWP模型可提供覆盖全球的气象数据，所以本方法可估计全球NWP模型覆盖下的任意一点的ZTD，适用范围很广。附图说明图1是本专利技术的工作流程图。图2是100个训练站点分布图。图3是550个测试站点分布图。图4是测试站点采用本专利技术方法前后的ZTD残差对比图。图5是测试站点采用本专利技术方法前后的ZTD残差的RMSD对比图。说明：图4中(a)和(b)分别为550个测试站点2005年采用本专利技术方法改正前后的2005年ZTD残差，这里的平均残差计算方法是测站一年的残差取绝对值再平均；图5中(a)和(b)分别为550个测试站点采用本专利技术方法改正前后的ZTD残差的RMSD。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。实施例为证明基于GRNN模型的区域NWP对流层延迟改正方法的有效性，下面选取日本地区650个测站2005年连续一年的采样率为2小时的NCAR对流层数据，以及测站区域相应的欧洲中尺度天气预报中心(EuropeanCentreforMedium-RangeWeatherForecasts,简称ECMWF)再分析资料中ERA-Interim产品的分层气象数据，其平面分辨率为0.125°×0.125°，垂直分辨率为37层(顶层高度约为47km)，时间分辨率为6小时。日本区域的矩形面积约300万平方公里，实验区域范围为32°N～40°N，130°E～142°E。从650个测站中选取100个测站的数据为GRNN训练数据，其余550个测站的数据为GRNN测试数据。训练测站分布如图2所示，测试测站分布如图3所示。第一步，参考数据获取。获取日本地区近1000个测站2005年连续一本文档来自技高网...

【技术保护点】
1.一种基于GRNN模型的区域NWP对流层延迟改正方法，其特征在于，包含以下具体步骤：第一步，参考数据获取。获取作业区域里所有CORS站提供的连续一年的高精度ZTD产品，记做IGS_ZTD，并剔除粗差及数据缺失严重的站点，该数据可在各省市的CORS网中心网站获得；第二步，NWP数据获取。获取NWP模型提供的该区域与第一步同年的连续一年的再分析资料，其中包括气压、温度、相对湿度及位势高度，该数据可从各数值预报资料中心获得；第三步，NWP_ZTD计算。利用NWP模型的气象资料通过分层积分法反演该区域上述CORS站的ZTD，记做NWP_ZTD。采用分层积分法反演ZTD的计算公式为：

【技术特征摘要】
1.一种基于GRNN模型的区域NWP对流层延迟改正方法，其特征在于，包含以下具体步骤：第一步，参考数据获取。获取作业区域里所有CORS站提供的连续一年的高精度ZTD产品，记做IGS_ZTD，并剔除粗差及数据缺失严重的站点，该数据可在各省市的CORS网中心网站获得；第二步，NWP数据获取。获取NWP模型提供的该区域与第一步同年的连续一年的再分析资料，其中包括气压、温度、相对湿度及位势高度，该数据可从各数值预报资料中心获得；第三步，NWP_ZTD计算。利用NWP模型的气象资料通过分层积分法反演该区域上述CORS站的ZTD，记做NWP_ZTD。采用分层积分法反演ZTD的计算公式为：上式中，ZTDgrid代表测站周围格网点在测站所在高度上的ZTD值，H_IGS代表C测站的高程，H_top代表NWP资料顶层高度，N代表大气折射率，Nn代表第n个积分区域范围内的大气折射率，△Hn代表第n个积分区域的高度；其中大气折射率N的表达公式如下：上式中，k1＝77.604K/mbar，k2＝64.79K/mbar，k3＝377...

【专利技术属性】
技术研发人员：徐莹，闫俐孜，李雷，刘凡，刘国林，
申请(专利权)人：山东科技大学，
类型：发明
国别省市：山东,37