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.     

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.     

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.     

The dissolution of benzene, toluene, m-xylene and naphthalene from a residually trapped non-aqueous phase liquid under mass transfer limited conditions

The results of dissolution experiments for benzene, toluene, m-xylene and naphthalene (BTXN) from a relatively insoluble oil phase (tridecane), residually trapped in a non-sorbing porous medium, are described. This mixture was chosen to simulate dissolution of soluble aromatic compounds from a petroleum hydrocarbon mixture, e.g., crude oil, for which a large fraction of the mixture is relatively insoluble. The experiments were carried out at a small source length to interstitial velocity ratio, L/v, so that dissolution would be mass transfer limited (MTL). When fitted to data for toluene, a multiregion mass transfer model was found to predict the experimental data satisfactorily for the other components without adjustment of the mass transfer rate parameters. These results indicate that the dissolution process can be generalized for various hydrophobic organic compounds present in a multicomponent non-aqueous phase liquid (NAPL) when mass transfer limitations are present. This also suggests that dissolution data obtained for one compound can be useful for predicting the dissolution histories for other compounds from petroleum hydrocarbon mixtures.     

Pulsed gas injection: A minimum effort approach for enhanced natural attenuation of chlorobenzene in contaminated groundwater

Chlorobenzene-contaminated groundwater was used to assess pulsed gas sparging as a minimum effort aeration strategy to enhance intrinsic natural attenuation. In contrast to existing biosparging operations, oxygen was supplied at minimum rate by reducing the gas injection frequency to 0.33 day^?1. Field tests in a model aquifer were conducted in a 12 m long reactor, filled with indigenous aquifer material and continuously recharged with polluted groundwater over 3 years. The closed arrangement allowed yield balances, cost accounting as well as the investigation of spatial distributions of parameters which are sensitive to the biodegradation process. Depending on the injection frequency and on the gas chosen for injection (pure oxygen or air) oxygen-deficient conditions prevailed in the aquifer. Despite the limiting availability of dissolved oxygen in the groundwater, chlorobenzene degradation under oxygen-deficient conditions proved to be more effective than under conditions with dissolved oxygen being available in high concentrations.     

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.     

Multiphase flow and transport caused by spontaneous gas phase growth in the presence of dense non-aqueous phase liquid

Disconnected bubbles or ganglia of trapped gas may occur below the top of the capillary fringe through a number of mechanisms. In the presence of dense non-aqueous phase liquid (DNAPL), the disconnected gas phase experiences mass transfer of dissolved gases, including volatile components from the DNAPL. The properties of the gas phase interface can also change. This work shows for the first time that when seed gas bubbles exist spontaneous gas phase growth can be expected to occur and can significantly affect water-gas-DNAPL distributions, fluid flow, and mass transfer. Source zone behaviour was observed in three different experiments performed in a 2-dimensional flow cell. In each case, a DNAPL pool was created in a zone of larger glass beads over smaller glass beads, which served as a capillary barrier. In one experiment effluent water samples were analyzed to determine the vertical concentration profile of the plume above the pool. The experiments effectively demonstrated a) a cycle of spontaneous gas phase expansion and vertical advective mobilization of gas bubbles and ganglia above the DNAPL source zone, b) DNAPL redistribution caused by gas phase growth and mobilization, and c) that these processes can significantly affect mass transport from a NAPL source zone.     

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.     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

Transport of a non-volatile contaminant in unsaturated porous media in the presence of colloids

Stable colloidal particles can travel long distances in subsurface environments and carry particle-reactive contaminants with them to locations further than predicted by the conventional advective-dispersive transport equation. When such carriers exist in a saturated porous medium, the system can be idealized as consisting of three phases: an aqueous phase, a carrier phase, and a stationary solid matrix phase. However, when colloids are present in an unsaturated porous medium, the system representation should include one more phase, i.e. the air phase. In the work reported, a mathematical model was developed to describe the transport and fate of the colloidal particles and a non-volatile contaminant in unsaturated porous media. The model is based on mass balance equations in a four-phase porous medium. Colloid mass transfer mechanisms among aqueous, solid matrix, and air phases, and contaminant mass transfer between aqueous and colloid phases are represented by kinetic expressions. Governing equations are non-dimensionalized and solved to investigate colloid and contaminant transport in an unsaturated porous medium. A sensitivity analysis of the transport model was utilized to assess the effects of several parameters on model behavior. The colloid transport model matches successfully with experimental data of Wan and Wilson. The presence of air-water interface retards the colloid transport significantly counterbalancing the facilitating effect of colloids. However, the retardation of contaminant transport by colloids is highly dependent on the properties of the contaminant and the colloidal surface.     

氮源对气相生物滴滤床净化效果的影响

氮源对于气相生物滴滤床的净化效果有显著的影响。由于喷淋液体在生物滴滤床中分布不均匀 ,导致床层中局部氮源不充足 ,从而使滤床的去除效果降低。提高喷淋液中氮源的浓度 ,可以较大幅度提高净化效率。采用NH-4 N会导致滤床中硝化细菌的生长和累积 ,使得异养菌和硝化细菌在溶解氧、氮源和填料空间上形成竞争 ,并致使氮源分布不均匀 ,而采用NO3 N作为氮源可以改善这种状况     

Landfill Gas Effects on Groundwater Samples at a Municipal Solid Waste Facility

A study was performed to determine the source of low concentrations of volatile organic compounds (VOCs) detected in groundwater samples at a solid waste management facility. The affected wells were identified as hydraulically upgradient of an old unlined facility, but downgradient of a new clay-lined landfill. These monitoring wells are close to both sites. Subsurface landfill gas migration was identified after a low permeability cap was installed on the older site. Subsurface gas pressure was monitored to identify horizontal landfill gas migration. Monitoring well headspace gases were evaluated to identify depressed oxygen concentrations and methane because of landfill gas migration into the well. Monitoring well headspace gas VOC concentrations were compared to groundwater VOC concentrations to determine the direction of phase transfer. A ratio above 1.0 of the observed well headspace gas concentration of a VOC to the concentration that would be in equilibrium with the groundwater concentration indicates gas-to-water phase transfer within the well. For the major gas-phase and aqueous-phase VOC, cis-1,2-dichloroethene, gas-to-water phase transfer is clearly indicated from the data for two of the four wells. Fifteen other VOCs were detected in monitoring well headspace gases but not in groundwater samples from the four wells studied. Only one compound in one well was detected in the groundwater sample but not in the headspace gases, and only one compound in one well was detected in both matrices at concentrations that suggested water-to-gas phase transfer. This study suggests that if landfill gas is suspected as the source of detected VOCs, monitoring well construction and stratigraphy are important considerations when attempting to differentiate between groundwater contamination by landfill gas and contamination from other sources.     

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.     

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.     

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m,p-xylene, and o-xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o-xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m,p-xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o-xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.     

Natural remobilization of multicomponent DNAPL pools due to dissolution

Mixtures of dense nonaqueous phase liquids (DNAPLs) trapped in the subsurface can act as long-term sources of contamination by dissolving into flowing groundwater. If the components have different solubilities then dissolution will alter the composition of the remaining DNAPL. We theorized that a multicomponent DNAPL pool may become mobile due to the natural dissolution process. In this study, we focused on two scenarios: (1) a DNAPL losing light component(s), with the potential for downward migration; and (2) a DNAPL losing dense component(s), with the potential for upward migration following transformation into a less dense than water nonaqueous phase liquid (LNAPL). We considered three binary mixtures of common groundwater contaminants: benzene and tetrachloroethylene (PCE), PCE and dichloromethane (DCM), and DCM and toluene. A number of physical properties that control the retention and transport of DNAPL in porous media were measured for the mixtures, namely: density, interfacial tension, effective solubility, and viscosity. All properties except density exhibited nonlinear relationships with changing molar ratio of the DNAPL. To illustrate the potential for natural remobilization, we modelled the following two primary mechanisms: the reduction in pool height as mass is lost by dissolution, and the changes in fluid properties with changing molar ratio of the DNAPL. The first mechanism always reduces the capillary pressure in the pool, while the second mechanism may increase the capillary pressure or alter the direction of the driving force. The difference between the rate of change of each determines whether the potential for remobilization increases or decreases. Static conditions and horizontal layering were assumed along with a one-dimensional, compositional modelling approach. Our results indicated that for initial benzene/PCE ratios greater than 25:75, the change in density was sufficiently faster than the decline in pool height to promote DNAPL breakthrough into the adjacent porous medium. In contrast, there was no potential for natural remobilization of a PCE-DCM mixture, primarily because the densities of the components are not sufficiently different. Dissolution of a DCM-toluene mixture decreased the density, reducing the tendency for downward displacement. However, the ultimate transformation from a DNAPL to an LNAPL may induce upward displacement. These results suggest that at sites with DNAPL pools containing a mix of components of sufficiently different densities and relative solubilities, natural remobilization may be an active mechanism, with implications for site evaluation and remediation.     

Determination of the Sodium Carbonate-Sulfur Dioxide Reaction Rate Coefficient at 450°F

Abstract

Gas-phase dispersion in granular biofilter materials with a wide range of particle sizes was investigated using atmospheric air and nitrogen as tracer gases. Two types of materials were used: (1) light extended clay aggregates (LECA), consisting of highly porous particles, and (2) gravel, consisting of solid particles. LECA is a commercial material that is used for insulation, as a soil conditioner, and as a carrier material in biofilters for air cleaning. These two materials were selected to have approximately the same particle shape. Column gas transport experiments were conducted for both materials using different mean particle diameters, different particle size ranges, and different gas flow velocities. Measured breakthrough curves were modeled using the advection-dispersion equation modified for mass transfer between mobile and immobile gas phases. The results showed that gas dispersivity increased with increasing mean particle diameter for LECA but was independent of mean particle diameter for gravel. Gas dispersivity also increased with increasing particle size range for both media. Dispersivities in LECA were generally higher than for gravel. The mobile gas content in both materials increased with increasing gas flow velocity but it did not show any strong dependency on mean particle diameter or particle size range. The relative fraction of mobile gas compared with total porosity was highest for gravel and lowest for LECA likely because of its high internal porosity.     

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).     

Effect of scale and dimensionality on the surfactant-enhanced solubilization of a residual DNAPL contamination

The mass transfer rate from residual dense non-aqueous phase liquids (DNAPLs) to the mobile aqueous phase is an important parameter for the efficiency of surfactant-enhanced remediation through solubilization of this type of contamination. The mass transfer kinetics are highly dependent on the dimensionality of the system. In this study, irregularly shaped residual TCE saturations in two-dimensional saturated flow fields were flushed with a 2% polyoxyethylene sorbitan (20) monooleate (POESMO) solution until complete removal had been achieved. A numerical model was developed and used for the simulation of the various surfactant-flushing experiments with different initial saturation patterns and flow rates. Through optimization against in situ concentration and saturation data, a phenomenological power-law model for the relationship between the mass transfer rate from the DNAPL to the mobile aqueous phase on the one hand and the residual DNAPL saturation and the flow velocity on the other hand was derived. The obtained mass transfer rate parameters provide a reasonable fit to the experimental data, predicting the cleanup time and the general saturation and concentration pattern quite well but failing to predict the concentration curves at every individual sampling port. The obtained mass transfer rate model gives smaller values for the predicted mass transfer rate but shows a comparable dependence on water flow and saturation as in earlier published one-dimensional column experiments with identical characteristics for porous medium, DNAPL and surfactant. Mass transfer rate predictions were about one order of magnitude lower in the 2-D flow cell experiment than in 1-D column experiments. These results give an indication for the importance of dimensionality during surfactant remediation.     

Modeling dissolution and volatilization of LNAPL sources migrating on the groundwater table

A vertically averaged two-dimensional model was developed to describe areal spreading and migration of light nonaqueous-phase liquids (LNAPLs) introduced into the subsurface by spills or leaks from underground storage tanks. The NAPL transport model was coupled with two-dimensional contaminant transport models to predict contamination of soil gas and groundwater resulting from a LNAPL migrating on the water table. Numerical solutions were obtained by using the finite-difference method. Simulations and sensitivity analyses were conducted with a LNAPL of pure benzene to study LNAPL migration and groundwater contamination. The model was applied to subsurface contamination by jet fuel. Results indicated that LNAPL migration were affected mostly by volatilization. The generation and movement of the dissolved plume was affected by the geology of the site and the free-product plume. Most of the spilled mass remained as a free LNAPL phase 20 years after the spill. The migration of LNAPL for such a long period resulted in the contamination of both groundwater and a large volume of soil.