CN115356776B - Sea tide mixing error estimation method in satellite gravity measurement - Google Patents
Sea tide mixing error estimation method in satellite gravity measurement Download PDFInfo
- Publication number
- CN115356776B CN115356776B CN202210985887.9A CN202210985887A CN115356776B CN 115356776 B CN115356776 B CN 115356776B CN 202210985887 A CN202210985887 A CN 202210985887A CN 115356776 B CN115356776 B CN 115356776B
- Authority
- CN
- China
- Prior art keywords
- frequency
- mixing
- gravity
- satellite
- signal
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000005484 gravity Effects 0.000 title claims abstract description 57
- 238000005259 measurement Methods 0.000 title claims abstract description 21
- 238000005070 sampling Methods 0.000 claims abstract description 39
- 230000000737 periodic effect Effects 0.000 claims description 8
- 230000003068 static effect Effects 0.000 claims description 7
- 230000003252 repetitive effect Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses a sea tide mixing error estimation method in satellite gravity measurement, which comprises the following steps: determining a satellite orbit sampling process of a gravity satellite task; determining inversion period sampling processes of different gravity inversion processes; calculating the mixing frequency of each sea tide component according to the satellite orbit sampling process and the inversion period sampling process; the mixing error of the high-frequency sea tide signal in the satellite gravity measurement is estimated. The gravity satellite task provided by the invention has two processes, namely, the first process is orbital sampling and the second process is gravity period inversion sampling, the first mixing frequency and the second mixing frequency can be estimated by determining the signal frequency and the sampling frequency of each process, and further, the mixing error can be estimated according to the mixing frequency, so that a practical method is provided for estimating the mixing error of the sea tide signal in the satellite gravity task.
Description
Technical Field
The invention relates to the technical field of satellite gravity, in particular to a method for estimating a mixing error of a high-frequency sea tide signal in a satellite gravity task.
Background
The phenomenon of mixing a high frequency sea-tide signal into a low frequency due to undersampling of the signal is called sea-tide mixing. During the process of moving around the earth and acquiring observation data, the gravity satellites record geophysical signals of different spatial scales and different time frequencies on the earth in a single observation value. The continuous satellite observation can effectively reflect a static gravitational field and a low-frequency time-varying gravitational field caused by seasonal changes and the like. For the high-frequency sea tide time-varying signal, the satellite is sampled into a low-frequency signal due to undersampling in the satellite observation and data processing processes, so that a mixing error is generated. Sea tide mixing errors severely restrict satellite gravity measurement inversion accuracy. However, there is currently no effective method to estimate the sea tide mixing error.
Disclosure of Invention
In order to solve the technical problems, the invention aims to estimate the sea tide mixing error according to a mixing mechanism of a high-frequency sea tide signal in satellite gravity measurement and provides an estimation method for estimating the sea tide mixing error in satellite gravity measurement observation.
The aim of the invention is achieved by the following technical scheme:
a method for estimating sea tide mixing errors in satellite gravity measurement, comprising:
A. determining a satellite orbit sampling process of a gravity satellite task;
B. determining inversion period sampling processes of different gravity inversion processes;
C. calculating the mixing frequency of each sea tide component according to the satellite orbit sampling process and the inversion period sampling process;
D. the mixing error of the high-frequency sea tide signal in the satellite gravity measurement is estimated.
One or more embodiments of the present invention may have the following advantages over the prior art:
the invention discloses an estimation method of a high-frequency signal mixing error in a satellite gravity task, which is characterized in that the satellite gravity task is decomposed into two processes of orbit sampling and period inversion sampling, the sampling frequency and the signal frequency of each process are respectively determined, the mixing frequency of each process is further calculated, and finally the mixing error is estimated according to the mixing frequency, so that a practical method is provided for estimating the mixing error.
Drawings
Fig. 1 is a flow chart of a method for estimating sea tide mixing errors in satellite gravity measurement.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the following examples and the accompanying drawings.
As shown in FIG. 1, the sea tide mixing error estimation method in satellite gravity measurement comprises
Step 1, determining a satellite orbit sampling process of a gravity satellite task, which specifically comprises the following steps:
step 1.1 determining a general sampling frequency of a general track (including a recurring periodic track and a non-recurring periodic track):
for general orbits, including repetitive periodic orbits and non-repetitive periodic orbits, the orbit sampling frequency can be extracted from the satellite orbit parameters, provided that the orbit plane precesses at a rate of inertiaThe rotation angular velocity of the earth is omega E The track-to-ground sampling frequency is the precession rate of the track surface relative to the earth's surface,
step 1.2, determining the signal frequency of satellite gravity mission orbit sampling;
different sea tide components have different frequencies, the sea tide component signal frequency is determined according to the signal range of the selected high-frequency sea tide signal, and is recorded as f i Wherein i=1, 2,3 … N, different sea tide components, denoted by fThere is a sea tide component frequency set.
Step 2, determining inversion period sampling processes of different gravity inversion processes, which specifically comprises the following steps:
step 2.1, determining the sampling signal frequency during satellite gravity inversion: the satellite gravity inversion process is to resample the orbit observation signal, when the orbit observation undersamples the sea tide signal, the resampled signal frequency is not the frequency of the geophysical signal itself, but the mixing frequency caused by the orbit undersampling
Step 2.2 determining the sampling frequency of gravity cycle inversion: when the gravity inversion process is strictly continuous over a certain time frame and per T r When the length of the sky data is one inversion period, the sampling frequency of gravity period inversion is 1/T of the reciprocal of the inversion period r 。
Step 3, estimating mixing frequency of heavy high-frequency sea tide signals measured by satellite gravity: according to the satellite orbit sampling process and the inversion period sampling process, the mixing frequency of each sea tide component is calculated, and the method specifically comprises the following steps:
step 3.1 according to the satellite orbit sampling procedure of mixing the high frequency signal, the first mixing frequency can be estimated
Wherein f is the frequency of the sea tide signal,for the precession rate of the orbit surface with respect to the earth, N is any integer that minimizes this equation.
Step 3.2 according to the gravity inversion sampling process of satellite gravity measurement, the second mixing frequency can be estimated
wherein ,for the first mixing frequency, T r For inversion period, N is any integer that minimizes this equation.
Step 4, utilizing the mixing error of the high-frequency sea tide signal in the estimated satellite gravity measurement, and specifically comprising the following steps:
and 4.1, extracting a gravitational field model coefficient time sequence. Extracting the time sequence of spherical harmonic coefficients of each gravity field according to the inversion result of the gravity field, and synchronously extracting time marks;
step 4.2, calculating a time sequence of the spherical harmonic coefficient residual error of the time-varying gravitational field; the method specifically comprises the following steps:
and 4.2.1, determining and subtracting a static gravitational field or a gravitational field mean value from the inverted gravitational field time sequence to obtain a time-varying gravitational field time sequence. The static gravitational field can use the existing static gravitational field model, and the gravitational field mean value can be obtained by averaging the gravitational field time sequence;
step 4.2.2 subtracting the dominant time-varying gravitational field signal from the time-varying gravitational field time sequence to obtain a residual time sequence expressed as ΔC in terms of spherical harmonic coefficients lm (f)、ΔS lm (t). The main time-varying signals include annual signals, semi-annual signals, and other time-varying signals of definite physical significance. If the quantity period and magnitude information of the time-varying signal are known, the quantity period and magnitude information can be directly subtracted from the time sequence of the time-varying gravitational field; if only the time-varying signal period is known, its magnitude can be estimated from the period using least squares and then subtracted from the time-varying gravitational field time series.
Step 4.3, taking the residual time sequence as input data, and estimating a mixing error by utilizing a first mixing frequency and a second mixing frequency of the high-frequency signal mixing;
specifically, residual time sequence ΔC is used lm (t)、ΔS lm (t) estimating a mixing error from the mixing frequency, the observation equation being
Wherein t is the time in which,for mixing frequencies, including a first mixing frequency->And a second mixing frequency->K is the total number of all mixing frequencies, +.>For each mixing frequency +.>Corresponding parameters to be estimated.
Will estimate the resultFrequency mixing error spherical harmonic coefficient time sequence obtained by bringing the right end of the frequency mixing error spherical harmonic coefficient time sequence into the right end of the frequency mixing error spherical harmonic coefficient time sequence>
Mixing error time seriesThe method is characterized in that the method is used for expressing spherical harmonic coefficients of sea tide mixing errors in each inversion period, and sea tide mixing errors at any position in the inversion period can be obtained through spherical harmonic synthesis.
Although the embodiments of the present invention are described above, the embodiments are only used for facilitating understanding of the present invention, and are not intended to limit the present invention. Any person skilled in the art can make any modification and variation in form and detail without departing from the spirit and scope of the present disclosure, but the scope of the present disclosure is still subject to the scope of the appended claims.
Claims (6)
1. A method for estimating sea tide mixing errors in satellite gravity measurement, comprising:
A. determining a satellite orbit sampling process of a gravity satellite task;
B. determining inversion period sampling processes of different gravity inversion processes;
C. calculating the mixing frequency of each sea tide component according to the satellite orbit sampling process and the inversion period sampling process;
D. estimating the mixing error of the high-frequency sea tide signal in satellite gravity measurement;
d1, extracting a gravity field model coefficient time sequence; extracting the time sequence of spherical harmonic coefficients of each gravity field according to the inversion result of the gravity field, and synchronously extracting time marks;
d2, calculating a time sequence of spherical harmonic coefficient residual errors of the time-varying gravitational field; the method specifically comprises the following steps:
d2.1, determining and subtracting a static gravitational field or a gravitational field mean value from the inverted gravitational field time sequence to obtain a time-varying gravitational field time sequence; the static gravitational field can use the existing static gravitational field model, and the gravitational field mean value can be obtained by averaging the gravitational field time sequence;
d2.2 subtracting the main time-varying gravitational field signal from the time-varying gravitational field time sequence to obtain a residual time sequence expressed as ΔC by spherical harmonic coefficients lm (t)、ΔS lm (t); the main time-varying signals comprise annual signals, semi-annual signals and other time-varying signals with definite physical significance; if the quantity period and magnitude information of the time-varying signal are known, the quantity period and magnitude information can be directly subtracted from the time sequence of the time-varying gravitational field; if only the time-varying signal period is known, then one can first rootEstimating the magnitude of the time-varying gravitational field by utilizing least square according to the period, and deducting the magnitude from the time-varying gravitational field time sequence;
d3, taking the residual time sequence as input data, and estimating a mixing error by utilizing a first mixing frequency and a second mixing frequency of high-frequency signal mixing;
specifically, residual time sequence ΔC is used lm (t)、ΔS lm (t) estimating a mixing error from the mixing frequency, the observation equation being
Wherein t is the time in which,for mixing frequencies, including a first mixing frequency->And a second mixing frequency->K is the total number of all mixing frequencies, +.>For each mixing frequency +.>Corresponding parameters to be estimated;
will estimate the resultFrequency mixing error spherical harmonic coefficient time sequence obtained by bringing the right end of the frequency mixing error spherical harmonic coefficient time sequence into the right end of the frequency mixing error spherical harmonic coefficient time sequence>
Mixing error time seriesThe method is characterized in that the method is used for expressing spherical harmonic coefficients of sea tide mixing errors in each inversion period, and sea tide mixing errors at any position in the inversion period can be obtained through spherical harmonic synthesis.
2. The method for estimating sea tide mixing error in satellite gravity measurement according to claim 1, wherein a comprises:
a1, determining a general sampling frequency containing a repeated periodic track and a non-repeated periodic track;
a2, determining the signal frequency of satellite gravity mission orbit sampling.
3. The method for estimating sea tide mixing error in satellite gravity measurement according to claim 2, wherein said A1 is provided with a track plane including a repetitive periodic track and a non-repetitive periodic track having a precession rate under an inertial system ofThe rotation angular velocity of the earth is omega E The sampling frequency of the orbit to the ground is the precession rate of the orbit plane relative to the earth's surface +.>
4. The method for estimating sea tide mixing error in satellite gravity measurement according to claim 2, wherein said A2 specifically comprises selecting a signal range of a high frequency sea tide signal, determining a sea tide component signal frequency f i Where i=1, 2, 3..n is the different sea tide components and f represents all sea tide component frequency sets.
5. The method for estimating sea tide mixing error in satellite gravity measurement according to claim 1, wherein B specifically comprises:
b1, determining the frequency of a sampling signal during satellite gravity inversion;
the satellite gravity inversion process is to resample the orbit observation signal, when the orbit observation undersamples the sea tide signal, the resampled signal frequency is not the frequency of the geophysical signal itself, but the mixing frequency caused by the orbit undersampling
B2, determining sampling frequency of gravity period inversion;
when the gravity inversion process is strictly continuous over a certain time frame and per T r When the length of the sky data is one inversion period, the sampling frequency of gravity period inversion is 1/T r 。
6. The method for estimating sea tide mixing error in satellite gravity measurement according to claim 1, wherein C specifically comprises:
c1 estimates a first mixing frequency according to a satellite orbit sampling process of mixing high-frequency signals
Wherein f is the frequency of the sea tide signal,n is any integer that minimizes the precession rate of the orbit surface relative to the earth;
c2 can estimate the second mixing frequency according to the gravity inversion sampling process of satellite gravity measurement
wherein ,for the first mixing frequency, T r For inversion period, N is any integer that minimizes this equation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210985887.9A CN115356776B (en) | 2022-08-16 | 2022-08-16 | Sea tide mixing error estimation method in satellite gravity measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210985887.9A CN115356776B (en) | 2022-08-16 | 2022-08-16 | Sea tide mixing error estimation method in satellite gravity measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115356776A CN115356776A (en) | 2022-11-18 |
| CN115356776B true CN115356776B (en) | 2023-10-10 |
Family
ID=84002532
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210985887.9A Active CN115356776B (en) | 2022-08-16 | 2022-08-16 | Sea tide mixing error estimation method in satellite gravity measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115356776B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011157686A1 (en) * | 2010-06-14 | 2011-12-22 | Universita' Degli Studi Di Roma "La Sapienza" | System for measuring coseismic movements or vibrations of structures based on global navigation satellite systems-gnss and/or pseudolites |
| CN105783884A (en) * | 2016-03-03 | 2016-07-20 | 南京信息工程大学 | Method for recognizing tidal mixing frontal surface with satellite altimeter along-track data |
| CN106997061A (en) * | 2017-04-05 | 2017-08-01 | 中国空间技术研究院 | A kind of method that relative velocity improves gravitational field inversion accuracy between star based on disturbance |
| CN110082787A (en) * | 2019-04-11 | 2019-08-02 | 华东师范大学 | A method of extracting half Sunday of Sunday tide signal from GNSS time sequence |
| CN112278322A (en) * | 2020-10-26 | 2021-01-29 | 中国空间技术研究院 | Marine satellite regression orbit determination method and device and storage medium |
-
2022
- 2022-08-16 CN CN202210985887.9A patent/CN115356776B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011157686A1 (en) * | 2010-06-14 | 2011-12-22 | Universita' Degli Studi Di Roma "La Sapienza" | System for measuring coseismic movements or vibrations of structures based on global navigation satellite systems-gnss and/or pseudolites |
| CN105783884A (en) * | 2016-03-03 | 2016-07-20 | 南京信息工程大学 | Method for recognizing tidal mixing frontal surface with satellite altimeter along-track data |
| CN106997061A (en) * | 2017-04-05 | 2017-08-01 | 中国空间技术研究院 | A kind of method that relative velocity improves gravitational field inversion accuracy between star based on disturbance |
| CN110082787A (en) * | 2019-04-11 | 2019-08-02 | 华东师范大学 | A method of extracting half Sunday of Sunday tide signal from GNSS time sequence |
| CN112278322A (en) * | 2020-10-26 | 2021-01-29 | 中国空间技术研究院 | Marine satellite regression orbit determination method and device and storage medium |
Non-Patent Citations (4)
| Title |
|---|
| Aliasing of ocean tides in satellite gravimetry: a two-step mechanism;Wei Liu，Nico Sneeuw;Journal of Geodesy;1-12 * |
| 卫星编队用于消除海潮模型混频误差影响的可行性研究;赵倩;姜卫平;徐新禹;邹贤才;;中国科学:地球科学(第02期);169-176 * |
| 王长青 ; 许厚泽 ; 钟敏 ; 冉将军 ; 周江存 ; .海潮误差对GRACE时变重力场解算的影响研究.地球物理学报.2015,(第09期),全文. * |
| 赵倩 ; 姜卫平 ; 徐新禹 ; 邹贤才 ; .卫星编队用于消除海潮模型混频误差影响的可行性研究.中国科学:地球科学.2015,(第02期),169-176. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115356776A (en) | 2022-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schwarz et al. | The use of FFT techniques in physical geodesy | |
| Pail et al. | First GOCE gravity field models derived by three different approaches | |
| Nassar et al. | Modeling inertial sensor errors using autoregressive (AR) models | |
| CN110058236A (en) | It is a kind of towards three-dimensional Ground Deformation estimation InSAR and GNSS determine Quan Fangfa | |
| EP3213033B1 (en) | Method of estimating a navigation state constrained in terms of observability | |
| KR20070059105A (en) | Improved gps accumulated delta range processing for navigation processing | |
| CN103743395A (en) | Time delay compensation method in inertia gravity matching combined navigation system | |
| WO2020140378A1 (en) | Method for post-compensation of motion error of rotating accelerometer gravity gradiometer | |
| CN113820003A (en) | Acceleration real-time reconstruction dynamic displacement method suitable for bridge vibration monitoring | |
| CN109471192B (en) | High-precision dynamic data processing method for full-automatic gravity tester | |
| CN115356776B (en) | Sea tide mixing error estimation method in satellite gravity measurement | |
| Nassar et al. | Improving positioning accuracy during kinematic DGPS outage periods using SINS/DGPS integration and SINS data de-noising | |
| Rapp et al. | Role of ocean variability and dynamic ocean topography in the recovery of the mean sea surface and gravity anomalies from satellite altimeter data | |
| Marsh et al. | Seasat altimeter timing bias estimation | |
| CN111504278A (en) | Sea wave detection method based on self-adaptive frequency domain integration | |
| CN103994775A (en) | Inertia measurement unit calibration method for low precision dual-axis rotation-dwell equipment having orientation reference | |
| Ben et al. | Real-time heave motion measurement by adaptive band-pass filter based on strapdown INS | |
| Alldredge | Varying geomagnetic anomalies and secular variation | |
| CN109115229B (en) | Method for measuring high-frequency attitude of spacecraft by using low-frequency attitude measurement sensor | |
| Sansò et al. | Geoid and Sea Surface Topography from satellite and ground data in the Mediterranean region-A review and new proposals | |
| CN110058324B (en) | Strapdown gravimeter horizontal component error correction method using gravity field model | |
| Sun et al. | Novel approach to GPS/SINS integration for IMU alignment | |
| CN107664499B (en) | On-line noise reduction method for accelerometer of ship strapdown inertial navigation system | |
| Zhao et al. | On the inference of supergranular flows by time-distance helioseismology | |
| CN111351483A (en) | Recursive multi-subsample large dynamic inertial navigation 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |