Fracture properties are important in carbonate reservoir characterization, as they are responsible for a large part of the fluid transfer properties at all scales. It is especially true in tight rocks where the matrix transfer properties only slightly contribute to the fluid flow. Open fractures are known to strongly affect seismic velocities, amplitudes and anisotropy. Here, we explore the impact of fracture evolution on the geophysical signature and directional Vp anisotropy of fractured carbonates through diagenesis. For that purpose, we studied a meter-scale, parallelepiped quarry block of limestone using a detailed structural and diagenetic characterization, and numerous Vp measurements. The block is affected by two en-échelon fracture clusters, both being formed in opening mode (mode 1) and cemented, but only one being reactivated in shear. We compared the diagenetic evolution of the fractures, which are almost all 100% filled with successive calcite cements, with the P-wave velocities measured across this meter-scale block of carbonate, which recorded the tectonic and diagenetic changes of a South Provence sedimentary basin. We found that a directional Vp anisotropy magnitude as high as 8–16% correlates with the reactivated fractures’ cluster dip angle, which is explained by the complex filling sequence and softer material present inside the fractures that have been reactivated during the basin's tectonic inversion. We show that although a late karstification phase preferentially affected these reactivated fractures, it only amplified the pre-existing anisotropy due to tectonic shear. We conclude that Vp anisotropy measurements may help to identify the fracture sealing/opening processes associated with polyphased tectonic history, the anisotropy being independent of the current stress-state. This case shows that velocity anisotropies induced by fractures resulted here from a cause that is different from how these features have often been interpreted: selective reactivation of sealed fractures clusters rather than direction of currently open ones.

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.

Upper Barremian – Lower Aptian inner platform “Urgonian” limestones in the Mont de Vaucluse region, SE France, consist of alternating metre-scale microporous and tight intervals. This paper focuses on the influence of structural deformation on the reservoir properties of the Urgonian limestone succession in a study area near the town of Rustrel. Petrographic, petrophysical and structural data were recovered from five fully-cored boreholes, from the walls of a 100 m long underground tunnel, and from a 50 m long transect at a nearby outcrop. The data allowed reservoir property variations in the Urgonian limestones to be studied from core to reservoir scale. Eleven Reservoir Rock Types (RRTs) were identified based on petrographic features (texture, grain size), reservoir properties (porosity, permeability), and the frequency of structural discontinuities such as fractures, faults and stylolites. Tight and microporous reservoir rock types were distinguished. Tight reservoir rock types were characterised by early cementation of intergranular pore spaces and by the presence of frequent structural discontinuities. By contrast microporous reservoir rock types contained preserved intragranular microporosity and matrix permeability, but had very few structural discontinuities. Observed vertical alternations of microporous and tight rock types are interpreted to have been controlled by the early diagenesis of the Urgonian carbonates. Deformation associated with regional-scale tectonic phases, including Albian – Cenomanian “Durancian” uplift (∼105 to 96 Ma) and Pyrenean compression (∼55 to 25 Ma), resulted in the modification of the initial petrophysical properties of the Urgonian limestones. An early diagenetic imprint conditioned both the intensity of structural deformations and the associated circulations of diagenetic and meteoric fluids. Evolution of the Reservoir Rock Types is therefore linked both to the depositional conditions and to subsequent phases of structural deformation.

Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far-reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture-dependent quantities of interest such as CO2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field- and laboratory-scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory-scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics-based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field-scale fracture systems. Finally, we review the use of machine learning-based emulators to rapidly investigate different fracture property scenarios and accelerate physics-based models by orders of magnitude to enable uncertainty quantification and near real-time analysis.

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

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments in the crystalline rock at the site to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, laboratory measurements and experiments, and we conducted an extended analysis of the local stress state using previously collected data. Some of these activities are still ongoing. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) experiments meet objectives in SubTER’s “stress” pillar and the “new subsurface signals” pillar. The kISMET site was established in the West Access Drift of SURF 4850 ft (1478 m) below ground (on the 4850L) in phyllite of the Precambrian Poorman Formation. We drilled and cored five near-vertical boreholes in a line on 3 m spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume at ~1528 m (5013 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. The load-displacement record on the core suggests that the elastic stiffness is anisotropic. Tensile strength ranges between 3‒7.5 MPa and 5‒12 MPa. Permeability measurements are planned, as are two types of laboratory miniature hydraulic fracturing experiments to investigate the importance of rock fabric (anisotropy and heterogeneity) on near-borehole hydraulic fracture generation. Pre-fracturing numerical simulations with INL’s FALCON discrete element code predicted a fracture radius of 1.2 m for a corresponding injection volume of 1.2 L for the planned fractures, and negligible microseismicity. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) pointing approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stress are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Monitoring using ERT and CASSM in the four monitoring boreholes, and passive seismic accelerometer-based measurements in the West Access Drift, was carried out during the generation of a larger fracture (so-called stimulation test) at a depth of 40 m below the invert. ERT was not able to detect the fracture created, nor were the accelerometers in the drift, but microseismicity was detected for first (deepest) hydraulic-fracturing stress measurement. The CASSM data have not yet been analyzed. Analytical solutions suggest fracture radius of the large fracture (stimulation test) was more than 6 m, depending on the unknown amount of leak-off. The kISMET results for stress state are consistent with large-scale mid-continent estimates of stress. Currently we are using the orientation of the stress field we determined to interpret a large number of borehole breakouts recorded in nearby boreholes at SURF to generate a more complete picture of the stress field and its variations at SURF. The efforts on the project have prompted a host of additional follow-on studies that we recommend be carried out at the kISMET site.

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.

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