CN108919371B - A kind of airborne gravity data downward continuation method and system for combining ground gravity station - Google Patents
A kind of airborne gravity data downward continuation method and system for combining ground gravity station Download PDFInfo
- Publication number
- CN108919371B CN108919371B CN201810816288.8A CN201810816288A CN108919371B CN 108919371 B CN108919371 B CN 108919371B CN 201810816288 A CN201810816288 A CN 201810816288A CN 108919371 B CN108919371 B CN 108919371B
- Authority
- CN
- China
- Prior art keywords
- data
- gravity
- ground
- airborne
- aerial
- 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
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
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
Abstract
The invention discloses a kind of airborne gravity data downward continuation method and system for combining ground gravity station.The described method includes: obtaining in Airborne Gravimetry System ground gravity data in aerial gravimetric data and ground gravity survey system；Determine that Airborne gravimetry influences data according to the ground gravity data；Aerial gravimetric data identical with Airborne gravimetry influence data in the aerial gravimetric data is filtered out, determines the aerial gravimetric data after filtering out；Point mass observation model is constructed according to the aerial gravimetric data after described filter out；The ground gravity data of ground white space are calculated according to the point mass observation model；The ground gravity data of the ground white space are assessed, determine the ground gravity data of non-geodetic face gravity station.It can be improved the ground gravity data assessment precision of non-geodetic face gravity station using method and system provided by the present invention, reduce amount of calculation.
Description
Technical field
The present invention relates to airborne gravity measurement fields, more particularly to a kind of airborne gravity data for combining ground gravity station
Downward continuation method and system.
Background technique
Airborne gravity measurement technology, because it can be in desert, marsh, glacier, hills, virgin forest, land-sea interface etc. one
The region for being difficult to carry out ground gravity survey a bit carries out operation, obtains to fast, economical in high-accuracy high-resolution earth gravitational field
High-frequency information obtains one of most effective technological means of earth gravitational field measurement data to become.
Airborne gravity measurement data is in use, generally require in continuation to earth surface or geoid, main use
It is big in the fusion of different type gravity survey data, the building in the whole world or region earth gravity field model, the whole world or region (seemingly)
Gravity datum map generalization etc. in the refining of ground-level, subaqueous gravity matching assisting navigation；Made with the gravimetric station on ground
For control, downward continuation is carried out to airborne gravity measurement data, the unstability of downward continuation can be weakened, inhibits boundary effect
Influence to continuation result improves continuation result precision；The airborne gravity data of currently used joint ground gravimetric station
The method of downward continuation is mainly least squares collocation.
The basic principle of least squares collocation is to utilize a small amount of high quality ground gravimetric station United Airlines gravity
Data estimate non-geodetic face gravity station, to realize the purpose of airborne gravity measurement data continuation to ground；Advantage is phase
For continuation more direct than single airborne gravity measurement data, continuation result precision can be improved, and continuation result can be provided
Accuracy extimate information；The disadvantage is that ground gravity survey point and the covariance function of airborne gravity data building difficulty are larger, and
And building effect will cause to directly affect to the precision of final continuation result (the ground gravity data of non-geodetic face gravity station)；By
This is it is found that is currently used combines the least squares collocation of the airborne gravity data downward continuation of ground gravimetric station
There are continuation result precision is low, the problem of heavy workload.
Summary of the invention
The object of the present invention is to provide it is a kind of combine ground gravity station airborne gravity data downward continuation method and system,
To solve the ground gravity number of the least squares collocation assessment non-geodetic face gravity station of existing airborne gravity data downward continuation
According to when Evaluation accuracy it is low, the problem of heavy workload.
To achieve the above object, the present invention provides following schemes:
A kind of airborne gravity data downward continuation method for combining ground gravity station, the airborne gravity data downward continuation
Method is based on Airborne Gravimetry System and ground gravity survey system, comprising:
Obtain in the Airborne Gravimetry System ground in aerial gravimetric data and the ground gravity survey system
Gravimetric data；The aerial gravimetric data is the gravity anomaly data that gravimeter acquires on course line, the ground gravity number
According to the gravity anomaly data obtained for ground gravity instrument measurement；
Determine that Airborne gravimetry influences data according to the ground gravity data；The aerial gravimetric data includes described
Airborne gravimetry influences data；It is to be calculated using ground gravity survey data that the Airborne gravimetry, which influences data,
The influence value to airborne gravity measurement data；
Aerial gravimetric data identical with Airborne gravimetry influence data in the aerial gravimetric data is filtered out, really
Surely the aerial gravimetric data after filtering out；
Point mass observation model is constructed according to the aerial gravimetric data after described filter out；
The ground gravity data of ground white space are calculated according to the point mass observation model；
The ground gravity data of the ground white space are assessed, determine the ground gravity of non-geodetic face gravity station
Data.
Optionally, described to determine that Airborne gravimetry influences data according to the ground gravity data, it specifically includes:
Using Poisson Poisson integral equation, according to formulaIt determines empty
Middle gravity measurement influences data；Wherein, Δ glandFor the ground gravity data of ground gravity point；Δg0For Airborne gravimetry shadow
Ring data；Ω is ground Mean height plane；R1It is the earth's core to diameter；(B', L') is the longitude and latitude of ground gravity point；R is aerial point
The earth's core to diameter；L is space length of the aerial measuring point to Ω upper integral face element.
Optionally, described to filter out in the aerial gravimetric data that data are influenced with the Airborne gravimetry is identical in the air
Gravimetric data determines the aerial gravimetric data after filtering out, specifically includes:
According to formulaDetermine the aerial gravimetric data after filtering out；Wherein, L is the aerial gravity after filtering out
Data；For i-th of gravity anomaly observed quantity, i
For i-th of gravity anomaly observed quantityGeocentric distance,For the radius of shellfish Ahab's horse Bjerhammar ball, Dj
For j-th of particle MjBuried depth, ρijFor i-th of observed quantity between j-th of particle at a distance from, ψijFor i-th of observed quantity with
Centre of sphere angle between j-th of particle, n are the number of observed quantity, and m is the number of particle.
Optionally, it is described filtered out according to after aerial gravimetric data building point mass observation model, specifically include:
Airborne gravity data according to the aerial gravimetric data building joint ground gravimetric station after described filter out are downward
The point mass observing matrix of continuation；
Point mass observation model is constructed according to the point mass observing matrix.
Optionally, before the building point mass observation model according to the point mass observing matrix, further includes:
Obtain current coverage of survey area；
Judge whether the current coverage of survey area is greater than coverage of survey area threshold value, obtains the first judging result；
If first judging result is expressed as current coverage of survey area greater than coverage of survey area threshold value, using Regularization Technique solution
It calculates observational equation and constructs point mass observation model；
If first judging result is expressed as current coverage of survey area no more than coverage of survey area threshold value, according to the point mass
Observing matrix constructs point mass observation model.
A kind of airborne gravity data downward continuation system for combining ground gravity station, comprising:
Aerial gravimetric data and ground gravity data acquisition module, it is empty in the Airborne Gravimetry System for obtaining
Ground gravity data in middle gravimetric data and the ground gravity survey system；The aerial gravimetric data is gravity on course line
The gravity anomaly data that instrument acquires, the ground gravity data are the gravity anomaly data that the measurement of ground gravity instrument obtains；
Airborne gravimetry influences data determining module, for determining Airborne gravimetry according to the ground gravity data
Influence data；The aerial gravimetric data includes that the Airborne gravimetry influences data；The Airborne gravimetry influences number
According to the influence value to airborne gravity measurement data to be calculated using ground gravity survey data；
Module is filtered out, for filtering out sky identical with Airborne gravimetry influence data in the aerial gravimetric data
Middle gravimetric data determines the aerial gravimetric data after filtering out；
Point mass observation model constructs module, for the aerial gravimetric data building point mass observation after filtering out according to
Model；
Ground gravity data computation module, for calculating the ground of ground white space according to the point mass observation model
Gravimetric data；
Evaluation module is assessed for the ground gravity data to the ground white space, determines non-geodetic face weight
The ground gravity data in force.
Optionally, the Airborne gravimetry influences data determining module, specifically includes:
Airborne gravimetry influences data determination unit, for utilizing Poisson Poisson integral equation, according to formulaDetermine that Airborne gravimetry influences data；Wherein, Δ glandFor ground gravity
The ground gravity data of point；Δg0Data are influenced for Airborne gravimetry；Ω is ground Mean height plane；R1It is the earth's core to diameter；
(B', L') is the longitude and latitude of ground gravity point；R is the earth's core for putting in the air to diameter；L is sky of the aerial measuring point to Ω upper integral face element
Between distance.
Optionally, described to filter out module, it specifically includes:
Unit is filtered out, for according to formulaDetermine the aerial gravimetric data after filtering out；Wherein, L is filter
Aerial gravimetric data after removing；For i-th of gravity anomaly observed quantity,riFor i-th of gravity anomaly observed quantityThe earth's core away from
From,For the radius of shellfish Ahab's horse Bjerhammar ball, DjFor j-th of particle MjBuried depth, ρijIt is i-th
Observed quantity between j-th of particle at a distance from, ψijFor the centre of sphere angle between i-th of observed quantity and j-th of particle, n is of observed quantity
Number, m are the number of particle.
Optionally, the point mass observation model constructs module, specifically includes:
Point mass observing matrix determination unit, for the aerial gravimetric data building joint ground weight after being filtered out according to
The point mass observing matrix of the airborne gravity data downward continuation of power measurement point；
Point mass observation model construction unit, for constructing point mass observation model according to the point mass observing matrix.
Optionally, before the point mass observation model construction unit, further includes:
Current coverage of survey area acquiring unit, for obtaining current coverage of survey area；
First judging unit obtains first and sentences for judging whether the current coverage of survey area is greater than coverage of survey area threshold value
Disconnected result；
First point mass observation model construction unit, if it is big to be expressed as current coverage of survey area for first judging result
In coverage of survey area threshold value, observational equation is resolved using Regularization Technique and constructs point mass observation model；
Second point quality observes model construction unit, if being expressed as current coverage of survey area not for first judging result
Greater than coverage of survey area threshold value, point mass observation model is constructed according to the point mass observing matrix.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: the invention proposes one kind
The airborne gravity data downward continuation method and system of joint ground gravity station, using the aviation weight of joint ground gravimetric station
The point mass method of force data downward continuation establishes point mass observation model, to calculate the ground gravity number of ground white space
According to, and then the ground gravity data are assessed, the ground gravity data of non-geodetic face gravity station are determined, due to not being related to
The building of covariance function can play an effective role in the control to continuation result, inhibit boundary effect to the shadow of continuation result
It rings, improves the effect of downward continuation, to improve the ground gravity data assessment precision of non-geodetic face gravity station, reduce meter
Operator workload.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is the airborne gravity data downward continuation system composition of joint ground gravimetric station provided by the present invention
Figure；
Fig. 2 is the airborne gravity data downward continuation method flow diagram of joint ground gravity station provided by the present invention；
Fig. 3 is provided by the present invention using Regularization Technique resolving observational equation building point mass observation model process
Figure；
Fig. 4 is the airborne gravity data downward continuation system construction drawing of joint ground gravity station provided by the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide it is a kind of combine ground gravity station airborne gravity data downward continuation method and system,
It can be improved the ground gravity data assessment precision of non-geodetic face gravity station, reduce amount of calculation.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is the airborne gravity data downward continuation system composition of joint ground gravimetric station provided by the present invention
Figure, as shown in Figure 1, Airborne Gravimetry System is by gravity sensor subsystem, worldwide navigation positioning system (Global
Navigation Satellite System, GNSS) positioning subsystem, posture subsystem, laser-measured height subsystem, data adopt
Collect the composition such as subsystem；The aerial grid point data that the present invention utilizes Airborne Gravimetry System to obtain, and then using jointly
The point mass method of the airborne gravity data downward continuation of face gravimetric station, by aerial grid point data downward continuation to earth table
On face or geoid, and the data are made to the application of next step.
Fig. 2 is the airborne gravity data downward continuation method flow diagram of joint ground gravity station provided by the present invention, such as
Shown in Fig. 2, the airborne gravity data downward continuation system based on above-mentioned joint ground gravimetric station, the invention proposes one kind
The airborne gravity data downward continuation method of joint ground gravity station, the airborne gravity data downward continuation method are based on aviation
Gravity measuring system and ground gravity survey system, comprising:
Step 201: obtaining aerial gravimetric data and the ground gravity survey system in the Airborne Gravimetry System
Ground gravity data in uniting；The aerial gravimetric data is the gravity anomaly data that gravimeter acquires on course line, describedly
Face gravimetric data is the gravity anomaly data that the measurement of ground gravity instrument obtains.
Step 202: determining that Airborne gravimetry influences data according to the ground gravity data；The aerial gravimetric data
Data are influenced including the Airborne gravimetry；The Airborne gravimetry influence data be using ground gravity survey data,
The influence value to airborne gravity measurement data being calculated.
According to imperial lattice theorem, outside of the earth disturbing potential can be with including that harmonic function outside one spherical surface of earth interior is forced
Closely, which is commonly referred to as shellfish Ahab's horse Bjerhammar ball.So-called virtual Point-mass Model, exactly to be distributed in multiple shellfish Ahabs
Virtual disturbing point masses system on horse Bjerhammar spherical surface equally to replace the practical gravitation field source of the earth, only requires these disturbances
Position that particle generates at the earth's surface and its derived quantity approach field member observed quantity with permissible accuracy.It is characterized with disturbing point masses external
The formula of gravity anomaly is
In formula (1), riFor i-th of gravity anomaly observed quantityGeocentric distance,For shellfish Ahab horse
The radius of Bjerhammar ball, DjFor j-th of particle MjBuried depth, ρij、ψijRespectively i-th of observed quantity and j-th of matter
Distance and centre of sphere angle between point, n, m are respectively the number of observed quantity and particle, and earth mean radius R, which is not used, here usually makes
6371km, and use the average the earth's core in required region to diameter:
In formula (2), HiIt is the ground elevation of observation point, RiFor corresponded on reference ellipsoid i point to diameter.
According to Poisson Poisson integral equation, by ground gravity point Δ glandCalculate aerial gravimetric data Δ g0Formula be
In formula (3), Ω is ground Mean height plane, and the earth's core is to diameter R1With indicating, (B', L') indicates ground gravity point
Longitude and latitude, r are the earth's core for putting in the air to diameter, and l is space length of the aerial measuring point to Ω upper integral face element:
Step 203: filtering out aerial weight identical with Airborne gravimetry influence data in the aerial gravimetric data
Force data determines the aerial gravimetric data after filtering out.
Ground gravity survey point is deducted from airborne gravity measurement data on the influential sky of Airborne gravimetry data tool
Middle gravity measurement influences data, i.e.,Then Point-mass Model is constructed with the data.
Step 204: point mass observation model is constructed according to the aerial gravimetric data after described filter out.
Write formula (1) as matrix form:
L=AM (5)
Wherein, L is treated airborne gravity dataThe vector of composition, M are the vector of point mass composition, A
For factor arrays.As Δ ghObservation number n it is equal with point mass number m when, formula (5) has only solution:
M=A-1L (6)
As n > m, formula (6) has least square solution:
M=(ATPA)-1ATPL (7)
P indicates Δ g hereinhWeight of observation battle array.
Step 205: the ground gravity data of ground white space are calculated according to the point mass observation model.
Step 206: the ground gravity data of the ground white space being assessed, determine non-geodetic face gravity station
Ground gravity data.
It is arranged on the internal ball of opposite ground survey blank area certain depth using airborne gravity data after processing corresponding
Point mass, as shown in formula (1).After acquiring point mass using formula (5), it can be calculated accordingly referring to formula (1)
The gravimetric data of face blank area, to realize the purpose of Downward Continuation of Airborne Gravity Data to ground.
Before the building point mass observation model according to the point mass observing matrix, further includes:
Obtain current coverage of survey area；
Judge whether current coverage of survey area is greater than coverage of survey area threshold value, if so, resolving observational equation using Regularization Technique
Construct point mass observation model；If it is not, constructing point mass observation model according to the point mass observing matrix.
Wherein, Fig. 3 is provided by the present invention using Regularization Technique resolving observational equation building point mass observation model
Flow chart, as shown in figure 3, establishing regularization according to the observational equation of treated airborne gravity data calculation Point-mass Model
Method observational equation；Spectral factorization is carried out to method of regularization observational equation, Spectral radius is obtained and decomposes equation, calculating speed can be improved
Degree；According to the methods of L-curve method, U curve method, the optimal regularization parameter of regularization equation is determined；By optimal regularization parameter
Method of regularization observational equation is substituted into, matrix calculation result is obtained, constructs point mass observation model.
Fig. 4 is the airborne gravity data downward continuation system construction drawing of joint ground gravity station provided by the present invention, such as
Shown in Fig. 4, a kind of airborne gravity data downward continuation system for combining ground gravity station, comprising:
Aerial gravimetric data and ground gravity data acquisition module 401, for obtaining the Airborne Gravimetry System
Ground gravity data in interior aerial gravimetric data and the ground gravity survey system；The aerial gravimetric data is on course line
The gravity anomaly data that gravimeter acquires, the ground gravity data are the gravity anomaly number that the measurement of ground gravity instrument obtains
According to.
Airborne gravimetry influences data determining module 402, for determining aerial gravity according to the ground gravity data
Measurement influences data；The aerial gravimetric data includes that the Airborne gravimetry influences data；The Airborne gravimetry shadow
Ringing data is the influence value to airborne gravity measurement data that is calculated using ground gravity survey data.
The Airborne gravimetry influences data determining module 402 and specifically includes:
Airborne gravimetry influences data determination unit, for utilizing Poisson Poisson integral equation, according to formulaDetermine that Airborne gravimetry influences data；Wherein, Δ glandFor ground gravity
The ground gravity data of point；Δg0Data are influenced for Airborne gravimetry；Ω is ground Mean height plane；R1It is the earth's core to diameter；
(B', L') is the longitude and latitude of ground gravity point；R is the earth's core for putting in the air to diameter；L is sky of the aerial measuring point to Ω upper integral face element
Between distance.
Module 403 is filtered out, it is identical as Airborne gravimetry influence data in the aerial gravimetric data for filtering out
Aerial gravimetric data, determine filter out after aerial gravimetric data.
The module 403 that filters out specifically includes:
Unit is filtered out, for according to formulaDetermine the aerial gravimetric data after filtering out；Wherein, L is filter
Aerial gravimetric data after removing；For i-th of gravity anomaly observed quantity,riFor i-th of gravity anomaly observed quantityThe earth's core away from
From,For the radius of shellfish Ahab's horse Bjerhammar ball, DjFor j-th of particle MjBuried depth, ρijIt is i-th
Observed quantity between j-th of particle at a distance from, ψijFor the centre of sphere angle between i-th of observed quantity and j-th of particle, n is of observed quantity
Number, m are the number of particle.
Point mass observation model constructs module 404, constructs point mass for the aerial gravimetric data after filtering out according to
Observation model.
The point mass observation model building module 404 specifically includes:
Point mass observing matrix determination unit, for the aerial gravimetric data building joint ground weight after being filtered out according to
The point mass observing matrix of the airborne gravity data downward continuation of power measurement point；
Point mass observation model construction unit, for constructing point mass observation model according to the point mass observing matrix.
Before the point mass observation model construction unit, further includes:
Current coverage of survey area acquiring unit, for obtaining current coverage of survey area；
First judging unit obtains first and sentences for judging whether the current coverage of survey area is greater than coverage of survey area threshold value
Disconnected result；
First point mass observation model construction unit, if it is big to be expressed as current coverage of survey area for first judging result
In coverage of survey area threshold value, observational equation is resolved using Regularization Technique and constructs point mass observation model；
Second point quality observes model construction unit, if being expressed as current coverage of survey area not for first judging result
Greater than coverage of survey area threshold value, point mass observation model is constructed according to the point mass observing matrix.
Ground gravity data computation module 405, for calculating ground white space according to the point mass observation model
Ground gravity data.
Evaluation module 406 is assessed for the ground gravity data to the ground white space, determines non-geodetic face
The ground gravity data of gravity station.
The present invention provides a kind of airborne gravity data downward continuation method and system for combining ground gravity station, it may be assumed that joint
The point mass method of the airborne gravity data downward continuation of ground gravity survey point, is calculated using Poisson Poisson integral equation
It to ground gravity survey point to the influence value of airborne gravity measurement data, and is deducted, is connect from airborne gravity measurement data
Using treated airborne gravity data calculation with respect to point mass on the internal ball of ground gravity survey blank area certain depth
Model is finally calculated the gravimetric data of corresponding ground blank area by Point-mass Model, to realize ground gravity survey point
The purpose of Combined Treatment is carried out with airborne gravity data.The present invention can get the airborne gravity on earth surface or geoid
The high quality continuation of measurement data is as a result, to fusion, the whole world or the region terrestrial gravitation for different type gravity survey data
The building of field model, the refining of the whole world or region (seemingly) geoid, gravity datum figure in subaqueous gravity matching assisting navigation
The application such as generation reliable basic data support is provided.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other
The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment
For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part
It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention；At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.
Claims (10)
1. a kind of airborne gravity data downward continuation method for combining ground gravity station, which is characterized in that the airborne gravity number
Airborne Gravimetry System and ground gravity survey system are based on according to downward continuation method, comprising:
Obtain in the Airborne Gravimetry System ground gravity in aerial gravimetric data and the ground gravity survey system
Data；The aerial gravimetric data is the gravity anomaly data that gravimeter acquires on course line, and the ground gravity data are
The gravity anomaly data that ground gravity instrument measurement obtains；
Determine that Airborne gravimetry influences data according to the ground gravity data；The aerial gravimetric data includes described aerial
Gravity measurement influences data；The Airborne gravimetry influences data for pair using ground gravity survey data, being calculated
The influence value of airborne gravity measurement data；
Aerial gravimetric data identical with Airborne gravimetry influence data in the aerial gravimetric data is filtered out, determines filter
Aerial gravimetric data after removing；
Point mass observation model is constructed according to the aerial gravimetric data after described filter out；
The ground gravity data of ground white space are calculated according to the point mass observation model；
The ground gravity data of the ground white space are assessed, determine the ground gravity number of non-geodetic face gravity station
According to.
2. airborne gravity data downward continuation method according to claim 1, which is characterized in that described according to ground weight
Force data determines that Airborne gravimetry influences data, specifically includes:
Using Poisson Poisson integral equation, according to formulaDetermine aerial weight
Power measurement influences data；Wherein, Δ glandFor the ground gravity data of ground gravity point；Δg0Number is influenced for Airborne gravimetry
According to；Ω is ground Mean height plane；R1It is the earth's core to diameter；(B', L') is the longitude and latitude of ground gravity point；R is the ground put in the air
Heart Xiang Jing；L is space length of the aerial measuring point to Ω upper integral face element.
3. airborne gravity data downward continuation method according to claim 2, which is characterized in that described to filter out the aerial weight
Aerial gravimetric data identical with Airborne gravimetry influence data, determines the aerial gravity number after filtering out in force data
According to specifically including:
According to formulaDetermine the aerial gravimetric data after filtering out；Wherein, L is the aerial gravity number after filtering out
According to；For i-th of gravity anomaly observed quantity,riFor
I-th of gravity anomaly observed quantityGeocentric distance,For the radius of shellfish Ahab's horse Bjerhammar ball, R is
Earth mean radius, DjFor j-th of particle MjBuried depth, ρijFor i-th of observed quantity between j-th of particle at a distance from, ψij
For the centre of sphere angle between i-th of observed quantity and j-th of particle, n is the number of observed quantity, and m is the number of particle.
4. airborne gravity data downward continuation method according to claim 3, which is characterized in that it is described filtered out according to after
Aerial gravimetric data construct point mass observation model, specifically include:
According to the airborne gravity data downward continuation of the aerial gravimetric data building joint ground gravimetric station after described filter out
Point mass observing matrix；
Point mass observation model is constructed according to the point mass observing matrix.
5. airborne gravity data downward continuation method according to claim 4, which is characterized in that described according to the point mass
Observing matrix constructs before point mass observation model, further includes:
Obtain current coverage of survey area；
Judge whether the current coverage of survey area is greater than coverage of survey area threshold value, obtains the first judging result；
If first judging result is expressed as current coverage of survey area greater than coverage of survey area threshold value, is resolved and seen using Regularization Technique
It surveys equation and constructs point mass observation model；
If first judging result is expressed as current coverage of survey area no more than coverage of survey area threshold value, observed according to the point mass
Matrix constructs point mass observation model.
6. a kind of airborne gravity data downward continuation system for combining ground gravity station characterized by comprising
Aerial gravimetric data and ground gravity data acquisition module, weigh in the air for obtaining in the Airborne Gravimetry System
Ground gravity data in force data and the ground gravity survey system；The aerial gravimetric data is that gravimeter obtains on course line
The gravity anomaly data obtained, the ground gravity data are the gravity anomaly data that the measurement of ground gravity instrument obtains；
Airborne gravimetry influences data determining module, for determining that Airborne gravimetry influences according to the ground gravity data
Data；The aerial gravimetric data includes that the Airborne gravimetry influences data；The Airborne gravimetry influences data
Using ground gravity survey data, the influence value to airborne gravity measurement data that is calculated；
Module is filtered out, for filtering out aerial weight identical with Airborne gravimetry influence data in the aerial gravimetric data
Force data determines the aerial gravimetric data after filtering out；
Point mass observation model constructs module, observes mould for the aerial gravimetric data building point mass after filtering out according to
Type；
Ground gravity data computation module, for calculating the ground gravity of ground white space according to the point mass observation model
Data；
Evaluation module is assessed for the ground gravity data to the ground white space, determines non-geodetic face gravity station
Ground gravity data.
7. airborne gravity data downward continuation system according to claim 6, which is characterized in that the Airborne gravimetry shadow
Data determining module is rung, is specifically included:
Airborne gravimetry influences data determination unit, for utilizing Poisson Poisson integral equation, according to formulaDetermine that Airborne gravimetry influences data；Wherein, Δ glandFor ground gravity
The ground gravity data of point；Δg0Data are influenced for Airborne gravimetry；Ω is ground Mean height plane；R1It is the earth's core to diameter；
(B', L') is the longitude and latitude of ground gravity point；R is the earth's core for putting in the air to diameter；L is sky of the aerial measuring point to Ω upper integral face element
Between distance.
8. airborne gravity data downward continuation system according to claim 7, which is characterized in that it is described to filter out module, specifically
Include:
Unit is filtered out, for according to formulaDetermine the aerial gravimetric data after filtering out；Wherein, L is after filtering out
Aerial gravimetric data；For i-th of gravity anomaly observed quantity,riFor i-th of gravity anomaly observed quantityThe earth's core
Distance,For the radius of shellfish Ahab's horse Bjerhammar ball, R is earth mean radius, DjFor j-th of particle Mj's
Buried depth, ρijFor i-th of observed quantity between j-th of particle at a distance from, ψijFor the ball between i-th of observed quantity and j-th of particle
Heart angle, n are the number of observed quantity, and m is the number of particle.
9. airborne gravity data downward continuation system according to claim 8, which is characterized in that the point mass observation model
Module is constructed, is specifically included:
Point mass observing matrix determination unit is surveyed for the aerial gravimetric data building joint ground gravity after filtering out according to
Measure the point mass observing matrix of the airborne gravity data downward continuation of point；
Point mass observation model construction unit, for constructing point mass observation model according to the point mass observing matrix.
10. airborne gravity data downward continuation system according to claim 9, which is characterized in that the point mass observes mould
Before type construction unit, further includes:
Current coverage of survey area acquiring unit, for obtaining current coverage of survey area；
First judging unit obtains the first judgement knot for judging whether the current coverage of survey area is greater than coverage of survey area threshold value
Fruit；
First point mass observation model construction unit is surveyed if being expressed as current coverage of survey area for first judging result and being greater than
Area's range threshold resolves observational equation using Regularization Technique and constructs point mass observation model；
Second point quality observes model construction unit, is not more than if being expressed as current coverage of survey area for first judging result
Coverage of survey area threshold value constructs point mass observation model according to the point mass observing matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810816288.8A CN108919371B (en) | 2018-07-24 | 2018-07-24 | A kind of airborne gravity data downward continuation method and system for combining ground gravity station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810816288.8A CN108919371B (en) | 2018-07-24 | 2018-07-24 | A kind of airborne gravity data downward continuation method and system for combining ground gravity station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108919371A CN108919371A (en) | 2018-11-30 |
| CN108919371B true CN108919371B (en) | 2019-10-08 |
Family
ID=64417239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810816288.8A Active CN108919371B (en) | 2018-07-24 | 2018-07-24 | A kind of airborne gravity data downward continuation method and system for combining ground gravity station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108919371B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113466959A (en) * | 2021-07-02 | 2021-10-01 | 中国地震局地球物理研究所 | Local gravity field modeling method and system based on ground gravity measurement data |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7069147B2 (en) * | 2004-05-28 | 2006-06-27 | Honeywell International Inc. | Airborne based monitoring |
| CN103399350B (en) * | 2013-07-29 | 2016-02-24 | 中国人民解放军国防科学技术大学 | A kind of airborne gravity downward continuation method based on integral iteration algorithm |
| CN103869376A (en) * | 2014-03-20 | 2014-06-18 | 中国石油大学(华东) | Three-dimensional gravity potential field regularized downward-extension method based on depth change and application thereof |
| US20160131794A1 (en) * | 2014-11-12 | 2016-05-12 | Cgg Services Sa | Systems and methods for a gravity survey using a free-fall gravity sensor |
| CN106125149B (en) * | 2016-06-18 | 2017-12-15 | 中国人民解放军信息工程大学 | The optimal buried depth of Point-mass Model middle-shallow layer high-resolution point mass determines method |
-
2018
- 2018-07-24 CN CN201810816288.8A patent/CN108919371B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108919371A (en) | 2018-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kong et al. | Development and application of UAV-SfM photogrammetry for quantitative characterization of rock mass discontinuities | |
| Mohammed et al. | Positional accuracy testing of Google Earth | |
| CN104931022B (en) | Satellite image stereoblock adjustment method based on spaceborne laser altimeter system data | |
| CN103674030B (en) | The deviation of plumb line dynamic measurement device kept based on astronomical attitude reference and method | |
| Cai et al. | The first airborne scalar gravimetry system based on SINS/DGPS in China | |
| CN110017812B (en) | Measuring device and method for vertical line deviation, gravitational acceleration and gravity gradient | |
| CN109283591A (en) | Using ground point as the airborne gravity data downward continuation method and system of control | |
| CN108594319A (en) | A kind of Downward Continuation of Airborne Gravity Data method and system | |
| Goetz et al. | Quantifying uncertainties in snow depth mapping from structure from motion photogrammetry in an alpine area | |
| Marek et al. | Monitoring of the shallow landslide using UAV photogrammetry and geodetic measurements | |
| Salvini et al. | Stability analysis of “Grotta delle Felci” Cliff (Capri Island, Italy): structural, engineering–geological, photogrammetric surveys and laser scanning | |
| CN104567802B (en) | The survey line formula land-sea height transfer method of integrated boat-carrying gravity and GNSS | |
| CN108919371B (en) | A kind of airborne gravity data downward continuation method and system for combining ground gravity station | |
| CN113900069A (en) | Vertical deviation calculation method and system based on interference imaging altimeter | |
| Zhang et al. | Surface topography around the summit of Dome A, Antarctica, from real-time kinematic GPS | |
| Schuler et al. | New data, new techniques and new challenges for updating the state of Svalbard glaciers (SvalGlac) | |
| Kosmatin Fras et al. | Assessment of the quality of digital terrain model produced from unmanned aerial system imagery | |
| Sildvee | Gravity measurements of Estonia | |
| Kosarev et al. | Determining deflections of the vertical in the Western Siberia region: The results of comparison | |
| Rizeei et al. | Extraction and accuracy assessment of DTMs derived from remotely sensed and field surveying approaches in GIS framework | |
| Wang et al. | Factors influencing measurement accuracy of unmanned aerial systems (UAS) and photogrammetry in construction earthwork | |
| Alexander et al. | Cost-effective aerial imagery and soil CO2 flux surveys for geothermal exploration | |
| Papp et al. | Effect of the difference between surface and terrain models on gravity field related quantities | |
| CN110110347B (en) | Aviation gravity vector downward continuation method and system based on point mass method | |
| Spiroiu et al. | Geo-gravimetric quasi-geoid determination over Romania |
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 |