Solute transport through fractured media: 1. The effect of matrix diffusion

Solute transport through fractured media is described by numerically combining advective-dispersive transport, which is dominant in the fractures, and diffusive transport, which is usually dominant in the unfractured matrix. Transport is considered in a manner conceptually similar to ‘double-porosity’ or ‘intra-aggregate’ transport models. A finite element model is developed for simulating nonreactive and reactive solute transport by advection, mechanical dispersion, and diffusion in a unidirectional flow field. The effect of the value of the solute diffusion coefficient in the matrix (termed the matrix diffusion coefficient) is illustrated by solute breakthrough curves and concentration profiles in the fracture as well as in the matrix. The illustrated conditions are similar to the laboratory tracer study on fractured till described in the accompanying paper (Grisak et al.). The effects on solute transport of fracture aperture size, water velocity in the fracture, matrix porosity, matrix distribution coefficient, and dispersivity in the fracture are illustrated with breakthrough curves and concentration profiles. The net effect of large matrix diffusion coefficients and/or large distribution coefficients in the matrix is to reduce significantly the effective solute velocity in the fracture. Reduction in the effective solute velocity can also be seen to be possible in materials with low matrix porosities such as crystalline rocks. The aperture size, matrix porosity, matrix diffusion coefficient, and distribution coefficient, all are important in determining the relative amounts of solute transported in the fracture and stored in the matrix. Implications with regard to aquifer recharge, groundwater chemistry, contaminant transport, tracer tests, and groundwater age dating in fractured media are discussed. The numerical model and the laboratory tracer test data provide considerable insight into the processes controlling solute transport in fractured media.

[2]  D. Norton,et al.  Transport phenomena in hydrothermal systems: the nature of porosity , 1977 .

[3]  R. B. Lantz Quantitative Evaluation of Numerical Diffusion (Truncation Error) , 1971 .

[4]  O. Zienkiewicz The Finite Element Method In Engineering Science , 1971 .

[5]  W. D. Kemper,et al.  Contribution of Adsorbed Cations to Diffusion in Clay-Water Systems , 1966 .

[6]  D. B. Oakes,et al.  Prediction of Future Nitrate Concentrations in Ground Watera , 1976 .

[7]  George F. Pinder,et al.  A Galerkin‐finite element simulation of groundwater contamination on Long Island, New York , 1973 .

[8]  W. A. Beetem,et al.  Porosity and Dispersion Constant Calculations for a Fractured Carbonate Aquifer Using the Two Well Tracer Method , 1971 .

[9]  R. E. Phillips,et al.  The diffusion of copper, manganese, and zinc as affected by concentration, clay mineralogy, and associated anions. , 1970 .

[10]  S. Foster The Chalk groundwater tritium anomaly — A possible explanation , 1975 .

[11]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[12]  W. D. Kemper,et al.  Movement of Nutrients to Plant Roots , 1968 .

[13]  R. K. Stoessell,et al.  A nonsteady state method for determining diffusion coefficients in porous media , 1975 .

[14]  Robert W. Gillham,et al.  Finite-element analysis of the transport of water and solutes in title-drained soils , 1979 .

[15]  R. Banks,et al.  A solution of the differential equation of longitudinal dispersion in porous media , 1961 .

[16]  R. Garrels,et al.  DIFFUSION OF IONS THROUGH INTERGRANULAR SPACES IN WATER-SATURATED ROCKS , 1949 .

[17]  William C. Lennox,et al.  Numerical simulation of waste movement in steady groundwater flow systems , 1976 .

[18]  A. Peck Mass transport in porous rocks , 1967 .

[19]  R. Krizek,et al.  DISPERSION OF A CONTAMINANT IN FISSURED ROCK. , 1972 .

[20]  John A. Cherry,et al.  Solute transport through fractured media: 2. Column study of fractured till , 1980 .