Uncertainty Quantification in Scale‐Dependent Models of Flow in Porous Media

Equations governing flow and transport in randomly heterogeneous porous media are stochastic and scale-dependent. In the Moment Equations (ME) method, exact deterministic equations for the leading moments of state variables are obtained at the same support scale as the governing equations. Computable approximations of the MEs can be derived via perturbation expansion in orders of the standard deviation of the random model parameters. As such, their convergence is guaranteed only for standard deviation smaller than one. Here we consider steady-state saturated flow in a porous medium with random second-order stationary conductivity field. We show it is possible to identify a support scale η*, where the typically employed approximate formulations of ME yield accurate (statistical) moments of a target state variable. Therefore, at support scale η* and larger, ME presents an attractive alternative to slowly convergent Monte Carlo (MC) methods whenever lead-order statistical moments of a target state variable are needed. We also demonstrate that a surrogate model for statistical moments could be constructed from MC simulations at larger support scales and be used to accurately estimate moments at smaller scales, where MC simulations are expensive and the ME method is not applicable.

[1]  E. C. Childs.  Dynamics of fluids in Porous Media , 1973 .

[2]  Sabine Attinger,et al.  Generalized Coarse Graining Procedures for Flow in Porous Media , 2003 .

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

[4]  Alexandre M. Tartakovsky,et al.  Numerical Studies of Three-dimensional Stochastic Darcy’s Equation and Stochastic Advection-Diffusion-Dispersion Equation , 2010, J. Sci. Comput..

[5]  C. Clauser.  Permeability of crystalline rocks , 1992 .

[6]  D. Xiu.  Fast numerical methods for stochastic computations: A review , 2009 .

[7]  Alberto Guadagnini,et al.  Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion , 2003 .

[8]  Alexandre M. Tartakovsky,et al.  Divergence of solutions to solute transport moment equations , 2008 .

[9]  Alexandre M. Tartakovsky,et al.  Immiscible front evolution in randomly heterogeneous porous media , 2003 .

[10]  S. Hassanizadeh,et al.  Upscaling Multiphase Flow in Porous Media , 2005 .

[11]  Erik H. Vanmarcke,et al.  Random Fields: Analysis and Synthesis. , 1985 .

[12]  Daniel M. Tartakovsky,et al.  Uncertainty quantification via random domain decomposition and probabilistic collocation on sparse grids , 2010, J. Comput. Phys..

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

[14]  Emil M. Constantinescu,et al.  Probabilistic Density Function Method for Stochastic ODEs of Power Systems with Uncertain Power Input , 2015, SIAM/ASA J. Uncertain. Quantification.

[15]  Alberto Guadagnini,et al.  EnKF coupled with groundwater flow moment equations applied to Lauswiesen aquifer, Germany , 2015 .

[16]  E. Vanmarcke,et al.  Random Fields: Analysis and Synthesis (Revised and Expanded New Edition) , 2010 .

[17]  S. P. Neuman,et al.  Inverse stochastic moment analysis of steady state flow in randomly heterogeneous media , 2006 .

[18]  S. P. Neuman,et al.  Generalized scaling of permeabilities: Validation and effect of support scale , 1994 .

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

[20]  Daniel M Tartakovsky,et al.  Probability density function method for Langevin equations with colored noise. , 2013, Physical review letters.

[21]  Clayton V. Deutsch,et al.  Geostatistical Software Library and User's Guide , 1998 .

[22]  Zhiming Lu,et al.  A simple kriging method incorporating multiscale measurements in geochemical survey , 2009 .

[23]  Daniel M Tartakovsky,et al.  PDF equations for advective-reactive transport in heterogeneous porous media with uncertain properties. , 2011, Journal of contaminant hydrology.

[24]  S. P. Neuman,et al.  Use of variable-scale pressure test data to estimate the log hydraulic conductivity covariance and dispersivity of fractured granites near Oracle, Arizona , 1988 .

[25]  Daniel M. Tartakovsky,et al.  Localization of Mean Flow and Apparent Transmissivity Tensor for Bounded Randomly Heterogeneous Aquifers , 2002 .

[26]  Guang Lin,et al.  An efficient, high-order probabilistic collocation method on sparse grids for three-dimensional flow and solute transport in randomly heterogeneous porous media , 2009 .

[27]  C. Newell,et al.  The New Potential for Understanding Groundwater Contaminant Transport , 2014, Ground water.

[28]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains: 2. Computational examples , 1999 .

[29]  Daniel M. Tartakovsky,et al.  Nonlocal and localized analyses of conditional mean transient flow in bounded, randomly heterogeneous porous media , 2004 .

[30]  Divergence of solutions to perturbation-based advection–dispersion moment equations , 2011 .

[31]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains: 1. Theory and computational approach , 1999 .

[32]  Alberto Guadagnini,et al.  Comparison of Ensemble Kalman Filter groundwater-data assimilation methods based on stochastic moment equations and Monte Carlo simulation , 2014 .

[33]  Alexandre M. Tartakovsky,et al.  Stochastic analysis of immiscible displacement of the fluids with arbitrary viscosities and its dependence on support scale of hydrological data , 2004 .

[34]  Ram Rajagopal,et al.  Towards optimal allocation of computer resources: Trade-offs between uncertainty quantification, discretization and model reduction , 2013, Environ. Model. Softw..

[35]  Vittorio Di Federico,et al.  Multifaceted nature of hydrogeologic scaling and its interpretation , 2003 .

[36]  S. P. Neuman,et al.  Correction to “Prediction of steady state flow in nonuniform geologic media by conditional moments: Exact nonlocal formalism, effective conductivities, and weak approximation” by Shlomo P. Neuman and Shlomo Orr , 1993 .

[37]  Alberto Guadagnini,et al.  Convergence assessment of numerical Monte Carlo simulations in groundwater hydrology , 2004 .

[38]  S. P. Neuman,et al.  Radial Flow in a Bounded Randomly Heterogeneous Aquifer , 2001 .

[39]  S. P. Neuman,et al.  Three‐dimensional steady state flow to a well in a randomly heterogeneous bounded aquifer , 2003 .

[40]  Alberto Guadagnini,et al.  Data assimilation and parameter estimation via ensemble Kalman filter coupled with stochastic moment equations of transient groundwater flow , 2013 .

[41]  Carolyn Conner Seepersad,et al.  Building Surrogate Models Based on Detailed and Approximate Simulations , 2004, DAC 2004.