Hamiltonian Monte Carlo inversion of seismic reflection data in the acoustic approximation

Open access
Date
2021Type
- Other Conference Item
ETH Bibliography
yes
Altmetrics
Abstract
The Hamiltonian Monte Carlo method (HMC) is gaining popularity in the geophysical community to fully address nonlinear inverse problems and related uncertainty quantification. We present here an application of HMC to invert seismic data in the acoustic approximation in the context of reflection seismology. We address a 2-D problem, in the form of a vertical cross section where both source and receivers are located near the surface of the model. To solve the forward problem we utilise the finite-difference method with PML absorbing boundary conditions. The observed data are represented by a set of shotgathers.
The crucial aspect for a successful application of the HMC lies in the capability of performing gradient computations in an efficient manner. To this end, we use the adjont state method to compute the gradient of the misfit functional, which has a computational cost of only about twice that of the forward computation, a very efficient strategy. From the collection of samples characterising the posterior distribution obtained with the HMC, we can derive quantities of interest using statistical analysis and assess uncertainties.
We illustrate an application of this methodology on a synthetic test mimicking the setup encountered in exploration problems. Show more
Permanent link
https://doi.org/10.3929/ethz-b-000518163Publication status
publishedExternal links
Journal / series
EGUspherePages / Article No.
Publisher
CopernicusEvent
Subject
seismic inversion; inverse problem; geophysicsOrganisational unit
03971 - Fichtner, Andreas / Fichtner, Andreas
More
Show all metadata
ETH Bibliography
yes
Altmetrics