Fast computation of uncertainty quantification measures in the geostatistical approach to solve inverse problems

Abstract

We consider the computational challenges associated with uncertainty quantification involved in parameter estimation such as seismic slowness and hydraulic transmissivity fields. The reconstruction of these parameters can be mathematically described as Inverse Problems which we tackle using the Geostatistical approach. The quantification of uncertainty in the Geostatistical approach involves computing the posterior covariance matrix which is prohibitively expensive to fully compute and store. We consider an efficient representation of the posterior covariance matrix at the maximum a posteriori (MAP) point as the sum of the prior covariance matrix and a low-rank update that contains information from the dominant generalized eigenmodes of the data misfit part of the Hessian and the inverse covariance matrix. The rank of the low-rank update is typically independent of the dimension of the unknown parameter. The cost of our method scales as $\bigO(m\log m)$ where $m$ dimension of unknown parameter vector space. Furthermore, we show how to efficiently compute measures of uncertainty that are based on scalar functions of the posterior covariance matrix. The performance of our algorithms is demonstrated by application to model problems in synthetic travel-time tomography and steady-state hydraulic tomography. We explore the accuracy of the posterior covariance on different experimental parameters and show that the cost of approximating the posterior covariance matrix depends on the problem size and is not sensitive to other experimental parameters.

Published in Advances in Water Resources

Cited by in corpus (9)

• Efficient generalized Golub-Kahan based methods for dynamic inverse problems
• Randomized Truncated SVD Levenberg-Marquardt Approach to Geothermal Natural State and History Matching
• A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design
• Computational methods for large-scale inverse problems: a survey on hybrid projection methods
• A Data-Scalable Randomized Misfit Approach for Solving Large-Scale PDE-Constrained Inverse Problems
• Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initializaiton
• Fast increased fidelity approximate Gibbs samplers for Bayesian Gaussian process regression
• Low Rank Independence Samplers in Bayesian Inverse Problems
• Uncertainty quantification in large Bayesian linear inverse problems using Krylov subspace methods