3D gravity interface inversion constrained by a few points and its GPU acceleration

•Present a fast procedure for the interface inversion of density contrast.•Apply GPU parallel computing with CUDA using NVIDIA video cards.•Analyze the optimization strategies of GPU for the interface inversion.•Introduce controlled-point constraints the in interface inversion.

sWe present a fast procedure for the interface inversion of density contrast on the platform of GPU. Firstly, based on the iterative method in spatial domain, we introduce soft constraints by some few control points, making the inversion results closer to reality. In this problem, we assume that the density contrast decays with depth according to a parabolic law. Secondly, in order to meet the interface inversion for a wide range large-scale gravity data, we develop GPU program of the 3D interface inversion of basement relief, and give the corresponding optimization strategies and the speedup tests for the calculation results. We demonstrate the applicability and efficacy of this technique by inverting gravity anomalies caused by a synthetic model of a density interface and a set of real field data.