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.     

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.     

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.     

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.     

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.     

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.     

Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.     

Refinement of the density-modified displacement method for efficient treatment of tetrachloroethene source zones

A novel method to remediate dense nonaqueous phase liquid (DNAPL) source zones that incorporates in situ density conversion of DNAPL via alcohol partitioning followed by displacement with a low interfacial tension (IFT) surfactant flood has been developed. Previous studies demonstrated the ability of the density-modified displacement (DMD) method to recover chlorobenzene (CB) and trichloroethene (TCE) from heterogeneous porous media without downward migration of the dissolved plume or free product. However, the extent of alcohol (n-butanol) partitioning required for in situ density conversion of high-density NAPLs, such as tetrachloroethene (PCE), could limit the utility of the DMD method. Hence, the objective of this study was to compare the efficacy of two n-butanol delivery approaches: an aqueous solution of 6% (wt) n-butanol and a surfactant-stabilized macroemulsion containing 15% (vol) n-butanol in water, to achieve density reduction of PCE-NAPL in two-dimensional (2-D) aquifer cells. Results of liquid-liquid equilibrium studies indicated that density conversion of PCE relative to water occurred at an n-butanol mole fraction of 0.56, equivalent to approximately 5 ml n-butanol per 1 ml of PCE when in equilibrium with an aqueous solution. In 2-D aquifer cell studies, density conversion of PCE was realized using both n-butanol preflood solutions, with effluent NAPL samples exhibiting density reductions ranging from 0.51 to 0.70 g/ml. Although the overall PCE mass recoveries were similar (91% and 93%) regardless of the n-butanol delivery method, the surfactant-stabilized macroemulsion preflood removed approximately 50% of the PCE mass. In addition, only 1.2 pore volumes of the macroemulsion solution were required to achieve in situ density conversion of PCE, compared to 6.4 pore volumes of the 6% (wt) n-butanol solution. These findings demonstrate that use of the DMD method with a surfactant-stabilized macroemulsion containing n-butanol holds promise as an effective source zone remediation technology, allowing for efficient recovery of PCE-DNAPL while mitigating downward migration of the dissolved plume and free product.     

Field demonstration of surfactant-enhanced solubilization of DNAPL at Dover Air Force Base, Delaware

This study reports on a surfactant-based flood for tetrachloroethylene (PCE) removal from a control test cell at the Dover National Test Site. The surfactant formulation (sodium dihexyl sulfosuccinate (Aerosol-MA or AMA), isopropanol and calcium chloride) was able to achieve a high concentration of PCE in swollen micelles (supersolubilization) without vertical PCE migration. The hydraulic system included eight screened wells that were operated in both vertical circulation and line drive configurations. After 10 pore volumes of flushing, the overall PCE removal was 68% (65% of which corresponded to the surfactant flooding alone). In addition, the residual PCE saturation was reduced from 0.7% to 0.2%, and the concentration of PCE in the groundwater was reduced from 37-190 mg/L before the flushing to 7.3 mg/L after flooding. Recycling the surfactant solution reduced the required surfactant mass (and thus cost, and waste) by 90%. Close to 80% of the total PCE removal was obtained during the first five pore volumes which were operated in an upward vertical circulation flow scheme. No free oil phase was observed during the test. Further analysis of multilevel sampler data suggests that most of the trapped oil remaining in the cell was likely localized in secluded regions of the aquifer, which helps explain the lower PCE groundwater concentration after remedial activities. In summary, this field study demonstrated the feasibility of surfactant-enhanced remediation to reduce the mass in the source zone and significantly reduce the PCE aqueous concentration and therefore the risk associated with the contaminant plume.     

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.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 1. Experimental studies

A matrix of batch, column and two-dimensional (2-D) box experiments was conducted to investigate the coupled effects of rate-limited solubilization and layering on the entrapment and subsequent recovery of a representative dense NAPL, tetrachloroethylene (PCE), during surfactant flushing. Batch experiments were performed to determine the equilibrium solubilization capacity of the surfactant, polyoxyethylene (20) sorbitan monooleate (Tween 80), and to measure fluid viscosity, density and interfacial tension. Results of one-dimensional column studies indicated that micellar solubilization of residual PCE was rate-limited at Darcy velocities ranging from 0.8 to 8.2 cm/h and during periods of flow interruption. Effluent concentration data were used to develop effective mass transfer coefficient (Ke) expressions that were dependent upon the Darcy velocity and duration of flow interruption. To simulate subsurface heterogeneity, 2-D boxes were packed with layers of F-70 Ottawa sand and Wurtsmith aquifer material within 20-30 mesh Ottawa sand. A 4% Tween 80 solution was then flushed through PCE-contaminated boxes at several flow velocities, with periods of flow interruption. Effluent concentration data and visual observations indicated that both rate-limited solubilization and pooling of PCE above the fine layers reduced PCE recovery to levels below those anticipated from batch and column measurements. These experimental results demonstrate the potential impact of both mass transfer limitations and subsurface layering on the recovery of PCE during surfactant enhanced aquifer remediation.     

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.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

Wettability hysteresis and its implications for DNAPL source zone distribution

Subsurface heterogeneity at sites contaminated with nonaqueous phase liquids (NAPLs) reduces the effectiveness of traditional remediation measures. One cause may be the increased proportion of NAPL that is hydraulically isolated due to capillary trapping in heterogeneously-wetted materials. This study examines the wettability of ten materials, ranging from minerals, such as calcite and dolomite, to carbonaceous materials, such shale and coal, in air and water, NAPL and air, and NAPL and water systems. The wettability differed depending on which phase the solid material was initially immersed in: the less crystalline solids, if initially contacted by water were water-wet, but if initially contacted by NAPL were NAPL-wet. This difference, termed here wettability hysteresis, was observed for a suite of halogenated NAPLs and was independent of equilibration time. The degree of wettability hysteresis was greatest in the NAPL and water systems, with the magnitude of the difference increasing with the carbonaceous materials. Since the degree of capillary trapping in subsurface materials is related to wettability, the phenomenon of wettability hysteresis suggests that system history is a factor that may increase the heterogeneity of NAPL source zones.     

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.     

Investigating the role of gas bubble formation and entrapment in contaminated aquifers: Reactive transport modelling

In many natural and contaminated aquifers, geochemical processes result in the production or consumption of dissolved gases. In cases where methanogenesis or denitrification occurs, the production of gases may result in the formation and growth of gas bubbles below the water table. Near the water table, entrapment of atmospheric gases during water table rise may provide a significant source of O(2) to waters otherwise depleted in O(2). Furthermore, the presence of bubbles will affect the hydraulic conductivity of an aquifer, resulting in changes to the groundwater flow regime. The interactions between physical transport, biogeochemical processes, and gas bubble formation, entrapment and release is complex and requires suitable analysis tools. The objective of the present work is the development of a numerical model capable of quantitatively assessing these processes. The multicomponent reactive transport code MIN3P has been enhanced to simulate bubble growth and contraction due to in-situ gas production or consumption, bubble entrapment due to water table rise and subsequent re-equilibration of the bubble with ambient groundwater, and permeability changes due to trapped gas phase saturation. The resulting formulation allows for the investigation of complex geochemical systems where microbially mediated redox reactions both produce and consume gases as well as affect solution chemistry, alkalinity, and pH. The enhanced model has been used to simulate processes in a petroleum hydrocarbon contaminated aquifer where methanogenesis is an important redox process. The simulations are constrained by data from a crude oil spill site near Bemidji, MN. Our results suggest that permeability reduction in the methanogenic zone due to in-situ formation of gas bubbles, and dissolution of entrapped atmospheric bubbles near the water table, both work to attenuate the dissolved gas plume emanating from the source zone. Furthermore, the simulations demonstrate that under the given conditions more than 50% of all produced CH(4) partitions to the gas phase or is aerobically oxidised near the water table, suggesting that these processes should be accounted for when assessing the rate and extent of methanogenic degradation of hydrocarbons.     

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.     

Laboratory investigation of flux reduction from dense non-aqueous phase liquid (DNAPL) partial source zone remediation by enhanced dissolution

This study investigated the benefits of partial removal of dense nonaqueous phase liquid (DNAPL) source zones using enhanced dissolution in eight laboratory scale experiments. The benefits were assessed by characterizing the relationship between reductions in DNAPL mass and the corresponding reduction in contaminant mass flux. Four flushing agents were evaluated in eight controlled laboratory experiments to examine the effects of displacement fluid property contrasts and associated override and underride on contaminant flux reduction (R(j)) vs. mass reduction (R(m)) relationships (R(j)(R(m))): 1) 50% ethanol/50% water (less dense than water), 2) 40% ethyl-lactate/60% water (more dense than water), 3) 18% ethanol/26% ethyl-lactate/56% water (neutrally buoyant), and 4) 2% Tween-80 surfactant (also neutrally buoyant). For each DNAPL architecture evaluated, replicate experiments were conducted where source zone dissolution was conducted with a single flushing event to remove most of the DNAPL from the system, and with multiple shorter-duration floods to determine the path of the R(j)(R(m)) relationship. All of the single-flushing experiments exhibited similar R(j)(R(m)) relationships indicating that override and underride effects associated with cosolvents did not significantly affect the remediation performance of the agents. The R(j)(R(m)) relationship of the multiple injection experiments for the cosolvents with a density contrast with water tended to be less desirable in the sense that there was less R(j) for a given R(m). UTCHEM simulations supported the observations from the laboratory experiments and demonstrated the capability of this model to predict R(j)(R(m)) relationships for non-uniformly distributed NAPL sources.     

Simultaneous optimization of dense non-aqueous phase liquid (DNAPL) source and contaminant plume remediation

A framework is developed for simultaneous, optimal design of groundwater contaminant source removal and plume remediation strategies. The framework allows for varying degrees of effort and cost to be dedicated to source removal versus plume remediation. We have accounted for the presence of physical heterogeneity in the DNAPL source, since source heterogeneity controls mass release into the plume and the efficiency of source removal efforts. We considered high and low estimates of capital and operating costs for chemical flushing removal of the source, since these are expected to vary form site to site. Using the lower chemical flushing cost estimates, it is found that the optimal allocation of funds to source removal or plume remediation is sensitive to the degree of heterogeneity in the source. When the time elapsed between the source release and the implementation of remediation was varied, it was found that, except for the longest elapsed time (50,000 days), a combination of partial source removal and plume remediation was most efficient. When first-order, dissolved contaminant degradation was allowed, source removal was found to be unnecessary for the cases where the degradation rate exceeded intermediate values of the first-order rate constant. Finally, it was found that source removal became more necessary as the degree of aquifer heterogeneity increased.     

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.