Search

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.

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 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.

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.

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.

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.

We investigate the influence of foliation orientation and fine-scale folding on the static and dynamic elastic properties and unconfined strength of the Poorman schist. Measurements from triaxial and uniaxial laboratory experiments reveal a significant amount of variability in the static and dynamic Young’s modulus depending on the sample orientation relative to the foliation plane. Dynamic P-wave modulus and S-wave modulus are stiffer in the direction parallel to the foliation plane as expected for transversely isotropic mediums with average Thomsen parameters values 0.133 and 0.119 for epsilon and gamma, respectively. Static Young’s modulus varies significantly between 21 and 117 GPa, and a peculiar trend is observed where some foliated sample groups show an anomalous decrease in the static Young’s modulus when the symmetry axis (x3-axis) is oriented obliquely to the direction of loading. Utilizing stress and strain relationships for transversely isotropic medium, we derive the analytical expression for Young’s modulus as a function of the elastic moduli E1, E3, ν31, and G13 and sample orientation to fit the static Young’s modulus measurements. Regression of the equation to the Young’s modulus data reveals that the decrease in static Young’s modulus at oblique symmetry axis orientations is directly influenced by a low shear modulus, G13, which we attribute to shear sliding along foliation planes during static deformation that occurs as soon as the foliation is subject to shear stress. We argue that such difference between dynamic and static anisotropy is a characteristic of near-zero porosity anisotropic rocks. The uniaxial compressive strength also shows significant variability ranging from 21.9 to 194.6 MPa across the five sample locations and is the lowest when the symmetry axis is oriented 45° or 60° from the direction of loading, also a result of shear sliding along foliation planes during static deformation.

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). The EGS Collab project is a collaborative research environment studying stimulation of crystalline rock at the 10-meter scale. High-quality characterization and monitoring data are collected during stimulation and flow tests to allow comparison to numerical coupled process models in an effort to build confidence in the codes and improve modeling techniques. The Experiment 1 test bed, located at the Sanford Underground Research Facility (SURF) in Lead, SD at a depth of approximately 1.5 km is being used to examine hydraulic fracturing. The testbed was characterized using numerous field-based geophysical and geologic techniques and laboratory testing, and a well-instrumented test bed was created to allow us to carefully monitor fracture stimulation events and flow tests. A number of hydraulic stimulation tests at several locations in one well were performed, creating new fractures that connect to existing fractures between the injection and production boreholes. Long-term ambient and chilled water injection tests have been performed as an analog to EGS, and system changes resulting from these water injections have been monitored using geophysical monitoring, flow, temperature, and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the scientific findings to date from Experiment 1. We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2).

Copyright 2019 ARMA, American Rock Mechanics Association. The EGS (Enhanced Geothermal Systems) Collab project is performing stimulation and flow experiments in highly-monitored and well-characterized intermediate-scale (approximately10 to 20 meter) field test beds at a depth of approximately 1,500 meters in the Sanford Underground Research Facility (SURF) in the Black Hills of South Dakota. Our fracture stimulation and interwell flow tests are performed to better understand processes that control formation of effective subsurface heat exchangers that are critical to the development and success of EGS. Different EGS Collab stimulations will be performed under dissimilar stress conditions to produce data for model comparisons that better differentiate stimulation mechanisms and the evolution of permeability enhancement in crystalline rock. EGS Collab experiments provide a means of testing tools, concepts, and strategies that could later be employed under geothermal reservoir conditions at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) and other enhanced geothermal systems. Key to the project is using numerical simulations in the experiment design and interpretation of

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful presite selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory‐scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time‐lapse 3‐dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear‐stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation. Plain Language Summary: Time‐lapse electrical geophysical sensing is used to image 3D changes in rock stress generated by an isolated and pressurized interval of a borehole in a deep, dense, fractured rock formation. Key Points: Remotely monitoring stress is challenging but important for relating geomechanical behavior to flow pathways during energy productionBulk electrical conductivity is sensitive to stress in crystalline rockTime‐lapse electrical resistivity tomography can be used to remotely monitor 3D changes in effective stress [ABSTRACT FROM AUTHOR]

With radar sounders, coherent backscattering simulations from global planetary digital elevation models (DEMs) typically display a deficit in diffuse clutter, which is mainly due to the implicit assumption that roughness at scales below the resolution of the DEM is absent. Indeed, while polynomial approximations of the phase evolution across the facet allow for fast and mathematically rigorous simulators, the coarse resolution of these planetary DEMs leads to a potentially significant portion of the backscattering response being neglected. In this paper, we derive the analytical phase response of a rough rectangular facet characterized by Gaussian roughness and a Gaussian isotropic correlation function under the linear phase approximation. Formulae for the coherent and incoherent power scattered by such an object are obtained for arbitrary bistatic scattering angles. Validation is done both in isolation and after inclusion in different Stratton‐Chu simulators. In order to illustrate the different uses of such a formulation, we reproduce two lunar radargrams acquired by the Lunar Radar Sounder instrument with a Stratton‐Chu simulator incorporating the proposed rough facet phase integral, and we show that the original radargrams are significantly better‐reproduced than with state‐of‐the‐art methods, at a similar computational cost. We also show how the rough facet integral formulation can be used in isolation to better characterize subglacial water bodies on Earth. Key Points: Planetary digital elevation models are often of coarse resolution and depict a surface that is smooth at scales below that resolutionPolynomial phase approximations can be used to simulate radar scattering rigorously but they overestimate the coherence of reflected signalsWe analytically derive the linear phase approximation formula on a rough rectangular facet, leading to much better clutter simulations [ABSTRACT FROM AUTHOR]