Adaptive multiscale MCMC algorithm for uncertainty quantification in seismic parameter estimation

Xiaosi Tan, Richard L. Gibson, Wing Tat Leung, Yalchin R. Efendiev

Research output: Contribution to journalArticlepeer-review

Abstract

Formulating an inverse problem in a Bayesian framework has several major advantages (Sen and Stoffa, 1996). It allows finding multiple solutions subject to flexible a priori information and performing uncertainty quantification in the inverse problem. In this paper, we consider Bayesian inversion for the parameter estimation in seismic wave propagation. The Bayes' theorem allows writing the posterior distribution via the likelihood function and the prior distribution where the latter represents our prior knowledge about physical properties. One of the popular algorithms for sampling this posterior distribution is Markov chain Monte Carlo (MCMC), which involves making proposals and calculating their acceptance probabilities. However, for large-scale problems, MCMC is prohibitevely expensive as it requires many forward runs. In this paper, we propose a multilevel MCMC algorithm that employs multilevel forward simulations. Multilevel forward simulations are derived using Generalized Multiscale Finite Element Methods that we have proposed earlier (Efendiev et al., 2013a; Chung et al., 2013). Our overall Bayesian inversion approach provides a substantial speed-up both in the process of the sampling via preconditioning using approximate posteriors and the computation of the forward problems for different proposals by using the adaptive nature of multiscale methods. These aspects of the method are discussed n the paper. This paper is motivated by earlier work of M. Sen and his collaborators (Hong and Sen, 2007; Hong, 2008) who proposed the development of efficient MCMC techniques for seismic applications. In the paper, we present some preliminary numerical results.
Original languageEnglish (US)
Pages (from-to)4665-4669
Number of pages5
JournalSEG Technical Program Expanded Abstracts 2014
Volume33
DOIs
StatePublished - Aug 5 2014
Externally publishedYes

Fingerprint Dive into the research topics of 'Adaptive multiscale MCMC algorithm for uncertainty quantification in seismic parameter estimation'. Together they form a unique fingerprint.

Cite this