Abstract:The iterative ensemble smoother (IES) has been widely used to estimate parameters and states of dynamic models where the data are collected at all observation steps simultaneously. A large number of IES ensemble samples may be required in the estimation. This implies that we need to repeatedly compute the forward model corresponding to the ensemble samples. This leads to slow efficiency for large-scale and strongly nonlinear models. To accelerate the posterior inference in the estimation, a low rank approximation using a variable-separation (VS) method is presented to reduce the cost of computing the forward model. It will be efficient to construct a surrogate model based on the low rank approximation, which gives a separated representation of the solution for the stochastic partial differential equations (SPDEs). The separated representation is the product of deterministic basis ftnctions and stochastic basis functions. For the anomalous diffusion reaction equations, the solution of the next moment depends on all of the previous moments, and this causes expensive computation for the Bayesian inverse problem. The presented VS can avoid this process through a few deterministic basis functions. The surrogate model can work well as the iteration moves on because the stochastic basis becomes more accurate when the uncertainty of random parameters decreases. To enhance the applicability in Bayesian inverse problems, we apply the VS-based IES method to complex structure patterns, which can be parameterized by discrete cosine transform (DCT). The post-processing technique based on a regularization method is employed after the iterations to improve the connectivity of the main features. In the paper, we focus on the time fractional diffusion reaction models in porous media and investigate their Bayesian inverse problems using the VS-based IES. A few numerical examples are presented to show the performance of the proposed IES method by taking account of structure inversion in permeability fields, parameters in permeability and reaction fields, and source ftinctions.

Bayesian Identification of Discontinuous Fields with an Ensemble-Based Variable Separation Multiscale Method

Sparse Variational Bayesian Inference for Water Pipeline Systems with Parameter Uncertainties

Robust Mixture Bayesian Latent Variable Regression with Structural Sparsity and Application to Inferential Sensing of Quality Variables

Variable-Separation Based Iterative Ensemble Smoother For Bayesian Inverse Problems In Anomalous Diffusion Reaction Models

A multi-resolution, non-parametric, Bayesian framework for identification of spatially-varying model parameters

Bayesian Multiscale Finite Element Methods. Modeling missing subgrid information probabilistically

A Sparse Grid Based Bayesian Method for Contaminant Source Identification

Multivariate Gaussian Process Regression for Multiscale Data Assimilation and Uncertainty Reduction

An iterative local updating ensemble smoother for high-dimensional inverse modeling with multimodal distributions

A multiscale strategy for Bayesian inference using transport maps

Multiscale Model Reduction for Stochastic Elasticity Problems Using Ensemble Variable-Separated Method

Multiscale model reduction method for Bayesian inverse problems of subsurface flow

A Segmentation Approach for Stochastic Geological Modeling Using Hidden Markov Random Fields

A hierarchical Bayesian modeling framework for identification of Non-Gaussian processes

A Hierarchical Bayesian Framework Embedded with an Improved Orthogonal Series Expansion for Gaussian Processes and Fields Identification

Parallel ensemble Kalman method with total variation regularization for large-scale field inversion

Bayesian data-driven discovery of partial differential equations with variable coefficients

Nonlocal multicontinua with Representative Volume Elements. Bridging separable and non-separable scales

Identification of nonlinear conservation laws for multiphase flow based on Bayesian inversion

Coupling of multiscale and multi-continuum approaches

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems