Aerial-and-ground data combined gravity conversion method and system

Abstract

An aerial-and-ground data combined gravity conversion method includes the following steps: calculate the first estimated ground gravity by the Runge-Kutta format 1, and calculate the first error between the first estimated ground gravity and the measured ground gravity; calculate the second estimated ground gravity by the Runge-Kutta format 2, and calculate the second error between the second estimated ground gravity and the measured ground gravity; and select the smaller one from the first and second errors, use the corresponding Runge-Kutta format as the Runge-Kutta format for gravity conversion, and finish the gravity data conversion using the mentioned Runge-Kutta format.

Claims

1. A method of mineral resource exploration using an aerial-and-ground data combined gravity conversion method, comprising the following steps: S1: calculating, by one or more processors, a first estimated ground gravity by inputting a first measured aerial gravity data into a first Runge-Kutta format, and calculating a first error between the first estimated ground gravity and a measured ground gravity; S2: calculating, by the one or more processors, a second estimated ground gravity by inputting the first measured aerial gravity data into a second Runge-Kutta format, and calculating a second error between the second estimated ground gravity and the measured ground gravity; and S3: selecting, by the one or more processors, a smaller error from the first error and the second error, and converting a second measured aerial gravity data into a third estimated ground gravity by using the Runge-Kutta format to the smaller error; S4: estimating, by the one or more processors, a spatial and physical property distribution of underground elements by utilizing the second measured aerial gravity data and the third estimated ground gravity; and S5: carrying out mineral resource exploration based on the spatial and physical property distribution.

2. The method according to claim 1, wherein the step S1 comprises the following steps: S11: calculating, by the one or more processors, an estimated aerial gravity u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h) by the aerial gravity data u(x,y,z.sub.0) and calculating a first estimated aerial gradient k.sub.1 by using the first Runge-Kutta format, wherein x, y, and z.sub.0 represents a three-dimensional spatial location of the first measured aerial gravity data, with z.sub.0 representing the initial altitude, a.sub.2 is a constant, and h is a positive number indicating the conversion altitude; S12: calculating, by the one or more processors, a second estimated aerial gravity gradient k.sub.2 by an integrated second vertical derivative; S13: calculating, by the one or more processors, a third estimated gravity gradient k.sub.3=u.sub.z(k.sub.3.sub.)(x,y,z.sub.0+a.sub.3h) and a fourth estimated gravity gradient k.sub.4=u.sub.z(k.sub.4.sub.)(x,y,z.sub.0+a.sub.4h), wherein a.sub.3 and a.sub.4 are constants; S14: calculating, by the one or more processors, the first estimated ground gravity u.sub.c(x,y,z.sub.0+h) by the first Runge-Kutta format; S15: calculating, by the one or more processors, the first error ε by a root-mean-square error method.

3. The method according to claim 2, wherein the step S11 comprises the following steps: calculating, by the one or more processors, the estimated aerial gravity according to the following formula:
u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h)=u(x,y,z.sub.0)+b.sub.21k.sub.1 wherein a.sub.2 and b.sub.21 are constants, u(x,y,z.sub.0) is the arial gravity data, k.sub.1=u.sub.z(x,y,z.sub.0) is the first aerial gravity gradient data, and u.sub.(k.sub.2)(x,y,z.sub.0+a.sub.2h) is the estimated aerial gravity.

4. The aerial and ground data combined gravity conversion method according to claim 2, wherein the first Runge-Kutta format is: $\begin{array}{cc}\{\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/2h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+1/2h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/6k_1+2/6k_2+2/6k_3+1/6k_4\right)\end{array}& \end{array}$ wherein k.sub.1, k.sub.2 and k.sub.3 are the first estimated gravity gradient, the second estimated gravity gradient and the third estimated gravity gradient for three air positions, k.sub.4 is the fourth estimated gravity gradient for a ground position, u(x,y,z.sub.0) is the aerial gravity data, and u.sub.c(x,y,z.sub.0+h) is the first estimated ground gravity, wherein u.sub.c denotes ground gravity, u.sub.z denotes aerial gravity, and u denotes gravity in general.

5. The method according to claim 1, wherein the second Runge-Kutta format is: $\{\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/3h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+2/3h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c^{\prime }\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/8k_1+3/8k_2+3/8k_3+1/8k_4\right)\end{array}$ wherein k.sub.1, k.sub.2 and k.sub.3 are the first estimated gravity gradient, the second estimated gravity gradient and the third estimated gravity gradient for different air positions, k.sub.4 is the estimated gravity gradient for a ground position, u(x,y,z.sub.0) is the aerial gravity data, and u.sub.c′(x,y,z.sub.0+h) is the second estimated ground gravity, wherein u.sub.c denotes ground gravity, u.sub.z denotes aerial gravity, and u denotes gravity in general.

6. An aerial-and-ground data combined gravity conversion system, comprising: one or more processors; a storage device, for storing one or more programs, when the one or more programs are executed by the one or more processors, the one or more processors execute the method according to claim 1.

7. The aerial-and-ground data combined gravity conversion system according to claim 6, wherein the step S1 comprises the following steps: S11: calculating, by the one or more processors, an estimated aerial gravity u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h) by the aerial gravity data u(x,y,z.sub.0) and calculating a first estimated aerial gradient k.sub.1 by using the first Runge-Kutta format, wherein x, y, and z.sub.0 represents a three-dimensional spatial location of the first measured aerial gravity data, with z.sub.0 representing the initial altitude, a.sub.2 is a constant, and h is a positive number indicating the conversion altitude; S12: calculating, by the one or more processors, a second estimated aerial gravity gradient k.sub.2 by an integrated second vertical derivative; S13: calculating, by the one or more processors, a third estimated gravity gradient k.sub.3=u.sub.z(k.sub.3.sub.)(x,y,z.sub.0+a.sub.3h) and a fourth estimated gravity gradient k.sub.4=u.sub.z(k.sub.4.sub.)(x,y,z.sub.0+a.sub.4h), wherein a.sub.3 and a.sub.4 are constants; S14: calculating, by the one or more processors, the first estimated ground gravity u.sub.c(x,y,z.sub.0+h) by the first Runge-Kutta format; S15: calculating, by the one or more processors, the first error ε by a root-mean-square error method.

8. The aerial-and-ground data combined gravity conversion system according to claim 7, wherein the step S11 comprises the following steps: calculating, by the one or more processors, the estimated aerial gravity according to the following formula:
u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h)=u(x,y,z.sub.0)+b.sub.21k.sub.1 wherein a.sub.2 and b.sub.21 are constants, u(x,y,z.sub.0) is the aerial gravity data, k.sub.1=u.sub.z(x,y,z.sub.0) is the first aerial gravity gradient data, and u.sub.(k.sub.2)(x,y,z.sub.0+a.sub.2h) is the estimated aerial gravity.

9. The aerial-and-ground data combined gravity conversion system according to claim 7, wherein the first Runge-Kutta format is: $\{\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/2h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+1/2h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/6k_1+2/6k_2+2/6k_3+1/6k_4\right)\end{array}$ wherein k.sub.1, k.sub.2 and k.sub.3 are the first estimated gravity gradient, the second estimated gravity gradient and the third estimated gravity gradient for three air positions, k.sub.4 is the fourth estimated gravity gradient for a ground position, u(x,y,z.sub.0) is the aerial gravity data, and u.sub.c(x,y,z.sub.0+h) is the first estimated ground gravity, wherein u.sub.c denotes ground gravity, u.sub.z denotes aerial gravity, and u denotes gravity in general.

10. The aerial-and-ground data combined gravity conversion system according to claim 6, wherein the second Runge-Kutta format is: $\{\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/3h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+2/3h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c^{\prime }\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/8k_1+3/8k_2+3/8k_3+1/8k_4\right)\end{array}$ wherein k.sub.1, k.sub.2 and k.sub.3 are the first estimated gravity gradient, the second estimated gravity gradient and the third estimated gravity gradient for different air positions, k.sub.4 is the estimated gravity gradient for a ground position, u(x,y,z.sub.0) is the aerial gravity data, and u.sub.c′(x,y,z.sub.0+h) is the second estimated ground gravity, wherein u.sub.c denotes ground gravity, u.sub.z denotes aerial gravity, and u denotes gravity in general.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) All the above is only an overview of the technical solutions for the invention. To better understand the technical means involved, the figures and specific embodiments are integrated to illustrate the invention in detail.

(2) FIG. 1 is a flow chart of the aerial-and-ground data combined gravity conversion method provided in the embodiments of the invention;

(3) FIG. 2 is a structural diagram of the aerial-and-ground data combined gravity conversion system provided in the embodiments of the invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(4) The preferred embodiments of the invention are described in conjunction with the figures. It is to be understood that the preferred embodiments herein are intended to describe and illustrate the invention only, but without any limit to the invention.

(5) The invention provides an aerial-and-ground data combined gravity conversion Runge-Kutta method. Compared with the conventional methods, the invention utilizes the fourth-order Runge-Kutta format to implement easy and straightforward computation. The fourth-order Runge-Kutta format hereof uses the vertical derivatives at four middle altitude points before and after the conversion for calculation, which improves the accuracy of the conversion result. In addition, the derivatives are calculated by the integrated second vertical derivative method, significantly improving the stability and accuracy of the conversion. The invention aims to provide an aerial-and-ground data combined gravity conversion Runge-Kutta method, to solve the problems in the background set forth above.

(6) FIG. 1 is a flow chart of the aerial-and-ground data combined gravity conversion method provided in the embodiments of the invention. As shown in FIG. 1, the aerial-and-ground data combined gravity conversion method provided in the embodiments of the invention comprises:

(7) S11, calculate the first estimated ground gravity by the Runge-Kutta format 1, and calculate the first error between the first estimated ground gravity and the measured ground gravity.

(8) S12, calculate the second estimated ground gravity by the Runge-Kutta format 2, and calculate the second error between the second estimated ground gravity and the measured ground gravity.

(9) S13, select the smaller one from the first and second errors, use the corresponding Runge-Kutta format as the Runge-Kutta format for gravity conversion, and finish the gravity data conversion using the mentioned Runge-Kutta format. The Runge-Kutta format mentioned herein may be either the Runge-Kutta format 1 or the Runge-Kutta format 2.

(10) In order to overcome the existing technical shortcomings and achieve the goals set forth above, the invention proposes an aerial-and-ground data combined gravity conversion Runge-Kutta method. The specific steps of its technical solution are as follows: 1. Calculate the estimated aerial gravity u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h) by the aerial gravity data u(x,y,z.sub.0) and gradient data k.sub.1 through fourth-order Runge-Kutta format 1; 2. Calculate the estimated gravity gradient k.sub.2 by the Integrated Second Vertical Derivative (ISVD) method; 3. Repeat steps 1-2 to calculate the estimated gravity gradients k.sub.3 and k.sub.4 in the air positions respectively; 4. Calculate the ground gravity data u.sub.c(x,y,z.sub.0+h) by the fourth-order Runge-Kutta format 1; 5. Compare and calculate the error ε between the gravity data u.sub.c(x,y,z.sub.0+h) and the measured data u.sub.e(x,y,z.sub.0+h) by the root-mean-square error method; 6. Repeat steps 1-5 by the fourth-order Runge-Kutta format 2; 7. Select the fourth-order Runge-Kutta format with the smaller error, and repeat steps 1-4 by the aerial gravity data u(x,y,z.sub.0) to obtain the aerial-and-ground data combined gravity Runge-Kutta conversion result u(x,y,z.sub.0+h).

(11) The invention is a breakthrough in the field of geophysical technology. While realizing the conversion of gravity data obtained from different observation surfaces, the technology enhances the accuracy of conversion, with practical significance to the accuracy improvement of gravity exploration as well as very important significance to the theoretical geological research and practical mineral exploration.

(12) In step 1, u(x,y,z.sub.0) and k.sub.1=u.sub.z(x,y,z.sub.0) are gravity data and gradient data at the altitude of z.sub.0, u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h) is the aerial gravity at a distance of z.sub.0+a.sub.2h above ground:
u.sub.(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h)=u(x,y,z.sub.0)+b.sub.21k.sub.1  (1)

(13) h is a positive number indicating the conversion altitude and a.sub.2 and b.sub.21 are constants obtained by the fourth-order Runge-Kutta format 1.

(14) In step 2, k.sub.2=u.sub.z(k.sub.2.sub.)(x,y,z.sub.0+a.sub.2h) is the estimated gravity gradient at an altitude of z.sub.0+a.sub.2h.

(15) The computational process of the Integrated Second Vertical Derivative (ISVD) method in step 2 above is as follows:

(16) Calculate the scalar potential {tilde over (v)}(x,y,z.sub.0+a.sub.2h) of ũ(x,y,z.sub.0+a.sub.2h) by the gravity data conversion method in the wavenumber domain:

(17) $\begin{array}{cc}\tilde{v}\left(x,y,z_0+a_{2}h\right)=F^{-1}\left[F\left[\tilde{u}\left(x,y,z_0+a_{2}h\right)\right]\frac{1}{2\pi \sqrt{k_x^2+k_y^2}}\right].& \left(2\right)\end{array}$

(18) where F is the Fourier transform, F.sup.−1 is the inverse Fourier transform, k.sub.x,k.sub.y is the coordinate in the wavenumber domain corresponding to the coordinate x,y in the spatial domain, ũ(k.sub.x,k.sub.x,z.sub.0) is the expression of u(x,y,z.sub.0) in the wavenumber domain, and

(19) $\frac{1}{2\pi \sqrt{k_x^2+k_y^2}}$
is the vertical integral operator in the wavenumber domain.

(20) The wavenumber domain mentioned above is obtained by Fourier transform, wherein the Fourier transform and the inverse transform are calculated as follows:

(21) $\begin{array}{cc}\{\begin{array}{c}\tilde{u}\left(k_x,k_x,z_0\right)=F\left[u\left(x,y,z_0\right)\right]\\ u\left(x,y,z_0\right)=F^{-1}\left[\tilde{u}\left(k_x,k_x,z_0\right)\right]\end{array}.& \left(3\right)\end{array}$

(22) Then calculate the second derivative of the scalar potential {tilde over (v)}(x,y,z.sub.0+a.sub.2h) in the horizontal direction by the finite difference method:

(23) $\begin{array}{cc}\{\begin{array}{c}{\tilde{v}}_{\mathrm{xx}}\left(x,y,z_0+a_{2}h\right)=\left[\tilde{v}\left(x+\Delta x,y,z_0+a_{2}h\right)+\tilde{v}\left(x-\Delta x,y,z_0+a_{2}h\right)-2\tilde{v}\left(x,y,z_0+a_{2}h\right)\right]/\Delta x^2\\ {\tilde{v}}_{\mathrm{yy}}\left(x,y,z_0+a_{2}h\right)=\left[\tilde{v}\left(x+\Delta y,y,z_0+a_{2}h\right)+\tilde{v}\left(x-\Delta y,y,z_0+a_{2}h\right)-2\tilde{v}\left(x,y,z_0+a_{2}h\right)\right]/\Delta y^2\end{array}.& \left(4\right)\end{array}$

(24) where Δx,Δy is the sampling interval in the horizontal direction.

(25) At last, obtain the second derivative of the scalar potential {tilde over (v)}(x,y,z.sub.0+a.sub.2h) in the vertical direction by the Laplace's equation, that is, the first derivative ũ.sub.z(x,y,z.sub.0+a.sub.2h) of the gravity data ũ(x,y,z.sub.0+a.sub.2h) data in the vertical direction:

(26) $\begin{array}{cc}\begin{array}{c}{\tilde{u}}_z\left(x,y,z_0+a_{2}h\right)={\tilde{v}}_{\mathrm{zz}}\left(x,y,z_0+a_{2}h\right)\\ =-\left[{\tilde{v}}_{\mathrm{xx}}\left(x,y,z_0+a_{2}h\right)+{\tilde{v}}_{\mathrm{yy}}\left(x,y,z_0+a_{2}h\right)\right].\end{array}& \left(5\right)\end{array}$

(27) Obtain the following by the inverse Fourier transform,

(28) $\begin{array}{cc}\begin{array}{c}{u}_{z\left(k_2\right)}\left(x,y,z_0+a_{2}h\right)=u_z\left(x,y,z_0+a_{2}h\right)\\ =F^{-1}\left[{\tilde{u}}_z\left(x,y,z_0+a_{2}h\right)\right].\end{array}& \left(6\right)\end{array}$

(29) In step 3, k.sub.3=u.sub.z(k.sub.3.sub.)(x,y,z.sub.0+a.sub.3h) and k.sub.4=u.sub.z(k.sub.4.sub.)(x,y,z.sub.0+a.sub.4h) are the estimated gravity gradients calculated by steps 1-2 for the altitudes z.sub.0+a.sub.3h and z.sub.0+a.sub.4h.

(30) The fourth-order Runge-Kutta format 1 in step 4 is:

(31) $\begin{array}{cc}\{\text{⁠}\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/2h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+1/2h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/6k_1+2/6k_2+2/6k_3+1/6k_4\right)\end{array}& \left(7\right)\end{array}$

(32) u.sub.c(x,y,z.sub.0+h) is calculated using the format (7).

(33) The root mean square error in step 5 is

(34) 0$\epsilon =\sqrt{\frac{1}{N}\underset{i=1,j=1}{\overset{N}{.Math.}}\left(u_c\left(x_i,y_j,z_0+h\right)-u_e\left(x_i,y_j,z_0+h\right)\right)}.$

(35) The fourth-order Runge-Kutta format 2 in step 6 is:

(36) $\begin{array}{cc}\{\text{⁠}\begin{array}{c}{k}_1=u_z\left(x,y,z_0\right)\\ k_2=u_{z\left(k_2\right)}\left(x,y,z_0+1/3h\right)\\ k_3=u_{z\left(k_3\right)}\left(x,y,z_0+2/3h\right)\\ k_4=u_{z\left(k_4\right)}\left(x,y,z_0+h\right)\\ u_c^{\prime }\left(x,y,z_0+h\right)=u\left(x,y,z_0\right)+\left(1/8k_1+3/8k_2+3/8k_3+1/8k_4\right)\end{array}& \left(8\right)\end{array}$

(37) FIG. 2 is a structural diagram of the aerial-and-ground data combined gravity conversion system of the invention. As shown in FIG. 2, the aerial-and-ground data combined gravity conversion system comprises: a Central Processing Unit (CPU) 201 that can execute appropriate actions and processes in accordance with the programs stored in the read-only memory (ROM) or those loaded into the random access memory (RAM) 203 from a memory 208. A variety of programs and data required for system operation are also stored in the RAM 203. The CPU 201, ROM 202, and RAM 203 are interconnected through a bus 204. The input/output (I/O) interface 205 is also coupled to the bus 204.

(38) The following components are connected to the I/O interface 205: an input section 206 comprising a keyboard and a mouse; an output section 207 comprising a cathode ray tube (CRT), a liquid crystal display (LCD) and a loudspeaker; a storage section 208 comprising a hard disk; and a communication section 209 comprising a LAN card, a modem and other network interfaces. The communication section 209 performs communication processing via a network, such as the Internet. A driver 210 is also connected to the I/O interface 205 as needed. A detachable medium 211, which may be a disk, CD a magneto-optical disk or a semiconductor memory, is mounted on the driver 210 as needed, so that the computer programs read from the medium can be installed into the storage section 208 as required.

(39) In particular, according to the embodiments of the invention, the processes described in the flow chart may be implemented as computer software programs. For example, a computer program product is involved in the embodiment of the invention, which comprises a computer program carried on a computer readable medium; the computer program contains the program codes for executing the method as shown in the flow chart. In such embodiment, the computer program can be downloaded from the network and installed via the communication section 209, and/or be installed through the detachable medium 211. When the computer program is executed by CPU 201, the functions defined in the invention are available. It is important to note that the computer-readable medium involved in the invention may be a computer readable signal medium, or a computer-readable storage medium, or any combination of both. The computer-readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus or device, or any combination of them. The specific examples of the computer readable storage medium may include, but are not limited to, an electric connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read only memory (ROM), an erasable programmable read only memory (EPROM or flash memory), an optical fiber, a portable compact disk read only memory (CD-ROM), an optical storage device, a magnetic storage device, or any appropriate combination of them. In the invention, the computer readable storage medium may be any tangible medium that can carry or store a program available for independent or combined application by the instruction execution system, apparatus or device. The computer readable signal medium in the invention may include some data signals communicated in the baseband or as part of a carrier, which carries the computer readable program codes. The communicated data signals can take a variety of forms including, but not limited to, electromagnetic signals, optical signals, or any appropriate combination of them. The computer readable signal medium can also be any other computer readable medium than a computer readable storage medium, which can transmit, communicate or transport the programs available for independent or combined application by the instruction execution system, apparatus or device. The program codes carried on the computer readable medium can be transmitted by any suitable medium, including but not limited to wireless, wire, fiber optic cable, RF, or any appropriate combination of them.

(40) The flow chart and block diagram illustrate the architecture, functionality and operation potentially accomplished by the systems, methods and computer program products involved in the embodiments of the invention. In this regard, each block of the flow chart or block diagram can represent a module, a program segment, or a section of a code that includes one or more executable instructions for implementing the specified logic functions. It should be noted that in some alternative embodiments, the functions indicated in the blocks may also be performed in a different order than that illustrated in the diagrams. For example, two consecutive blocks may be executed substantially in parallel, which may sometimes be executed in the reverse order as required by the functionality involved. Furthermore, every block in the block diagram and/or flow chart, as well as the combinations of the blocks in the block diagram and/or flow chart can be implemented by a dedicated hardware-based system performing the specified function or operation, or by a combination of dedicated hardware and computer instructions.

(41) The units involved in the embodiments of the invention may be implemented by either software or hardware.

(42) The invention provides an aerial-and-ground data combined gravity conversion method. Compared with the conventional methods, the invention utilizes the fourth-order Runge-Kutta format to implement easy and straightforward computation. The fourth-order Runge-Kutta format hereof uses the vertical derivatives at four middle altitude points before and after the conversion for calculation, which improves the accuracy of the conversion result. In addition, the derivatives are calculated by the integrated second vertical derivative method, significantly improving the stability and accuracy of the conversion. The fourth-order Runge-Kutta format has different expression formulas, making it possible to obtain the conversion results of different ground gravity data. The accuracy of the final conversion results can be optimized by the comparing different converted ground gravity data with the measured ground gravity data of some known regions and selecting the fourth-order Runge-Kutta format with smaller calculation errors.

(43) All the above is nothing but a preferred embodiment of the invention, and by no means shall it impose any limitation against any form of the invention. Any person skilled in the art who makes simple modifications, equivalent changes or modifications to the disclosed technical contents will be held accountable within the protection scope of the invention.