Gravitational field numerical simulation method and device based on complex terrain and computer equipment

Abstract

The invention relates to a gravitational field numerical simulation method and device based on complex terrains, computer equipment and a storage medium. The method comprises the following steps: constructing a three-dimensional cuboid model containing a complex terrain and an exploration target, performing mesh generation to obtain a plurality of small cuboids, and setting spatial domain density values of the small cuboids; obtaining an offset wave number according to grid parameters of grid subdivision and preset Gaussian parameters, dividing the three-dimensional cuboid model into a lower part model and an upper part model according to topographic information, performing block superposition to calculate a frequency domain gravity anomaly field of a plurality of observation surfaces, and performing two-dimensional Gaussian Fourier inversion to obtain spatial domain gravity anomaly fields of multiple observation surfaces, and obtaining a spatial domain gravity anomaly field of the complex topographic relief surface through a spline interpolation algorithm. The vertical direction is reserved in a spatial domain, mesh generation is flexible, the method and device are suitable for simulating complex terrains and complex geologic bodies, multiple offset wave numbers are convenient for parallel calculation, and the gravity anomaly calculation efficiency 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