[发明专利]一种超短期地球自转参数预报方法

权利要求书

1.一种超短期地球自转参数预报方法，其特征在于，该方法包括以下步骤：

S1、通过文件或者网络获取IERS提供的UT1-UTC基础数据和LOD基础数据，同时获取GNSS技术快速低延迟观测导出的最新的LOD基础数据；

S2、移除UT1-UTC基础数据的跳秒和潮汐，获得UT1R-TAI基础数据，并移除LOD基础数据的潮汐，获得LODR基础数据，UT1R-TAI基础数据和LODR基础数据的时间序列均分为周期项、线性项和残差项；

S3、计算UT1R-TAI基础数据在时间t+1的预报值具体包括以下步骤：

S31、计算UT1R-TAI基础数据的一阶差分序列：

式中，Δf_t^ut1为UT1R-TAI基础数据在时间t的一阶差分序列，f_t^ut1为时间t的UT1R-TAI基础数据，为时间t-1的UT1R-TAI基础数据；

S32、利用最小二乘模型拟合并预报Δf_t^ut1的周期项和线性项、LODR的周期项和线性项，同时，计算Δf_t^ut1的残差项dΔf_t^ut1和LODR的残差项dLODRt；

S33、利用多元自回归模型拟合Δf_t^ut1的残差项dΔf_t^ut1和LODR的残差项dLODR_t，并预报Δf_t^ut1的残差项；

S34、综合预报Δf_t^ut1的周期项、线性项和残差项，计算Δf_t^ut1在t+1时刻的预报值

S35、利用差分恢复算子计算UT1R-TAI基础数据在时间t+1的预报值

S4、增加的跳秒和潮汐，获取UT1-UTC基础数据在时间t+1的预报值。

2.根据权利要求1所述的一种超短期地球自转参数预报方法，其特征在于：步骤S1中，最新的LOD基础数据通过以下方法获取：下载延迟3～9小时的国际GNSS服务的LOD基础数据或者下载延迟3～9小时的国内iGMAS的LOD基础数据或者利用GNSS技术解算导出的最新超快速产品的LOD基础数据。

3.根据权利要求1所述的一种超短期地球自转参数预报方法，其特征在于：步骤S1中，最新的LOD基础数据根据GNSS原始观测数据进行超快速LOD参数解算获取：

a、实时获取GNSS卫星最新的24小时观测值文件、导航文件、先验的地球自转参数文件和GNSS地球自转参数精密确定的通用配置文件；

b、根据GNSS观测数据建立GNSS观测方程：

式中，为理论距离，ρ为观测距离，Δρ为观测误差，X_s为卫星在地固坐标系下的坐标，X_r为测站在地固坐标系下的坐标，P、N、S、U分别为岁差、章动、自转和极移从地固坐标系到惯性坐标系的旋转矩阵；

GNSS观测方程线性化为：

式中，LOD₀为先验LOD参数，θ为格林泥治真恒星时，γ≈1.0027379093，t为观测时刻的儒略日，t₀为所选参考时刻的儒略日；

c、线性化后的GNSS观测方程，根据最小二乘法解算先验LOD参数LOD₀的改正数，最终获取GNSS技术导出的最新的快速低延迟的LOD基础数据。

4.根据权利要求1所述的一种超短期地球自转参数预报方法，其特征在于：步骤S2中，当预报时长大于2天时，采用最新的LOD基础数据。