Stochastic Analysis of Reactive Solute Transport in Heterogeneous, Fractured Porous Media: A Dual-Permeability Approach

A nonlocal, first-order, Eulerian stochastic theory is developed for reactive chemical transport in a heterogeneous, fractured porous medium. A dual-permeability model is adopted to describe the flow and transport in the medium, where the solute convection and dispersion in the matrix are considered. The chemical is under linear nonequilibrium sorption and first-order degradation. The hydraulic conductivities, sorption coefficients, degradation rates in both fracture and matrix regions, and interregional mass transfer coefficient are all assumed to be random variables. The resultant theory for mean concentrations in both regions is nonlocal in space and time. Under spatial Fourier and temporal Laplace transforms, the mean concentrations are explicitly expressed. The transformed results are then numerically inverted to the real space via Fast Fourier Transform method. The theory developed in this study generalizes the stochastic studies for a reactive chemical transport in a one-domain flow field (Hu et al., 1997a) and in a mobile/immobile flow field (Huang and Hu, 2001). In comparison with one-domain transport, the dual-permeability model predicts a larger second moment in the longitudinal direction, but smaller one in the transverse direction. In addition, various simplification assumptions have been made based on the general solution. The validity of these assumptions has been tested via the spatial moments of the mean concentration in both fracture and matrix regions.

[1]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[2]  M. Brusseau,et al.  Nonideal transport of reactive solutes in heterogeneous porous media: 6. Microscopic and macroscopic approaches for incorporating heterogeneous rate‐limited mass transfer , 2000 .

[3]  K. H. Coats,et al.  Dead-End Pore Volume and Dispersion in Porous Media , 1964 .

[4]  Emil O. Frind,et al.  Contaminant Transport in Fractured Porous Media: Analytical Solution for a Two‐Member Decay Chain in a Single Fracture , 1984 .

[5]  J. H. Cushman,et al.  Hierarchical Approaches to Transport in Heterogeneous Porous Media , 1991 .

[6]  Ralph E. Showalter,et al.  Diffusion models for fractured media , 1990 .

[7]  J. H. Cushman,et al.  A FAST FOURIER TRANSFORM STOCHASTIC ANALYSIS OF THE CONTAMINANT TRANSPORT PROBLEM , 1993 .

[8]  K. Novakowski,et al.  2. A Semianalytical model for the determination of effective diffusion coefficients, porosity, and adsorption , 1996 .

[9]  Alexander Y. Sun,et al.  Stochastic analysis of transient saturated flow through heterogeneous fractured porous media: A double‐permeability approach , 2000 .

[10]  Eva Hakami,et al.  Aperture measurements and flow experiments on a single natural fracture , 1996 .

[11]  Vladimir Cvetkovic,et al.  Transport of reactive tracers in rock fractures , 1999, Journal of Fluid Mechanics.

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

[13]  James O. Duguid,et al.  Flow in fractured porous media , 1977 .

[14]  Mark N. Goltz,et al.  A natural gradient experiment on solute transport in a sand aquifer: 3. Retardation estimates and mass balances for organic solutes , 1986 .

[15]  R. Gillham,et al.  Laboratory and field measurements of non-equilibrium transport in the Borden aquifer, Ontario, Canada , 1992 .

[16]  Peter K. Kitanidis,et al.  Generalized Taylor-Aris moment analysis of the transport of sorbing solutes through porous media with spatially-periodic retardation factor , 1992 .

[17]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[18]  Luis Moreno,et al.  Analysis of a long-term pumping and tracer test using the channel network model , 1998 .

[19]  Bill X. Hu,et al.  Nonlocal nonreactive transport in heterogeneous porous media with interregional mass diffusion , 2000 .

[20]  P. Witherspoon,et al.  Porous media equivalents for networks of discontinuous fractures , 1982 .

[21]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[22]  Todd Arbogast,et al.  Derivation of the double porosity model of single phase flow via homogenization theory , 1990 .

[23]  John F. Pickens,et al.  Solute transport through fractured media: 1. The effect of matrix diffusion , 1980 .

[24]  Emil O. Frind,et al.  Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures , 1982 .

[25]  Z. J. Kabala,et al.  A stochastic model of reactive solute transport with time‐varying velocity in a heterogeneous aquifer , 1991 .

[26]  Stephen R. Brown.  Simple mathematical model of a rough fracture , 1995 .

[27]  Stephen R. Brown,et al.  Correlation between the surfaces of natural rock joints , 1986 .

[28]  Allen F. Moench,et al.  Double‐Porosity Models for a Fissured Groundwater Reservoir With Fracture Skin , 1984 .

[29]  E. Sudicky,et al.  The interpretation of a tracer experiment conducted in a single fracture under conditions of natural groundwater flow , 1999 .

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

[31]  P. Maloszewski,et al.  Solute Transport Processes in a Highly Permeable Fault Zone of Lindau Fractured Rock Test Site (Germany) , 1998 .

[32]  J. V. Tracy,et al.  Flow through fractures , 1981 .

[33]  Piotr Maloszewski,et al.  On the theory of tracer experiments in fissured rocks with a porous matrix , 1985 .

[34]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[35]  Horst H. Gerke,et al.  Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models , 1993 .

[36]  John H. Cushman,et al.  Nonlocal Reactive Transport with Physical and Chemical Heterogeneity: Linear Nonequilibrium Sorption with Random Kd , 1995 .

[37]  J. H. Cushman,et al.  Nonlocal reactive transport with physical, chemical, and biological heterogeneity , 1997 .

[38]  J. Noorishad,et al.  An upstream finite element method for solution of transient transport equation in fractured porous media , 1982 .

[39]  John H. Cushman,et al.  Nonlocal Reactive Transport with Physical and Chemical Heterogeneity: Localization Errors , 1995 .

[40]  P. A. Witherspoon,et al.  Hydromechanical behavior of a deformable rock fracture subject to normal stress , 1981 .

[41]  Cass T. Miller,et al.  Sorption of hydrophobic organic pollutants in saturated soil systems , 1986 .

[42]  J. Gale,et al.  Comparison of coupled fracture deformation and fluid flow models with direct measurements of fracture pore structure and stress-flow properties , 1987 .

[43]  W. W. Wood,et al.  Large-Scale Natural Gradient Tracer Test in Sand and Gravel, , 1991 .

[44]  L. Gelhar,et al.  Stochastic analysis of sorption macrokinetics in heterogeneous aquifers , 1996 .

[45]  G. Dagan,et al.  Spatial moments of a kinetically sorbing solute plume in a heterogeneous aquifer , 1993 .

[46]  S. P. Neuman,et al.  Reconstruction and geostatistical analysis of multiscale fracture apertures in a large block of welded tuff , 1992 .

[47]  M. Celia,et al.  Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts, 2, Analysis of spatial moments for a nonreactive tracer , 1991 .

[48]  Robert Bibby,et al.  Mass transport of solutes in dual‐porosity media , 1981 .

[49]  E. Sudicky,et al.  Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media , 1996 .

[50]  Bill X. Hu,et al.  Nonlocal reactive transport in heterogeneous dual‐porosity media with rate‐limited sorption and interregional mass diffusion , 2001 .

[51]  M. V. Genuchten,et al.  A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media , 1993 .

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