Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
International Journal for Uncertainty Quantification
IF: 0.967 5-Year IF: 1.301 SJR: 0.531 SNIP: 0.8 CiteScore™: 1.52

ISSN Print: 2152-5080
ISSN Online: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2019028759
Forthcoming Article

Variable-separation based iterative ensemble smoother for Bayesian inverse problems in anomalous diffusion reaction models

Lijian Jiang
Tongji University
Yuming Ba
Hunan University
Na Ou
Hunan University

ABSTRACT

Iterative ensemble smoother (IES) has been widely used to estimate parameters and states of dynamic models where the data is 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 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 functions 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 uncertainty quantification. 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 o