CN103018783A - Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method - Google Patents
Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method Download PDFInfo
- Publication number
- CN103018783A CN103018783A CN2012105779732A CN201210577973A CN103018783A CN 103018783 A CN103018783 A CN 103018783A CN 2012105779732 A CN2012105779732 A CN 2012105779732A CN 201210577973 A CN201210577973 A CN 201210577973A CN 103018783 A CN103018783 A CN 103018783A
- Authority
- CN
- China
- Prior art keywords
- satellite
- earth
- expression
- orbital
- gravity
- 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.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
The invention relates to a gravity satellite formation orbital stability optimization design and an earth gravity field precision inversion method on the basis of disturbing inter-satellite distance principle, in particular to an orbital stability optimization design method for a four-satellite system (FSS). In order to guarantee the stability of the four-satellite system, the quantity of satellite orbits is optimally designed, orbital semi-major axes, orbital eccentricities, orbital inclinations and right ascensions of ascending nodes keep unchanged, the difference of arguments of perigees each pair of satellites and the difference of mean anomalies of the pair of satellites are 180 degrees respectively, an initial argument of perigee of each satellite is arranged at the equator, an initial mean anomaly of each satellite is arranged at a pole, and the ratio of the semi-major axis of each elliptical orbit of the four-satellite system to a semi-minor axis of the elliptical orbit of the four-satellite system is 2:1. The gravity satellite formation orbital stability optimization design and the earth gravity field precision inversion method have the advantages that an earth gravity field is precisely and quickly inverted on the basis of a disturbing inter-satellite distance process; and the orbits are high in stability owing to the method, the earth gravity field computation precision is effectively improved, the gravity field inversion speed is increased to a great extent, and requirements on the performance of a computer are low.
Description
Technical field
The present invention relates to the interleaving techniques fields such as space geodesy, satellite gravimetry, geophysics, cosmonautics, particularly relate to a kind of four stars formation system (FSS:Four-Satellite System) orbital stability Optimization Design, and then utilize the interstellar distance of satellite borne laser interference measuring instrument, the satellite orbital position of GPS receiver and the satellite nonconservative force observation data of accelerometer, based on disturbance interstellar distance method of inversion Effective Raise earth gravity field precision.
Background technology
Be to utilize Satellite Tracking satellite height/low low technical (SST-HL/LL:Satellite-to-Satellite Tracking in the High-Low/Low-Low mode) to promote the New Times to " the wisdom earth " cognitive ability this century.As shown in Figure 1, international Gravity Satellite CHAMP(Challenging Minisatellite Payload), GRACE(Gravity Recovery and Climate Experiment) and GOCE(Gravity Field and Steady-State Ocean Circulation Explorer) in succession succeed in sending up and more the emission that is about to of high precision GRACE Follow-On satellite indicating that the mankind will welcome unprecedented high precision Satellite gravity Detection Techniques epoch.Although the more traditional gravimetry technology of GRACE satellite gravity survey technology (vehicle-mounted, boat-carrying and airborne) can be efficiently, high precision and high spatial resolution are surveyed the earth static (middle long wave) and the time change (long wave) gravity field, but the inherent defect of GRACE Gravity Satellite system (the first, can't reduce satellite orbital altitude; The second, can't improve the load measurement precision; The 3rd, can't obtain vertical gravity gradient; The 4th, can't deduct high-frequency signal and obscure) can't eliminate by self-regulation.Therefore, in order effectively to remedy the weak point of GRACE, and then high precision and high spatial resolution set up earth static state (medium short wave) and the time become (middle long wave) gravity field model, only have the New Satellite of enforcement gravimetry plan just can efficient solution to determine (as shown in Figure 2).
Massonnet was applied to the passive radar interferometry with the four stars Formation Technology first in 1998; Sneeuw and Schaub have proposed the new technology based on four stars formation system precision detection earth gravity field in 2004.As shown in Figure 2, as follows based on the main thought of four stars Formation Technology high-acruracy survey earth gravity field: (1) every single star orbits the earth along elliptical orbit separately; (2) barycenter of four stars formation system is with the earth movements of circular orbit pattern ring; (3) every single star rotates (semi-major axis around formation system barycenter with the elliptical orbit form
:Semi-minor axis=2
:1).Gravitation potential of earth
Right respectively
Second derivative be expressed as follows
Wherein, the gravitation potential of earth second derivative is symmetric tensor, satisfies the Laplace equation simultaneously under vacuum condition and shows as without mark
V Xx +
V Yy +
V Zz =0, therefore in 9 gravity gradient components, there are 5 to be independently.GRACE is the string type Formation Technology of rail double star employing front and back tracking altogether, and being equivalent to base length is interstellar distance
The horizontal gravity gradient instrument, therefore only can measure the horizontal gravity gradient component of direction of visual lines
The four stars Formation Technology can be measured horizontal and vertical gravity gradient component simultaneously
Because the power spectrum of vertical gravity gradient is about 2 times of the horizontal gravity gradient power spectrum
Therefore, owing to not only having increased vertical gravity gradient signal, and reduced largely satellite orbital altitude and improved the combined influence of crucial load measurement precision, improve at least 10 times based on the precision of four stars Formation Technology inverting earth gravity field than the measuring accuracy of GRACE string type Formation Technology.
Because asymmetry and heterogeneity and the compression of the earth item of earth gravity field
J 2Combined influence and effect, the orbital stability of four stars formation system will occur sharply and fast drift.Although the stability that can keep total system by track and the attitude thruster of every satellite self, a large amount of jet fuel consumption will cause the rapidly shortening of four stars formation lifetime of system.Therefore, the optimal design of four stars formation system track radical is to set up the gordian technique of high precision, high spatial resolution and high order earth gravity field model.
Summary of the invention
The objective of the invention is: proposed the orbital stability optimal design of four stars formation system, and further improved the earth gravity field precision based on the disturbance interstellar distance method of inversion.
For achieving the above object, the present invention has adopted following technical scheme:
1) optimal design of orbit elements of satellite
Obtain four Gravity Satellites (FSS-1/2/3/4) formation systematic survey data, every single star edge elliptical orbit separately orbits the earth in described four Gravity Satellites, the barycenter of four stars formation system is with the earth movements of circular orbit pattern ring, and every single star rotates with the elliptical orbit form around formation system barycenter.
Four stars formation system can be by setting 6 Kepler orbit elements (semi-major axis of orbits
a, orbital eccentricity
e, orbit inclination
i, right ascension of ascending node
Ω, argument of perigee
ωWith initial mean anomaly
M) realize, specific rules is as follows: (1) semi-major axis of orbit
a, orbital eccentricity
e, orbit inclination
iAnd right ascension of ascending node
ΩRemain unchanged; (2) argument of perigee of every pair of satellite (FSS-1/2 and FSS-3/4)
ω 1,
ω 2,
ω 3With
ω 4And mean anomaly
M 1,
M 2,
M 3With
M 4Differ respectively 180
o(
ω 1=
ω 2+ 180
o,
ω 3=
ω 4+ 180
o M 1=
M 2+ 180
o,
M 3=
M 4+ 180
o), the initial argument of perigee of every satellite is arranged at the place, equator and initial mean anomaly designs in the limit place; (3) setting the semi-major axis of four stars formation system elliptical orbit and the ratio of semi-minor axis is
ρ Max : ρ Min=2
:1.The initial argument of perigee of every satellite is arranged at the place, equator and but initial mean anomaly is arranged at the drift of limit place establishment four stars formation system, and then guarantees earth gravity field Measurement sensibility and the accuracy of the system of forming into columns.
The measurement data of described four Gravity Satellite formation systems comprises interstellar distance, satellite orbital position and satellite nonconservative force; By described measurement data inverting earth gravity field.
2) the accurate inverting of earth gravity field
Based on Newton's interpolation formula, single star actual measurement orbital position
Taylor expansion be expressed as follows
Wherein,
The expression binomial coefficient,
,
The moment of expression calculation level,
The initial time of expression interpolation point,
The expression sampling interval,
The quantity of expression interpolation point.
By calculating the poor of actual measurement track and reference orbit
, can obtain the Taylor expansion of single star disturbance orbital position
Based on the second derivative of formula (4), single star disturbance orbital acceleration is expressed as follows
Based on formula (5), double star disturbance orbital acceleration differential representation is as follows
Wherein,
With
Represent respectively double star Relative Perturbation orbital position vector and Relative Perturbation orbital acceleration vector,
With
Represent respectively double star disturbance orbital position vector separately,
With
Represent respectively double star disturbance orbital acceleration vector separately.
Wherein,
,
The disturbance interstellar distance observed quantity of expression laser interferometry instrument.
Wherein,
Expression other Relative Perturbation conservative force except terrestrial gravitation mainly comprises: lunisolar attraction, solid tide power, oceanic tide power, atmospheric tide power and utmost point tidal force, and it can calculate acquisition by setting up model;
Expression Relative Perturbation nonconservative force mainly comprises: atmospherical drag, sun optical pressure, terrestrial radiation are pressed and orbit altitude and attitude control, and it can measure acquisition in real time by star accelerometer;
Expression Relative Perturbation center gravitation
Wherein,
The expression earth quality
And universal gravitational constant
Long-pending;
Disturbance the earth's core radius of expression double star,
Expression disturbance orbital position vector
3 components;
Represent relative earth disturbance gravitation,
,
The expression gradient operator,
With
The earth disturbing potential that represents respectively double star
With
Gradient, earth disturbing potential
Be expressed as follows
(12)
Wherein,
Represent respectively the earth's core radius, geocentric colatitude degree and geocentric longitude,
Expression earth mean radius,
LThe maximum order that expression earth disturbing potential is pressed spherical function expansion;
Represent normalized association Legendre function,
lThe expression exponent number,
mThe expression number of times;
With
Represent Geopotential coefficient to be estimated.
By with formula (10) substitution formula (9), disturbance interstellar distance observation equation is expressed as follows
The earth gravity field inversion step of the four stars formation system principle among the present invention is as follows: the FSS-1/2 Double Satellites Observation data of the system of at first four stars being formed into columns and FSS-3/4 Double Satellites Observation data are substitution disturbance interstellar distance observation equation (13) respectively, then unites and finds the solution Geopotential coefficient
With
Because disturbance moonscope amount (actual measured amount-reference quantity) is more responsive to the earth gravity field precision, therefore can improve than before 3~5 times based on disturbance interstellar distance observation equation (13) inverting earth gravity field precision.
The orbital stability optimal design that the present invention is based on four stars formation system is conducive to accurately and the characteristics of fast inversion earth gravity field design, and its advantage is: 1) orbital stability of four stars formation system is better; 2) the earth gravity field computational accuracy is higher; 3) orbital stability design process physical meaning is clear and definite; 4) improve largely the earth gravity field inversion speed; 5) require lower to computing power.Be different from domestic and international prior art achievement, the present invention has carried out the stability Design research of four stars formation system track radical again, adopts the initial argument of perigee with every satellite
ωBe arranged at place, equator and initial mean anomaly
MDesign in north and south limit place, be conducive to keep four stars formation Systems balanth, the drift of establishment four stars formation system; Based on disturbance interstellar distance method, accurately and fast inversion the FSS-1/2/3/4 earth gravity field.
Four, description of drawings
Fig. 1Be the international three phase Gravity Satellites that succeeded in sending up, comprise the single star of CHAMP, NASA (NASA) and the GRACE double star of German space agency (DLR) and the GOCE gravity gradient satellite of European Space Agency (ESA) at German Potsdam earth science research center (GFZ).
Fig. 2Four stars formation systematic survey schematic diagram, wherein satellite orbit semi-major axis
a, orbital eccentricity
e, orbit inclination
iAnd right ascension of ascending node
ΩRemain unchanged the argument of perigee of every pair of satellite
ωAnd mean anomaly
MDiffer respectively 180
o, the initial argument of perigee of every satellite
ωBe arranged at place, equator and initial mean anomaly
MDesign in north and south limit place, and the ratio of the semi-major axis of satellites formation system elliptical orbit and semi-minor axis is 2:1.
Fig. 3It is the optimal design parameter table of four stars formation system track radical.
Fig. 4 aIt is the argument of perigee (north and south limit) of FSS-1/2 satellite.
Fig. 4 bIt is the mean anomaly (equator) of FSS-1/2 satellite.
Fig. 4 cIt is the interstellar distance of FSS-1/2 in the four stars formation system.
Fig. 5 aIt is the argument of perigee (equator) of FSS-3/4 satellite.
Fig. 5 bIt is the mean anomaly (north and south limit) of FSS-3/4 satellite.
Fig. 5 cIt is the interstellar distance of FSS-3/4 in the four stars formation system.
Fig. 6 aIt is the interstellar distance coloured noise of laser interferometry instrument.
Fig. 6 bIt is the orbital position coloured noise of GPS receiver.
Fig. 6 cIt is the nonconservative force coloured noise of accelerometer.
Fig. 7Be based on disturbance interstellar distance method inverting 120 rank FSS-1/2/3/4 accumulative total geoid surface precision.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is further described.
The design of Gravity Satellite Formation Orbit optimizing stability and accurate inverting terrestrial gravitation field method are:
1) optimal design of orbit elements of satellite
Obtain four Gravity Satellites (FSS-1/2/3/4) formation systematic survey data, every single star edge elliptical orbit separately orbits the earth in described four Gravity Satellites, the barycenter of four stars formation system is with the earth movements of circular orbit pattern ring, and every single star rotates with the elliptical orbit form around formation system barycenter.
As shown in Figure 2, four stars formation system can be by setting 6 Kepler orbit elements (semi-major axis of orbits
a, orbital eccentricity
e, orbit inclination
i, right ascension of ascending node
Ω, argument of perigee
ωWith initial mean anomaly
M) realize, specific rules is as follows: (1) semi-major axis of orbit
a, orbital eccentricity
e, orbit inclination
iAnd right ascension of ascending node
ΩRemain unchanged; (2) argument of perigee of every pair of satellite (FSS-1/2 and FSS-3/4)
ω 1,
ω 2,
ω 3With
ω 4And mean anomaly
M 1,
M 2,
M 3With
M 4Differ respectively 180
o(
ω 1=
ω 2+ 180
o,
ω 3=
ω 4+ 180
o M 1=
M 2+ 180
o,
M 3=
M 4+ 180
o), the initial argument of perigee of every satellite is arranged at the place, equator and initial mean anomaly designs in north and south limit place; (3) setting the semi-major axis of four stars formation system elliptical orbit and the ratio of semi-minor axis is
ρ Max : ρ Min=2
:1.The measurement data of described four Gravity Satellite formation systems comprises the interstellar distance of satellite borne laser interference measuring instrument, the satellite orbital position of GPS receiver and the satellite nonconservative force of accelerometer; By described measurement data inverting earth gravity field.
The present invention utilize the linear single-step process of Runge-Kutta in conjunction with 12 rank Adams-Cowell linear multistep method Numerical Integral Formulas analogue simulations the FSS-1/2/3/4 satellite orbit, wherein orbit altitude 250 km, observation time 30 days, sampling interval 10 s, reference to gravitational model E GM2008, Kepler orbit elements as shown in Figure 3.When the satellite argument of perigee
ω(Fig. 4 a) and mean anomaly to be arranged at limit place, north and south
MWhen being arranged at Chi Daochu (Fig. 4 b), four stars formation system drifting will increasingly sharpen; Shown in Fig. 4 c, the interstellar distance of FSS-1/2 increases by 50 km every day at least in the four stars formation system, and then can't realize the stable mutually tracking of FSS-1/2 satellite.The present invention adopts the argument of perigee of every satellite
ω(Fig. 5 a) and mean anomaly to be arranged at Chi Daochu
MBe arranged at limit place, north and south (Fig. 5 b), so four stars formation system is highly stable; Shown in Fig. 5 c, the interstellar distance stable maintenance of FSS-3/4 does not increase and unlimited the increase in time at 50~110 km in the four stars formation system.Contrast as can be known by Fig. 4 c and 5c, the initial argument of perigee of every satellite is arranged at the place, equator and but initial mean anomaly is arranged at the drift of limit place, north and south establishment four stars formation system, and then guarantees earth gravity field Measurement sensibility and the accuracy of the system of forming into columns.
2) the accurate inverting of earth gravity field
Based on Newton's interpolation formula, the Taylor expansion of single star actual measurement orbital position is expressed as follows
(14)
Wherein,
The expression binomial coefficient,
,
The moment of expression calculation level,
The initial time of expression interpolation point,
The expression sampling interval,
The quantity of expression interpolation point.
Based on formula (14)-(15), single star disturbance orbital position
Taylor expansion be expressed as follows
(16)
Based on the second derivative of formula (16), single star disturbance orbital acceleration is expressed as follows
Based on formula (17), double star disturbance orbital acceleration differential representation is as follows
(18)
Wherein,
With
Represent respectively double star Relative Perturbation orbital position vector and Relative Perturbation orbital acceleration vector,
With
Represent respectively double star disturbance orbital position vector separately,
With
Represent respectively double star disturbance orbital acceleration vector separately.
In formula (18),
The sight line representation in components as follows
In formula (19),
Can be rewritten as
Wherein,
,
The disturbance interstellar distance observed quantity of expression laser interferometry instrument.
(22)
Wherein,
Expression other Relative Perturbation conservative force (established model calculating) except terrestrial gravitation mainly comprises: lunisolar attraction, solid tide power, oceanic tide power, atmospheric tide power, utmost point tidal force etc.;
Expression Relative Perturbation nonconservative force (star accelerometer measurement) mainly comprises: atmospherical drag, sun optical pressure, terrestrial radiation pressure, orbit altitude and attitude control etc.;
Expression Relative Perturbation center gravitation
Wherein,
The expression earth quality
And universal gravitational constant
Long-pending;
Disturbance the earth's core radius of expression double star,
Expression disturbance orbital position vector
3 components;
Represent relative earth disturbance gravitation,
,
The expression gradient operator,
With
The earth disturbing potential that represents respectively double star
With
Gradient, earth disturbing potential
Be expressed as follows
(24)
Wherein,
Represent respectively the earth's core radius, geocentric colatitude degree and geocentric longitude,
Expression earth mean radius,
LThe maximum order that expression earth disturbing potential is pressed spherical function expansion;
Represent normalized association Legendre function,
lThe expression exponent number,
mThe expression number of times;
With
Represent Geopotential coefficient to be estimated.
By with formula (22) substitution formula (21), disturbance interstellar distance observation equation is expressed as follows
The earth gravity field inversion step of the four stars formation system principle among the present invention is as follows: at first, and with the four stars respectively substitution disturbance interstellar distance observation equation (25) of the FSS-1/2 Double Satellites Observation data of system and FSS-3/4 Double Satellites Observation data of forming into columns; Then, two prescription journeys are merged into a large-scale overdetermined equation group, namely the equation number is greater than the unknown number number; At last, by least square method, solve Geopotential coefficient
With
It will be understood by those skilled in the art that based on above-mentioned about Geopotential coefficient
With
System of equations, also can adopt other data processing method combined calculation to obtain earth gravity field.Because disturbance moonscope amount (actual measured amount-reference quantity) is more responsive to the earth gravity field precision, therefore can improve than before 3~5 times based on disturbance interstellar distance observation equation (25) inverting earth gravity field precision.
Based on the Gauss-Markov model, the coloured noise of analog satellite observed reading is expressed as follows
Wherein,
The expression related coefficient;
The random white noise of expression normal distribution,
jThe number of expression observation station;
The coloured noise that expression has correlativity.
Fig. 6 a~6c represents respectively based on Gauss-Markov coloured noise model, utilizes interstellar distance coloured noise and orbital position and the axial coloured noise of nonconservative force x of related coefficient (nonconservative force 0.90 of the interstellar distance 0.85 of laser interferometry instrument, the orbital position 0.95 of GPS receiver, star accelerometer) and sampling interval 10 s simulation.By adding coloured noise, make the nonconservative force observation data of the orbital position of interstellar distance, GPS receiver of the four stars formation system laser interferometry instrument of numerical simulation and star accelerometer closer to measured data, improve the checking precision of earth gravity field inverting of the present invention.
As shown in Figure 7, dotted line represents the measured precision of the 120 terrace gravity field model E IGEN-GRACE02S that announce at German Potsdam earth science research center (GFZ), and accumulative total geoid surface precision is 1.894 * 10 at place, 120 rank
-1M; Asterisk line and solid line represent respectively based on double star string type (GRACE-II) and four stars formation system of the present invention (FSS-1/2/3/4), by related coefficient (nonconservative force 0.90 of the interstellar distance 0.85 of laser interferometry instrument, the orbital position 0.95 of GPS receiver and star accelerometer), observation time 30 days and sampling interval 10 s, based on the simulation precision of disturbance interstellar distance method inverting GRACE-II and FSS-1/2/3/4 earth gravity field, accumulative total geoid surface precision is respectively 4.785 * 10 at place, 120 rank
-4M and 1.162 * 10
-4M.Result of study shows: the first, with respect to GRACE-II double star string type formation system because FSS-1/2/3/4 four stars formation of the present invention system has increased the measurement of vertical gravity gradient signal, so Effective Raise the sensing precision of middle long wave earth gravity field.At 120 rank places, improve approximately 4 times based on the precision of FSS-1/2/3/4 four stars formation parity inversion earth gravity field of the present invention than the inversion accuracy of GRACE-II double star string type formation system.The second, improve at least 10 times based on the inversion accuracy of the more current GRACE satellite of precision of FSS-1/2/3/4 four stars formation parity inversion earth gravity field of the present invention.
3) characteristics of four stars formation system
The first, can reduce largely satellite orbital altitude.Because GRACE double star string type formation system adopts star accelerometer to measure in real time nonconservative force at rail, then deducts the nonconservative force effect in Data Post again.Therefore, because the negative interference of nonconservative force, the orbit altitude of GRACE satellite (500 km) can't effectively reduce.Studies show that both at home and abroad: every reduction by 100 km of Gravity Satellite orbit altitude, the nonconservative force (take atmospherical drag as main) that acts on satelloid approximately increases by 10 times.For quiet and stable working environment and then its measuring accuracy of Effective Raise (suitably shorten and measure dynamic range) being provided for the crucial load of Gravity Satellite, and by shielding action in the nonconservative force (atmospherical drag, sun optical pressure, terrestrial radiation pressure, orbit altitude and attitude control etc.) of satelloid and then effectively prolong serviceable life of satellite, nonconservative force bucking-out system (DFCS:Drag-Free Control System) will be carried by four stars formation system.Therefore, four stars formation system can substantively reduce satellite orbital altitude (250 km), and then establishment terrestrial gravitation field signal is with the exponential damping effect of satellite orbital altitude increase.
The second, can substantively improve crucial load measurement precision.GRACE double star string type formation system adopts K wave band measuring instrument to measure speed (1 μ m/s) between interstellar distance (10 μ m) and star, acts on simultaneously the nonconservative force (10 of satelloid based on the SuperSTAR accelerometer measures
-10M/s
2).In order further to reduce crucial load error to the negative effect of earth gravity field precision, four stars formation system will carry more high-precision laser interferometry instrument and accurately measure speed (1 nm/s) between interstellar distance (10 nm) and star, act on simultaneously the nonconservative force (10 of satelloid by nonconservative force bucking-out system high precision balance
-13M/s
2).Because the crucial load measurement precision of four stars formation system is able to significant lifting, therefore can improve at least 10 times based on the more current model accuracy based on GRACE double star string type formation system made of the earth gravity field model precision of four stars formation system made.
The 3rd, can measure simultaneously three-dimensional gravity gradient observed reading.Because GRACE double star string type formation system is equivalent to the horizontal gravity gradient instrument that base length is interstellar distance (220 km), so GRACE can't obtain vertical gravity gradient signal, and then lost largely the earth gravity field precision.Four stars formation system can accurately obtain three-dimensional gravity gradient (level and vertical) observation signal simultaneously, not only can increase substantially precision and the spatial resolution of earth gravity field, and can effectively remove the warp-wise band error (isotropy is more excellent) that terrestrial time that the disappearance owing to vertical gravity gradient signal causes becomes gravity field.
The 4th, but the confusion effect of establishment high-frequency signal.The confusion effect of the high frequency errors such as atmosphere and oceanic tide variation is the most critical factor of loss GRACE earth gravity field inversion accuracy.Because the high frequency errors such as GRACE atmosphere and oceanic tide variation are close to atmosphere and oceanic tide computation model error, thus the high frequency errors such as atmosphere and oceanic tide variation difficult from the GRACE earth gravity field model accurate deduction.Four stars formation system, not only can be similar to the gravity gradient component of GRACE double star string type formation systematic survey direction of visual lines, to increase simultaneously the gravity gradient signal of measuring perpendicular to direction of visual lines, be intended to cut down the high frequency confusion effect to the negative interference of earth gravity field inversion accuracy.
The 5th, become the gravity field signal in the detectable earth during long wave.Become gravity field signal (spatial resolution is 400 km approximately) when GRACE double star string type formation system only can survey earth long wave, become gravity field signal (spatial resolution is better than 200 km) when can't high precision determining in the earth long wave.Four stars formation system becomes the sensing precision of gravity field signal when adopting measuring instrument and nonconservative force bucking-out system between novel and accurate laser interference star will increase substantially long wave in the earth, be intended to for the science such as seismology, thalassography, glaciology, hydrology and national defence crossing research field provide more high precision and more the terrestrial time of high spatial resolution become gravity field information.
Claims (2)
1. one kind is adopted Gravity Satellite formation inverting terrestrial gravitation field method, it is characterized in that:
Obtain four Gravity Satellite FSS-1/2/3/4 formation systematic survey data, every single star of described four Gravity Satellites orbits the earth along elliptical orbit separately, the barycenter of four stars formation system is with the earth movements of circular orbit pattern ring, and every single star rotates with the elliptical orbit form around formation system barycenter;
Wherein, 6 Kepler orbit elements of four Gravity Satellite formation defaults comprise semi-major axis of orbit
a, orbital eccentricity
e, orbit inclination
i, right ascension of ascending node
Ω, argument of perigee
ωWith initial mean anomaly
M, wherein,
The semi-major axis of orbit of the first, four Gravity Satellite
a, orbital eccentricity
e, orbit inclination
iAnd right ascension of ascending node
ΩIn the satellite transit process, remain unchanged;
The second, the argument of perigee of the every couple of satellite FSS-1/2 and FSS-3/4
ω 1,
ω 2,
ω 3With
ω 4And mean anomaly
M 1,
M 2,
M 3With
M 4Differ respectively 180
o,
ω 1=
ω 2+ 180
o,
ω 3=
ω 4+ 180
o,
M 1=
M 2+ 180
o,
M 3=
M 4+ 180
o, the initial argument of perigee of every satellite is arranged at the place, equator and initial mean anomaly designs in the limit place;
The semi-major axis of the three, four Gravity Satellite formation system elliptical orbit
ρ MaxAnd semi-minor axis
ρ MinRatio be
ρ Max : ρ Min=2
:1;
The measurement data of described four Gravity Satellite formation systems comprises the interstellar distance of satellite borne laser interference measuring instrument, the satellite orbital position of GPS receiver and the satellite nonconservative force of accelerometer; By described measurement data inverting earth gravity field.
2. according to the method for the described Gravity Satellite formation of claim 1 inverting earth gravity field, it is characterized in that,
Described method by measurement data inverting earth gravity field is,
Calculate the Taylor expansion of single star actual measurement orbital position based on Newton's interpolation formula
Wherein,
The expression binomial coefficient,
,
The moment of expression calculation level,
The initial time of expression interpolation point,
The expression sampling interval,
The quantity of expression interpolation point;
The Taylor expansion of single star reference orbit position is expressed as follows
Second derivative based on formula (3) can obtain single star disturbance orbital acceleration
Can obtain double star disturbance orbital acceleration difference based on formula (4)
Wherein,
With
Represent respectively double star Relative Perturbation orbital position vector and Relative Perturbation orbital acceleration vector,
With
Represent respectively double star disturbance orbital position vector separately,
With
Represent respectively double star disturbance orbital acceleration vector separately;
In formula (5),
Concrete form be expressed as follows
Wherein,
Expression other Relative Perturbation conservative force except terrestrial gravitation mainly comprises: lunisolar attraction, solid tide power, oceanic tide power, atmospheric tide power and utmost point tidal force;
Expression Relative Perturbation nonconservative force mainly comprises: atmospherical drag, sun optical pressure, terrestrial radiation are pressed and orbit altitude and attitude control;
Expression Relative Perturbation center gravitation
Wherein,
The expression earth quality
And universal gravitational constant
It is long-pending,
Disturbance the earth's core radius of expression double star,
Expression disturbance orbital position vector
3 components;
Represent relative earth disturbance gravitation,
,
The expression gradient operator,
With
The earth disturbing potential that represents respectively double star
With
Gradient, earth disturbing potential
Be expressed as follows
Wherein,
Represent respectively the earth's core radius, geocentric colatitude degree and geocentric longitude,
Expression earth mean radius,
LThe maximum order that expression earth disturbing potential is pressed spherical function expansion;
Represent normalized association Legendre function,
lThe expression exponent number,
mThe expression number of times;
With
Represent Geopotential coefficient to be estimated;
By with formula (6) substitution formula (5), disturbance interstellar distance observation equation is expressed as follows
Wherein,
,
The disturbance interstellar distance observed quantity of expression laser interferometry instrument,
The unit vector of second satellite is pointed in expression by first satellite;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210577973.2A CN103018783B (en) | 2012-12-27 | 2012-12-27 | Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210577973.2A CN103018783B (en) | 2012-12-27 | 2012-12-27 | Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103018783A true CN103018783A (en) | 2013-04-03 |
CN103018783B CN103018783B (en) | 2015-04-08 |
Family
ID=47967595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210577973.2A Expired - Fee Related CN103018783B (en) | 2012-12-27 | 2012-12-27 | Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103018783B (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353767A (en) * | 2015-12-03 | 2016-02-24 | 上海新跃仪表厂 | Satellite formation keeping control method capable of improving payload observation efficiency |
CN106468787A (en) * | 2016-08-31 | 2017-03-01 | 昆明理工大学 | GRAVITY EARTH TIDE signal independence harmonic component extraction method based on PBIL |
CN106840159A (en) * | 2016-12-12 | 2017-06-13 | 西安空间无线电技术研究所 | A kind of coordinate transformation method based on ellipticity angle compensation |
CN108020866A (en) * | 2017-11-20 | 2018-05-11 | 中国空间技术研究院 | A kind of method and system and processor of the inverting of celestial body gravitational field |
CN108762285A (en) * | 2018-05-25 | 2018-11-06 | 北京控制工程研究所 | A kind of targeted attitude collaborative planning method and system of spacecraft multistage composite control |
CN110830103A (en) * | 2019-08-30 | 2020-02-21 | 中国西安卫星测控中心 | Centralized thrust type deployment method for space satellite constellation |
CN111324851A (en) * | 2020-02-19 | 2020-06-23 | 中国科学院遥感与数字地球研究所 | Method for calculating atmospheric boundary observation area of equatorial orbit satellite |
CN111366984A (en) * | 2020-03-23 | 2020-07-03 | 东华理工大学 | Method for determining gravitational field model based on gravity satellite inter-satellite laser ranging system |
CN111439394A (en) * | 2020-03-20 | 2020-07-24 | 上海航天控制技术研究所 | High-precision formation control method based on accelerometer combination |
CN111552003A (en) * | 2020-05-11 | 2020-08-18 | 中国人民解放军军事科学院国防科技创新研究院 | Asteroid gravitational field full-autonomous measurement system and method based on ball satellite formation |
CN113310486A (en) * | 2021-05-31 | 2021-08-27 | 武汉大学 | Universal gravitation positioning method and system |
CN113608244A (en) * | 2021-07-27 | 2021-11-05 | 中国科学院微小卫星创新研究院 | Space gravitational wave detection satellite constellation ground demonstration verification system |
CN115774928A (en) * | 2022-11-20 | 2023-03-10 | 中国人民解放军63921部队 | Improved Laplace model-based initial orbit optimization method for space debris short arc angle measurement only |
CN117270557A (en) * | 2023-09-14 | 2023-12-22 | 中国西安卫星测控中心 | Optimal satellite formation control method for inclination angle and semi-long axis combined bias |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909381A (en) * | 1997-02-19 | 1999-06-01 | Itt Manufacturing Enterprises, Inc. | System of on board prediction of trajectories for autonomous navigation of GPS satellites |
CN102262248B (en) * | 2011-06-03 | 2012-12-26 | 中国科学院测量与地球物理研究所 | Satellite gravity inversion method based on double-satellite spatial three-dimensional interpolation principle |
-
2012
- 2012-12-27 CN CN201210577973.2A patent/CN103018783B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909381A (en) * | 1997-02-19 | 1999-06-01 | Itt Manufacturing Enterprises, Inc. | System of on board prediction of trajectories for autonomous navigation of GPS satellites |
CN102262248B (en) * | 2011-06-03 | 2012-12-26 | 中国科学院测量与地球物理研究所 | Satellite gravity inversion method based on double-satellite spatial three-dimensional interpolation principle |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353767B (en) * | 2015-12-03 | 2019-06-14 | 上海新跃仪表厂 | A kind of satellites formation configuration holding control method improving payload observed efficiency |
CN105353767A (en) * | 2015-12-03 | 2016-02-24 | 上海新跃仪表厂 | Satellite formation keeping control method capable of improving payload observation efficiency |
CN106468787A (en) * | 2016-08-31 | 2017-03-01 | 昆明理工大学 | GRAVITY EARTH TIDE signal independence harmonic component extraction method based on PBIL |
CN106468787B (en) * | 2016-08-31 | 2019-04-09 | 昆明理工大学 | GRAVITY EARTH TIDE signal independence harmonic component extraction method based on PBIL |
CN106840159A (en) * | 2016-12-12 | 2017-06-13 | 西安空间无线电技术研究所 | A kind of coordinate transformation method based on ellipticity angle compensation |
CN106840159B (en) * | 2016-12-12 | 2019-06-18 | 西安空间无线电技术研究所 | A kind of coordinate transformation method based on ellipticity angle compensation |
CN108020866A (en) * | 2017-11-20 | 2018-05-11 | 中国空间技术研究院 | A kind of method and system and processor of the inverting of celestial body gravitational field |
CN108762285A (en) * | 2018-05-25 | 2018-11-06 | 北京控制工程研究所 | A kind of targeted attitude collaborative planning method and system of spacecraft multistage composite control |
CN108762285B (en) * | 2018-05-25 | 2021-08-10 | 北京控制工程研究所 | Target attitude collaborative planning method and system for multi-stage composite control of spacecraft |
CN110830103A (en) * | 2019-08-30 | 2020-02-21 | 中国西安卫星测控中心 | Centralized thrust type deployment method for space satellite constellation |
CN111324851A (en) * | 2020-02-19 | 2020-06-23 | 中国科学院遥感与数字地球研究所 | Method for calculating atmospheric boundary observation area of equatorial orbit satellite |
CN111324851B (en) * | 2020-02-19 | 2023-11-24 | 中国科学院遥感与数字地球研究所 | Method for calculating atmospheric boundary observation area of equatorial orbit satellite |
CN111439394A (en) * | 2020-03-20 | 2020-07-24 | 上海航天控制技术研究所 | High-precision formation control method based on accelerometer combination |
CN111366984A (en) * | 2020-03-23 | 2020-07-03 | 东华理工大学 | Method for determining gravitational field model based on gravity satellite inter-satellite laser ranging system |
CN111366984B (en) * | 2020-03-23 | 2022-10-14 | 东华理工大学 | Method for determining gravitational field model based on gravity satellite inter-satellite laser ranging system |
CN111552003A (en) * | 2020-05-11 | 2020-08-18 | 中国人民解放军军事科学院国防科技创新研究院 | Asteroid gravitational field full-autonomous measurement system and method based on ball satellite formation |
CN113310486A (en) * | 2021-05-31 | 2021-08-27 | 武汉大学 | Universal gravitation positioning method and system |
CN113310486B (en) * | 2021-05-31 | 2023-03-21 | 武汉大学 | Universal gravitation positioning method and system |
CN113608244A (en) * | 2021-07-27 | 2021-11-05 | 中国科学院微小卫星创新研究院 | Space gravitational wave detection satellite constellation ground demonstration verification system |
CN113608244B (en) * | 2021-07-27 | 2023-12-29 | 中国科学院微小卫星创新研究院 | Space gravitational wave detection satellite constellation ground demonstration verification system |
CN115774928A (en) * | 2022-11-20 | 2023-03-10 | 中国人民解放军63921部队 | Improved Laplace model-based initial orbit optimization method for space debris short arc angle measurement only |
CN117270557A (en) * | 2023-09-14 | 2023-12-22 | 中国西安卫星测控中心 | Optimal satellite formation control method for inclination angle and semi-long axis combined bias |
CN117270557B (en) * | 2023-09-14 | 2024-05-03 | 中国西安卫星测控中心 | Optimal satellite formation control method for inclination angle and semi-long axis combined bias |
Also Published As
Publication number | Publication date |
---|---|
CN103018783B (en) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103018783A (en) | Gravity satellite formation orbital stability optimization design and earth gravity field precision inversion method | |
Jekeli | A review of gravity gradiometer survey system data analyses | |
CN102305949B (en) | Method for building global gravitational field model by utilizing inter-satellite distance interpolation | |
CN102313905B (en) | Satellite gravity inversion method based on inter-satellite velocity interpolation principle | |
Lu et al. | Geodesy | |
CN102393535B (en) | Satellite gravity inversion method based on two-star energy interpolation principle | |
CN104374388A (en) | Flight attitude determining method based on polarized light sensor | |
CN102262248A (en) | Satellite gravity inversion method based on double-satellite spatial three-dimensional interpolation principle | |
CN106997061B (en) | A method of gravitational field inversion accuracy is improved based on relative velocity between disturbance star | |
Sun et al. | Autonomous orbit determination via Kalman filtering of gravity gradients | |
CN103076640A (en) | Method for inverting earth gravitational field by using variance-covariance diagonal tensor principle | |
Sneeuw et al. | Science requirements on future missions and simulated mission scenarios | |
CN103645489A (en) | A spacecraft GNSS single antenna attitude determination method | |
CN102519485A (en) | Gyro information-introduced double-position strapdown inertial navigation system initial alignment method | |
Baumann et al. | Absolute airborne gravimetry: a feasibility study | |
CN103076639B (en) | Method for inverting earth gravity field of residual inter-star velocity | |
CN105487405B (en) | Low tracking Gravisat semi-physical system | |
CN103091722A (en) | Satellite gravity inversion method based on load error analysis theory | |
CN103852812A (en) | Method for forward simulation of occultation | |
CN103091721B (en) | Satellite joint inversion earth gravitational field method using different orbit inclination angles | |
CN103093101A (en) | Satellite gravity inversion method base on gravity gradient error model principle | |
Qin et al. | Relativistic tidal effects on clock-comparison experiments | |
CN103901496A (en) | Gravity measuring method based on fiber-optic gyroscope SINS and Big Dipper | |
CN103091723A (en) | Method of reducing influences of gravity satellite centroid adjustment errors to earth gravitational field accuracy | |
US20240142660A1 (en) | Method and system for indirect measurement of gravity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150408 Termination date: 20151227 |
|
EXPY | Termination of patent right or utility model |