A Modeling and experimental study of light nonaqueous phase liquid (LNAPL) accumulation in wells and LNAPL recovery from wells

A formulation is presented for numerical modeling of light nonaqueous phase liquid (LNAPL) accumulation in observation wells and LNAPL removal from recovery wells. The formulation includes fluid gravity segregation, changing water and LNAPL levels in the well, and the corresponding changes in fluid storage in the well bore. The method was added to a three-dimensional finite difference model for hysteretic multiphase flow. The model predictions are compared to three-dimensional pilot scale experiments of LNAPL infiltration and recovery. Model predictions matched experimental results well, demonstrating the validity of the model formulation. Model results indicated that soil layering and lack of mechanical equilibrium had a significant influence on LNAPL thicknesses in wells. Thus empirical approaches based on assumptions of soil homogeneity and mechanical equilibrium would not have been appropriate for prediction of LNAPL volume from toluene thicknesses measured in wells.

[1]  David B. McWhorter,et al.  Volume Estimation of Light Nonaqueous Phase Liquids in Porous Media , 1990 .

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

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

[4]  Jonathan F. Sykes,et al.  Compositional simulation of groundwater contamination by organic compounds: 1. Model development and verification , 1993 .

[5]  D. W. Peaceman.  Interpretation of Well-Block Pressures in Numerical Reservoir Simulation(includes associated paper 6988 ) , 1978 .

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

[7]  J. C. Parker,et al.  Estimation of Free Hydrocarbon Volume from Fluid Levels in Monitoring Wells , 1990 .

[8]  Brent E. Sleep,et al.  Simulation of Bioventing for Soil and Ground-Water Remediation , 1996 .

[9]  D. W. Peaceman.  Interpretation of wellblock pressures in numerical reservoir simulation: Part 3; Off-center and multiple wells within a wellblock , 1990 .

[10]  B. Sleep,et al.  The effect of temperature on capillary pressure‐saturation relationships for air‐water and perchloroethylene‐water systems , 1998 .

[11]  E. A. Sudicky,et al.  Influence of Leaky Boreholes on Cross‐Formational Groundwater Flow and Contaminant Transport , 1995 .

[12]  David A. Collins,et al.  Fully Coupled Multiblock Wells in Oil Simulation , 1985 .

[13]  Jagath J. Kaluarachchi,et al.  Effect of subsurface heterogeneity on free-product recovery from unconfined aquifers , 1996 .

[14]  L. Nghiem,et al.  Field-Scale Simulation Of Horizontal Wells , 1992 .

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

[16]  P. Forsyth,et al.  Discrete wellbore simulations of pump-and-treat strategies for remediation of LNAPL-contaminated aquifers , 1998 .

[17]  Yu-Shu Wu,et al.  A consistent approach for applying numerical boundary conditions for multiphase subsurface flow , 1996 .

[18]  Yu-Shu Wu,et al.  A virtual node method for handling well bore boundary conditions in modeling multiphase flow in porous and fractured media. , 2000 .

[19]  Calvin C. Chien,et al.  Modeling wells in variably saturated soil with wellbore fluid gravity segregation , 2000 .

[20]  D. W. Peaceman.  Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability , 1983 .