Modeling in situ ozonation for the remediation of nonvolatile PAH-contaminated unsaturated soils

Mathematical models were developed to investigate the characteristics of gaseous ozone transport under various soil conditions and the feasibility of in situ ozone venting for the remediation of unsaturated soils contaminated with phenanthrene. On the basis of assumptions for the mass transfer and reactions of ozone, three approaches were considered: equilibrium, kinetic, and lump models. Water-saturation-dependent reactions of gaseous ozone with soil organic matter (SOM) and phenanthrene were employed. The models were solved numerically by using the finite-difference method, and the model parameters were determined by using the experimental data of Hsu [The use of gaseous ozone to remediate the organic contaminants in the unsaturated soils, PhD Thesis, Michigan State Univ., East Lansing, MI, 1995]. The transport of gas-phase ozone is significantly retarded by ozone consumption due to reactions with SOM and phenanthrene, in addition to dissolution. An operation time of 156 h was required to completely remove phenanthrene in a 5-m natural soil column. In actual situations, however, the operation time is likely to be longer than the ideal time because of unknown factors including heterogeneity of the porous medium and the distribution of SOM and contaminant. The ozone transport front length was found to be very limited (< 1 m). The sensitivity analysis indicated that SOM is the single most important factor affecting in situ ozonation for the remediation of unsaturated soil contaminated with phenanthrene. Models were found to be insensitive to the reaction mechanisms of phenathrene with either gas-phase ozone or dissolved ozone. More study is required to quantify the effect of OH* formation on the removal of contaminant and on ozone transport in the subsurface.     

Investigations on mobility of carbon colloid supported nanoscale zero-valent iron (nZVI) in a column experiment and a laboratory 2D-aquifer test system

Nanoscale zero-valent iron (nZVI) has recently gained great interest in the scientific community as in situ reagent for installation of permeable reactive barriers in aquifer systems, since nZVI is highly reactive with chlorinated compounds and may render them to harmless substances. However, nZVI has a high tendency to agglomerate and sediment; therefore it shows very limited transport ranges. One new approach to overcome the limited transport of nZVI in porous media is using a suited carrier colloid. In this study we tested mobility of a carbon colloid supported nZVI particle “Carbo-Iron Colloids” (CIC) with a mean size of 0.63 μm in a column experiment of 40 cm length and an experiment in a two-dimensional (2D) aquifer test system with dimensions of 110?×?40?×?5 cm. Results show a breakthrough maximum of 82 % of the input concentration in the column experiment and 58 % in the 2D-aquifer test system. Detected residuals in porous media suggest a strong particle deposition in the first centimeters and few depositions in the porous media in the further travel path. Overall, this suggests a high mobility in porous media which might be a significant enhancement compared to bare or polyanionic stabilized nZVI.     

Nitrogen as an indicator of mass transfer during in-situ gas sparging

Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.     

Empirical correlations to estimate agglomerate size and deposition during injection of a polyelectrolyte-modified Fe0 nanoparticle at high particle concentration in saturated sand

Controlled emplacement of polyelectrolyte-modified nanoscale zerovalent iron (NZVI) particles at high particle concentration (1-10 g/L) is needed for effective in situ subsurface remediation using NZVI. Deep bed filtration theory cannot be used to estimate the transport and deposition of concentrated polyelectrolyte-modified NZVI dispersions (>0.03 g/L) because particles agglomerate during transport which violates a fundamental assumption of the theory. Here we develop two empirical correlations for estimating the deposition and transport of concentrated polyelectrolyte-modified NZVI dispersions in saturated porous media when NZVI agglomeration in porous media is assumed to reach steady state quickly. The first correlation determines the apparent stable agglomerate size formed during NZVI transport in porous media for a fixed hydrogeochemical condition. The second correlation estimates the attachment efficiency (sticking coefficient) of the stable agglomerates. Both correlations are described using dimensionless numbers derived from parameters affecting deposition and agglomeration in porous media. The exponents for the dimensionless numbers are determined from statistical analysis of breakthrough data for polyelectrolyte-modified NZVI dispersions collected in laboratory scale column experiments for a range of ionic strength (1, 10, and 50mM Na(+) and 0.25, 1, and 1.25 mM Ca(2+)), approach velocity (0.8 to 55 × 10(-4)m/s), average collector sizes (d(50)=99 μm, 300 μm, and 880 μm), and polyelectrolyte surface modifier properties. Attachment efficiency depended on approach velocity and was inversely related to collector size, which is contrary to that predicted from classic filtration models. High ionic strength, the presence of divalent cations, lower extended adsorbed polyelectrolyte layer thickness, decreased approach velocity, and a larger collector size promoted NZVI agglomeration and deposition and thus limited its mobility in porous media. These effects are captured quantitatively in the two correlations developed. The application and limitations of using the correlations for preliminary design of in situ NZVI emplacement strategies is discussed.     

Investigation of the vadose zone using barometric pressure cycles

This paper documents a technique for investigating one-dimensional airflow in the vadose zone. Variations in pore gas pressures resulting from barometric cycles were measured at depths as great as 180 m in several gas monitoring wells. The data were transformed to the frequency domain, enabling comparison with closed-form analytic expressions of one-dimensional pressure transport in layered porous media. The data reveal evidence for vertical fracture flow that was not apparent from in situ measurements of permeability. The data also reveal that the basalt underlying the site at depths greater than 100 m has permeability exceeding 1000 darcies, and is vented to the atmosphere at an estimated distance of a few kilometers from the site.     

Effects of pH and ionic strength on sulfamethoxazole and ciprofloxacin transport in saturated porous media

Many antibiotics regarded as emerging contaminants have been frequently detected in soils and groundwater; however, their transport behaviors in soils remain largely unknown. This study examined the transport of two antibiotics, sulfamethoxazole (SMZ) and ciprofloxacin (CIP), in saturated porous media. Laboratory columns packed with quartz sand was used to test the effects of solution pH and ionic strength (IS) on their retention and transport. The results showed that these two antibiotics behaved differently in the saturated sand columns. In general, SMZ manifested a much higher mobility than CIP for all experimental conditions tested. Almost all SMZ transported through the columns within one pore volume in deionized water (i.e., pH=5.6, IS=0), but no CIP was detected in the effluents under the same condition after extended column flushing. Perturbations in solution pH (5.6 and 9.5) and IS (0 and 0.1M) showed no effect on SMZ transport in the saturated columns. When pH increased to 9.5, however, ~93% of CIP was eluted from the sand columns. Increase of IS from 0 to 0.1M also slightly changed the distribution of adsorbed CIP within the sand column at pH 5.6, but still no CIP was detected in the effluents. A mathematical model based on advection-dispersion equation coupled with equilibrium and kinetic reactions successfully simulated the transport of the antibiotics in water-saturated porous media with R(2)=0.99.     

臭氧水中传质模型的研究

在传统鼓泡塔中对臭氧在水中的传质过程进行实验研究。利用实验改变进水流量、臭氧进气流量以及臭氧进气浓度等，得到臭氧传质系数、臭氧传质效率和臭氧消耗等有关参数。建立一个臭氧传质模型，可以预测不同操作条件下臭氧传质效率，从而优化控制臭氧化反应动力学。     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

The velocity of DNAPL fingering in water-saturated porous media: laboratory experiments and a mobile-immobile-zone model.

Dense nonaqueous phase liquids (DNAPLs) are immiscible with water and can give rise to highly fingered fluid distributions when infiltrating through water-saturated porous media. In this paper, a conceptual mobile-immobile-zone (MIZ) model is presented to describe the structure of a DNAPL finger in water-saturated porous media and the velocity of finger propagation. A finger is composed of a finger body and a tip. The finger body has a mobile core and an immobile sheath. All the DNAPL within the tip of a finger is mobile. Lab experiments utilizing image analyses of a DNAPL (PCE) penetrating into water-saturated homogeneous glass beads were carried out in a two-dimensional transparent chamber. The results show that the fingers elongated almost linearly with time. The fingers did not grow laterally after the tip of the finger had passed. The average finger diameters were between 3.9 and 5.4 mm for PCE propagation in water-saturated glass bead porous media with mean particle diameters from 0.32 to 1.36 mm. The estimated mobile core diameters were 51-60% of the average finger diameters.     

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.     

The velocity of DNAPL fingering in water-saturated porous media: laboratory experiments and a mobile–immobile–zone model

Dense nonaqueous phase liquids (DNAPLs) are immiscible with water and can give rise to highly fingered fluid distributions when infiltrating through water-saturated porous media. In this paper, a conceptual mobile–immobile–zone (MIZ) model is presented to describe the structure of a DNAPL finger in water-saturated porous media and the velocity of finger propagation. A finger is composed of a finger body and a tip. The finger body has a mobile core and an immobile sheath. All the DNAPL within the tip of a finger is mobile. Lab experiments utilizing image analyses of a DNAPL (PCE) penetrating into water-saturated homogeneous glass beads were carried out in a two-dimensional transparent chamber. The results show that the fingers elongated almost linearly with time. The fingers did not grow laterally after the tip of the finger had passed. The average finger diameters were between 3.9 and 5.4 mm for PCE propagation in water-saturated glass bead porous media with mean particle diameters from 0.32 to 1.36 mm. The estimated mobile core diameters were 51–60% of the average finger diameters.     

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

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

堆肥 + 零价铁可渗透反应墙修复黄土高原地下水中铬铅复合污染简

为了研究堆肥+零价铁混合可渗透反应墙（PRB）修复黄土高原地下水中铬铅复合污染的可行性,分别用堆肥、零价铁、堆肥+ 零价铁、堆肥+ 零价铁+活性炭为反应介质,通过模拟柱实验考察PRB修复铬铅复合污染黄土高原地下水的效果。结果表明,在实验进行30 d后当反应柱1和2对六价铬的去除率接近于零,而且对二价铅的去除率迅速下降时,反应柱3对2种污染物仍保持较高的去除率;反应介质质量比为10：2：1的反应柱4和质量比为10：1：2的反应柱5对污染物的去除效果均优于质量比为10：1：1的反应柱3;反应50 d后,添加活性炭的反应柱6对2种污染物的去除率仍在90%。这说明使用堆肥+零价铁混合可渗透反应墙修复黄土高原地下水中铬铅复合污染是可行的;且以堆肥+零价铁作为介质的反应柱去除效果优于单独以堆肥或铁粉为介质的反应柱;增加铁粉或堆肥的用量有利于铬铅复合污染的去除;且同时添加活性炭更有助于污染物的去除。     

填料塔中Na2CO3吸收高浓度H2S的传质特性

在填料吸收塔中考察了Na2CO3溶液吸收高浓度H2S气体的气液传质特性。通过测量填料塔进出口气体中H2S浓度计算了Na2CO3溶液吸收高浓度H2S气体的总体积传质系数（KGa）,并研究了进气流速、吸收液流量、吸收温度和吸收液浓度对KGa的影响。结果表明,KGa随Na2CO3浓度、吸收液流量的增加而增加,随吸收温度、进气流速的升高而降低;在高浓度H2S吸收过程中液相传质阻力不能忽略。     

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.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

Persulfate oxidation of trichloroethylene with and without iron activation in porous media

In situ chemical oxidation with persulfate anion (S2O82*) is a viable technique for remediation of groundwater contaminants such as trichloroethylene (TCE). An accelerated reaction using S2O82* to destroy TCE can be achieved via chemical activation with ferrous ion to generate sulfate radicals (SO4*)(E degrees =2.6 V). The column study presented here simulates persulfate oxidation of TCE in porous media (glass beads and a sandy soil). Initial experiments were conducted to investigate persulfate transport in the absence of TCE in the column. The persulfate flushing exhibited a longer residence time and revealed a moderate persulfate interaction with soils. In TCE treatment experiments, the results indicate that the water or persulfate solution would push dissolved TCE from the column. Therefore, the effluent TCE concentration gradually increased to a maximum when about one pore volume was replaced with the flushing solution in the column. The presence of Fe2+ concentration within the column caused a quick drop in effluent TCE concentration and more TCE degradation was observed. When a TCE solution was flushing through the soil column, breakthrough of TCE concentration in the effluent was relatively slow. In contrast, when the soil column was flushed with a mixed solution of persulfate and TCE, persulfate appeared to preferentially oxidize soil oxidizable matter rather than TCE during transport. Hence, persulfate oxidation of soil organics may possibly reduce the interaction between TCE and soil (e.g., adsorption) and facilitate the transport of TCE through soil columns resulting in faster breakthrough.     

Closed-Loop Recovery of Site Remediation Gas Phase Contaminants

A novel gas phase treatment system (contaminant absorption and recovery [CAR]) for removal and subsequent recycling of gas phase VOCs from soil vapor extraction/gas stripping systems has been developed. Gas phase removal efficiencies using a packed column contactor exceed 99 percent The VOC-laden absorption fluid is subsequently vacuum-stripped of the VOCs, allowing potential condensation into liquid solvent concentrates. Partition coefficients for trichloroethylene (TCE) in triethylene glycol (TEG) ranged to ca. 5.0 mole fraction gas/mole fraction liquid, indicating a significant capacity for removal from the gas phase. Results of pilot-scale operation indicate favorable removal efficiencies and cost-effective performance in comparison to GAC or thermal destruction processes. System mass transfer coefficient predictions were done, using a variety of mathematical models and compared to experimental results. A modified Mangers and Ponten correlation was found to describe system mass transfer coefficients well. The impact of water carry-over on TCE/TEG partitioning was found to be significant. The standard change in entholpy (ΔH°) and the standard change in volume (ΔV°) values were also calculated, and predictions of temperature and pressure on system performance were evaluated.     

Locally-calibrated light transmission visualization methods to quantify nonaqueous phase liquid mass in porous media

Five locally-calibrated light transmission visualization (LTV) methods were tested to quantify nonaqueous phase liquid (NAPL) mass and mass reduction in porous media. Tetrachloroethylene (PCE) was released into a two-dimensional laboratory flow chamber packed with water-saturated sand which was then flushed with a surfactant solution (2% Tween 80) until all of the PCE had been dissolved. In all the LTV methods employed here, the water phase was dyed, rather than the more common approach of dyeing the NAPL phase, such that the light adsorption characteristics of NAPL did not change as dissolution progressed. Also, none of the methods used here required the use of external calibration chambers. The five visualization approaches evaluated included three methods developed from previously published models, a binary method, and a novel multiple wavelength method that has the advantage of not requiring any assumptions about the intra-pore interface structure between the various phases (sand/water/NAPL). The new multiple wavelength method is also expected to be applicable to any translucent porous media containing two immiscible fluids (e.g., water-air, water-NAPL). Results from the sand-water-PCE system evaluated here showed that the model that assumes wetting media of uniform pore size (Model C of Niemet and Selker, 2001) and the multiple wavelength model with no interface structure assumptions were able to accurately quantify PCE mass reduction during surfactant flushing. The average mass recoveries from these two imaging methods were greater than 95% for domain-average NAPL saturations of approximately 2.6x10(-2), and were approximately 90% during seven cycles of surfactant flushing that sequentially reduced the average NAPL saturation to 7.5x10(-4).     

Modeling the Reduction of Vapor Phase Emissions from Surface Soils Due to Soil Matrix Effects: Porosity/Tortuosity Concepts

ABSTRACT

A major route for transport of volatile organic compounds within porous media is vapor phase diffusion. The diffusion rate through a porous medium is less than that through free-air due to the decreased cross-sectional area available for gas movement and the increased path length due to pore tortuosity. Numerous empirical expressions have been published that relate the diffusion coefficient in porous media to the diffusion coefficient in free-air (unobstructed gas phase). Published measurements of relative diffusivity and air-filled porosity were combined into a database. Empirical expressions available in the literature, including the popular Millington-Quirk equation, were evaluated along with a fourth-degree polynomial expression developed by the authors to determine the best type of equation to predict relative diffusivity as a function of air-filled porosity over the domain of values for porosity ranging from 0.071 to 1 for different types of materials. Mean square deviations were used as the statistical test to compare equations. The polynomial expression developed in this project produced a significantly different effective diffusion coefficient (1.3 x 10^-6 m²/sec) compared to values of 9.2 x 10^-6 m²/sec and 3.1 x 10^-6 m²/ sec predicted by forms of the Millington-Quirk equation for a specific case.