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.     

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.     

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.     

DNAPL remediation with in situ chemical oxidation using potassium permanganate. II. Increasing removal efficiency by dissolving Mn oxide precipitates

In situ chemical oxidation (ISCO) schemes using MnO4- have been effective in destroying chlorinated organic solvents dissolved in ground water. Laboratory experiments and field pilot tests reveal that the precipitation of Mn oxide, one of the reaction products, causes a reduction of permeability, which can lead to flow bypassing and inefficiency of the scheme. Without a solution to this problem of plugging, it is difficult to remove DNAPL from the subsurface completely. In a companion paper, we showed with batch experiments that Mn oxide can be dissolved rapidly with certain organic acids. This study utilizes 2-D flow-tank experiments to examine the possibility of nearly complete DNAPL removal by ISCO with MnO4-, when organic acids are used to remove Mn oxide. The experiments were conducted in a small 2-D glass flow tank containing a lenticular silica-sand medium. Blue-dyed trichloroethylene (TCE) provided residual, the perched and pooled DNAPL. KMnO4 at 200 mg/l was flushed through the DNAPL horizontally. Once plugging reduced permeability and prevented further delivery of the oxidant, citric or oxalic acids were pumped into the flow tank to dissolve the Mn oxide precipitates. Organic ligands removed the Mn oxide precipitates relatively quickly, and permitted another cycle of MnO4- flooding. Cycles of MnO4-/acid flooding continued until all of the visible DNAPL was removed. The experiments were monitored with chemical analysis and visualization. A mass-balance calculation indicated that by the end of the experiments, all the DNAPL was removed. The results show also how heterogeneity adds complexity to initial redistribution of DNAPL, and to the efficiency of the chemical flooding.     

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.     

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.     

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.     

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

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

Characterization of DNAPL from the U.S. DOE Savannah River Site

The composition of chlorinated hydrocarbon DNAPLs (dense non-aqueous phase liquids) from field sites can be substantially different than the material originally purchased for use as a solvent. Waste management practices at the U.S. Department of Energy's (DOE) Savannah River Site (SRS) included co-disposal of a wide range of organic and inorganic wastes. In 1991, a clear, orange-colored DNAPL was found in two wells near the SRS M-area settling basin. Waste effluent from the fuel and target fabrication facilities that were discharged to this settling basin included acids, caustics, metals and chlorinated solvents. The characterization of the SRS DNAPL suggests that numerous constituents partitioned into the DNAPL during its use as a solvent, co-disposal and ultimate migration through the subsurface. Trace constituents in the DNAPL include metals, from processing operations or co-disposal practices and subsurface minerals, high molecular weight hydrocarbons and alkyl esters, and acids. This complex mixture results in DNAPL-water interfacial properties that are substantially different than would be expected from a simple mixture of PCE and TCE. Under conditions when there is a high DNAPL to water volume ratio, a semi-rigid film accumulates on water droplets suspended in the DNAPL. It is concluded that the array of precipitated metal species comprising this film contributes to the interfacial tension that is over an order of magnitude lower than expected for a "clean" PCE/TCE mixture.     

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.     

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.     

Multidimensional validation of a numerical model for simulating a DNAPL release in heterogeneous porous media

A fixed-volume release of 1,2-DCE, tracked in space and time with a light transmission/image analysis system, provided a data set for the infiltration, redistribution, and immobilisation of a dense non-aqueous phase liquid (DNAPL) in a heterogeneous porous medium. The two-dimensional bench scale flow cell was packed with a spatially correlated, random heterogeneous distribution of six sand types. In order to provide the necessary modelling parameters, detailed constitutive relationships were measured at the local scale for the six sands. These experiments revealed that nonwetting phase (NWP) relative permeability-saturation (k(rN)-S(W)) relationships are strongly correlated to sand type. Trends in the best-fit k(rN)-S(W) parameters reflected a positive correlation between mean grain diameter and the maximum NWP relative permeability, k(rN)(max). Multiphase flow simulations of the bench scale experiment best reproduced the experimental observations, producing excellent matches in both time and space, when the measured, correlated local scale k(rN)-S(W) relationships were employed.     

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.     

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.     

Variations in expression of carbon isotope fractionation of chlorinated ethenes during biologically enhanced PCE dissolution close to a source zone

The stable carbon isotope values of tetrachloroethene (PCE) and its degradation products were monitored during studies of biologically enhanced dissolution of PCE dense nonaqueous phase liquid (DNAPL) to determine the effect of PCE dissolution on observed isotope values. The degradation of PCE was monitored in a 2-dimensional model aquifer and in a pilot test cell (PTC) at Dover Air Force Base, both with emplaced PCE DNAPL sources. Within the plume down gradient from the source, the isotopic fractionation of dissolved PCE and its degradation products were consistent with those observed in biodegradation laboratory studies. However, close to the source zone significant shifts in the isotope values of dissolved PCE were not observed in either the model aquifer or PTC due to the constant input of newly dissolved, non fractionated PCE, and the small isotopic fractionation associated with PCE reductive dechlorination by the mixed microbial culture used. Therefore the identification of reductive dechlorination in the presence of PCE DNAPL was based upon the appearance of daughter products and the isotope values of those daughter products. An isotope model was developed to simulate isotope values of PCE during the dissolution and degradation of PCE adjacent to a DNAPL source zone. With the exception of very high degradation rate constants (>1/day) stable carbon isotope values of PCE estimated by the model remained within error of the isotope value of the PCE DNAPL, consistent with measured isotope values in the model aquifer and in the PTC.     

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.     

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.     

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.     

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.