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.     

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.     

表面活性剂强化抽出处理含水层中DNAPL污染物的去除特征

为明确表面活性剂强化抽出处理含水层中DNAPL污染物过程中表面活性剂的增强修复效果,及DNAPL自身理化性质和介质孔径的影响,利用数码图像分析技术对1,2-二氯乙烷和四氯乙烯2种DNAPL在石英砂填充的二维砂箱中的抽取迁移过程进行了实验模拟研究,并对抽出水样中DNAPL的浓度进行了测试分析。结果表明,实验条件下加入低浓度(0.18%)的十二烷基苯磺酸钠(SDBS)大幅提高了对弱透水层截留的2种DNAPL聚集体的抽出处理效率。1,2-二氯乙烷在该表面活性剂溶液中的表观溶解度远高于四氯乙烯,因此其短时间内的绝对去除率更高。SDBS强化抽出处理DNAPL的作用机理以增溶作用为主,而其增流作用使DNAPL迁移流动后分布面积增大,增加了与表面活性剂溶液接触的面积,对增溶作用起到促进效果。细粒介质中DNAPL迁移后的最大分布面积较大,因此体系中DNAPL的溶解速率较高。在DNAPL聚集体质量与水力梯度固定的条件下,油水界面张力越低,DNAPL的密度越大,DNAPL垂向迁移的风险就越大。本研究为修复工程中如何依据DNAPL种类与场地多孔介质的情况选择表面活性剂提供了参考。     

Removal of DNAPL contamination from the saturated zone by the combined effect of vertical upward flushing and density reduction

A common aspect of innovative remediation techniques is that they tend to reduce the interfacial tension between the aqueous and non-aqueous phase liquids, resulting in mobilization of the organic contaminant. This complicates the remediation of aquifers, contaminated with Dense Non-Aqueous Phase Liquids (DNAPLs), as they are likely to migrate downwards, deeper into the aquifer and into finer layers. A possible solution is the use of swelling alcohols, which tend to reduce the density difference between the aqueous phase and the DNAPL. To avoid premature mobilization upon the initial contact between the DNAPL and the alcohol, several researchers have proposed the use of vertical upward flow of the alcohol. In this paper, we present an equation, which describes the upward mobilization of both continuous and discontinuous DNAPLs and so the important parameters governing the upward controlled mobilization of the DNAPL. The need and required magnitude of this specific discharge was investigated by conducting four column experiments in which the initial density of the DNAPL and the permeability was varied. It was shown that the required flow velocities increase with the permeability of the porous medium and the initial density difference between the aqueous phase and the DNAPL. Whenever the specific discharge falls below the critical value, the DNAPL moves downward. A second set of column experiments looked at the impact of permeability of porous medium on the solubilization and mobilization of DNAPL during alcohol flooding. Columns, packed with coarse or fine sand, containing a residual trichloroethylene (TCE) or perchloroethylene (PCE) saturation were flushed with the alcohol mixture at a fixed specific discharge rate. The induced pressure gradients in the aqueous phase, which were higher in the fine sand, resulted for this porous medium in extensive mobilization of the DNAPL against the direction of the buoyancy force. The density of the first NAPL coming out of the top of the fine sand was close to that of the pure DNAPL. In the coarser sand, the pressure gradients were sufficient to prevent downward migration of the DNAPL, but upward mobilization was minimal. The predominant removal mechanism in this case was the much slower solubilization.     

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.     

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.     

Relative velocities of DNAPL and aqueous phase plume migration

Numerical simulation is used to examine the relative velocities of DNAPL and aqueous phase plumes in sandy aquifers where lateral spreading of DNAPL has occurred at the base of the aquifer. The scenario being modeled is one where a permeable aquifer is underlain by a sloping aquitard, which results in lateral migration of the DNAPL down the slope, in addition to lateral migration of an aqueous phase plume subject to a specified hydraulic gradient. A sensitivity analysis is presented to the impacts of both DNAPL properties and geologic properties. The most important chemical properties governing the relative velocities of the DNAPL and the shallow aqueous phase plume are the DNAPL viscosity and the aqueous component soil-water partition coefficient (Kd). The dip of the underlying aquitard was found to be relatively unimportant, at least for the range of values studied. The scenario under consideration can be important in conceptual model development and remedial design, as in certain cases DNAPL could be migrating in areas without the evidence of a well-developed aqueous phase plume. The implication of this work is that the absence of a shallow aqueous phase plume directly downgradient of a DNAPL source zone does not rule out the possibility of deep occurrences of DNAPL beyond the shallow monitoring well network. A further finding of this study is that the occurrence of a highly sorbing compound in groundwater at virtually any concentration may indicate the immediate upgradient presence of residual or pooled DNAPL.     

Mass removal of chlorinated ethenes from rough-walled fractures using permanganate

In situ chemical oxidation (ISCO) employing permanganate is an emerging technology that has been successful at enhancing mass removal from DNAPL source zones in unconsolidated media at the pilot-scale. The focus of this study was to evaluate the applicability of flushing a permanganate solution across two single vertical fractures in a laboratory environment to remove free phase DNAPL. The fracture experiments were designed to represent a portion of a larger fractured aquifer system impacted by a near-surface DNAPL spill over a shallow fractured rock aquifer. Each fracture was characterized by hydraulic and tracer tests, and the aperture field for one of the fractures was mapped using a co-ordinate measurement machine. Following DNAPL emplacement, a series of water and permanganate flushes were performed. To support observations from the fracture experiments, a set of batch experiments was conducted. The data from both fracture experiments showed that the post-oxidation effluent concentration was not impacted by the oxidant flush; however, changes in the aperture distribution, flow field, and flow rate were observed. These changes resulted in a significant decrease to the mass loading from the fractures, and were attributed to the build-up of oxidation by-products (manganese oxides and carbon dioxide) within the fracture which was corroborated by the batch experiment data and visual examination of the walls of one fracture. These results provide insight into the potential impact that a permanganate solution and oxidation by-products can have on the aperture distribution within a fracture and on DNAPL mass transfer rates. A permanganate flush or injection completed within a fractured rock aquifer may lead to the development of an insoluble product adjacent to the DNAPL which results in the reduction or complete elimination of advective regions near the DNAPL and reduces mass transfer rates. This outcome would have significant implications on the plume generating potential of the remaining DNAPL.     

Removal of PCB-DNAPL from a rough-walled fracture using alcohol/polymer flooding

Phase behaviour experiments employing PCB (Aroclor 1242)/alcohol/water systems were conducted with ethanol (EtOH) and n-propanol (nPA). Both exhibited an affinity for the aqueous phase within the entire two-phase region. As much as 88% by volume (88% vol.) EtOH and 80% vol. nPA were necessary to achieve full miscibility of the PCB in the aqueous phase. DNAPL-water interfacial tension (IFT) was reduced from 38.9 dyn/cm to 4.7 dyn/cm and 2.4 dyn/cm with 80% vol. EtOH and 76% vol. nPA. The addition of alcohol brought about 41% and 54% reductions in DNAPL viscosity at maximal concentrations of EtOH and nPA. Density of the PCB-DNAPL was relatively unaffected by the presence of alcohol. A series of seven experiments were conducted where successive slugs of nPA and xanthan gum polymer solutions were injected into a fractured shale sample. A 30% vol. nPA solution injected under a hydraulic gradient of 0.36 allowed enhanced PCB removal primarily through reduction of IFT and resulted in 72% DNAPL recovery. Several pore volumes of alcohol solution were necessary to displace all the potentially mobile non-wetting phase since the high-viscosity DNAPL was mobilized at a lower flow rate than the overall fluid velocity, illustrating non-piston displacement. The injection of a 95% vol. nPA alcohol solution, theoretically at a sufficient concentration to produce fully miscible displacement of the residual DNAPL at equilibrium, resulted in non-equilibrium partitioning of the PCB into the flushing solution, likely due to the high fluid velocities in the fracture. The injection of 200 pore volumes of 95% vol. nPA solution resulted in 94% DNAPL recovery. Alcohol floods operated below the miscibility envelope appear to be a valuable source zone remedial alternative where the objective is to reduce DNAPL mobility to zero, but it should be noted that DNAPL mobility is increased during the application of the technology and steps may need to be taken to prevent unwanted vertical mobilization.     

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.     

A new concept: the use of neutrally-buoyant microemulsions for DNAPL remediation

Even in the absence of mobilization of dense nonaqueous phase liquid (DNAPL), the microemulsion that forms when the surfactant solubilizes a dense contaminant such as trichloroethylene will be more dense than water and tends to migrate downward. This paper addresses the issue of migration with a new concept: surfactant enhanced aquifer remediation at neutral buoyancy. Laboratory results of surfactant remediation in two-dimensional model aquifers show that downward migration of microemulsion containing solubilized dense contaminants can be reduced to an acceptable level, even in the absence of capillary barriers in the aquifer. One model experiment was designed to exhibit a small degree of vertical migration and full capture of the microemulsion at the extraction well. The second experiment was designed to demonstrate the effect of large buoyancy forces that lead to excessive downward migration of the microemulsion. Density measurements of aqueous solutions containing sodium dihexyl sulfosuccinate surfactant, isopropanol, trichloroethylene, and sodium chloride are presented. A companion paper presents the results of the flow and transport calculations needed for this approach to surfactant flooding.     

Microbial physiology-based model of ethanol metabolism in subsurface sediments

Dense non-aqueous phase liquids (DNAPLs) present in the subsurface may contain surface active compounds that impact DNAPL migration and distribution. While a number of studies have revealed the role surface active compounds play in altering the wettability of quartz sand, few have considered the implications for other minerals common to contaminated sites. This study extends understanding of DNAPL/surfactant wettability to iron oxide surfaces. Specifically, quartz and iron oxide-coated sands in a tetrachloroethene (PCE)/water system containing the organic base (an organic molecule that acts as a base) dodecylamine (DDA) were compared at a variety of scales. Wettability of the minerals' surfaces, and the impact of wettability on capillary resistance to DNAPL entry, were assessed as a function of pH through: (i) advancing and receding contact angles, (ii) primary drainage capillary pressure-saturation experiments, and (iii) small, two-dimensional, flow cell experiments. The work revealed that, at neutral pH and under identical boundary capillary pressures, DNAPL invaded quartz sand but not iron oxide-coated sand; however, at low pH, DNAPL invaded both sands equally. These differences were demonstrated to be due to wettability alterations associated with the strength of attractive forces between DDA and the mineral surface, dictated by the isolectric point of the minerals and system pH. Observed differences in DNAPL invasion behavior were consistent with measured intrinsic contact angles and P(c)-S relationships, the latter requiring scaling by the operative contact angle inside the porous medium for a meaningful comparison. This study suggests that the distribution of minerals (and, more specifically, their isoelectric points), as well as the aqueous phase pH at a given site, may have a significant impact on the DNAPL source zone architecture.     

Using polymer solutions to enhance recovery of mobile coal tar and creosote DNAPLs

Direct pumping and enhanced recovery of coal tar and creosote dense, non-aqueous phase liquids (DNAPLs) from the subsurface have had mixed results because these DNAPLs are viscous fluids that can potentially alter aquifer wettability. To improve the inefficiencies associated with waterflooding, the research presented here considered the use of a polymer solution that can be added to the injected flood solution to increase the viscosity and decrease the velocity of the flooding solution. Results from one-dimensional, vertically oriented laboratory column experiments that evaluate the recovery of coal-derived DNAPL with both water and polymer flooding solutions are presented. The final DNAPL saturation remaining in the column was assessed in water and oil-wet systems for three viscous DNAPLs. Adding polymer to increase the aqueous solution viscosity did not have a significant impact in water-wet systems. A final DNAPL saturation of approximately 19% was achieved for both water and polymer floods. In contrast, the addition of polymer significantly improved recovery in oil-wet systems. The final saturation was over 40% in oil-wet systems after waterflooding, but approximately 19% with a polymer flushing solution. Although the final saturation produced with polymer flooding was similar between the oil- and water-wet systems, differences in the relative permeability and distribution of DNAPL in the porous matrix caused the DNAPL recovery to be much slower in the oil-wet system.     

Noninvasive monitoring of DNAPL migration through a saturated porous medium using electrical impedance tomography

Electrical impedance tomography (EIT) was used to monitor the movement of a fluorinated hydrocarbon dense nonaqueous phase liquid (DNAPL) through a saturated porous medium within a laboratory column. Impedance measurements were made using a horizontal plane of 12 electrodes positioned at regular intervals around the centre of the column. A 2D inversion algorithm, which incorporated the cylindrical geometry of the column, was used to reconstruct resistivity and phase images from the measured data. Differential time-lapse images of DNAPL movement past the plane of electrodes were generated by the cell-by-cell subtraction of resistivity and phase baseline models from those associated with the DNAPL release stage of the experiment.The DNAPL pulse was clearly delineated as resistive anomalies in the differential time-lapse resistivity images. The spatial extent of the resistive anomalies indicated that in addition to vertical migration, some lateral spreading of the DNAPL had occurred. Residual contamination could be detected after quasi-static conditions were reestablished. Residual DNAPL saturation was estimated from the resistivity model data by applying Archie's second equation.Despite significant measured phase changes due to DNAPL contamination, the differential phase images revealed only weak anomalies associated with DNAPL flow; these anomalies could be seen only in the initial stages of the experiment during peak flow through the plane of electrodes.     

PCE dissolution and simultaneous dechlorination by nanoscale zero-valent iron particles in a DNAPL source zone

While the capability of nanoscale zero-valent iron (NZVI) to dechlorinate organic compounds in aqueous solutions has been demonstrated, the ability of NZVI to remove dense non-aqueous phase liquid (DNAPL) from source zones under flow-through conditions similar to a field scale application has not yet been thoroughly investigated. To gain insight on simultaneous DNAPL dissolution and NZVI-mediated dechlorination reactions after direct placement of NZVI into a DNAPL source zone, a combined experimental and modeling study was performed. First, a DNAPL tetrachloroethene (PCE) source zone with emplaced NZVI was built inside a small custom-made flow cell and the effluent PCE and dechlorination byproducts were monitored over time. Second, a model for rate-limited DNAPL dissolution and NZVI-mediated dechlorination of PCE to its three main reaction byproducts with a possibility for partitioning of these byproducts back into the DNAPL was formulated. The coupled processes occurring in the flow cell were simulated and analyzed using a detailed three-dimensional numerical model. It was found that subsurface emplacement of NZVI did not markedly accelerate DNAPL dissolution or the DNAPL mass-depletion rate, when NZVI at a particle concentration of 10g/L was directly emplaced in the DNAPL source zone. To react with NZVI the DNAPL PCE must first dissolve into the groundwater and the rate of dissolution controls the longevity of the DNAPL source. The modeling study further indicated that faster reacting particles would decrease aqueous contaminant concentrations but there is a limit to how much the mass removal rate can be increased by increasing the dechlorination reaction rate. To ensure reduction of aqueous contaminant concentrations, remediation of DNAPL contaminants with NZVI should include emplacement in a capture zone down-gradient of the DNAPL source.     

Infiltration of PCE in a system containing spatial wettability variations

A two-dimensional infiltration experiment was conducted to investigate and quantify the effect of spatial wettability variations on DNAPL migration and entrapment in saturated sands. Experimental observations of tetrachloroethylene (PCE) infiltration showed that organic-wet sand lenses acted as very effective capillary barriers, retaining PCE and inhibiting its downward migration. A multiphase numerical simulator was used to model this sand box experiment. The simulator incorporates wettability-modified van Genuchten and Brooks-Corey capillary pressure/saturation relationships as well as Burdine and Mualem relative permeability relationships. PCE mass distributions, estimated by image analysis of digital photographs taken during the infiltration event, were compared to simulation results. Although both relative permeability models were qualitatively able to predict the PCE retention in the organic-wet layers, simulations with the Mualem model failed to capture the observed rate of PCE migration. A traditional multiphase simulator, incorporating water-wet capillary retention relations, failed to predict both PCE pathways and retention behavior. This study illustrates the potential influence of subsurface wettability variations on DNAPL migration and entrapment and supports the use of modified capillary relations in conjunction with the Burdine model in multiphase flow simulators.     

Aquifer washing by micellar solutions: 1: Optimization of alcohol–surfactant–solvent solutions

Phase diagrams were used for the formulation of alcohol–surfactant–solvent and to identify the DNAPL (Dense Non Aqueous Phase Liquid) extraction zones. Four potential extraction zones of Mercier DNAPL, a mixture of heavy aliphatics, aromatics and chlorinated hydrocarbons, were identified but only one microemulsion zone showed satisfactory DNAPL recovery in sand columns. More than 90 sand column experiments were performed and demonstrate that: (1) neither surfactant in water, alcohol–surfactant solutions, nor pure solvent can effectively recover Mercier DNAPL and that only alcohol–surfactant–solvent solutions are efficient; (2) adding salts to alcohol–surfactant or to alcohol–surfactant–solvent solutions does not have a beneficial effect on DNAPL recovery; (3) washing solution formulations are site specific and must be modified if the surface properties of the solids (mineralogy) change locally, or if the interfacial behavior of liquids (type of oil) changes; (4) high solvent concentrations in washing solutions increase DNAPL extraction but also increase their cost and decrease their density dramatically; (5) maximum DNAPL recovery is observed with alcohol–surfactant–solvent formulations which correspond to the maximum solubilization in Zone C of the phase diagram; (6) replacing part of surfactant SAS by the alcohol n-butanol increases washing solution efficiency and decreases the density and the cost of solutions; (7) replacing part of n-butanol by the nonionic surfactant HOES decreases DNAPL recovery and increases the cost of solutions; (8) toluene is a better solvent than D-limonene because it increases DNAPL recovery and decreases the cost of solutions; (9) optimal alcohol–surfactant–solvent solutions contain a mixture of solvents in a mass ratio of toluene to D-limonene of one or two. Injection of 1.5 pore volumes of the optimal washing solution of n-butanol–SAS–toluene–D-limonene in water can recover up to 95% of Mercier DNAPL in sand columns. In the first pore volume of the washing solution recovered in the sand column effluent, the DNAPL is in a water-in-oil microemulsion lighter than the excess aqueous phase (Winsor Type II system), which indicates that part of the DNAPL was mobilized. In the next pore volumes, DNAPL is dissolved in a oil-in-water microemulsion phase and is mobilized in an excess oil phase lighter than the microemulsion (Winsor Type I system). The main drawback of this oil extraction process is the high concentration of ingredients necessary for DNAPL dissolution, which makes the process expensive. Because mobilization of oil seems to occur at the washing solution front, an injection strategy must be developed if there is no impermeable limit at the aquifer base. DNAPL recovery in the field could be less than observed in sand columns because of a smaller sweep efficiency related to field sand heterogeneities. The role of each component in the extraction processes in sand column as well as the Winsor system type have to be better defined for modeling purposes. Injection strategies must be developed to recover ingredients of the washing solution that can remain in the soil at the end of the washing process. ©1997 Elsevier Science B.V.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Prediction and minimization of vertical migration of DNAPLS using surfactant enhanced aquifer remediation at neutral buoyancy

Surfactant enhanced aquifer remediation (SEAR) had previously been thought to require a capillary barrier below the contaminated zone to prevent the vertical migration of a microemulsion containing solubilized DNAPL. This paper shows the vertical migration of the dense microemulsion is described and predicted by the value of three dimensionless scaling groups. Embedded within these scaling groups are four design parameters. The value of these parameters can be manipulated in order to reduce the amount of vertical migration anticipated for a given remediation design. Plots have been constructed that illustrate the relationship between vertical migration and the value of the scaling groups; such plots can be used to predict vertical migration and to determine appropriate screen intervals of extraction wells to ensure full capture of the contaminants. This predictive capability has been verified in laboratory experiments. Predicted migration of the microemulsion agreed within about 2% of that observed. Development of the scaling groups is presented, remediation design implications are discussed, and laboratory verification is shown. Additional discussion of the laboratory work is given in a companion paper [Kostarelos, K., Pope, G.A., Rouse, B.A., Shook, G.M., 1998. A new concept: the use of neutrally-buoyant microemulsions for DNAPL remediation, J. Contam. Hydro., this edition].     

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.