Informed proposal Monte Carlo. Issue 2 (29th April 2021)

• Record Type: Journal Article
• Authors: Khoshkholgh, Sarouyeh; Zunino, Andrea; Mosegaard, Klaus
• Is Part Of: Geophysical journal international. Volume 226:Issue 2(2021)
• Page Start: 1239
• Page End: 1248
• Publication Date: 2021-04-29
• Subjects: Inverse theory -- Statistical methods -- Waveform inversion -- Computational seismology; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggab173
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Any search or sampling algorithm for solution of inverse problems needs guidance to be efficient. Many algorithms collect and apply information about the problem on the fly, and much improvement has been made in this way. However, as a consequence of the No-Free-Lunch Theorem, the only way we can ensure a significantly better performance of search and sampling algorithms is to build in as much external information about the problem as possible. In the special case of Markov Chain Monte Carlo (MCMC) sampling we review how this is done through the choice of proposal distribution, and we show how this way of adding more information about the problem can be made particularly efficient when based on an approximate physics model of the problem. A highly non-linear inverse scattering problem with a high-dimensional model space serves as an illustration of the gain of efficiency through this approach.