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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

A laboratory study of cryogenic fracturing was performed to test its ability to improve oil/gas recovery from low-permeability reservoirs. Our objective is to develop well-stimulation technologies using cryogenic fluids [e.g., liquid nitrogen (LN)] to increase permeability in a large reservoir volume surrounding wells. The new technology has the potential to reduce formation damage caused by current stimulation methods and minimize or eliminate water usage. The concept of cryogenic fracturing is that a sharp thermal gradient (thermal shock) created at the surfaces of formation rocks by applying cryogenic fluid can cause strong local tensile stress and start fractures. We developed a laboratory system for cryogenic fracturing under true-triaxial loading, with LN-delivery/control and -measurement systems. The loading system simulates confining stresses by independently loading each axis up to approximately 5,000 psi on 888-in. cubes. Temperature in boreholes and at block surfaces and fluid pressure in boreholes were continuously monitored. Acoustic and pressure-decay measurements were obtained before and at various stages of stimulations. Cubic blocks (8 88-in.) of Niobrara shale, concrete, and sandstones were tested, and stress levels and anisotropies varied. Three schemes were considered: Gas fracturing without cryo-stimulation, gas fracturing after low-pressure cryogen flow-through, and gas fracturing after high-pressure cryogen flow-through. Results from pressure-decay tests show that LN stimulation clearly increases permeability, and repeated stimulations further increase the permeability. Acoustic velocities and amplitudes decreased significantly after cryo-stimulation, indicating fracture creation. In the gas fracturing without the stimulation, breakdown (complete fracturing) occurs suddenly without any initial leaking, and major fracture planes form along the plane containing principal-stress and intermediate-stress directions, as expected theoretically. However, in the gas fracturing after cryogenic stimulations, breakdown occurred gradually and with massive leaking because of thermal fractures created during stimulation. In addition, the major fracture direction does not necessarily follow the plane containing the principal-stress direction, especially at low confining-stress levels. In tests, we observed that cryogenic stimulation seems to disrupt the internal stress field. The increase in borehole temperature after stimulation affects the permeability of the specimen. When a stimulated specimen is still cold, it maintains high permeability because fractures remain open and local thermal tension is maintained near the borehole. When the rock warms back, fractures close and permeability decreases. In these tests, we have not used proppants. Overall, fractures are clearly generated by low- and high-pressure thermal shocks. The added pressure of the high-pressure thermal shocks helps to further propagate cryogenic fractures generated by thermal shock. Breakdown pressure is significantly lowered by LN stimulation, with observed breakdown-pressure reductions up to approximately 40%.

A laboratory study of cryogenic fracturing was performed to test its ability to improve oil/gas recovery from low-permeability reservoirs. Our objective is to develop well-stimulation technologies using cryogenic fluids [e.g., liquid nitrogen (LN)] to increase permeability in a large reservoir volume surrounding wells. The new technology has the potential to reduce formation damage caused by current stimulation methods and minimize or eliminate water usage. The concept of cryogenic fracturing is that a sharp thermal gradient (thermal shock) created at the surfaces of formation rocks by applying cryogenic fluid can cause strong local tensile stress and start fractures. We developed a laboratory system for cryogenic fracturing under true-triaxial loading, with LN-delivery/control and -measurement systems. The loading system simulates confining stresses by independently loading each axis up to approximately 5,000 psi on 888-in. cubes. Temperature in boreholes and at block surfaces and fluid pressure in boreholes were continuously monitored. Acoustic and pressure-decay measurements were obtained before and at various stages of stimulations. Cubic blocks (8 88-in.) of Niobrara shale, concrete, and sandstones were tested, and stress levels and anisotropies varied. Three schemes were considered: Gas fracturing without cryo-stimulation, gas fracturing after low-pressure cryogen flow-through, and gas fracturing after high-pressure cryogen flow-through. Results from pressure-decay tests show that LN stimulation clearly increases permeability, and repeated stimulations further increase the permeability. Acoustic velocities and amplitudes decreased significantly after cryo-stimulation, indicating fracture creation. In the gas fracturing without the stimulation, breakdown (complete fracturing) occurs suddenly without any initial leaking, and major fracture planes form along the plane containing principal-stress and intermediate-stress directions, as expected theoretically. However, in the gas fracturing after cryogenic stimulations, breakdown occurred gradually and with massive leaking because of thermal fractures created during stimulation. In addition, the major fracture direction does not necessarily follow the plane containing the principal-stress direction, especially at low confining-stress levels. In tests, we observed that cryogenic stimulation seems to disrupt the internal stress field. The increase in borehole temperature after stimulation affects the permeability of the specimen. When a stimulated specimen is still cold, it maintains high permeability because fractures remain open and local thermal tension is maintained near the borehole. When the rock warms back, fractures close and permeability decreases. In these tests, we have not used proppants. Overall, fractures are clearly generated by low- and high-pressure thermal shocks. The added pressure of the high-pressure thermal shocks helps to further propagate cryogenic fractures generated by thermal shock. Breakdown pressure is significantly lowered by LN stimulation, with observed breakdown-pressure reductions up to approximately 40%.

Niobrara shale cubes of 20 cm from Colorado were employed to investigate gas and supercritical CO2 injection-induced fracturing in naturally fractured caprocks of deep aquifers/depleted reservoirs and fractured shale reservoirs. Under tri-axial stresses, gas or supercritical CO2 was injected into the center of the cubes to induce fracturing. Real-time pressure and temperature, acoustic wave, pressure decay, fracture coloring, and gas fracturing were used to characterize the fracturing process and fracture morphology. Without pore pressure, CO2 injection-induced fracturing occurred and completed instantly, accompanied by an evident temperature drop. Strongly bonded fractures barely affected transverse fracture propagation, whereas weakly bonded or open fractures arrested the injected fluid first and then allowed it to generate new fractures perpendicular to the minimum horizontal stress. Breakdown pressures for cubes with preexisting fractures using gas and supercritical CO2 are much lower than both poroelastic predictions and slick-water fracturing pressure, and some are even lower than the minimum horizontal stress. This is attributed to unconformable preexisting fractures and the low viscosity of CO2. Moreover, decreasing tri-axial stress levels and increasing stress differences tend to lower the breakdown pressure. This study is instructive for understanding and tackling geomechanical issues related to CO2 geological storage and fracturing of shale reservoirs.

Niobrara shale cubes of 20 cm from Colorado were employed to investigate gas and supercritical CO2 injection-induced fracturing in naturally fractured caprocks of deep aquifers/depleted reservoirs and fractured shale reservoirs. Under tri-axial stresses, gas or supercritical CO2 was injected into the center of the cubes to induce fracturing. Real-time pressure and temperature, acoustic wave, pressure decay, fracture coloring, and gas fracturing were used to characterize the fracturing process and fracture morphology. Without pore pressure, CO2 injection-induced fracturing occurred and completed instantly, accompanied by an evident temperature drop. Strongly bonded fractures barely affected transverse fracture propagation, whereas weakly bonded or open fractures arrested the injected fluid first and then allowed it to generate new fractures perpendicular to the minimum horizontal stress. Breakdown pressures for cubes with preexisting fractures using gas and supercritical CO2 are much lower than both poroelastic predictions and slick-water fracturing pressure, and some are even lower than the minimum horizontal stress. This is attributed to unconformable preexisting fractures and the low viscosity of CO2. Moreover, decreasing tri-axial stress levels and increasing stress differences tend to lower the breakdown pressure. This study is instructive for understanding and tackling geomechanical issues related to CO2 geological storage and fracturing of shale reservoirs.

Leakage risk assessment is an inevitable procedure in permanent sequestration and storage of CO2 in deep saline aquifers and depleted oil and gas reservoirs, where the integrity of caprock is most critical. Low porosity and low permeability concrete cubes were employed as caprock analogs to investigate the supercritical CO2 injection-induced fracturing processes under true tri-axial stress conditions. A systematic experimental procedure, consisting of active acoustic emission measurement, pressure decay, injection pressure and temperature monitoring, fracture coloring, and gas fracturing, is formulated to qualitatively and quantitatively characterize the injection-induced fracturing processes and fracture morphology. Occurrence of injection-induced fracturing can be directly identified from peaks of borehole pressure profiles as well as sharp drops on temperature profiles due to CO2 expansion, but generally there was no fracture propagation plateau appearing for the 20 cm rock cubes with zero pore pressure. Acoustic wave signatures, including both waveform change and arrival time delay, can effectively capture the extension of induced fractures inside the opaque rock. Initiation and propagation of supercritical CO2 injection-induced fractures are highly dominated by the tri-axial stresses, following the general trend of continuum mechanics at high stress levels with large stress difference. As the stress difference decreases, induced fractures branch off in relatively arbitrary directions. For these supercritical CO2 injection-induced fracturing experiments, poroelastic mechanics model makes a decent fit with the measured breakdown pressure of the caprock analogs. Findings in this study are valuable for risk analysis and operation optimization of geological CO2 sequestration and storage as well as for CO2 fracturing design and implementation in shale and tight reservoirs.

The concept of cryogenic fracturing is that a sharp thermal gradient developed by applying a cryogenic fluid on a rock surface causes a strong local tensile stress that initiates fractures. Prior field tests suggest that field application with special equipment rated for cryogenic temperatures may deliver potential benefits. The tests did not, however, identify the fracture mechanisms at work in downhole conditions. In this study, we present our laboratory designs and procedures developed for studying cryogenic fracturing mechanisms in a well environment, and examine typical data indicative of the performance of the system. The experimental apparatus and procedures were specifically designed to conduct cryogenic fracturing tests in specimens under confining stress, with integrated cryogen transport, measurements, and fracture characterization. A true-triaxial loading system was built to simulate reservoir stress levels and anisotropic stress application, and was designed to avoid thermal stresses in tests involving cryogen by arranging discrete components in an open chamber. To maximize thermal shock on the wellbore surface, tubing and wellhead are configured so that liquid nitrogen enters the borehole and flow out after contributing to thermal shock. The temperature at boreholes and along flow lines, borehole pressure, and liquid nitrogen consumption are monitored throughout treatments. Acoustic transmission and pressure-decay measurements are used to characterize fractures before and after the experiments. Breakdown tests performed on un-stimulated specimens and stimulated specimens compare breakdown pressures of the two groups and thus evaluate the performance of thermal fracturing. The laboratory design was able to effectively apply cryogenic stimulations to laboratory rock specimens. The characterization methods were able to capture fracture creation and rock property changes due to cryogenic fracturing.

The concept of cryogenic fracturing is that a sharp thermal gradient developed by applying a cryogenic fluid on a rock surface causes a strong local tensile stress that initiates fractures. Prior field tests suggest that field application with special equipment rated for cryogenic temperatures may deliver potential benefits. The tests did not, however, identify the fracture mechanisms at work in downhole conditions. In this study, we present our laboratory designs and procedures developed for studying cryogenic fracturing mechanisms in a well environment, and examine typical data indicative of the performance of the system. The experimental apparatus and procedures were specifically designed to conduct cryogenic fracturing tests in specimens under confining stress, with integrated cryogen transport, measurements, and fracture characterization. A true-triaxial loading system was built to simulate reservoir stress levels and anisotropic stress application, and was designed to avoid thermal stresses in tests involving cryogen by arranging discrete components in an open chamber. To maximize thermal shock on the wellbore surface, tubing and wellhead are configured so that liquid nitrogen enters the borehole and flow out after contributing to thermal shock. The temperature at boreholes and along flow lines, borehole pressure, and liquid nitrogen consumption are monitored throughout treatments. Acoustic transmission and pressure-decay measurements are used to characterize fractures before and after the experiments. Breakdown tests performed on un-stimulated specimens and stimulated specimens compare breakdown pressures of the two groups and thus evaluate the performance of thermal fracturing. The laboratory design was able to effectively apply cryogenic stimulations to laboratory rock specimens. The characterization methods were able to capture fracture creation and rock property changes due to cryogenic fracturing.

CO2 flooding is a proven enhanced oil recovery (EOR) technique and is also considered as a potential method for CO2 sequestration. Despite having successful field trials on CO2 EOR, the effects of reservoir wettability and permeability heterogeneity on the efficiency of miscible CO2 flooding are not well understood. In this work, laboratory investigations have been carried out to evaluate the influence of these properties on the miscible CO2 EOR performance. The wettability of hydrophilic Berea core samples was altered to be oil-wet by vacuum saturation of the clean and dry core samples with n-hexadecane. The permeability heterogeneity was obtained by combining two half pieces of axially split water-wet core samples of different permeabilities. Core flooding experiments were conducted for n-hexadecane - synthetic brine - CO2 systems at 1400 psig backpressure to achieve minimum miscibility pressure (MMP) of CO2 in n-hexadecane at the test temperature (24 ± 1 °C). It was found that wettability strongly influences CO2 EOR. For the alternate cases of previously brine flooded (to remaining oil saturation) oil-wet and water-wet core samples, five pore volumes (PVs) of CO2 recovered 100% and only 43% of remaining oil in place (ROIP) respectively. Three PVs of CO2 could recover only about 0-5% ROIP from the split core samples. The mechanisms underlying these results are discussed. This study sheds light on the significant influence of reservoir wettability and permeability heterogeneity on the performance of miscible CO2 EOR.