Seismic tomography using variational inference methods

Abstract

Seismic tomography is a methodology to image the interior of solid or fluid media, and is often used to map properties in the subsurface of the Earth. In order to better interpret the resulting images it is important to assess imaging uncertainties. Since tomography is significantly nonlinear, Monte Carlo sampling methods are often used for this purpose, but they are generally computationally intractable for large datasets and high-dimensional parameter spaces. To extend uncertainty analysis to larger systems we use variational inference methods to conduct seismic tomography. In contrast to Monte Carlo sampling, variational methods solve the Bayesian inference problem as an optimization problem, yet still provide probabilistic results. In this study, we applied two variational methods, automatic differential variational inference (ADVI) and Stein variational gradient descent (SVGD), to 2D seismic tomography problems using both synthetic and real data and we compare the results to those from two different Monte Carlo sampling methods. The results show that variational inference methods can produce accurate approximations to the results of Monte Carlo sampling methods at significantly lower computational cost, provided that gradients of parameters with respect to data can be calculated efficiently. We expect that the methods can be applied fruitfully to many other types of geophysical inverse problems.

26 pages, 14 figures

References in corpus (6)

• Variational Inference: A Review for Statisticians
• Markov Chain Monte Carlo and Variational Inference: Bridging the Gap
• Automatic Differentiation Variational Inference
• An introduction to sampling via measure transport
• On the precision of noise correlation interferometry
• Structured Stochastic Variational Inference

Cited by in corpus (10)

• 3D Bayesian Variational Full Waveform Inversion
• Learned multiphysics inversion with differentiable programming and machine learning
• Solving multiphysics-based inverse problems with learned surrogates and constraints
• Variational Prior Replacement in Bayesian Inference and Inversion
• Introducing Conceptual Geological Information into Bayesian Tomographic Imaging
• Annealed Stein Variational Gradient Descent for Improved Uncertainty Estimation in Full-Waveform Inversion
• Surface wave dispersion inversion using an energy likelihood function
• VIP -- Variational Inversion Package with example implementations of Bayesian tomographic imaging
• Bayesian Time-Lapse Full Waveform Inversion using Hamiltonian Monte Carlo
• Interrogating probabilistic inversion results for subsurface structural information