A Gravity Forward Modeling Acceleration Method

Abstract

The invention discloses a novel gravity forward acceleration method. The novel gravity forward acceleration method comprises steps that step1, a three-dimensional geological body is divided into a plurality of tetrahedron units; corresponding octree structures are built for observation points and the three-dimensional geological body; step2, the three-dimensional geological body is divided into a near region and a far region; step3, a near region gravity parameter and a far region gravity parameter are respectively calculated; the integral of the near region is directly calculated based on the gravity integral singularity-free analytical expressions of the tetrahedron units included by the near region, and then the near region gravity parameter is acquired; based on the octrees of the observation points and a source octree, the integral of the far region is calculated by adopting a self-adaptive rapid multipole expansion algorithm, and then the far region gravity parameter is acquired; step 4, the sum of the acquired near region gravity parameter and the acquired far region gravity parameter from the step3 is calculated, and an observation point gravity parameter is acquired. Based on unstructured mesh generation, the gravity forward calculation of the observations distributed anywhere is carried out, and the efficiency of the gravity forward calculation is greatly improved.

Description

technical field [0001] The invention relates to a gravity forward modeling method in the field of geophysics, in particular to a high-precision and fast forward modeling method for large-scale aviation gravity and satellite gravity. Background technique [0002] Gravity prospecting is an important method of geophysical prospecting, and is widely used in mineral resource prospecting, crustal structure and Moho interface detection, and environmental geophysical prospecting. The measured gravity data has a direct relationship with the subsurface density distribution body. Therefore, the inversion of subsurface density can be used to discriminate subsurface density anomalies. With the development of gravity exploration technology and the gradual maturity of airborne gravity technology, for a target with a given density distribution, gravity exploration often contains more observation points, and the time-consuming characteristics of traditional gravity forward modeling methods b...

AI Technical Summary

Problems solved by technology

The above two methods have fast computational efficiency and high precision, but both require uniform hexahedral meshing for 3D geological bodies, which makes the meshing difficult to fit complex models well, and the method is not suitable for dealing with 3D geology with terrain At the same time, this method can only calculate the gravity forward modeling of the observation point on the same plane

Images