Main content

Date created: | Last Updated:

: DOI | ARK

Creating DOI. Please wait...

Create DOI

Category: Project

Description: We present an introduction to Hamiltonian Monte Carlo (HMC) sampling of high-dimensional model spaces with focus on linear and weakly nonlinear inverse problems. This includes the theoretical foundations, intuitive examples, as well as applications to linear and nonlinear, adjoint-based traveltime tomography. HMC rests on the construction of an artificial Hamiltonian system where a model is treated as a high-dimensional particle moving along a trajectory in an extended model space. Using derivatives of the forward equations, HMC is able to make long-distance moves from the current towards a new model, thereby promoting sample independence while maintaining high acceptance rates. Though HMC is well-known and appreciated in many branches of the physical and computational sciences, it is only starting to be explored for the solution of geophysical inverse problems. The first part of this article is therefore dedicated to an introduction to HMC using common geophysical terminology. Based on a brief review of equilibrated random walks, we provide a self-contained proof of the HMC algorithm. This leads us to consider simple examples with linear forward problems and Gaussian priors. Having analytical solutions, these problems are not the primary target of Monte Carlo methods. However, they allow us to gain intuition and to study critical tuning issues, the independence of successive samples, the tractable number of dimensions, and limitations of the algorithm. The second part covers applications to linear traveltime tomography, reflector imaging, and nonlinear traveltime tomography. The latter rests on a fast-sweeping method combined with adjoint techniques for the computation of derivatives. Using HMC allows us to solve the tomographic problem fully probabilistically for several thousand unknowns, and without the requirement of supercomputing resources.

License: Academic Free License (AFL) 3.0

Files

Loading files...

Citation

Tags

Recent Activity

Loading logs...

OSF does not support the use of Internet Explorer. For optimal performance, please switch to another browser.
Accept
This website relies on cookies to help provide a better user experience. By clicking Accept or continuing to use the site, you agree. For more information, see our Privacy Policy and information on cookie use.
Accept
×

Start managing your projects on the OSF today.

Free and easy to use, the Open Science Framework supports the entire research lifecycle: planning, execution, reporting, archiving, and discovery.