Flow and travel time statistics in highly heterogeneous porous media

Abstract

[1] In this paper we present flow and travel time ensemble statistics based on a new simulation methodology, the adaptive Fup Monte Carlo method (AFMCM). As a benchmark case, we considered two-dimensional steady flow in a rectangular domain characterized by multi-Gaussian heterogeneity structure with an isotropic exponential correlation and lnK variance s Y 2 up to 8. Advective transport is investigated using the travel time framework where Lagrangian variables (e.g., velocity, transverse displacement, or travel time) depend on space rather than on time. We find that Eulerian and Lagrangian velocity distributions diverge for increasing lnK variance due to enhanced channeling. Transverse displacement is a nonnormal for all s Y 2 and control planes close to the injection area, but after xI Y = 20 was found to be nearly normal even for high s Y 2 . Travel time distribution deviates from the Fickian model for large lnK variance and exhibits increasing skewness and a power law tail for large lnK variance, the slope of which decreases for increasing distance from the source; no anomalous features are found. Second moment of advective transport is analyzed with respect to the covariance of two Lagrangian velocity variables: slowness and slope which are directly related to the travel time and transverse displacement variance, which are subsequently related to the longitudinal and transverse dispersion. We provide simple estimators for the Eulerian velocity variance, travel time variance, slowness, and longitudinal dispersivity as a practical contribution of this analysis. Both two-parameter models considered (the advection-dispersion equation and the lognormal model) provide relatively poor representations of the initial part of the travel time probability density function in highly heterogeneous porous media. We identify the need for further theoretical and experimental scrutiny of early arrival times, and the need for computing higher-order moments for a more accurate characterization of the travel time probability density function. A brief discussion is presented on the challenges and extensions for which AFMCM is suggested as a suitable approach.Citation: Gotovac, H., V. Cvetkovic, and R. Andricevic (2009), Flow and travel time statistics in highly heterogeneous porous media, Water Resour. Res., 45, W07402,

Concentration statistics of solute transport for the near field zone of an estuary

“…The presented methodology provides a potentially useful tool for the screening assessment, enabling one to estimate the concentration statistics of any conservative substances continuously discharged in the sea by the river's surface water. The developed concept follows the similar methodology already applied in other environmental flows, such as in the atmosphere ( Chatwin and Sullivan, 1990a;Sawford and Sullivan, 1995;Sullivan, 2004 ), and in the groundwater ( Andricevic, 1998;2008;Gotovac et al, 2009 ).…”

Heterogeneity of Infiltration Rates in Alluvial Floodplains as Measured with a Berm Infiltration Technique

“…Linear deposits of coarse-grained sediments with high infiltration rates such as gravel outcrops and macropores create preferential flow paths (Gotovac et al, 2009;Najm et al, 2010). These flow paths can link distant floodplains areas directly to streams.…”

“…For water movement through soil, macropores have been shown to have a large impact on flow and solute transport (Thomas and Phillips, 1979;Fox et al, 2004;Djodjic et al, 2004;Akay and Fox, 2007;Gotovac et al, 2009). Gold and Kellogg (1997) specifically called for the development of unique sampling schemes and simulation models for situations where substantial infiltration occurs through preferential flow pathways, but limited work in monitoring, theoretical model development, and application has been achieved to date.…”

Tracer travel and residence time distributions in highly heterogeneous aquifers: Coupled effect of flow variability and mass transfer

“…Although two-dimensional cases have been extensively studied numerically [e.g., 1,8,48,29,33,16,25,26], solutions for three-dimensional cases are still difficult and rare. Furthermore, mass 5 transfer processes can affect transport significantly; these are known to be ratelimiting due to diffusion, yet their modelling is not straightforward.…”