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.     

Diffusive partitioning tracer test for the quantification of nonaqueous phase liquid (NAPL) in the vadose zone: Performance evaluation for heterogeneous NAPL distribution

A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers.     

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.     

Radon as a naturally occurring tracer for the assessment of residual NAPL contamination of aquifers

The noble gas radon has a strong affinity to non-aqueous phase-liquids (NAPLs). That property makes it applicable as naturally occurring partitioning tracer for assessing residual NAPL contamination of aquifers. In a NAPL contaminated aquifer, radon dissolved in the groundwater partitions preferably into the NAPL. The magnitude of the resulting radon deficit in the groundwater depends on the NAPL-specific radon partition coefficient and on the NAPL saturation of the pore space. Hence, if the partition coefficient is known, the NAPL saturation is attainable by determination of the radon deficit. After a concise discussion of theoretical aspects regarding radon partitioning into NAPL, related experimental data and results of a field investigation are presented. Aim of the laboratory experiments was the determination of radon partition coefficients of multi-component NAPLs of environmental concern. The on-site activities were carried out in order to confirm the applicability of the "radon method" under field conditions.     

Tracer test for the measurement of gas diffusion and non-aqueous phase liquid (NAPL) saturation in soil

During soil bioremediation, the diffusion of oxygen into the soil is an important prerequisite for aerobic biodegradation, and the decrease of petroleum products is the ultimate goal. Both processes need to be monitored. The aim of this work was to develop a gas tracer test that yields information on both, gas diffusion and residual saturation with non-aqueous phase liquids (NAPLs) in unsaturated soil heaps. One conservative tracer (methane) and 4 partitioning gas tracers (diethylether, methyl tert-butyl ether, chloroform and n-heptane) were injected as vapors into laboratory columns filled with unsaturated sand with increasing NAPL saturation. Breakthrough curves of gaseous compounds were measured at two points and compared to analytical solutions of an analytical diffusive-reactive transport equation. By fitting of methane data, robust results for effective diffusivity (tortuosity) were obtained. NAPL saturation was most accurately measured by the moderately water soluble tracers (ethers and chloroform). The hydrophobic tracer n-heptane did not partition into water-immersed NAPL. An easy and accurate way to assess air-NAPL partitioning constants from gas chromatography retention times is furthermore reported. It is concluded that gas tracer tests have the potential for measuring two important properties in soil bioremediation systems easily and quickly.     

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.     

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.     

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.     

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.     

Performance assessment of NAPL remediation in heterogeneous alluvium

Over the last few years, more than 40 partitioning interwell tracer tests (PITTs) have been conducted at many different sites to measure nonaqueous phase liquid (NAPL) saturations in the subsurface. While the main goal of these PITTs was to estimate the NAPL volume in the subsurface, some were specifically conducted to assess the performance of remedial actions involving NAPL removal. In this paper, we present a quantitative approach to assess the performance of remedial actions to recover NAPL that can be used to assess any NAPL removal technology. It combines the use of PITTs (to estimate the NAPL volume in the swept pore volume between injection and extraction wells of a test area) with the use of several cores to determine the vertical NAPL distribution in the subsurface. We illustrate the effectiveness of such an approach by assessing the performance of a surfactant/foam flood conducted at Hill Air Force Base, UT, to remove a TCE-rich NAPL from alluvium with permeability contrasts as high as one order of magnitude. In addition, we compare the NAPL volumes determined by the PITTs with volumes estimated through geostatistical interpolation of aquifer sediment core data collected with a vertical frequency of 5-10 cm and a lateral borehole spacing of 0.15 m. We demonstrate the use of several innovations including the explicit estimation of not only the errors associated with NAPL volumes and saturations derived from PITTs but also the heterogeneity of the aquifer sediments based upon permeability estimates. Most importantly, we demonstrate the reliability of the     

Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.     

A set of constitutive relationships accounting for residual NAPL in the unsaturated zone

Although laboratory experiments show that non-aqueous phase liquid (NAPL) is retained in the unsaturated zone, no existing multiphase flow model has been developed to account for residual NAPL after NAPL drainage in the unsaturated zone. We developed a static constitutive set of saturation-capillary pressure relationships for water, NAPL and air that accounts for both this residual NAPL and entrapped NAPL. The set of constitutive relationships is formulated similarly to the set of scaled relationships that is frequently applied in continuum models. The new set consists of three fluid-phase systems: a three-phase system and a two-phase system, that both comply with the original constitutive model, and a newly introduced residual NAPL system. The new system can be added relatively easily to the original two- and three-phase systems. Entrapment is included in the model. The constitutive relationships of the non-drainable residual NAPL system are based on qualitative fluid behavior derived from a pore scale model. The pore scale model reveals that the amount of residual NAPL depends on the spreading coefficient and the water saturation. Furthermore, residual NAPL is history-dependent. At the continuum scale, a critical NAPL pressure head defines the transition from free, mobile NAPL to residual NAPL. Although the Pc-S relationships for water and total liquid are not independent in case of residual NAPL, two two-phase Pc-S relations can represent a three-phase residual system of Pc-S relations. A newly introduced parameter, referred to as the residual oil pressure head, reflects the mutual dependency of water and oil. Example calculations show consistent behavior of the constitutive model. Entrapment and retention in the unsaturated zone cooperate to retain NAPL. Moreover, the results of our constitutive model are in agreement with experimental observations.     

Design of aquifer remediation systems: (1) Describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.     

Biosurfactant-enhanced solubilization of NAPL mixtures

Remediation of nonaqueous phase liquids (NAPLs) by conventional pump-and-treat methods (i.e., water flushing) is generally considered to be ineffective due to low water solubilities of NAPLs and to mass-transfer constraints. Chemical flushing techniques, such as surfactant flushing, can greatly improve NAPL remediation primarily by increasing the apparent solubility of NAPL contaminants. NAPLs at hazardous waste sites are often complex mixtures. However, the equilibrium and nonequilibrium mass-transfer characteristics between NAPL mixtures and aqueous surfactant solutions are not well understood. This research investigates the equilibrium solubilization behavior of two- and three-component NAPL mixtures (containing akylbenzenes) in biosurfactant solutions. NAPL solubilization is found to be ideal in water (i.e., obeys Raoult's Law), while solubilization in biosurfactant solutions was observed to be nonideal. Specifically, the relatively hydrophobic compounds in the mixture experienced solubility enhancements that were greater than those predicted by ideal enhanced solubilization theory, while the solubility enhancements for the relatively hydrophilic compounds were less than predicted. The degree of nonideality is shown to be a nonlinear function of the NAPL-phase mole fraction. Empirical relationships based on the NAPL-phase mole fraction and/or micelle-aqueous partition coefficients measured in single-component NAPL systems are developed to estimate values for the multicomponent partition coefficients. Empirical relationships that incorporate both the NAPL-phase mole fraction and single-component partition coefficients yield much improved estimates for the multicomponent partition coefficient.     

Magnetic resonance imaging of nonaqueous phase liquid during soil vapor extraction in heterogeneous porous media

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.     

Numerical simulations of radon as an in situ partitioning tracer for quantifying NAPL contamination using push-pull tests

Presented here is a reanalysis of results previously presented by [Davis, B.M., Istok, J.D., Semprini, L., 2002. Push-pull partitioning tracer tests using radon-222 to quantify non-aqueous phase liquid contamination. J. Contam. Hydrol. 58, 129-146] of push-pull tests using radon as a naturally occurring partitioning tracer for evaluating NAPL contamination. In a push-pull test where radon-free water and bromide are injected, the presence of NAPL is manifested in greater dispersion of the radon breakthrough curve (BTC) relative to the bromide BTC during the extraction phase as a result of radon partitioning into the NAPL. Laboratory push-pull tests in a dense or DNAPL-contaminated physical aquifer model (PAM) indicated that the previously used modeling approach resulted in an overestimation of the DNAPL (trichloroethene) saturation (S(n)). The numerical simulations presented here investigated the influence of (1) initial radon concentrations, which vary as a function of S(n), and (2) heterogeneity in S(n) distribution within the radius of influence of the push-pull test. The simulations showed that these factors influence radon BTCs and resulting estimates of S(n). A revised method of interpreting radon BTCs is presented here, which takes into account initial radon concentrations and uses non-normalized radon BTCs. This revised method produces greater radon BTC sensitivity at small values of S(n) and was used to re-analyze the results from the PAM push-pull tests reported by Davis et al. The re-analysis resulted in a more accurate estimate of S(n) (1.8%) compared with the previously estimated value (7.4%). The revised method was then applied to results from a push-pull test conducted in a light or LNAPL-contaminated aquifer at a field site, resulting in a more accurate estimate of S(n) (4.1%) compared with a previously estimated value (13.6%). The revised method improves upon the efficacy of the radon push-pull test to estimate NAPL saturations. A limitation of the revised method is that 'background' radon concentrations from a non-contaminated well in the NAPL-contaminated aquifer are needed to accurately estimate NAPL saturation. The method has potential as a means of monitoring the progress of NAPL remediation.     

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.     

A constitutive model for air-NAPL-water flow in the vadose zone accounting for immobile, non-occluded (residual) NAPL in strongly water-wet porous media

A hysteretic constitutive model describing relations among relative permeabilities, saturations, and pressures in fluid systems consisting of air, nonaqueous-phase liquid (NAPL), and water is modified to account for NAPL that is postulated to be immobile in small pores and pore wedges and as films or lenses on water surfaces. A direct outcome of the model is prediction of the NAPL saturation that remains in the vadose zone after long drainage periods (residual NAPL). Using the modified model, water and NAPL (free, entrapped by water, and residual) saturations can be predicted from the capillary pressures and the water and total-liquid saturation-path histories. Relations between relative permeabilities and saturations are modified to account for the residual NAPL by adjusting the limits of integration in the integral expression used for predicting the NAPL relative permeability. When all of the NAPL is either residual or entrapped (i.e., no free NAPL), then the NAPL relative permeability will be zero. We model residual NAPL using concepts similar to those used to model residual water. As an initial test of the constitutive model, we compare predictions to published measurements of residual NAPL. Furthermore, we present results using the modified constitutive theory for a scenario involving NAPL imbibition and drainage.     

Bioenhancement of NAPL pool dissolution: experimental evaluation

Experiments were conducted to quantify nonaqueous phase liquid (NAPL) pool dissolution and its enhancement by in situ biodegradation. The experiments were performed using square cross-section, glass-bead packed column reactors with a small pool of a toluene-in-dodecane mixture (toluene mole fraction, X(tol) approximately 0.02 or 0.09). Experimental quasi-steady-state toluene dissolution fluxes were determined using a 14C-carbon mass-balance approach during water flushing with and without biodegradation. The experiments demonstrated a statistically significant bioenhancement of the toluene dissolution flux of up to roughly twofold at average pore water velocities of approximately 0.1 and 1 m/day when the toluene mole fraction was low ( approximately 0.02); however, little or no bioenhancement was observed with the higher mole fraction ( approximately 0.09). Although it cannot be determined conclusively, the weight of evidence based on biomass measurements and model analyses suggests that the reduced bioenhancement for the high mole fraction was due to higher dissolved toluene concentrations, which may have caused toxicity effects. Importantly, even though NAPL dissolution was not bioenhanced in every case, the biodegradation reduced toluene concentrations to low levels in the reactor effluents.     

A proposed model to include a residual NAPL saturation in a hysteretic capillary pressure-saturation relationship

A residual non-aqueous phase liquid (NAPL) present in the vadose zone can act as a contaminant source for many years as the compounds of concern partition to infiltrating groundwater and air contained in the soil voids. Current pressure-saturation-relative permeability relationships do not include a residual NAPL saturation term in their formulation. This paper presents the results of series of two- and three-phase pressure cell experiments conducted to evaluate the residual NAPL saturation and its impact on the pressure-saturation relationship. A model was proposed to incorporate a residual NAPL saturation term into an existing hysteretic three-phase parametric model developed by Parker and Lenhard [Water Resour. Res. 23(12) (1987) 2187], Lenhard and Parker [Water Resour. Res. 23(12) (1987) 2197] and Lenhard [J. Contam. Hydrol. 9 (1992) 243]. The experimental results indicated that the magnitude of the residual NAPL saturation was a function of the maximum total liquid saturation reached and the water saturation. The proposed model to incorporate a residual NAPL saturation term is similar in form to the entrapment model proposed by Parker and Lenhard, which was based on an expression presented by Land [Soc. Pet. Eng. J. (June 1968) 149].