Stochastic uncertainty analysis for solute transport in randomly heterogeneous media using a Karhunen‐Loève‐based moment equation approach

[1] A new approach has been developed for solving solute transport problems in randomly heterogeneous media using the Karhunen-Loeve-based moment equation (KLME) technique proposed by Zhang and Lu (2004). The KLME approach combines the Karhunen-Loeve decomposition of the underlying random conductivity field and the perturbative and polynomial expansions of dependent variables including the hydraulic head, flow velocity, dispersion coefficient, and solute concentration. The equations obtained in this approach are sequential, and their structure is formulated in the same form as the original governing equations such that any existing simulator, such as Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems (MT3DMS), can be directly applied as the solver. Through a series of two-dimensional examples, the validity of the KLME approach is evaluated against the classical Monte Carlo simulations. Results indicate that under the flow and transport conditions examined in this work, the KLME approach provides an accurate representation of the mean concentration. For the concentration variance, the accuracy of the KLME approach is good when the conductivity variance is 0.5. As the conductivity variance increases up to 1.0, the mismatch on the concentration variance becomes large, although the mean concentration can still be accurately reproduced by the KLME approach. Our results also indicate that when the conductivity variance is relatively large, neglecting the effects of the cross terms between velocity fluctuations and local dispersivities, as done in some previous studies, can produce noticeable errors, and a rigorous treatment of the dispersion terms becomes more appropriate.

[1]  L. Gelhar,et al.  Transport in three-dimensionally heterogeneous aquifers: 2. Predictions and observations of concentration fluctuations , 1994 .

[2]  Ahmed E. Hassan,et al.  A Monte Carlo assessment of Eulerian flow and transport perturbation models , 1998 .

[3]  Arturo A. Keller,et al.  A stochastic analysis of steady state two‐phase flow in heterogeneous media , 2005 .

[4]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[5]  S. P. Neuman,et al.  Eulerian‐Lagrangian Analysis of Transport Conditioned on Hydraulic Data: 1. Analytical‐Numerical Approach , 1995 .

[6]  Giovanna Darvini,et al.  Nonstationary flow and nonergodic transport in random porous media , 2006 .

[7]  John H. Cushman,et al.  The Physics of Fluids in Hierarchical Porous Media: Angstroms to Miles , 1997 .

[8]  G. Dagan Flow and transport in porous formations , 1989 .

[9]  Dongxiao Zhang,et al.  Conditional simulations of flow in randomly heterogeneous porous media using a KL-based moment-equation approach , 2004 .

[10]  G. Dagan,et al.  A solute flux approach to transport in heterogeneous formations: 1. The general framework , 1992 .

[11]  S. P. Neuman,et al.  Effect of Local Dispersion on Solute Transport in Randomly Heterogeneous Media , 1996 .

[12]  P. Kitanidis,et al.  Advection-diffusion in spatially random flows: Formulation of concentration covariance , 1997 .

[13]  Y. Rubin Applied Stochastic Hydrogeology , 2003 .

[14]  Dongxiao Zhang,et al.  Stochastic uncertainty analysis for unconfined flow systems , 2006 .

[15]  V. Fiorotto,et al.  Dispersion tensor evaluation in heterogeneous media for finite Peclet values , 2000 .

[16]  Roger Ghanem,et al.  Stochastic Finite Element Analysis for Multiphase Flow in Heterogeneous Porous Media , 1998 .

[17]  S. P. Neuman,et al.  A quasi-linear theory of non-Fickian and Fickian subsurface dispersion , 1990 .

[18]  G. Dagan Solute transport in heterogeneous porous formations , 1984, Journal of Fluid Mechanics.

[19]  Dennis McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 1. Unconditional moments , 1989 .

[20]  You‐Kuan Zhang Stochastic Methods for Flow in Porous Media: Coping with Uncertainties , 2001 .

[21]  Shlomo P. Neuman Blueprint for perturbative solution of flow and transport in strongly heterogeneous composite media using fractal and variational multiscale decomposition , 2006 .

[22]  Ahmed E. Hassan,et al.  Eulerian solutions of O(σNv) for the stochastic transport problem for conservative tracers coupled with O(σ4ƒ) solutions for the flow problem in an infinite domain , 1999 .

[23]  E. Frind,et al.  Simulation of contaminant transport in three dimensions: 2. Dimensionality effects , 1987 .

[24]  Dongxiao Zhang,et al.  An efficient, high-order perturbation approach for flow in random porous media via Karhunen-Loève and polynomial expansions , 2004 .

[25]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media : Computational analysis , 2006 .

[26]  S. P. Neuman,et al.  Eulerian‐Lagrangian Theory of transport in space‐time nonstationary velocity fields: Exact nonlocal formalism by conditional moments and weak approximation , 1993 .

[27]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[28]  Vivek Kapoor,et al.  Transport in three-dimensionally heterogeneous aquifers: 1. Dynamics of concentration fluctuations , 1994 .

[29]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[30]  Alberto Guadagnini,et al.  Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media : Theoretical framework , 2006 .

[31]  L. Gelhar Stochastic Subsurface Hydrology , 1992 .

[32]  Dongxiao Zhang,et al.  Stochastic analysis of saturated–unsaturated flow in heterogeneous media by combining Karhunen-Loeve expansion and perturbation method , 2004 .

[33]  Allan L. Gutjahr,et al.  Stochastic analysis of macrodispersion in a stratified aquifer , 1979 .

[34]  E. Frind,et al.  Simulation of contaminant transport in three dimensions: 1. The alternating direction Galerkin Technique , 1987 .

[35]  G. Dagan,et al.  The influence of pore‐scale dispersion on concentration statistical moments in transport through heterogeneous aquifers , 1997 .