Bayesian seismic inversion: Measuring Langevin MCMC sample quality with kernels
KAUST DepartmentComputer, Electrical and Mathematical Science and Engineering (CEMSE) Division
Earth Science and Engineering Program
Extreme Computing Research Center
Physical Science and Engineering (PSE) Division
Permanent link to this recordhttp://hdl.handle.net/10754/667572
MetadataShow full item record
AbstractThe Bayesian framework is commonly used to quantify uncertainty in seismic inversion. To perform Bayesian inference, Markov chain Monte Carlo (MCMC) algorithms are regarded as the gold standard technique for sampling from the posterior probability distribution. Consistent MCMC methods have trouble for complex, high-dimensional models, and most methods scale poorly to large datasets, such as those arising in seismic inversion. As an alternative, approximate MCMC methods based on unadjusted Langevin dynamics offer scalability and more rapid sampling at the cost of biased inference. However, when assessing the quality of approximate MCMC samples for characterizing the posterior distribution, most diagnostics fail to account for these biases. In this work, we introduce the kernel Stein discrepancy (KSD) as a diagnostic tool to determine the convergence of MCMC samples for Bayesian seismic inversion. We demonstrate the use of the KSD for measuring sample quality and selecting the optimal Langevin MCMC algorithm for two Gaussian Bayesian inference problems.
CitationIzzatullah, M., Baptista, R., Mackey, L., Marzouk, Y., & Peter, D. (2020). Bayesian seismic inversion: Measuring Langevin MCMC sample quality with kernels. SEG Technical Program Expanded Abstracts 2020. doi:10.1190/segam2020-3422419.1
SponsorsThis publication is based on work supported by King Abdul-lah University of Science and Technology (KAUST) and the Seismic Modeling and Inversion (SMI) group at KAUST.
PublisherSociety of Exploration Geophysicists
Conference/Event nameSEG International Exposition and 90th Annual Meeting