Bayesian approach to quantify parameter uncertainty and impacts on predictive flow and mass transport in heterogeneous aquifer

Groundwater flow and mass transport predictions are subjected to uncertainty due to heterogeneity of hydraulic conductivity, whose variability in space is considerably higher than that of other hydraulic properties relevant to groundwater flow. To characterize the distribution of hydraulic conductivity, random space function (RSF) is often used. The Bayesian approach was applied to quantitatively study the effect of parameter uncertainty in RSF on a hypothetical two-dimensional uniform groundwater flow and mass transport. Specifically, the parameter uncertainty transmitted to macrodispersion in mass transport model was also inferred. The results showed that the posterior probability distributions of parameters were updated after Bayesian inference. The numerical experiments indicated that the overall predictive uncertainty was increased with simulating time along the flow direction. As to the relative contribution of the two types of uncertainty, it indicated that parametric uncertainty was a little more important than stochastic uncertainty for the predictive uncertainty of hydraulic head. When the uncertainty of hydraulic head as well as macrodispersion was transported to mass transport model, a much bigger contribution of stochastic uncertainty was observed. Therefore, parametric uncertainty should not be neglected during the process of subsurface simulation.

[1]  John L. Wilson,et al.  Two-Dimensional Plume in Uniform Ground-Water Flow , 1978 .

[2]  Hesham M. Bekhit,et al.  Using Markov Chain Monte Carlo to quantify parameter uncertainty and its effect on predictions of a groundwater flow model , 2009, Environ. Model. Softw..

[3]  William W.-G. Yeh,et al.  Aquifer parameter identification with optimum dimension in parameterization , 1981 .

[4]  J. Gómez-Hernández,et al.  Uncertainty assessment and data worth in groundwater flow and mass transport modeling using a blocking Markov chain Monte Carlo method. , 2009 .

[5]  Philippe Renard,et al.  Issues in characterizing heterogeneity and connectivity in non-multiGaussian media , 2008 .

[6]  Luk Peeters,et al.  Application of a multimodel approach to account for conceptual model and scenario uncertainties in groundwater modelling , 2010 .

[7]  G. Zeng,et al.  Uncertainty Analysis of Stochastic Solute Transport in a Heterogeneous Aquifer , 2009 .

[8]  E. Sudicky A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process , 1986 .

[9]  Stochastic delineation of capture zones: classical versus Bayesian approach , 2003 .

[10]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[11]  N. Sun Inverse problems in groundwater modeling , 1994 .

[12]  George E. P. Box,et al.  Bayesian Inference in Statistical Analysis: Box/Bayesian , 1992 .

[13]  G. Zeng,et al.  Spatial analysis of human health risk associated with ingesting manganese in Huangxing Town, Middle China. , 2009, Chemosphere.

[14]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[15]  Yoram Rubin,et al.  A Geostatistical Approach to the Conditional Estimation of Spatially Distributed Solute Concentration and Notes on the Use of Tracer Data in the Inverse Problem , 1996 .

[16]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[17]  G. Zeng,et al.  Optimal solute transport in heterogeneous aquifer: coupled inverse modelling , 2010 .

[18]  W. Yeh Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .