Search

With the goal of better understanding stimulation in crystalline rock for improving enhanced geothermal systems (EGS), the EGS Collab Project performed a series of stimulations and flow tests at 1.25 and 1.5 km depths. The tests were performed in two well-instrumented testbeds in the Sanford Underground Research Facility in Lead, South Dakota, United States. The testbed for Experiment 1 at 1.5 km depth contained two open wells for injection and production and six instrumented monitoring wells surrounding the targeted stimulation zone. Four multi-step stimulation tests targeting hydraulic fracturing and nearly year-long ambient temperature and chilled water flow tests were performed in Experiment 1. The testbed for Experiments 2 and 3 was at 1.25 km depth and contained five open wells in an outwardly fanning five-spot pattern and two fans of well-instrumented monitoring wells surrounding the targeted stimulation zone. Experiment 2 targeted shear stimulation, and Experiment 3 targeted low-flow, high-flow, and oscillating pressure stimulation strategies. Hydraulic fracturing was successful in Experiments 1 and 3 in generating a connected system wherein injected water could be collected. However, the resulting flow was distributed dynamically, and not entirely collected at the anticipated production well. Thermal breakthrough was not observed in the production well, but that could have been masked by the Joule-Thomson effect. Shear stimulation in Experiment 2 did not occur – despite attempting to pressurize the fractures most likely to shear – because of the inability to inject water into a mostly-healed fracture, and the low shear-to-normal stress ratio. The EGS Collab experiments are described to provide a background for lessons learned on topics including induced seismicity, the correlation between seismicity and permeability, distributed and dynamic flow systems, thermoelastic and pressure effects, shear stimulation, local geology, thermal breakthrough, monitoring stimulation, grouting boreholes, modeling, and system management.

Undocumented Orphaned Wells (UOWs) are wells without an operator that have limited or no documentation with regulatory authorities. An estimated 310,000 to 800,000 UOWs exist in the United States (US), whose locations are largely unknown. These wells can potentially leak methane and other volatile organic compounds to the atmosphere, and contaminate groundwater. In this study, we developed a novel framework utilizing a state-of-the-art computer vision neural network model to identify the precise locations of potential UOWs. The U-Net model is trained to detect oil and gas well symbols in georeferenced historical topographic maps, and potential UOWs are identified as symbols that are further than 100 m from any documented well. A custom tool was developed to rapidly validate the potential UOW locations. We applied this framework to four counties in California and Oklahoma, leading to the discovery of 1301 potential UOWs across >40,000 km2. We confirmed the presence of 29 UOWs from satellite images and 15 UOWs from magnetic surveys in the field with a spatial accuracy on the order of 10 m. This framework can be scaled to identify potential UOWs across the US since the historical maps are available for the entire nation.

In the United States, hundreds of thousands of undocumented orphan wells have been abandoned, leaving the burden of managing environmental hazards to governmental agencies or the public. These wells, a result of over a century of fossil fuel extraction without adequate regulation, lack basic information like location and depth, emit greenhouse gases, and leak toxic substances into groundwater. For most of these wells, basic information such as well location and depth is unknown or unverified. Addressing this issue necessitates innovative and interdisciplinary approaches for locating, characterizing, and mitigating their environmental impacts. Our survey of the United States revealed the need for tools to identify well locations and assess conditions, prompting the development of technologies including machine learning to automatically extract information from old records (95%+ accuracy), remote sensing technologies like aero-magnetometers to find buried wells, and cost-effective methods for estimating methane emissions. Notably, fixed-wing drones equipped with magnetometers have emerged as cost-effective and efficient for discovering unknown wells, offering advantages over helicopters and quadcopters. Efforts also involved leveraging local knowledge through outreach to state and tribal governments as well as citizen science initiatives. These initiatives aim to significantly contribute to environmental sustainability by reducing greenhouse gases and improving air and water quality.

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. We use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

Continuous microseismic monitoring using three-component (3C) accelerometers deployed in multiple boreholes allows for tracking the detailed evaluation of mesoscale (∼10 m scale) fracture growth during the fracture stimulation experiments at the first Enhanced Geothermal Systems (EGS) Collab testbed. Building on a well-constrained microseismic event catalog, we invert for moment tensor of the events to better understand the fracture geometry and stress orientations. However, it is challenging because of the unknown orientation of 3C accelerometers and low signal-to-noise-ratio nature of high-frequency (several kHz) monitoring. To address these challenges, we first perform the hodogram analysis on the continuous active-source seismic monitoring (CASSM) data to determine the orientations of the 18 3C accelerometers. We then apply the principal component analysis (PCA) to the observed microseismic waveforms to improve the signal-to-noise ratios. We perform a grid search for the full moment tensor by fitting the PCA-denoised waveforms at a frequency range of 5 to 8 kHz. The moment tensor results show both the creation of hydraulic fractures and the reactivation of natural fractures during the hydraulic stimulations. Our stress inversion based on the inverted moment tensors reveals the alteration of stress regime caused by hydraulic fracture stimulations.

Fracture injection/shut-in tests are commonly used to measure the state of stress in the subsurface. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of (Hayashi and Haimson, 1991). In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. The calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

Enhanced geothermal systems (EGS) require cost-effective monitoring of fracture networks. We validate the capability of using borehole distributed acoustic sensing (DAS) ambient noise for fracture monitoring using core photos and core logs. The EGS Collab project has conducted 10 m scale field experiments of hydraulic fracture stimulation using 50-60 m deep experimental wells at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The first EGS Collab testbed is located at 1616.67 m (4850 ft) depth at SURF and consists of one injection well, one production well, and six monitoring wells. All wells are drilled subhorizontally from an access tunnel called a drift. The project uses a single continuous fiber-optic cable installed sequentially in the six monitoring wells to record DAS data for monitoring hydraulic fracturing during stimulation. We analyze 60 s time records of the borehole DAS ambient noise data and compute the noise root-mean-square (rms) amplitude on each channel (points along the fiber cable) to obtain DAS ambient noise rms amplitude depth profiles along the monitoring wellbore. Our noise rms amplitude profiles indicate amplitude peaks at distinct depths. We compare the DAS noise rms amplitude profiles with borehole core photos and core logs and find that the DAS noise rms amplitude peaks correspond to the locations of fractures or lithologic changes indicated in the core photos or core logs. We then compute the hourly DAS noise rms amplitude profiles in two monitoring wells during three stimulation cycles in 72 h and find that the DAS noise rms amplitude profiles vary with time, indicating the fracture opening/growth or closing during the hydraulic stimulation. Our results demonstrate that borehole DAS passive ambient noise can be used to detect fractures and monitor fracturing processes in EGS reservoirs.

The EGS Collab project performed well-monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to inform challenges in implementing enhanced geothermal system (EGS) technology. This project, supported by the US Department of Energy, gathered data and observations from the field tests and compared these to simulation results to understand processes and to build confidence in numerical modeling of the processes. The project consisted of 3 Experiments, each comprising test and testbed design, many individual tests, numerical simulation, and analysis. The Experiments were performed in two deep underground testbeds at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Field experiments are now complete, significant data sets have been collected and analyzed, and some analysis continues. Experiments using underground test facilities have many advantages in that they allow:•Three-dimensional characterization of the stimulated volume by complementary geophysical methods surrounding the experiment•Using techniques that are currently not applicable under geothermal condition to provide processes insight •Comprehensive tracer testing and detailed characterization of complex fluid movements•Understanding the geometry of the stimulated network at the meso-scale and its implications for effective fracture surface area, rock block size, and heat exchange.Underground testing has its own set of complications, however, which affect the ability to perform tests as desired and affect the experiment results. Included here are the inability to flow at desired injection rates due to stress gradients caused by drift cooling, and the need to strongly limit induced seismicity because the distance to people and equipment.Experiment 1 examined hydraulic fracturing at a depth of 1.5 km in a well-characterized phyllite. Eight subhorizontal boreholes were used in this Experiment. Geophysical monitoring instrumentation was deployed in six boreholes to monitor stimulation events and flow tests. The other two boreholes were used to perform and carefully measure water injection and production. More than a dozen stimulations and nearly one year of flow tests in the testbed were performed. Detailed observations of processes occurring during stimulation and dynamic flow tests were collected and analyzed. Flow tests using ambient-temperature and chilled water were performed with intermittent tracer tests to examine system behavior. We achieved adaptive control of the tests using close monitoring of rapidly disseminated data and near-real-time simulation. Numerical simulation was critical in answering key experimental design questions, forecasting fracture behavior, and analyzing results. We were successful in performing many simulations in near-real-time in conjunction with the field experiments, with more detailed simulations performed later. Experiment 2 was intended to examine hydraulic shearing of natural fractures at a depth of 1.25 km in amphibolite. The stresses, rock type, and fracture conditions are different than in Experiment 1. The testbed consists of 9 boreholes, in addition to 2 exploratory characterization boreholes. Four boreholes drilled as two fans of 2 monitoring holes contained grouted-in monitoring sensors. The remaining five open boreholes drilled in a five-spot pattern were adaptively used for injection, production, and monitoring. Approximately five fracture set orientations were encountered in the testbed along with a low-stress rhyolite sill at 35 m below the access drift in exploratory well TV4100. The testbed was designed to optimize the potential for shear stimulation while also avoiding the low-stress rhyolite. Experiment 2 focused on stimulating a fracture in the most likely orientation to shear, however shear stimulation did not occur probably due to cementation from natural secondary mineralization. Other fracture sets encountered were also cemented and had orientations less likely to shear.Experiment 3 consisted of several stimulations in the same testbed as Experiment 2 allowing different stimulation approaches including ramped-rate injection, rapid injection, and oscillating-pressure injection. Ultimately these methods created hydraulic fractures, one of which was used for a medium-duration cold water injection test. The major findings of the EGS Collab Project include: 1.Significant shear stimulation did not occur during our stimulation attempts. Shear stimulation may occur but under a limited set of conditions not encountered.2.Our stimulations resulted in hydraulic fractures that required hydraulic propping. Pumping at pressures exceeding the minimum principal stress may not be feasible in an enhanced geothermal system.3.The systems we generated were complex hydraulic fracture/natural fracture systems, and these systems changed over time in response to applied pressures and flowrates and to unknown stimuli.4.The project attempted alternative stimulation methods, which did not provide significant flow improvement. 5.Thermal breakthrough was not achieved as designed, most likely because flow to production boreholes was not adequate. 6.The combination of geophysical tools used provided excellent understanding of many important processes. 7.Microearthquakes (MEQs) didn’t necessarily identify flow paths.8.Engineering tools bounding expected seismicity are needed.This report provides a summary of tests and analyses performed for EGS Collab Experiment 2 (Shear Stimulation in Testbed 2) and Experiment 3 (Alternative Stimulation methods). Much of the EGS Collab work has been published in journals and conference papers, presented in conferences, included in written reports, and submitted in data sets to the Geothermal Data Repository (GDR). The entirety of these written works is included as an appendix to this report, and this report serves as a summary and framework pointing to these published papers, presentations, and reports.

Fundamental mechanical and creep behavior of minerals at elevated temperatures is relevant to many subsurface exploration and resource recovery processes, yet remains largely elusive. This study uses instrumented indentation to measure microscale deformation and nanoscale creep of biotite, albite, and quartz, which are the major mineral phases of granitic rocks, at temperatures up to 400 °C. The results show that the elastic modulus of biotite is not only the lowest at room temperature but also decays the most at elevated temperature (i.e., by 25% reduction at 400 °C) compared to that of albite and quartz. The creep deformation within 20 s doubles when the temperature increases from 20 °C to 400 °C for albite, quartz, and biotite, and biotite showed at least three times the overall creep deformation compared to quartz and albite under the same temperature. The indentation creep deformation can be well characterized using a logarithmic time model, showing both transient and constant-rate creep deformation. The secondary creep rate for biotite is about three times that of albite and quartz at the same temperature, and as the temperature increases from 20 °C to 400 °C, the rate of secondary creep becomes approximately 10 times faster for the three tested minerals. Both transient creep deformation and the rate of secondary creep can be empirically correlated to the elastic moduli, which allows quick estimates of the temperature-dependent creep behavior of granitic minerals using their elastic properties, especially when those creep constants are not readily available at elevated temperatures. These results enhance the understanding of temperature-dependent creep deformation at small scales and provide insight into mineral-level damage in granitic rocks due to temperature and long-term deformation.

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

Lost circulation is one of the most common and expensive problems facing geothermal energy development, representing up to 30% of drilling costs. We examined drilling records from four geothermal fields—McGinness Hills in central Nevada, Don A. Campbell and Steamboat Hills in western Nevada, and Puna Geothermal Venture on the Big Island of Hawai'i— to identify geologies most prone to lost circulation, as well as common mitigation strategies. Depths of lost circulation events varied, but their frequency often increased in the production interval. Lost circulation commonly occurred near fault intersections, and heavily faulted fields like McGinness Hills and Don A. Campbell showed secondary mineralization within approximately 100 m (328 ft) or less of where circulation was lost. Lost circulation mitigation strategies included using locally available materials (e.g., cotton seed hulls) as well as more expensive proprietary lost circulation materials, cement plugs above the reservoir, and drilling blind with aerated, polymer-based mud in the production zone. Addressing lost circulation using a well thought out decision-making approach and materials above the reservoir will save time and cost, and provide needed well integrity. Mitigation often requires a series of steps, typically applied from perceived least expensive to most, and are dependent on the severity and location in the well where circulation was lost and availability of materials. Placing cement plugs can cure lost circulation events, however these plugs are often expensive, time-consuming, and may not be successful.

Many challenges related to carbon-dioxide ([Formula: see text]) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90[Formula: see text] (vertical fracture plane) to 30[Formula: see text]. The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0[Formula: see text]). However, for fracture inclinations greater than 10[Formula: see text], the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester [Formula: see text] through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

The EGS Collab project, supported by the US Department of Energy, is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to address challenges in implementing enhanced geothermal systems (EGS). Data and observations from the field tests are compared to simulations to understand processes and build confidence in numerical modeling of the processes. We have completed Experiment 1 (of 3), which examined hydraulic fracturing in a well-characterized underground fractured phyllite test bed at a depth of approximately 1.5 km at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Testbed characterization included fracture mapping, borehole acoustic and optical televiewers, full waveform sonic, conductivity, resistivity, temperature, campaign p- and s-wave investigations and electrical resistance tomography. Borehole geophysical techniques including passive seismic, continuous active source seismic monitoring, electrical resistance tomography, fiber-based distributed strain, distributed temperature, and distributed acoustic monitoring, were used to carefully monitor stimulation events and flow tests. More than a dozen stimulations and nearly one year of flow tests were performed. Quality data and detailed observations were collected and analyzed during stimulation and water flow tests using ambient temperature and chilled water. We achieved adaptive control of the tests using real-time monitoring and rapid dissemination of data and near-real-time simulation. More detailed numerical simulation was performed to answer key experimental design questions, forecast fracture propagation trajectories and extents, and analyze and evaluate results. Data are freely available from the Geothermal Data Repository. Experiment 2 examines the potential for hydraulic shearing in amphibolite at a depth of about 1.25 km at SURF. This site has a different set of stress and fracture conditions than Experiment 1. The Experiment 2 testbed consists of nine subhorizontal boreholes configured in two fans of two boreholes which surround the testbed and contain grouted-in electrical resistance tomography, seismic sensors, active seismic sources and distributed fiber sensors. A “five-spot” set of test wells that extends from a custom mined alcove includes an injection well and four production/monitoring wells. The testbed was characterized geophysically and hydrologically, and three stimulations have been performed using the Step-Rate Injection Method for Fracture In-Situ Properties (SIMFIP) tool to measure strains, and a new strain quantifying tool (downhole robotic strain analysis tool -DORSA) was deployed in a monitoring hole during stimulation. Real-time data were broadcast during stimulations to allow real-time response to arising issues.

The water-gas relative permeability is an important parameter to characterize multiphase flow in sediments. To study the water-gas relative permeability of hydrate-bearing porous media, multiphase flow simulations were carried out at the pore scale using the lattice Boltzmann method. The effects of hydrate saturation and hydrate-growth habits on the water-gas relative permeability, which is scaled by the relative permeability considering the hydrate only, were evaluated in a two-dimensional porous medium. Results show that the increase of hydrate saturation causes the decrease of water-gas effective permeability as expected. However, the effect of hydrate saturation on the water-gas relative permeability is different from that of hydrate saturation on the water-gas effective permeability. The water-gas relative permeability increases with the increase of hydrate saturation in the pore-filling case. The water-gas relative permeability decreases with the increase of hydrate saturation in the grain-coating case. The wettability of solid phase has a different effect on the relative permeability of wetting phase and nonwetting phase. The Jamin effect (phase blocking) was observed and may exist in the production of gas from natural gas hydrate reservoirs. This seriously affects the multiphase flow characteristics. The changes of microscale fluid distribution effect the changes of water-gas relative permeability. The relationship between the water-gas relative permeability and the characterization parameters of microscale fluid distribution was analyzed.

Feeding problems are common in autistic children. They include eating a limited repertoire of foods and often impact on health and quality of life. To examine the experiences of parents of autistic children who have feeding problems, and their service needs. A pragmatic, qualitative research approach was adopted. Six participants were recruited using purposive and snowball sampling. Data was collected via semi-structured interviews. Participants also submitted photographs with captions. These illustrated their child's feeding problems and the effect on them and their family. Interviews were transcribed. They were analyzed alongside the photos and captions using inductive thematic analysis. Rigor was ensured via member checking, peer review of themes, and reflexive journaling. Five themes emerged. Parents and children found 'mealtimes difficult.' Parents were 'unsure when or whether to pursue help.' They felt 'let down by a lack of support.' 'Lowered expectations of parents and professionals' prevailed about their child's feeding problems. Parents desired 'services that are easily accessible and supportive.' This is the first qualitative study of parents' experiences of feeding problems in autistic children in the UK. It has been the first study to incorporate participatory photography. These problems have a negative impact on the child's learning and the child's, and family's wellbeing. Support needs are largely unmet. This study has provided rich insight into autistic children's feeding problems. The findings will inform the design of better support for autistic children and their families.

We present hydro-mechanical measurements that characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens from the enhanced geothermal system (EGS) Collab Project's Experiment 1 and 2 sites (E1 and E2) at the Sanford Underground Research Facility. We employed a triaxial direct shear method augmented with X-ray imaging to perform hydroshearing (injection-induced shearing) and mechanical shearing on naturally fractured specimens at in situ stress conditions. Measurements included fracture permeability, strength, stress-dependent aperture, shear dilation, and frictional strength. Results reveal that in situ natural fractures must be permeable, weak, and shear-oriented to be hydrosheared, and only a subset of the observable in situ fractures were suitable. When sheared, the fracture permeability typically increased by a factor of 10 or more and this increase was retained over time. However, shear slip did not always result in permeability increase. High phyllosilicate content associated with exceptionally weak fractures exhibited poor or even decreased permeability after stimulation. These measurements in combination with site data were used to conduct a slip-tendency analysis for different fracture sets, and we selected the top candidate natural fractures for hydroshearing at the EGS Collab sites. We also found that the lower in situ shear stress and stronger fractures at the E2 site make hydroshearing more challenging than at the E1 site. Overall, the methods and analysis used in our work can be applied to any geothermal project to identify in situ joint sets that are best suited for hydroshearing, which in turn can help to optimize well placement and energy production.

The EGS Collab project is performing well-monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to inform challenges in implementing enhanced geothermal systems (EGS). This project, supported by the US Department of Energy, is gathering data and observations from the field tests and comparing these to simulation results to understand processes and to build confidence in numerical modeling of the processes. Experiment 1 (now complete) examined hydraulic fracturing in an underground test bed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km in a well-characterized phyllite. Geophysical monitoring instrumentation in six of eight sub-horizontal boreholes monitored stimulation events and flow tests. The other two boreholes were used to perform and carefully measure water injection and production. More than a dozen stimulations and nearly one year of flow tests in the testbed were performed. Detailed observations of processes occurring during stimulation and dynamic flow tests were collected and analyzed. Flow tests using ambient-temperature and chilled water were performed with intermittent tracer tests to examine system behavior. We achieved adaptive control of the tests using close monitoring of rapidly disseminated data and near-real-time simulation. Numerical simulation was critical in answering key experimental design questions, forecasting fracture behavior, and analyzing results. We were successful in performing many simulations in near-real-time in conjunction with the field experiments, with more detailed simulations performed later. The primary objective of Experiment 2 is to examine hydraulic shearing of natural fractures at a depth of 1.25 km in amphibolite at SURF. The stresses, rock type, and fracture conditions are different than in Experiment 1. The testbed consists of 9 boreholes, in addition to two earlier-drilled characterization boreholes. One borehole is used for injection, two fans of 2 monitoring wells have several geophysical measurement tools grouted in, and four open boreholes surrounding the injection hole are adaptively used for production and monitoring. We have encountered approximately five fracture set orientations in the testbed, and designed our testbed accordingly to maximize the potential for shear stimulation. Three stimulations have been performed to date from the injection borehole, each intersecting at least one production borehole. Different methods have been used for each stimulation, including a ramped flow, a high flow rate, and oscillating pressure.

4D subsurface visualization using a web platform can provide improved communication, education and outreach to non-experts and stakeholders. It allows for the improved understanding of complex relationships and interactions that occur in inaccessible locations where many experimental testbeds exist. The use of a web-based visualization tool lessens the cognitive impact further by reducing the need of specialized software. Users can define their own exploration of a 3D or 4D scene, adding and removing data as needed, and they are able to make informed decisions based on accurate model visualizations, while providing a collaborative tool that can be accessed on any platform using a browser window. Here we show the tool with example geophysical data collected from the EGS Collab project that includes borehole logs and preliminary conductivity model results. The tool is also able to point to real time streams to display up-to-the-second data as an experiment is ongoing. This level of latency can also provide operators with essential information that can direct an experiment's progress. Additionally, this tool can leverage augmented and virtual reality (AR, VR) capabilities of certain mobile devices and head-mounted displays, providing further engaging visualization possibilities.

Heat recovery from an enhanced geothermal system (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS commercialization but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Heat recovery from an enhanced geothermal system (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS commercialization but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

This document provides a distillation of learnings of the EGS Collab (Collab) project for several topics of interest requested by the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) project. This document was written following the completion of Collab Experiment 1, but prior to beginning Collab Experiment 2. Recommendations offered herein are perspectives offered by Collab scientists in support of FORGE.The United States has an enormous indigenous renewable energy potential from enhanced geothermal systems (EGS). To realize this potential, the US Department of Energy (DOE) Geothermal Technologies Office (GTO) has made significant investments into research to eliminate impediments to developing EGS. Two major current projects are Collab initiated in 2017, and the FORGE project initiated with a site selection process in 2015. Collab is a collaborative multi-national-lab, university, and commercial entity research endeavor bringing together a team of skilled and experienced subsurface process modeling, monitoring, and experimentation researchers and engineers to focus on intermediate-scale EGS reservoir creation processes and related model validation in crystalline rock. Collab is utilizing readily accessible underground facilities to refine the understanding of rock mass response to stimulation using experiments on the order of 10 m scale under EGS-relevant stress conditions. Experimental results from Collab stimulation, flow, tracer, and thermal tests are being used to validate coupled thermal-hydrological-mechanical-chemical (THMC) modeling approaches applicable to EGS. Collab is also testing and improving conventional and novel field monitoring tools. The project focuses on understanding and predicting permeability enhancement and evolution in crystalline rock. This focused research includes creating sustained and distributed permeability for heat extraction by generating new fractures that complement existing natural fractures. FORGE has the mission of establishing an EGS field test site that enables cutting-edge research and testing for EGS technology to identify a replicable, commercial pathway to EGS. The FORGE team is developing the EGS field test site near Milford, Utah. The two projects differ in some attributes, including spatial scale (Collab / FORGE - 10-m / reservoir), access to the rock (short boreholes and nearby instruments / deep wells and standard field geophysical equipment), environmental conditions (cool rock at reasonable stress / hot rock at reasonable stress), focus (direct investigation / development of a testbed and management of a research program), and project structure (integrated team / science and engineering framework supporting many individual research teams). The specific topics of this report requested by FORGE are: 1) data processing, annotation, and integration; 2) processing seismic data, 3) connecting geophysics, fractures, and flow systems, and 4) challenges of modeling fractured flow systems. Each of these topics is briefly discussed in this executive summary and covered in substantially greater detail in the main document. Numerous conference and journal papers are cited in the report that provide additional information. Many lessons have been learned. Some are not directly transferable to FORGE, but the underpinnings of the lessons may be applicable.

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried outin situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stressfield, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project alsoincluded pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability(k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the SanfordUnderground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the PrecambrianPoorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the twooutermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert(floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocityheterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Fieldmeasurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontalmaximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controllingfracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20ft), depending on the unknown amount of leak-off.

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk involved with developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated rocks, bulk electrical conductivity is sensitive to both primary and secondary porosity (i.e., matrix porosity and fractures), and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a dense fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150-250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: This increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: A thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. These experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

Oil and gas shales are a system where understanding the mechanical properties at the microscale is of paramount importance, e.g. to better understand the behavior of proppant-shale contacts and their role in the evolution of propped fractures in unconventional reservoirs. This work shows for the first time an experiment coupling indentation testing with in situ X-ray imaging in a Green River shale sample. A full compliance curve has been measured with the sample in water, allowing to visualize the indentation of the sample in function of axial load, in a purpose-built system for combined in situ indentation and X-ray imaging. A series of 3D datasets were used for a digital volume correlation study to obtain local strain fields. This analysis has been complemented with the analysis of cracks. Finally, geomechanical modeling has been carried out to replicate and generalize the observed behavior in the shale. This study validated this experimental approach, providing a breakthrough in understanding micro-mechanics in shales, and demonstrates how this class of experiments can be important for studies involving the prediction of the evolution of propped fractures in shale reservoirs, with possible applications in a much larger number of application fields (geothermal, materials science, etc.)

In cryogenic fracturing, a rock surface exposed to cryogenic fluids undergoes a large thermal gradient, and the resultant local tensile stress overcomes rock strength and initiates fractures. This study investigates the development of cracks generated from the cryogenic treatment of a borehole under no external confining stress on specimens. The experiments were performed on transparent PMMA specimens to observe fracture proliferation around boreholes. Liquid nitrogen was flowed through the boreholes to cool the borehole surface. The results show that initial fracture growth is characterized by abrupt starts and stops, and as the fracture propagates outward, the growth appears more continuous. In an early stage, horizontal/radial fractures and vertical fractures are the defining patterns. Horizontal fractures tend to be separated by a specific exclusion distance (i.e., spacing between cracks). While distinct horizontal/vertical fractures and exclusion distance manifest themselves at an early stage, fractures resulting from fracture interactions and curvatures can develop into complex shapes at later stages. Cryogenic thermal loading induces distinctively curved fractures. The tendency of curvature may prevent greater penetration. An increase in the borehole pressure during liquid nitrogen flow, however, can lessen fracture tortu-osity and facilitate radial propagation. A high flow pressure and rate are also advantageous in that they accelerate cooling and fracture propagation.

Oil and gas shales are a system where understanding the mechanical properties at the microscale is of paramount importance, e.g. to better understand the behavior of proppant-shale contacts and their role in the evolution of propped fractures in unconventional reservoirs. This work shows for the first time an experiment coupling indentation testing with in situ X-ray imaging in a Green River shale sample. A full compliance curve has been measured with the sample in water, allowing to visualize the indentation of the sample in function of axial load, in a purpose-built system for combined in situ indentation and X-ray imaging. A series of 3D datasets were used for a digital volume correlation study to obtain local strain fields. This analysis has been complemented with the analysis of cracks. Finally, geomechanical modeling has been carried out to replicate and generalize the observed behavior in the shale. This study validated this experimental approach, providing a breakthrough in understanding micro-mechanics in shales, and demonstrates how this class of experiments can be important for studies involving the prediction of the evolution of propped fractures in shale reservoirs, with possible applications in a much larger number of application fields (geothermal, materials science, etc.)

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, we show the classic fingering and transitions in flow regimes at (Formula presented.) (capillary number) from −7.2 to −3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO2 (scCO2) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

Estimation of CO2 storage capacity in fractured porous reservoirs requires a better understanding of the CO2-water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO2-water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. The core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO2-water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO2 saturations in the matrix increaze with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO2 storage efficiency in fractured porous media.

The final version of the above article was posted prematurely on 16 July 2021, owing to a technical error. The final, corrected version of record will be made fully available at a later date.

We provide direct observations of the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ∼1.5 km depth in the crystalline rock at the Sanford Underground Research Facility (USA). Micro-shearing of the borehole initiates on a foliation plane at 61% of the minimum principal stress σ3. As the fluid pressure increases further up to 112% of σ3, 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 important it is to consider monitoring zonal deformation in conjunction with pressure and flow to better manage the complex hydromechanical evolution of the growing fracture.

This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rear-rangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.

Energy extraction from subsurface reservoirs is important for addressing the increasing energy demand and environmental concerns such as global warming. However, the characterization of subsurface reservoirs, particularly reservoirs dominated by fracture networks remains a challenge due to the lack of means to directly observe subsurface processes. This study explores the feasibility and efficacy of characterizing fracture flow and transport processes in an enhanced geothermal system (EGS) testbed through stochastic tracer modeling. There are two enabling factors that allow application of stochastic modeling to characterize a subsurface reservoir. First, an abundance of geological and geophysical measurements enables the development of a high-fidelity and well-constrained fracture network model. Second, high-performance computing (HPC) allows running massive realizations efficiently. Six conservative tracer tests were stochastically modeled and produced satisfactory realizations that successfully reproduce field tracer recovery data from each tracer test. The evolution of flow and transport processes in the fracture network was then analyzed from these satisfactory realizations. The present study demonstrates that stochastic tracer modeling on a high-fidelity fracture network model is feasible and can provide important insights regarding flow and transport characteristics in subsurface fractured reservoirs.

Energy extraction from subsurface reservoirs is important for addressing the increasing energy demand and environmental concerns such as global warming. However, the characterization of subsurface reservoirs, particularly reservoirs dominated by fracture networks remains a challenge due to the lack of means to directly observe subsurface processes. This study explores the feasibility and efficacy of characterizing fracture flow and transport processes in an enhanced geothermal system (EGS) testbed through stochastic tracer modeling. There are two enabling factors that allow application of stochastic modeling to characterize a subsurface reservoir. First, an abundance of geological and geophysical measurements enables the development of a high-fidelity and well-constrained fracture network model. Second, high-performance computing (HPC) allows running massive realizations efficiently. Six conservative tracer tests were stochastically modeled and produced satisfactory realizations that successfully reproduce field tracer recovery data from each tracer test. The evolution of flow and transport processes in the fracture network was then analyzed from these satisfactory realizations. The present study demonstrates that stochastic tracer modeling on a high-fidelity fracture network model is feasible and can provide important insights regarding flow and transport characteristics in subsurface fractured reservoirs.

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. Here, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

The EGS Collab project, supported by the US Department of Energy, is addressing challenges in implementing enhanced geothermal systems (EGS). This includes improving understanding of the stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The project is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory. Data and observations from the field test are compared to simulations to understand processes and to build confidence in numerical modeling of the processes. In Experiment 1, we examined hydraulic fracturing an underground test bed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km. We drilled eight sub-horizontal boreholes in a well-characterized phyllite. Six of the boreholes were instrumented with many sensor types to allow careful monitoring of stimulation events and flow tests, and the other two boreholes were used for water injection and production. We performed a number of stimulations and flow tests in the testbed. Our monitoring systems allowed detailed observations and collection of numerous data sets of processes occurring during stimulation and during dynamic flow tests. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Data were rapidly analyzed, allowing adaptive control of the tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results. Many simulations were performed in near-real-time in conjunction with the field experiments, with more detailed process study simulations performed on a longer timeframe. Experiment 2 will examine hydraulic shearing in a test bed being built at the SURF at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1. Five sets of fracture orientations were considered in design, and three orientations seem to be consistently observed.

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. Here we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

There is an interest to ensure sub-saturated water content in lines containing carbon dioxide in applications such as enhanced oil recovery and carbon sequestration, to reduce risks of hydrate blockage and corrosion. The water content of carbon dioxide at various temperatures and pressures has been measured in the past, but there is no consistent set of measurements that could be used for carbon dioxide storage and transportation design work. The solubility of water in a carbon dioxide rich gas phase at hydrate forming conditions was measured in this work. Pressures ranged from 12.06 to 29.30 bar along two isotherms, 1 °C and −7 °C, all within the gaseous carbon dioxide and hydrate stability zone. For the first time in these types of measurements, the solid phase was also characterized and confirmed to be carbon dioxide hydrate via X-ray computed tomography, simultaneous with water content measurements of the gas phase. Once carbon dioxide hydrate conversion had reached a maximum value (65% estimated by X-ray computed tomography), the equilibrium water content was measured. Prior to reaching this maximum carbon dioxide hydrate conversion, the water content in carbon dioxide was observed to decrease as liquid water converted to carbon dioxide hydrate. This slow conversion to hydrate, metastability of the hydrate phase, or unexpected phases may be responsible for the large discrepancy between prior data sets for similar carbon dioxide water content measurements.

There is an interest to ensure sub-saturated water content in lines containing carbon dioxide in applications such as enhanced oil recovery and carbon sequestration, to reduce risks of hydrate blockage and corrosion. The water content of carbon dioxide at various temperatures and pressures has been measured in the past, but there is no consistent set of measurements that could be used for carbon dioxide storage and transportation design work. The solubility of water in a carbon dioxide rich gas phase at hydrate forming conditions was measured in this work. Pressures ranged from 12.06 to 29.30 bar along two isotherms, 1 °C and −7 °C, all within the gaseous carbon dioxide and hydrate stability zone. For the first time in these types of measurements, the solid phase was also characterized and confirmed to be carbon dioxide hydrate via X-ray computed tomography, simultaneous with water content measurements of the gas phase. Once carbon dioxide hydrate conversion had reached a maximum value (65% estimated by X-ray computed tomography), the equilibrium water content was measured. Prior to reaching this maximum carbon dioxide hydrate conversion, the water content in carbon dioxide was observed to decrease as liquid water converted to carbon dioxide hydrate. This slow conversion to hydrate, metastability of the hydrate phase, or unexpected phases may be responsible for the large discrepancy between prior data sets for similar carbon dioxide water content measurements.

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.