Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models

Variably saturated water flow in a dual-porosity medium may be described using two separate flow equations which are coupled by means of a sink source term Γw, to account for the transfer of water between the macropore (or fracture) and soil (or rock) matrix pore systems. In this study we propose a first-order rate expression for Γw, which assumes that water transfer is proportional to the difference in pressure head between the two pore systems. A general expression for the transfer coefficient αw was derived using Laplace transforms of the linearized horizontal flow equation. The value of αw could be related to the size and shape of the matrix blocks (or soil aggregates) and to the hydraulic conductivity Ka of the matrix at the fracture/matrix interface. The transfer term Γw, was evaluated by comparing simulation results with those obtained with equivalent one- and two-dimensional single-porosity flow models. Accurate results were obtained when Ka was evaluated using a simple arithmetic average of the interface conductivities associated with the fracture and matrix pressure heads. Results improved when an empirical scaling coefficient γw was included in αw. A single value of 0.4 for γw was found to be applicable, irrespective of the hydraulic properties or the initial pressure head of the simulated system.

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

[2]  A. Valocchi,et al.  Constraints on the Validity of Equilibrium and First-Order Kinetic Transport Models in Structured Soils , 1986 .

[3]  Anders Rasmuson,et al.  Exact solution of a model for diffusion in particles and longitudinal dispersion in packed beds , 1980 .

[4]  Ronald C. Dykhuizen,et al.  Transport of solutes through unsaturated fractured media , 1987 .

[5]  R. Wagenet,et al.  Simulation of water and chemicals in macropore soils Part 1. Representation of the equivalent macropore influence and its effect on soilwater flow , 1992 .

[6]  A. Valocchi.  Use of temporal moment analysis to study reactive solute transport in aggregated porous media , 1990 .

[7]  E. A. Klavetter,et al.  A continuum model for water movement in an unsaturated fractured rock mass , 1988 .

[8]  E. J. Tyler,et al.  Solute Movement in Structured Soils: Two-Region Model with Small Interaction1 , 1981 .

[9]  M. Brusseau,et al.  Modeling solute transport in structured soils: a review , 1990 .

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

[11]  G. Vachaud,et al.  Solute Transfer, with Exchange between Mobile and Stagnant Water, through Unsaturated Sand1 , 1977 .

[12]  A Green‐Ampt Model of infiltration in a cracked soil , 1984 .

[13]  Gudmundur S. Bodvarsson,et al.  Absorption of water into porous blocks of various shapes and sizes , 1990 .

[14]  S. R. Workman,et al.  PREFLO: A WATER MANAGEMENT MODEL CAPABLE OF SIMULATING PREFERENTIAL FLOW , 1990 .

[15]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[16]  E. A. Sudicky,et al.  The Laplace Transform Galerkin technique for efficient time-continuous solution of solute transport in double-porosity media , 1990 .

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

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

[19]  Keith Beven,et al.  On the Variation of Infiltration Into a Homogeneous Soil Matrix Containing a Population of Macropores , 1986 .

[20]  Malcolm R. Davidson.  Numerical calculation of saturated-unsaturated infiltration in a cracked soil , 1985 .

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

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

[23]  G. Bodvarsson,et al.  An approximate solution for one‐dimensional absorption in unsaturated porous media , 1989 .

[24]  C. Rappoldt.  THE APPLICATION OF DIFFUSION MODELS TO AN AGGREGATED SOIL , 1990 .

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

[26]  M. V. Genuchten,et al.  Models for simulating salt movement in aggregated field soils , 1986 .

[27]  Per-Erik Jansson,et al.  Modelling water and solute transport in macroporous soil. I. Model description and sensitivity analysis , 1991 .

[28]  R. C. Dykhuizen,et al.  A new coupling term for dual‐porosity models , 1990 .

[29]  A. W. Warrick,et al.  Numerical approximations of darcian flow through unsaturated soil , 1991 .

[30]  J. S. Wang,et al.  Flow and Transport in Fractured Rocks , 1991 .

[31]  Gudmundur S. Bodvarsson,et al.  Integral Method Solution for Diffusion into a Spherical Block , 1989 .

[32]  Bernd Diekkrüger,et al.  BIMODAL POROSITY AND UNSATURATED HYDRAULIC CONDUCTIVITY , 1991 .

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

[34]  Steven G. Thoma,et al.  Impact of fracture coatings on fracture/matrix flow interactions in unsaturated, porous media , 1992 .

[35]  J. Bouma,et al.  A simulation model for infiltration into cracked dry soil , 1980 .

[36]  R. Wagenet,et al.  Simulation of water and chemicals in macropore soils Part 2. Application of linear filter theory , 1992 .

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

[38]  Tomas Vogel.  SWMII-Numerical model of two-dimensional flow in a variably saturated porous medium. , 1988 .