Relative permeability characteristics necessary for simulating DNAPL infiltration, redistribution, and immobilization in saturated porous media

[1] This study presents a relative permeability-saturation (kr-S) constitutive model that incorporates the critical phenomena necessary for simulating the rates of nonwetting fluid infiltration, redistribution, and immobilization in a saturated porous medium. To develop a model validation data set, the migration of a dense, nonaqueous phase liquid (DNAPL) pool within a one-dimensional, 1 m tall, saturated sandpack was monitored under alternating drainage and imbibition conditions. A light transmission/image analysis system, able to distinguish between connected-phase and residual nonwetting phase (NWP) in the apparatus, measured the elevation of the top of the connected-phase DNAPL pool as a function of time. Light transmission calibration curves, correlating fluid saturation to transmitted color at the macroscopic scale, were found to exhibit a functional dependence on saturation history that must be taken into account. Applying the calibration curves to captured images of the experiment provided a continuous sequence of fluid saturation profiles. Numerical simulations of the bench-scale experiment, using model parameters measured independently at the macroscopic scale, predict within measurement uncertainty the observed timescales of DNAPL migration and immobilization. Additional simulations reveal that model validation for imbibition processes depends on properly accounting for NWP kr-S hysteresis, including imbibition function curvature and the abrupt extinction of NWP relative permeability.

[1]  J. Gerhard,et al.  Capillary pressure characteristics necessary for simulating DNAPL infiltration, redistribution, and immobilization in saturated porous media , 2003 .

[2]  J. Høst-Madsen,et al.  Laboratory and numerical investigations of immiscible multiphase flow in soil , 1992 .

[3]  J. H. Dane,et al.  Light Nonaqueous‐Phase Liquid Movement in a Variably Saturated Sand , 1997 .

[4]  William G. Gray,et al.  Thermodynamic basis of capillary pressure in porous media , 1993 .

[5]  Jacob H. Dane,et al.  Simultaneous measurement of capillary pressure, saturation, and effective permeability of immiscible liquids in porous media , 1998 .

[6]  The influence of capillarity in numerical modeling of organic liquid redistribution in two-phase systems , 1998 .

[7]  M. Oostrom,et al.  Comparison of relative permeability-saturation-pressure parametric models for infiltration and redistribution of a light nonaqueous-phase liquid in sandy porous media , 1998 .

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

[9]  J. Pankow,et al.  Dense chlorinated solvents in porous and fractured media : model experiments , 1988 .

[10]  Markus Hilpert,et al.  Investigation of the residual–funicular nonwetting-phase-saturation relation , 2000 .

[11]  Gary R. Olhoeft,et al.  Observed Migration of a Controlled DNAPL Release by Geophysical Methods , 1995 .

[12]  Ulrich Fischer,et al.  Experimental and Numerical Investigation of Soil Vapor Extraction , 1996 .

[13]  F. Dullien.  Porous Media: Fluid Transport and Pore Structure , 1979 .

[14]  D. Stonestrom,et al.  Air permeability and trapped‐air content in two soils , 1989 .

[15]  J. Gerhard,et al.  Influence of constitutive model parameters on the predicted migration of DNAPL in heterogeneous porous media , 2003 .

[16]  F. Brissaud,et al.  Evolution of Capillarity and Relative Permeability Hysteresis , 1972 .

[17]  K. Rathfelder,et al.  Flow and entrapment of dense nonaqueous phase liquids in physically and chemically heterogeneous aquifer formations , 1998 .

[18]  J. H. Dane,et al.  Measurement and simulation of one-dimensional transient three phase flow for monotonic liquid drainage , 1988 .

[19]  J. H. Dane,et al.  Movement and Remediation of Trichloroethylene in a Saturated Heterogeneous Porous Medium 1. Spill Behavior and initial dissolution , 1999 .

[20]  Jason I. Gerhard,et al.  Site‐Specific Design for Dual Phase Recovery and Stabilization of Pooled DNAPL , 2001 .

[21]  Linda M. Abriola,et al.  The influence of field-scale heterogeneity on the infiltration and entrapment of dense nonaqueous phase liquids in saturated formations , 2000 .

[22]  Jack C. Parker,et al.  A model for hysteretic constitutive relations governing multiphase flow: 1. Saturation-pressure relations , 1987 .

[23]  A. Turgeon.  An application of parametric mixed‐integer linear programming to hydropower development , 1987 .

[24]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[25]  David Redman,et al.  A Field Experiment to Study the Behavior of Tetrachloroethylene Below the Water Table: Spatial Distribution of Residual and Pooled DNAPL , 1993 .

[26]  N. T. Burdine.  Relative Permeability Calculations From Pore Size Distribution Data , 1953 .

[27]  David B. McWhorter,et al.  Exact integral solutions for two‐phase flow , 1990 .

[28]  N. Wardlaw,et al.  Mercury Capillary Pressure Curves and the Intepretation of Pore Structure and Capillary Behaviour in Reservoir Rocks , 1976 .

[29]  A. Corey.  Mechanics of Immiscible Fluids in Porous Media , 1986 .

[30]  Jason I. Gerhard,et al.  The Influence of Waterflood Design on the Recovery of Mobile DNAPLs , 1998 .

[31]  Mark D. White,et al.  Modeling fluid flow and transport in variably saturated porous media with the STOMP simulator. 2. Verification and validation exercises , 1995 .

[32]  U. Fischer,et al.  Modeling nonwetting-phase relative permeability accounting for a discontinuous nonwetting phase , 1997 .

[33]  Tammo S. Steenhuis,et al.  Visualization by light transmission of oil and water contents in transient two-phase flow fields , 1998 .

[34]  L. E. Scriven,et al.  Displacement of residual nonwetting fluid from porous media , 1981 .

[35]  Avery H. Demond,et al.  AN EXAMINATION OF RELATIVE PERMEABILITY RELATIONS FOR TWO-PHASE FLOW IN POROUS MEDIA , 1987 .

[36]  J. Stedinger,et al.  Algorithms for Optimizing Hydropower System Operation , 1985 .

[37]  R. A. Raschke,et al.  Estimation of Absolute Permeability From Capillary Pressure Measurements , 1987 .

[38]  Avery H. Demond,et al.  Estimation of two‐phase relative permeability relationships for organic liquid contaminants , 1993 .

[39]  N. T. Hoa.  A new method allowing the measurement of rapid variations of the water content in sandy porous media , 1981 .

[40]  Bernard H. Kueper,et al.  Experimental observations of multiphase flow in heterogeneous porous media , 1989 .

[41]  Tammo S. Steenhuis,et al.  Wetting front instability: 2. Experimental determination of relationships between system parameters and two‐dimensional unstable flow field behavior in initially dry porous media , 1989 .

[42]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[43]  Mark D. White,et al.  Modeling surfactant-enhanced nonaqueous-phase liquid remediation of porous media , 1998 .

[44]  F. Dullien.  3 – Pore Structure , 1979 .

[45]  Ulrich Fischer,et al.  A comparison of relative nonwetting‐phase permeability models , 1999 .

[46]  Vincent C. Tidwell,et al.  X ray and visible light transmission for laboratory measurement of two‐dimensional saturation fields in thin‐slab systems , 1994 .

[47]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 2. Numerical Simulation , 1985 .

[48]  D. R. Nielsen,et al.  A consistent set of parametric models for the two‐phase flow of immiscible fluids in the subsurface , 1989 .

[49]  W. Anderson.  Wettability literature survey - Part 5: The effects of wettability on relative permeability , 1987 .

[50]  Avery H. Demond,et al.  Effect of interfacial forces on two-phase capillary pressure-saturation relationships , 1991 .

[51]  R. A. Morse,et al.  Experimental Investigation of Factors Affecting Laboratory Relative Permeability Measurements , 1951 .

[52]  Tammo S. Steenhuis,et al.  MECHANISM FOR FINGER PERSISTENCE IN HOMOGENEOUS, UNSATURATED, POROUS MEDIA: THEORY AND VERIFICATION , 1989 .

[53]  J. Parker,et al.  Comparing Simulated and Experimental Hysteretic Two‐Phase Transient Fluid Flow Phenomena , 1991 .

[54]  J. H. Dane,et al.  Infiltration and redistribution of perchloroethylene in partially saturated, stratified porous media , 1998 .

[55]  Jack C. Parker,et al.  A model for hysteretic constitutive relations governing multiphase flow: 2. Permeability‐saturation relations , 1987 .

[56]  Jonathan F. Sykes,et al.  Laboratory and model simulations of a LNAPL spill in a variably-saturated sand, 1. Laboratory experiment and image analysis techniques , 1994 .

[57]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 1. Model development , 1991 .

[58]  Carlon S. Land,et al.  Calculation of Imbibition Relative Permeability for Two- and Three-Phase Flow From Rock Properties , 1968 .

[59]  Cass T. Miller,et al.  Multiphase flow and transport modeling in heterogeneous porous media: challenges and approaches , 1998 .

[60]  Karsten H. Jensen,et al.  Experimental study of movement and distribution of dense organic contaminants in heterogeneous aquifers , 1995 .

[61]  T. Ertekin,et al.  The effect of the capillary number and its constituents on two-phase relative permeability curves , 1983 .

[62]  Alain Bourgeat,et al.  TWO-PHASE FLOW IN HETEROGENEOUS POROUS MEDIA , 1991 .

[63]  J. H. Henderson,et al.  Imbibition Relative Permeability in Unconsolidated Porous Media , 1962 .