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

Geothermal Play Fairway Analysis, Part 2: GIS Methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development

Geothermal Play Fairway Analysis, Part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales

Implications of the Drift Scale Heater Test at Yucca Mountain for Epithermal Mineralization

An 8-year long, drift scale heater test (DST) is currently underway at the underground Exploratory Studies Facility at Yucca Mountain in Nevada. The host rock for the DST is a highly fractured, welded tuff. The rock has {approx}10% matrix porosity 90% filled with water. After a little more than two years of heating, the temperature at the drift wall reached {approx}200 C and has been maintained at that temperature for the past {approx}1.5 years. Gas and water (both vapor and liquid) have been collected from monitoring boreholes since the test began. The CO{sub 2} concentration of the gas and the isotopic compositions of the water and CO{sub 2} are measured. These data are used to constrain numerical models of coupled thermal, hydrological, and chemical processes occurring in the system. Despite obvious differences from epithermal systems (e.g., the DST is being conducted in an unsaturated system), the trends observed in the isotopic compositions of the water and CO{sub 2} have interesting implications for natural systems. In areas below boiling, the isotope ratios of the water are near that of the ambient pore water ({delta}{sup 18}O about -12{per_thousand}). Where significant amounts of vapor condensate occur (above the boiling front above the drift and in fracture zones to the sides of the drift), the {delta}{sup 18}O values of the water are lower than the pore water, reflecting addition of low-{delta}{sup 18}O steam condensate. Conversely, in boiling zones, the {delta}{sup 18}O values of the water become progressively higher, representing Rayleigh fractionation of the pore water as it is vaporized. As the temperature approaches boiling, the gas phase becomes dominated by water vapor. The remainder of the gas phase consists of air with elevated CO{sub 2} (up to 15%). The source of the CO, is primarily dissolved inorganic carbon (DIC) in the pore water. As the temperature increases, the {delta}{sup 13}C values of the CO{sub 2} shift from approximate equilibrium with the pore water DIC (-15{per_thousand}) to much higher values (>0{per_thousand}). Dissolution of calcite in fractures is also a significant source of CO{sub 2} in regions with drainage of vapor condensate. Isotopic data from several Mexican epithermal vein systems will be discussed in light of these findings

Modeling of coupled heat transfer and reactive transport processes in porous media: Application to seepage studies at Yucca Mountain, Nevad a

When hot radioactive waste is placed in subsurface tunnels, a series of complex changes occurs in the surrounding medium. The water in the pore space of the medium undergoes vaporization and boiling. Subsequently, vapor migrates out of the matrix pore space, moving away from the tunnel through the permeable fracture network. This migration is propelled by buoyancy, by the increased vapor pressure caused by heating and boiling, and through local convection. In cooler regions, the vapor condenses on fracture walls, where it drains through the fracture network. Slow imbibition of water thereafter leads to gradual rewetting of the rock matrix. These thermal and hydrological processes also bring about chemical changes in the medium. Amorphous silica precipitates from boiling and evaporation, and calcite from heating and CO2 volatilization. The precipitation of amorphous silica, and to a much lesser extent calcite, results in long-term permeability reduction. Evaporative concentration also results in the precipitation of gypsum (or anhydrite), halite, fluorite and other salts. These evaporative minerals eventually redissolve after the boiling period is over, however, their precipitation results in a significant temporary decrease in permeability. Reduction of permeability is also associated with changes in fracture capillary characteristics. In short, the coupled thermal-hydrological-chemical (THC) processes dynamically alter the hydrological properties of the rock. A model based on the TOUGHREACT reactive transport software is presented here to investigate the impact of THC processes on flow near an emplacement tunnel at Yucca Mountain, Nevada. We show how transient changes in hydrological properties caused by THC processes often lead to local flow channeling and saturation increases above the tunnel. For models that include only permeability changes to fractures, such local flow channeling may lead to seepage relative to models where THC effects are ignored. However, coupled THC seepage models that include both permeability and capillary changes to fractures may not show this additional seepage