CN114928428A - Precise time transfer method and system based on enhancement means - Google Patents
Precise time transfer method and system based on enhancement means Download PDFInfo
- Publication number
- CN114928428A CN114928428A CN202210585728.XA CN202210585728A CN114928428A CN 114928428 A CN114928428 A CN 114928428A CN 202210585728 A CN202210585728 A CN 202210585728A CN 114928428 A CN114928428 A CN 114928428A
- Authority
- CN
- China
- Prior art keywords
- data
- satellite
- reference station
- information
- user
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0635—Clock or time synchronisation in a network
- H04J3/0638—Clock or time synchronisation among nodes; Internode synchronisation
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to a precision time transmission method and a precision time transmission system based on an enhancement means. According to the precision time transmission method based on the enhancement means, firstly, based on the observation values of intensive regional enhancement stations, products such as a wide-area orbit, a clock error and a UPD are combined, after regional enhancement information is generated, a user can perform fine correction of a non-error after receiving the enhancement information, non-error PPP-RTK calculation is performed by combining the collected observation data and the received ephemeris and UPD products, and the time difference of each measurement station can be rapidly obtained, so that rapid precision time transmission among the users is realized, and the timeliness of the user is improved while the time transmission precision is improved.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a precision time transmission method and a precision time transmission system based on an enhancement means.
Background
The precise time transfer has wide application requirements in industries such as power networking, 5G communication, financial securities and the like. The time transfer method based on the satellite technology is one of the important means because of low cost, good continuity and easy implementation. The method of carrier phase Precise Point Positioning (PPP) time transfer in satellite time transfer is currently the mainstream technology due to high precision and has been widely used. However, the current PPP time transfer method has disadvantages that many errors in the non-differential observed value are difficult to process efficiently, and the time transfer precision can be further improved; secondly, the convergence time of the phase ambiguity in PPP is long, which affects the timeliness of the user.
How to further improve the precision and convergence speed of PPP time transfer has important value for popularizing the application of the technology and improving the user service experience.
Disclosure of Invention
The invention aims to provide a precision time transfer method and a precision time transfer system based on an enhancement means, which can improve the time transfer precision and the timeliness of the use of a user.
In order to achieve the purpose, the invention provides the following scheme:
a precision time transfer method based on an enhancement means, comprising:
acquiring wide area enhancement information, regional reference station data and resolving auxiliary data; the wide area augmentation information includes: wide-area satellite ephemeris, clock error and phase fractional deviation products; the regional reference station data includes: pseudo range and phase original observation data of a satellite on a reference station; the resolving assistance data includes: satellite antenna phase center parameters, earth rotation parameters and reference station coordinates;
extracting information of each satellite on a reference station based on the wide area augmentation information, the regional reference station data and the resolving auxiliary data to obtain regional augmentation information; the region enhancement information includes: pseudo range, phase, troposphere information and ionosphere information for each satellite;
acquiring user observation data; the user observation data comprises: user station measurement data and resolving auxiliary data of the user station measurement; the user station data comprises: pseudo-range and phase original observation data of a satellite on a user observation station;
and carrying out precise time transmission by adopting the user observation data and the region enhancement information.
Preferably, the extracting information of each satellite on the reference station based on the wide-area augmentation information, the regional reference station data, and the solution assistance data to obtain the regional augmentation information specifically includes:
performing quality inspection, cycle slip detection and gross error elimination on the regional reference station data to obtain first processed data;
and based on the satellite ephemeris and the reference station coordinates, correcting relativity, tide, an antenna phase center, a troposphere and earth rotation errors of the first processed data to obtain region enhancement information.
Preferably, the correcting relativity, tide, antenna phase center, troposphere and earth rotation error of the processed data based on the satellite ephemeris and the reference station coordinates to obtain the regional enhancement information specifically includes:
amending the relativity and the tide using the model specified in the IERS convections 2010;
correcting the antenna phase center by using an igs14.atx model;
correcting the troposphere by adopting a Saastamoinen model;
correcting the earth rotation error by adopting an IERS EOP C04 model;
and combining the phase decimal deviation product to carry out PPP solution and ambiguity fixing on the reference station, obtaining the clock error of the receiver after the PPP ambiguity is fixed, and further extracting the enhancement information of each satellite on each reference station based on the phase observation value and the satellite ephemeris to obtain the regional enhancement information.
Preferably, the performing precise time transfer by using the user observation data and the region enhancement information specifically includes:
obtaining user region enhancement information based on the user observation data;
determining an observation equation and a random model based on the user area enhancement information, and resolving the observation equation and the random model by adopting least square to obtain a floating solution;
and on the basis of the floating solution, performing single difference projection between the stars, recovering the whole cycle characteristic of the ambiguity, and fixing the ambiguity by adopting a Lambda search algorithm to realize precise time transfer of the ambiguity as a fixed solution.
Preferably, the corrected non-differential observation value of the atmosphere delay correction interpolated by the user station is expressed as:
in the formula (I), the compound is shown in the specification,
for the pseudoranges corrected for the atmospheric delay,
for the phase observation after the atmospheric delay correction,
for the tropospheric model value to be,
for the interpolated tropospheric correction values,
for the interpolated ionospheric corrections,
is the observation of the original pseudorange,
for raw phase observations, λ i A carrier wavelength at an ith frequency;
the observation equation after the enhanced information is corrected is as follows:
in the formula (I), the compound is shown in the specification,
for the floating-point ambiguities at the end of the streamer,
for the corrected receiver clock error, c is the speed of light,
for the pseudo-range to observe the noise,
for phase observation noise, ρ s Is the standing star geometric distance; i represents frequency and s represents satellite;
the stochastic model is determined according to the satellite altitude as:
in the formula, theta s Is the altitude angle of the s-th satellite (σ) s ) 2 Alpha is the cutoff satellite elevation angle, which is the prior variance of the observations.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
according to the precision time transmission method based on the enhancement means, firstly, based on the observation value of the intensive regional enhancement station, products such as a wide-area orbit, a clock error and a UPD are combined, after regional enhancement information is generated, a user can perform fine correction of a non-error after receiving the enhancement information, non-error PPP-RTK calculation is performed by combining the collected observation data and the received ephemeris and UPD products, time difference of each measurement station can be rapidly obtained, rapid precision time transmission among the users is achieved, time transmission precision is improved, and meanwhile timeliness of the user is improved.
Corresponding to the precision time transmission method based on the enhancement means, the invention also provides a precision time transmission system based on the enhancement means, which comprises the following steps:
the first data acquisition module is used for acquiring wide area enhancement information, regional reference station data and resolving auxiliary data; the wide area augmentation information includes: wide area satellite ephemeris, clock error and phase fractional deviation products; the regional reference station data includes: pseudo range and phase original observation data of a satellite on a reference station; the resolving assistance data includes: satellite antenna phase center parameters, earth rotation parameters and reference station coordinates;
the enhancement information generation module is used for extracting information of each satellite on a reference station based on the wide-area enhancement information, the regional reference station data and the resolving auxiliary data to obtain regional enhancement information; the region enhancement information includes: pseudo range, phase, troposphere information and ionosphere information for each satellite;
the second data acquisition module is used for acquiring user observation data; the user observation data comprises: user station measurement data and resolving auxiliary data of the user station measurement; the user station data comprises: pseudo-range and phase original observation data of a satellite on a user observation station;
and the time transmission module is used for carrying out precise time transmission by adopting the user observation data and the region enhancement information.
Preferably, the enhancement information generating module includes:
the first processing unit is used for carrying out quality inspection, cycle slip detection and gross error rejection processing on the regional reference station data to obtain first processing data;
and the first enhancement information generation unit is used for correcting relativity, tide, antenna phase center, troposphere and earth rotation error of the first processed data based on the satellite ephemeris and the reference station coordinates to obtain regional enhancement information.
Preferably, the enhancement information generating unit includes:
a first revision subunit for revising the relatives and the tides using the model specified in the IERS convections 2010;
a second correction subunit, configured to correct the antenna phase center using an igs14.atx model;
a third correction subunit, configured to correct the troposphere using a Saastamoinen model;
the fourth correcting subunit is used for correcting the earth rotation error by adopting an IERS EOP C04 model;
and the enhancement information generation subunit is used for combining the phase decimal deviation product to perform PPP solution and ambiguity fixing of the reference station, obtaining the clock error of the receiver after the PPP ambiguity is fixed, and further extracting the enhancement information of each satellite on each reference station based on the phase observation value and the satellite ephemeris to obtain the regional enhancement information.
Preferably, the time transfer module includes:
the second enhancement information generation unit is used for obtaining user area enhancement information based on the user observation data;
the least square resolving unit is used for determining an observation equation and a random model based on the user area enhancement information and resolving the observation equation and the random model by adopting least square to obtain a floating point solution;
and the time transfer unit is used for performing single difference projection between the stars on the basis of the floating solution, recovering the whole cycle characteristic of the ambiguity, and fixing the ambiguity by adopting a Lambda search algorithm to realize precise time transfer with the ambiguity as a fixed solution.
The technical effect achieved by the precision time transmission system based on the enhancement means provided by the invention is the same as that achieved by the precision time transmission method based on the enhancement means provided by the invention, so that the details are not repeated herein.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a flow chart of a precision time transfer method based on an enhancement means provided by the present invention;
FIG. 2 is a block diagram of an embodiment of a precision time transfer method based on an enhancement method provided by the present invention;
fig. 3 is a schematic structural diagram of a precision time transfer system based on an enhancement means provided in the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a precision time transfer method and a precision time transfer system based on an enhancement means, which can improve the time transfer precision and the timeliness of the use of a user.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the method for precision time transfer based on enhancement means provided by the present invention comprises:
step 100: and acquiring wide area enhancement information, regional reference station data and resolving auxiliary data. The wide area augmentation information includes: wide area satellite ephemeris, clock error and phase fractional deviation product (UPD). The invention mainly obtains real-time orbit, clock error and phase decimal deviation correction (UPD) products provided by the GNSS product center of the national time service center of the Chinese academy of sciences in real time through the Internet. The regional reference station data includes: pseudoranges and phase raw observations of satellites at a reference station. Resolving the assistance data includes: the satellite antenna phase center parameters, the earth rotation parameters and the reference station coordinates.
Step 101: and extracting information of each satellite on the reference station based on the wide-area augmentation information, the regional reference station data and the resolving auxiliary data to obtain regional augmentation information. The region enhancement information includes: pseudoranges, phases, tropospheric information and ionospheric information for each satellite.
For example, the implementation process of this step is: firstly, performing quality inspection, cycle slip detection and gross error elimination on original data (pseudo range and phase observation values) of a reference station, and deleting data without satellite ephemeris or incomplete observation values to obtain clean data. Based on precise ephemeris and strong constraint coordinates, the preprocessed clean data is corrected by relativity, tide, antenna phase center, troposphere and earth rotation error, wherein the relativity and tide correction uses the model correction specified in IERS convections 2010, the antenna phase center correction uses igs14.atx model correction, the troposphere correction uses Saastamoinen model correction, the earth rotation error correction uses IERS EOP C04 model correction, and then PPP solution and ambiguity fixing of the reference station are carried out by combining with a UPD product, after the PPP ambiguity fixing, the clock error of the receiver can be obtained, and further the enhanced information of each satellite on each reference station can be extracted based on the phase observation value and the ephemeris product.
The tropospheric correction Delta T on the inclined path is the calculated tropospheric value T s Minus the model quantity T model As shown in formula (1):
ΔT＝T s -T model (1)
here, T s And T model The calculation method of (c) is as follows:
T model ＝d dry F dry +d wet F wet (3)
wherein, lambda is carrier wave wavelength, phi is phase observation value, rho is geometric distance between satellite and station, t is receiver clock error, t is phase difference s Is the satellite clock error, c is the speed of light,
for phase ambiguity, s denotes satellite, IF denotes ionospheric-free combination, d dry And d wet Delay of dry and wet components of troposphere in zenith direction, F dry And F wet Respectively, corresponding projection coefficients.
The ionospheric delay corrections of pseudoranges and phases on the ramp path can be calculated by equation (4):
wherein, P and phi represent pseudo range and phase observation value, rho is station star geometric distance, delta P,i And
for receiver-side and satellite-side hardware delays, delta Φ,i And
respectively, the fractional phase offset at the receiver and satellite. And the hardware delay deviation and the phase decimal deviation of the satellite end are corrected by adopting a product, and the hardware delay and the phase decimal deviation of the receiver end are not considered and are fused in the enhancement information.
After the regional enhancement information is obtained, the pseudo-range, troposphere and ionosphere enhancement information of the phase of each satellite is broadcasted in a set format through a network.
Step 102: and acquiring user observation data. Specifically, pseudo range and phase observation data of a satellite on a user observation station, a wide-area satellite orbit, clock error and UPD (universal up-draft) product, enhancement information of each reference station and auxiliary resolving parameters including an antenna phase center parameter, an earth rotation parameter and an accurate coordinate of the observation station are obtained.
Step 103: and carrying out precise time transfer by adopting user observation data and regional enhancement information.
For example, in implementation, the ambiguity based regional enhancement approach is a time-transfer solution of the floating solution. Firstly, the data of the user station is subjected to quality inspection, cycle slip detection and gross error elimination, and the data without satellite ephemeris or incomplete observation values are deleted to obtain clean data. Next, error correction was performed, where relativistic and tidal corrections were corrected using the model specified in the IERS convections 2010, antenna phase center correction was corrected using the igs14.atx model, troposphere correction was corrected using the Saastamoinen model, earth rotation error correction was corrected using the IERS EOP C04 model, and interpolation and fine correction of the enhancement information were performed. And finally, determining an observation equation and a stochastic model, and performing least square solution to realize precise time transfer of the ambiguity of the user side as a floating solution.
The non-difference observation value after the correction of the atmospheric delay correction interpolated at the user terminal can be expressed as:
in the formula (I), the compound is shown in the specification,
and
respectively, the atmospheric delay corrected pseudorange and phase observations,
for the tropospheric model value to be,
for the raw pseudorange observations,
for raw phase observations, λ i Is the carrier wavelength of the ith frequency,
values are corrected for the interpolated tropospheric slopes,
for the interpolated ionospheric corrections, the calculation formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,
for the extracted tropospheric delays of the respective satellites,
is the extracted ionospheric delay for each satellite. a is a r For the interpolation coefficient, the calculation formula is as follows:
in the formula, x ru Is the difference in coordinates of the reference station and the subscriber station, a r For interpolating coefficients, (X) r ,Y r ) And (X) u ,Y u ) Respectively, the plane coordinates of the reference station and the user in the area coordinate system.
The observation equation after the enhancement information is corrected can be written as:
in the formula (I), the compound is shown in the specification,
floating ambiguity is the streaming station end.
And
representing the corrected pseudorange and phase observations,
for the corrected receiver clock error, c is the speed of light,
for the pseudo-range to observe the noise,
for phase observation noise, p s Is the standing star geometric distance; i represents frequency and s represents satellite.
The stochastic model is determined according to the satellite altitude as follows:
here, θ s Is the altitude angle, σ, of the s-th satellite s For standard deviation of observations, the observation is typically set to 0.3 meters for pseudorange, the observation is typically set to 0.002 meters for phase, and α is the cutoff satellite elevation angle, typically set to 30 degrees.
And finally, carrying out precise time transfer calculation by taking the ambiguity as a fixed solution. Specifically, on the basis of a floating solution, single difference projection between stars is carried out, the whole-cycle characteristic of the ambiguity is recovered, the ambiguity is fixed by adopting a Lambda search algorithm, and precise time transfer with the ambiguity as a fixed solution is realized.
Based on the above description, as shown in fig. 2, the precision time transfer method based on the enhancement means provided by the present invention mainly includes eight parts, namely, wide-area enhanced product acquisition, regional reference station data acquisition, calculation auxiliary product acquisition, enhanced information extraction, enhanced information dissemination, user data acquisition, time transfer calculation based on the ambiguity of the enhanced information as a floating solution, and time transfer calculation based on the ambiguity as a fixed solution.
Compared with the prior art, the precision time transfer method based on the enhancement means also has the following advantages:
first, the present invention eliminates or attenuates common errors, improves the observed value error correction level, and improves the time transfer accuracy based on the enhancement information. And the high-precision enhanced information provided by the reference station, including some unmodeled errors, can be used for improving the correction level of the observed value errors of the user station.
Second, the phase ambiguity is fixed, allowing for fast convergence of time transfer services. Compared with the conventional PPP time transfer method, the ambiguity adopts a floating solution, the method and the device rely on the enhancement information to carry out the normalization and the fixation of the ambiguity, accelerate the parameter convergence speed and realize the quick and precise time transfer of the user.
In addition, corresponding to the above-mentioned precision time transmission method based on enhancement means, the present invention also provides a precision time transmission system based on enhancement means, as shown in fig. 3, the system includes: the device comprises a first data acquisition module 1, an enhanced information generation module 2, a second data acquisition module 3 and a time transmission module 4.
The first data acquisition module 1 is used for acquiring wide area augmentation information, regional reference station data and resolving auxiliary data. The wide area augmentation information includes: wide area satellite ephemeris, clock error and phase fractional deviation products. The regional reference station data includes: the pseudoranges and phases of the satellites at the reference station are the raw observations. Resolving the assistance data includes: satellite antenna phase center parameters, earth rotation parameters and reference station coordinates.
The enhanced information generation module 2 is configured to extract information of each satellite on the reference station based on the wide-area enhanced information, the regional reference station data, and the solution assistance data to obtain regional enhanced information. The region enhancement information includes: pseudoranges, phases, tropospheric information and ionospheric information for each satellite.
The second data acquiring module 3 is used for acquiring user observation data. The user observation data includes: user station data and resolving auxiliary data of the user station. The user station data comprises: pseudoranges and phase raw observations of satellites at user stations.
The time transmission module 4 is used for performing precise time transmission by adopting user observation data and regional enhancement information.
In order to further improve the accuracy of obtaining the regional enhancement information, as an embodiment of the present invention, the enhancement information generating module 2 provided above may include: the device comprises a first processing unit and a first enhancement information generating unit.
The first processing unit is used for carrying out quality inspection, cycle slip detection and gross error rejection processing on the regional reference station data to obtain first processing data.
The first enhancement information generation unit is used for correcting relativity, tide, antenna phase center, troposphere and earth rotation error of the first processed data based on the satellite ephemeris and the reference station coordinates to obtain regional enhancement information.
Wherein the enhancement information generation unit includes: the device comprises a first correction subunit, a second correction subunit, a third correction subunit, a fourth correction subunit and an enhancement information generation subunit.
The first correction subunit is used to correct relativity and tides using the model specified in the IERS considerations 2010.
The second correction subunit is used to correct the antenna phase center using the igs14.atx model.
The third correction subunit is used for correcting the troposphere by adopting a Saastamoinen model.
The fourth correction subunit is used for correcting the earth rotation error by adopting an IERS EOP C04 model.
The enhancement information generation subunit is configured to perform PPP solution and ambiguity fixing for the reference station in combination with the phase decimal deviation product, obtain a clock error of the receiver after the PPP ambiguity fixing, and further extract enhancement information of each satellite on each reference station based on the phase observation value and the satellite ephemeris to obtain regional enhancement information.
Further to improve the precision of the precise time transmission, as another embodiment of the present invention, the time transmission module 4 adopted above may include: the device comprises a second enhancement information generation unit, a least square calculation unit and a time transmission unit.
The second enhancement information generation unit is used for obtaining user area enhancement information based on the user observation data.
The least square solving unit is used for determining an observation equation and a random model based on the user area enhancement information and adopting least square to solve the observation equation and the random model to obtain a floating point solution.
The time transfer unit is used for performing single difference projection between the satellites on the basis of the floating point solution, recovering the whole cycle characteristic of the ambiguity, fixing the ambiguity by adopting a Lambda search algorithm, and realizing precise time transfer of the ambiguity which is a fixed solution.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the description of the method part.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the foregoing, the description is not to be taken in a limiting sense.
Claims (9)
1. A precision time transfer method based on an enhancement means is characterized by comprising the following steps:
acquiring wide area enhancement information, regional reference station data and resolving auxiliary data; the wide area augmentation information includes: wide area satellite ephemeris, clock error and phase fractional deviation products; the regional reference station data includes: pseudo range and phase original observation data of a satellite on a reference station; the resolving assistance data includes: satellite antenna phase center parameters, earth rotation parameters and reference station coordinates;
extracting information of each satellite on a reference station based on the wide area augmentation information, the regional reference station data and the resolving auxiliary data to obtain regional augmentation information; the region enhancement information includes: pseudo range, phase, troposphere information and ionosphere information for each satellite;
acquiring user observation data; the user observation data comprises: user station measurement data and resolving auxiliary data of the user station measurement; the user station data comprises: pseudo-range and phase original observation data of a satellite on a user observation station;
and carrying out precise time transmission by adopting the user observation data and the region enhancement information.
2. The precision time transfer method based on the augmentation instrument of claim 1, wherein the extracting information of each satellite on a reference station based on the wide-area augmentation information, the regional reference station data and the solution assistance data to obtain regional augmentation information specifically comprises:
performing quality inspection, cycle slip detection and gross error elimination on the regional reference station data to obtain first processed data;
and based on the satellite ephemeris and the reference station coordinates, correcting relativity, tide, an antenna phase center, a troposphere and earth rotation errors of the first processed data to obtain region enhancement information.
3. The method for precision time transfer based on the augmentation instrument as claimed in claim 2, wherein the correction of relativity, tide, antenna phase center, troposphere and earth rotation error is performed on the processed data based on the satellite ephemeris and the reference station coordinates to obtain the regional augmentation information, specifically comprising:
modifying the relatives and the tides using the model specified in the IERS convections 2010;
correcting the antenna phase center by using an igs14.atx model;
correcting the troposphere by adopting a Saastamoinen model;
correcting the earth rotation error by adopting an IERS EOP C04 model;
and combining the phase decimal deviation product to carry out PPP solution and ambiguity fixing on the reference station, obtaining the clock error of the receiver after the PPP ambiguity is fixed, and further extracting the enhancement information of each satellite on each reference station based on the phase observation value and the satellite ephemeris to obtain the regional enhancement information.
4. The precision time transfer method based on the enhancement means as claimed in claim 1, wherein the precision time transfer using the user observation data and the region enhancement information specifically comprises:
obtaining user region enhancement information based on the user observation data;
determining an observation equation and a random model based on the user area enhancement information, and resolving the observation equation and the random model by adopting least square to obtain a floating solution;
and on the basis of the floating solution, performing single difference projection between the stars, recovering the whole cycle characteristic of the ambiguity, and fixing the ambiguity by adopting a Lambda search algorithm to realize precise time transfer of the ambiguity as a fixed solution.
5. The enhanced means-based precision time transfer method of claim 4, wherein the user station interpolated atmosphere delay correction corrected non-differential observation value is represented as:
in the formula (I), the compound is shown in the specification,
for the pseudoranges corrected for the atmospheric delay,
for the phase observation after the atmospheric delay correction,
for the tropospheric model value to be,
for the interpolated tropospheric correction values,
for the interpolated ionospheric corrections,
for the raw pseudorange observations,
for raw phase observations, λ i A carrier wavelength at an ith frequency;
the observation equation after the enhanced information is corrected is as follows:
in the formula (I), the compound is shown in the specification,
for the floating-point ambiguities at the end of the streamer,
for the corrected receiver clock error, c is the speed of light,
for the pseudo-range to observe the noise,
for phase observation noise, p s Is the standing star geometric distance; i represents frequency and s represents satellite;
the stochastic model is determined as follows according to the satellite altitude:
in the formula, theta s Is the altitude angle of the s-th satellite, (σ) s ) 2 Alpha is the cutoff satellite elevation angle, which is the prior variance of the observations.
6. A precision time transfer system based on an enhancement means, comprising:
the first data acquisition module is used for acquiring wide area enhancement information, regional reference station data and resolving auxiliary data; the wide area augmentation information includes: wide-area satellite ephemeris, clock error and phase fractional deviation products; the regional reference station data includes: pseudo range and phase original observation data of a satellite on a reference station; the resolving assistance data includes: satellite antenna phase center parameters, earth rotation parameters and reference station coordinates;
the enhancement information generation module is used for extracting information of each satellite on a reference station based on the wide-area enhancement information, the regional reference station data and the resolving auxiliary data to obtain regional enhancement information; the region enhancement information includes: pseudo range, phase, troposphere information and ionosphere information of each satellite;
the second data acquisition module is used for acquiring user observation data; the user observation data comprises: user station measurement data and resolving auxiliary data of the user station measurement; the user station data comprises: pseudo-range and phase original observation data of a satellite on a user observation station;
and the time transmission module is used for carrying out precise time transmission by adopting the user observation data and the region enhancement information.
7. The enhanced instrument based precision time delivery system of claim 6, wherein the enhanced information generation module comprises:
the first processing unit is used for carrying out quality inspection, cycle slip detection and gross error rejection processing on the regional reference station data to obtain first processing data;
and the first enhancement information generation unit is used for correcting relativity, tide, antenna phase center, troposphere and earth rotation error of the first processed data based on the satellite ephemeris and the reference station coordinates to obtain regional enhancement information.
8. The enhanced instrument-based precision time transfer system according to claim 7, wherein the enhancement information generating unit includes:
a first amendment subunit for amending the relativity and the tide using the model specified in the IERS convections 2010;
a second correction subunit, configured to correct the antenna phase center using an igs14.atx model;
a third correction subunit, configured to correct the troposphere using a Saastamoinen model;
the fourth correcting subunit is used for correcting the earth rotation error by adopting an IERS EOP C04 model;
and the enhancement information generation subunit is used for combining the phase decimal deviation product to perform PPP solution and ambiguity fixing of the reference station, obtaining the clock error of the receiver after the PPP ambiguity is fixed, and further extracting the enhancement information of each satellite on each reference station based on the phase observation value and the satellite ephemeris to obtain the regional enhancement information.
9. The enhanced utility based precision time transfer system of claim 6, wherein said time transfer module comprises:
the second enhancement information generation unit is used for obtaining user area enhancement information based on the user observation data;
the least square resolving unit is used for determining an observation equation and a random model based on the user area enhancement information and resolving the observation equation and the random model by adopting least square to obtain a floating point solution;
and the time transfer unit is used for performing single difference projection between the stars on the basis of the floating solution, recovering the whole cycle characteristic of the ambiguity, and fixing the ambiguity by adopting a Lambda search algorithm to realize precise time transfer with the ambiguity as a fixed solution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210585728.XA CN114928428A (en) | 2022-05-26 | 2022-05-26 | Precise time transfer method and system based on enhancement means |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210585728.XA CN114928428A (en) | 2022-05-26 | 2022-05-26 | Precise time transfer method and system based on enhancement means |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114928428A true CN114928428A (en) | 2022-08-19 |
Family
ID=82809741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210585728.XA Pending CN114928428A (en) | 2022-05-26 | 2022-05-26 | Precise time transfer method and system based on enhancement means |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114928428A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115225245A (en) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | Non-differential non-combination PPP-RTK long baseline real-time transfer method |
| CN116540282A (en) * | 2023-05-08 | 2023-08-04 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
-
2022
- 2022-05-26 CN CN202210585728.XA patent/CN114928428A/en active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115225245A (en) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | Non-differential non-combination PPP-RTK long baseline real-time transfer method |
| CN116540282A (en) * | 2023-05-08 | 2023-08-04 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
| CN116540282B (en) * | 2023-05-08 | 2023-09-29 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111045034B (en) | GNSS multi-system real-time precise time transfer method and system based on broadcast ephemeris | |
| CN109581452B (en) | GNSS reference station carrier phase integer ambiguity resolution method | |
| CN114928428A (en) | Precise time transfer method and system based on enhancement means | |
| CN111580132B (en) | Beidou local precise time transfer method based on time laboratory enhancement information | |
| CN111856534B (en) | Dual-mode GNSS carrier precise single-point positioning method and system of intelligent terminal | |
| CN109001781A (en) | A kind of tri- frequency Ambiguity Solution Methods of BDS for taking ionosphere constraint into account | |
| CN113848572A (en) | Multi-frequency PPP sequential single epoch positioning method based on atmospheric error enhancement | |
| CN110824505B (en) | Deviation estimation method and system, positioning method and terminal of GNSS satellite receiver | |
| CN113325446B (en) | Multimode common-frequency GNSS carrier phase time transfer method and system | |
| CN116243591B (en) | Subnanosecond time service method integrating UTC (k) and Beidou broadcast ephemeris | |
| CN108535749B (en) | Positioning enhancement method and system based on CORS and positioning system | |
| CN114966760A (en) | Ionosphere weighted non-differential non-combination PPP-RTK technology implementation method | |
| CN116203598A (en) | Ionosphere modeling method, device and medium based on foundation and star-based enhanced fusion | |
| CN114879239B (en) | Regional three-frequency integer clock error estimation method for enhancing instantaneous PPP fixed solution | |
| CN114721012B (en) | Ionized layer delay effect correction method and device | |
| CN110988917B (en) | Real-time monitoring method for satellite orbit maneuvering state | |
| CN111562594A (en) | Beidou precision time transfer method based on interstation single difference and ambiguity fixation | |
| CN115327593B (en) | Positioning method, system and storage medium based on unmanned aerial vehicle | |
| CN117665871A (en) | GNSS non-difference fixed solution precision clock difference determination method under carrier distance support | |
| CN115933356B (en) | High-precision time synchronization system and method for virtual atomic clock | |
| CN115856945B (en) | Precise single-point time service method, device and medium suitable for GNSS multi-frequency points | |
| CN114355420B (en) | PPP product positioning method and device for distributed Beidou position service center | |
| CN115144878B (en) | PPP-based short-distance large-height-difference NRTK troposphere delay correction method | |
| CN111551975A (en) | BDS/GPS reference station low-altitude angle satellite whole-cycle ambiguity determination method | |
| CN114545470B (en) | Cross-network continuous network RTK positioning service method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |