In this paper, a Pareto inversion based global optimization approach, to obtain results of joint inversion of two types of geophysical data sets, is formulated. 2D magnetotelluric and gravity data were used for tests, but presented solution is flexible enough to be used for combination of any kind of two or more target functions, as long as misfits can be calculated and forward problems solved. To minimize dimensionality of the solution, space and introduce straightforward regularization Sharp Boundary Interface (SBI) method was applied. As a main optimization engine, Particle Swarm Optimization (PSO) was used. Synthetic examples based on a real geological model were used to test proposed approach and show its usefulness in practical applications.
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.
Joint inversion is a widely used geophysical method that allows model parameters to be obtained from the observed data. Pareto inversion results are a set of solutions that include the Pareto front, which consists of non-dominated solutions. All solutions from the Pareto front are considered the most feasible models from which a particular one can be chosen as the final solution. In this paper, it is shown that models represented by points on the Pareto front do not reflect the shape of the real model. In this contribution, a collective approach is proposed to interpret the geometry of models retrieved in inversion. Instead of choosing single solutions from the Pareto front, all obtained solutions were combined in one “heat map”, which is a plot representing the frequency of points belonging to all returned objects from the solution set. The conducted experiment showed that this approach limits the problem of equivalence and is a promising way of representing the geometry of the model that was retrieved in the inversion process.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.