Temperature Effects on seepage Fluid Compositions at Yucca Mountain

This project investigated the effect of two repository operating temperature modes on coupled thermal, hydrological, and chemical processes around potential nuclear waste-emplacement tunnels (drifts) at Yucca Mountain, Nevada. The main objective of this study was to evaluate the composition of fluids (water and gas) that could enter the drifts, because these data directly relate to the performance of waste canisters and other in-drift engineered systems over the life of the potential repository. Multicomponent reactive transport simulations were performed using TOUGHREACT, initially written by T. Xu and K. Pruess at LBNL and modified here to handle high-temperature and boiling environments. Two repository operating temperature modes were investigated: (1) a ''high-temperature'' mode, which considered a short preclosure ventilation period (50 years) and gave rise to above-boiling temperatures in rocks around the drift for hundreds of years, and (2) a ''low-temperature'' mode with a smaller heat load and longer preclosure ventilation (300 years), yielding temperatures at the surface of the waste package below 85 C (a design threshold) and thus below boiling conditions. Simulations under ambient conditions (no heat load) were also conducted to serve as a baseline for comparing results of thermal-loading simulations

TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic Media, V1.2.1

Coupled modeling of subsurface multiphase fluid and heat flow, solute transport, and chemical reactions can be applied to many geologic systems and environmental problems, including geothermal systems, diagenetic and weathering processes, subsurface waste disposal, acid mine drainage remediation, contaminant transport, and groundwater quality. TOUGHREACT has been developed as a comprehensive non-isothermal multi-component reactive fluid flow and geochemical transport simulator to investigate these and other problems. A number of subsurface thermo-physical-chemical processes are considered under various thermohydrological and geochemical conditions of pressure, temperature, water saturation, and ionic strength. TOUGHREACT can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity. The code can accommodate any number of chemical species present in liquid, gas and solid phases. A variety of equilibrium chemical reactions are considered, such as aqueous complexation, gas dissolution/exsolution, and cation exchange. Mineral dissolution/precipitation can take place subject to either local equilibrium or kinetic controls, with coupling to changes in porosity and permeability and capillary pressure in unsaturated systems. Chemical components can also be treated by linear adsorption and radioactive decay. The first version of the non-isothermal reactive geochemical transport code TOUGHREACT was developed (Xu and Pruess, 1998) by introducing reactive geochemistry into the framework of the existing multi-phase fluid and heat flow code TOUGH2 (Pruess, 1991). TOUGHREACT was further enhanced with the addition of (1) treatment of mineral-water-gas reactive-transport under boiling conditions, (2) an improved HKF activity model for aqueous species, (3) gas species diffusion coefficients calculated as a function of pressure, temperature, and molecular properties, (4) mineral reactive surface area formulations for fractured and porous media, and (5) porosity, permeability, and capillary pressure changes owing to mineral precipitation/dissolution (Sonnenthal et al., 1998, 2000, 2001; Spycher et al., 2003a). Subsequently, TOUGH2 V2 was released with additional EOS modules and features (Pruess et al., 1999). The present version of TOUGHREACT includes all of the previous extensions to the original version, along with the replacement of the original TOUGH2 (Pruess, 1991) by TOUGH2 V2 (Pruess et al., 1999). TOUGHREACT has been applied to a wide variety of problems, some of which are included as examples, such as: (1) Supergene copper enrichment (Xu et al., 2001); (2) Mineral alteration in hydrothermal systems (Xu and Pruess, 2001a; Xu et al., 2004b; Dobson et al., 2004); (3) Mineral trapping for CO{sub 2} disposal in deep saline aquifers (Xu et al., 2003b and 2004a); (4) Coupled thermal, hydrological, and chemical processes in boiling unsaturated tuff for the proposed nuclear waste emplacement site at Yucca Mountain, Nevada (Sonnenthal et al., 1998, 2001; Sonnenthal and Spycher, 2000; Spycher et al., 2003a, b; Xu et al., 2001); (5) Modeling of mineral precipitation/dissolution in plug-flow and fracture-flow experiments under boiling conditions (Dobson et al., 2003); (6) Calcite precipitation in the vadose zone as a function of net infiltration (Xu et al., 2003); and (7) Stable isotope fractionation in unsaturated zone pore water and vapor (Singleton et al., 2004). The TOUGHREACT program makes use of 'self-documenting' features. It is distributed with a number of input data files for sample problems. Besides providing benchmarks for proper code installation, these can serve as a self-teaching tutorial in the use of TOUGHREACT, and they provide templates to help jump-start new applications. The fluid and heat flow part of TOUGHREACT is derived from TOUGH2 V2, so in addition to the current manual, users must have the manual of the TOUGH2 V2 (Pruess et al., 1999). The present version of TOUGHREACT provides the following TOUGH2 fluid property or 'EOS' (equation-of-state) modules: (1) EOS1 for water, or two waters with typical applications to hydrothermal problems, (2) EOS2 for multiphase mixtures of water and CO{sub 2} also with typical applications to hydrothermal problems, (3) EOS3 for multiphase mixtures of water and air with typical applications to vadose zone and nuclear waste disposal problems, (4) EOS4 that has the same capabilities as EOS3 but with vapor pressure lowering effects due to capillary pressure, (5) EOS9 for single phase water (Richards equation) with typical applications to ambient temperature and pressure reactive geochemical transport problems, and (6) ECO2N for multiphase mixtures of water, CO{sub 2} and NaCl with typical applications to CO{sub 2} disposal in deep brine aquifers

Modeling coupled thermal-hydrological-chemical processes in theunsaturated fractured rock of Yucca Mountain, Nevada: Heterogeneity andseepage

An understanding of processes affecting seepage intoemplacement tunnels is needed for correctly predicting the performance ofunderground radioactive waste repositories. It has been previouslyestimated that the capillary and vaporization barriers in the unsaturatedfractured rock of Yucca Mountain are enough to prevent seepage underpresent day infiltration conditions. It has also been thought that asubstantially elevated infiltration flux will be required to causeseepage after the thermal period is over. While coupledthermal-hydrological-chemical (THC) changes in Yucca Mountain host rockdue to repository heating has been previously investigated, those THCmodels did not incorporate elements of the seepage model. In this paper,we combine the THC processes in unsaturated fractured rock with theprocesses affecting seepage. We observe that the THC processes alter thehydrological properties of the fractured rock through mineralprecipitation and dissolution. We show that such alteration in thehydrological properties of the rock often leads to local flow channeling.We conclude that such local flow channeling may result in seepage undercertain conditions, even with nonelevated infiltrationfluxes

Advances in Hydrogeochemical Indicators for the Discovery of New Geothermal Resources in the Great Basin, USA

This report summarizes the results of Phase I work for a go/no go decision on Phase II funding. In the first objective, we assessed the extent to which fluid-mineral equilibria controlled deep water compositions in geothermal systems across the Great Basin. Six systems were evaluated: Beowawe; Desert Peak; Dixie Valley; Mammoth; Raft River; Roosevelt. These represent a geographic spread of geothermal resources, in different geological settings and with a wide range of fluid compositions. The results were used for calibration/reformulation of chemical geothermometers that reflect the reservoir temperatures in producing reservoirs. In the second objective, we developed a reactive -transport model of the Desert Peak hydrothermal system to evaluate the processes that affect reservoir fluid geochemistry and its effect on solute geothermometry. This included testing geothermometry on “reacted” thermal water originating from different lithologies and from near-surface locations where the temperature is known from the simulation. The integrated multi-component geothermometer (GeoT, relying on computed mineral saturation indices) was tested against the model results and also on the systems studied in the first objective

Temporal changes in noble gas compositions within the Aidlinsector ofThe Geysers geothermal system

The use of nonreactive isotopic tracers coupled to a full thermal-hydrological reservoir simulation allows for an improved method of investigating how reservoir fluids contained within matrix and fractures contribute over time to fluids produced from geothermal systems. A combined field and modeling study has been initiated to evaluate the effects of injection, production, and fracture-matrix interaction on produced noble gas contents and isotopic ratios. Gas samples collected periodically from the Aidlin steam field at The Geysers, California, between 1997 and 2006 have been analyzed for their noble gas compositions, and reveal systematic shifts in abundance and isotopic ratios over time. Because of the low concentrations of helium dissolved in the injection waters, the injectate itself has little impact on the helium isotopic composition of the reservoir fluids over time. However, the injection process may lead to fracturing of reservoir rocks and an increase in diffusion-controlled variations in noble gas compositions, related to gases derived from fluids within the rock matrix

Microbial U isotope fractionation depends on U(VI) reduction rate

U isotope fractionation may serve as an accurate proxy for U(VI) reduction in both modern and ancient environments, if the systematic controls on the magnitude of fractionation (ε) are known. We model the effect of U(VI) reduction kinetics on U isotopic fractionation during U(VI) reduction by a novel Shewanella isolate, Shewanella sp. (NR), in batch incubations. The measured ε values range from 0.96 ± 0.16 to 0.36 ± 0.07‰ and are strongly dependent on the U(VI) reduction rate. The ε decreases with increasing reduction rate constants normalized by cell density and initial U(VI). Reactive transport simulations suggest that the rate dependence of ε is due to a two-step process, where diffusive transport of U(VI) from the bulk solution across a boundary layer is followed by enzymatic reduction. Our results imply that the spatial decoupling of bulk U(VI) solution and enzymatic reduction should be taken into account for interpreting U isotope data from the environment

Modeling Coupled Processes in Clay Formations for Radioactive Waste Disposal

As a result of the termination of the Yucca Mountain Project, the United States Department of Energy (DOE) has started to explore various alternative avenues for the disposition of used nuclear fuel and nuclear waste. The overall scope of the investigation includes temporary storage, transportation issues, permanent disposal, various nuclear fuel types, processing alternatives, and resulting waste streams. Although geologic disposal is not the only alternative, it is still the leading candidate for permanent disposal. The realm of geologic disposal also offers a range of geologic environments that may be considered, among those clay shale formations. Figure 1-1 presents the distribution of clay/shale formations within the USA. Clay rock/shale has been considered as potential host rock for geological disposal of high-level nuclear waste throughout the world, because of its low permeability, low diffusion coefficient, high retention capacity for radionuclides, and capability to self-seal fractures induced by tunnel excavation. For example, Callovo-Oxfordian argillites at the Bure site, France (Fouche et al., 2004), Toarcian argillites at the Tournemire site, France (Patriarche et al., 2004), Opalinus clay at the Mont Terri site, Switzerland (Meier et al., 2000), and Boom clay at Mol site, Belgium (Barnichon et al., 2005) have all been under intensive scientific investigations (at both field and laboratory scales) for understanding a variety of rock properties and their relations with flow and transport processes associated with geological disposal of nuclear waste. Clay/shale formations may be generally classified as indurated and plastic clays (Tsang et al., 2005). The latter (including Boom clay) is a softer material without high cohesion; its deformation is dominantly plastic. For both clay rocks, coupled thermal, hydrological, mechanical and chemical (THMC) processes are expected to have a significant impact on the long-term safety of a clay repository. For example, the excavation-damaged zone (EDZ) near repository tunnels can modify local permeability (resulting from induced fractures), potentially leading to less confinement capability (Tsang et al., 2005). Because of clay's swelling and shrinkage behavior (depending on whether the clay is in imbibition or drainage processes), fracture properties in the EDZ are quite dynamic and evolve over time as hydromechanical conditions change. To understand and model the coupled processes and their impact on repository performance is critical for the defensible performance assessment of a clay repository. Within the Natural Barrier System (NBS) group of the Used Fuel Disposition (UFD) Campaign at DOE's Office of Nuclear Energy, LBNL's research activities have focused on understanding and modeling such coupled processes. LBNL provided a report in this April on literature survey of studies on coupled processes in clay repositories and identification of technical issues and knowledge gaps (Tsang et al., 2010). This report will document other LBNL research activities within the natural system work package, including the development of constitutive relationships for elastic deformation of clay rock (Section 2), a THM modeling study (Section 3) and a THC modeling study (Section 4). The purpose of the THM and THC modeling studies is to demonstrate the current modeling capabilities in dealing with coupled processes in a potential clay repository. In Section 5, we discuss potential future R&D work based on the identified knowledge gaps. The linkage between these activities and related FEPs is presented in Section 6

Reactive transport model of sulfur cycling as impacted by perchlorate and nitrate treatments

Microbial souring in oil reservoirs produces toxic, corrosive hydrogen sulfide through microbial sulfate reduction, often accompanying (sea)­water flooding during secondary oil recovery. With data from column experiments as constraints, we developed the first reactive-transport model of a new candidate inhibitor, perchlorate, and compared it with the commonly used inhibitor, nitrate. Our model provided a good fit to the data, which suggest that perchlorate is more effective than nitrate on a per mole of inhibitor basis. Critically, we used our model to gain insight into the underlying competing mechanisms controlling the action of each inhibitor. This analysis suggested that competition by heterotrophic perchlorate reducers and direct inhibition by nitrite produced from heterotrophic nitrate reduction were the most important mechanisms for the perchlorate and nitrate treatments, respectively, in the modeled column experiments. This work demonstrates modeling to be a powerful tool for increasing and testing our understanding of reservoir-souring generation, prevention, and remediation processes, allowing us to incorporate insights derived from laboratory experiments into a framework that can potentially be used to assess risk and design optimal treatment schemes