4.6 Article

Full-waveform inversion by informed-proposal Monte Carlo

Journal

GEOPHYSICAL JOURNAL INTERNATIONAL
Volume 230, Issue 3, Pages 1824-1833

Publisher

OXFORD UNIV PRESS
DOI: 10.1093/gji/ggac150

Keywords

Inverse theory; Statistical methods; Waveform inversion; Computational seismology

Funding

  1. Innovation Fund Denmark through the OPTION Project [5184-00025B]

Ask authors/readers for more resources

Utilizing a global proposal distribution informed by the physics of the problem can significantly enhance the performance of the MCMC algorithm, especially when solving highly nonlinear inverse problems with vast model spaces. This improvement mainly manifests in a dramatic reduction in burn-in time and a better ability to explore high-probability regions.
Markov chain Monte Carlo (MCMC) sampling of solutions to large-scale inverse problems is, by many, regarded as being unfeasible due to the large number of model parameters. This statement, however, is only true if arbitrary, local proposal distributions are used. If we instead use a global proposal, informed by the physics of the problem, we may dramatically improve the performance of MCMC and even solve highly nonlinear inverse problems with vast model spaces. We illustrate this by a seismic full-waveform inverse problem in the acoustic approximation, involving close to 10(6) parameters. The improved performance is mainly seen as a dramatic shortening of the burn-in time (the time used to reach at least local equilibrium), but also the algorithm's ability to explore high-probability regions (through more accepted perturbations) is potentially better. The sampling distribution of the algorithm asymptotically converges to the posterior probability distribution, but as with all other inverse methods used to solve highly nonlinear inverse problems we have no guarantee that we have seen all high-probability solutions in a finite number of iterations. On the other hand, with the proposed method it is possible to sample more high-probability solutions in a shorter time, without sacrificing asymptotic convergence. This may be a practical advantage for problems with many parameters and computer-intensive forward calculations.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.6
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available