TY - JOUR
T1 - A Bayesian approach to estimate uncertainty for full-waveform inversion using a priori information from depth migration
AU - Zhu, Hejun
AU - Li, Siwei
AU - Fomel, Sergey
AU - Stadler, Georg
AU - Ghattas, Omar
N1 - Publisher Copyright:
© 2016 Society of Exploration Geophysicists.
PY - 2016/9/1
Y1 - 2016/9/1
N2 - Full-waveform inversion (FWI) enables us to obtain highresolution subsurface images; however, estimating model uncertainties associated with this technique is still a challenging problem. We have used a Bayesian inference framework to estimate model uncertainties associated with FWI. The uncertainties were assessed based on an a posteriori covariance operator, evaluated at the maximum a posteriori model. For the prior distribution, we have used a spatially nonstationary covariance operator based on a plane-wave construction with local dips measured from migrated images. Preconditioned frequencydomain FWI was used to estimate the maximum a posteriori model. Efficient manipulation of the posterior covariance was based on a low-rank approximation of the data misfit Hessian preconditioned by the prior covariance operator. The strong decay of the singular values indicated that data were mostly informative about a low-dimensional subspace of model parameters. To reduce computational cost of the randomized singular value decomposition, we have used a Hessian approximation based on point-spread functions. The 2D numerical examples with synthetic data confirmed that the method can effectively estimate uncertainties for FWI. Visual comparisons of random samples drawn from the prior and posterior distributions have allowed us to evaluate model uncertainties. Furthermore, we found out how statistical quantities, such as means and pointwise standard deviation fields, can be efficiently extracted from the prior and posterior distributions. These fields helped us to objectively assess subsurface images provided by FWI.
AB - Full-waveform inversion (FWI) enables us to obtain highresolution subsurface images; however, estimating model uncertainties associated with this technique is still a challenging problem. We have used a Bayesian inference framework to estimate model uncertainties associated with FWI. The uncertainties were assessed based on an a posteriori covariance operator, evaluated at the maximum a posteriori model. For the prior distribution, we have used a spatially nonstationary covariance operator based on a plane-wave construction with local dips measured from migrated images. Preconditioned frequencydomain FWI was used to estimate the maximum a posteriori model. Efficient manipulation of the posterior covariance was based on a low-rank approximation of the data misfit Hessian preconditioned by the prior covariance operator. The strong decay of the singular values indicated that data were mostly informative about a low-dimensional subspace of model parameters. To reduce computational cost of the randomized singular value decomposition, we have used a Hessian approximation based on point-spread functions. The 2D numerical examples with synthetic data confirmed that the method can effectively estimate uncertainties for FWI. Visual comparisons of random samples drawn from the prior and posterior distributions have allowed us to evaluate model uncertainties. Furthermore, we found out how statistical quantities, such as means and pointwise standard deviation fields, can be efficiently extracted from the prior and posterior distributions. These fields helped us to objectively assess subsurface images provided by FWI.
UR - http://www.scopus.com/inward/record.url?scp=84984623467&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84984623467&partnerID=8YFLogxK
U2 - 10.1190/GEO2015-0641.1
DO - 10.1190/GEO2015-0641.1
M3 - Article
AN - SCOPUS:84984623467
SN - 0016-8033
VL - 81
SP - R307-R323
JO - Geophysics
JF - Geophysics
IS - 5
ER -