Newberry Volcano in Central Oregon is the site of a Department of Energy funded Enhanced Geothermal System (EGS) Demonstration Project. Stimulation and production of an EGS is a strong perturbation to the physical and chemical environment, giving rise to coupled Thermal-Hydrological-Mechanical-Chemical (THMC) processes leading to permeability changes as a result of mineral dissolution and precipitation, rock deformation, and fracture reactivation. To evaluate these processes, and to help guide EGS stimulation and reservoir development strategies, a combined native-state and reservoir model of the west flank of Newberry Volcano was created that encompasses the planned stimulation zone and a several km region of the west flank from the surface down to the supercritical region, likely close to a postulated cooling intrusive body. Temperature and pressure distributions were first modeled using TOUGHREACT with boundary conditions estimated from nearby drill holes, and compared to measurements made in the over 3 km deep NWG 55-29 drill hole. With estimates of the porosity and heat capacities for the major hydrogeologic units, thermal conductivities were calibrated by matching to the measured temperature profile. To simulate the development of the observed hydrothermal mineralogy, a reaction-transport model (THC) was developed using the pre-alteration mineralogy and shallow groundwater chemistry as the initial geochemical conditions, assuming that modeled temperature and pressure distributions were relatively constant over several thousand years. Close correspondence of modeled and observed epidote distributions support the observation that past hydrothermal activity took place under thermal gradients similar to current values, whereas calcite and sulfide abundances at depth likely require a magmatic gas component. Multicomponent geothermometry was used to estimate potential temperatures of equilibration of waters, and to evaluate the effects of kinetics on calculated mineral equilibration temperatures. The ultimate goal will be to capture both the local chemical and mechanical changes in the rock owing to stimulation as well as the potential long-term response and sustainability of the larger-scale geothermal reservoir.

ASCII text files containing grid-block name, X-Y-Z location, and multiple parameters from TOUGH2 simulation output of CO2 injection into an idealized single fault representing a dipping normal fault at the Desert Peak geothermal field (readable by GMS). The fault is composed of a damage zone, a fault gouge and a slip plane. The runs are described in detail in the following: Borgia A., Oldenburg C.M., Zhang R., Jung Y., Lee K.J., Doughty C., Daley T.M., Chugunov N., Altundas B, Ramakrishnan T.S., 2017. Carbon Dioxide Injection for Enhanced Characterization of Faults and Fractures in Geothermal Systems. Proceedings of the 42st Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, February 13-17.

Directional Cooling-Induced Fracturing (DCIF) experiments were conducted on a short, cylindrical sample of Westerly granite (diameter = 4 inches, height ~ 2 inches). Liquid nitrogen was poured in a copper cup attached to the top of the sample, and the resulting acoustic emissions (AEs) and temperature changes on the surface of the sample were monitored. The obtained AEs were used to determine the microcracking source locations and amplitude, and the associated moment tensors. Included in this submission is an animation of the AEs, a graphic displaying the temperature changes, and the measured data.

This package includes data and footage from two rounds of downhole camera surveys performed at the Sanford Underground Research Facility (SURF) on the 4850 level. The exercise was performed once on 25 May 2018 and once on 21 December 2018. On May 25th, the first round was done during fluid injection at the 164-ft stimulation zone in the injection well (E1-I). On December 21st, the second round was carried out during fluid injection at the 142-ft stimulation zone. Prior to the injections, downhole instrumentation was removed from the production well (E1-P) to allow room for the downhole camera system. The water within E1-P was then lifted out by the application of air pressure and the downhole camera system was conveyed into the production well. Finally, the water was injected into E1-I and the camera was used to scan for jetting points, or fluid entry, in E1-P. There is a survey description in this package that further describes the procedure of the survey and the overall results. Additionally, there is a detailed analysis of the surveys in the form of a PowerPoint, which includes animations/visualizations from the camera footage, presents interpretations in detail, and provides some general conclusions. Three animations, along with the two video segments that show the jetting into E1-P, are also provided. The video footage was collected using a GeoVISION Dual-Scan Micro Video Camera, the specs of which are also included in this package as a resource.

These data and test descriptions comprise a chilled circulation test conducted at the 164' fracture in the EGS Collab Experiment 1 testbed on the 4850 ft level of the Sanford Underground Research Facility. Descriptions of the meta data, design drawings for the flow testing system, and evaluation of the thermistor data are provided here. The test ran from April 2019 through early March of 2020, when testing was concluded at the experiment 1 site. These data are are complementary to the stimulation data provided in another submission which is linked below (i.e. stimulation at the 164' notch). More information about the test itself as well as the rationale and process of data processing is available on the EGS Collab Experiment 1 Long Term Circulation Test wiki page which is also linked below.

Two broadband seismometers were installed on the 4100 level and recorded for the duration of EGS Collab Experiment #2. Inspired by published data from similar instruments installed in the Aspo Hard Rock Lab, these long-period instruments aimed to measure the tilting of the drift in response to the injection of fluid into the testbed. One instrument was installed underneath the wellheads in Site A (aka the "battery" alcove) and the other was installed along the east wall of the drift, south of Site B. Due to the feet of gravel (ballast) laid along the floor of the drift, we were unable to anchor the sensors directly to the rock. As a result, the coupling of the sensors to the experiment rock volume is likely poor. In addition, there are a number of noise sources that complicate the interpretation of the data. For example, sensor BBB is installed adjacent (within 3 ft) to the rail line that runs towards the Ross shaft. Trains (motors) run along this line almost daily and produce a large signal in these data. Careful extraction of periods of interest, as well as filtering for specific signals, is necessary. The sensors are Nanometrics Trillium Compact Posthole seismometers, sensitive down to 120 seconds period. They were installed as close to the drift wall and as deep as we could manually excavate (only about 1 ft or so). The holes were leveled with sand and the sensors were placed on a paver before backfilling with sand. The hole was then covered by a bucket filled with insulation to improve the sensor's isolation from daily temperature variations, which are minor but present due to drift ventilation from the surface. Data were recorded on Nanometrics Centaur digitizers at 100 Hz. The full response information is available in the StationXML file provided here, or by querying the sensors through the IRIS DMC (see links below). These instruments were provided free of charge through the IRIS PASSCAL instrument center. The network code is XP and the station codes are BBA and BBB. The waveform data can be queried through the IRIS FDSN server using any method the user likes. One convenient option is to use the Obspy python package: https://docs.obspy.org/packages/obspy.clients.fdsn.html

Distributed fiber optic sensing was an important part of the monitoring system for EGS Collab Experiment #2. A single loop of custom fiber package was grouted into the four monitoring boreholes that bracketed the experiment volume. This fiber package contained two multi-mode fibers and four single-mode fibers. These fibers were connected to an array of fiber optic interrogator units, each targeting a different measurement. The distributed temperature system (DTS) consisted of a Silixa XT-DTS unit, connected to both ends of one of the two multi-mode fibers. This system measured absolute temperature along the entire length of fiber for the duration of the experiment at a sampling rate of approximately 10 minutes. This dataset includes both raw data in XML format from the XT-DTS, as well as a processed dataset with the sections of data pertaining only to the boreholes are extracted. We have also included a report that provides all of the relevant details necessary for users to process and interpret the data for themselves. Please read this accompanying report. If, after reading it, there are still outstanding questions, please do not hesitate to contact us. Happy processing.

Stimulation data from Experiment 1 of EGS Collab, which occurred on the 4850 ft level of the Sanford Underground Research Facility (SURF). A detailed description of the stimulation data is provided in the StimulationDataNotes.docx and is also available on the EGS Collab Wiki. A Meta Data Cheat Sheet, which describes all of the channels in the Raw CSV files, is available as well. Note that this cheat sheet is a comprehensive meta data descriptor and channels were added as the experiment evolved. This means that some columns may not be populated in early data. Additionally, we have included the chat logs from these experiments. The experiments were broadcast over teleconferencing software and real-time data displays were available to remote observers. The logs contain important observations from those personnel performing the experiment and the remote contributors. Finally, we have included summary and individual plots of all of the data for the user to compare to.

The U.S. Department of Energy's Enhanced Geothermal System (EGS) Collab project aims to improve our understanding of hydraulic stimulations in crystalline rock for enhanced geothermal energy production through execution of intensely monitored meso-scale experiments. The first experiment was performed at the 4850 ft level of the Sanford Underground Research Facility (SURF), approximately 1.5 km below the surface at Lead, South Dakota. The data reported here were collected by the continuous active-source seismic monitoring (CASSM) system (Ajo-Franklin et al., 2011). This system was permanently installed in the testbed and consisted of 17 piezoelectric sources that were recorded by 2-12 channel hydrophone arrays, 18 3-C accelerometers, and 4 3-C geophones at a Nyquist frequency of 24kHz. The source array was activated in a repeated sequence of shots (each source fired 16 times and stacked into resultant waveforms) for the duration of the experiment (April 25, 2018 - March 7, 2019) with few exceptions. Please see the attached documents describing the source / receiver geometry. The data are available in both seg2 (.dat extension) and segy (.sgy extension) format. Each segy file contains multiple seg2 files.

Core logs from the EGS Collab project Experiment 1 for the stimulation (Injection) well (E1-I), the Production well (E1-P), and monitoring wells (E1-OT, E1-OB, E1-PST, E1-PSB, E1-PDT, and E1-PDB) on the 4850 Level of SURF (the Sanford Underground Research Facility), single PDF file, 5-ft run intervals. In the monitoring well IDs, "O" indicates that the well is orthogonal to the anticipated fracture plane, "P" indicates that the well is parallel to the anticipated fracture plane, "S" indicates a shallow well, "D" indicates a deep well, "T" refers to top, and "B" refers to bottom. Logs include: experiment number; borehole ID; depth interval; run number; final packed core box number; scribe line (yes/no; red-on-right convention); logging dates; logger initials; as well as sketches of core foliation, folding, and fracturing with additional details and notes on other features of interest.

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and a transient afterslip phase. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment. This submission includes the pressure and deformation data recorded by the SIMFIP probe during the first injection at the 164 ft (50 m) notch of borehole E1-I. The injection was performed on on 05/22/2018 as part of Experiment 1 of the EGS Collab project. This data accompanies a manuscript submitted to GRL, linked in this submission.

Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).

Core logs and photos from the EGS Collab project Experiment 2 for the Top Vertical well (TV4100) and the Top Horizontal well (TV 4100) on the 4100 Level of SURF (the Sanford Underground Research Facility). The core logs are stored in a single PDF file with 5-ft run intervals. In the monitoring well IDs, "O" indicates that the well is orthogonal to the anticipated fracture plane, "T" refers to top, and "H" refers to horizontal. A core log CT scan for TV4100 and a layout image of the 4100 wells are included as well. Logs include: experiment number; borehole ID; depth interval; run number; final packed core box number; scribe line (yes/no; red-on-right convention); logging dates; logger initials; as well as sketches of core foliation, folding, and fracturing with additional details and notes on other features of interest. Shift reports include: date, location, personnel, summary of site activity, and field notes.

This submission contains the presentation slides and recordings from the first 98 EGS Collab Modeling and Simulation Working Group teleconferences. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.

This submission contains the presentation slides and recordings from EGS Collab Modeling and Simulation Working Group (MSWG) teleconferences number 99 through 128. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.

This package includes data and models that support hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF). A paper by Schwering et al. (2020) describes the deterministic basis for developing a "common" discrete fracture network (CDFN) model of significant natural fractures in EGS Collab Testbed 1 on the 4850-Level of SURF. The ReadMe for this model shows drift, wells, scanlines, fracture data, interpreted fractures, and geophysical visualizations. There is also a summary of the data that was used in this experiment and includes results from reviewing core, televiewer (TV) logs, core-TV depth/feature registration, and from mapping weeps in the 4850-Level drift. The CDFN is intended to be a baseline model of the pre-stimulated testbed (though some observations from stimulation helped inform the model).

As part of the geophysical characterization suite for the first EGS Collab tesbed, here are the baseline cross-well seismic data and resultant models. The campaign seismic data have been organized, concatenated with geometry and compressional (P-) & and shear (S-) wave picks, and submitted as SGY files. P-wave data were collected and analyzed in both 2D and 3D, while S-wave data were collected and analyzed in 2D only. Inversion models are provided as point volumes; the volumes have been culled to include only the points within source/receiver array coverage. The full models space volumes are also included, if relevant. An AGU 2018 poster by Linneman et al. is included that provides visualizations/descriptions of the cross-well seismic characterization method, elastic moduli calculations, and images of model inversion results.

This file contains the first set of tracer data for the EGS Collab testbed. The first set of tracer tests were conducted during October-November, 2018. We have included tracer data for C-dots, chloride, fluorescein, and rhodamine-B. The details about the tracer test can be found in Background and Methods of Tracer Tests (Mattson et al. (2019)) (also included in this package). References Mattson, E.D., Neupane, G., Plummer, M.A., Hawkins, A., Zhang, Y. and the EGS Collab Team 2019. Preliminary Collab fracture characterization results from flow and tracer testing efforts. In Proceedings 44th Workshop on Geothermal Reservoir Engineering, edited, Stanford University, Stanford, California.

The submission contains a .xls files consisting of 10 excel sheets, which contain combined list of pressure, saturation, salinity, temperature profiles from the simulation of CO2 push-pull using Brady reservoir model and the corresponding effective compressional and shear velocity, bulk density, and fluid and time-lapse neutron capture cross section profiles of rock at times 0 day (baseline) through 14 days. First 9 sheets (each named after the corresponding CO2 push-pull simulation time) contains simulated pressure, saturation, temperature, salinity profiles and the corresponding effective elastic and neutron capture cross section profiles of rock matrix at the time of CO2 injection. Each sheet contains two sets of effective compressional velocity profiles of the rock, one based on Gassmann and the other based on Patchy saturation model. Effective neutron capture cross section calculations are done using a proprietary neutron cross-section simulator (SNUPAR) whereas for the thermodynamic properties of CO2 and bulk density of rock matrix filled with fluid, a standalone fluid substitution tool by Schlumberger is used. Last sheet in the file contains the bulk modulus of solid rock, which is inverted from the rock properties (porosity, sound speed etc) based on Gassmann model. Bulk modulus of solid rock in turn is used in the fluid substitution.

The link points to a website at NCEDC to download the full moment tensors inversion software The moment tensor analysis conducted in the current project is based on the full moment tensor model described in Minson and Dreger (2008). The software including source, examples and tutorial can be obtained from ftp://ncedc.org/outgoing/dreger (download file pasi-nov282012.tar.gz). Performance criteria, mathematics and test results are provided by Minson and Dreger (2008), Ford et al. (2008, 2009, 2010, 2012) and Saikia (1994). References: Ford, S., D. Dreger and W. Walter (2008). Source Characterization of the August 6, 2007 Crandall Canyon Mine Seismic Event in Central Utah, Seism. Res. Lett., 79, 637-644. Ford, S. R., D. S. Dreger and W. R. Walter (2009). Identifying isotropic events using a regional moment tensor inversion, J. Geophys. Res., 114, B01306, doi:10.1029/2008JB005743. Ford, S. R., D. S. Dreger and W. R. Walter (2010). Network sensitivity solutions for regional moment tensor inversions, Bull. Seism. Soc. Am., 100, p. 1962-1970. Ford, S. R., W. R. Walter, and D. S. Dreger (2012). Event discrimination using regional moment 665 tensors with teleseismic-P constraints, Bull. Seism. Soc. Am. 102, 867-872. Minson, S. and D. Dreger (2008), Stable Inversions for Complete Moment Tensors, Geophys. J. Int., 174, 585-592. Saikia, C.K. (1994), Modified Frequency-Wavenumber Algorithm for Regional Seismograms using Filons Quadrature: Modeling of Lg Waves in Eastern North America. Geophys. J. Int., 118, 142-158.

This data file is for test 2. In this test a sample of granite with a pre-cut (man made fracture) is confined, heated and differential stress is applied. The max temperature in this this system development test is 95C. More test details can be found on the spreadsheets--note that there are 2 spreadsheets

This submission contains a link to the EDX Collaborative Workspace where the MT data collected in support of the DOE GTO 4D EGS monitoring project is stored. Daily production reports-- Oregon State University (OSU) had 6 stations running continuously. --Dynamic survey map, KML file with MT locations on the west flank. Read off at location, created excel file for locations of each NBL. Zonge has two N-S lines of MT stations, 1-x and 2-x. Created excel file for locations of each 1-x and 2-x. --In stations, each station has day file with calibration of magnetometers, 6 channels. .Z3d are proprietary data files (refer to Zonge Int'l) --MT Section: has four channels that went into it, in edi format are given frequenices, coordinate system. Tensor-- four elements of this tensor at each frequency. The tensor is complex-valued-- it has a real part and an imaginary part at each frequency. In .zxr the impedance sensor relates N-S to E-W. "r" is imaginary. "ZXYVAR" is variance, error on each impedance tensor. Transmuted into apparent resistivity "ro". phase.

The Newberry Volcano EGS Demonstration in central Oregon, a 5 year project begun in 2010, tests recent technological advances designed to reduce the cost of power generated by EGS in a hot, dry well (NWG 55-29) drilled in 2008. First, the stimulation pumps used were designed to run for weeks and deliver large volumes of water at moderate well-head pressure. Second, to stimulate multiple zones, AltaRock developed thermo-degradable zonal isolation materials (TZIMs) to seal off fractures in a geothermal well to stimulate secondary and tertiary fracture zones. The TZIMs degrade within weeks, resulting in an optimized injection/ production profile of the entire well. Third, the project followed a project-specific Induced Seismicity Mitigation Plan (ISMP) to evaluate, monitor for, and mitigate felt induced seismicity. An initial stimulation was conducted in 2012 and continued for 7 weeks, with over 41,000 m3 of water injected. Further analysis indicated a shallow casing leak and an unstable formation in the open hole. The well was repaired with a shallow casing tieback and perforated liner in the open hole and re-stimulated in 2014. The second stimulation started September 23rd, 2014 and continued for 3 weeks with over 9,500 m3 of water injected. The well was treated with several batches of newly tested TZIM diverter materials and a newly designed Diverter Injection Vessel Assembly (DIVA), which was the main modification to the original injection system design used in 2012. A second round of stimulation that included two perforation shots and additional batches of TZIM was conducted on November 11th, 2014 for 9 days with an additional 4,000 m3 of water injected. The stimulations resulted in a 3-4 fold increase in injectivity, and PTS data indicates partial blocking and creation of flow zones near the bottom of the well.

The Newberry Volcano EGS Demonstration in central Oregon, a 3 year project started in 2010, tests recent technological advances designed to reduce the cost of power generated by EGS in a hot, dry well (NWG 55-29) drilled in 2008. First, the stimulation pumps used were designed to run for weeks and deliver large volumes of water at moderate well-head pressure. Second, to stimulate multiple zones, AltaRock developed thermo-degradable zonal isolation materials (TZIMs) to seal off fractures in a geothermal well to stimulate secondary and tertiary fracture zones. The TZIMs degrade within weeks, resulting in an optimized injection/ production profile of the entire well. Third, the project followed a project-specific Induced Seismicity Mitigation Plan (ISMP) to evaluate, monitor for, and mitigate felt induced seismicity. Stimulation started October 17, 2012 and continued for 7 weeks, with over 41,000 m3 of water injected. Two TZIM treatments successfully shifted the depth of stimulation. Injectivity, DTS, and seismic analysis indicate that fracture permeability in well NWG 55-29 was enhanced by two orders of magnitude.

This submission contains 167 deviatoric moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration. Also included is a statistical representation of the properties of 751 fractures derived from the analysis of seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration Project. The locations of the fractures are taken from microseismic hypocenters, the fracture areas are derived from moment magnitudes via scaling relationships, and the azimuths (sigma 1) and dips (sigma 3) are derived from the results of stress analyses.

This submission contains 167 full moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities. Also included are the azimuth and plunge angles for the three main stress directions sigma1, sigma 2 and sigma 3 at the Prati32 EGS demonstration site in the northwest Geysers geothermal reservoir. The data are divided into 15 time periods spanning a range of five years, including two years prior to start of injection until three years post start of injection activities.

This dataset includes earthquake catalogues for the three stages of the 2022 well 16A(78)-32 stimulation provided by Geo Energie Suisse. Events in these catalogues have been visually inspected. There are additional events of lower signal to noise that were automatically detected. Those events will require additional analysis and processing.

Input and output files used for fault characterization through numerical simulation using iTOUGH2. The synthetic data for the push period are generated by running a forward simulation (input parameters are provided in iTOUGH2 Brady GF6 Input Parameters.txt [InvExt6i.txt]). In general, the permeability of the fault gouge, damage zone, and matrix are assumed to be unknown. The input and output files are for the inversion scenario where only pressure transients are available at the monitoring well located 200 m above the injection well and only the fault gouge permeability is estimated. The input files are named InvExt6i, INPUT.tpl, FOFT.ins, CO2TAB, and the output files are InvExt6i.out, pest.fof, and pest.sav (names below are display names). The table graphic in the data files below summarizes the inversion results, and indicates the fault gouge permeability can be estimated even if imperfect guesses are used for matrix and damage zone permeabilities, and permeability anisotropy is not taken into account.

Compressibility values were obtained in a range of pressures at 250degC by employing a fixed volume view cell completely filled with PAA aqueous solution and injecting CO2 at constant flow rate (0.3mL/min). Pressure increase as a function of supercritical CO2 (scCO2) mass fraction in the mixture was monitored. The plot shows the apparent compressibility of Stimufrac as a function of scCO2 mass fraction obtained in a pressure range between 2100-7000 psi at 250degC. At small mass fractions of scCO2 the compressibility increases probably due to the dissolution/reaction of CO2 in aqueous PAA and reaches a maximum at mCO2/mH2O = 0.06. Then, compressibility decreases showing a linear relationship with scCO2 mass fraction due to the continuous increase in density of the binary fluid associated to the pressure increase.

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4. Using history matching with tracers and available production data, the model should be tuned to represent the subsurface reservoir as accurately as possible 5. A large number of simulations is run using the history-matched reservoir model. Each simulation assumes a different wellbore flow rate allocation across the injection and production wells, where the individual selected flow rates do not violate the practical constraints for the corresponding wells. 6. ML models are trained using the simulation data. The code in our GitHub repository demonstrates how these models can be trained and evaluated. 7. The trained ML models can be used to evaluate a large set of candidate flow allocations with the goal of selecting the most optimal allocations, i.e., producing the largest amounts of thermal energy over the modeled period of time. The referenced paper provides more details about this optimization process

Lab-scale stimulation was performed on Coso samples obtained from a single core (1623 feet TVD, reservoir Coso CGC 18-27) using StimuFrac and control fluid in the absence of stimuli-responsive polymer.

Borehole W1 is a NQ core hole drilled at our test site in Socorro. The rock is rhyolite. Borehole W1 which was used to test gas-gas explosive mixtures is 55 feet deep with casing (pinkish in the drawing) set to 35 feet. The model is a representation of the borehole and the holes we cored around the central borehole after the test. The brown colored core holes showed dye when we filled W1 with water and slightly pressurized it. This indicates there was some path between W1 and the colored core hole. The core holes are shown to their TD in the drawing. The green plane is a fracture plane which we believe is the result of the explosions of the gas mixture in W1. Data resource is a 2D .pdf Solid Works Drawing of borehole W1.

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.

CO2 breakthrough experiments in a column packed with sand and filled with StimuFrac fluid or water to learn about transport of the stimuli (CO2) on environments where either water or the stimuli-responsive polymer aqueous solution (StimuFrac) is present. Results suggest co-injection of StimuFrac and CO2 as the potentially best alternative to deploy this novel fracking fluid.

List of triggered events recorded on LBNL's permanent EGS seismic array at Brady's geothermal field. This submission also includes links to the NCEDC EGS Earthquake Catalog Search page and to the metadata for the seismic array installed at Brady's Geothermal Field.

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 is a link to the website where DAS seismic data, collected from wells 78-32 and 78B-32 during the Utah FORGE 2022 stimulation, is now available for download. The data can be accessed at "Well 16A78-32 2022 Stimulation Seismicity Data" link in the submission under the Silixa heading. The page includes surface acquisition nodal datasets, downhole geophone data, and Silixa fiber data. Raw seismic stimulation data and the scripts to process this data is under the Silixa heading.

This data includes a document that describes the effort to collect and analyze water and gas samples from deep Utah FORGE wells 16A(78)-32, 58-32, 56-32 and 78B-32 along with additional pdf files showing ThermoChem's analyses attached as an appendix.

This is a project description video by Dr. William W. Fleckenstein related to their "Development of Multi-Stage Fracturing System and Wellbore Tractor to Enable Zonal Isolation During Stimulation and EGS Operations in Horizontal Wellbores" R&D project at Utah FORGE which is linked bellow.

These are revised catalogs, related to the April, 2022 well 16A(78)-32 stimulation (phases 1,2, & 3), provided by Geo Energie Suisse (GES) that include additional events at the start of Stage 1 and some tidying up of some locations. These catalogs also include events for additional events that were auto-located to provide a larger dataset for statistical analyses, like b-value calculations. The actual auto-locations have been removed to prevent spurious location plots being created.

This is an analysis of the pressure falloff in stage 1 fracture stimulation of FORGE well 16A(78)-32. The objective of this research is to understand the information content of the well stimulation data of FORGE Well 16A(78)-32. The Stage 1 step-rate test, a variant of the classic diagnostic fracture injection test (DFIT), contains valuable information about the success of well fracturing and the nature of resulting formation stimulation in the drainage volume of Well 16A(78). The analysis we have provided is based on the classic pressure transient analysis in petroleum reservoirs. The next step in the analysis is to use the information we have discovered in the analysis of tracer flowback data. This set of slides we have provided includes the pressure falloff analysis of the data recorded during stimulation of Stage 1 in injection Well 16A(78)-32 conducted in April of 2022. To honor multiple rate a superposition approach for linear flow regime was applied. The analysis yielded a permeability two orders of magnitude larger than permeability from cores. Our calculated permeability is essentially the effective permeability of micro- and macro-fracture system in the stimulated volume of the Well 16A(78)-32. Another observation is that after using the classic G-function plot, no closure stress was observed. This could suggest that pre-existing natural fractures were reopened during stimulation and yet had no propensity to close in accordance to the poroelastic properties.

This dataset includes files used to fit planar fractures through the preliminary earthquake catalogs of the three stages of the April 2022 well 16A(78)-32 stimulation which is linked bellow. These planar features have been used to update the FORGE reference Discrete Fracture Network (DFN) model. The files are provided to encourage other modelers to use additional workflows to find additional/alternative features. To this end, the dataset includes the cleaned earthquake catalog data translated to the FORGE reference model global reference frame, the well trajectory of 16A(78)-32 in those same coordinates, the fit 15 planar features in csv format, and a pdf file with slides illustrating the process used to fit the features. A recorded presentation of this material is available from the October 2022 FORGE Modeling and Simulation Forum which is also linked below.

This is a set of data related to the stimulation program at Utah FORGE well 16A(78)-32 during April, 2022. This includes daily reports, 1 second Pason data, tracer data, and shear stimulation data and information including a report of an evolving prognosis for the stimulation operations.

This paper discusses the progress on a project funded by the DOE Utah FORGE (Frontier Observatory for Research in Geothermal Energy) for the development of a subsurface heat exchanger for Enhanced Geothermal Systems (EGS) using unique casing sleeves cemented in place and are used first as a system for rapid and inexpensive multi-stage stimulations and second to perform conformance control functions at 225 oC. The proposed sleeves will use a single-sized dissolvable ball to open for fracture stimulation. After stimulation, and once the balls dissolve, the sleeves are open for immediate fluid injection. A separately designed wellbore tractor specific for both fluid detection and valve manipulation is then deployed to detect and control the injection entry points to create an effective EGS through paired horizontal injectors and open hole producers. The wells will be connected through multiple networks of induced and natural fractures that can be controlled throughout the field life.

This regards a high-resolution DAS microseismic dataset produced by Silixa from Utah FORGE Phase 2C seismic monitoring well 78-32 during stimulation testing of well 58-32. It is a very large dataset and as such it is currently not directly available on GDR. However, it is available for download from the Center for High Performance Computing (CHPC) at the University of Utah using the shell script below. Additional survey information and tips for running the script are available in the attached downloadable Word Document file.

Data from two Tensor Optical Fiber Strainmeters that were operational during Stages 1, 2, and 3 of the April, 2022 stimulation of well 16A(78)-32. Each csv file contains data from each stimulation stage (stage1, stage2, stage3) for both Phase 1a strainmeter installations (FS01, formerly FS-C, and FS02, formerly FS1-2) in human-readable comma-separated value text files. There are two header lines in each file describing the data contained in that column along with their units, respectively. Data have been decimated from 2 to 1 Hz to match the Pason data found in the linked GDR dataset below (16A78-32 Stimulation Pason Data). These files contain the time series spanning the same time interval as the Pason data as well as ambient data for 5 hours before the stimulation and 5 hours following shut in. The station locations were chosen based on their proximity to the borehole seismometers owned and operated by the University of Utah. See README.txt for more information.

This is the Utah FORGE well 16A(78)-32 stimulation microseismic detection and location report from Silixa LLC. The stimulation was accomplished during April, 2022.

This archive contains data from the tracer test performed during the Utah FORGE well 16A(78)-32 stimulation. This includes: 1. the raw data file from Pason in the csv format; 2. a one-page Word document from Pason that explains the headers and units used in their data file; 3. an Excel spreadsheet with concentration data and relevant Pason flow data; 4. a Word document with jpg versions of concentration charts (Figs. 1-3) and a tracer mass recovery chart (Fig. 4).

The U.S. Department of Energy's (U.S. DOE) Frontier Observatory for Research in Geothermal Energy (FORGE) is a field laboratory that provides a unique opportunity to develop and test new technologies for characterizing, creating and sustaining Enhanced Geothermal Systems (EGS) in a controlled environment. In 2018, the U.S. DOE selected a site in south-central Utah for the FORGE laboratory. Numerous geoscientific studies have been conducted in the region since the 1970s in support of geothermal development at Roosevelt Hot Springs. A vertical scientific well, 58-32, was drilled and tested to a depth of 2290 m (7515 ft) GL in 2017 on the FORGE site to provide additional characterization of the reservoir rocks. The well encountered a conductive thermal regime and a bottom hole temperature of 199degC (390degF). More than 2000 natural fractures were identified, but measured permeabilities are low, less than 30 micro-darcies. Induced fractures indicate that the maximum horizontal stress trends NNE-SSW, consistent with geologic and well observations from the surrounding area. Approximately 45 m (147 ft) at the base of the well was left uncased. A maximum wellhead pressure of 27.6 MPa (4000 psig) at an injection rate of ~1431 L/min (~9 bpm) was measured during stimulation testing in September 2017. Conventional diagnostic evaluations of the data suggest that hydraulic fracturing and shearing occurred. Estimates of the stress gradient for delta_h_min range from of 16.7 to 17.6 kPa/m (0.74 to 0.78 psi/ft). A gradient of 25.6 kPa/m (1.13psi/ft) was calculated for delta_V. In 2019, the 2017 open-hole stimulation in well 58-32 was repeated with injection rates up to 2385 L/min (15 bpm). Two additional stimulations were conducted in the cased portion of the well; one to stimulate critically stressed fractures and the second to test noncritically stressed fractures. Breakdown of the zone spanning critically-stressed fractures occurred at a surface pressure of approximately 29.0 MPa (4200 psig). Although stimulation of the noncritically stressed fractures was interrupted by failure of the bridge plug beneath the perforated interval, micro-seismic data suggests stimulation of the fractures may have been initiated at a surface pressure of 45.5 MPa (6600 psig). These stimulation results support the conclusion the Mineral Mountains granitoid is an appropriate host for EGS development. Micro-seismicity was monitored during the stimulations using surface and downhole instrumentation. Five seismometers and a nodal array of 150 seismic sensors were deployed on the surface. A Distributed Acoustic Sensing (DAS) cable and a string of 12 geophones were deployed in well 78-32, drilled to a depth of 998 m (3274 ft) GL. A broadband sensor and a high-temperature geophone were deployed in well 68-32, drilled to a depth of 303 m (994 ft) GL. More than 420 micro-seismic events were detected by the geophone string. Other instruments detected fewer events.

Pressure, temperature, and flow data from open-hole, upper perforation, and lower perforation well stimulations gathered from various tools collected at well 58-32 during phase 2C.