Affiliation:
1. School of GeoSciences, University of Edinburgh, Edinburgh, United Kingdom
Abstract
SUMMARY
We test a fully non-linear method to solve Bayesian seismic tomographic problems using data consisting of observed traveltimes of first-arriving waves. Rather than using Monte Carlo methods to sample the posterior probability distribution that embodies the solution of the tomographic inverse problem, we use variational inference. Variational methods solve the Bayesian inference problem under an optimization framework by seeking the best approximation to the posterior distribution from a family of distributions, while still providing fully probabilistic results. We introduce a new variational method for geophysics—normalizing flows. The method models the posterior distribution by using a series of invertible and differentiable transforms—the flows. By optimizing the parameters of these transforms the flows are designed to convert a simple and analytically known probability distribution into a good approximation of the posterior distribution. Numerical examples show that normalizing flows can provide an accurate tomographic result including full uncertainty information while significantly decreasing the computational cost compared to Monte Carlo and other variational methods. In addition, this method provides analytic solutions for the posterior distribution rather than an ensemble of posterior samples. This opens the possibility that subsequent calculations that use the posterior distribution might be performed analytically.
Publisher
Oxford University Press (OUP)
Subject
Geochemistry and Petrology,Geophysics
Cited by
29 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献