SciTech Connect

Title: An adaptive sparse-grid high-order stochastic collocation method for Bayesian inference in groundwater reactive transport modeling

An adaptive sparse-grid high-order stochastic collocation method for Bayesian inference in groundwater reactive transport modeling Although Bayesian analysis has become vital to the quantification of prediction uncertainty in groundwater modeling, its application has been hindered due to the computational cost associated with numerous model executions needed for exploring the posterior probability density function (PPDF) of model parameters. This is particularly the case when the PPDF is estimated using Markov Chain Monte Carlo (MCMC) sampling. In this study, we develop a new approach that improves computational efficiency of Bayesian inference by constructing a surrogate system based on an adaptive sparse-grid high-order stochastic collocation (aSG-hSC) method. Unlike previous works using first-order hierarchical basis, we utilize a compactly supported higher-order hierar- chical basis to construct the surrogate system, resulting in a significant reduction in the number of computational simulations required. In addition, we use hierarchical surplus as an error indi- cator to determine adaptive sparse grids. This allows local refinement in the uncertain domain and/or anisotropic detection with respect to the random model parameters, which further improves computational efficiency. Finally, we incorporate a global optimization technique and propose an iterative algorithm for building the surrogate system for the PPDF with multiple significant modes. Once the surrogate system is determined, the PPDF can be evaluated by sampling the surrogate more » system directly with very little computational cost. The developed method is evaluated first using a simple analytical density function with multiple modes and then using two synthetic groundwater reactive transport models. The groundwater models represent different levels of complexity; the first example involves coupled linear reactions and the second example simulates nonlinear ura- nium surface complexation. The results show that the aSG-hSC is an effective and efficient tool for Bayesian inference in groundwater modeling in comparison with conventional MCMC sim- ulations. The computational efficiency is expected to be more beneficial to more computational expensive groundwater problems. « less
Authors: ; ;
Publication Date:
OSTI Identifier:1055118
Report Number(s):ORNL/TM-2012/499
KJ0401000; ERKJE45
DOE Contract Number:DE-AC05-00OR22725
Resource Type:Technical Report
Data Type:
Research Org:Oak Ridge National Laboratory (ORNL)
Sponsoring Org:ORNL LDRD Director's R&D; SC USDOE - Office of Science (SC)
Country of Publication:United States
Language:English
Subject: