Predictive modelling of NAPL injection tests in variable aperture spatially correlated fractures.

In preparation for a field experiment where a NAPL will be injected into a fractured sandstone aquifer, a 2D invasion percolation model of DNAPL migration in a single horizontal fracture with varying aperture has been developed. This simulation investigated the effect of spatially correlated fracture aperture fields on pressure-saturation relationships as a function of variable aperture mean, standard deviation, and spatial correlation statistics under hydrostatic conditions. Results from spatially correlated variable aperture fields can be significantly different from random fields. Longer ranges decreased entry pressures and higher standard deviations decreased nonwetting phase saturations. Mean aperture is the major control on displacement pressure, entry pressure and the form of the pressure-saturation curve. Simulation results using statistical parameters for a variable aperture natural sandstone fracture as described by Yeo et al. [International Journal of Rock Mechanics and Mining Sciences 35 (1998) 1051] closely resemble a Brooks-Corey pressure-saturation function, and exhibit a distinct entry pressure followed by a rapid increase in nonwetting phase saturation. Graphical estimates of entry pressure provide a good approximation of the critical aperture controlling the formation of a continuous nonwetting phase pathway in a variable aperture fracture. Consequently, we show that multiphase flow concepts developed for porous media can successfully be applied to variable aperture fractures. Entry pressure correlates well to the mean aperture in these simulations when using a Gaussian distribution of fracture aperture. Interpreted aperture distributions from NAPL injection experiments do not fit the true distribution well at low nonwetting phase saturations, but do at higher saturations above the entry pressure. Consequently, true, mechanical aperture variation within a fracture plane cannot be determined from NAPL injection experiments either in the field or laboratory.     

Laboratory-scale in situ chemical oxidation of a perchloroethylene pool using permanganate

In situ chemical oxidation (ISCO) is an emerging technology for the destruction of some chlorinated solvents present in subsurface environments. A laboratory investigation using a physical model was designed to assess the effectiveness of using permanganate as an oxidant to reduce the mass of a perchloroethylene (PCE) pool. The physical model was filled with silica sand overlying a silica flour base, simulating a two-dimensional saturated sand zone overlying a capillary barrier. PCE was introduced into the model so that it rested on top of the silica flour base, forming a dense nonaqueous phase liquid pool. The experimental methodology involved flushing the model with a permanganate solution for 146 days. During this period, measurements of chloride were used to assess the extent of pool oxidation. Before and after the oxidant flush, the quasi-steady state dissolution from the PCE pool was evaluated. Additionally, tracer studies were completed to assess changes in the flow field due to the oxidation process. At the termination of the experiment nine soil cores extracted from the model were used to detect the presence of MnO2 deposits and to quantify the mass of PCE remaining in the system. Excavation of the remaining material in the model revealed that the MnO2 distribution throughout the model was consistent with that observed in the cores. The oxidant flush was concluded before all of the pure phase PCE had been completely oxidized; however, approximately 45% of the PCE mass was removed, resulting in a fourfold decrease in the quasi-steady state aqueous phase mass loading of PCE from the pool. Measurements of chloride during the oxidant flush and of PCE in the soil cores suggested that the oxidation reaction occurred primarily at the upgradient edge of the PCE pool. MnO2 deposits within the model aquifer decreased the velocity of water directly above the pool, and the overall mass transfer from the remaining PCE pool. The results of this experimental study indicate that ISCO using permanganate is capable of removing substantial mass from a DNAPL pool; however, the performance of ISCO as a pool removal technology will be limited by the formation and precipitation of hydrous MnO2 that occurs during the oxidation process.     

Interphase mass transfer during chemical oxidation of TCE DNAPL in an aqueous system

The chemical oxidation of trichloroethene dense non-aqueous phase liquid by permanganate was studied in an aqueous system using micro-reaction/extraction vessels in a novel approach. Experiments were conducted at ambient temperature ( approximately 20 degrees C) under static and mixed conditions to evaluate the rate of TCE(DNAPL) dissolution as a function of permanganate concentration. Chemical oxidation by permanganate was shown to increase the rate of TCE(DNAPL) dissolution under static conditions and decrease the rate of dissolution under mixed conditions. The apparent inconsistency in results appears to result from the local deposition of a film at the DNAPL interface composed of manganese oxide solids as discovered through visual observation with the aid of a Goniometer. Data from interfacial deposition experiments suggest that the film formed rapidly and reached maturation within approximately 2 h with little or no growth occurring thereafter. A conceptual model of the reaction and mass transfer processes occurring at the DNAPL interface was proposed based on the experimental results.     

Multiphase flow and transport in fractured clay/sand sequences.

A numerical model (Queen's University Multi-Phase Flow Simulator, QUMPFS) was used to assess the rate of trichloroethylene (TCE) dense, non-aqueous phase liquid (DNAPL) migration through fractured clay, with special attention focused on the influence of interbedded sand lenses. The presence of these sand lenses was found to increase the time required for the non-wetting phase to migrate through the full 30 m vertical extent of the clay sequence from a few days to several years. Applied vertical hydraulic gradients were found to be moderately influential in systems consisting solely of fractured clays, yet one of the dominant factors controlling speed of vertical migration when sand lenses were present. Larger displacement pressure of the sands relative to that of the fractures leads to slower DNAPL migration rates, due to the delays that occur during build-up of capillary pressures. Dissolution of DNAPL and subsequent matrix diffusion of the aqueous phase has little effect on the rate of DNAPL migration through systems consisting of fractured clay only, yet slows the rate of migration in systems containing sand lenses. In all cases examined, the rate of DNAPL loading to the lower aquifer far exceeded the rate of aqueous phase mass loading. It was also found that DNAPL reaches the lower aquifer at approximately the same time as the aqueous phase plumes even for systems experiencing downward groundwater flow due to the attenuation of the aqueous phase through matrix diffusion.     

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.     

Using constant head step tests to determine hydraulic apertures in fractured rock

The initial step in the analysis of contaminant transport in fractured rock requires the consideration of groundwater velocity. Practical methods for estimating the average linear groundwater velocity (vˉ) in fractured rock require determination of hydraulic apertures which are commonly calculated by applying the cubic law using transmissivity (T) values and the number of hydraulically active fractures in the test interval. High-resolution, constant-head step injection testing of cored boreholes in a 100 m thick fractured dolostone aquifer was conducted using inflatable packers to isolate specific test intervals from the rest of the borehole. The steps in each test interval were gradually increased from very low to much higher injection rates. At smaller injection rates, the flow rate vs. applied pressure graph projects through the origin and indicates Darcian flow; non Darcian flow is evident at higher injection rates. Non-Darcian flow results in significantly lower calculated T values, which translates to smaller hydraulic aperture values. Further error in the calculated hydraulic aperture stems from uncertainty in the number of hydraulically active fractures in each test interval. This estimate can be inferred from borehole image and core logs, however, all of the fractures identified are not necessarily hydraulically active. This study proposes a method based on Reynolds number calculations aimed at improving confidence in the selection of the number of active fractures in each test interval.     

Two-phase flow and transport in a single fracture-porous medium system

The prognosis for the remediation of contaminated fractured media is much worse than that for more homogeneous units. Fractures act as conduits for the flow of dense non-aqueous phase liquids (DNAPLs), while diffusion is responsible for the storage of dissolved mass in the surrounding matrix. A numerical model incorporating aqueous phase transport in a variable-aperture fracture and its surrounding matrix is developed and coupled with an existing two-phase flow model. The processes of transient two-phase flow, non-equilibrium dissolution, advective–dispersive transport in the fracture, and three-dimensional matrix diffusion are included in the model. Results from various investigations show that the DNAPL distribution is very sensitive to variations in aperture within a single fracture. Diffusion-controlled mass removal from both the matrix and from the hydraulically inaccessible zones within the fracture itself result in extremely large time frames for significant mass removal from these systems. Success in aqueous phase mass removal from the matrix is very sensitive to the effective fracture spacing. The hydraulic gradient in the fracture only affects the amount of water removed from the system, and does not greatly affect the amount of time required to remove the contaminant mass from the matrix. The ability to remove mass is somewhat sensitive to the porosity and effective matrix diffusion coefficient over the range of expected values.     

Intermediate-scale 2D experimental investigation of in situ chemical oxidation using potassium permanganate for remediation of complex DNAPL source zones

In situ chemical oxidation is a technology that has been applied to speed up remediation of a contaminant source zone by inducing increased mass transfer from DNAPL sources into the aqueous phase for subsequent destruction. The DNAPL source zone can consist of one or more individual sources that may be present as an interconnected pool of high saturation, as a region of disconnected ganglia at residual saturation, or as combinations of these two morphologies. Potassium permanganate (KMnO(4)) is a commonly employed oxidant that has been shown to rapidly destroy DNAPL compounds like PCE and TCE following second-order kinetics in an aqueous system. During the oxidation of a target DNAPL compound, or naturally occurring reduced species in the subsurface, manganese oxide (MnO(2)) solids are produced. Research has shown that these manganese oxide solids may result in permeability reductions in the porous media thus reducing the ability for oxidant to be transported to individual DNAPL sources. It can also occur at the DNAPL-water interface, decreasing contact of the oxidant with the DNAPL. Additionally, MnO(2) formation at the DNAPL-water interface, and/or flow-bypassing as a result of permeability reductions around the source, may alter the mass transfer from the DNAPL into the aqueous phase, potentially diminishing the magnitude of any DNAPL mass depletion rate increase induced by oxidation. An experiment was performed in a two-dimensional (2D) sand-filled tank that included several discrete DNAPL source zones. Spatial and temporal monitoring of aqueous PCE, chloride, and permanganate concentrations was used to relate changes in mass depletion of, and mass flux, from DNAPL residual and pool source zones to chemical oxidation performance and MnO(2) formation. During the experiment, permeability changes were monitored throughout the 2D tank and these were related to MnO(2) deposition as measured through post-oxidation soil coring. Under the conditions of this experiment, MnO(2) formation was found to reduce permeability in and around DNAPL source zones resulting in changes to the overall flow pattern, with the effects depending on source zone configuration. A pool with little or no residual around it, in a relatively homogeneous flow field, appeared to benefit from resulting MnO(2) pore-blocking that substantially reduced mass transfer from the pool even though there was relatively little PCE mass removed from the pool. In contrast, a pool with residual around it (in a more typical heterogeneous flow field) appeared to undergo increased mass transfer as MnO(2) reduced permeability, altering the water flow and increasing the mixing at the DNAPL-water interface. Further, the magnitude of increased PCE mass depletion during oxidation appeared to depend on the PCE source configuration (pool versus ganglia) and decreased as MnO(2) was formed and deposited at the DNAPL-water interface. Overall, the oxidation of PCE mass appeared to be rate-limited by the mass transfer from the DNAPL to aqueous phase.     

Bench-scale visualization of DNAPL remediation processes in analog heterogeneous aquifers: surfactant floods and in situ oxidation using permanganate

We have conducted well-controlled DNAPL remediation experiments within a 2-D, glass-walled, sand-filled chamber using surfactants (Aerosol MA and Tween 80) to increase solubility and an oxidant (permanganate) to chemically degrade the DNAPL. Initial conditions for each remediation experiment were created by injecting DNAPL as a point source at the top of the chamber and allowing the DNAPL to migrate downward through a water-filled, heterogeneous, sand-pack designed to be evocative of a fluvial depositional environment. This migration process resulted in the DNAPL residing as a series of descending pools. Lateral advection across the chamber was used to introduce the remedial fluids. Photographs and digital image analysis illustrate interactions between the introduced fluids and the DNAPL. In the surfactant experiments, we found that DNAPL configured in a series of pools was easily mobilized. Extreme reductions in DNAPL/water interfacial tension occurred when using the Aerosol MA surfactant, resulting in mobilization into low permeability regions and thus confounding the remediation process. More modest reductions in interfacial tension occurred when using the Tween 80 surfactant resulting in modest mobilization. In this experiment, capillary forces remained sufficient to exclude DNAPL migration into low permeability regions allowing the excellent solubilizing properties of the surfactant to recover almost 90% of the DNAPL within 8.6 pore volumes. Injection of a potassium permanganate solution resulted in precipitation of MnO2, a reaction product, creating a low-permeability rind surrounding the DNAPL pools. Formation of this rind hindered contact between the permanganate and the DNAPL, limiting the effectiveness of the remediation. From these experiments, we see the value of performing visualization experiments to evaluate the performance of proposed techniques for DNAPL remediation.     

Measurements of groundwater velocity in discrete rock fractures

Estimating groundwater velocity in fracture networks using a Darcy or cubic law calculation is complicated by the wide distribution of fracture aperture often found in these systems and by the difficulty in measuring hydraulic head in discrete fracture features. Although difficult to conduct in a fractured rock setting, the point dilution method can be utilized to collect direct measurements of groundwater velocity in individual fractures. To compare measured against calculated velocities, more than 100 point dilution experiments were conducted within a 35 x 35 m area of a single fracture and in discrete fracture features within a fracture network at a larger scale. The dilution experiments were conducted by isolating a fracture feature in a borehole, measuring the hydraulic aperture, and measuring the decay of an injected tracer due to the advective groundwater flux across the fracture. Groundwater velocity was estimated using the hydraulic aperture and the rate of decay of the injected tracer. Estimates of the local hydraulic gradient were calculated via the cubic law using the velocity estimate and the hydraulic aperture. The results of the tests conducted in the single fracture show variable (1 to 33 m/day) but on average higher velocities in comparison to that measured during a natural gradient tracer experiment conducted previously (in which the effects of matrix diffusion were accounted for) and to that which would be calculated using the cubic law. Based on these results, it was determined that the best estimate of the average groundwater velocity, at the scale of the measurement area used for the cubic law calculations, could only be obtained using the largest apertures in the aperture distribution. Variability of the velocity measurements was also observed over time. Increases in velocity were attributed to the effect of rainfall although concurrent increases in hydraulic gradient were not detected (likely within the tolerance of the measuring devices). The groundwater velocities measured in the fracture network varied over a wider range than at the scale of the single fracture (from 2 to 388 m/day). No correlation, however, was observed between the size of the fracture aperture and measured velocity.     

Multiphase flow and transport through fractured heterogeneous porous media

The migration of Dense, Non-Aqueous Phase Liquid (DNAPL) and dissolved phase contamination through a fractured heterogeneous porous medium has been investigated through the use of a multiphase compositional model. The sensitivity of the timescales of migration and the distribution of contaminant in the subsurface to the mean permeability, the variance of the permeability, and the degree of fracturing of the domain were examined. It was found that increasing the mean permeability of the domain allowed the DNAPL to penetrate deeper into the subsurface, while decreasing the mean permeability caused the DNAPL to pool at shallower depths. The presence of fractures within the system was found to control the infiltration only in the most fractured domain. Moment analysis of the nonwetting phase showed that large-scale movement had ceased after approximately 9 years (maximum duration of the source-on condition was approximately 4.5 years). This tended to be due to a redistribution of the DNAPL towards a residual configuration, as was evidenced by the gradual trending of average nonwetting phase saturations within the domain to a static value. The dissolved phase plume was found to migrate at essentially the same rate as the nonwetting phase, due to the reduced relative permeability of lenses containing DNAPL, and due to diffusive losses of mass to the matrix of fractured clay and silty-clay lenses. Some exceptions to this were found when the DNAPL could not overcome the displacement pressure of a lens, and could not by-pass the lens due to the lack of available driving force after the source had been shut off.     

Cosolvent-enhanced chemical oxidation of perchloroethylene by potassium permanganate

A laboratory study was conducted to examine cosolvent-enhanced in-situ chemical oxidation (ISCO) of perchloroethylene (PCE) using potassium permanganate (KMnO4). The conceptual basis for this new technique is to enhance permanganate oxidation of dense non-aqueous phase liquids (DNAPLs) with the addition of a cosolvent, thereby increasing DNAPL solubility while avoiding mobilization. Among 17 cosolvent candidates screened, tertiary butyl alcohol (TBA) and acetone were the most stable in the presence of KMnO4, both of which increased PCE aqueous solubility significantly, and therefore are suitable to be used as cosolvent in this study. Batch experiments indicated that the second-order rate constant for PCE oxidation by potassium permanganate was 0.043+/-0.002 M(-1) s(-1) in the purely aqueous (no cosolvent) solution. In the presence of 20% cosolvent (volume fraction=fc=0.2), the rate constant decreased to 0.036+/-0.003 M(-1) s(-1) with TBA and to 0.031+/-0.002 M(-1) s(-1) with acetone. However, in the presence of free-phase PCE, chloride ion concentration from PCE oxidation in acetone/water solutions (fc=0.2) was about twice that in aqueous solutions, indicating that the increase in PCE solubility more than compensated for the decrease in reaction rate constant, such that the oxidation efficiency of PCE was increased with cosolvent. A complete chlorine mass balance was observed in the aqueous system, whereas approximately 70% was obtained in TBA/water or acetone/water (fc=0.2). In soil columns containing residual DNAPL and subjected to isocratic flushing with step-wise increases in f(c) cosolvent, TBA at fc=0.2 resulted in PCE mobilization, whereas acetone at fc相似文献   

Radionuclide transport and retardation in rock fracture and crushed rock column experiments

Transport and retardation of non-sorbing tritiated water and chloride and slightly sorbing sodium was studied in Syyry area SY-KR7 mica gneiss, in altered porous tonalite and in fresh tonalite. Experiments were performed using dynamic fracture and crushed rock column methods. Static batch method for sodium was introduced to compare retardation values from static and dynamic experiments. The ¹⁴C-PMMA method was used to study the pore structure of matrices. The pore aperture distribution was evaluated from Hg-porosimetry determinations and the surface areas were determined using the B.E.T. method. The flow characteristics and transport behavior of tracers were interpreted using a numerical compartment model for dispersion. The effect of matrix diffusion was calculated using an analytical solution to the advection-matrix diffusion problem in which surface retardation was taken into account. Radionuclide transport behavior in rock fractures was explained on the basis of rock structure.     

Long-term ground penetrating radar monitoring of a small volume DNAPL release in a natural groundwater flow field

An earlier field experiment at Canadian Forces Base Borden by Brewster and Annan [Geophysics 59 (1994) 1211] clearly demonstrated the capability of ground penetrating radar (GPR) reflection profiling to detect and monitor the formation of DNAPL layers in the subsurface. Their experiment involved a large volume release (770 L) of tetrachloroethylene into a portion of the sand aquifer that was hydraulically isolated from groundwater flow by sheet pile walls. In this study, we evaluated the ability of GPR profiling to detect and monitor much smaller volume releases (50 L). No subsurface confining structure was used in this experiment; hence, the DNAPL impacted zone was subjected to the natural groundwater flow regime. This condition allowed us to geophysically monitor the DNAPL mass loss over a 66 month period. Reflectivity variations on the GPR profiles were used to infer the presence and evolution of the solvent layers. GPR imaging found significant reflectivity increases due to solvent layer formation during the two week period immediately after the release. These results demonstrated the capacity of GPR profiling for the detection and monitoring of lesser volume DNAPL releases that are more representative of small-scale industrial spills. The GPR imaged solvent layers subsequently reduced in both areal extent and reflectivity after 29 months and almost completely disappeared by the end of the 66 month monitoring period. Total DNAPL mass estimates based on GPR profiling data indicated that the solvent mass was reduced to 34%-36% of its maximum value after 29 months; only 4%-9% of the solvent mass remained in the study area after 66 months. These results are consistent with independent hydrogeological estimates of remaining DNAPL mass based on the downgradient monitoring of the dissolved solvent phase. Hence, we have concluded that the long-term GPR reflectivity changes of the DNAPL layers are likely the result from the dissolution of chlorinated solvents residing in those layers. The long-term monitoring results demonstrated that GPR profiling is a promising non-invasive method for use at DNAPL contaminated sites in sandy aquifers where temporal information about immiscible contaminant mass depletion due to either natural flow or remediation is needed. However, our results also indicated that the GPR signature of older DNAPL impacted zones may not differ greatly from the uncontaminated background if significant mass reduction due to dissolution has occurred.     

Field study of TCE diffusion profiles below DNAPL to assess aquitard integrity

An area where a free-product accumulation of trichloroethylene (TCE) dense non-aqueous phase liquid (DNAPL) occurs at the bottom of a 10-m-thick surficial sand aquifer was studied to determine the integrity of the underlying, 20-m-thick, clayey silt aquitard formed of glaciolacustrine sediment. TCE concentration-versus-depth profiles determined from aquitard cores collected at five locations indicated penetration of detectable TCE 2.5 to 3.0 m into the aquitard. Two of the profiles show persistent DNAPL at the aquitard interface, while two others indicate that DNAPL, present initially, was completely dissolved away producing concentration declines at the aquitard interface. The fifth profile suggests shallow DNAPL penetration (<0.5 m) into the aquitard, however, this penetration, which was likely caused by cross-contamination during core collection or cone penetrometry (CPT) of the aquitard interface, did not increase the maximum depth of TCE penetration. Combining the field profiles with one-dimensional model simulations, downward migration of the aqueous TCE front, defined as the EPA MCL of 5 microg/l, which was below the analytical detection limit, was projected to a distance between 4 and 5 m below the top of the aquitard. Using a single set of estimated aquitard parameter values, simulations of aqueous TCE migration into the aquitard provided a good fit to four of the field profiles with a migration time of 35 to 45 years, consistent with the history of TCE use at the site. These simulations indicate aqueous TCE migration is diffusion-dominated with only small advective influence by the downward groundwater velocity of 2 to 3 cm/year or less in the aquitard due to pumping of the underlying aquifer to supply water to the facility in the past 50 years. The applicability of the parameter values was confirmed by in situ diffusion experiments of 1-year duration, in which stainless steel cylinders containing DNAPL were inserted into the aquitard. The diffusion-dominated nature of the profiles indicates that the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells, concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard provides long-term protection of the underlying aquifer from contamination from this DNAPL zone. Simulations of long-term migration of the TCE solute front indicate breakthrough to the lower aquifer at 1200 years for the no advection scenario and at 500 years if the strong downward hydraulic gradient persists. However, even after the breakthrough, the mass flux through the aquitard to the underlying aquifer remains relatively low, and when considered in terms of potential impacts to pumping wells , concentrations are not expected to increase significantly above present-day MCLs. The use of contaminant profiles of different time and distance scales, in addition to hydraulic data, dramatically improves the ability to assess aquitard integrity, and provides improved transport parameter values for estimating contaminant arrival times and fluxes. The apparent lack of deep preferential pathways for TCE migration, such as open fractures, is probably due to the softness of the silty aquitard deposit and minimal physical or chemical weathering of the aquitard.     

Evaluation of trichloroethene recovery processes in heterogeneous aquifer cells flushed with biodegradable surfactants

The ability of two biodegradable surfactants, polyoxyethylene (20) sorbitan monooleate (Tween 80) and sodium dihexyl sulfosuccinate (Aerosol MA), to recover a representative dense non-aqueous-phase liquid (DNAPL), trichloroethene (TCE), from heterogeneous porous media was evaluated through a combination of batch and aquifer cell experiments. An aqueous solution containing 3.3% Aerosol MA, 8% 2-propanol and 6 g/l CaCl(2) yielded a weight solubilization ratio (WSR) of 1.21 g TCE/g surfactant, with a corresponding liquid-liquid interfacial tension (IFT) of 0.19 dyn/cm. Flushing of aquifer cells containing a TCE-DNAPL source zone with approximately two pore volumes of the AMA formulation resulted in substantial (>30%) mobilization of TCE-DNAPL. However, a TCE mass recovery of 81% was achieved when the aqueous-phase flow rate was sufficient to displace the mobile TCE-DNAPL toward the effluent well. Aqueous solutions of Tween 80 exhibited a greater capacity to solubilize TCE (WSR=1.74 g TCE/g surfactant) and exerted markedly less reduction in IFT (10.4 dyn/cm). These data contradict an accepted empirical correlation used to estimate IFT values from solubilization capacity, and indicate a unique capacity of T80 to form concentrated TCE emulsions. Flushing of aquifer cells with less than 2.5 pore volumes of a 4% T80 solution achieved TCE mass recoveries ranging from 66 to 85%, with only slight TCE-DNAPL mobilization (<5%) occurring when the total trapping number exceeded 2 x 10(-5). These findings demonstrate the ability of Tween 80 and Aerosol MA solutions to efficiently recover TCE from a heterogeneous DNAPL source zone, and the utility of the total trapping number as a design parameter for a priori prediction of DNAPL mobilization and bank angle formation when flushing with low-IFT solutions. Given their potential to stimulate microbial reductive dechlorination at low concentrations, these surfactants are well-suited for remedial action plans that couple aggressive mass removal followed by enhanced bioremediation to treat chlorinated solvent source zones.     

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.     

Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation

This paper concludes that back diffusion from one or a few thin clayey beds in a sand aquifer can cause contaminant persistence above MCLs in a sand aquifer long after the source zone initially causing the plume is isolated or removed. This conclusion is based on an intensive case study of a TCE contaminated site in Florida, with the processes evaluated using numerical modeling. At this site, the TCE DNAPL zone formed decades ago, and was hydraulically isolated by means of an innovative system performing groundwater extraction, treatment and re-injection. Treated water is re-injected in a row of injection wells situated a short distance downgradient of the extraction wells, creating a clean-water displacement front to efficiently flush the downgradient plume. This scheme avoids the creation of stagnation zones typical of most groundwater pump-and-treat systems, thereby minimizing the time for aquifer flushing and therefore downgradient cleanup. The system began operation in August 2002 and although the performance monitoring shows substantial declines in concentrations, detectable levels of TCE and degradation products persist downgradient of the re-injection wells, long after the TCE should have disappeared based on calculations assuming a nearly homogenous sand aquifer. Three hypotheses were assessed for this plume persistence: 1) incomplete source-zone capture, 2) DNAPL occurrence downgradient of the re-injection wells, and 3) back diffusion from one or more thin clay beds in the aquifer. After careful consideration, the first two hypotheses were eliminated, leaving back diffusion as the only plausible hypothesis, supported by detailed measurements of VOC concentrations within and near the clay beds and also by numerical model simulations that closely represent the field site hydrogeologic conditions. The model was also used to simulate a more generalized, hypothetical situation where more thin clayey beds occur in a sand aquifer with an underlying aquitard. While there is no doubt that DNAPL source mass reduction can eventually improve downgradient groundwater quality, the magnitude and time scale over which the improvement occurs is the major uncertainty given current characterization approaches. This study shows that even one thin clay bed, less than 0.2 m thick, can cause plume persistence due to back diffusion for several years or even decades after the flux from the source is completely isolated. Thin clay beds, which have a large storage capacity for dissolved and sorbed contaminant mass, are common in many types of sandy aquifers. However, without careful inspection of continuous cores and sampling, such thin clay beds, and their potential for causing long-term back-diffusion effects, can easily go unnoticed during site characterization.     

Identification of key parameters controlling dissolved oxygen migration and attenuation in fractured crystalline rocks

In the crystalline rocks of the Canadian Shield, geochemical conditions are currently reducing at depths of 500-1000 m. However, during future glacial periods, altered hydrologic conditions could potentially result in enhanced recharge of glacial melt water containing a relatively high concentration of dissolved oxygen (O2). It is therefore of interest to investigate the physical and geochemical processes, including naturally-occurring redox reactions, that may control O2 ingress. In this study, the reactive transport code MIN3P is used in combination with 2k factorial analyses to identify the most important parameters controlling oxygen migration and attenuation in fractured crystalline rocks. Scenarios considered are based on simplified conceptual models that include a single vertical fracture, or a fracture zone, contained within a rock matrix that extends from the ground surface to a depth of 500 m. Consistent with field observations, Fe(II)-bearing minerals are present in the fractures (i.e. chlorite) and the rock matrix (biotite and small quantities of pyrite). For the parameter ranges investigated, results indicate that for the single fracture case, the most influential factors controlling dissolved O2 ingress are flow velocity in the fracture, fracture aperture, and the biotite reaction rate in the rock matrix. The most important parameters for the fracture zone simulations are flow velocity in the individual fractures, pO2 in the recharge water, biotite reaction rate, and to a lesser degree the abundance and reactivity of chlorite in the fracture zone, and the fracture zone width. These parameters should therefore receive increased consideration during site characterization, and in the formulation of site-specific models intended to predict O2 behavior in crystalline rocks.     

Attempt to model laboratory-scale diffusion and retardation data

Different approaches for measuring the interaction between radionuclides and rock matrix are needed to test the compatibility of experimental retardation parameters and transport models used in assessing the safety of the underground repositories for the spent nuclear fuel. In this work, the retardation of sodium, calcium and strontium was studied on mica gneiss, unaltered, moderately altered and strongly altered tonalite using dynamic fracture column method. In-diffusion of calcium into rock cubes was determined to predict retardation in columns. In-diffusion of calcium into moderately and strongly altered tonalite was interpreted using a numerical code FTRANS. The code was able to interprete in-diffusion of weakly sorbing calcium into the saturated porous matrix. Elution curves of calcium for the moderately and strongly altered tonalite fracture columns were explained adequately using FTRANS code and parameters obtained from in-diffusion calculations. In this paper, mass distribution ratio values of sodium, calcium and strontium for intact rock are compared to values, previously obtained for crushed rock from batch and crushed rock column experiments. Kd values obtained from fracture column experiments were one order of magnitude lower than Kd values from batch experiments.