Patrick F. Dobson

Porosity, permeability and fluid flow in the Yellowstone Geothermal System, Wyoming

Cores from two of 13 U.S. Geological Survey research holes at Yellowstone National Park (Y-5 and Y-8) were evaluated to characterize lithology, texture, alteration, and the degree and nature of fracturing and veining. Porosity and matrix permeability measurements and petrographic examination of the cores were used to evaluate the effects of lithology and hydrothermal alteration on porosity and permeability. The intervals studied in these two core holes span the conductive zone and the upper portion of the convective geothermal reservoir. Variations in porosity and matrix permeability observed in the Y-5 and Y-8 cores are primarily controlled by lithology. Y-8 intersects three distinct lithologies: volcaniclastic sandstone, perlitic rhyolitic lava, and non-welded pumiceous ash-flow tuff. The sandstone typically has high permeability and porosity, and the tuff has very high porosity and moderate permeability, while the perlitic lava has very low porosity and is essentially impermeable. Hydrothermal self-sealing appears to have generated localized permeability barriers within the reservoir. Changes in pressure and temperature in Y-8 correspond to a zone of silicification in the volcaniclastic sandstone just above the contact with the perlitic rhyolite; this silicification has significantly reduced porosity and permeability. In rocks with inherently low matrix permeability (such as densely welded ash-flow tuff), fluid flow is controlled by the fracture network. The Y-5 core hole penetrates a thick intracaldera section of the 0.6-Ma Lava Creek ash-flow tuff. In this core, the degree of welding appears to be responsible for most of the variations in porosity, matrix permeability, and the frequency of fractures and veins. Fractures are most abundant within the more densely welded sections of the tuff. However, the most prominent zones of fracturing and mineralization are associated with hydrothermal breccias within densely welded portions of the tuff. These breccia zones represent transient conduits of high fluid flow that formed by the explosive release of overpressure in the underlying geothermal reservoir and that were subsequently sealed by supersaturated geothermal fluids. In addition to this fracture sealing, hydrothermal alteration at Yellowstone appears generally to reduce matrix permeability and focus flow along fractures, where multiple pulses of fluid flow and self-sealing have occurred.

Experimental and numerical simulation of dissolution and precipitation: implications for fracture sealing at Yucca Mountain, Nevada

Plugging of flow paths caused by mineral precipitation in fractures above the potential repository at Yucca Mountain, Nevada could reduce the probability of water seeping into the repository. As part of an ongoing effort to evaluate thermal-hydrological-chemical (THC) effects on flow in fractured media, we performed a laboratory experiment and numerical simulations to investigate mineral dissolution and precipitation under anticipated temperature and pressure conditions in the repository. To replicate mineral dissolution by vapor condensate in fractured tuff, water was flowed through crushed Yucca Mountain tuff at 94 degrees C. The resulting steady-state fluid composition had a total dissolved solids content of about 140 mg/l; silica was the dominant dissolved constituent. A portion of the steady-state mineralized water was flowed into a vertically oriented planar fracture in a block of welded Topopah Spring Tuff that was maintained at 80 degrees C at the top and 130 degrees C at the bottom. The fracture began to seal with amorphous silica within 5 days.A 1-D plug-flow numerical model was used to simulate mineral dissolution, and a similar model was developed to simulate the flow of mineralized water through a planar fracture, where boiling conditions led to mineral precipitation. Predicted concentrations of the major dissolved constituents for the tuff dissolution were within a factor of 2 of the measured average steady-state compositions. The mineral precipitation simulations predicted the precipitation of amorphous silica at the base of the boiling front, leading to a greater than 50-fold decrease in fracture permeability in 5 days, consistent with the laboratory experiment.These results help validate the use of a numerical model to simulate THC processes at Yucca Mountain. The experiment and simulations indicated that boiling and concomitant precipitation of amorphous silica could cause significant reductions in fracture porosity and permeability on a local scale. However, differences in fluid flow rates and thermal gradients between the experimental setup and anticipated conditions at Yucca Mountain need to be factored into scaling the results of the dissolution/precipitation experiments and associated simulations to THC models for the potential Yucca Mountain repository.

Seismic and aseismic deformations and impact on reservoir permeability: The case of EGS stimulation at The Geysers, California, USA

Author(s): Jeanne, P; Rutqvist, J; Rinaldi, AP; Dobson, PF; Walters, M; Hartline, C; Garcia, J | Abstract: ©2015. American Geophysical Union. All Rights Reserved. In this paper, we use the Seismicity-Based Reservoir Characterization approach to study the spatiotemporal dynamics of an injection-induced microseismic cloud, monitored during the stimulation of an enhanced geothermal system, and associated with the Northwest Geysers Enhanced Geothermal System (EGS) Demonstration project (California). We identified the development of a seismically quiet domain around the injection well surrounded by a seismically active domain. Then we compare these observations with the results of 3-D Thermo-Hydro-Mechanical simulations of the EGS, which accounts for changes in permeability as a function of the effective normal stress and the plastic strain. The results of our modeling show that (1) the aseismic domain is caused by both the presence of the injected cold water and by thermal processes. These thermal processes cause a cooling-stress reduction, which prevent shear reactivation and favors fracture opening by reducing effective normal stress and locally increasing the permeability. This process is accompanied by aseismic plastic shear strain. (2) In the seismic domain, microseismicity is caused by the reactivation of the preexisting fractures, resulting from an increase in injection-induced pore pressure. Our modeling indicates that in this domain, permeability evolves according to the effective normal stress acting on the shear zones, whereas shearing of preexisting fractures may have a low impact on permeability. We attribute this lack of permeability gain to the fact that the initial permeabilities of these preexisting fractures are already high (up to 2 orders of magnitude higher than the host rock) and may already be fully dilated by past tectonic straining.

Geomechanical simulation of the stress tensor rotation caused by injection of cold water in a deep geothermal reservoir

We present a three-dimensional thermohydromechanical numerical study of the evolution and distribution of the stress tensor within the northwest part of The Geysers geothermal reservoir (in California), including a detailed study of the region around one injection well from 2003 to 2012. Initially, after imposing a normal faulting stress regime, we calculated local changes in the stress regime around injection wells. Our results were compared with previously published studies in which the stress state was inferred from inverting the focal plane mechanism of seismic events. Our main finding is that changes in stress tensor orientation are caused by injection-induced progressive cooling of the reservoir, as well as by the seasonal variations in injection rate. Because of the gravity flow and cooling around a liquid zone formed by the injection, the vertical stress reduction is larger and propagates far below the injection well. At the same time, the horizontal stress increases, mostly because of stress redistribution below and above the cooling area. These two phenomena cause the rotation of the stress tensor and the appearance of a strike-slip regime above, inside, and below the cooling area. The cooling and the associated rotation of the stress regime can play a significant role in the observed long-term deepening of the microseismicity below active injection wells.

Petrology and Geochemistry of Boninite Series Volcanic Rocks,Chichi-jima, Bonin Islands, Japan

An Eocene submarine boninite-series volcanic center is exposed on the island of Chichi-jima, Bonin Islands, Japan. Five rock types—boninite, bronzite andesite, dacite, quartz dacite, and rhyolite—were distinguished within the boninite volcanic sequence on the basis of petrography and geochemistry. Boninites contain high Mg, Ni, and Cr contents indicative of primitive melts, but have high Si contents relative to other mantle-derived magmas. All boninite-series rocks contain very low incompatible-element concentrations, and concentrations of high-field-strength elements in primitive boninites are less than half of those found in depleted mid-ocean ridge basalts. Abundances of large-ion lithophile elements are relatively high in boninite-series rocks, similar to the enrichments observed in many island-arc lavas. Trends for both major- and trace-element data suggest that the more evolved members of the boninite magma series were derived primarily through high-level fractional crystallization of boninite. Textural features, such as resorption and glomeroporphyrocrysts, and reverse chemical zonations suggest that magma mixing contributed to the development of the quartz dacites.

Uranium-Series Constraints on Radionuclide Transport and Groundwater Flow at the Nopal I Uranium Deposit, Sierra Peña Blanca, Mexico

Uranium-series data for groundwater samples from the Nopal I uranium ore deposit were obtained to place constraints on radionuclide transport and hydrologic processes for a nuclear waste repository located in fractured, unsaturated volcanic tuff. Decreasing uranium concentrations for wells drilled in 2003 are consistent with a simple physical mixing model that indicates that groundwater velocities are low ( approximately 10 m/y). Uranium isotopic constraints, well productivities, and radon systematics also suggest limited groundwater mixing and slow flow in the saturated zone. Uranium isotopic systematics for seepage water collected in the mine adit show a spatial dependence which is consistent with longer water-rock interaction times and higher uranium dissolution inputs at the front adit where the deposit is located. Uranium-series disequilibria measurements for mostly unsaturated zone samples indicate that (230)Th/(238)U activity ratios range from 0.005 to 0.48 and (226)Ra/(238)U activity ratios range from 0.006 to 113. (239)Pu/(238)U mass ratios for the saturated zone are <2 x 10(-14), and Pu mobility in the saturated zone is >1000 times lower than the U mobility. Saturated zone mobility decreases in the order (238)U approximately (226)Ra > (230)Th approximately (239)Pu. Radium and thorium appear to have higher mobility in the unsaturated zone based on U-series data from fractures and seepage water near the deposit.

Stratigraphy of the PB-1 well, nopal I Uranium deposit, Sierra Peña Blanca, Chihuahua, Mexico

The Nopal I site in the Peña Blanca uranium district has a number of geologic and hydrologic similarities to the proposed high-level radioactive waste repository at Yucca Mountain, making it a useful analogue to evaluate process models for radionuclide transport. The PB-1 well was drilled in 2003 at the Nopal I uranium deposit as part of a DOE-sponsored natural analogue study to constrain processes affecting radionuclide transport. The well penetrates through the Tertiary volcanic section down to Cretaceous limestone and intersects the regional aquifer system. The well, drilled along the margin of the Nopal I ore body, was continuously cored to a depth of 250 m, thus providing an opportunity to document the local stratigraphy. Detailed observations of these units were afforded through petrographic description and rockproperty measurements of the core, together with geophysical logs of the well. The uppermost unit encountered in the PB-1 well is the Nopal Formation, a densely welded, crystal-rich, rhyolitic ashflow tuff. This cored section is highly altered and devitrified, with kaolinite, quartz, chlorite, and montmorillonite replacing feldspars and much of the groundmass. Breccia zones within the tuff contain fracture fillings of hematite, limonite, goethite, jarosite, and opal. A zone of intense clay alteration, encountered in the depth interval 17.45-22.30 m, was interpreted to represent the basal vitrophyre of this unit. Underlying the Nopal Formation is the Coloradas Formation, which consists of a welded lithic-rich rhyolitic ash-flow tuff. The cored section of this unit has undergone devitrification and oxidation, and has a similar alteration mineralogy to that observed in the Nopal tuff. A sharp contact between the Coloradas tuff and the underlying Pozos Formation was observed at a depth of 136.38 m. The Pozos Formation consists of poorly sorted conglomerate containing clasts of subangular to subrounded fragments of volcanic rocks, limestone, and chert. Three thin (2-6 m) intervals of intercalated pumiceous tuffs are present within this unit. The contact between the Pozos Formation and the underlying Cretaceous limestone basement was encountered at a depth of 244.40 m. The water table is located at a depth of ~223 m. Several zones with elevated radioactivity in the PB-1 core occur above the current water table. These zones may be associated with changes in redox conditions that could have resulted in the precipitation of uraninite from downward-flowing waters transporting U from the overlying Nopal deposit. All of the intersected units have low (typically submillidarcy) matrix permeability, thus fluid flow in this area is dominated by fracture flow. These stratigraphic and rock-property observations can be used to constrain flow and transport models for the Peña Blanca natural analogue.

Summary of Rock-Property Measurements for Hong Kong TuffSamples

A series of rock-property measurements was performed on a suite of rhyolitic tuff samples from the area above the Aberdeen Tunnel of Hong Kong. The goal of this study was to determine the mechanical properties of these samples after weathering. This report contains petrographic descriptions, porosity, bulk and grain density, as well as ultrasonic measurements, elastic modulii calculations, and rock-strength determinations. Variations in rock properties are related to alteration and the presence of fractures in the tuff. Granitic rocks located adjacent to the altered tuffs would be better candidates for underground excavations.

Infiltration and Seepage Through Fractured Welded Tuff

The Nopal I mine in Pena Blanca, Chihuahua, Mexico, contains a uranium ore deposit within fractured tuff. Previous mining activities exposed a level ground surface 8 m above an excavated mining adit. In this paper, we report results of ongoing research to understand and model percolation through the fractured tuff and seepage into a mined adit both of which are important processes for the performance of the proposed nuclear waste repository at Yucca Mountain. Travel of water plumes was modeled using one-dimensional numerical and analytical approaches. Most of the hydrologic properly estimates were calculated from mean fracture apertures and fracture density. Based on the modeling results, we presented constraints for the arrival time and temporal pattern of seepage at the adit.

Deep Borehole Field Test Research Activities at LBNL

Author(s): Dobson, Patrick; Tsang, Chin-Fu; Kneafsey, Timothy; Borglin, Sharon; Piceno, Yvette; Andersen, Gary; Nakagawa, Seiji; Nihei, Kurt; Rutqvist, Jonny; Doughty, Christine; Reagan, Matthew