Three-dimensional gravity anomaly forward modeling method and device

Abstract

The invention discloses a three-dimensional gravity anomaly forward modeling method and device. The method comprises the following steps: constructing a gravity anomaly model comprising a simulation area; determining a gravity anomaly three-dimensional propagation formula of the simulation area based on the gravity anomaly model; performing cubic finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula to obtain a one-dimensional gravity wavenumber domain spectral integral of the gravity anomaly three-dimensional propagation formula; performing subdivision interpolation on the one-dimensional gravity wave number spectrum integral of the gravity anomaly three-dimensional propagation formula by adopting a quadratic finite element method, and then performing accumulation to obtain an accumulated one-dimensional gravity wave number spectrum integral; the accumulated one-dimensional gravity wave number spectrum integral is converted into a spatial domain through cubic finite element method Fourier transform, spatial domain gravity anomaly is obtained, interpolation is carried out on a spectrum coefficient part in a discrete unit through the cubic finite element method, spatial domain gravity anomaly can be deduced through a spectrum integral interpolation result in the unit, and the spatial domain gravity anomaly is obtained. The calculation precision of gravity anomaly three-dimensional forward calculation is improved, and the algorithm can be suitable for gravity anomaly three-dimensional forward calculation under complex conditions.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of gravity anomaly exploration, more particularly, to a three-dimensional gravity anomaly forward simulation method and device. BACKGROUND

[0002] In the numerical simulation of gravity anomaly frequency domain, the calculation precision of Fourier transform is one of the key factors affecting the simulation precision. Fourier transform is mathematically an integral of oscillatory function. At present, the methods for solving the integral of oscillatory function can be roughly divided into three categories: Filon method, Levin method and asymptotic expansion method. For Fourier transform, since the matrix of operator e -ikx is known, the Filon method is mostly used to realize Fourier transform, such as Clendenin (1966), Piessens and Poleunis (1971), Iserles (2004). The difference between the specific methods lies in that different forms of polynomial functions are used to approximate the kernel function in the subdivision interval. Clendenin (1966) uses linear function, Piessens and Poleunis (1971) use more complex Chebyshev polynomials, and Iserles (2004) makes in-depth analysis on the Filon method for realizing Fourier transform. Wu (2014, 2016) applies Gaussian integral to the integral of the kernel function, eliminating certain boundary effects.

[0003] The above-mentioned Fourier transform methods can be well applied to the calculation of low oscillation geophysical field, but with the increase of oscillation frequency, the calculation precision of these methods is reduced or even invalid, and they are not universal.

[0004] Therefore, a better solution is needed. SUMMARY

[0005] The present application provides a three-dimensional gravity anomaly forward simulation method and device to solve the technical problem of low calculation precision of three-dimensional forward calculation of gravity anomaly under complex conditions in the prior art. The method comprises:

[0006] constructing a gravity anomaly model comprising a simulation region;

[0007] determining a gravity anomaly three-dimensional propagation formula of the simulation region based on the gravity anomaly model;

[0008] performing three-dimensional finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula to obtain one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula;

[0009] performing subdivision interpolation on the one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula using quadratic finite element method and then accumulating to obtain the one-dimensional gravity spectrum integral after accumulation;

[0010] The one-dimensional gravity spectrum integral after accumulation is converted to the spatial domain by three finite element method Fourier transform to obtain the spatial domain gravity anomaly.

[0011] In some embodiments, a gravity anomaly model of a simulation region is constructed, specifically:

[0012] The simulation region, the anomaly region, and the residual density are determined.

[0013] The gravity anomaly model is constructed based on the simulation region, the anomaly region, and the residual density.

[0014] In some embodiments, a three-dimensional propagation formula of the gravity anomaly of the simulation region is determined based on the gravity anomaly model, specifically:

[0015] The three-dimensional propagation formula of the gravity anomaly of the simulation region is determined based on the simulation region, the anomaly region, and the anomaly body density.

[0016] The three-dimensional propagation formula of the gravity anomaly is specifically:

[0017]

[0018] Wherein, x, y, z respectively represent the observation point position of the simulation region, x', y', z' respectively represent the anomaly point position of the anomaly region; The distance from the observation point of the simulation region to the point of the anomaly region is represented; v represents the volume of the anomaly; u(r') represents the residual density; γ = 6.67 × 10 -11 N·m 2 ·kg -2 Is a constant.

[0019] In some embodiments, the three-dimensional propagation formula of the gravity anomaly is subjected to three finite element method Fourier transform to obtain a one-dimensional gravity spectrum integral of the three-dimensional propagation formula of the gravity anomaly, specifically:

[0020] The spatial domain of the simulation region is determined based on the three-dimensional propagation formula of the gravity anomaly.

[0021] The spatial domain is discretely profiled to obtain a plurality of element integrals.

[0022] The plurality of element integrals are subjected to three finite element method interpolation transform to obtain a plurality of transformed spectral element integrals.

[0023] The plurality of transformed spectral element integrals are subjected to finite element quadratic interpolation to obtain interpolated spectral elements.

[0024] All element accumulations are performed on the interpolated spectral elements to obtain spectral element accumulation results.

[0025] The one-dimensional gravity spectrum integral of the three-dimensional propagation formula of the gravity anomaly is obtained based on the spectrum unit accumulation result.

[0026] In some embodiments, the one-dimensional gravity spectrum integral of the three-dimensional propagation formula of the gravity anomaly is specifically:

[0027]

[0028] wherein k x , k y respectively represent the gravity spectrum wave number corresponding to the x and y directions of the simulation area, represent the sum of the gravity spectrum wave number; ～ u(k x , k y , z') represents the Fourier transform result of the residual density wave number of the gravity anomaly residual density u(r').

[0029] In some embodiments, the one-dimensional gravity spectrum integral of the three-dimensional propagation formula of the gravity anomaly is accumulated after being divided and interpolated by using the quadratic finite element method, to obtain the accumulated one-dimensional gravity spectrum integral, which is specifically:

[0030] The one-dimensional gravity spectrum integral of the three-dimensional propagation formula of the gravity anomaly is discretized into a plurality of one-dimensional integral unit integrals in the z direction;

[0031] In each one-dimensional integral unit integral, the residual density is represented as a quadratic interpolation function of the node residual density, to obtain the residual density converted unit integral;

[0032] The density converted unit integral is accumulated to obtain the accumulated one-dimensional spectrum integral.

[0033] In some embodiments, the spatial domain gravity anomaly is specifically:

[0034]

[0035] wherein w(x, y, z) is the spatial domain gravity anomaly.

[0036] Correspondingly, the application further provides a three-dimensional gravity anomaly forward modeling device, which comprises:

[0037] A model construction module is configured to construct a gravity anomaly model comprising a simulation area.

[0038] A determination module is configured to determine a three-dimensional propagation formula of the gravity anomaly of the simulation area based on the gravity anomaly model.

[0039] A one-dimensional gravity spectrum integral obtaining module is configured to perform cubic finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula to obtain one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula;

[0040] An accumulation module is configured to accumulate the one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula after subdivision interpolation by quadratic finite element method to obtain accumulated one-dimensional gravity spectrum integral.

[0041] A spatial domain gravity anomaly module is configured to convert the accumulated one-dimensional gravity spectrum integral to the spatial domain by cubic finite element method Fourier transform to obtain the spatial domain gravity anomaly.

[0042] One of the embodiments of the present application further provides a computing device, comprising a memory and a processor, wherein the memory is configured to store computer executable instructions, and the processor is configured to execute the computer executable instructions, and the computer executable instructions, when executed by the processor, implement the steps of the three-dimensional gravity anomaly forward modeling method according to any one of the above.

[0043] One of the embodiments of the present application further provides a computer readable storage medium, which stores computer executable instructions, and the computer executable instructions, when executed by the processor, implement the steps of the three-dimensional gravity anomaly forward modeling method according to any one of the above.

[0044] By applying the above technical solution, a three-dimensional gravity anomaly forward modeling method is provided, which comprises the following steps: constructing a gravity anomaly model comprising a simulation region; determining a gravity anomaly three-dimensional propagation formula of the simulation region based on the gravity anomaly model; performing cubic finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula to obtain one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula; accumulating the one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula after subdivision interpolation by quadratic finite element method to obtain accumulated one-dimensional gravity spectrum integral; and converting the accumulated one-dimensional gravity spectrum integral to the spatial domain by cubic finite element method Fourier transform to obtain the spatial domain gravity anomaly. The coefficients in the discrete unit are interpolated by cubic finite element method, and the interpolation results in the unit can be deduced to the spatial domain gravity anomaly, thereby improving the calculation precision of the three-dimensional gravity anomaly calculation, and the algorithm can be applied to the three-dimensional gravity anomaly calculation under complex conditions. BRIEF DESCRIPTION OF DRAWINGS

[0045] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0046] Figure 1 is a flow chart of a three-dimensional gravity anomaly forward simulation method provided by an embodiment of the present application;

[0047] Figure 2 is a Fourier transform numerical solution and analytical solution contour comparison and relative error graph provided by an embodiment of the present application;

[0048] Figure 3 is a Fourier transform numerical solution and analytical solution contour comparison and relative error graph provided by an embodiment of the present application;

[0049] Figure 4 is a structure schematic diagram of a gravity anomaly model provided by an embodiment of the present application;

[0050] Figure 5 is an effect schematic diagram of a numerical solution, an analytical solution and an absolute error of a Fourier transform gravity anomaly numerical simulation based on a cubic finite element method on a z=0m plane provided by an embodiment of the present application;

[0051] Figure 6 is a structure schematic diagram of a three-dimensional gravity anomaly forward simulation device provided by an embodiment of the present application;

[0052] Figure 7 is a structure block diagram of a computing device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0053] In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present specification. However, the present specification can be practiced without the specific details, other than in the examples, and it is understood that the scope of the present specification is not limited to the details below.

[0054] The terminology used in one or more embodiments of the present specification is for the purpose of describing particular embodiments only and is not intended to be limiting of one or more embodiments of the present specification. As used in one or more embodiments of the present specification and the accompanying claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and / or "comprising," when used in one or more embodiments of the present specification, specify the presence of stated features, integers, steps, operations, elements, and / or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0055] It should be understood that, although the terms first, second, etc. can be employed in describing various information in one or more embodiments, the information should not be limited to these terms. These terms are only used to differentiate one piece of information from another. For example, a first can be termed a second, and, similarly, a second can be termed a first, without departing from the scope of one or more embodiments. The word "if' can be interpreted to mean "when" or "upon" or "in response to determining" depending on the context.

[0056] As shown in Figure 1 The present application proposes a three-dimensional gravity anomaly forward simulation method, which comprises the following steps:

[0057] Step S101, constructing a gravity anomaly model comprising a simulation region.

[0058] In a possible implementation, the gravity anomaly model comprising the simulation region is constructed, specifically:

[0059] determining the simulation region, the anomaly region and the residual density;

[0060] constructing the gravity anomaly model based on the simulation region, the anomaly region and the residual density.

[0061] In this embodiment, the gravity anomaly model comprising the target region is constructed, any observation point (x, y, z) in the simulation region, any anomaly point (x', y', z') in the anomaly region and the residual density u(r') are given.

[0062] Step S102, determining a three-dimensional propagation formula of the gravity anomaly of the simulation region based on the gravity anomaly model.

[0063] In a possible implementation, the three-dimensional propagation formula of the gravity anomaly of the simulation region is determined based on the gravity anomaly model, specifically:

[0064] determining the three-dimensional propagation formula of the gravity anomaly of the simulation region based on the simulation region, the anomaly region and the anomaly body density;

[0065] wherein the three-dimensional propagation formula of the gravity anomaly is specifically:

[0066]

[0067] wherein x, y and z respectively represent the observation point position of the simulation region, x', y' and z' respectively represent the anomaly point position of the anomaly region; represents the distance from the observation point of the simulation region to the point of the anomaly region; v represents the volume of the anomaly; u(r') represents the residual density; γ = 6.67 x 10 -11N·m 2 ·kg -2 is a constant.

[0068] In this embodiment, according to the simulation area, a gravity anomaly three-dimensional propagation formula of the simulation area is determined:

[0069]

[0070] wherein x, y, z respectively represent the observation point position of the simulation area, x', y', z' respectively represent the anomaly point position of the anomaly area; represents the distance from the observation point of the simulation area to the point of the anomaly area; v represents the volume of the anomaly; u(r') represents the residual density; γ = 6.67 x 10 -11 N·m 2 ·kg -2 is a constant.

[0071] In step S103, a one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula is obtained by performing a cubic finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula.

[0072] In one possible implementation, the one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula is obtained by performing a cubic finite element method Fourier transform on the gravity anomaly three-dimensional propagation formula, and specifically:

[0073] The spatial domain of the simulation area is determined based on the gravity anomaly three-dimensional propagation formula;

[0074] The spatial domain is discretely profiled to obtain a plurality of element integrals;

[0075] The parameters of the plurality of element integrals are subjected to a cubic finite element method interpolation transform to obtain a plurality of spectrum element integrals after transformation;

[0076] The plurality of spectrum element integrals after transformation are subjected to finite element quadratic interpolation to obtain interpolated spectrum elements;

[0077] The interpolated spectrum elements are subjected to all element accumulation to obtain spectrum element accumulation results;

[0078] The one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula is obtained based on the spectrum element accumulation results.

[0079] In one possible implementation, the one-dimensional gravity spectrum integral of the gravity anomaly three-dimensional propagation formula is specifically:

[0080]

[0081] wherein k x , k yrespectively represent the corresponding spectral wave numbers in x, y directions of the simulation area, represent the sum of the spectral wave numbers; ～ u(k x , k y , z') represents the Fourier transform result of the residual density u(r') of the gravity anomaly - residual density wave number spectrum.

[0082] In this embodiment, the cubic finite element method Fourier transform is used for the gravity anomaly spatial domain integral expression of step one, and the three-dimensional spatial domain is converted to a one-dimensional wave number domain:

[0083]

[0084] In formula (2), k x , k y respectively represent the corresponding spectral wave numbers in x, y directions of the simulation area, represent the sum of the spectral wave numbers; ～ u(k x , k y , z') represents the Fourier transform result of the residual density u(r') of the gravity anomaly - residual density wave number spectrum.

[0085] The cubic finite element method Fourier transform used in the present application has high calculation accuracy, and the specific idea is as follows:

[0086] For the three-dimensional residual density u(r') in formula (1), the one-dimensional (x direction) Fourier transform calculation wave number spectrum expression in the gravity anomaly field simulation is as follows:

[0087] (3)

[0088] Wherein, i represents a complex number, x represents the position of the simulation area measuring point, k x represents the corresponding spectral wave number in the x direction of the simulation area, u(x) represents the residual density in the x direction, ～ u(k x ) represents the one-dimensional Fourier transform result of the residual density u(x) of the gravity anomaly - residual density wave number spectrum.

[0089] Discretely divide (3):

[0090]

[0091] Wherein, ～ u(k x ) is the residual density wave number spectrum unit integral, and M is the number of units.

[0092] The cubic finite element method interpolation is used for the integral unit coefficient in formula (4):

[0093] The parameters are interpolated by cubic finite element, which can be expressed as:

[0094]

[0095] Where, u1, u2, u3, u4 and x1, x2, x3, x4 represent the interpolation point residual density value and simulation area coordinates respectively, and N1, N2, N3, N4 represent the interpolation shape function of residual density, and the specific expression is as follows:

[0096]

[0097] Where, ξ represents the new space position of conversion. Bringing (6) into (5) gives

[0098]

[0099] According to formulas (5), (6) and (7) and the derivation rule, we can get:

[0100]

[0101] The residual density wave spectrum integral formula (4) is converted from x domain to ξ domain, and the one-dimensional residual density wave spectrum integral expression becomes:

[0102]

[0103] Bringing (8) into formula (9), the one-dimensional residual density wave spectrum calculation formula can be obtained:

[0104]

[0105] Where,

[0106]

[0107] c = (x1-27x2+27x3-x4)

[0108] When the wave number k x is equal to zero, the residual density wave spectrum formula (10) is simplified as:

[0109]

[0110] After solving each integral element in the above step three ～ u e (k x ), the final one-dimensional spectrum result is obtained by element accumulation for different discrete elements:

[0111]

[0112] The two-dimensional Fourier transform expression for calculating the residual density wave spectrum in the frequency domain simulation of the gravity anomaly field is:

[0113]

[0114] wherein k x , k y respectively represent the spectral wave numbers corresponding to the x and y directions of the simulation area, ～ u(k x , k y ) is the two-dimensional residual density wave spectrum.

[0115] The strategy adopted by the present application is to sequentially perform one-dimensional integral calculation on the x and y directions respectively, and then combine them to obtain the two-dimensional Fourier transform result. First, one-dimensional integral is performed on the x direction:

[0116]

[0117] The above formula (9) is the same as the one-dimensional Fourier transform formula (1), and can be solved by using the one-dimensional Fourier transform calculation method. On this basis, one-dimensional integral is performed on the y direction:

[0118]

[0119] The above formula (15) is the same as the one-dimensional Fourier transform formula (1), and the residual density wave spectrum result ～ u(k x , k y ) can be obtained by using the one-dimensional Fourier transform algorithm.

[0120] In step S104, the one-dimensional spectral integral of the gravity anomaly three-dimensional propagation formula is divided and interpolated by using the quadratic finite element method, and then accumulated to obtain the accumulated one-dimensional spectral integral.

[0121] In one possible implementation, the one-dimensional spectral integral of the gravity anomaly three-dimensional propagation formula is divided and interpolated by using the quadratic finite element method, and then accumulated to obtain the accumulated one-dimensional integral, specifically:

[0122] The one-dimensional spectral integral of the gravity anomaly three-dimensional propagation formula is discretized into multiple one-dimensional integral unit integrals according to the z direction;

[0123] In each one-dimensional spectral integral unit integral, the residual density is represented as a quadratic difference function of the node residual density to obtain the density-converted unit spectral integral;

[0124] The density-converted unit spectral integral is accumulated to obtain the accumulated one-dimensional spectral integral.

[0125] In this embodiment, the transformed one-dimensional spectral integral is divided and interpolated by using the quadratic finite element method, and then accumulated.

[0126] Discretizing (2) along z direction into N elements, we can get

[0127]

[0128] where I denotes the element spectral integral.

[0129] Within each element, the residual density is expressed as a quadratic interpolation function of the nodal residual densities:

[0130]

[0131] In (17), the are the wave-number domain residual density values at the three nodes of the element, N j , N p , N m are quadratic shape functions, and their explicit expressions are:

[0132]

[0133] Substituting (17) into the element integral in (16), we can get

[0134]

[0135] where z′ j , z′ m , z p are the three interpolation points of the element integral, and the specific sampling interval of the vertical element can be flexibly selected according to the residual density.

[0136] In step S105, the one-dimensional integral after accumulation is converted to the spatial domain by a cubic finite element method Fourier transform to obtain a spatial domain gravity anomaly.

[0137] In a possible implementation, the spatial domain gravity anomaly is specifically:

[0138]

[0139] where w(x, y, z) is the spatial domain gravity anomaly.

[0140] The following gives a specific example of the application to verify the correctness of the method proposed by the application and the effectiveness of the method for numerical simulation of gravity anomaly field.

[0141] 1. Correctness verification example

[0142] The Fourier forward transform is performed on a one-dimensional spatial domain Gaussian function to obtain a one-dimensional wave-number domain Gaussian function The Fourier inverse transform is performed, and the transform results are compared with the analytical solution to verify the correctness of the method, wherein a = 0.001. The positive transform space sampling range is (-100m, 100m), the space sampling points are 101, and the sampling interval is 2m; the inverse transform wave number k x The sampling range is (-1.5, 1.5), and the wave number sampling number is 121.

[0143] Figure 2 The Fourier transform numerical solution and the analytical solution are compared, and the relative error graph is shown in the following table: Figure 3 The Fourier transform numerical solution and the analytical solution are compared, and the relative error graph is shown in the following table: Figure 2 a and Figure 3 It can be seen that the numerical solution of the positive and inverse transform is basically coincident with the analytical solution, and the coincidence degree is high; from Figure 2 b and Figure 3 It can be seen that the relative error of the positive and inverse transform is below 0.06%, and the calculation precision is high. The positive and inverse transform calculation results show that the Fourier transform theoretical method proposed in the application is correct and has high calculation precision.

[0144] 2. Gravity anomaly numerical simulation example

[0145] The model is shown in the following table: Figure 4 The simulation area range is 10km*10km*5km, the prismatic anomaly body is located in the center of the simulation area, the anomaly body size is 4km*4km*2km, and the top surface is buried at a depth of 2km. The simulation area grid is divided into 101*101*101, and the horizontal sampling interval is 100m. The vertical sampling interval is 50m. The residual density p = 1000kg*m -3 The wave number sampling range is (-3*1.0-2, 3*1.0-2), and the observation plane z = 0m o

[0146] Figure 5 The numerical solution, the analytical solution and the absolute error of the Fourier transform gravity anomaly numerical simulation based on the cubic finite element method in the z = 0m plane are shown in the following table: Figure 6 It can be seen that the absolute error of the gravity anomaly field is about 0.01mGal, which is lower than the analytical solution of the gravity anomaly field by three orders of magnitude. The calculation result shows that the method has high calculation precision in the numerical simulation of the gravity anomaly field.

[0147] By applying the above technical scheme, a three-dimensional gravity anomaly forward simulation method is provided, which comprises: constructing a gravity anomaly model comprising a simulation region; determining a three-dimensional propagation formula of gravity anomaly of the simulation region based on the gravity anomaly model; performing three-order finite element method Fourier transform on the three-dimensional propagation formula of gravity anomaly to obtain one-dimensional integral of the three-dimensional propagation formula of gravity anomaly; performing subdivision and interpolation on the one-dimensional integral of the three-dimensional propagation formula of gravity anomaly by using quadratic finite element method and then accumulating to obtain one-dimensional integral after accumulation; converting the one-dimensional integral after accumulation to a spatial domain by three-order finite element method Fourier transform to obtain spatial domain gravity anomaly, and deriving the spatial domain gravity anomaly by interpolating the coefficients in the discrete unit by using three-order finite element method, which improves the calculation precision of three-dimensional gravity anomaly forward calculation, and the algorithm can be applied to three-dimensional gravity anomaly forward calculation under complex conditions.

[0148] The embodiment of the present application further provides a three-dimensional gravity anomaly forward simulation device, as shown in the accompanying drawings, which comprises: Figure 6

[0149] A preprocessing module 601 is configured to collect historical seismic exploration data, and preprocess single-channel seismic waveform data in the historical seismic exploration data to obtain preprocessed single-channel seismic waveform data and amplitude-related auxiliary information.

[0150] A division module 602 is configured to divide the preprocessed single-channel seismic waveform data into a training set and a test set.

[0151] A model construction module 603 is configured to construct a deep learning convolutional neural network, train the deep learning convolutional neural network by taking the training set and the auxiliary information as input, and obtain a trained neural network model.

[0152] A test module 604 is configured to test the trained neural network model by using the test set, and obtain a seismic data classification model.

[0153] A type determination module 605 is configured to preprocess to-be-classified seismic waveform data to obtain preprocessed to-be-classified seismic waveform data, and determine the type of the preprocessed to-be-classified seismic waveform data based on the seismic data classification model and a cosine similarity algorithm.

[0154] Figure 7 A structural block diagram of a computing device 400 provided by one embodiment of the present specification is shown. The components of the computing device 400 include but are not limited to a memory 410 and a processor 420. The processor 420 is connected with the memory 410 through a bus 430, and a database 450 is used to save data.

[0155] ​The computing device 400 also includes an access device 440 that enables the computing device 400 to communicate via one or more networks 460. Examples of these networks include a public switched telephone network (PSTN), a local area network (LAN), a wide area network (WAN), a personal area network (PAN), or combinations of these and / or other types of networks. The access device 440 can include one or more of any type of network interface (for example, a network interface card (NIC)), such as an IEEE 802.11 wireless local area network (WLAN) wireless interface, a Worldwide Interoperability for Microwave Access (Wi-MAX) interface, an Ethernet interface, a Universal Serial Bus (USB) interface, a cellular network interface, a Bluetooth interface, a Near Field Communication (NFC).

[0156] In one embodiment of the present specification, the above-described components of the computing device 400, as well as other components not shown in FIG. 4, can be connected to each other by, for example, a bus. It should be understood that the computing device structure diagram shown is for the purpose of example only and is not a limitation on the scope of the present specification. Other components can be added or replaced as needed by those skilled in the art. Figure 7 In one embodiment of the present specification, the above-described components of the computing device 400, as well as other components not shown in FIG. 4, can be connected to each other by, for example, a bus. It should be understood that the computing device structure diagram shown is for the purpose of example only and is not a limitation on the scope of the present specification. Other components can be added or replaced as needed by those skilled in the art. Figure 7 The computing device structure diagram shown is for the purpose of example only and is not a limitation on the scope of the present specification. Other components can be added or replaced as needed by those skilled in the art.

[0157] The computing device 400 can be any type of stationary or mobile computing device, including a mobile computer or mobile computing device (for example, a tablet computer, a personal digital assistant, a laptop computer, a notebook computer, a netbook, and the like), a mobile phone (for example, a smartphone), a wearable computing device (for example, a smartwatch, smartglasses, and the like), or other types of mobile devices, or a stationary computing device such as a desktop computer or a personal computer (PC). The computing device 400 can also be a mobile or stationary server.

[0158] The processor 420 is configured to execute computer-executable instructions to implement the steps of the method for simulating three-dimensional gravity anomaly forward modeling. The above is a schematic solution of the computing device of the embodiment. It should be noted that the technical solution of the computing device and the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling belong to the same concept, and details of the technical solution of the computing device not described in detail can be referred to the description of the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling.

[0159] The embodiment of the present specification also provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the steps of the method for simulating three-dimensional gravity anomaly forward modeling.

[0160] The above is a schematic solution of the computer-readable storage medium of the embodiment. It should be noted that the technical solution of the storage medium and the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling belong to the same concept, and details of the technical solution of the storage medium not described in detail can be referred to the description of the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling.

[0161] The embodiment of the present specification also provides a computer program, which, when executed in a computer, causes the computer to perform the steps of the method for simulating three-dimensional gravity anomaly forward modeling.

[0162] The above is a schematic solution of the computer program of the embodiment. It should be noted that the technical solution of the computer program and the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling belong to the same concept, and details of the technical solution of the computer program not described in detail can be referred to the description of the technical solution of the method for simulating three-dimensional gravity anomaly forward modeling.

[0163] The above describes specific embodiments of the present specification. Other embodiments are within the scope of the appended claims. In some cases, the acts or steps recited in the claims can be performed in a different order than the order in which they are recited and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order to achieve the desired results. In some implementations, multitasking and parallel processing can be advantageous.

[0164] The computer readable medium can include any entity or apparatus capable of carrying the computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, software distribution medium, etc. It should be noted that the computer readable medium can include appropriate additions or subtractions according to the requirements of legislation and patent practice in the jurisdiction, for example, in some jurisdictions, according to legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunication signals.

[0165] It should be noted that for the foregoing method embodiments, the descriptions are each simply a combination of a series of acts for the sake of brevity, but those skilled in the art should know that the present application is not limited by the order of the acts described, because some steps can be performed in other orders or at the same time in accordance with the present application. In addition, those skilled in the art should know that the embodiments described in the specification are all preferred embodiments, and the acts and modules involved are not necessarily essential to the present application.

[0166] In the above embodiments, the description of each embodiment has its own focus, and the parts not described in detail in a certain embodiment can be referred to the relevant description of other embodiments.

[0167] The preferred embodiments of the present application disclosed above are only used to help explain the present application. The alternative embodiments do not describe all the details and limit the present application to the specific embodiments described. Obviously, according to the content of the present application, many modifications and changes can be made. The present application selects and specifically describes these embodiments in order to better explain the principles and practical applications of the present application, so that those skilled in the art can well understand and use the present application. The present application is limited only by the claims and their full scope and equivalents.

Claims

1. A three-dimensional forward modeling method for gravity anomalies, characterized in that, The method includes: Construct a gravity anomaly model including the simulated region; Based on the gravity anomaly model, the three-dimensional propagation formula of gravity anomaly in the simulated region is determined; A three-dimensional finite element Fourier transform is performed on the three-dimensional propagation formula of gravity anomaly to obtain the one-dimensional gravity spectrum integral of the three-dimensional propagation formula of gravity anomaly. The one-dimensional integral of the three-dimensional propagation formula of gravity anomaly is obtained by interpolating and subdividing it using the quadratic finite element method and then summing the results. The accumulated one-dimensional gravity spectrum integral is converted to the spatial domain by using the three-dimensional finite element method Fourier transform, thus obtaining the spatial domain gravity anomaly.

2. The method according to claim 1, characterized in that, Construct a gravity anomaly model including the simulated region, specifically as follows: Determine the simulated region, the abnormal region, and the remaining density; The gravity anomaly model is constructed based on the simulated region, the anomaly region, and the remaining density.

3. The method according to claim 1, characterized in that, Based on the aforementioned gravity anomaly model, the three-dimensional propagation formula for gravity anomalies in the simulated region is determined as follows: The three-dimensional propagation formula of gravity anomaly in the simulated region is determined based on the simulated region, the abnormal region, and the density of the abnormal body. The specific formula for the three-dimensional propagation of gravity anomalies is as follows: Where x, y, and z represent the locations of observation points in the simulated region, and x′, y′, and z′ represent the locations of anomaly points in the anomaly region, respectively; γ represents the distance from the observation point in the simulated area to the point in the anomalous area; v represents the volume of the anomalous area; u(r′) represents the residual density; γ = 6.67 × 10⁻⁶. -11 N·m 2 ·kg -2 It is a constant.

4. The method according to claim 1, characterized in that, Performing a cubic finite element Fourier transform on the three-dimensional propagation formula of the gravity anomaly yields the one-dimensional gravity spectrum integral of the formula, specifically: The spatial domain of the simulated region is determined based on the three-dimensional propagation formula of gravity anomaly. The spatial domain is discretized to obtain multiple unit integrals; The parameters of the multiple element integrals are subjected to a cubic finite element method interpolation transformation to obtain the transformed multiple spectral element integrals; The transformed spectral elements are integrated and then subjected to finite element quadratic interpolation to obtain the interpolated spectral elements. The interpolated spectral elements are summed to obtain the spectral element summation result. The one-dimensional gravity spectrum integral of the three-dimensional propagation formula of gravity anomaly is obtained based on the accumulated result of the spectral units.

5. The method according to claim 4, characterized in that, The one-dimensional gravity spectrum integral of the three-dimensional propagation formula for gravity anomalies is specifically as follows: Where, k x ,k y These represent the gravitational wavenumbers in the x and y directions of the simulated region, respectively. This represents the sum of the wavenumbers of the gravity spectrum; The Fourier transform result of the residual density u(r′) of the gravity anomaly is represented as the residual density wavenumber spectrum. This represents the gravity wavenumber spectrum.

6. The method according to claim 5, characterized in that, The one-dimensional gravity spectrum integral of the three-dimensional propagation formula of gravity anomaly is obtained by subdividing and interpolating using the quadratic finite element method and then summing the results. Specifically: The one-dimensional gravity spectrum integral of the three-dimensional propagation formula of gravity anomaly is discretized into multiple unit integrals of one-dimensional integrals along the z-direction; Within each element integral of the one-dimensional gravity spectrum integral, the residual density is expressed as a quadratic difference function of the nodal residual density, thus obtaining the element integral after residual density transformation. The unit integrals after the residual density conversion are accumulated to obtain the accumulated one-dimensional gravity spectrum.

7. The method according to claim 6, characterized in that, The aforementioned spatial gravity anomaly is specifically as follows: Where w(x,y,z) represents the gravity anomaly in the spatial domain.

8. A three-dimensional gravity anomaly forward modeling simulation device, characterized in that, The device includes: The model building module is used to construct a gravity anomaly model that includes the simulated region. The determination module is used to determine the three-dimensional propagation formula of gravity anomaly in the simulated region based on the gravity anomaly model; A one-dimensional gravity spectrum integral acquisition module is used to perform a cubic finite element Fourier transform on the three-dimensional propagation formula of gravity anomaly to obtain the one-dimensional gravity spectrum integral of the three-dimensional propagation formula of gravity anomaly. The accumulation module is used to accumulate the one-dimensional integral of the three-dimensional propagation formula of gravity anomaly after performing subdivision and interpolation using the quadratic finite element method, so as to obtain the accumulated one-dimensional gravity spectrum integral. The spatial domain gravity anomaly module is used to convert the accumulated one-dimensional gravity spectrum integral to the spatial domain through the three-dimensional finite element method Fourier transform, thereby obtaining the spatial domain gravity anomaly.

9. A computing device, characterized in that, include: Memory and processor; The memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions, which, when executed by the processor, implement the steps of the three-dimensional gravity anomaly forward modeling method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing computer-executable instructions that, when executed by a processor, implement the steps of the three-dimensional gravity anomaly forward modeling method according to any one of claims 1 to 7.

Images