DNAPL source depletion: linking architecture and flux response

The relationship between dense non-aqueous phase liquid (DNAPL) mass reduction and contaminant mass flux was investigated experimentally in four model source zones. The flow cell design for the experiments featured a segmented extraction well that allowed for analysis of spatially resolved flux information. This flux information was coupled with image analysis of the NAPL spatial distribution to investigate the relationship between flux and the up-gradient NAPL architecture. Results indicate that in the systems studied, the relationship between DNAPL mass reduction and contaminant mass flux was primarily controlled by the NAPL architecture. A specific definition of NAPL architecture was employed where the source zone is resolved into a collection of streamtubes with spatial variability in NAPL saturation along each streamtube integrated and transformed into an effective NAPL content for each streamtube. The distribution of NAPL contents among the streamtubes (NAPL architecture) controlled dissolution dynamics. Two simplified models, a streamtube model and an effective Damkohler number model, were investigated for their ability to simulate dissolution dynamics.     

Characterization of pore scale NAPL morphology in homogeneous sands as a function of grain size and NAPL dissolution

In this study, we investigate pore scale morphology of nonaqueous phase liquids (NAPLs) trapped in different pore sizes using tracer techniques. Specific interfacial area and saturation of NAPL trapped in homogeneous sands were measured using the interfacial and partitioning tracer techniques. The observed NAPL-water interfacial areas increased in a log-linear fashion with decreasing sand grain size, but showed no clear trend with residual NAPL saturation formed in the various grain sizes. The measured values were used to calculate the NAPL morphology index, which characterizes the spatial NAPL distribution within the pore space. The NAPL morphology indices, increased exponentially with decreasing grain size, indicating that the NAPL becomes smaller, but more blobs. For a fixed grain size, the specific interfacial area and saturation of the NAPL were measured following changes caused by dissolution using alcohol. The observed interfacial areas showed a decrease linearly as a function of the NAPL saturation.     

Effect of nonionic surfactant partitioning on the dissolution kinetics of residual perchloroethylene in a model porous medium

At concentrations above the critical micelle concentration, surfactants can significantly enhance the solubilization of residual nonaqueous phase liquids (NAPL) and, for this reason, are the focus of research on surfactant-enhanced aquifer remediation (SEAR). As a consequence of their amphiphilic nature, surfactants may also partition to various extents between the organic and aqueous phases, thereby affecting SEAR performance. We report here on the observation and analysis of the effect of surfactant partitioning on the dissolution kinetics of residual perchloroethylene (PCE) by aqueous solutions (1000 mg/L) of the non-ionic surfactant Triton X-100 in a model porous medium. For this fluid system, batch equilibration experiments showed that the surfactant partitions strongly into the NAPL (NAPL-water partition coefficient equal to 12.5). Dynamic interfacial tension (IFT) measurements were employed to study surfactant diffusion and interfacial adsorption. The dynamic IFT measurements were consistent with partitioning of the surfactant between the two liquid phases. PCE dissolution experiments, conducted in a transparent glass micromodel using an aqueous surfactant solution, were contrasted to experiments using clean water. Surfactant partitioning was observed to delay significantly the onset of micellar solubilization of PCE, an observation reproduced by a numerical model. This effect is attributed to the reduction of surfactant concentration in the immediate vicinity of the NAPL-water interface, which accompanies transport of the surfactant into the NAPL. Accordingly, it is suggested that both the rate and the extent of diffusion of the surfactant into the NAPL affect the onset of and the driving force for micellar solubilization. While many surfactants do not readily partition in NAPL, this possibility must be considered when selecting non-ionic surfactants for the enhanced solubilization of residual chlorinated solvents in porous media.     

Modelling dissolution of multi-component non-aqueous phase liquids using Pareto distributions

A distributed-site non-aqueous phase liquid (NAPL) dissolution model that can capture the physical process of dissolution is presented. The model uses a Pareto distribution of effective interfacial area and a correlation for the mass transfer. Model parameters are obtained using a formal constrained non-linear optimisation technique. The model is applied on effluent data from two soil column experiments. The results of the Pareto distributed site model, with only two parameters to fit, were compared to two-site and multi-site models as well as the beta distributed site model. The Pareto distributed-site model produces equal or smaller error estimates than the other models.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

Secondary imbibition in NAPL-invaded mixed-wet sediments

A previously developed pore network model is used here to study the spontaneous and forced secondary imbibition of a NAPL-invaded sediment, as in the displacement of NAPL by waterflooding a mixed-wet soil. We use a 3D disordered pore network with a realistic representation of pore geometry and connectivity, and a quasi-static displacement model that fully describes the pore-scale physics. After primary drainage (NAPL displacing water) up to a maximum capillary pressure, we simulate secondary imbibition (water displacing NAPL). We conduct a parametric study of imbibition by varying systematically the controlling parameters: the advancing contact angles, the fraction of NAPL-wet pores, the interfacial tension, and the initial water saturation. Once the secondary imbibition is completed, the controlling displacement mechanisms, capillary pressures, relative permeabilities, and trapped NAPL saturations are reported. It is assumed that NAPL migrates into an initially strongly water-wet sediment, i.e., the receding contact angles are very small. However, depending on the surface mineralogy and chemical compositions of the immiscible fluid phases, the wettability of pore interiors is altered while the neighborhoods of pore corners remain strongly water-wet-resulting in a mixed-wet sediment. Here, we compare three different levels of wettability alteration: water-wet (advancing contact angles (20 degrees to 55 degrees), intermediate-wet (55 degrees to 120 degrees), and NAPL-wet (120 degrees to 155 degrees). The range of advancing contact angles and the fraction of NAPL-wet pores have dramatic effects on the NAPL-water capillary pressures and relative permeabilities. The spatially inhomogeneous interfacial tension has a minor impact on the trapped NAPL saturation and relative permeability to NAPL, and a slight effect on the relative permeability to water. The initial water saturation has a slight effect on the two-phase flow characteristics of water-wet sediments; however, with more NAPL-wet pores in the sediment, it starts to have a profound effect on the water and NAPL relative permeabilities.     

A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model

A pore network model with cubic chambers and rectangular tubes was used to estimate the nonaqueous phase liquid (NAPL) dissolution rate coefficient, Kdissai, and NAPL/water total specific interfacial area, ai. Kdissai was computed as a function of modified Peclet number (Pe') for various NAPL saturations (SN) and ai during drainage and imbibition and during dissolution without displacement. The largest contributor to ai was the interfacial area in the water-filled corners of chambers and tubes containing NAPL. When Kdissai was divided by ai, the resulting curves of dissolution coefficient, Kdiss versus Pe' suggested that an approximate value of Kdiss could be obtained as a weak function of hysteresis or SN. Spatially and temporally variable maps of Kdissai calculated using the network model were used in field-scale simulations of NAPL dissolution. These simulations were compared to simulations using a constant value of Kdissai and the empirical correlation of Powers et al. [Water Resour. Res. 30(2) (1994b) 321]. Overall, a methodology was developed for incorporating pore-scale processes into field-scale prediction of NAPL dissolution.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 2. Numerical simulation

A numerical model of surfactant enhanced solubilization was developed and applied to the simulation of nonaqueous phase liquid recovery in two-dimensional heterogeneous laboratory sand tank systems. Model parameters were derived from independent, small-scale, batch and column experiments. These parameters included viscosity, density, solubilization capacity, surfactant sorption, interfacial tension, permeability, capillary retention functions, and interphase mass transfer correlations. Model predictive capability was assessed for the evaluation of the micellar solubilization of tetrachloroethylene (PCE) in the two-dimensional systems. Predicted effluent concentrations and mass recovery agreed reasonably well with measured values. Accurate prediction of enhanced solubilization behavior in the sand tanks was found to require the incorporation of pore-scale, system-dependent, interphase mass transfer limitations, including an explicit representation of specific interfacial contact area. Predicted effluent concentrations and mass recovery were also found to depend strongly upon the initial NAPL entrapment configuration. Numerical results collectively indicate that enhanced solubilization processes in heterogeneous, laboratory sand tank systems can be successfully simulated using independently measured soil parameters and column-measured mass transfer coefficients, provided that permeability and NAPL distributions are accurately known. This implies that the accuracy of model predictions at the field scale will be constrained by our ability to quantify soil heterogeneity and NAPL distribution.     

The effect of multicomponent diffusion on NAPL dissolution from spherical ternary mixtures

This paper investigates the dissolution characteristics of ternary nonaqueous phase liquid (NAPL) mixtures with the goal of comparing the relative contributions of multicomponent (intra-NAPL) diffusion, film transfer and thermodynamic nonideality. These contributions are compared at the pore scale and intermediate scale (several centimeters downstream from the source). Trichloroethene (TCE), tetrachloroethene (PCE) and 1,1,1-trichloroethane (TCA) were selected to model a reasonably ideal mixture; TCE, PCE and octanol were selected as a relevant nonideal mixture. A multicomponent diffusion-based dissolution model incorporating hydrodynamic theory was formulated to estimate intra-NAPL concentration gradients and associated aqueous interfacial concentrations for ideally shaped (spherical) NAPL blobs. Pore scale dissolution times for this model were compared to those generated using the conventional well-mixed NAPL dissolution model, applying the same film transfer boundary condition in both cases. Activity coefficients (spatially and temporally variable for the diffusion model, temporally variable for the well-mixed model) were estimated using UNIFAC. NAPL interfacial concentration histories generated using the pore scale models were used as input in a three-dimensional groundwater transport model (MT3DMS) to compare downstream concentration distributions. For the relatively large NAPL bodies simulated (r=0.6 cm), intra-NAPL diffusion effects were found to be significant at the pore scale and were strongly impacted by the mixture's thermodynamic ideality. At the intermediate scale, and for the conditions tested, modest differences in the simulations suggested that intra-NAPL diffusion effects would be negligible compared to those associated with mixture composition uncertainty, dissolution rate processes related to NAPL-induced permeability effects and hydrodynamic issues associated with flow field heterogeneity.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

Investigation of surfactant-enhanced dissolution of entrapped nonaqueous phase liquid chemicals in a two-dimensional groundwater flow field.

Because of their low solubility, waste chemicals in the form of nonaqueous phase liquids (NAPLs) that are entrapped in subsurface formations act as long-term sources of groundwater contamination. In the design of remediation schemes that use surfactants, it is necessary to estimate the mass transfer rate coefficients under multi-dimensional flow fields that exit at field sites. In this study, we investigate mass transfer under a two-dimensional flow field to obtain an understanding of the basic mechanisms of surfactant-enhanced dissolution and to quantify the mass transfer rates. Enhanced dissolution experiments in a two-dimensional test cell were conducted to measure rates of mass depletion from entrapped NAPLs to a flowing aqueous phase containing a surfactant. In situ measurement of transient saturation changes using a gamma attenuation system revealed dissolution patterns that are affected by the dimensionality of the groundwater flow field. Numerical modeling of local flow fields that changed with time, due to depletion of NAPL sources, enabled the examination of the basic mechanisms of NAPL dissolution in complex groundwater systems. Through nonlinear regression analysis, mass transfer rates were correlated to porous media properties, NAPL saturation and aqueous phase velocity. Results from the experiments and numerical analyses were used to identify deficiencies in existing methods of analysis that uses assumptions of one-dimensional flow, homogeneity of aquifer properties, local equilibrium and idealized transient mass transfer.     

Mass flux from a non-aqueous phase liquid pool considering spontaneous expansion of a discontinuous gas phase

The partitioning of non-aqueous phase liquid (NAPL) compounds to a discontinuous gas phase results in the repeated spontaneous expansion, snap-off, and vertical mobilization of the gas phase. This mechanism has the potential to significantly affect the mass transfer processes that control the dissolution of NAPL pools by increasing the vertical transport of NAPL mass and increasing the total mass transfer rate from the surface of the pool. The extent to which this mechanism affects mass transfer from a NAPL pool depends on the rate of expansion and the mass of NAPL compound in the gas phase. This study used well-controlled bench-scale experiments under no-flow conditions to quantify for the first time the expansion of a discontinuous gas phase in the presence of NAPL. Air bubbles placed in glass vials containing NAPL increased significantly in volume, from a radius of 1.0 mm to 2.0 mm over 215 days in the presence of tetrachloroethene (PCE), and from a radius of 1.2 mm to 2.3 mm over 22 days in the presence of trans-1,2-dichloroethene (tDCE). A one-dimensional mass transfer model, fit to the experimental data, showed that this expansion could result in a mass flux from the NAPL pool that was similar in magnitude to the mass flux expected for the dissolution of a NAPL pool in a two-fluid (NAPL and water) system. Conditions favouring the significant effect of a discontinuous gas phase on mass transfer were identified as groundwater velocities less than ~0.01 m/day, and a gas phase that covers greater than ~10% of the pool surface area and is located within ~0.01 m of the pool surface. Under these conditions the mass transfer via a discontinuous gas phase is expected to affect, for example, efforts to locate NAPL source zones using aqueous concentration data, and predict the lifetime and risk associated with NAPL source zones in a way that is not currently included in the common conceptual models used to assess NAPL-contaminated sites.     

Locally-calibrated light transmission visualization methods to quantify nonaqueous phase liquid mass in porous media

Five locally-calibrated light transmission visualization (LTV) methods were tested to quantify nonaqueous phase liquid (NAPL) mass and mass reduction in porous media. Tetrachloroethylene (PCE) was released into a two-dimensional laboratory flow chamber packed with water-saturated sand which was then flushed with a surfactant solution (2% Tween 80) until all of the PCE had been dissolved. In all the LTV methods employed here, the water phase was dyed, rather than the more common approach of dyeing the NAPL phase, such that the light adsorption characteristics of NAPL did not change as dissolution progressed. Also, none of the methods used here required the use of external calibration chambers. The five visualization approaches evaluated included three methods developed from previously published models, a binary method, and a novel multiple wavelength method that has the advantage of not requiring any assumptions about the intra-pore interface structure between the various phases (sand/water/NAPL). The new multiple wavelength method is also expected to be applicable to any translucent porous media containing two immiscible fluids (e.g., water-air, water-NAPL). Results from the sand-water-PCE system evaluated here showed that the model that assumes wetting media of uniform pore size (Model C of Niemet and Selker, 2001) and the multiple wavelength model with no interface structure assumptions were able to accurately quantify PCE mass reduction during surfactant flushing. The average mass recoveries from these two imaging methods were greater than 95% for domain-average NAPL saturations of approximately 2.6x10(-2), and were approximately 90% during seven cycles of surfactant flushing that sequentially reduced the average NAPL saturation to 7.5x10(-4).     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland-Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.     

Dissolution kinetics of high explosives particles in a saturated sandy soil

Solid phase high explosive (HE) residues from munitions detonation may be a persistent source of soil and groundwater contamination at military training ranges. Saturated soil column tests were conducted to observe the dissolution behavior of individual components (RDX, HMX, and TNT) from two HE formulations (Comp B and C4). HE particles dissolved readily, with higher velocities yielding higher dissolution rates, higher mass transfer coefficients, and lower effluent concentrations. Effluent concentrations were below solubility limits for all components at superficial velocities of 10-50 cm day(-1). Under continuous flow at 50 cm day(-1), RDX dissolution rates from Comp B and C4 were 34.6 and 97.6 microg h(-1) cm(-2) (based on initial RDX surface area), respectively, significantly lower than previously reported dissolution rates. Cycling between flow and no-flow conditions had a small effect on the dissolution rates and effluent concentrations; however, TNT dissolution from Comp B was enhanced under intermittent-flow conditions. A model that includes advection, dispersion, and film transfer resistance was developed to estimate the steady-state effluent concentrations.     

The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation

Understanding the process of mass transfer from source zones of aquifers contaminated with organic chemicals in the form of dense non-aqueous phase liquids (DNAPL) is of importance in site management and remediation. A series of intermediate-scale tank experiments was conducted to examine the influence of aquifer heterogeneity on DNAPL mass transfer contributing to dissolved mass emission from source zone into groundwater under natural flow before and after remediation. A Tetrachloroethylene (PCE) spill was performed into six source zone models of increasing heterogeneity, and both the spatial distribution of the dissolution behavior and the net effluent mass flux were examined. Experimentally created initial PCE entrapment architecture resulting from the PCE migration was largely influenced by the coarser sand lenses and the PCE occupied between 30 and 60% of the model aquifer depth. The presence of DNAPL had no apparent effect on the bulk hydraulic conductivity of the porous media. Up to 71% of PCE mass in each of the tested source zone was removed during a series of surfactant flushes, with associated induced PCE mobilization responsible for increasing vertical DNAPL distributions. Effluent mass flux due to water dissolution was also found to increase progressively due to the increase in NAPL-water contact area even though the PCE mass was reduced. Doubling of local groundwater flow velocities showed negligible rate-limited effects at the scale of these experiments. Thus, mass transfer behavior was directly controlled by the morphology of DNAPL within each source zone. Effluent mass flux values were normalized by the up-gradient DNAPL distributions. For the suite of aquifer heterogeneities and all remedial stages, normalized flux values fell within a narrow band with mean of 0.39 and showed insensitivity to average source zone saturations.     

Simple screening models of NAPL dissolution in the subsurface

Three simple screening models of nonaqueous phase liquid (NAPL) dissolution in the subsurface are proposed based on the NAPL mass conservation and the assumption of proportionality between the residual NAPL source zone concentration and the remaining residual NAPL mass. The purpose of the proposed models is to predict the solute concentration in the zone of the residual NAPL as a result of dissolution. The predicted source zone concentration decrease is used to simulate and account for the decrease of dissolution rate with time. The proposed simple NAPL dissolution models enable the pseudo-equilibrium formulation to be used and therefore the numerical simulations for field application problems can be simplified compared to the non-equilibrium counterpart. With proper choice of empirical parameters, the proposed simple screening models can work as well as more complex dissolution rate correlation models, such as that of Imhoff et al. [Water Resour. Res. 30 (1994) 307-320]. It is found that the proposed models are very good for quantifying non-equilibrium dissolution, which is characterized by tailing of breakthrough curves. The models are especially useful for situations of small residual NAPL saturation, which are typical for many field applications.     

Hot water flushing for immiscible displacement of a viscous NAPL

Thermal remediation techniques, such as hot water flooding, are emerging technologies that have been proposed for the removal of nonaqueous phase liquids (NAPLs) from the subsurface. In this study a combined laboratory and modeling investigation was conducted to determine if hot water flooding techniques would improve NAPL mass removal compared to ambient temperature water flushing. Two experiments were conducted in a bench scale two-dimensional sandbox (55 cmx45 cmx1.3 cm) and NAPL saturations were quantified using a light transmission apparatus. In these immiscible displacement experiments the aqueous phase, at 22 degrees C and 50 degrees C, displaced a zone with initial NAPL saturations on the order of 85%. The interfacial tension and viscosity of the selected light NAPL, Voltesso 35, are strongly temperature-dependent. Experimental results suggest that hot water flooding reduced the size of the high NAPL saturation zone, in comparison to the cold water flood, and yielded greater NAPL mass recovery (75% NAPL removal vs. 64%). Hot water flooding did not, however, result in lower residual NAPL saturations. A numerical simulator was modified to include simultaneous flow of water and organic phases, energy transport, temperature and pressure. Model predictions of mass removal and NAPL saturation profiles compared well with observed behavior. A sensitivity analysis indicates that the utility of hot water flooding improves with the increasing temperature dependence of NAPL hydraulic properties.     

Parameters that control the cleanup of fractured permeable aquifers

This study develops a modeling approach for simulating and evaluating entrapped light nonaqueous-phase liquid (light NAPL-LNAPL) dissolution and transport of the solute in a fractured permeable aquifer (FPA). The term FPA refers to an aquifer made of porous blocks of high permeability that embed fractures. The fracture network is part of the domain characterized by high permeability and negligible storage. Previous studies show that sandstone aquifers often represent FPAs. The basic model developed in this study is a two-dimensional (2-D) model of permeable blocks that embed oblique equidistant fractures with constant aperture and orientation. According to this model, two major parameters govern NAPL dissolution and transport of the solute. These parameters are: 1) the dimensionless interphase mass transfer coefficient, K(f0), and 2) the mobility number, N(M0). These parameters represent measures of heterogeneity affecting flow, NAPL dissolution, and transport of the solute in the domain. The parameter K(f0) refers to the rate at which organic mass is transferred from the NAPL into the water phase. The parameter N(M0) represents the ratio of flow through the porous blocks to flow through the fracture network in regions free of entrapped NAPL. It also provides a measure of groundwater flow bypassing regions contaminated by entrapped NAPL. In regions contaminated by entrapped NAPL our simulations have often indicated very low permeability of the porous blocks, enabling a significant increase of the fracture flow at the expense of the permeable block flow. Two types of constitutive relationships also affect the rate of FPA cleanup: 1) the relationship between the saturation of the entrapped NAPL and the permeability of the porous blocks, and 2) the relationships representing effects of the entrapped NAPL saturation and the permeable block flow velocity on rates of interphase mass transfer. This study provides basic tools for evaluating the characteristics of pump-and-treat cleanup of FPAs by referring to sets of parameters and constitutive relationships typical of FPAs. The numerical simulations carried out in this study show that at high initial saturation of the entrapped NAPL, during initial stages of the FPA cleanup the contaminant concentration increases, but later it decreases. This phenomenon originates from significant groundwater bypassing the NAPL entrapped in the permeable blocks via the fracture network.