No description available.

No description available.

No description available.

<p> Matrix diffusion from planar fractures was studied both mathematically and through physical model experiments. A conceptual model was developed based on previous work by Parker (1994) and Crank (1956). Mathematical models were developed to simulate diffusion from 2D and 3D instantaneous disk sources and a 3 D continuous disk source. The models were based on analytical solutions previously developed by Carslaw and Jaeger (1959). Analytical solution is not available for the total mass diffused into the porous matrix for a 3D continuous disk source, and it was therefore calculated through the summation of the iso-concentration lines, which were assumed to be a semi-spherical shape.</p> <p> The mathematical simulations indicated that the 2D scenario produces significantly different results from the 3D scenario, the time for mass disappearance is significantly larger for continuous sources than for instantaneous sources, the normalized concentration generally decreased over time for instantaneous sources while it increased over time for continuous sources, diffusion rates decrease significantly over time and space, and the normalized mass loss from the source zone never reaches 1 for continuous sources due to the semi-infinite integral. The simulations also showed that disappearance times increase exponentially with increasing source radii and matrix porosity, and decrease with increasing aqueous-phase NAPL solubilities.</p> <p> The observations from the physical model experiments were very close to the simulated data at z = 0, validating the 3D mathematical models for this elevation. A plot of the observed vs simulated data did not reveal any trends, indicating that the majority of the differences can be attributed to experimental error. The experimental concentrations were below the method detection limit at depths of 3 and 6 cm however, indicating that either the experiments should have been conducted over a longer time period or a more sensitive analytical method should have been employed, to enable model validation at these depths.</p> / Thesis / Master of Applied Science (MASc)

Creosote, containing many high molecular weight hydrophobic polycyclic aromatic hydrocarbons (PAH's), is present in the subsurface environment at the Oneida Tie-Yard in Oneida, Tennessee. Phytoremediation using hybrid poplar trees was chosen as the remedial technology on-site. Since monitoring began, the contaminant plume has been shrinking consistently and evidence has shown that remediation is taking place. However, remediation may be rate-limited by the desorption and dissolution kinetics of the PAH's on-site. The objectives of this research are to: (1) estimate the desorption and dissolution rates of 10 PAH's found in the subsurface and (2) estimate the amount of each PAH and total mass of contaminant that is irreversibly sorbed to the soil. Three laboratory desorption and dissolution experiments were performed using contaminated soil samples from the Oneida Tie-Yard site. The first experiment was a batch desorption equilibrium experiment, the second was a batch desorption kinetics experiment, and the third was a soil column dissolution kinetics experiment. The target compounds in this study were: naphthalene, acenaphthylene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, pyrene, chrysene, and benzo(b)fluoranthene. The resulting data for the desorption equilibrium experiment revealed that rates of equilibrium were truly not instantaneous in the systems studied. However, because approximately 76% of PAH's desorbed by the first sampling event (3 days), an equilibrium isotherm was considered appropriate. Results showed that there is a sorbed reversible concentration that readily desorbs to the aqueous phase for each PAH. Additionally, it was determined that the percent removal of sorbed PAH's decreases with increasing molecular weight. Desorption curves based on experimental data were found to exhibit linear behavior over large variations in aqueous concentration, but showed exponential behavior as concentrations approached zero. Freundlich sorption equilibrium isotherms for the 10 monitored PAH's on-site were generally found to have N coefficient values over 1, especially over large variations in solution phase concentration, indicating a non-uniform sorbent. Dissolution of resistant PAH's under field-like conditions was determined to occur over long periods of time. Dissolution rates calculated from experimental data were shown to generally decrease with increasing molecular weight. Overall, desorption and dissolution kinetics of PAH's were shown to be rate-limiting factors to remediation at the Oneida Tie-Yard. / Master of Science

In Situ Thermal Treatment (ISTT) is a candidate remediation technology for dense non-aqueous phase liquids (DNAPLs). However, the relationships between gas production, gas flow, and contaminant mass removal during ISTT are not fully understood. A laboratory study was conducted to assess the degree of mass removal, as well as the gas generation rate and the composition of the gas phase as a function of different heating times and initial DNAPL saturations. The temperature of the contaminated soil was measured continuously using a thermocouple to identify periods of heating, co-boiling and boiling. Samples were collected from the aqueous and DNAPL phase of the condensate, as well as from the source soil, at different heating times, and analyzed by gas chromatography/mass spectrometry. In addition to laboratory experiments, a mathematical model was developed to predict the co-boiling temperature and transient composition of the gas phase during heating of a uniform source. Predictions for single-component sources matched the experiments well, with a co-boiling plateau at 88°C ± 1°C for experiments with tetrachloroethene (PCE) and water. A comparison of predicted and observed boiling behaviour showed a discrepancy at the end of the co-boiling period, with earlier temperature increases occurring in the experiments. The results of this study suggest that temperature observations related to the co-boiling period during ISTT applications may not provide a clear indication of complete NAPL mass removal, and that multi-compartment modeling associated with various NAPL saturation zones is required to consider mass-transfer limitations within the heated zone. Predictions for multi-component DNAPL, containing 1,2-Dichloroethane (1,2-DCA), PCE and Chlorobenzene, showed no co-boiling plateau. CB is the least volatile component and dominates in the vapour phase at the end of the co-boiling process, and it can be used as an indicator of the end of the co-boiling stage. Two field NAPL mixtures were simulated using the screening-level analytical model to demonstrate its potential application on ISTT. The two mixtures with similar composition but different mass fractions result in distinct co-boiling temperature and mass transfer behaviour. The non-volatile component in the NAPL mixture results in larger amounts of water consumption and longer ISTT operation time. / Thesis (Master, Civil Engineering) -- Queen's University, 2013-09-30 09:26:00.857

Fractured rock formations represent a valuable source of groundwater and can be highly susceptible to contamination by dense, non-aqueous phase liquids (DNAPLs). The goal of this research is to evaluate the effectiveness of three accepted remediation technologies for addressing DNAPL contamination in fractured rock environments. The technologies under investigation in this study are chemical oxidation, bioremediation, and surfactant flushing. Numerical simulations were employed to examine the performance of each of these technologies at the field scale. The numerical model DNAPL3D-RX, a finite difference multiphase flow-dissolution-aqueous transport code that incorporates RT3D for multiple species reactions, was modified to simulate fractured rock environments. A gridding routine was developed to allow the model to accurately capture DNAPL migration in fractures and aqueous phase diffusion gradients in the matrix while retaining overall model efficiency. Reaction kinetics code subroutines were developed for each technology so as to ensure the key processes were accounted for in the simulations. The three remedial approaches were systematically evaluated via simulations in two-dimensional domains characterized by heterogeneous orthogonal fracture networks parameterized to be representative of sandstone, granite, and shale. Each simulation included a DNAPL release at the water table, redistribution to pools and residual, followed by 20 years of ‘ageing’ under ambient gradient conditions. Suites of simulations for each technology examined a variety of operational issues including the influence of DNAPL type and remedial fluid injection protocol. Performance metrics included changes in mass flux exiting, mass destruction in the matrix versus the fractures, and percentage of injected remedial fluid interacting with the target contaminant. The effectiveness of the three remediation technologies covered a wide range; the mass of contaminants destroyed were found to range from 15% to 99.5% of the initial mass present. Effectiveness of each technology was found to depend on a variety of critical factors particular to each approach. For example, in-situ chemical oxidation was found to be limited by the organic material present in the matrix of the rocks, while the efficiency of enhanced bioremediation was found to be related to factors such as the location of indigenous bacteria present in the domain and rate of bioremediation. In the chemical oxidation study, the efficiency of oxidant consumption was observed to be poor across the suite of scenarios, with greater than 90% of the injected permanganate consumed by natural oxidant demand. This study further revealed that the same factors that contributed to forward diffusion of contaminants prior to treatment are critical to this remediation method as they can determine the extent of contaminant destruction during the injection period. Bioremediation in fractured rock was demonstrated to produce relatively good results under robust first-order decay rates and active microorganisms throughout the fractures and matrix. It was demonstrated that under ideal conditions, of the total initial mass present, up to 3/4 could be reduced to ethene, indicating bioremediation may be a promising treatment approach due to the effective penetration of electron donor into the matrix during the treatment period and the ongoing treatment that occurs after injection ceases. However, when indigenous bacteria was assumed to exist only within the fractured walls of sandstone, it was found that under the same conditions, the rate of dechlorination was 200 times less than the Base Case. Since the majority of the mass resided in the matrix, lack of bioremediation in the matrix significantly reduced the effectiveness of treatment. Surfactant treatment with Tween-80 was proven to be a relatively effective technique in enhanced solubilisation of DNAPL from the fractures within the domain. However, by comparing the aqueous and sorbed mass at the start and end of the Treatment stage, it is revealed that surfactant treatment is not efficient in removing these masses that reside within the matrix. Furthermore, DNAPLs identified in dead end vertical fractures were found to remain in the domain by the end of the simulations across all scenarios studied; indicating that the injected surfactant experiences difficulty in accessing DNAPLs entrapped in dead end fractures. Altogether, the results underscore the challenge of restoring fractured rock aquifers due to the field scale limitations on sufficient contact between remedial fluids and in situ contaminants in all but the most ideal circumstances.

Degassing and in situ formation of a mobile gas phase takes place when an aqueous phase equilibrated with a gas at a pressure higher than the subsurface pressure is injected in water-saturated porous media. This process, which has been termed supersaturated water injection (SWI), is a novel and hitherto unexplored means of introducing a gas phase into the subsurface. Herein is a first macroscopic account of the SWI process on the basis of continuum scale simulations and column experiments with CO2 as the dissolved gas. A published empirical mass transfer correlation (Nambi and Powers, Water Resour Res, 2003) is found to adequately describe the non-equilibrium transfer of CO2 between the aqueous and gas phases. Remarkably, the dynamics of gas-water two-phase flow, observed in a series of SWI experiments in homogeneous columns packed with silica sand or glass beads, are accurately predicted by traditional two-phase flow theory which allows the corresponding gas phase relative permeability to be determined. A key consequence of the finding, that the displacement of the aqueous phase by gas is compact at the macroscopic scale, is consistent with pore scale simulations of repeated mobilization, fragmentation and coalescence of large gas clusters (i.e., large ganglion dynamics) driven entirely by mass transfer. The significance of this finding for the efficient delivery of a gas phase below the water table in relation to the alternative process of in-situ air sparging and the potential advantages of SWI are discussed. SWI has been shown to mobilize a previously immobile oil phase in the subsurface of 3-phase systems (oil, water and gas). A macroscopic account of the SWI process is given on the basis of continuum-scale simulations and column experiments using CO2 as the dissolved gas and kerosene as the trapped oil phase. Experimental observations show that the presence of oil ganglia in the subsurface alters gas phase mobility from 2-phase predictions. A corresponding 3-phase gas relative permeability function is determined, whereas a published 3-phase relative permeability correlation (Stone, Journal of Cana Petro Tech, 1973) is found to be inadequate for describing oil phase flow during SWI. A function to predict oil phase relative permeability is developed for use during SWI at high aqueous phase saturations with a disconnected oil phase and quasi-disconnected gas phase. Remarkably, the dynamics of gas-water-oil 3-phase flow, observed in a series of SWI experiments in homogeneous columns packed with silica sand or glass beads, are accurately predicted by traditional continuum-scale flow theory. The developed relative permeability function is compared to Stone’s Method and shown to approximate it in all regions while accurately describing oil flow during SWI. A published validation of Stone’s Method (Fayers and Matthews, Soc of Petro Eng Journal, 1984) is cited to validate this approximation of Stone’s Method.

Dimethyl Sulfide (CH3SCH3, DMS) is the most important natural sulfur compound emitted from the ocean and its oxidation in the atmosphere has been proposed to play an important role in climate modification because some products from DMS oxidation become non-volatile and could participate in particle formation and growth processes. Although it has been demonstrated that aqueous phase reactions are potentially important for understanding DMS oxidation, the kinetics database for aqueous phase transformations is rather limited. In this work, a laser flash photolysis (LFP) ??ng path UV-visible absorption (LPA) technique was employed to investigate the kinetics of the aqueous phase reactions of four organic sulfur compounds produced from DMS oxidation, i.e., dimethylsulfoxide (DMSO), dimethyl-sulfone (DMSO2), methanesulfinate (MSI) and methanesulfonate (MS), with four important aqueous phase radicals, OH, SO4 and #8722;, Cl and Cl2 and #8722;. The temperature-dependent kinetics of the OH and SO4 and #8722; reactions with DMSO, DMSO2 and MS were studied for the first time. OH is found to be the most reactive, while Cl2 and #8722; is the least reactive toward all the sulfur species studied. The less oxidized DMSO and MSI are found to be more reactive than the more oxidized DMSO2 and MS for each radical. The kinetic data have been employed in a Trajectory Ensemble Model to simulate DMS oxidation in the marine atmosphere as a means of assessing the contribution of aqueous phase reactions to the growth of particulate matter. For the first time, oxidation of organic sulfur compounds by SO4 and #8722;, Cl and Cl2 and #8722; are included in the model to simulate DMS chemistry. Our simulations suggest that aqueous phase reactions contribute >97% of MS and ~90% of NSS (Non-Seasalt Sulfate) production, and aqueous phase reactions of the organic sulfur compounds contribute 30% of total particle mass growth. When our kinetic data for the MS + OH reaction were used in the model, it was found that MS + OH could consume ~20% of MS and produce ~8% of NSS, within 3 days under typical marine atmospheric conditions.

Groundwater at many existing and former industrial sites and disposal areas is contaminated by halogenated organic compounds that were released into the environment. The chlorinated solvent trichloroethylene (TCE) is one of the most ubiquitous of these compounds. In situ chemical oxidation (ISCO) has been successfully used for the removal of TCE. The objective of this study was to apply the ISCO technology to remediate TCE-contaminated groundwater. In this study, potassium permanganate (KMnO4) was used as the oxidant during the ISCO process. The study consisted bench-scale and pilot-scale experiments. In the laboratory experiments, the major controlling factors included oxidant concentrations, effects of soil oxidant demand (SOD) on oxidation efficiency, and addition of dibasic sodium phosphate on the inhibition of production of manganese dioxide (MnO2). Results show that higher molar ratios of KMnO4 to TCE corresponded with higher TCE oxidation rate under the same initial TCE concentration condition. Moreover, higher TCE concentration corresponded with higher TCE oxidation rate under the same molar ratios of KMnO4 to TCE condition. Results reveal that KMnO4 is a more stable and dispersive oxidant, which is able to disperse into the soil materials and react with organic contaminants effectively. Significant amount of MnO2 production can be effectively inhibited with the addition of Na2HPO4. Results show that the increase in the first-order decay rate was observed when the oxidant concentration was increased, and the half-life was approximately 24.3 to 251 min. However, the opposite situation was observed when the second-order decay rate was used to describe the reaction. Results from the column experiment show that the breakthrough volumes were approximately 50.4 to 5.06 pore volume (PV). Injection of KMnO4 would cause the decrease in TCE concentration through oxidation. Results also indicate that the addition of Na2HPO4 would not inhibit the TCE removal rate. In the second part of this study, a TCE-contaminated site was selected for the conduction of pilot-scale study. A total of eight remediation wells were installed for this pilot-scale study. The initial TCE concentrations of the eight wells were as follows: C1 = 0.59 mg/L, C1-E = 0.64 mg/L, C1-W = 0.61 mg/L, EW-1 = 0.65 mg/L, EW-1E = 0.62 mg/L, EW-1W = 0.57 mg/L, C2 = 0.62 mg/L, C3 = 0.35 mg/L. C1, EW-1, C2, and C3 were located along the groundwater flow direction from the upgradient (C1) to the downgradient location (C3), and the distance between each well was 3 m. C1-E and C1-W were located in lateral to C1 with a distance of 3 m to C1. EW-1E and EW-1W were in lateral to EW-1 with a distance of 3 m to EW-1. In the first test, 2,700 L of KMnO4 solution was injected into each of the three injection wells (C1, C1-E, and C1-W) with concentration of 5,000 mg/L. Three injections were performed with an interval of 6 hr between each injection. After injection, the TCE concentrations in those three wells dropped down to below detection limit (<0.0025 mg/L). However, no significant variations in TCE concentrations were observed in other wells. In the second test, 2,700 L of KMnO4 solution was injected into injection well (EW-1) with concentration of 5,000 mg/L. Six injections were performed with an interval of 6 hr between each injection. After injection, the TCE concentrations in the injection well dropped down to below detection limit (<0.0025 mg/L). TCE concentrations in (C1, C1-E, C1-W, EW-1E, EW-1W, C2, and C3) dropped to 0.35-0.49 mg/L. After injection, no significant temperature and pH variation was observed. However, increase in conductivity and oxidation-reduction potential (ORP) was observed. This indicates that the KMnO4 oxidation process is a potential method for TCE-contaminate site remediation. The groundwater conductivity increased from 500 £gS/cm to 1,000 £gS/cm, and ORP increased from 200 to 600 mv. Increase in KMnO4, MnO2, and total Mn was also observed in wells. Results from the slug tests show that the hydraulic conductivity remained in the range from 10-4 to 10-5 m/sec before and after the KMnO4 injection.