Nonideal transport of reactive contaminants in heterogeneous porous media: 7. distributed-domain model incorporating immiscible-liquid dissolution and rate-limited sorption/desorption

The purpose of this work is to present a distributed-domain mathematical model incorporating the primary mass-transfer processes that mediate the transport of immiscible organic liquid constituents in water-saturated, locally heterogeneous porous media. Specifically, the impact of grain/pore-scale heterogeneity on immiscible-liquid dissolution and sorption/desorption is represented in the model by describing the system as comprising a continuous distribution of mass-transfer domains. With this conceptualization, the distributions of the initial dissolution rate coefficient and the sorption/desorption rate coefficient are represented as probability density functions. Several sets of numerical experiments are conducted to examine the effects of heterogeneous dissolution and sorption/desorption on contaminant transport and elution. Four scenarios with different combinations of uniform/heterogeneous rate-limited dissolution and uniform/heterogeneous rate-limited sorption/desorption are evaluated. The results show that both heterogeneous rate-limited sorption/desorption and heterogeneous rate-limited dissolution can significantly increase the time or pore volumes required to elute immiscible-liquid constituents from a contaminated porous medium. However, sorption/desorption has minimal influence on elution behavior until essentially all of the immiscible liquid has been removed. For typical immiscible-liquid constituents that have relatively low sorption, the asymptotic elution tailing produced by heterogeneous rate-limited sorption/desorption begins at effluent concentrations that are several orders of magnitude below the initial steady-state concentrations associated with dissolution of the immiscible liquid. Conversely, the enhanced elution tailing associated with heterogeneous rate-limited dissolution begins at concentrations that are approximately one-tenth of the initial steady-state concentrations. Hence, dissolution may generally control elution behavior of immiscible-liquid constituents in cases wherein grain/pore-scale heterogeneity significantly influences both dissolution and sorption/desorption.     

Mass-removal and mass-flux-reduction behavior for idealized source zones with hydraulically poorly-accessible immiscible liquid

A series of flow-cell experiments was conducted to investigate aqueous dissolution and mass-removal behavior for systems wherein immiscible liquid was non-uniformly distributed in physically heterogeneous source zones. The study focused specifically on characterizing the relationship between mass flux reduction and mass removal for systems for which immiscible liquid is poorly accessible to flowing water. Two idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. The results showed that significant reductions in mass flux occurred at relatively moderate mass-removal fractions for all systems. Conversely, minimal mass flux reduction occurred until a relatively large fraction of mass (>80%) was removed for the control experiment, which was designed to exhibit ideal mass removal. In general, mass flux reduction was observed to follow an approximately one-to-one relationship with mass removal. Two methods for estimating mass-flux-reduction/mass-removal behavior, one based on system-indicator parameters (ganglia-to-pool ratio) and the other a simple mass-removal function, were used to evaluate the measured data. The results of this study illustrate the impact of poorly accessible immiscible liquid on mass-removal and mass-flux processes, and the difficulties posed for estimating mass-flux-reduction/mass-removal behavior.     

Relationship between mass-flux reduction and source-zone mass removal: analysis of field data

The magnitude of contaminant mass-flux reduction associated with a specific amount of contaminant mass removed is a key consideration for evaluating the effectiveness of a source-zone remediation effort. Thus, there is great interest in characterizing, estimating, and predicting relationships between mass-flux reduction and mass removal. Published data collected for several field studies were examined to evaluate relationships between mass-flux reduction and source-zone mass removal. The studies analyzed herein represent a variety of source-zone architectures, immiscible-liquid compositions, and implemented remediation technologies. There are two general approaches to characterizing the mass-flux-reduction/mass-removal relationship, end-point analysis and time-continuous analysis. End-point analysis, based on comparing masses and mass fluxes measured before and after a source-zone remediation effort, was conducted for 21 remediation projects. Mass removals were greater than 60% for all but three of the studies. Mass-flux reductions ranging from slightly less than to slightly greater than one-to-one were observed for the majority of the sites. However, these single-snapshot characterizations are limited in that the antecedent behavior is indeterminate. Time-continuous analysis, based on continuous monitoring of mass removal and mass flux, was performed for two sites, both for which data were obtained under water-flushing conditions. The reductions in mass flux were significantly different for the two sites (90% vs. approximately 8%) for similar mass removals ( approximately 40%). These results illustrate the dependence of the mass-flux-reduction/mass-removal relationship on source-zone architecture and associated mass-transfer processes. Minimal mass-flux reduction was observed for a system wherein mass removal was relatively efficient (ideal mass-transfer and displacement). Conversely, a significant degree of mass-flux reduction was observed for a site wherein mass removal was inefficient (non-ideal mass-transfer and displacement). The mass-flux-reduction/mass-removal relationship for the latter site exhibited a multi-step behavior, which cannot be predicted using some of the available simple estimation functions.     

Source-zone characterization of a chlorinated-solvent contaminated Superfund site in Tucson, AZ

An extensive site-characterization project was conducted at a large chlorinated-solvent contaminated Superfund site in Tucson, AZ. The project consisted of several components, including traditional site-characterization activities, tracer tests, laboratory experiments conducted with core material collected from the site, and mathematical modeling. The primary focus of the work presented herein is the analysis of induced-gradient contaminant elution tests conducted in a source zone at the site, investigation of the potential occurrence of immiscible liquid in the saturated zone, characterization of the relationship between mass flux reduction and mass removal, and evaluation of the impact of source-zone management on site remediation. The results of the present study, along with those of prior work, indicate that immiscible liquid is likely present in the saturated zone at the site source zones. Extensive tailing and rebound was observed for the contaminant-elution tests, indicating nonideal transport and mass-transfer behavior. The elution data were analyzed with a source-zone-scale mathematical model, and the results indicated that nonideal immiscible-liquid dissolution was the primary cause of the observed behavior. The time-continuous relationship between mass flux reduction and mass removal associated with the plume-scale pump-and-treat operation exhibited an initial large drop in mass flux with minimal mass removed, followed by a period of minimal mass flux reduction and a second period of large reduction. This behavior reflects the impact of both source-zone and aqueous-plume mass removal dynamics. Ultimately, a greater than 90% reduction in mass flux was achieved for a mass removal of approximately 50%. The influence of source-zone management on site remediation was evaluated by conducting two predictive simulations, one for which the source zones were controlled and one for which they were not. A plume-scale model was used to simulate the composite contaminant concentrations associated with groundwater extracted with the pump-and-treat system, which were compared to measured data. The information generated from this study was used to enhance the site conceptual model, help optimize operation of the pump-and-treat system, and evaluate the utility of source-zone remediation.     

The effects of matrix diffusion on radionuclide migration in rock column experiments

Rock column experiments were performed to examine the effects of matrix diffusion and hydrodynamic dispersion on the migration of radionuclides at the laboratory scale. Tritiated water and chloride transportation was studied in intact mica gneiss and in altered more porous tonalite columns with narrow flow channels. The column diffusion properties were estimated prior to water flow experiments using the gas diffusion method with helium as the tracer gas. The numerical compartment model for advection and dispersion, with and without matrix diffusion, was used to interpret the tracer transport in the columns. Matrix diffusion behavior was also distinguished from dominating hydrodynamic dispersion in rock column experiments at the slowest water flow rates.     

A micromodel analysis of factors influencing NAPL removal by surfactant foam flooding

A methodology to study the trichloroethylene (TCE) and dodecane removal in porous media by surfactant foams (SF) was presented by using etched-glass micromodels. The purpose of this work was to systematically evaluate the impact of various physicochemical factors such as gas fraction (GF), surfactant concentration, pore structure and nonaqueous phase liquid (NAPL) types on NAPL removal during SF flooding. The TCE displacement by SF was dependent on the gas fraction of SF. Low GFs (50% and 66%) were more efficient for TCE removal and sweep efficiencies than a high GF (85%). An increase in TCE removal was observed with increasing surfactant concentration at a fixed GF. TCE removal by SF flooding appeared to be dependent more to the value of Capillary number rather than to the concentration of surfactant solution. The effect of the pore heterogeneity was evaluated by employing two different types of micromodels. The Capillary number is an important parameter in the determination of sweep efficiency or gas saturation of SF in a nonhomogeneous porous medium. However, the TCE removal from a nonhomogeneous porous medium may not be associated with sweep efficiency. The initial configuration of residual TCE blobs in a nonhomogeneous porous medium would also be influential in displacing TCE. Sweep efficiencies and pressure responses of two NAPL systems (TCE and dodecane) were monitored to evaluate foam stability when the foam contacts the NAPLs. Stable foam contacting with TCE is implied, while it appears that dodecane cause the SF to collapse. All results indicate that the Capillary number (a ratio of viscous forces to capillary forces) is the most important parameter for TCE removal by SF flooding. Micromodel visualizations of water, surfactant and SF floods were showed and also discussed.     

Infiltration and redistribution of perchloroethylene in partially saturated, stratified porous media

Contamination of the subsurface by nonaqueous phase liquids (NAPLs) is a widespread problem. To investigate the behavior of a nonspreading, dense NAPL (DNAPL) in the vadose zone, we conducted perchloroethylene (PCE) infiltration experiments in nominally 1- and 2-dimensional (D), stratified porous media. In addition, the usefulness and limitations of a multifluid flow simulator to describe PCE infiltration and redistribution under the experimental conditions were tested. The physical simulations were conducted in a column (1-D) and a flow container (2-D) which were packed with two distinct layers of coarse-grained sand and a fine-grained sand layer in between. Volumetric water and PCE contents were determined with a fully automated dual-energy gamma radiation system. While migrating through the drier parts of the coarse-grained sand layers, PCE appeared to wet the water–air interface rather than displacing any water. In the wetter parts of the porous medium, PCE displaced water and behaved as a true nonwetting fluid. PCE showed a limited response to gradients in capillary pressure and rather high values for the volumetric PCE content were measured in the fine-grained sand layers. This was attributed to the nonspreading nature of PCE. The multifluid flow simulator appeared to predict the initial PCE movement in the vadose zone reasonably well. However, the model was not capable of predicting the final amounts of PCE retained in either the unsaturated or saturated part of the flow domain, mainly because the simulator does not consider the nonspreading flow behavior of NAPLs.     

In-situ oxidation of trichloroethene by permanganate: effects on porous medium hydraulic properties

In-situ oxidation of dense nonaqueous-phase liquids (DNAPLs) by strong oxidants such as potassium permanganate (KMnO4) has been proposed as a possible DNAPL remediation strategy. In this study, we investigated the effects of in-situ trichloroethene (TCE) oxidation by KMnO4 on porous medium hydraulic properties. In particular, we wanted to determine the overall effects of concurrent solid phase (MnO2) precipitation, gas (CO2) evolution and TCE dissolution resulting from the oxidation reaction on the porous medium's aqueous-phase relative permeability, krw. Three TCE removal experiments were conducted in a 95-cm long, 5.1-cm i.d. glass column, which was homogeneously packed with well-characterized 30/40-mesh silica sand. TCE was emplaced in the sand-pack in residual, entrapped form through a sequence of water/TCE imbibition and drainage steps. The column was then flushed under constant aqueous flux conditions for up to 104 h with either deionized water (reference experiment), deionized water containing 5 mM KMnO4 or deionized water containing 5 mM KMnO4 and 300 mM Na2HPO4. Aqueous-phase relative permeabilities were computed from measured flow rates and measurements of aqueous-phase pressure head, h obtained using pressure transducers connected to tensiometers distributed along the column length. A dual-energy gamma radiation system was used to monitor changes in fluid saturation that occurred during each experiment. In addition, column effluent samples were collected for chemical analyses. Dissolution of TCE during deionized water flushing led to an increase in krw by approximately 22% and a local reduction in h. On the other hand, vigorous CO2 gas production and precipitation of MnO2 was visually observed during flushing with deionized water that contained 5 mM KMnO4. As a consequence, krw declined by approximately 96% and h increased locally by more than 1000 cm H2O during the first 24 h of the experiment, causing sand-pack ruptures and pump failure. Conversely, less CO2 gas production and MnO2 precipitation was visually observed during flushing with deionized water that contained 5 mM KMnO4 and 300 mM Na2HPO4. Consequently, only small increases in h (< 15 cm H2O) were observed in this experiment due to a reduction in krw of approximately 53%. While we must attribute changes in h due to variations in krw to our specific experimental design (constant aqueous flux, one-dimensional flow experiments), these experiments nevertheless confirm that successful application of in situ chemical oxidation of TCE requires consideration of detrimental processes such as MnO2 precipitation and CO2 gas formation. In addition, our results indicate that utilization of a buffered oxidant solution may improve the effectiveness of in-situ oxidation of TCE by KMnO4 in otherwise weakly buffered porous media.     

Immiscible fluid flow in porous media: dielectric properties

The movement and retention of immiscible fluids in porous media was determined using a microwave dielectric measuring technique. The technique was used to monitor the composite dielectric constant of columns containing model soil materials to observe the effect of displacement of various saturating fluids by water and determine the retained immiscible components. Various combinations of displacing fluids, both miscible and immiscible, were studied, yielding time-dependent flow characteristics as well as the total retention factor. The flow circumstances were similar to what would be encountered in an accidental spill of chlorocarbons or crude oil.The retention of chlorocarbons and crude oil was determined for model soils consisting of sand or glass spheres. The storage capacity of the sand for immiscible fluids was dependent on the fluid initially saturating the porous medium. For chlorocarbons, it was increased in water-wet compared to dry sands while the inverse was true for crude oil. The retention of crude oil in the sand was dependent on the water flow velocity and was much larger than that for chlorocarbons. The crude oil could be displaced by solvents such as cyclohexane in the sand but the net volume occupied by the two immiscible fluids after water flushing did not change appreciably, i.e. the crude oil component was replaced by the solvent. Application of a miscible fluid-water mixture or surfactants improved the ability of water to displace the immiscible fluid from the sand column.     

The effect of local-scale physical heterogeneity and nonlinear,rate-limited sorption/desorption on contaminant transport in porous media

Nonideal transport of contaminants in porous media has often been observed in laboratory characterization studies. It has long been recognized that multiple processes associated with both physical and chemical factors can contribute to this nonideal transport behavior. To fully understand system behavior, it is important to determine the relative contributions of these multiple factors when conducting contaminant transport and fate studies. In this study, the relative contribution of physical-heterogeneity-related processes versus those of nonlinear, rate-limited sorption/desorption to the observed nonideal transport of trichloroethene in an undisturbed aquifer core was determined through a series of miscible-displacement experiments. The results of experiments conducted using the undisturbed core, collected from a Superfund site in Tucson, AZ, were compared to those obtained from experiments conducted using the same aquifer material packed homogeneously. The results indicate that both physical and chemical factors, specifically preferential flow and associated rate-limited diffusive mass-transfer and rate-limited sorption/desorption, respectively, contributed to the nonideal behavior observed for trichloroethene transport in the undisturbed core. A successful prediction of trichloroethene transport in the undisturbed core was made employing a mathematical model incorporating multiple sources of nonideal transport, using independently determined model parameters to account for the multiple factors contributing to the nonideal transport behavior. The simulation results indicate that local-scale physical heterogeneity controlled the nonideal transport behavior of trichloroethene in the undisturbed core, and that nonlinear, rate-limited sorption/desorption were of secondary importance.     

活性炭吸附有机蒸气过程中的初始浓度测定

以PTA生产尾气为实验体系，讨论了在活性炭对有机气体动态吸附过程中，有机气体初始浓度的三种测定方法，即气相色谱法、吸附称量法和气体方程计算法。结果得出，吸附称量法误差最小，其次是气体方程计算法，误差最大的是气相色谱法。     

Assessment of human exposure to PCDDs,PCDFs and Co-PCBs using hair as a human pollution indicator sample I: Development of analytical method for human hair and evaluation for exposure assessment

Dioxins including polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs) and co-planar polychlorinated biphenyls (Co-PCBs) are highly toxic. Even at environmental pollution levels, they cause hormonal damage in women, and they have been shown to induce immunosuppression and genital function damage in humans. In this study, a new method using isotope dilution was established to detect PCDDs, PCDFs and Co-PCBs in human hair. This method, comprised of washing and cutting of hair, alkaline decomposition, hexane extraction, multilayer silica gel column chromatography, high performance liquid chromatography with a porous graphite carbon column and analysis by high resolution gas chromatography/high resolution mass spectrometry, enabled us to analyze PCDDs, PCDFs and Co-PCBs at trace levels of less than pg/g with good reproducibility. In addition, there was a correlation between some isomers in human hair and blood collected from identical donors. Human hair analysis is useful to evaluate human risk assessment including that due to environmental pollution.     

Biodegradation during contaminant transport in porous media: 4. Impact of microbial lag and bacterial cell growth

Miscible-displacement experiments were conducted to examine the impact of microbial lag and bacterial cell growth on the transport of salicylate, a model hydrocarbon compound. The impacts of these processes were examined separately, as well as jointly, to determine their relative effects on biodegradation dynamics. For each experiment, a column was packed with porous medium that was first inoculated with bacteria that contained the NAH plasmid encoding genes for the degradation of naphthalene and salicylate, and then subjected to a step input of salicylate solution. The transport behavior of salicylate was non-steady for all cases examined, and was clearly influenced by a delay (lag) in the onset of biodegradation. This microbial lag, which was consistent with the results of batch experiments, is attributed to the induction and synthesis of the enzymes required for biodegradation of salicylate. The effect of microbial lag on salicylate transport was eliminated by exposing the column to two successive pulses of salicylate, thereby allowing the cells to acclimate to the carbon source during the first pulse. Elimination of microbial lag effects allowed the impact of bacterial growth on salicylate transport to be quantified, which was accomplished by determining a cell mass balance. Conversely, the impact of microbial lag was further investigated by performing a similar double-pulse experiment under no-growth conditions. Significant cell elution was observed and quantified for all conditions/systems. The results of these experiments allowed us to differentiate the effects associated with microbial lag and growth, two coupled processes whose impacts on the biodegradation and transport of contaminants can be difficult to distinguish.     

Laboratory-scale application of fiber optic transflection dip probe (FOTDP) for in situ monitoring of gas phase ozone in unsaturated porous media

A fiber optic transflection dip probe (FOTDP) system was developed for in situ and real-time monitoring of the transport of gas phase ozone in unsaturated porous media. A unique property of this system is the employment of a dip probe, which is inserted within the porous media. At the probe's tip, incoming light interacts with gas phase ozone and is partially reflected back into the probe by a mirror attached to the tip. Calibration of the FOTDP system was successfully carried out with various ozone concentrations using a column packed with glass beads. The ozone breakthrough curves (BTCs) were obtained by converting normalized UV intensities into gas phase ozone concentrations. The FOTDP system worked well for in situ monitoring of gas phase ozone using a column packed with sand under various water saturations in the presence of SOM and reflected the ideal transport phenomena of gas phase ozone for various flow rates.     

Influence of surfactant-facilitated interfacial tension reduction on chlorinated solvent migration in porous media: observations and numerical simulation

The ability of a multiphase flow model to capture the migration behavior of chlorinated solvents under conditions of surfactant-facilitated interfacial tension (IFT) reduction is assessed through comparison of model predictions with observations from controlled laboratory experiments. Tetrachloroethene (PCE) was released in two-dimensional saturated systems, packed with sandy media that incorporated rectangular lenses of capillary contrast. Spatially uniform interfacial tension conditions were created in the tanks by pre-flushing the porous medium with either Milli Q water or an aqueous surfactant solution. Experimental observations showed that surfactant-facilitated IFT reductions substantially lowered capillary resistance to the vertical downward migration of PCE and enabled PCE to enter finer grained, less permeable lenses that were not penetrated in the absence of surfactant. An immiscible flow model was used to simulate the conditions of the laboratory experiments. Under higher IFT conditions (47.5 and 5 dyn/cm), the model could successfully predict the general migration behavior of the organic liquid. Model predictions, however, exhibited poorer agreement with observed migration pathways under low IFT conditions (0.5 dyn/cm). In all cases, the predicted PCE distributions were influenced by selection of the parametric model for capillary retention and relative permeability. Simulated migration rates were more consistent with observed behavior when the Brooks-Corey/Burdine model was employed. For low interfacial tensions, improved predictions of migration pathways were obtained through grid refinement and incorporation of small-scale packing variability. Simulations highlight the substantial sensitivity of model predictions to the capillary pressure-scaling factor, grid resolution, and small-scale porosity variations at interfaces of permeability contrast under reduced IFT conditions.     

Effect of nonionic surfactant partitioning on the dissolution kinetics of residual perchloroethylene in a model porous medium

At concentrations above the critical micelle concentration, surfactants can significantly enhance the solubilization of residual nonaqueous phase liquids (NAPL) and, for this reason, are the focus of research on surfactant-enhanced aquifer remediation (SEAR). As a consequence of their amphiphilic nature, surfactants may also partition to various extents between the organic and aqueous phases, thereby affecting SEAR performance. We report here on the observation and analysis of the effect of surfactant partitioning on the dissolution kinetics of residual perchloroethylene (PCE) by aqueous solutions (1000 mg/L) of the non-ionic surfactant Triton X-100 in a model porous medium. For this fluid system, batch equilibration experiments showed that the surfactant partitions strongly into the NAPL (NAPL-water partition coefficient equal to 12.5). Dynamic interfacial tension (IFT) measurements were employed to study surfactant diffusion and interfacial adsorption. The dynamic IFT measurements were consistent with partitioning of the surfactant between the two liquid phases. PCE dissolution experiments, conducted in a transparent glass micromodel using an aqueous surfactant solution, were contrasted to experiments using clean water. Surfactant partitioning was observed to delay significantly the onset of micellar solubilization of PCE, an observation reproduced by a numerical model. This effect is attributed to the reduction of surfactant concentration in the immediate vicinity of the NAPL-water interface, which accompanies transport of the surfactant into the NAPL. Accordingly, it is suggested that both the rate and the extent of diffusion of the surfactant into the NAPL affect the onset of and the driving force for micellar solubilization. While many surfactants do not readily partition in NAPL, this possibility must be considered when selecting non-ionic surfactants for the enhanced solubilization of residual chlorinated solvents in porous media.     

Transport of volatile compounds in porous media in the presence of a trapped gas phase

The presence of an immobile gaseous phase in an otherwise-saturated porous medium affects the transport of volatile compounds. The linear theory of partitioning tracers suggests that a volatile tracer introduced into such a system should be retarded with a constant retardation factor. Using high concentrations, however, the saturation of the gaseous phase will change as an effect of the tracer test itself. Competitive gas transfer among all volatile compounds and the change of saturation may lead to tracer concentrations that are temporarily higher than those injected. We analyze the system in the framework of the coherence theory by Helfferich [Soc. Pet. Eng. J. 21 (1) (1981) 51]. The governing equations are formulated as functions of total concentration, i.e., the mass of solutes in all phases per pore volume. Neglecting dispersion and mass-transfer kinetics, we derive the characteristic form of the resulting system of hyperbolic equations. In a system with N volatile compounds, a variation of the concentrations splits up into N waves, each traveling with its own characteristic velocity. If the presence of a gaseous phase is sustained, one wave will be a standing one. We perform numerical model calculations for tracers with various Henry's law coefficients and show that the results agree with the semi-analytical solution obtained by coherence theory.     

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.