On the optimization of electromagnetic geophysical data: Application of the PSO algorithm

Abstract Particle Swarm optimization (PSO) algorithm resolves constrained multi-parameter problems and is suitable for simultaneous optimization of linear and nonlinear problems, with the assumption that forward modeling is based on good understanding of ill-posed problem for geophysical inversion. We apply PSO for solving the geophysical inverse problem to infer an Earth model, i.e. the electrical resistivity at depth, consistent with the observed geophysical data. The method doesn't require an initial model and can be easily constrained, according to external information for each single sounding. The optimization process to estimate the model parameters from the electromagnetic soundings focuses on the discussion of the objective function to be minimized. We discuss the possibility to introduce in the objective function vertical and lateral constraints, with an Occam-like regularization. A sensitivity analysis allowed us to check the performance of the algorithm. The reliability of the approach is tested on synthetic, real Audio-Magnetotelluric (AMT) and Long Period MT data. The method appears able to solve complex problems and allows us to estimate the a posteriori distribution of the model parameters.

[1]  Robert L. Parker,et al.  Optimal one-dimensional inversion and bounding of magnetotelluric apparent resistivity and phase measurements , 1996 .

[2]  S. Yuan,et al.  Swarm intelligence optimization and its application in geophysical data inversion , 2009 .

[3]  D. Oldenburg,et al.  A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems , 2004 .

[4]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[5]  James G. Berryman,et al.  Optimization and geophysical inverse problems , 2000 .

[6]  Junying Chen,et al.  Particle Swarm Optimization with Local Search , 2005, 2005 International Conference on Neural Networks and Brain.

[7]  L. Cagniard.  Basic theory of the magneto-telluric method of geophysical prospecting , 1953 .

[8]  Esben Auken,et al.  Layered and laterally constrained 2D inversion of resistivity data , 2004 .

[9]  Alan G. Jones,et al.  The magnetotelluric method : theory and practice , 2012 .

[10]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[11]  J. Wait.  On the relation between telluric currents and the Earth's magnetic field , 1954 .

[12]  Yue Shi,et al.  A modified particle swarm optimizer , 1998, 1998 IEEE International Conference on Evolutionary Computation Proceedings. IEEE World Congress on Computational Intelligence (Cat. No.98TH8360).

[13]  M. N. Berdichevsky,et al.  Basic principles of interpretation of magnetotelluric sounding curves , 1976 .

[14]  Esperanza García Gonzalo,et al.  PSO: A powerful algorithm to solve geophysical inverse problems: Application to a 1D-DC resistivity case , 2010 .

[15]  J. T. Smith,et al.  Magnetotelluric inversion for minimum structure , 1988 .

[16]  Alberto Godio,et al.  Multi Population Genetic Algorithm to estimate snow properties from GPR data , 2016 .

[17]  James Kennedy,et al.  Particle swarm optimization , 2002, Proceedings of ICNN'95 - International Conference on Neural Networks.

[18]  Andries P. Engelbrecht,et al.  Computational Intelligence: An Introduction , 2002 .

[19]  Alberto Godio,et al.  Improved Monte Carlo 1D-Inversion of vertical electrical sounding and time-domain electromagnetic data , 2010 .

[20]  Mrinal K. Sen,et al.  Global Optimization Methods in Geophysical Inversion , 1995 .

[21]  R. Shaw,et al.  Particle swarm optimization : A new tool to invert geophysical data , 2007 .

[22]  Alan G. Jones,et al.  A multi-station magnetotelluric study in southern Scotland – I. Fieldwork, data analysis and results , 1979 .

[23]  J. Fernández-Martínez,et al.  Particle swarm optimization applied to solving and appraising the streaming-potential inverse problem , 2010 .

[24]  Adele Manzella,et al.  Electromagnetic and DC methods for geothermal exploration in Italy – case studies and future developments , 2015 .

[25]  Michele Pipan,et al.  Rayleigh wave dispersion curve inversion via genetic algorithms and Marginal Posterior Probability Density estimation , 2007 .

[26]  Esperanza García Gonzalo,et al.  The PSO family: deduction, stochastic analysis and comparison , 2009, Swarm Intelligence.

[27]  A. K. Agarwal,et al.  Characterization of the magnetotelluric tensor in terms of its invariants , 2000 .