Towards large‐scale stochastic refraction tomography: a comparison of three evolutionary algorithms

ABSTRACT The main goal of this study is to assess the potential of evolutionary algorithms to solve highly non‐linear and multi‐modal tomography problems (such as first arrival traveltime tomography) and their abilities to estimate reliable uncertainties. Classical tomography methods apply derivative‐based optimization algorithms that require the user to determine the value of several parameters (such as regularization level and initial model) prior to the inversion as they strongly affect the final inverted model. In addition, derivative‐based methods only perform a local search dependent on the chosen starting model. Global optimization methods based on Markov chain Monte Carlo that thoroughly sample the model parameter space are theoretically insensitive to the initial model but turn out to be computationally expensive. Evolutionary algorithms are population‐based global optimization methods and are thus intrinsically parallel, allowing these algorithms to fully handle available computer resources. We apply three evolutionary algorithms to solve a refraction traveltime tomography problem, namely the differential evolution, the competitive particle swarm optimization and the covariance matrix adaptation–evolution strategy. We apply these methodologies on a smoothed version of the Marmousi velocity model and compare their performances in terms of optimization and estimates of uncertainty. By performing scalability and statistical analysis over the results obtained with several runs, we assess the benefits and shortcomings of each algorithm.

[1]  Anne Auger,et al.  Comparison-based natural gradient optimization in high dimension , 2014, GECCO.

[2]  Petros Koumoutsakos,et al.  Reducing the Time Complexity of the Derandomized Evolution Strategy with Covariance Matrix Adaptation (CMA-ES) , 2003, Evolutionary Computation.

[3]  M. Nafi Toksöz,et al.  Nonlinear refraction traveltime tomography , 1998 .

[4]  Henri Calandra,et al.  High-performance 3D first-arrival traveltime tomography , 2010 .

[5]  Michael Andrew Christie,et al.  Application of Particle Swarms for History Matching in the Brugge Reservoir , 2010 .

[6]  S. Billings.  Simulated annealing for earthquake location , 1994 .

[7]  Jean Virieux,et al.  Seismic imaging of complex structures by non-linear traveltime inversion of dense wide-angle data: application to a thrust belt , 2002 .

[8]  Mark Noble,et al.  Stochastic seismic tomography by interacting Markov chains , 2016 .

[9]  Rini Akmeliawati,et al.  Performance Comparison of Differential Evolution and Particle Swarm Optimization in Constrained Optimization , 2012 .

[10]  Antonio Bolufé Röhler,et al.  Measuring the curse of dimensionality and its effects on particle swarm optimization and differential evolution , 2014, Applied Intelligence.

[11]  Alexandrine Gesret,et al.  New parameterizations for Bayesian seismic tomography , 2018 .

[12]  Andries Petrus Engelbrecht,et al.  Particle swarm optimization: Velocity initialization , 2012, 2012 IEEE Congress on Evolutionary Computation.

[13]  M. Sambridge,et al.  Seismic tomography with the reversible jump algorithm , 2009 .

[14]  Colin A. Zelt,et al.  Three‐dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin , 1998 .

[15]  Andries Petrus Engelbrecht,et al.  A study of particle swarm optimization particle trajectories , 2006, Inf. Sci..

[16]  C. Haberland,et al.  Bayesian inversion of refraction seismic traveltime data , 2018 .

[17]  Peter J. Angeline,et al.  Evolutionary Optimization Versus Particle Swarm Optimization: Philosophy and Performance Differences , 1998, Evolutionary Programming.

[18]  Esperanza García Gonzalo,et al.  Reservoir characterization and inversion uncertainty via a family of particle swarm optimizers , 2012 .

[19]  Rainer Storn.  Real-world applications in the communications industry - when do we resort to Differential Evolution? , 2017, 2017 IEEE Congress on Evolutionary Computation (CEC).

[20]  F. Boschetti,et al.  Inversion of seismic refraction data using genetic algorithms , 1996 .

[21]  Thomas Bodin,et al.  Probabilistic surface reconstruction from multiple data sets: An example for the Australian Moho , 2012 .

[22]  D. J. White,et al.  Two-Dimensional Seismic Refraction Tomography , 1989 .

[23]  J. Tronicke,et al.  Assessing uncertainty in refraction seismic traveltime inversion using a global inversion strategy , 2015 .

[24]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

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

[26]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[27]  M. Noble,et al.  Accurate 3-D finite difference computation of traveltimes in strongly heterogeneous media , 2014 .

[28]  Alexandrine Gesret,et al.  A parallel competitive Particle Swarm Optimization for non-linear first arrival traveltime tomography and uncertainty quantification , 2018, Comput. Geosci..

[29]  Roelof Versteeg,et al.  The Marmousi experience; velocity model determination on a synthetic complex data set , 1994 .

[30]  R. Eberhart,et al.  Comparing inertia weights and constriction factors in particle swarm optimization , 2000, Proceedings of the 2000 Congress on Evolutionary Computation. CEC00 (Cat. No.00TH8512).

[31]  Malcolm Sambridge,et al.  A Parallel Tempering algorithm for probabilistic sampling and multimodal optimization , 2014 .

[32]  W. Rabbel,et al.  On the application of Particle Swarm Optimization strategies on Scholte-wave inversion , 2012 .

[33]  A. Montanari,et al.  Syn‐convergent extension observed using the RETREAT GPS network, northern Apennines, Italy , 2012 .

[34]  C. D. Boor,et al.  On Calculating B-splines , 1972 .

[35]  A. Tarantola,et al.  Inverse problems = Quest for information , 1982 .

[36]  Daniele Boiero,et al.  Improved Monte Carlo inversion of surface wave data , 2008 .

[37]  Xianhai Song,et al.  Application of particle swarm optimization to interpret Rayleigh wave dispersion curves , 2012 .

[38]  G. Amdhal,et al.  Validity of the single processor approach to achieving large scale computing capabilities , 1967, AFIPS '67 (Spring).

[39]  Michael Andrew Christie,et al.  Comparison of Stochastic Sampling Algorithms for Uncertainty Quantification , 2010 .

[40]  B. Růžek,et al.  Differential Evolution Algorithm in the Earthquake Hypocenter Location , 2001 .

[41]  Jens Tronicke,et al.  Crosshole traveltime tomography using particle swarm optimization: A near-surface field example , 2012 .

[42]  Michael Andrew Christie,et al.  Population MCMC methods for history matching and uncertainty quantification , 2010, Computational Geosciences.

[43]  Malcolm Sambridge,et al.  Earthquake hypocenter location using genetic algorithms , 1993, Bulletin of the Seismological Society of America.

[44]  Henri Calandra,et al.  First-arrival traveltime tomography based on the adjoint-state method , 2009 .

[45]  Mrinal K. Sen,et al.  Bayesian inference, Gibbs' sampler and uncertainty estimation in geophysical inversion , 1996 .