The method of distributions for dispersive transport in porous media with uncertain hydraulic properties

Predictions of solute transport in subsurface environments are notoriously unreliable due to aquifer heterogeneity and uncertainty about the values of hydraulic parameters. Probabilistic framework, which treats the relevant parameters and solute concentrations as random fields, allows for quantification of this predictive uncertainty. By providing deterministic equations for either probability density function or cumulative distribution function (CDF) of predicted concentrations, the method of distributions enables one to estimate, e.g., the probability of a contaminant's concentration exceeding a safe dose. We derive a deterministic equation for the CDF of solute concentration, which accounts for uncertainty in flow velocity and initial conditions. The coefficients in this equation are expressed in terms of the mean and variance of concentration. The accuracy and robustness of the CDF equations are analyzed by comparing their predictions with those obtained with Monte Carlo simulations and an assumed beta CDF.

[1]  Daniel M. Tartakovsky,et al.  Assessment and management of risk in subsurface hydrology: A review and perspective , 2013 .

[2]  Daniel M. Tartakovsky,et al.  Perspective on theories of non-Fickian transport in heterogeneous media , 2009 .

[3]  Dennis McLaughlin,et al.  Asymptotic properties of the Eulerian truncation approximation: Analysis of the perfectly stratified transport problem , 2002 .

[4]  D M Tartakovsky,et al.  Cumulative distribution function solutions of advection–reaction equations with uncertain parameters , 2014, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[5]  Daniel M. Tartakovsky,et al.  Probability density functions for passive scalars dispersed in random velocity fields , 2010 .

[6]  John H. Cushman,et al.  Fickian dispersion is anomalous , 2015 .

[7]  Marco Dentz,et al.  Concentration statistics for transport in heterogeneous media due to stochastic fluctuations of the center of mass velocity , 2012 .

[8]  D. Haworth.  Progress in probability density function methods for turbulent reacting flows , 2010 .

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

[10]  T. Lundgren.  Distribution Functions in the Statistical Theory of Turbulence , 1967 .

[11]  Aldo Fiori,et al.  First‐order based cumulative distribution function for solute concentration in heterogeneous aquifers: Theoretical analysis and implications for human health risk assessment , 2014 .

[12]  A. Fiori.  The Lagrangian concentration approach for determining dilution in aquifer transport: Theoretical analysis and comparison with field experiments , 2001 .

[13]  Daniel M. Tartakovsky,et al.  Stochastic analysis of effective rate constant for heterogeneous reactions , 2003 .

[14]  M. Dentz,et al.  Modeling non‐Fickian transport in geological formations as a continuous time random walk , 2006 .

[15]  Kenzi Karasaki,et al.  Probability Density Functions for Solute Transport in Random Field , 2003 .

[16]  Heinz Pitsch,et al.  Hybrid large-eddy simulation/Lagrangian filtered-density-function approach for simulating turbulent combustion , 2005 .

[17]  Andrew J. Majda,et al.  An explicit example with non‐Gaussian probability distribution for nontrivial scalar mean and fluctuation , 1996 .

[18]  Daniel M. Tartakovsky,et al.  Exact PDF equations and closure approximations for advective-reactive transport , 2013, J. Comput. Phys..

[19]  A. Bellin,et al.  On the local concentration probability density function of solutes reacting upon mixing , 2011 .

[20]  Peter K. Kitanidis,et al.  Concentration fluctuations and dilution in aquifers , 1998 .

[21]  Alberto Bellin,et al.  Probability density function of non-reactive solute concentration in heterogeneous porous formations. , 2007, Journal of contaminant hydrology.

[22]  Daniel M. Tartakovsky,et al.  Probability density functions for advective‐reactive transport with uncertain reaction rates , 2009 .

[23]  Daniel M. Tartakovsky,et al.  Extension of “Transient flow in bounded randomly heterogeneous domains: 1, Exact conditional moment equations and recursive approximations” , 1999 .

[24]  Laurent Falk,et al.  A generalized mixing model for initial contacting of reactive fluids , 1994 .

[25]  John H. Cushman,et al.  A primer on upscaling tools for porous media , 2002 .

[26]  C. L. Winter,et al.  A Perturbation Expansion for Diffusion in a Random Velocity Field , 1984 .

[27]  Alexandre M. Tartakovsky,et al.  A Comparison of Closures for Stochastic Advection-Diffusion Equations , 2013, SIAM/ASA J. Uncertain. Quantification.

[28]  M. Dentz,et al.  Mixing and spreading in stratified flow , 2007 .

[29]  Gedeon Dagan,et al.  Theory of Solute Transport by Groundwater , 1987 .