1887

Abstract

Summary

In this abstract, we review the gradient-based Markov Chain Monte Carlo (MCMC) and demonstrate its applicability in inferring the uncertainty in seismic inversion. There are many flavours of gradient-based MCMC; here, we will only focus on the Unadjusted Langevin algorithm (ULA) and Metropolis-Adjusted Langevin algorithm (MALA). We propose an adaptive step-length based on the Lipschitz condition within ULA to automate the tuning of step-length and suppress the Metropolis-Hastings acceptance step in MALA. We consider the linear seismic travel-time tomography problem as a numerical example to demonstrate the applicability of both methods.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202010496
2021-10-18
2024-03-29
Loading full text...

Full text loading...

References

  1. Brooks, S., Gelman, A., Jones, G. and Meng, X.L.
    [2011] Handbook of Markov Chain Monte Carlo. CRC press.
    [Google Scholar]
  2. Fichtner, A., Zunino, A. and Gebraad, L.
    [2018] Hamiltonian Monte Carlo solution of tomographic inverse problems.Geophysical Journal International, 216(2), 1344–1363.
    [Google Scholar]
  3. Gazzola, S., Hansen, P.C. and Nagy, J.G.
    [2019] IR Tools: a MATLAB package of iterative regularization methods and large-scale test problems.Numerical Algorithms, 81(3), 773–811.
    [Google Scholar]
  4. Girolami, M. and Calderhead, B.
    [2011] Riemann manifold Langevin and Hamiltonian Monte Carlo methods.Journal of the Royal Statistical Society: Series B (Statistical Methodology), 73(2), 123–214.
    [Google Scholar]
  5. Martin, J., Wilcox, L., Burstedde, C. and Ghattas, O.
    [2012] A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion.SIAM Journal on Scientific Computing, 34(3), A1460–A1487.
    [Google Scholar]
  6. Neal, P. and Roberts, G.O.
    [2006] Optimal Scaling For Partially Updating Mcmc Algorithms.
    [Google Scholar]
  7. Roberts, G.O. and Rosenthal, J.S.
    [1998] Optimal scaling of discrete approximations to Langevin diffusions.
    [Google Scholar]
  8. Simsekli, U., Badeau, R., Cemgil, A.T. and Richard, G.
    [2016] Stochastic Quasi-Newton Langevin Monte Carlo. In: Proceedings of the 33rd International Conference on International Conference on Machine Learning - Volume 48, ICML’16. JMLR.org, 642–651.
    [Google Scholar]
  9. Tarantola, A.
    [2005] Inverse Problem Theory and Methods for Model Parameter Estimation. Society for Industrial and Applied Mathematics.
    [Google Scholar]
  10. Tong, X.T., Morzfeld, M. and Marzouk, Y.M.
    [2019] MALA-within-Gibbs samplers for high-dimensional distributions with sparse conditional structure.
    [Google Scholar]
  11. Welling, M. and Teh, Y.W.
    [2011] Bayesian Learning via Stochastic Gradient Langevin Dynamics. In: Proceedings of the 28th International Conference on International Conference on Machine Learning, ICML’11. Omnipress, USA, 681–688.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.202010496
Loading
/content/papers/10.3997/2214-4609.202010496
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error