Two‐phase flow in heterogeneous porous media: 2. Model application

The migration of a dense, nonaqueous phase liquid through heterogeneous porous media is examined using numerical simulation. Laboratory measurements of capillary pressure-saturation curves were performed on samples obtained from a sand aquifer and scaled to permeability to provide a data base of model input parameters. Numerical simulations incorporating 25,200 finite difference cells are carried out in a spatially correlated, random permeability field to illustrate the influence of fluid properties on the migration of a nonwetting liquid below the water table. The simulation results are characterized by spatial moments to reflect the relative degrees of lateral spreading exhibited by the migrating nonwetting body in the presence of lenses of differing permeability. In addition, numerical simulations were performed in a solution domain containing a single lens of lower permeability material in order to examine the local-scale sensitivity to porous media and fluid properties. The results of the study show the migration of a nonwetting liquid to be extremely sensitive to subtle variations in the capillary properties of the porous medium and to be influenced strongly by the fluid physical properties.

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

[2]  John A. Cherry,et al.  Migration of contaminants in groundwater at a landfill: A case study: 1. Groundwater flow and plume delineation , 1983 .

[3]  D. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 1. Approach and overview of plume movement , 1986 .

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

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

[6]  Jack C. Parker,et al.  Finite‐element analysis of multiphase immiscible flow through soils , 1987 .

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

[8]  J. Parker,et al.  An efficient finite element method for modeling multiphase flow , 1989 .

[9]  Donald L. Katz,et al.  Threshold pressure phenomena in porous media , 1968 .

[10]  J. M. Dumoré,et al.  Drainage Capillary-Pressure Functions and the Influence of Connate Water , 1974 .

[11]  David L. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 2. Spatial moments and the advection and dispersion of nonreactive tracers , 1986 .

[12]  R. R. Berg,et al.  Capillary Pressures in Stratigraphic Traps , 1975 .

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

[14]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[15]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 1. Equation Development , 1985 .

[16]  Charles R. Faust,et al.  Simulation of three‐dimensional flow of immiscible fluids within and below the unsaturated zone , 1989 .

[17]  Charles R. Faust,et al.  Transport of Immiscible Fluids Within and Below the Unsaturated Zone: A Numerical Model , 1985 .

[18]  Peter A. Forsyth Simulation of nonaqueous phase groundwater contamination , 1988 .

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

[20]  R. Lenhard,et al.  Measurement and prediction of saturation-pressure relationships in three-phase porous media systems , 1987 .

[21]  J. F. Sykes,et al.  Numerical Modeling of Immiscible Organic Transport at the Hyde Park Landfill , 1986 .