Characterization of the velocity field organization in heterogeneous media by conditional correlation

The purpose of the present work is to quantify the correlation structure of simulated velocity fields in heterogeneous permeability fields and to discuss how to represent it in upscaled transport models. We investigate the velocity field correlation structure for multinormal log permeability fields. The simulated velocity distributions are analyzed in a Lagrangian framework, i.e., along the particles' paths.To quantify the different spatial organization of low- and high-velocity zones, we condition the estimated velocity correlation length and time on the initial particle velocity. The velocity correlation length is found to increase with the initial particle velocity, following a power law. Such an effect is likely due to the channeling of high-velocity zones, which implies that particles keep memory of their initial velocity over longer distances for high initial velocities than for low initial velocities. Two distinct regimes are identified for the velocity correlation time. For low initial particle velocity the correlation time is controlled by the large time needed to escape from the low- velocity zones. For high initial particle velocity it is controlled by the large time needed for particles to sample the whole velocity field, in particular low- velocity zones. One of the consequences of these results is that for such velocity fields the nonlinear dependence of both the correlation length and time on the particle initial velocity restricts the se of spatial or temporal Markovian assumptions for modeling velocity transitions in effective transport models.

[1]  J. Klafter,et al.  Transport aspects in anomalous diffusion: Lévy walks. , 1989, Physical review. A, General physics.

[2]  W. Kinzelbach,et al.  The Random Walk Method in Pollutant Transport Simulation , 1988 .

[3]  Brian Berkowitz,et al.  Theory of anomalous chemical transport in random fracture networks , 1998 .

[4]  C. Tsang,et al.  Flow channeling in heterogeneous fractured rocks , 1998 .

[5]  D. W. Pollock Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[6]  J. Bouchaud,et al.  Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications , 1990 .

[7]  D. McLaughlin,et al.  An efficient multivariate random field generator using the fast Fourier transform , 1998 .

[8]  C.A.J. Appelo,et al.  FLOW AND TRANSPORT , 2004 .

[9]  Jesús Carrera,et al.  On the relationship between indicators of geostatistical, flow and transport connectivity , 2005 .

[10]  Virgilio Fiorotto,et al.  Solute transport in highly heterogeneous aquifers , 1998 .

[11]  L. Gelhar,et al.  Field study of dispersion in a heterogeneous aquifer: 2. Spatial moments analysis , 1992 .

[12]  G. Matheron,et al.  Is transport in porous media always diffusive? A counterexample , 1980 .

[13]  Hans-Peter Scheffler,et al.  Governing equations and solutions of anomalous random walk limits. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[15]  D. Benson,et al.  Application of a fractional advection‐dispersion equation , 2000 .

[16]  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 .

[17]  Sean Andrew McKenna,et al.  On the late‐time behavior of tracer test breakthrough curves , 2000 .

[18]  J. Klafter,et al.  The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics , 2004 .

[19]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[20]  M. Shlesinger,et al.  Random walks with infinite spatial and temporal moments , 1982 .

[21]  Ralf Metzler,et al.  Physical pictures of transport in heterogeneous media: Advection‐dispersion, random‐walk, and fractional derivative formulations , 2002, cond-mat/0202327.

[22]  A. Rinaldo,et al.  Simulation of dispersion in heterogeneous porous formations: Statistics, first‐order theories, convergence of computations , 1992 .

[23]  S. Gorelick,et al.  Rate‐limited mass transfer or macrodispersion: Which dominates plume evolution at the macrodispersion experiment (MADE) site? , 2000 .

[24]  Allen M. Shapiro,et al.  Interpreting tracer breakthrough tailing from different forced‐gradient tracer experiment configurations in fractured bedrock , 2002 .

[25]  Gaisheng Liu,et al.  Limits of applicability of the advection‐dispersion model in aquifers containing connected high‐conductivity channels , 2004 .

[26]  J. Gómez-Hernández,et al.  To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology , 1998 .

[27]  Brian Berkowitz,et al.  Time behavior of solute transport in heterogeneous media: transition from anomalous to normal transport , 2003 .

[28]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[29]  J. Klafter,et al.  The random walk's guide to anomalous diffusion: a fractional dynamics approach , 2000 .

[30]  R. Ababou,et al.  Numerical simulation of three-dimensional saturated flow in randomly heterogeneous porous media , 1989 .

[31]  Y. Rubin Stochastic modeling of macrodispersion in heterogeneous porous media , 1990 .

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

[33]  Graham E. Fogg,et al.  Role of Molecular Diffusion in Contaminant Migration and Recovery in an Alluvial Aquifer System , 2001 .

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