High-dimensional Stochastic Inversion via Adjoint Models and Machine Learning

03/16/2018
by   Charanraj A. Thimmisetty, et al.
0

Performing stochastic inversion on a computationally expensive forward simulation model with a high-dimensional uncertain parameter space (e.g. a spatial random field) is computationally prohibitive even with gradient information provided. Moreover, the `nonlinear' mapping from parameters to observables generally gives rise to non-Gaussian posteriors even with Gaussian priors, thus hampering the use of efficient inversion algorithms designed for models with Gaussian assumptions. In this paper, we propose a novel Bayesian stochastic inversion methodology, characterized by a tight coupling between a gradient-based Langevin Markov Chain Monte Carlo (LMCMC) method and a kernel principal component analysis (KPCA). This approach addresses the `curse-of-dimensionality' via KPCA to identify a low-dimensional feature space within the high-dimensional and nonlinearly correlated spatial random field. Moreover, non-Gaussian full posterior probability distribution functions are estimated via an efficient LMCMC method on both the projected low-dimensional feature space and the recovered high-dimensional parameter space. We demonstrate this computational framework by integrating and adapting recent developments such as data-driven statistics-on-manifolds constructions and reduction-through-projection techniques to solve inverse problems in linear elasticity.

READ FULL TEXT

page 28

page 29

page 30

page 32

page 33

page 36

research
11/08/2022

Domain-decomposed Bayesian inversion based on local Karhunen-Loève expansions

In many Bayesian inverse problems the goal is to recover a spatially var...
research
05/26/2021

Cauchy Markov Random Field Priors for Bayesian Inversion

The use of Cauchy Markov random field priors in statistical inverse prob...
research
08/25/2018

Stochastic Collocation with Non-Gaussian Correlated Parameters via a New Quadrature Rule

This paper generalizes stochastic collocation methods to handle correlat...
research
09/19/2021

On the curvatures of Gaussian random field manifolds

Information geometry is concerned with the application of differential g...
research
09/30/2014

Hyper-Spectral Image Analysis with Partially-Latent Regression and Spatial Markov Dependencies

Hyper-spectral data can be analyzed to recover physical properties at la...
research
05/09/2013

Inferring Team Strengths Using a Discrete Markov Random Field

We propose an original model for inferring team strengths using a Markov...

Please sign up or login with your details

Forgot password? Click here to reset