Pareto Joint Inversion of 2D magnetometric and gravity data- synthetic study

Abstract

Pareto joint inversion for two or more data sets is an attractive and promising tool which eliminates target functions weighing and scaling, providing a set of acceptable solutions composing a Pareto front. In former author’s study MARIA (Modular Approach Robust Inversion Algorithm) was created as a flexible software based on global optimization engine (PSO) to obtain model parameters in process of Pareto joint inversion of two geophysical data sets. 2D magnetotelluric and gravity data were used for preliminary tests, but the software is ready to handle data from more than two geophysical methods. In this contribution, the authors’ magnetometric forward solver was implemented and integrated with MARIA. The gravity and magnetometry forward solver was verified on synthetic models. The tests were performed for different models of a dyke and showed, that even when the starting model is a homogeneous area without anomaly, it is possible to recover the shape of a small detail of the real model. Results showed that the group analysis of models on the Pareto front gives more information than the single best model. The final stage of interpretation is the raster map of Pareto front solutions analysis.

Keywords

• pareto joint inversion
• pso
• global optimization
• scientific computing
• magnetometry
• gravimetry