Abstract

We present a 3D high-resolution modeling methodology based on the interpretation of gravity gradient data and its joint inversion with the simulated annealing (SA) global optimization method. The geometry of the model, used as computational domain in the solution of the forward and inverse problems, is defined with an irregular ensemble of cubic prisms that recreates the interpreted shape of the target, derived from the results of applying different interpretation methods to the gravity gradient data. In our inversion approach, the linear inverse problem resulting from the domain discretization is not solved. Instead, the cost function is explored with the SA algorithm, its low misfit region is identified, and models belonging to it are selected for obtaining the mean model, which represents the most likely model among them, as well as for estimating its uncertainty. The SA inversion algorithm we applied was numerically optimized to reduce the computational burden required by the forward problem, and it was driven by optimal tuning parameters, determined by a parametric analysis. Tests on synthetic data show the efficiency of our methodology to obtain a model that approximates the synthetic target and the usefulness of the estimated uncertainty to complement the interpretation. Finally, by applying our methodology to gravity gradient data acquired over the Vinton dome located in Louisiana, USA, we obtained results that are in agreement with geological information and previous studies.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call