Bayesian Inversion, Uncertainty Analysis and Interrogation Using Boosting Variational Inference
-
Published:2024-01
Issue:1
Volume:129
Page:
-
ISSN:2169-9313
-
Container-title:Journal of Geophysical Research: Solid Earth
-
language:en
-
Short-container-title:JGR Solid Earth
Author:
Zhao Xuebin1ORCID,
Curtis Andrew1ORCID
Affiliation:
1. School of Geosciences University of Edinburgh Edinburgh UK
Abstract
AbstractGeoscientists use observed data to estimate properties of the Earth's interior. This often requires non‐linear inverse problems to be solved and uncertainties to be estimated. Bayesian inference solves inverse problems under a probabilistic framework, in which uncertainty is represented by a so‐called posterior probability distribution. Recently, variational inference has emerged as an efficient method to estimate Bayesian solutions. By seeking the closest approximation to the posterior distribution within any chosen family of distributions, variational inference yields a fully probabilistic solution. It is important to define expressive variational families so that the posterior distribution can be represented accurately. We introduce boosting variational inference (BVI) as a computationally efficient means to construct a flexible approximating family comprising all possible finite mixtures of simpler component distributions. We use Gaussian mixture components due to their fully parametric nature and the ease with which they can be optimized. We apply BVI to seismic travel time tomography and full waveform inversion, comparing its performance with other methods of solution. The results demonstrate that BVI achieves reasonable efficiency and accuracy while enabling the construction of a fully analytic expression for the posterior distribution. Samples that represent major components of uncertainty in the solution can be obtained analytically from each mixture component. We demonstrate that these samples can be used to solve an interrogation problem: to assess the size of a subsurface target structure. To the best of our knowledge, this is the first method in geophysics that provides both analytic and reasonably accurate probabilistic solutions to fully non‐linear, high‐dimensional Bayesian full waveform inversion problems.
Publisher
American Geophysical Union (AGU)
Subject
Space and Planetary Science,Earth and Planetary Sciences (miscellaneous),Geochemistry and Petrology,Geophysics
Reference129 articles.
1. Abadi M. Agarwal A. Barham P. Brevdo E. Chen Z. Citro C. et al. (2016).Tensorflow: Large‐scale machine learning on heterogeneous distributed systems. arXiv preprint arXiv:1603.04467.
2. A simple statistical estimation procedure for Monte Carlo Inversion in geophysics
3. A tutorial on adaptive MCMC
4. Interrogation theory
5. On adaptive Markov chain Monte Carlo algorithms