CN-122017887-A - Combined correction method, system and storage medium for GNSS non-structural deformation

CN122017887ACN 122017887 ACN122017887 ACN 122017887ACN-122017887-A

Abstract

The invention discloses a GNSS non-structural deformation joint correction method, a GNSS non-structural deformation joint correction system and a storage medium. According to the method, GNSS pseudo-range and carrier phase observation data, precise products such as precise orbit and clock error are obtained, a three-component displacement time sequence and a temperature-driven thermoelastic displacement time sequence which are caused by non-tidal atmospheric load, non-tidal ocean load and land hydrologic load are constructed, and the three-component displacement time sequence and the temperature-driven thermoelastic displacement time sequence are unified into a reference frame and a time system which are consistent with GNSS calculation. The method comprises the steps of interpolating the unstructured displacement vector to an observation epoch according to time, converting the unstructured displacement vector to a geocentric rectangular coordinate system by a site local coordinate system ENU, projecting the coordinate system to a sight line direction from a satellite to a site reference point to form an equivalent geometric distance correction term, applying the correction term to pseudo-range observation and carrier phase equivalent meter domain observation respectively in the same domain with the observation, further executing precise resolving of GNSS daily arc segments to output a coordinate time sequence, and preferably being used for a precise single-point positioning whole-cycle ambiguity fixing technology.

Inventors

LU RAN
LI ZHAO
WANG JIAN
ZHANG MINGYUAN
CHEN SHAN

Assignees

武汉大学

Dates

Publication Date: 20260512
Application Date: 20260416

Claims (10)

1. A method for joint correction of GNSS non-structural deformations, comprising: acquiring GNSS observation data of a station to be processed and a precision product, wherein the GNSS observation data comprises pseudo-range observation and carrier phase observation of at least one satellite system, and the precision product comprises an orbit, a clock error and a deviation product for fixing ambiguity; The construction of the non-structural displacement prior correction comprises the steps of obtaining or calculating a three-component displacement time sequence caused by non-tidal atmospheric load, non-tidal ocean load and land hydrologic load and a thermoelastic displacement time sequence forcedly generated by temperature, summing the three-component displacement time sequence to obtain a non-tidal load displacement time sequence, unifying the non-tidal load displacement time sequence and the thermoelastic displacement time sequence to a reference frame and a time system consistent with GNSS (Global navigation satellite System) calculation, and forming a non-structural displacement vector under a site local coordinate system; Performing time interpolation on the unstructured displacement vector to obtain a displacement increment of each observation epoch, and converting the displacement increment under the local coordinate system into a displacement increment under a geocentric rectangular coordinate system; Calculating the line-of-sight unit vector from the satellite to the reference point of the measuring station for each observation epoch and each satellite, and projecting the displacement increment under the geocentric rectangular coordinate system to the line-of-sight direction to form an equivalent geometric distance correction term; applying the equivalent geometric distance correction term to pseudo-range observation and carrier phase observation in the same domain as the observation to obtain a corrected observation value; and executing GNSS arc-segment precise calculation based on the corrected observation value, and outputting a coordinate time sequence.
2. The joint correction method of GNSS non-structural deformations according to claim 1, wherein said non-structural displacement vectors are expressed as: , Wherein, the Is the sum of three-component displacement caused by non-tidal atmospheric load, non-tidal ocean load and land hydrologic load, Is a thermo-elastic displacement comprising at least a vertical component.
3. The method of claim 2, wherein when the thermoelastic displacement provides only a vertical component, the horizontal component of the non-structural displacement vector is zero and the vertical component is the sum of the non-tidal load vertical displacement and the thermoelastic vertical displacement.
4. The method for joint correction of GNSS non-structural deformation according to claim 1, wherein the time interpolation of the non-structural displacement vectors comprises calculating the non-structural displacement vectors with sampling rate of full hour by time interpolation to obtain the displacement increment of each observation epoch.
5. The joint correction method of GNSS non-structural deformation according to claim 1, wherein the local coordinate system is an east-north-day coordinate system, and the conversion of the displacement increment in the local coordinate system to the displacement increment in the geocentric rectangular coordinate system is achieved by the following formula: , Wherein, the The unit vectors of the stations in the local east-north-day coordinate system are respectively, Respectively in the observation epoch The displacement of the station in the local east direction, north direction and vertical direction.
6. The joint correction method of GNSS non-structural deformations according to claim 1, wherein said equivalent geometric distance correction term is calculated by: , Wherein, the Is the line of sight unit vector of the satellite to the station reference point, Is the displacement increment under the rectangular coordinate system of the earth center.
7. The method of claim 1, wherein applying the equivalent geometric distance correction term to the pseudorange observations and the carrier phase observations in the same domain as the observations comprises correcting the pseudorange observations directly in the rice domain, and multiplying the carrier phase observations by the wavelength converted equivalent rice domain amounts and performing the same domain correction.
8. The method of claim 1, wherein performing precise GNSS arc segment resolution includes further performing whole-cycle ambiguity fixing using a precise single-point positioning ambiguity fixing technique, and outputting a time sequence of day coordinates.
9. A joint correction system for GNSS non-structural deformations, comprising: The system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for acquiring GNSS observation data of a station to be processed and a precision product, the GNSS observation data comprise pseudo-range observation and carrier phase observation of at least one satellite system, and the precision product comprises a track, a clock error and a deviation product for fixing ambiguity; The non-structural displacement construction module is used for acquiring or calculating a three-component displacement time sequence caused by non-tidal atmospheric load, non-tidal ocean load and land hydrologic load and a thermoelastic displacement time sequence forcedly generated by temperature, summing the three-component displacement time sequence to obtain a non-tidal load displacement time sequence, unifying the non-tidal load displacement time sequence and the thermoelastic displacement time sequence into a reference frame and a time system consistent with GNSS (Global navigation satellite System) calculation, and forming a non-structural displacement vector under a site local coordinate system; The time interpolation and coordinate conversion module is used for performing time interpolation on the unstructured displacement vector to obtain a displacement increment of each observation epoch, and converting the displacement increment under the local coordinate system into a displacement increment under a geocentric rectangular coordinate system; The sight projection module is used for calculating the sight unit vector from the satellite to the reference point of the measuring station for each observation epoch and each satellite, and projecting the displacement increment under the geocentric rectangular coordinate system to the sight direction to form an equivalent geometric distance correction term; The observation correction module is used for applying the equivalent geometric distance correction term to pseudo-range observation and carrier phase observation in the same domain as the observation to obtain a corrected observation value; And the precise GNSS data processing module is used for executing precise GNSS arc segment calculation based on the corrected observation value and outputting a coordinate time sequence.
10. A non-transitory computer readable storage medium storing computer instructions that cause a computer to perform the steps of the method of any one of claims 1 to 8.

Description

Combined correction method, system and storage medium for GNSS non-structural deformation Technical Field The invention belongs to the field of satellite geodetic measurement and high-precision GNSS data processing, and relates to non-structural deformation modeling and correction in a GNSS coordinate time sequence, in particular to a processing method and a processing system for jointly introducing non-tidal load displacement and temperature-driven thermoelastic displacement into a GNSS observation equation for solving in an observation layer mode. Background Continuous GNSS reference stations are capable of providing millimeter-scale three-dimensional deformation monitoring capabilities and have been widely used for reference frame maintenance, structural deformation monitoring, and geophysical process research. However, the mass redistribution of the earth's surface such as the atmosphere, sea, land hydrology and the like can cause elastic load displacement, and meanwhile, the near-surface temperature change can drive the observation piers of the measuring station and the surrounding medium thereof to generate thermal expansion and contraction response. The non-structural deformation signals may reach a considerable level on an hourly to annual scale, and particularly in a regional dense station network, the common mode errors among stations and the spatially-correlated noise are more likely to appear, so that the repeatability of vertical coordinates, the speed estimation and the uncertainty thereof are affected. In the existing engineering application, the non-tidal load and the thermoelastic displacement are more commonly subjected to a mode of coordinate post-processing deduction, namely, coordinate estimation is finished according to a set observation model, and then the displacement of an external model is deducted on a daily coordinate time sequence layer. The method is simple and convenient to realize, but when the non-structural signal is not explicitly modeled in the parameter estimation process, the non-structural signal can be coupled with parameters such as station coordinates, troposphere delay, clock error, ambiguity and the like and partially absorbed, so that low-frequency system errors and common-mode components still remain after correction, and related noise and signal leakage are difficult to be fundamentally inhibited. Furthermore, in recent years, while the site-shift sequence of an hour resolution has been obtainable by a layered thermal conduction-elastic coupling model, its consistent introduction and joint correction path in the solution link still lacks a reusable implementation. Disclosure of Invention To overcome the above-mentioned shortcomings of the prior art, the present invention provides a method, system and storage medium for joint correction of GNSS non-structural deformations. The invention aims to provide a GNSS-oriented non-structural displacement observation layer joint correction method and a GNSS-oriented non-structural displacement observation layer joint correction system, wherein geometric distance correction is applied to epoch-by-epoch pseudo-range observation and carrier phase observation in an observation equation layer, so that non-structural deformation is explicitly interpreted before parameter estimation, coupling with parameters to be estimated and signal partial absorption are reduced, and vertical repeatability and spatial consistency of a dense station network calculation result are improved. According to an aspect of the present disclosure, there is provided a joint correction method for GNSS non-structural deformations, including: acquiring GNSS observation data of a station to be processed and a precision product, wherein the GNSS observation data comprises pseudo-range observation and carrier phase observation of at least one satellite system, and the precision product comprises an orbit, a clock error and a deviation product for fixing ambiguity; The construction of the non-structural displacement prior correction comprises the steps of obtaining or calculating a three-component displacement time sequence caused by non-tidal atmospheric load, non-tidal ocean load and land hydrologic load and a thermoelastic displacement time sequence forcedly generated by temperature, summing the three-component displacement time sequence to obtain a non-tidal load displacement time sequence, unifying the non-tidal load displacement time sequence and the thermoelastic displacement time sequence to a reference frame and a time system consistent with GNSS (Global navigation satellite System) calculation, and forming a non-structural displacement vector under a site local coordinate system; Performing time interpolation on the unstructured displacement vector to obtain a displacement increment of each observation epoch, and converting the displacement increment under the local coordinate system into a displacement increment under a geocentric rectangular coordi