Search

CN-121541223-B - Wide-area troposphere modeling method integrating navigation satellite and meteorological model

CN121541223BCN 121541223 BCN121541223 BCN 121541223BCN-121541223-B

Abstract

The invention discloses a wide-area troposphere modeling method integrating a navigation satellite and a meteorological model, which relates to the technical field of satellite positioning, and comprises the steps of acquiring first troposphere zenith wet delay through meteorological data provided by a global numerical meteorological prediction model in a global navigation satellite ground observation station distribution missing area, interpolating through an Akima interpolation algorithm to generate second troposphere zenith wet delay synchronous with the global navigation satellite ground observation station, acquiring corresponding third troposphere zenith wet delay through the second troposphere zenith wet delay after constructing a virtual reference station meeting requirements according to a preset principle, integrating fourth troposphere zenith wet delay and third troposphere zenith wet delay estimated by the global navigation satellite ground observation station, utilizing the integrated troposphere zenith wet delay and an improved troposphere polynomial fitting model, and determining a target coefficient of the improved troposphere polynomial fitting model suitable for satellite end broadcasting for use.

Inventors

  • Cui Bobin
  • WANG LE
  • HUANG GUANWEN
  • DU SHI

Assignees

  • 长安大学

Dates

Publication Date
20260512
Application Date
20260120

Claims (4)

  1. 1. A wide-area troposphere modeling method integrating navigation satellites and meteorological models is applied to a server and is characterized by comprising the following steps: Acquiring first tropospheric zenith humidity delay through meteorological data provided by a global numerical meteorological prediction model in a distributed missing area of a global navigation satellite ground observation station, and interpolating the first tropospheric zenith humidity delay through an Akima interpolation algorithm to generate second tropospheric zenith humidity delay synchronous with the global navigation satellite ground observation station; Constructing a virtual reference station meeting the requirements according to a preset principle, and acquiring a third troposphere zenith humidity delay corresponding to the virtual reference station through the second troposphere zenith humidity delay; The fourth troposphere zenith humidity delay estimated by the global navigation satellite ground observation station and the third troposphere zenith humidity delay are fused, and the fused troposphere zenith humidity delay is obtained; modeling by utilizing the fused troposphere zenith wet delay and the improved troposphere polynomial fitting model, and determining a target coefficient of the improved troposphere polynomial fitting model suitable for satellite end broadcasting; broadcasting the target coefficient through a satellite end for the user end to use; Modeling by using the fused troposphere zenith moisture delay and the improved troposphere polynomial fitting model, determining target coefficients of the improved troposphere polynomial fitting model suitable for satellite end broadcasting, wherein the method comprises the following steps: Based on the fused troposphere zenith wet delay, carrying out parameter estimation by adopting a least square method, and removing abnormal values in the estimation process by a median detection method to obtain an initial troposphere polynomial fitting model; when the initial troposphere polynomial fitting model meets a preset condition, taking the coefficient of the initial troposphere polynomial fitting model as a target coefficient of the current period, otherwise taking the coefficient of the initial troposphere polynomial fitting model of the previous period as a target coefficient; The interpolation is performed on the first troposphere zenith humidity delay through an Akima interpolation algorithm to generate a second troposphere zenith humidity delay synchronous with a global navigation satellite ground observation station, and the method comprises the following steps: interpolation in space dimension and time dimension is carried out on the first troposphere zenith moisture delay at discrete time through an Akima interpolation algorithm, so that a second troposphere zenith moisture delay matched with the observation time of the global navigation satellite ground observation station is obtained; the preset principle comprises the steps of setting virtual reference stations according to the distribution and the topographic features of the global navigation satellite ground observation stations, wherein the distance between the virtual reference stations and the global navigation satellite ground observation stations does not exceed a fixed threshold value, and the height difference between different global navigation satellite ground observation stations is kept within a preset height.
  2. 2. The method of claim 1, wherein the improved tropospheric polynomial fitting model is represented by the following formula: ; In the above-mentioned method, the step of, For the fused troposphere zenith wet delay, As the scale factor of the water vapor, 、 、 、 、 And All of the parameters are the parameters of the change, Is the height of the object to be processed, For the latitude deviation corresponding to the global navigation satellite ground observation station, And the longitude deviation corresponding to the global navigation satellite ground observation station.
  3. 3. The method of claim 1, wherein the fourth tropospheric zenith wet delay and the third tropospheric zenith wet delay estimated by the global navigation satellite terrestrial survey station are fused, the method further comprising, prior to obtaining the fused tropospheric zenith wet delay: And estimating the zenith wet delay of the fourth troposphere of the global navigation satellite ground observation station in real time by using a non-differential and non-combined precise single-point positioning ambiguity fixing technology.
  4. 4. A wide-area troposphere modeling method integrating navigation satellites and meteorological models is applied to a user side, and is characterized by comprising the following steps: receiving a target coefficient broadcasted from a satellite terminal; Calculating the value of local zenith wet delay by using the received target coefficient according to the longitude and latitude of the position of the user terminal and the user elevation; The target coefficient is obtained by distributing a missing area in a global navigation satellite ground observation station, acquiring first troposphere zenith humidity delay through meteorological data provided by a global numerical weather forecast model, and carrying out interpolation on the first troposphere zenith humidity delay at discrete time in space dimension and time dimension through an Akima interpolation algorithm to obtain second troposphere zenith humidity delay matched with the observation time of the global navigation satellite ground observation station; Constructing a virtual reference station meeting the requirements according to a preset principle, and acquiring a third troposphere zenith humidity delay corresponding to the virtual reference station through the second troposphere zenith humidity delay; The fourth troposphere zenith humidity delay estimated by the global navigation satellite ground observation station and the third troposphere zenith humidity delay are fused, and the fused troposphere zenith humidity delay is obtained; based on the fused troposphere zenith wet delay, carrying out parameter estimation by adopting a least square method, and removing abnormal values in the estimation process by a median detection method to obtain an initial troposphere polynomial fitting model; when the initial troposphere polynomial fitting model meets a preset condition, taking the coefficient of the initial troposphere polynomial fitting model as a target coefficient of the current period, otherwise taking the coefficient of the initial troposphere polynomial fitting model of the previous period as a target coefficient; the preset principle comprises the steps of setting virtual reference stations according to the distribution and the topographic features of the global navigation satellite ground observation stations, wherein the distance between the virtual reference stations and the global navigation satellite ground observation stations does not exceed a fixed threshold value, and the height difference between different global navigation satellite ground observation stations is kept within a preset height.

Description

Wide-area troposphere modeling method integrating navigation satellite and meteorological model Technical Field The invention relates to the technical field of satellite positioning, in particular to a wide-area troposphere modeling method for fusing a navigation satellite and a meteorological model. Background Tropospheric delay accurate modeling is critical to global navigation satellite system (Global Navigation SATELLITE SYSTEM, GNSS) in high-precision applications such as geodetic positioning, navigation, and timing. Troposphere, and particularly its wet delay portion, introduces signal delays into the precise single point positioning (Precise Point Positioning, PPP) of GNSS, which, if not effectively attenuated, would significantly reduce positioning accuracy. However, the sparse spatial distribution of global navigation satellite ground observation stations limits the construction of high resolution tropospheric delay fields, especially in areas of complex terrain, near coastlines or lacking sufficient observation infrastructure, which highlights the importance of spatial resolution and accuracy of tropospheric delay estimation. In the related technology, the total zenith delay is calculated through the observation data of the global navigation satellite ground observation station, and the zenith wet delay is obtained by subtracting the zenith dry delay calculated by the meteorological parameters and the model, but the method depends on the quality of GNSS observation data and a data processing strategy, can not obtain results when the global navigation satellite ground observation station is sparse or the data is interrupted, and is difficult to independently and accurately calculate in a region without a meteorological sensor. Therefore, how to calculate the zenith wet delay of different areas more accurately and effectively becomes a problem to be solved urgently. Disclosure of Invention Therefore, the implementation of the invention provides a wide-area troposphere modeling method integrating navigation satellites and meteorological models, which solves the problem that the related technology cannot accurately and effectively calculate zenith wet delay of different areas. According to a first aspect of an embodiment of the present invention, there is provided a wide-area troposphere modeling method for fusing a navigation satellite and a meteorological model, applied to a server, including: Acquiring first tropospheric zenith humidity delay through meteorological data provided by a global numerical meteorological prediction model in a distributed missing area of a global navigation satellite ground observation station, and interpolating the first tropospheric zenith humidity delay through an Akima interpolation algorithm to generate second tropospheric zenith humidity delay synchronous with the global navigation satellite ground observation station; Constructing a virtual reference station meeting the requirements according to a preset principle, and acquiring a third troposphere zenith humidity delay corresponding to the virtual reference station through the second troposphere zenith humidity delay; The fourth troposphere zenith humidity delay estimated by the global navigation satellite ground observation station and the third troposphere zenith humidity delay are fused, and the fused troposphere zenith humidity delay is obtained; modeling by utilizing the fused troposphere zenith wet delay and the improved troposphere polynomial fitting model, and determining a target coefficient of the improved troposphere polynomial fitting model suitable for satellite end broadcasting; And broadcasting the target coefficient through the satellite end for the user end to use. Optionally, the improved tropospheric polynomial fitting model is expressed by the following formula: ; In the above-mentioned method, the step of, For the fused troposphere zenith wet delay,As the scale factor of the water vapor,、、、、AndAll of the parameters are the parameters of the change,Is the height of the object to be processed,For the latitude deviation corresponding to the global navigation satellite ground observation station,And the longitude deviation corresponding to the global navigation satellite ground observation station. Optionally, the modeling using the fused tropospheric zenith wet delay and the improved tropospheric polynomial fitting model, determining target coefficients of the improved tropospheric polynomial fitting model suitable for satellite-side broadcasting includes: Based on the fused troposphere zenith wet delay, carrying out parameter estimation by adopting a least square method, and removing abnormal values in the estimation process by a median detection method to obtain an initial troposphere polynomial fitting model; And when the initial tropospheric polynomial fitting model meets a preset condition, taking the coefficient of the initial tropospheric polynomial fitting model as a target coefficient of