There is an ongoing effort at Oak Ridge National Laboratory to develop a unique experimental capability for investigating flow through porous and fractured media using neutron imaging techniques. This capability is expected to support numerous areas of investigation associated with flow processes relevant to EGS including, but not limited to: experimental visualization and measurement of velocity profiles and other flow characteristics to better inform reduced-order modeling of flow through fractures; laboratory scale validation of flow models and simulators; and a 'real-time' tool for studying geochemical rock/fluid interactions by noninvasively measuring material effects such as precipitation and dissolution in EGS-representative conditions. Demonstrating the ability of the technique to generate useful quantitative data is the primary focus at this stage of the effort. Details of the experimental setup and neutron imaging technique will be discussed in this communication, including the description of a custom designed, high pressure, neutron scattering and imaging compatible triaxial flow cell.

Substantial unexploited opportunity exists for the US, and the world, in Enhanced Geothermal Systems (EGS). This cost analysis study provides a baseline cost for a 50 MW Geothermal Power Plant in a challenging New England environment. The study then assesses how project costs would change as a function of (1) geothermal working fluid, that is, CO2 vs. water, (2) drilling technology, (3) CO2 generation technology, and (4) other geographic locations.

This submission supersedes pressure data from March 2017 which can be found as a link in the submission resources. This submission contains 3 .csv files with time series pressure data in 3 observation wells at Brady Geothermal Field as part of the PoroTomo project. These pressure files correct a time stamp issue that was in older data which did not correct for daylight savings time which occurred 13 Mar 2016 at 0900 UTC. The data here provides borehole pressures at different temperatures and times. The timeframe each resource was taken in varies between each resource and can be found in the resource descriptions.

*This submission provides corrections to GDR Submissions 844 and 845* Poroelastic Tomography (PoroTomo) by Adjoint Inverse Modeling of Data from Hydrology. The 3 *csv files containing pressure data are the corrected versions of the pressure dataset found in Submission 844. The dataset has been corrected in the sense that the atmospheric pressure has been subtracted from the total pressure measured in the well. Also, the transducers used at wells 56A-1 and SP-2 are sensitive to surface temperature fluctuations. These temperature effects have been removed from the corrected datasets. The 4th *csv file contains corrected version of the pumping data found in Submission 845. The data has been corrected in the sense that the data from several wells that were used during the PoroTomo deployment pumping tests that were not included in the original dataset has been added. In addition, several other minor changes have been made to the pumping records due to flow rate instrument calibration issues that were discovered.

Numerical simulations have shown that the use of supercritical CO2 instead of water as a heat transfer fluid yields significantly greater heat extraction rates for geothermal energy. If this technology is implemented successfully, it could increase geothermal energy production and offset atmospheric emissions of greenhouse gases. However, the impact of geochemical reactions between acidic waters in equilibrium with supercritical CO2 and the reservoir rock have not been evaluated. At issue are enhanced rock-water interactions that may reduce reservoir porosity and permeability and may exacerbate downstream scaling. The publications included in this submission aim to assess the geochemical impact of CO2 on geothermal energy production by analyzing the geochemistry of existing geothermal fields with elevated natural CO2, to measure realistic rock-water rates for geothermal systems using laboratory and field-based experiments, and to develop reactive transport models using the filed-based rates to simulate production scale impacts.

Chemical reactions pose an important but poorly understood threat to EGS long-term success because of their impact on fracture permeability. This report summarizes the dissolution rate equations for layered silicates where data were lacking for geothermal systems. Here we report updated rate laws for chlorite (Carroll and Smith 2013), biotite (Carroll and Smith, 2015), illite (Carroll and Smith, 2014), and for muscovite. Also included is a spreadsheet with rate data and rate equations for use in reactive transport simulators.

The project provides an updated Energy Return on Investment (EROI) for Enhanced Geothermal Systems (EGS). Results incorporate Argonne National Laboratory's Life Cycle Assessment and base case assumptions consistent with other projects in the Analysis subprogram. EROI is a ratio of the energy delivered to the consumer to the energy consumed to build, operate, and decommission the facility. EROI is important in assessing the viability of energy alternatives. Currently EROI analyses of geothermal energy are either out-of-date, of uncertain methodology, or presented online with little supporting documentation. This data set is a collection of files documenting data used to calculate the Energy Return On Investment (EROI) of Engineered Geothermal Systems (EGS) and erratum to publications prior to the final report. Final report is available below, or from the OSTI web site (http://www.osti.gov/geothermal/). Data in this collections includes the well designs used, input parameters for GETEM, a discussion of the energy needed to haul materials to the drill site, the baseline mud program, and a summary of the energy needed to drill each of the well designs. EROI is the ratio of the energy delivered to the customer to the energy consumed to construct, operate, and decommission the facility. Whereas efficiency is the ratio of the energy delivered to the customer to the energy extracted from the reservoir.

The primary objective of our current research is to develop a computational test bed for evaluating borehole techniques to enhance fluid flow and heat transfer in enhanced geothermal systems (EGS). Simulating processes resulting in hydraulic fracturing and/or the remobilization of existing fractures, especially the interaction between propagating fractures and existing fractures, represents a critical goal of our project. This paper details the basic methodology of our approach. Two numerical examples showing the capability and effectiveness of our simulator are also presented.

Wells for Engineered Geothermal Systems (EGS) typically occur in conditions presenting significant challenges for conventional rotary and percussive drilling technologies: granitic rocks that reduce drilling speeds and cause substantial equipment wear. Thermal spallation drilling, in which rock is fragmented by high temperature rather than mechanical means, offers a potential solution to these problems. However, much of the knowledge surrounding this drilling technique is empirical - based on laboratory experiments that may or may not represent field conditions. This paper outlines a new numerical modeling effort investigating the grain-scale processes governing thermal spallation drilling. Several factors affect spall production at the mesoscale, including grain size and size distribution, surface temperatures and material heterogeneity. To investigate the relative influence of these factors, we have conducted a series of simulations using GEODYN - a parallel Eulerian solid and fluid dynamics code. In this paper, we describe a two-dimensional model used to simulate the grain-scale processes and present preliminary results from this modeling effort.

Combined Gravity Station Locations Initial 3D gravity results from Zonge Int'l recorded for the 4D EGS Monitoring project at Newberry, during stimulation of Well 55-29 by AltaRock Energy

Hydraulic studies of drilling micropores at various depths and with various hole sizes, tubing, fluids and rates to show theoretical feasibility. WELLFLO Simulations Report separated into three parts: Step 4: Drilling 10,000 ft Wells with Supercritical Steam, Nitrogen, and Carbon Dioxide Step 5: Drilling 20,000 ft Wells with Supercritical Steam, Nitrogen, and Carbon Dioxide Step 6: Drilling 30,000 ft Wells with Supercritical Steam, Nitrogen, and Carbon Dioxide

This dataset contains maps of deformation covering Raft River, Idaho from 2007 to 2010 calculated from interferometric synthetic aperture radar data. This dataset is used in the study entitled "Inferring geothermal reservoir processes at the Raft River Geothermal Field, Idaho, USA through modeling InSAR-measured surface deformation" by F. Liu, et al. This dataset was derived from raw SAR data from the Envisat satellite missions operated by the European Space Agency (ESA) that are copyrighted by ESA and were provided through the WInSAR consortium at the UNAVCO facility. All pair directories use the image acquired on 3/11/2007 as a reference image. To view specific information for each grd file, please use the GMT command "grdinfo" - e.g., for grd file In20070311_20071111/drho_utm.grd, use terminal command: grdinfo In20070311_20071111/drho_utm.grd

Hydraulic fracturing is currently the primary method for stimulating low-permeability geothermal reservoirs and creating Enhanced (or Engineered) Geothermal Systems (EGS) with improved permeability and heat production efficiency. Complex natural fracture systems usually exist in the formations to be stimulated and it is therefore critical to understand the interactions between existing fractures and newly created fractures before optimal stimulation strategies can be developed. Our study aims to improve the understanding of EGS stimulation-response relationships by developing and applying computer-based models that can effectively reflect the key mechanisms governing interactions between complex existing fracture networks and newly created hydraulic fractures. In this paper, we first briefly describe the key modules of our methodology, namely a geomechanics solver, a discrete fracture flow solver, a rock joint response model, an adaptive remeshing module, and most importantly their effective coupling. After verifying the numerical model against classical closed-form solutions, we investigate responses of reservoirs with different preexisting natural fractures to a variety of stimulation strategies. The factors investigated include: the in situ stress states (orientation of the principal stresses and the degree of stress anisotropy), pumping pressure, and stimulation sequences of multiple wells.

Lab-scale stimulation work on non-porous fused silica (similar mechanical properties to igneous rock) was performed using pure water, pure CO2 and water/CO2 mixtures to compare back to back fracturing performance of these fluids with PNNL's StimuFrac.

A map showing location of wells permitted, drilled and seismic test, as part of validation of innovative exploration technologies done for the Newberry Volcano project in 2012.

Hydraulic shearing is an appealing reservoir stimulation strategy for Enhanced Geothermal Systems. It is believed that hydro-shearing is likely to simulate a fracture network that covers a relatively large volume of the reservoir whereas hydro-fracturing tends to create a small number of fractures. In this paper, we examine the geomechanical and hydraulic behaviors of natural fracture systems subjected to hydro-shearing stimulation and develop a coupled numerical model within the framework of discrete fracture network modeling. We found that in the low pressure hydro-shearing regime, the coupling between the fluid phase and the rock solid phase is relatively simple, and the numerical model is computationally efficient. Using this modified model, we study the behavior of a random fracture network subjected to hydro-shearing stimulation.

Project status report as of end of 2011. Key milestones and accomplishments in the year 2011 include: Proof of feasibility of analog and digital components with their lifetime tested for more than 500 hours at 300 degrees C, Demonstration of ceramic-based packaging capable of surviving 300 degrees C, Demonstration of a temperature sensor system, Completion of the design of a pressure sensor system, and Demonstration of a prototype of the digital telemetry system.

This entry contains raw data files from experiments performed on the Vulcan beamline at the Spallation Neutron Source at Oak Ridge National Laboratory using a pressure cell. Cylindrical granite and marble samples were subjected to confining pressures of either 0 psi or approximately 2500 psi and internal pressures of either 0 psi, 1500 psi or 2500 psi through a blind axial hole at the center of one end of the sample. The sample diameters were 1.5" and the sample lengths were 6". The blind hole was 0.25" in diameter and 3" deep. One set of experiments measured strains at points located circumferentially around the center of the sample with identical radii to determine if there was strain variability (this would not be expected for a homogeneous material based on the symmetry of loading). Another set of experiments measured load variation across the radius of the sample at a fixed axial and circumferential location. Raw neutron diffraction intensity files and experimental parameter descriptions are included.

This folder contains X-ray CT images and an explanation related to the shear induced permeability testing of Stripa granite

Fabricated SiC diodes are tested in the temperature range of 300 degrees C to 600 degrees C.

Low permeability geothermal reservoirs can be stimulated by hydraulic fracturing to create Enhanced (or Engineered) Geothermal Systems (EGS) with higher permeability and improved heat transfer to increase heat production. In this paper, we document our effort to develop a numerical simulator with explicit geomechanics-discrete flow network coupling by utilizing and further advancing the simulation capabilities of the Livermore Distinct Element Code (LDEC). The important modules of the simulator include an explicit finite element solid solver, a finite volume method flow solver, a joint model using the combined FEM-DEM capability of LDEC, and an adaptive remeshing module. The numerical implementation is verified against the classical KGD model. The interaction between two fractures with simple geometry and the stimulation of a relatively complex existing fracture network under different in-situ stress conditions are studied with the simulator.

Widespread adoption of geothermal energy will require access to deeply buried resources in granitic basement rocks at high temperatures and pressures. Exploiting these resources necessitates novel methods for drilling, stimulation, and maintenance, under operating conditions that are often difficult or impossible to reproduce in laboratory settings. Physically rigorous numerical modeling tools are vital to highlight potential risks, guide process optimization and reduce the uncertainties involved in developing new technologies for these environments. Lawrence Livermore National Laboratory has developed, and is constantly improving, several multi-physics solid/structural mechanics, fluid dynamics, chemistry, and discrete element codes. Integration of the LLNL simulation tools into a coherent simulation environment will provide a predictive capability for the thermomechanical response - in particular the spall and fracture - of basement rocks at high temperatures and pressures useful for drilling and other geothermal applications. This paper outlines a modeling effort investigating the processes involved in hydrothermal spallation drilling. These include interconnected phenomena on several length and time scales: from system-scale fluid dynamics and heat transfer of the high temperature jet to the rock face to the grain-scale thermomechanics of spallation. Three models are described to capture these different scale processes: a grain-scale model to investigate the onset of spallation; a particulate fluids model to simulate the transport of the produced spalls; and a borehole-scale model to represent the integrated system behavior.

Polycrystalline diamond compact (PDC) bits are routinely used in the oil and gas industry for drilling medium to hard rock but have not been adopted for geothermal drilling, largely due to past reliability issues and higher purchase costs. The Sandia Geothermal Research Department has recently completed a field demonstration of the applicability of advanced synthetic diamond drill bits for production geothermal drilling. Two commercially-available PDC bits were tested in a geothermal drilling program in the Chocolate Mountains in Southern California. These bits drilled the granitic formations with significantly better Rate of Penetration (ROP) and bit life than the roller cone bit they are compared with. Drilling records and bit performance data along with associated drilling cost savings are presented herein. The drilling trials have demonstrated PDC bit drilling technology has matured for applicability and improvements to geothermal drilling. This will be especially beneficial for development of Enhanced Geothermal Systems whereby resources can be accessed anywhere within the continental US by drilling to deep, hot resources in hard, basement rock formations.

We investigate the flow-channeling phenomenon caused by thermal drawdown in fractured geothermal reservoirs. A discrete fracture network-based, fully coupled thermal "hydrological" mechanical simulator is used to study the interactions between fluid flow, temperature change, and the associated rock deformation. The responses of a number of randomly generated 2D fracture networks that represent a variety of reservoir characteristics are simulated with various injection-production well distances. We find that flow channeling, namely flow concentration in cooled zones, is the inevitable fate of all the scenarios evaluated. We also identify a secondary geomechanical mechanism caused by the anisotropy in thermal stress that counteracts the primary mechanism of flow channeling. This new mechanism tends, to some extent, to result in a more diffuse flow distribution, although it is generally not strong enough to completely reverse flow channeling. We find that fracture intensity substantially affects the overall hydraulic impedance of the reservoir but increasing fracture intensity generally does not improve heat production performance. Increasing the injection-production well separation appears to be an effective means to prolong the production life of a reservoir. DOI: 10.1007/s00603-015-0776-0.

Compiled by MacLeod and Sherrod, 1995. Includes the geologic map of the Newberry area and two cross-sections.

A corrigendum was submitted to the journal of Geothermics on our article "Environmentally friendly, rheoreversible, hydraulic-fracturing fluids for enhanced geothermal systems" Shao et al Geothermics 58 (2015) 22-31. In the original article some permeability values were underestimated, in particular, for rock samples fractured by the stimuli-responsive fracking fluid (PAA-CO2). In addition, effective pressures were determined to be lower for three control experiments (deionized water-carbon dioxide, DIW-CO2). Therefore, we revised values of permeability and effective pressure as well as performed additional lab-scale stimulation experiments under identical conditions to further verify/update the deductions presented in the discussion section. This is the reason for the additional data introduced in the below Table 1 (grey color). The authors regret the following inadvertent errors and corresponding modifications. These modifications do not change the scientific conclusions of the article.

This submission includes a shapefile of the Opal Mound Fault, and multiple datasets of lineaments mapped in the Mineral Mountains which overlook the Utah FORGE site, hyperlinked to rose diagrams in a polygon grid shapefile.

This submission includes the geographic extent shapefile of the Milford FORGE site located in Utah, along with a shapefile of seismometer positions throughout the area, and models of basin depth and potentiometric contours.

Six samples were evaluated in unconfined and triaxial compression, their data are included in separate excel spreadsheets, and summarized in the word document. Three samples were plugged along the axis of the core (presumed to be nominally vertical) and three samples were plugged perpendicular to the axis of the core. A designation of "V"indicates vertical or the long axis of the plugged sample is aligned with the axis of the core. Similarly, "H" indicates a sample that is nominally horizontal and cut orthogonal to the axis of the core. Stress-strain curves were made before and after the testing, and are included in the word doc. The confining pressure for this test was 2800 psi. A series of tests are being carried out on to define a failure envelope, to provide representative hydraulic fracture design parameters and for future geomechanical assessments. The samples are from well 52-21, which reaches a maximum depth of 3581 ft +/- 2 ft into a gneiss complex.

Three shapefiles in this submission show the position of proposed seismic line surveys. The mid-crustal velocity anomaly file shows the extent of an anomalously low P-wave velocity zone in the subsurface. Two other files show the extent of known hydrothermal systems in the Roosevelt Hot Springs area. Another file shows the location of the proposed water pipeline to pump water from the supply wells to the deep drill site.

The individual shapefiles in this dataset delineate estimated temperature contours (20, 40, 60, and 80 deg C) at a depth of 200 m in the Milford, Utah FORGE area. Contours were derived from 86 geothermal, gradient, and other wells drilled in the area since the mid-1970s with depths greater than 50 m. Conductive temperature profiles for wells less than 200 m were extrapolated to determine the temperature at the desired depth. Because 11 wells in the eastern section of the study area (in and around the Mineral Mountains) are at higher elevations compared to those closer to the center of the basin, temperature profiles were extrapolated to a constant elevation of 200 m below the 1830 m (6000 ft) a.s.l. datum (approximate elevation of alluvial fans at the base of the Mineral Mountains) to smooth the contours across the ridges and valleys.

Newberry seeks to explore "blind" (no surface evidence) convective hydrothermal systems associated with a young silicic pluton on the flanks of Newberry Volcano. This project will employ a combination of innovative and conventional techniques to identify the location of subsurface geothermal fluids associated with the hot pluton. Newberry project drill site location map 2010. This submission contains a topographic drill site location map.

Drilling summary from validation of innovative exploration technologies for Newberry Volcano, including depths, dates, and drilling statistics from 2012

Amplitude images (both vertical and depth slices) extracted from 3D seismic reflection survey over area of Walker Ranch area (adjacent to Raft River). Crossline spacing of 660 feet and inline of 165 feet using a Vibroseis source. Processing included depth migration. Micro-earthquake hypocenters on images. Stratigraphic information and nearby well tracks added to images. Images are embedded in a Microsoft Word document with additional information. Exact location and depth restricted for proprietary reasons. Data collection and processing funded by Agua Caliente. Original data remains property of Agua Caliente.