A method, system, device and medium for determining artificial terrain regional gravity anomalies

Abstract

The application discloses a kind of artificial terrain area gravity anomaly determination method, system, equipment and medium, it is related to geological exploration field.The method comprises the following steps: determining the center position coordinates of ground control point signboard;Determine the image data in artificial terrain area;According to the center position coordinates, correct image data, obtain topographic data;Determine the plane coordinates, geographic latitude and elevation of all observation points in artificial terrain area;Determine the observed gravity value of all observation points in artificial terrain area;According to the plane coordinates, geographic latitude, elevation, observed gravity value and topographic data of all observation points in artificial terrain area, calculate the gravity anomaly value of all observation points in artificial terrain area;The gravity anomaly value of all observation points in artificial terrain area is gridded, to determine the gravity anomaly variation characteristics of artificial terrain area, complete the missing data of gravity anomaly of artificial terrain area.The application can solve the problem of missing gravity anomaly data in artificial terrain area.

Description

Technical Field

[0001] This invention relates to the field of geological exploration, and in particular to a method, system, equipment and medium for determining gravity anomalies in artificial terrain areas. Background Technology

[0002] Surface gravity exploration is a geophysical method that studies gravity variations (called gravity anomalies) caused by uneven distribution of underground material density to understand and infer the Earth's structure, crustal formation, and to explore for mineral resources. The first step in surface gravity exploration is to measure gravity values ​​on the ground using a high-precision gravimeter. These measurements are influenced by numerous external factors and are not entirely caused by underground geological bodies. They include topographic relief, elevation changes at the measuring point, the fact that the Earth is not a sphere, and gravity variations caused by the Earth's rotation. Only by removing these influences can the gravity anomalies caused by uneven distribution of underground material density be obtained. Topographic relief has the most significant impact on gravity measurements and is the most important component of gravity anomaly errors. Therefore, the primary problem to be solved in obtaining high-quality gravity anomalies is to properly handle the influence of topographic relief.

[0003] In surface gravity exploration, gravity anomalies are often derived from collected raw topographic data. However, in actual measurements, artificial terrain areas such as mines and quarries are frequently encountered. The actual terrain undulations differ significantly from the raw topographic data, leading to distorted gravity anomalies measured in these artificial terrain areas, necessitating the omission of these points. This results in incomplete data, missing anomaly values ​​from the mine area, which inevitably affects the reasonable interpretation of the target. Studying the deep resource potential of known mineral deposits is currently an important direction for mineral exploration. Surface gravity exploration is a crucial means of exploring deep geological structures. Gravity anomalies can effectively reflect the distribution characteristics of deep fault structures and intrusive rocks, which are essential geological conditions for basic geological research and mineralization.

[0004] Therefore, existing technologies suffer from incomplete gravity anomaly data, making it impossible to accurately study geological phenomena. Summary of the Invention

[0005] The purpose of this invention is to provide a method, system, equipment, and medium for determining gravity anomalies in artificial terrain areas, so as to solve the problem of missing gravity anomaly data in artificial terrain areas and achieve accurate research on geological phenomena.

[0006] To achieve the above objectives, the present invention provides the following solution:

[0007] A method for determining gravity anomalies in an artificial terrain area, wherein the method is applied to a gravity anomaly observation device for the artificial terrain area, the observation device comprising: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers, the ground control point markers being placed within the artificial terrain area; the method for determining gravity anomalies in the artificial terrain area includes:

[0008] Using the GPS measuring instrument, the center coordinates of the ground control point sign were determined;

[0009] The drone was used to measure image data within the artificial terrain area;

[0010] The image data is corrected based on the center location coordinates to obtain terrain data;

[0011] Using the GPS measuring instrument, the plane coordinates, geographical latitude, and altitude of all observation points within the artificial terrain area were determined;

[0012] Using the gravimeter, the observed gravity values ​​at all observation points within the artificial terrain area were measured;

[0013] Based on the plane coordinates, geographical latitude, altitude, observed gravity values, and topographic data of all observation points within the artificial terrain area, calculate the gravity anomaly values ​​of all observation points within the artificial terrain area.

[0014] The gravity anomaly values ​​of all observation points within the artificial terrain area are gridded based on the planar coordinates of all observation points within the artificial terrain area to determine the gravity anomaly variation characteristics of the artificial terrain area and to supplement the missing gravity anomaly data of the artificial terrain area.

[0015] Optionally, based on the planar coordinates, geographical latitude, elevation, observed gravity values, and the terrain data of all observation points within the artificial terrain area, the gravity anomaly values ​​of all observation points within the artificial terrain area are calculated, specifically including:

[0016] Calculate the Bouguer correction values ​​for all observation points within the artificial terrain area based on their geographical latitude and elevation.

[0017] Based on the plane coordinates, elevation, and topographic data of all observation points within the artificial terrain area, calculate the terrain correction values ​​for all observation points within the artificial terrain area.

[0018] Calculate the normal gravity values ​​of all observation points within the artificial terrain area based on their geographical latitude.

[0019] Based on the observed gravity values, normal gravity values, Bouguer correction values, and terrain correction values ​​of all observation points within the artificial terrain area, calculate the gravity anomaly values ​​of all observation points within the artificial terrain area.

[0020] Optionally, based on the geographical latitude and elevation of all observation points within the artificial terrain area, the Bouguer correction value for all observation points within the artificial terrain area is calculated, using the following formula:

[0021]

[0022] Where: δ gB This is the Bouguer correction value, in units of 10. -5 m / s 2 ; ρ1 is the geographical latitude of the observation point; h is the altitude of the observation point; ρ1 is the density of the intermediate layer; a is the correction radius of the intermediate layer of the circular domain.

[0023] Optionally, based on the plane coordinates, elevation, and topographic data of all observation points within the artificial terrain area, the topographic correction value for all observation points within the artificial terrain area is calculated using the following formula:

[0024]

[0025] Among them, g gT ρ is the topographic correction value; G is the gravitational constant; ρ is the average density of the Earth's crust; l is the integral grid spacing; C ij r is the integration constant, specifically the trapezoidal coefficient. ij h is the distance between the integration node (i, j) and the observation point; ij The elevation difference between the integration node (i, j) and the observation point is denoted by ; i is the row number of the integration node; j is the column number of the integration node; the integration node is determined based on the terrain data; the distance between the integration node and the observation point is determined based on the plane coordinates of the observation point; the elevation difference between the integration node and the observation point is determined based on the altitude of the observation point.

[0026] Optionally, the normal gravity values ​​of all observation points within the artificial terrain area are calculated based on their geographical latitude. The specific formula is as follows:

[0027]

[0028] Where g0 is the normal gravity value, in units of 10. -5 m / s 2 ; The latitude of the observation point.

[0029] Optionally, based on the observed gravity values, normal gravity values, Bouguer corrections, and topographic corrections of all observation points within the artificial terrain area, the gravity anomaly values ​​of all observation points within the artificial terrain area are calculated using the following formula:

[0030] Δg B =g - g0 + δ gB +g gT

[0031] Where: Δg B δ represents the gravity anomaly; g represents the observed gravity value; g0 represents the normal gravity value; δ gB g is the Bouguer correction value. gT This is the terrain correction value.

[0032] Optionally, the image data is corrected based on the center location coordinates to obtain terrain data, specifically including:

[0033] Based on the center position coordinates, aerial triangulation is performed on the image data to obtain a three-dimensional elevation data volume;

[0034] The three-dimensional elevation data volume is processed by data synthesis and format conversion to obtain terrain data.

[0035] A gravity anomaly determination system for artificial terrain areas, wherein the system is applied to a gravity anomaly observation device for artificial terrain areas, the observation device comprising: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers, the ground control point markers being placed within the artificial terrain area; the gravity anomaly determination system for artificial terrain areas includes:

[0036] The center position coordinate acquisition module is used to determine the center position coordinates of the ground control point sign using the GPS measuring instrument;

[0037] The image data acquisition module is used to measure image data within the artificial terrain area using the UAV;

[0038] The terrain data acquisition module is used to correct the image data according to the center location coordinates to obtain terrain data;

[0039] The observation point information acquisition module is used to determine the plane coordinates, geographical latitude, and altitude of all observation points within the artificial terrain area using the GPS measuring instrument.

[0040] The gravity value acquisition module is used to measure the gravity values ​​of all observation points within the artificial terrain area using the gravimeter.

[0041] The gravity anomaly determination module is used to calculate the gravity anomaly value of all observation points within the artificial terrain area based on the plane coordinates, geographical latitude, altitude, observed gravity value, and the terrain data of all observation points within the artificial terrain area.

[0042] The gridding module is used to grid the gravity anomaly values ​​of all observation points within the artificial terrain area based on the planar coordinates of all observation points within the artificial terrain area, determine the gravity anomaly variation characteristics of the artificial terrain area, and supplement the missing gravity anomaly data of the artificial terrain area.

[0043] An electronic device includes a memory and a processor, the memory storing a computer program, and the processor running the computer program to cause the electronic device to perform the above-described method for determining gravity anomalies in artificial terrain areas.

[0044] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for determining gravity anomalies in artificial terrain areas.

[0045] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0046] The present invention provides a method for determining gravity anomalies in artificial terrain areas. This method utilizes a GPS measuring instrument to determine the center coordinates of ground control point markers and the planar coordinates, latitude, and elevation of all observation points within the artificial terrain area. It also utilizes a drone to measure image data within the artificial terrain area and corrects the image data based on the center coordinates to obtain terrain data. Finally, it uses a gravimeter to measure the observed gravity values ​​at all observation points within the artificial terrain area. Based on this data, the method calculates the gravity anomaly values ​​at all observation points within the artificial terrain area and grids these values ​​to determine the gravity anomaly variation characteristics of the artificial terrain area, thus supplementing missing gravity anomaly data. Therefore, this invention solves the problem of missing gravity anomaly data in artificial terrain areas, which is beneficial for geological structure research and mineral exploration. Attached Figure Description

[0047] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0048] Figure 1 A flowchart of the method for determining gravity anomalies in artificial terrain areas provided by the present invention;

[0049] Figure 2 A schematic diagram of a device for marking ground control points;

[0050] Figure 3 A schematic diagram illustrating the relationship between flight altitude and ground resolution;

[0051] Figure 4 A flowchart of the aerial triangulation process;

[0052] Figure 5 This is a spurious gravity anomaly map of the artificial terrain area in the experimental zone;

[0053] Figure 6 This is a gravity anomaly map of the artificial terrain area in the test zone obtained using the present invention;

[0054] Figure 7 A block diagram of the gravity anomaly determination system for artificial terrain areas provided by the present invention.

[0055] Symbol explanation:

[0056] Center location coordinate acquisition module—1, image data acquisition module—2, terrain data acquisition module—3, observation point information acquisition module—4, observation gravity value acquisition module—5, gravity anomaly value determination module—6, gridding module—7. Detailed Implementation

[0057] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0058] The purpose of this invention is to provide a method, system, equipment, and medium for determining gravity anomalies in artificial terrain areas, so as to solve the problem of missing gravity anomaly data in artificial terrain areas and achieve accurate research on geological phenomena.

[0059] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0060] Example 1

[0061] like Figure 1 As shown, the present invention provides a method for determining gravity anomalies in artificial terrain areas. The method is applied to a gravity anomaly observation device for artificial terrain areas. The gravity anomaly observation device for artificial terrain areas includes: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers. The ground control point markers are placed within the artificial terrain area to form ground control points.

[0062] The method for determining gravity anomalies in the artificial terrain area includes:

[0063] Step S1: Use the GPS measuring instrument to determine the center coordinates of the ground control point sign.

[0064] Specifically, in areas with artificial terrain, after the ground control point sign is fixed in place, the coordinates of the center of the ground control point sign are determined using the GPS measuring instrument; that is, the center position coordinates. The placement of the ground control point sign is mainly determined based on the designed point coordinates. The GPS measuring instrument is used to navigate to the location of the point, and the ground control point sign is fixed according to the actual situation. The coordinates of the actual center position of the ground control point sign are then determined.

[0065] When setting up ground control point signs, the following should be noted:

[0066] 1. Ground control point signs should be placed in locations that are easy to interpret, convenient for joint measurement, and have minimal elevation changes.

[0067] 2. Ground control point signs should be placed on fixed ground features with open views.

[0068] 3. Ground control point signs are roughly evenly distributed in the artificial terrain area.

[0069] As one specific implementation method, Figure 2 This is a schematic diagram of a ground control point sign device provided by the present invention. Figure 2 As shown, the ground control point sign is 100cm long and 100cm wide. To ensure more accurate data measurement by the drone, a fixing device for the ground control point sign is also installed to prevent displacement. The frame of the ground control point sign is made of metal to ensure stability, specifically aluminum alloy; the sign face is made of brightly colored, high-contrast material, specifically red and white cloth.

[0070] Step S2: Using the aforementioned UAV, measure the image data within the artificial terrain area. The area of ​​terrain data measured by the UAV should be larger than the area of ​​the artificial terrain to facilitate stitching with surrounding terrain data. In this step, attention needs to be paid to aspects such as flight path layout, flight altitude design, photography time, supplementary and repeated photography, and UAV image data acquisition, as detailed below:

[0071] 1. Route layout

[0072] Based on the survey area and elevation data, flight path design software is used, combined with the elevations of the lowest and highest points in the survey area, to design flight paths. Flight paths are designed based on the DEM of the survey area, with the lowest point within the area selected as the reference plane. Aerial survey factors such as absolute flight altitude, heading and lateral overlap, flight path spacing, number of flight paths, flight path length, and number of images within the survey area are designed. Flight paths generally follow a straight east-west direction. Under specific conditions, adjustments can be made according to the actual terrain to ensure properly overlapping stereo pairs. The heading overlap of the flight paths is set between 60-80%, and the lateral overlap is between 50-80%.

[0073] 2. Flight Altitude Design

[0074] Ground resolution (GSD) in aerial photography depends on flight altitude; the flight altitude should be designed appropriately according to surveying requirements. The relationship between ground resolution and flight altitude is as follows: Figure 3 As shown, we have:

[0075]

[0076] Therefore, we have:

[0077]

[0078] In the formula: H—relative flight altitude; f—lens focal length; A—pixel size; GSD—ground resolution. The flight altitude corresponding to the ground resolution GSD can be calculated using the above formula.

[0079] 3. Photography time

[0080] Aerial photography should be conducted under the most favorable weather conditions in the area, ensuring sufficient illumination while avoiding excessive shadows. The selection is generally based on the solar altitude angle and shadow ratio of the area.

[0081] 4. Reshoots and retakes

[0082] Absolute gaps, relative gaps, and other serious defects discovered during aerial photography must be promptly corrected. Correction photography should generally follow the original flight path, with both ends of the correction path extending beyond the gap by two points.

[0083] 5. UAV image data acquisition

[0084] Unmanned aerial vehicle (UAV) image data acquisition utilizes a multi-rotor UAV equipped with an oblique photography camera, following a pre-planned flight path to collect image data. The acquired image data is then standardized in naming, storage, and backup. To ensure data acquisition quality during UAV image acquisition, the following points should be noted:

[0085] (1) Image overlap control

[0086] Because aerial imagery employs multi-view matching algorithms for aerial triangulation, a high degree of image overlap is required to match more corresponding points. Generally, the forward overlap of the imagery from the downward-looking camera should be no less than 60%, but it should not be too large. Excessive overlap will not only reduce acquisition efficiency but also significantly impact post-processing speed. The lateral overlap of the downward-looking imagery should be designed to be 50%–80%.

[0087] (2) Coverage guarantee control of the shooting area

[0088] To ensure images can be captured even outside the photographic area, the coverage of the photographic area boundary is larger than that of perpendicular photography. The theoretical and actual values ​​for the number of baselines and flight paths whose coverage extends beyond the boundary line are calculated using the following formulas:

[0089]

[0090] In the formula: p is the heading or lateral overlap; θ is the tilt angle; β is the field of view angle.

[0091] In actual flight, due to the influence of various factors such as the atmosphere, the actual value of the heading or lateral coverage exceeding the boundary line is generally calculated according to the following formula:

[0092] Baseline number = Theoretical value + 2;

[0093] Number of routes = theoretical value + 1;

[0094] In actual aerial photography, 4 to 5 additional flight paths can be added outside the boundary of the survey area in the direction of the flight path to ensure that the images from the left and right view lenses cover the entire survey area. Outside the boundary of the survey area in the direction of the flight path, each flight path can be extended by more than 1.5 km to ensure that the images from the forward and backward view lenses cover the entire survey area.

[0095] (3) Altitude maintenance control

[0096] Regarding altitude maintenance, the altitude difference between adjacent images on the same route should not exceed 3m, the difference between the maximum and minimum altitude should not exceed 5m, and the difference between the actual altitude and the designed altitude within the area should not exceed 5m (except for manually encrypted shooting).

[0097] (4) Image quality control

[0098] Image quality inspection mainly checks for phenomena such as clouds, cloud shadows, smoke, fog, and reflections in the images.

[0099] (5) Mapping quality control

[0100] The ground resolution of the DOM in the image is set at 0.5m. The positional error of a feature point relative to its corresponding feature point, the horizontal error and the elevation error of the check point shall not exceed the provisions in the table below. For special areas, the allowance may be relaxed by 0.5 times, and twice the standard error shall be taken as the maximum error.

[0101] Step S3: Correct the image data according to the center position coordinates to obtain terrain data. This step specifically includes:

[0102] Step S3.1: Based on the center position coordinates, perform aerial triangulation on the image data to obtain a three-dimensional elevation data volume.

[0103] Step S3.2: Perform data synthesis and format conversion on the three-dimensional elevation data volume to obtain terrain data.

[0104] As a specific implementation method, such as Figure 4 As shown, aerial triangulation calculations are performed using image data from UAV photography, POS data (i.e., the three-dimensional coordinates and flight attitude of the UAV at the moment of taking the photo, where the three-dimensional coordinates include longitude, latitude, and flight altitude, and the flight attitude includes heading angle, pitch angle, and roll angle), and image control point measurement results, employing 3D mapping software. The overall aerial triangulation workflow is as follows: First, data analysis and area network division are performed. Then, new projects are created, automatic interior orientation is performed, densification point measurements are conducted, and free network adjustment is performed. Densification points are repeatedly edited until they are qualified. Next, image control points and checkpoint measurements are performed, and area network adjustment is performed. Image control points are repeatedly edited until they are qualified. Finally, bundle adjustment and area network edge joining are performed. Common points are repeatedly edited until the scale accuracy requirements are met, and then the results are output.

[0105] The specific steps of the three-dimensional encryption process include:

[0106] (1) Correct the distortion of the original image according to the camera calibration data.

[0107] (2) Relative orientation. Establish a project, perform densification point matching on the distortion-corrected image, perform free network adjustment after matching, remove gross errors until the relative orientation requirement is met.

[0108] (3) Absolute orientation. After adding control points, perform bundle adjustment of the area network until the scale accuracy requirements are met.

[0109] The specific process of outputting and organizing the encrypted aerial triangulation results includes: outputting the encrypted aerial triangulation results according to the requirements of subsequent work, and organizing the output results to prepare for subsequent modeling work.

[0110] The terrain data volume output process specifically includes: performing aerial triangulation on the collected image data to restore the true position and attitude of the camera; generating a high-precision three-dimensional elevation data volume under automatic software calculation; synthesizing and outputting DEM data; and then converting it into the required data format.

[0111] Step S4: Using the GPS measuring instrument, determine the plane coordinates, geographical latitude, and altitude (i.e., the location information of the observation points) of all observation points within the artificial terrain area.

[0112] Step S5: Using the gravimeter, measure the observed gravity values ​​(i.e., gravity field data) at all observation points within the artificial terrain area. The gravimeter used in this invention is briefly described below:

[0113] 1. Instrument name: CG-5, which is currently the mainstream gravity instrument in China. Other models include LCR-G.

[0114] 2. Manufacturer: Scintrex, Canada.

[0115] 3. Instrument Composition: The gravimeter consists of a main unit and auxiliary parts such as a tripod.

[0116] 4. Gravity field measurement process:

[0117] (1) Place the tripod: On a flat ground at the measuring point, press down on the tripod so that the steel tips of each leg are inserted into the ground in a "T" shape. This makes it easier to adjust the bottom screw leveling instrument.

[0118] (2) Instrument placement: Place the instrument stably on the tripod, ensuring that the V-shaped groove on the bottom of the instrument is in contact with the ball end of the tripod foot screw. This will fix the instrument on the tripod. Turn on the power and level the instrument.

[0119] (3) Data acquisition: After determining the observation point information, data acquisition begins, and the measurement results are stored in the memory along with the date and time.

[0120] Step S6: Based on the planar coordinates, latitude, elevation, observed gravity values, and topographic data of all observation points within the artificial terrain area, calculate the gravity anomaly values ​​of all observation points within the artificial terrain area. Further, this step specifically includes:

[0121] Step S6.1: Calculate the Bouguer correction value for all observation points within the artificial terrain area based on their geographical latitude and elevation.

[0122] Step S6.2: Calculate the terrain correction value for all observation points within the artificial terrain area based on the plane coordinates, elevation, and terrain data of all observation points within the artificial terrain area.

[0123] Step S6.3: Calculate the normal gravity value of all observation points within the artificial terrain area based on their geographical latitude.

[0124] Step S6.4: Calculate the gravity anomaly value of all observation points within the artificial terrain area based on the observed gravity value, normal gravity value, Bouguer correction value, and terrain correction value of all observation points within the artificial terrain area.

[0125] Specifically, the formula for calculating gravity anomaly values ​​is as follows:

[0126] Δg B =g - g0 + δ gB +g gT ;

[0127] In the formula:

[0128] Δg B —Gravity anomaly;

[0129] g0—Normal gravity value;

[0130] g — observed gravity value;

[0131] δ gB —Bougue correction;

[0132] g gT — Terrain correction value.

[0133] The formula for calculating the normal gravity value is:

[0134]

[0135] In the formula:

[0136] —The geographical latitude of the observation point.

[0137] The formula for calculating the Bouguer correction is:

[0138]

[0139] In the formula:

[0140] h——Elevation of the observation point (negative value when the observation point is below the elevation datum, unit: m);

[0141] ρ1 — density of the intermediate layer (2.67 g / cm³) 3 );

[0142] a——Corrected radius of the intermediate layer of the circular domain (20000m).

[0143] When calculating the terrain correction value of the gravity observation point, the terrain data in the artificial terrain area is first divided into quadrangular prisms. The terrain correction values ​​of the four nodes near the observation point are calculated first. When calculating the terrain correction value of the nodes, the elevation value of the observation point is used to replace the elevation value of the four nodes. Then, the terrain correction values ​​of the four nodes are interpolated to the location of the observation point as the terrain correction value of the observation point.

[0144] The formula for calculating terrain correction values ​​is:

[0145]

[0146] In the formula:

[0147] G – Gravitational constant, taken as 6.67 × 10⁻⁶ -8 cm 3 / (g·s 2 );

[0148] ρ—average density of the Earth's crust, taken as 2.67 g / cm³ 3 ;

[0149] l—Integral grid spacing, selected as 1000m;

[0150] C ij —Integral constant, specifically the trapezoidal coefficient;

[0151] r ij —The distance between the integration node (i, j) and the calculation point (i.e., the observation point) is determined based on the planar coordinates of the observation point;

[0152] h ij —The elevation difference between the integration node (i, j) and the calculation point (i.e., the observation point) is determined based on the elevation of the observation point;

[0153] i — row number of the integration node;

[0154] j — the column number of the integration node.

[0155] The specific values ​​of the integration constant are:

[0156]

[0157] That is: for the internal nodes of the trapezoid, the integral constant is 1; for the edge points of the trapezoid, the integral constant is 0.5; for the external corner points of the trapezoid, the integral constant is 0.25; and for the internal corner points of the trapezoid, the integral constant is 0.75.

[0158] Step S7: Based on the planar coordinates of all observation points within the artificial terrain area, the gravity anomaly values ​​of all observation points within the artificial terrain area are gridded (that is, the gravity anomaly values ​​of observation points that are not uniformly distributed in space are converted into representative values ​​in a regular grid), to determine the gravity anomaly variation characteristics of the artificial terrain area and to complete the missing gravity anomaly data of the artificial terrain area.

[0159] Figure 5 This is a spurious gravity anomaly map of the artificial terrain area in the experimental zone, such as... Figure 5 As shown, in artificial terrain areas, gravity anomalies are false anomalies, so the true gravity anomaly changes in artificial terrain areas cannot be known. Figure 6 This is a gravity anomaly map of the artificial terrain area in the test zone obtained using this invention. (See image below.) Figure 6 As shown, after using drones to acquire gravity anomalies in artificial terrain areas, complete gravity anomaly data is obtained, and the true gravity anomalies in the artificial terrain areas are discovered. This allows us to obtain the characteristics of the anomalies, analyze their geological origins, and conduct basic geological research and mineral exploration prediction.

[0160] Example 2

[0161] To implement the method corresponding to Embodiment 1 above and achieve the corresponding functions and technical effects, a gravity anomaly determination system for artificial terrain areas is provided below. This system is applied to a gravity anomaly observation device for artificial terrain areas. The observation device includes: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers, which are placed within the artificial terrain area. Figure 7 As shown, the gravity anomaly determination system for the artificial terrain area includes:

[0162] The center position coordinate acquisition module 1 is used to determine the center position coordinates of the ground control point sign using the GPS measuring instrument.

[0163] Image data acquisition module 2 is used to measure image data within the artificial terrain area using the UAV.

[0164] The terrain data acquisition module 3 is used to correct the image data according to the center position coordinates to obtain terrain data.

[0165] The observation point information acquisition module 4 is used to determine the plane coordinates, geographical latitude and altitude of all observation points within the artificial terrain area using the GPS measuring instrument.

[0166] The gravity value acquisition module 5 is used to measure the gravity values ​​of all observation points within the artificial terrain area using the gravimeter.

[0167] The gravity anomaly determination module 6 is used to calculate the gravity anomaly value of all observation points within the artificial terrain area based on the plane coordinates, geographical latitude, altitude, observed gravity value, and the terrain data of all observation points within the artificial terrain area.

[0168] The gridding module 7 is used to grid the gravity anomaly values ​​of all observation points in the artificial terrain area based on the planar coordinates of all observation points in the artificial terrain area, determine the gravity anomaly variation characteristics of the artificial terrain area, and complete the missing gravity anomaly data of the artificial terrain area.

[0169] Example 3

[0170] This invention also provides an electronic device, including a memory and a processor. The memory stores a computer program, and the processor runs the computer program to enable the electronic device to perform the gravity anomaly determination method for artificial terrain areas described in Embodiment 1. The electronic device may be a server.

[0171] In addition, the present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for determining gravity anomalies in artificial terrain areas in Embodiment 1.

[0172] In summary, this invention provides a method, system, equipment, and medium for determining gravity anomalies in artificial terrain areas. It involves using a GPS measuring instrument to determine the coordinates of ground control point markers; using a drone to measure terrain data within the artificial terrain area; correcting the terrain data based on the coordinates of the control point markers to obtain the corrected terrain data; using a GPS measuring instrument to determine the location information (including planar coordinates, geographical latitude and longitude, and altitude) of all observation points within the artificial terrain area; using a gravimeter to measure the gravity field data (including observed gravity values) of all observation points within the artificial terrain area; calculating the gravity anomaly of all observation points within the artificial terrain area based on the corrected terrain data, the location information of all observation points, and the gravity field data of all observation points; and gridding the gravity anomaly of all observation points within the artificial terrain area based on their planar coordinates to determine the characteristics of gravity anomaly changes and obtain missing data for the artificial terrain area. This solves the problem of missing gravity anomaly data in artificial terrain areas, which is beneficial for geological structure research and mineral exploration.

[0173] Compared with the prior art, the present invention has the following advantages:

[0174] 1. This invention provides complete gravity anomaly data for research in basic geology, mineral exploration breakthroughs and other related fields.

[0175] 2. This invention improves the rationality, accuracy, and reliability of inference and interpretation of geological, mineral resource, or other detection targets.

[0176] 3. This invention improves the accuracy and efficiency of gravity anomaly observation.

[0177] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.

[0178] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for determining gravity anomalies in artificial terrain areas, characterized in that, The method for determining gravity anomalies in artificial terrain areas is applied to a gravity anomaly observation device for artificial terrain areas. The observation device includes: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers, with the ground control point markers placed within the artificial terrain area. The method for determining gravity anomalies in artificial terrain areas includes: Using the GPS measuring instrument, the center coordinates of the ground control point sign were determined; The drone was used to measure image data within the artificial terrain area; The image data is corrected based on the center location coordinates to obtain terrain data; Using the GPS measuring instrument, the plane coordinates, geographical latitude, and altitude of all observation points within the artificial terrain area were determined; Using the gravimeter, the observed gravity values ​​at all observation points within the artificial terrain area were measured; Based on the plane coordinates, geographical latitude, altitude, observed gravity values, and topographic data of all observation points within the artificial terrain area, calculate the gravity anomaly values ​​of all observation points within the artificial terrain area. The gravity anomaly values ​​of all observation points within the artificial terrain area are gridded based on the planar coordinates of all observation points within the artificial terrain area to determine the gravity anomaly variation characteristics of the artificial terrain area and to supplement the missing gravity anomaly data of the artificial terrain area.

2. The method for determining gravity anomalies in artificial terrain areas according to claim 1, characterized in that, Based on the plane coordinates, geographical latitude, elevation, observed gravity values, and topographic data of all observation points within the artificial terrain area, the gravity anomaly values ​​of all observation points within the artificial terrain area are calculated, specifically including: Calculate the Bouguer correction values ​​for all observation points within the artificial terrain area based on their geographical latitude and elevation. Based on the plane coordinates, elevation, and topographic data of all observation points within the artificial terrain area, calculate the terrain correction values ​​for all observation points within the artificial terrain area. Calculate the normal gravity values ​​of all observation points within the artificial terrain area based on their geographical latitude. Based on the observed gravity values, normal gravity values, Bouguer correction values, and terrain correction values ​​of all observation points within the artificial terrain area, calculate the gravity anomaly values ​​of all observation points within the artificial terrain area.

3. The method for determining gravity anomalies in artificial terrain areas according to claim 2, characterized in that, Based on the geographical latitude and elevation of all observation points within the artificial terrain area, calculate the Bouguer correction value for all observation points within the artificial terrain area using the following formula: Where: δ gB This is the Bouguer correction value, in units of 10. -5 m / s 2 ; ρ1 is the geographical latitude of the observation point; h is the altitude of the observation point; ρ1 is the density of the intermediate layer; a is the correction radius of the intermediate layer of the circular domain.

4. The method for determining gravity anomalies in artificial terrain areas according to claim 2, characterized in that, Based on the plane coordinates, elevation, and topographic data of all observation points within the artificial terrain area, the topographic correction value for all observation points within the artificial terrain area is calculated using the following formula: Among them, g gT ρ is the topographic correction value; G is the gravitational constant; ρ is the average density of the Earth's crust; l is the integral grid spacing; C ij r is the integration constant, specifically the trapezoidal coefficient. ij h is the distance between the integration node (i, j) and the observation point; ij The elevation difference between the integration node (i, j) and the observation point is denoted as i; i is the row number of the integration node; j is the column number of the integration node; the integration node is determined based on the terrain data; the distance between the integration node and the observation point is determined based on the plane coordinates of the observation point; the elevation difference between the integration node and the observation point is determined based on the altitude of the observation point.

5. The method for determining gravity anomalies in artificial terrain areas according to claim 2, characterized in that, Based on the geographical latitude of all observation points within the artificial terrain area, the normal gravity values ​​of all observation points within the artificial terrain area are calculated using the following formula: Where g0 is the normal gravity value, in units of 10. -5 m / s 2 ; The latitude of the observation point.

6. The method for determining gravity anomalies in artificial terrain areas according to claim 2, characterized in that, Based on the observed gravity values, normal gravity values, Bouguer correction values, and topographic correction values ​​of all observation points within the artificial terrain area, the gravity anomaly values ​​of all observation points within the artificial terrain area are calculated using the following formula: Δg B =g-g0+δ gB +g gT Where: Δg B δ represents the gravity anomaly; g represents the observed gravity value; g0 represents the normal gravity value; δ gB g is the Bouguer correction value. gT This is the terrain correction value.

7. The method for determining gravity anomalies in artificial terrain areas according to claim 2, characterized in that, The image data is corrected based on the center location coordinates to obtain terrain data, specifically including: Based on the center position coordinates, aerial triangulation is performed on the image data to obtain a three-dimensional elevation data volume; The three-dimensional elevation data volume is processed by data synthesis and format conversion to obtain terrain data.

8. A system for determining gravity anomalies in artificial terrain areas, characterized in that, The gravity anomaly determination system for artificial terrain areas is applied to a gravity anomaly observation device for artificial terrain areas. The artificial terrain area gravity anomaly observation device includes: a drone, a gravimeter, a GPS measuring instrument, and ground control point markers, with the ground control point markers placed within the artificial terrain area. The gravity anomaly determination system for artificial terrain areas includes: The center position coordinate acquisition module is used to determine the center position coordinates of the ground control point sign using the GPS measuring instrument; The image data acquisition module is used to measure image data within the artificial terrain area using the UAV; The terrain data acquisition module is used to correct the image data according to the center location coordinates to obtain terrain data; The observation point information acquisition module is used to determine the plane coordinates, geographical latitude, and altitude of all observation points within the artificial terrain area using the GPS measuring instrument. The gravity value acquisition module is used to measure the gravity values ​​of all observation points within the artificial terrain area using the gravimeter. The gravity anomaly determination module is used to calculate the gravity anomaly value of all observation points within the artificial terrain area based on the plane coordinates, geographical latitude, altitude, observed gravity value, and the terrain data of all observation points within the artificial terrain area. The gridding module is used to grid the gravity anomaly values ​​of all observation points within the artificial terrain area based on the planar coordinates of all observation points within the artificial terrain area, determine the gravity anomaly variation characteristics of the artificial terrain area, and supplement the missing gravity anomaly data of the artificial terrain area.

9. An electronic device, characterized in that, The device includes a memory and a processor, the memory being used to store a computer program, and the processor running the computer program to cause the electronic device to perform the gravity anomaly determination method for artificial terrain areas as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, It stores a computer program that, when executed by a processor, implements the method for determining gravity anomalies in artificial terrain areas as described in any one of claims 1 to 7.

Images