Gravity anomaly simulation method, device, equipment, medium and product

Abstract

The invention discloses a gravity anomaly simulation method and device, equipment, a medium and a product, and relates to the technical field of geophysical exploration, and the method comprises the steps: constructing a two-dimensional body model of a target exploration region, carrying out the grid division of the two-dimensional body model, carrying out the abnormal body residual density assignment of each grid node, and obtaining a geologic body density model; determining a Poisson equation satisfied by a wave number domain gravity potential based on the geologic body density model; solving the Poisson equation satisfied by the wave number domain gravitational potential by adopting one-dimensional non-uniform fast Fourier transform to obtain the wave number domain gravitational potential; determining a wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential; and performing one-dimensional non-uniform fast Fourier inversion on the wavenumber domain gravity anomaly field to obtain a spatial domain gravity anomaly field of the target exploration area. According to the method, the calculation efficiency can be improved by utilizing the non-uniform sampling Fourier transform.

Description

Technical Field

[0001] The present application relates to the field of geophysical exploration technology, and in particular to a gravity anomaly simulation method, device, equipment, medium and product. Background Art

[0002] Gravity exploration is a geophysical exploration method that uses the density difference of underground media as the physical property basis to study the deep structure of the earth, detect hidden rock bodies or rock layers and faults, and search for mineral resources. The actual observed gravity data is directly related to the underground density distribution. Therefore, underground density anomalies can be detected by inverting the underground density. With the advancement of gravity exploration technology, the refined exploration of large-scale and complex geological bodies has gradually attracted the attention of scholars, and the traditional gravity anomaly forward modeling method has become more and more time-consuming. Therefore, efficient gravity anomaly forward modeling becomes more and more important. Summary of the invention

[0003] The purpose of this application is to provide a gravity anomaly simulation method, device, equipment, medium and product that can improve computing efficiency.

[0004] To achieve the above objectives, this application provides the following solutions:

[0005] In a first aspect, the present application provides a gravity anomaly simulation method, comprising:

[0006] Constructing a two-dimensional body model of the target exploration area, gridding the two-dimensional body model, assigning anomaly body residual density to each grid node, and obtaining a geological body density model;

[0007] Determine the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model;

[0008] A one-dimensional non-uniform fast Fourier transform is used to solve the Poisson equation satisfied by the wave number domain gravity potential to obtain the wave number domain gravity potential;

[0009] determining a wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential;

[0010] A one-dimensional non-uniform inverse fast Fourier transform is performed on the wavenumber domain gravity anomaly field to obtain the space domain gravity anomaly field of the target exploration area.

[0011] Optionally, determining the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model specifically includes:

[0012] Based on the geological body density model, a Poisson equation satisfied by the spatial domain gravity potential of the target exploration area is obtained;

[0013] A one-dimensional fast Fourier transform is performed on the Poisson equation satisfied by the gravity potential in the space domain to obtain the Poisson equation satisfied by the gravity potential in the wavenumber domain.

[0014] Optionally, the Poisson equation satisfied by the spatial domain gravity potential is expressed as: 2 U(x,z)=-4πγρ(x,z);

[0015] Among them, 2 represents the Laplace operator, U(x,z) represents the gravity potential in the spatial domain, ρ(x,z) represents the residual density of the abnormal body in the spatial domain, (x,z) is the spatial coordinate of any point in the geological body density model, x is the x-axis coordinate, z is the z-axis coordinate, and γ is the gravitational constant.

[0016] Optionally, the Poisson equation satisfied by the wavenumber domain gravity potential is expressed as:

[0017] in, represents the gravity potential in the wavenumber domain, It represents the residual density of the anomaly body in the wavenumber domain, and k is the wavenumber.

[0018] Optionally, using a one-dimensional non-uniform fast Fourier transform to solve the Poisson equation satisfied by the wave number domain gravity potential to obtain the wave number domain gravity potential specifically includes:

[0019] The one-dimensional non-uniform fast Fourier transform is realized by combining an interpolation algorithm with an express Fourier transform, and the Poisson equation satisfied by the wavenumber domain gravity potential is solved by the one-dimensional non-uniform fast Fourier transform to obtain the wavenumber domain gravity potential.

[0020] Optionally, determining the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential specifically includes:

[0021] According to the formula Determine the gravity anomaly field in the wavenumber domain;

[0022] in, represents the horizontal component of the gravity anomaly field in the wavenumber domain, represents the vertical component of the gravity anomaly field in the wavenumber domain, i is the imaginary unit, k is the wavenumber, Represents the gravity potential in the wavenumber domain.

[0023] In a second aspect, the present application provides a gravity anomaly simulation device, the gravity anomaly simulation device comprising:

[0024] A geological body density model construction module is used to construct a two-dimensional body model of the target exploration area, grid the two-dimensional body model, assign an abnormal body residual density to each grid node, and obtain a geological body density model;

[0025] A module for determining the Poisson equation satisfied by the wavenumber domain gravity potential, used for determining the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model;

[0026] A wave number domain gravity potential solving module is used to solve the Poisson equation satisfied by the wave number domain gravity potential by using a one-dimensional non-uniform fast Fourier transform to obtain the wave number domain gravity potential;

[0027] A wavenumber domain gravity anomaly field determination module, used to determine the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential;

[0028] The spatial domain gravity anomaly field determination module is used to perform a one-dimensional non-uniform fast Fourier inverse transform on the wave number domain gravity anomaly field to obtain the spatial domain gravity anomaly field of the target exploration area.

[0029] In a third aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any one of the gravity anomaly simulation methods described above.

[0030] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any one of the gravity anomaly simulation methods described above.

[0031] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of any one of the gravity anomaly simulation methods described above.

[0032] According to the specific embodiments provided in this application, this application discloses the following technical effects:

[0033] The present application provides a gravity anomaly simulation method, device, equipment, medium and product, which adopts a one-dimensional non-uniform fast Fourier transform to solve the Poisson equation satisfied by the wavenumber domain gravity potential to obtain the wavenumber domain gravity potential, and performs a one-dimensional non-uniform fast Fourier inverse transform on the wavenumber domain gravity anomaly field to obtain the space domain gravity anomaly field of the target exploration area. The observation points are selected by one-dimensional non-uniform fast Fourier transform and inverse transform, and the gravity anomaly forward calculation can be performed on the observation points arbitrarily distributed in the wavenumber domain, thereby reducing the calculation workload and improving the calculation efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0035] Figure 1 A schematic diagram of a flow chart of a gravity anomaly simulation method provided in one embodiment of the present application;

[0036] Figure 2 A schematic diagram of a two-dimensional body with a matrix cross section provided in an embodiment of the present application;

[0037] Figure 3 A schematic diagram of the effect of the horizontal component of the gravity anomaly field provided in one embodiment of the present application;

[0038] Figure 4 A schematic diagram of the effect of the vertical component of the gravity anomaly field provided in one embodiment of the present application;

[0039] Figure 5 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0040] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0041] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0042] This application provides a gravity anomaly simulation method, such as Figure 1 As shown, the gravity anomaly simulation method includes:

[0043] Step 101: construct a two-dimensional volume model of the target exploration area, grid the two-dimensional volume model, assign anomaly residual density to each grid node, and obtain a geological volume density model.

[0044] Step 102: Determine the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model.

[0045] Step 103: using a one-dimensional non-uniform fast Fourier transform to solve the Poisson equation satisfied by the wave number domain gravity potential, and obtain the wave number domain gravity potential.

[0046] Step 104: Determine the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential.

[0047] Step 105: Perform a one-dimensional non-uniform inverse fast Fourier transform on the wavenumber domain gravity anomaly field to obtain the space domain gravity anomaly field of the target exploration area.

[0048] In an exemplary embodiment, step 101 specifically includes: establishing a spatial model of the target exploration area and determining the starting position of the two-dimensional body model in the x and z directions; then discretizing the two-dimensional body model space along parallel to the x and z axes respectively, and the grid size can be arbitrarily divided according to the size and density changes of the exploration target; finally, assigning the residual density of the abnormal body to the grid nodes in the two-dimensional body model space to characterize the density model of the complex geological body with arbitrary density distribution.

[0049] For grid division, at least three points are required to pass above the anomaly during gravity exploration, so the grid interval is less than one-third of the anomaly size, and the smaller the grid interval, the better the approximation of the anomaly. In addition, the denser the grid is where the density of the underground area changes rapidly, the higher the accuracy of the simulation results.

[0050] In an exemplary embodiment, step 102 specifically includes: obtaining a Poisson's equation satisfied by the spatial domain gravity potential of the target exploration area based on the geological body density model; performing a one-dimensional fast Fourier transform on the Poisson's equation satisfied by the spatial domain gravity potential to obtain a Poisson's equation satisfied by the wavenumber domain gravity potential.

[0051] The Poisson equation satisfied by the gravity potential in the space domain is expressed as: 2 U(x,z)=-4πγρ(x,z)(1)

[0052] Among them, 2 represents the Laplace operator, U(x,z) represents the gravity potential in the spatial domain, ρ(x,z) represents the residual density of the abnormal body in the spatial domain, (x,z) is the spatial coordinate of any point in the geological body density model, x is the x-axis coordinate, z is the z-axis coordinate, and γ is the gravitational constant.

[0053] A one-dimensional fast Fourier transform is performed on formula (1) along the horizontal direction (x-axis direction) to transform the two-dimensional Poisson's equation into a one-dimensional partial differential equation problem that is independent of each other at different wave numbers, that is, the Poisson's equation satisfied by the gravity potential in the wavenumber domain.

[0054] The Poisson equation satisfied by the gravity potential in the wavenumber domain is expressed as:

[0055] in, represents the gravity potential in the wavenumber domain, It represents the residual density of the anomaly body in the wavenumber domain, and k is the wavenumber.

[0056] In an exemplary embodiment, step 103 specifically includes: using an interpolation algorithm combined with a fast Fourier transform to implement the one-dimensional non-uniform fast Fourier transform, using a one-dimensional non-uniform fast Fourier transform to solve the Poisson equation satisfied by the wavenumber domain gravity potential to obtain the wavenumber domain gravity potential.

[0057] For the wave number domain gravity potential, that is, the solution of formula (2), the accuracy and efficiency of Fourier transform directly affect the efficiency and accuracy of the final calculated gravity anomaly. This application uses a one-dimensional non-uniform fast Fourier transform to solve formula (2).

[0058] The one-dimensional non-uniform fast Fourier transform used in this application is:

[0059]

[0060] Among them, i is the imaginary unit, α n For a given discrete sampling point x n ∈[-N / 2,N / 2] corresponds to the value of the sampling point, that is, the residual density value of the abnormal body, F(α) j Each discrete sampling point {α n}, N represents the total number of sampling points.

[0061] The implementation of the non-uniform fast Fourier transform algorithm mainly includes the following three steps:

[0062] (1) Using the adjacent x n The non-uniform Fourier transform basis is approximated by the Fourier transform basis of q uniform sampling points

[0063]

[0064] Where: m represents the oversampling factor, ω l represents the weight factor, s j is the precision factor, [mx n ] indicates that mx n Round up.

[0065] (2) Calculate the Fourier transform coefficients of the new Fourier transform basis:

[0066]

[0067] Among them, η p are Fourier transform coefficients.

[0068] (3) Using uniform fast Fourier transform, we can get:

[0069]

[0070] Among them, F j The values ​​of each discrete sample point calculated for the jth sample point.

[0071] According to formulas (3), (4), (5) and (6), it can be seen that the non-uniform sampling fast Fourier transform is a combination of the standard fast Fourier transform and local interpolation to realize the non-uniform fast Fourier transform algorithm. The essence of the non-uniform Fourier transform of this application is to interpolate non-uniform points to uniform points through sampling, and then use uniform fast Fourier transform for calculation.

[0072] Determining the wave number domain gravity anomaly field according to the wave number domain gravity potential specifically includes: According to the simple characteristic of derivative in the wave number domain, a calculation formula of the wave number domain gravity anomaly field can be given, that is,

[0073]

[0074] in, represents the horizontal component of the gravity anomaly field in the wavenumber domain, represents the vertical component of the gravity anomaly field in the wavenumber domain, i is the imaginary unit, k is the wavenumber, Represents the gravity potential in the wavenumber domain.

[0075] This application uses a one-dimensional non-uniform Fourier transform to transform the two-dimensional partial differential equation satisfied by the gravity potential into one-dimensional partial differential equations independent of different wave numbers. Solving the one-dimensional partial differential equations corresponding to different wave numbers can obtain the gravity potential in the wave number domain The wave number here is the wave number selected in the wave number domain, which uses the non-uniform Fourier transform.

[0076] In an exemplary embodiment, during the one-dimensional non-uniform inverse fast Fourier transform of the wavenumber domain gravity anomaly field, the selection of wavenumbers directly affects the accuracy of the Fourier transform and the final calculation effect. If too many wavenumbers are selected, the calculation efficiency is low but the accuracy is high, while if too few wavenumbers are selected, the efficiency is high but the accuracy is low. Therefore, the wavenumbers are selected at equal intervals on the logarithmic axis from 0 to the cutoff frequency, so that the unity of accuracy and efficiency is achieved with the least wavenumber selection. This wavenumber selection method is introduced into Formula 6 to perform a one-dimensional non-uniform fast Fourier transform on the wavenumber domain gravity anomaly field in Formula 7 to obtain the gravity anomaly field, and the calculation result is output. The cutoff frequency is calculated by the Nyquist sampling theorem, that is, the cutoff frequency k max =π / Δx min , Δx min Indicates the minimum sampling interval in the x direction.

[0077] In an exemplary embodiment, a gravity anomaly simulation method of the present application is verified.

[0078] The target area has a two-dimensional volume with a cross section of a matrix such as Figure 2 As shown in the figure, the target area ranges from -400m to 400m in the x direction and from 0m to 400m in the z direction. The horizontal and vertical small rectangle sizes are both 2m and 1m. The number of nodes for dividing the two-dimensional model area is 401×401. The anomaly body is buried at a depth of 200m and has a size of 100×80m. 2 , residual density is 1000kg / m 3 , the calculated survey line is the gravity anomaly on the horizontal ground (the wave number is selected as 401).

[0079] The present application method is implemented by programming in Fortran language, and the configuration of the personal computer used to run the program is: CPU-Inter Core i7-8700, main frequency is 3.2GHz, and running memory is 8.00GB. Figure 3 Part a) is the result diagram of the horizontal component of the gravity anomaly field calculated by the method of this application and the analytical solution. Figure 3 In part a), it can be seen that the analytical solution agrees well with the numerical solution. Figure 3 Part b) is a relative error diagram between the analytical solution and the numerical solution of the horizontal component of the gravity field. It can be seen that the calculation accuracy of the method of the present application is relatively high, and the relative error is less than 0.04%. Figure 4 Part a) is the result diagram of the vertical component of the gravity anomaly field calculated by the method of this application and the analytical solution. Figure 4 In part a), it can be seen that the analytical solution agrees well with the numerical solution. Figure 4 Part b) is the relative error diagram between the analytical solution and the numerical solution of the vertical component of the gravity field. It can be seen that this method has high calculation accuracy and the relative error is less than 0.04%.

[0080] The present application can realize the forward modeling of gravity anomalies at observation points arbitrarily distributed in the wavenumber domain without the aid of any interpolation method. The gravity forward modeling process in the wavenumber domain is affected by the boundary truncation effect of the fast Fourier transform. The Fourier edge expansion method and Gauss-FFT are usually used to reduce the error caused by the boundary truncation effect of the fast Fourier transform by increasing the number of Gaussian points, while increasing the amount of calculation and memory requirements. The non-uniform Fourier transform of the present application can be sampled according to the spectral transformation, and the problem of suppressing the boundary truncation effect can be solved with the least wavenumber selection. When performing large-scale numerical simulation of gravity anomalies, this method not only greatly improves the calculation efficiency, but also takes up less memory during calculation, and is easy to parallelize.

[0081] Based on the same inventive concept, the embodiment of the present application also provides a gravity anomaly simulation device for implementing the gravity anomaly simulation method involved above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations of one or more gravity anomaly simulation device embodiments provided below can refer to the limitations of the gravity anomaly simulation method above, and will not be repeated here.

[0082] In an exemplary embodiment, the present application provides a gravity anomaly simulation device comprising:

[0083] A geological body density model construction module is used to construct a two-dimensional body model of the target exploration area, grid the two-dimensional body model, assign an abnormal body residual density to each grid node, and obtain a geological body density model;

[0084] A module for determining the Poisson equation satisfied by the wavenumber domain gravity potential, used for determining the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model;

[0085] A wave number domain gravity potential solving module is used to solve the Poisson equation satisfied by the wave number domain gravity potential by using a one-dimensional non-uniform fast Fourier transform to obtain the wave number domain gravity potential;

[0086] A wavenumber domain gravity anomaly field determination module, used to determine the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential;

[0087] The spatial domain gravity anomaly field determination module is used to perform a one-dimensional non-uniform fast Fourier inverse transform on the wave number domain gravity anomaly field to obtain the spatial domain gravity anomaly field of the target exploration area.

[0088] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 5As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store gravity anomaly simulation data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a gravity anomaly simulation method is implemented.

[0089] Those skilled in the art will understand that Figure 5 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components. In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0090] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0091] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0092] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0093] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0094] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a data processing logic of a programmable logic device, etc., but is not limited thereto.

[0095] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0096] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A gravity anomaly simulation method, characterized in that: The gravity anomaly simulation method comprises: Constructing a two-dimensional body model of the target exploration area, gridding the two-dimensional body model, assigning anomaly body residual density to each grid node, and obtaining a geological body density model; Determine the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model; A one-dimensional non-uniform fast Fourier transform is used to solve the Poisson equation satisfied by the wave number domain gravity potential to obtain the wave number domain gravity potential; determining a wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential; A one-dimensional non-uniform inverse fast Fourier transform is performed on the wavenumber domain gravity anomaly field to obtain the space domain gravity anomaly field of the target exploration area.

2. The gravity anomaly simulation method according to claim 1, characterized in that: Determining the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model specifically includes: Based on the geological body density model, a Poisson equation satisfied by the spatial domain gravity potential of the target exploration area is obtained; A one-dimensional fast Fourier transform is performed on the Poisson equation satisfied by the gravity potential in the space domain to obtain the Poisson equation satisfied by the gravity potential in the wavenumber domain.

3. The gravity anomaly simulation method according to claim 2, characterized in that: The Poisson equation satisfied by the spatial domain gravity potential is expressed as: in, represents the Laplace operator, U(x,z) represents the gravity potential in the space domain, ρ(x,z) represents the residual density of the abnormal body in the space domain, (x,z) is the spatial coordinate of any point in the geological body density model, x is the x-axis coordinate, z is the z-axis coordinate, and γ is the gravitational constant.

4. The gravity anomaly simulation method according to claim 3, characterized in that: The Poisson equation satisfied by the wavenumber domain gravity potential is expressed as: in, represents the gravity potential in the wavenumber domain, It represents the residual density of the anomaly body in the wavenumber domain, and k is the wavenumber.

5. The gravity anomaly simulation method according to claim 1, characterized in that: The Poisson equation satisfied by the wave number domain gravity potential is solved by using a one-dimensional non-uniform fast Fourier transform to obtain the wave number domain gravity potential, which specifically includes: The one-dimensional non-uniform fast Fourier transform is realized by combining an interpolation algorithm with an express Fourier transform, and the Poisson equation satisfied by the wavenumber domain gravity potential is solved by the one-dimensional non-uniform fast Fourier transform to obtain the wavenumber domain gravity potential.

6. The gravity anomaly simulation method according to claim 1, characterized in that: Determining the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential specifically includes: According to the formula Determine the gravity anomaly field in the wavenumber domain; in, represents the horizontal component of the gravity anomaly field in the wavenumber domain, represents the vertical component of the gravity anomaly field in the wavenumber domain, i is the imaginary unit, k is the wavenumber, Represents the gravity potential in the wavenumber domain.

7. A gravity anomaly simulation device, characterized in that: The gravity anomaly simulation device comprises: A geological body density model construction module is used to construct a two-dimensional body model of the target exploration area, grid the two-dimensional body model, assign an abnormal body residual density to each grid node, and obtain a geological body density model; A module for determining the Poisson equation satisfied by the wavenumber domain gravity potential, used for determining the Poisson equation satisfied by the wavenumber domain gravity potential based on the geological body density model; A wave number domain gravity potential solving module is used to solve the Poisson equation satisfied by the wave number domain gravity potential by using a one-dimensional non-uniform fast Fourier transform to obtain the wave number domain gravity potential; A wavenumber domain gravity anomaly field determination module, used to determine the wavenumber domain gravity anomaly field according to the wavenumber domain gravity potential; The spatial domain gravity anomaly field determination module is used to perform a one-dimensional non-uniform fast Fourier inverse transform on the wave number domain gravity anomaly field to obtain the spatial domain gravity anomaly field of the target exploration area.

8. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the gravity anomaly simulation method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the gravity anomaly simulation method described in any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the gravity anomaly simulation method described in any one of claims 1 to 6 is implemented.

Images