CN112526625B - Computing device for abnormal value of Bragg gravity of aviation gravity measurement point - Google Patents
Computing device for abnormal value of Bragg gravity of aviation gravity measurement point Download PDFInfo
- Publication number
- CN112526625B CN112526625B CN202011329110.4A CN202011329110A CN112526625B CN 112526625 B CN112526625 B CN 112526625B CN 202011329110 A CN202011329110 A CN 202011329110A CN 112526625 B CN112526625 B CN 112526625B
- Authority
- CN
- China
- Prior art keywords
- measurement point
- cuboid
- elevation
- gravity
- aviation
- 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
- 230000005484 gravity Effects 0.000 title claims abstract description 203
- 238000005259 measurement Methods 0.000 title claims abstract description 121
- 230000002159 abnormal effect Effects 0.000 title abstract description 9
- 238000012937 correction Methods 0.000 claims abstract description 55
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 238000010276 construction Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 238000012876 topography Methods 0.000 description 11
- 238000009826 distribution Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 8
- 238000003860 storage Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses a calculation device for a Bragg gravity abnormal value of an aviation gravity measurement point, which comprises the following components: the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity anomaly value; the projection module is used for converting the GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographical position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point; constructing a module; a terrain correction value calculation module; and the Bragg gravity anomaly calculation module is used for obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction. The calculation device for the abnormal value of the Bragg gravity of the aviation gravity measurement point can obtain the abnormal value of the aviation Bragg gravity.
Description
Technical Field
The invention relates to the technical field of aviation gravity measurement, in particular to a device for calculating a Bragg gravity abnormal value of an aviation gravity measurement point.
Background
With the rapid development of technology, the precision requirement of gravity measurement in the gravity exploration field is higher and higher.
The data obtained by gravity measurement are gravity observation data. The gravity observation data includes anomalies caused by geologic bodies with uneven underground density, and also includes influences of different surrounding terrains of each measuring point, different latitude, and the like, so that the gravity observation data needs to be corrected. The airborne gravity measurement is a gravity measurement method, and particularly, an airborne gravity measurement system is installed on an airplane to perform continuous measurement.
At present, the Bragg gravity anomaly (Bouguer gravity anomaly) has important significance for processing gravity data. Ground-based Bragg gravity anomalies are typically obtained by ground-based gravity measurements. The gravity obtained by subtracting the normal gravity value after the observation result of the ground gravity instrument is corrected by latitude and altitude, the middle layer is corrected and the topography is corrected is called ground Bragg gravity anomaly.
Based on this, the inventor of the present application found that the calculation method for acquiring the bragg gravity anomaly through the ground gravity measurement value in the prior art has a certain inapplicability in the calculation for acquiring the aviation bragg gravity anomaly through the aviation gravity anomaly.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person of ordinary skill in the art.
Disclosure of Invention
In order to solve the above-mentioned problems, an object of an embodiment of the present invention is to provide a computing device for a bragg gravity anomaly value of an aeronautical gravity measurement point.
In a first aspect, an embodiment of the present invention provides a method for calculating a bragg gravity anomaly value of an airborne gravity measurement point, including: acquiring geographical position information of an aviation gravity measurement point, and GNSS elevation of the aviation gravity measurement point and a gravity anomaly value; converting the GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographic position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point; in a terrain elevation data grid, constructing a first cuboid according to four nodes closest to the elevation projection point as a bottom surface, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises elevation values of the plurality of nodes; respectively taking the four nodes as starting points, taking the step length of a grid in a terrain elevation data grid as a side length, extending towards the direction deviating from a projection point, constructing a second cuboid, wherein the height of the second cuboid is the average value of elevation values corresponding to the four nodes on the bottom surface, and repeating the step of constructing the second cuboid until the coordinate value of the node exceeds a preset range; calculating a terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid; summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point; and obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction.
In one possible implementation, the terrain elevation data grid includes a local terrain elevation data grid and a regional terrain elevation data grid; the preset range of the node coordinate values in the local terrain elevation data grid is 0-20km, and the preset range of the node coordinate values in the regional terrain data grid is 20-166.7 km.
In one possible implementation, the coordinate values of the nodes are beyond a preset range: the coordinate value of the node exceeds 166.7km.
In one possible implementation manner, before the building the first cuboid according to the four nodes closest to the elevation projection point in the terrain elevation data grid as the bottom surface, the method further includes: projecting geographical position information of a measuring point into a rock density distribution data grid to obtain a density projection point, wherein the rock density distribution data grid comprises grid data formed by coordinate values of a plurality of nodes, and the rock density distribution data grid comprises density values of the plurality of nodes; the density of the first cuboid is a density value corresponding to the density projection point.
In one possible implementation, the density value of the second cuboid is an average value of four nodes on the bottom surface of the second cuboid projected into a rock density distribution data grid, respectively.
In one possible implementation, the division rule of the local terrain elevation data grid is high-precision and large-scale; the dividing rule of the regional terrain elevation data grids is low-grade precision and small scale.
In one possible implementation, the computing method further includes: obtaining geographical position information of a next aviation gravity measurement point, GNSS elevation and gravity measurement value of the next aviation gravity measurement point, and repeating the steps of calculating the terrain correction corresponding to the next aviation gravity measurement point until the calculation of the terrain correction of all aviation gravity measurement points is completed; carrying out point-by-point low-pass filtering on the terrain correction values of all the aviation gravity measurement points, wherein the filtering scale of the low-pass filtering corresponds to the filtering scale of the gravity abnormality of the aviation gravity measurement points; the step of obtaining the Bragg gravity anomaly value of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction comprises the following steps: and obtaining the Bragg gravity anomaly value of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the filtered terrain correction.
In a second aspect, an embodiment of the present invention further provides a device for calculating a bragg gravity anomaly value of an airborne gravity measurement point, including: the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity anomaly value; the projection module is used for converting the GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographical position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point; the construction module is used for constructing a first cuboid according to four nodes closest to the elevation projection point in a terrain elevation data grid, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises elevation values of the plurality of nodes; the construction module is further used for constructing a second cuboid by taking the four nodes as starting points and taking the step length of the grid in the terrain elevation data grid as the side length and extending towards the direction deviating from the projection point, the height of the second cuboid is the average value of elevation values corresponding to the four nodes on the bottom surface, and the step of constructing the second cuboid is repeated until the coordinate value of the node exceeds a preset range; the terrain correction value calculation module is used for calculating the terrain correction value of each cuboid according to the height difference between the aviation gravity measurement point and each cuboid, the distance difference between the aviation gravity measurement point and each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid; the terrain correction computing module is used for summing all the terrain correction values to obtain the terrain correction of the aviation gravity measurement point; and the Bragg gravity anomaly calculation module is used for obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction.
In a third aspect, embodiments of the present invention further provide a storage medium storing computer-executable instructions for performing the above-described calculation of a bragg gravity anomaly value for an airborne gravity measurement point.
According to the computing device for the Bragg gravity anomaly value of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point is converted into the elevation of the ground level, the geographic position information of the aviation gravity measurement point is projected into the terrain elevation data grid, a first cuboid is built according to four nodes closest to the elevation projection point as the bottom surface, and the four nodes are respectively used as the starting points to repeatedly build the cuboid; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point; and obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction, so that the calculation of the aviation Bragg gravity anomaly can be realized.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for calculating a Bragg gravity anomaly value of an airborne gravity measurement point provided by an embodiment of the invention;
FIG. 2 illustrates a schematic diagram of a terrain elevation data grid provided by an embodiment of the present invention;
fig. 3 is a schematic diagram of a cuboid calculation method according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a calculation device for a bragg gravity anomaly value of an aeronautical gravity measurement point according to an embodiment of the present invention.
Detailed Description
The following detailed description of embodiments of the invention is, therefore, to be taken in conjunction with the accompanying drawings, and it is to be understood that the scope of the invention is not limited to the specific embodiments.
In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings are merely for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention.
Throughout the specification and claims, unless explicitly stated otherwise, the term "comprise" or variations thereof such as "comprises" or "comprising", etc. will be understood to include the stated element or component without excluding other elements or components.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present invention, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
Because the actual topography fluctuation of the aviation gravity measurement is complex, the embodiment divides the topography around the aviation gravity measurement point into a plurality of cuboids, calculates the gravity influence value of the topography mass of each cuboid on the aviation gravity measurement point, and finally adds up the gravity influence values to obtain the topography influence value of the point. Because the topography around the measuring point is a very irregular curved surface, the topography can not be expressed by a function of a space coordinate, the digital simulation of the topography of the ground is realized by adopting the topography elevation data, and the grid data of the regular topography elevation are established. When the aerogravity terrain correction calculation is carried out, the elevation data are respectively read according to the grid nodes in the local terrain elevation grid data and the regional terrain elevation grid data within a certain range.
The embodiment of the invention provides a flow chart of a calculation method of a Bragg gravity anomaly value of an aviation gravity measurement point, which is shown in fig. 1 and comprises the following steps: step 1-step 7.
Step 1, obtaining geographical position information of an aviation gravity measurement point, GNSS elevation of the aviation gravity measurement point and a gravity anomaly value;
and 2, converting the GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographical position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point.
In one implementation, the terrain elevation data grid includes a local terrain elevation data grid and a regional terrain elevation data grid; the preset range of the node coordinate values in the local terrain elevation data grid is 0-20km, and the preset range of the node coordinate values in the regional terrain data grid is 20-166.7 km. Referring to fig. 2, an exemplary elevation data grid is shown, where EFGH is a local terrain, IJKL is a regional terrain, and the middle black region is the range of the airborne gravity measurement points.
The dividing rule of the local terrain elevation data grid is high precision and large scale; the dividing rule of the regional terrain elevation data grids is low-grade precision and small scale.
Step 3, constructing a first cuboid in a terrain elevation data grid according to four nodes closest to the elevation projection point as the bottom surface, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises elevation values of the plurality of nodes;
step 4, respectively taking the four nodes as starting points, taking the grid length in the terrain elevation data grid as side length, extending towards the direction deviating from the projection point, constructing a second cuboid, wherein the height of the second cuboid is the average value of elevation values corresponding to the four nodes on the bottom surface, and repeating the step of constructing the second cuboid until the coordinate value of the node exceeds a preset range;
the coordinate value of the node exceeds the preset range, namely the coordinate value of the node exceeds 166.7km.
Step 5, calculating a terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;
step 6, summing all terrain correction values to obtain the terrain correction value of the aviation gravity measurement point;
and 7, obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction.
According to the method for calculating the Bragg gravity anomaly value of the aviation gravity measurement point, GNSS elevation of the aviation gravity measurement point is converted into the elevation of the ground level, geographic position information of the aviation gravity measurement point is projected into a terrain elevation data grid, a first cuboid is constructed according to four nodes closest to the elevation projection point as a bottom surface, and the four nodes are respectively used as starting points to repeatedly construct the cuboid; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point; and obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction, so that the calculation of the aviation Bragg gravity anomaly can be realized.
In one implementation, step 3 is preceded by: projecting geographical position information of a measuring point into a rock density distribution data grid to obtain a density projection point, wherein the rock density distribution data grid comprises grid data formed by coordinate values of a plurality of nodes, and the rock density distribution data grid comprises density values of the plurality of nodes; the density of the first cuboid is a density value corresponding to the density projection point.
Correspondingly, the density value of the second cuboid is an average value of four nodes on the bottom surface of the second cuboid projected into the rock density distribution data grid respectively.
In one implementation, step 5 may specifically include:
for each cuboid, calculating an aviation gravity terrain correction value delta g of the cuboid according to a formula I, wherein the formula I comprises:
wherein G is a universal gravitation constant, ρ is the density of each cuboid, X, Y, Z are observation point space positions, and (X1, X2), (Y1, Y2), and (Z1, Z2) are cuboid space positions. R is the distance from the observation point to each vertex of the cuboid.
Further, the implementation manners of step 6 and step 7 are as follows. Referring to fig. 3, which is a schematic diagram of a cuboid calculation method provided in this embodiment, a position of an aviation gravity observation point projected to the terrain height Cheng Wangge data is a point P (X0, Y0), wherein a point A, B, C, D is a node closest to the point P in the terrain grid data. Taking A, B, C, D points as sides of a cuboid, taking the height of the sides as an average value of elevation values of the nodes A, B, C, D, and calculating a terrain influence value delta g from the cuboid to an aviation gravity measurement point 1 . Expanding to two directions of the grid, and calculating the terrain influence value delta g of the fitting cuboid corresponding to each four nodes to the aviation gravity measurement point A by using the method i 。
The sum of the terrain corrections for the aero gravity measurement points is:
the abnormal values of the Bragg gravity of the aviation gravity measurement points are as follows:
G bragg cell ＝G Self-air -G Ground (floor) (2)
Wherein G is Self-air The spatial gravity anomaly value for the aero gravity measurement point a is obtained by aero gravity measurement.
In one implementation, step 7 may further include: and obtaining the geographical position information of the next aviation gravity measurement point, the GNSS elevation of the next aviation gravity measurement point and the gravity measurement value, and repeating the steps of calculating the terrain correction corresponding to the next aviation gravity measurement point until the calculation of the terrain correction of all aviation gravity measurement points is completed.
Specifically, the calculation steps of the steps 1 to 6 are repeated for each airborne gravity measurement point until the topography correction of all the airborne gravity measurement points is calculated.
In one implementation, the method may further include: and carrying out point-by-point low-pass filtering on the terrain correction values of all the aviation gravity measurement points, wherein the filtering scale of the low-pass filtering corresponds to the filtering scale of the aviation space gravity anomaly.
Accordingly, step 7 may be: and obtaining the Bragg gravity anomaly value of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the filtered terrain correction.
The embodiment of the invention provides a schematic structural diagram of a device for calculating the abnormal value of the Bragg gravity of the aviation gravity measurement point, which is shown in fig. 4, and comprises the following steps: the system comprises an acquisition module 1, a projection module 2, a construction module 3, a terrain correction value calculation module 4, a terrain correction value calculation module 5 and a Bragg gravity anomaly calculation module 6.
The acquisition module 1 is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity anomaly value;
the projection module 2 is used for converting GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographical position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point;
the construction module 3 is configured to construct a first cuboid according to four nodes closest to the elevation projection point as a bottom surface in a terrain elevation data grid, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises elevation values of the plurality of nodes;
the building module 3 is further configured to build a second cuboid by using the four nodes as starting points and using a grid length in a terrain elevation data grid as a side length, extending in a direction away from the projection point, where a height of the second cuboid is an average value of elevation values corresponding to the four nodes on the bottom surface, and repeating the step of building the second cuboid until coordinate values of the nodes exceed a preset range;
the terrain correction value calculating module 4 is used for calculating the terrain correction value of each cuboid according to the height difference from the aviation gravity measuring point to each cuboid, the distance difference from the aviation gravity measuring point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;
the terrain correction computing module 5 is used for summing all the terrain correction values to obtain the terrain correction of the aviation gravity measurement point;
and the Bragg gravity anomaly calculation module 6 is used for obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction.
According to the computing device for the Bragg gravity anomaly value of the aviation gravity measurement point, GNSS elevation of the aviation gravity measurement point is converted into the elevation of the ground level, geographic position information of the aviation gravity measurement point is projected into a terrain elevation data grid, a first cuboid is built according to four nodes closest to the elevation projection point as the bottom surface, and the four nodes are respectively used as the starting points to repeatedly build the cuboid; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point; and obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction, so that the calculation of the aviation Bragg gravity anomaly can be realized.
The embodiment of the invention also provides a storage medium, which stores computer executable instructions, including a program for executing the method for calculating the abnormal value of the Bragg gravity of the aviation gravity measurement point, wherein the computer executable instructions can execute the method in any method embodiment.
The storage medium may be any available medium or data storage device that can be accessed by a computer, including, but not limited to, magnetic storage (e.g., floppy disks, hard disks, magnetic tape, magneto-optical disks (MOs), etc.), optical storage (e.g., CD, DVD, BD, HVD, etc.), and semiconductor storage (e.g., ROM, EPROM, EEPROM, nonvolatile storage (NAND FLASH), solid State Disk (SSD)), etc.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing descriptions of specific exemplary embodiments of the present invention are presented for purposes of illustration and description. It is not intended to limit the invention to the precise form disclosed, and obviously many modifications and variations are possible in light of the above teaching. The exemplary embodiments were chosen and described in order to explain the specific principles of the invention and its practical application to thereby enable one skilled in the art to make and utilize the invention in various exemplary embodiments and with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims and their equivalents.
Claims (1)
1. A computing device for a bragg gravity anomaly value for an airborne gravity measurement point, comprising:
the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity anomaly value;
the projection module is used for converting the GNSS elevation of the aviation gravity measurement point into the elevation of the ground level, and projecting the geographical position information of the aviation gravity measurement point into a terrain elevation data grid to obtain an elevation projection point;
the construction module is used for constructing a first cuboid according to four nodes closest to the elevation projection point in a terrain elevation data grid, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises elevation values of the plurality of nodes;
the construction module is further used for constructing a second cuboid by taking the four nodes as starting points and taking the grid step length in the terrain elevation data grid as side length and extending in the direction deviating from the elevation projection point, wherein the height of the second cuboid is the average value of elevation values corresponding to the four nodes on the bottom surface; repeating the step of constructing the second cuboid until the coordinate value of the node exceeds a preset range;
the terrain correction value calculation module is used for calculating the terrain correction value of each cuboid according to the height difference between the aviation gravity measurement point and each cuboid, the distance difference between the aviation gravity measurement point and each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;
the terrain correction computing module is used for summing all the terrain correction values to obtain the terrain correction of the aviation gravity measurement point;
and the Bragg gravity anomaly calculation module is used for obtaining the Bragg gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011329110.4A CN112526625B (en) | 2020-06-08 | 2020-06-08 | Computing device for abnormal value of Bragg gravity of aviation gravity measurement point |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011329110.4A CN112526625B (en) | 2020-06-08 | 2020-06-08 | Computing device for abnormal value of Bragg gravity of aviation gravity measurement point |
| CN202010511402.3A CN111679336B (en) | 2020-06-08 | 2020-06-08 | Calculation method and device for bump gravity abnormal value of aviation gravity measurement point |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010511402.3A Division CN111679336B (en) | 2020-06-08 | 2020-06-08 | Calculation method and device for bump gravity abnormal value of aviation gravity measurement point |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112526625A CN112526625A (en) | 2021-03-19 |
| CN112526625B true CN112526625B (en) | 2024-01-30 |
Family
ID=72454001
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011329110.4A Active CN112526625B (en) | 2020-06-08 | 2020-06-08 | Computing device for abnormal value of Bragg gravity of aviation gravity measurement point |
| CN202010511402.3A Expired - Fee Related CN111679336B (en) | 2020-06-08 | 2020-06-08 | Calculation method and device for bump gravity abnormal value of aviation gravity measurement point |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010511402.3A Expired - Fee Related CN111679336B (en) | 2020-06-08 | 2020-06-08 | Calculation method and device for bump gravity abnormal value of aviation gravity measurement point |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN112526625B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112462443B (en) * | 2020-11-13 | 2021-07-20 | 中国自然资源航空物探遥感中心 | Synchronous terrain correction method and device for aerial gravity measurement |
| RU2764108C1 (en) * | 2020-12-03 | 2022-01-13 | Иван Владимирович Бригадин | Method for masking underground objects |
| CN115629426A (en) * | 2022-10-24 | 2023-01-20 | 中国自然资源航空物探遥感中心 | Method, system, device and medium for extracting local gravity anomaly |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166170A (en) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | Method for positioning and detecting low-density concealed ore bodies in gallery gravity total spatial domain |
| WO2015104633A2 (en) * | 2014-01-13 | 2015-07-16 | Cgg Services Sa | System and method for two dimensional gravity modeling with variable densities |
| CN110646858A (en) * | 2019-10-12 | 2020-01-03 | 山东省物化探勘查院 | Submarine gravity measurement middle-far one-region terrain correction calculation method |
| CN110706275A (en) * | 2019-10-16 | 2020-01-17 | 中国石油大学(华东) | Local gravity anomaly extraction method based on satellite altimetry gravity data |
| CN110967778A (en) * | 2019-10-24 | 2020-04-07 | 西北大学 | Dynamic coordinate system polyhedral subdivision gravity grid distribution correction method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170108612A1 (en) * | 2015-10-15 | 2017-04-20 | King Saud University | Inertial system for gravity difference measurement |
| CN108008459B (en) * | 2017-11-28 | 2019-11-01 | 北京中科地物能源技术有限公司 | A kind of method and device obtaining residual gravity anomaly |
-
2020
- 2020-06-08 CN CN202011329110.4A patent/CN112526625B/en active Active
- 2020-06-08 CN CN202010511402.3A patent/CN111679336B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015104633A2 (en) * | 2014-01-13 | 2015-07-16 | Cgg Services Sa | System and method for two dimensional gravity modeling with variable densities |
| CN104166170A (en) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | Method for positioning and detecting low-density concealed ore bodies in gallery gravity total spatial domain |
| CN110646858A (en) * | 2019-10-12 | 2020-01-03 | 山东省物化探勘查院 | Submarine gravity measurement middle-far one-region terrain correction calculation method |
| CN110706275A (en) * | 2019-10-16 | 2020-01-17 | 中国石油大学(华东) | Local gravity anomaly extraction method based on satellite altimetry gravity data |
| CN110967778A (en) * | 2019-10-24 | 2020-04-07 | 西北大学 | Dynamic coordinate system polyhedral subdivision gravity grid distribution correction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111679336A (en) | 2020-09-18 |
| CN112526625A (en) | 2021-03-19 |
| CN111679336B (en) | 2020-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112526625B (en) | Computing device for abnormal value of Bragg gravity of aviation gravity measurement point | |
| AU2016305571B2 (en) | System and method for gravity and/or gravity gradient terrain corrections | |
| CN101839710B (en) | Method for optimizing quasi-geoid calculation | |
| CN102645209B (en) | Joint positioning method for spatial points by means of onboard LiDAR point cloud and high resolution images | |
| CN112363236B (en) | Gravity field data equivalent source continuation and data type conversion method based on PDE | |
| CN103217688B (en) | Airborne laser radar point cloud adjustment computing method based on triangular irregular network | |
| CN109934911B (en) | OpenGL-based three-dimensional modeling method for high-precision oblique photography of mobile terminal | |
| CN103256914B (en) | A kind of method and system calculating silt arrester inundated area based on DEM | |
| CN103278115A (en) | Method and system for calculating deposition volume of check dam based on DEM (digital elevation model) | |
| CN111257956A (en) | Matlab-based regional quasi-geoid surface refinement method | |
| CN113311494A (en) | Satellite gravity field inversion method | |
| Yildiz | A study of regional gravity field recovery from GOCE vertical gravity gradient data in the Auvergne test area using collocation | |
| CN113671570B (en) | Seismic surface wave travel time and gravity anomaly joint inversion method and system | |
| CN113109883B (en) | Method for forward modeling satellite gravity field based on isoparametric transformation global discrete grid spherical coordinates | |
| CN102650703B (en) | CPM-CGCS2000 building method for CGCS2000 plate model | |
| Jekeli et al. | Geoid determination in South Korea from a combination of terrestrial and airborne gravity anomaly data | |
| CN104240301A (en) | Geological surface reconstruction method and device | |
| CN115112090A (en) | Mountain area digital elevation datum establishing method based on few GNSS leveling points | |
| CN113552644B (en) | Density determination method and device and electronic equipment | |
| Jekeli | Deflections of the vertical from full-tensor and single-instrument gravity gradiometry | |
| CN110968656B (en) | Topographic map processing method and device | |
| US20240085182A1 (en) | Methods for lifting and lowering magnetic anomalies for magnetic navigation | |
| Grzyb et al. | The effect of topography and quality of a digital terrain model on the accuracy of terrain corrections for centimetre quasigeoid modelling | |
| CN114964169B (en) | Remote sensing image adjustment method for image space object space cooperative correction | |
| CN104614730B (en) | A kind of method using lunar surface laser reflection battle array coordinate calibration moon elevation model |
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 |