Numerical simulation method, device and computer equipment for gravity field based on complex terrain

Abstract

The present application relates to a method, device, computer equipment and storage medium for numerical simulation of gravity field based on complex terrain. The method includes: constructing a three-dimensional cuboid model containing complex terrain and exploration targets, performing grid division to obtain a plurality of small cuboids, and setting the spatial domain density value of the small cuboids; Gaussian parameters are set to obtain the offset wavenumber. According to the terrain information, the three-dimensional cuboid model is divided into the lower part model and the upper part model, and the frequency domain gravity anomaly field of multiple observation surfaces is superimposed to calculate the two-dimensional Gaussian Fourier transform. The inverse transformation is used to obtain the gravity anomaly field in the space domain of multiple observation surfaces, and the gravity anomaly field in the space domain of the complex terrain undulation surface is obtained through the spline interpolation algorithm. The invention retains the vertical direction in the space domain, has flexible grid division, and is suitable for simulating complex terrain and complex geological bodies. Multiple offset wave numbers are convenient for parallel calculation, and the calculation efficiency of gravity anomaly is improved.

Description

technical field [0001] The present application relates to the technical field of computer numerical simulation, in particular to a method, device, computer equipment and storage medium for numerical simulation of gravity field based on complex terrain. Background technique [0002] At this stage, due to the increase of gravity observation data, the memory and calculation required for gravity field forward modeling have increased sharply, which is difficult for ordinary computers to realize, and has become a difficult problem for rapid three-dimensional density imaging. [0003] At present, the space domain method is still dominant due to its accuracy, but with the application of the standard FFT edge expansion method and NUFFT in the forward numerical simulation of gravity anomaly, the frequency domain method has become the processing method due to its simplicity, accuracy and high efficiency. The preferred method for large-scale forward modeling. Literature (Chai Y , Hinze...

AI Technical Summary

Problems solved by technology

[0004] In the existing technology, large-scale and large-scale airborne gravity surveys have emerged, but the terrain in actual aerial surveys is not a plane, but complex undulating terrain. Traditional space domain methods and frequency domain methods are difficult to meet the demands of massive high resolution data. Rapid inversion requires poor adaptability, so it is of practical significance to study a high-efficiency and high-precision numerical simulation method for 3D gravity anomalies based on complex terrain

Images