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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Geologic reservoirs containing gas hydrate occur beneath permafrost environments and within marine continental slope sediments, representing a potentially vast natural gas source. Numerical simulators provide scientists and engineers with tools for understanding how production efficiency depends on the numerous, interdependent (coupled) processes associated with potential production strategies for these gas hydrate reservoirs. Confidence in the modeling and forecasting abilities of these gas hydrate reservoir simulators (GHRSs) grows with successful comparisons against laboratory and field test results, but such results are rare, particularly in natural settings. The hydrate community recognized another approach to building confidence in the GHRS: comparing simulation results between independently developed and executed computer codes on structured problems specifically tailored to the interdependent processes relevant for gas hydrate-bearing systems. The United States Department of Energy, National Energy Technology Laboratory, (DOE/NETL), sponsored the first international gas hydrate code comparison study, IGHCCS1, in the early 2000s. IGHCCS1 focused on coupled thermal and hydrologic processes associated with producing gas hydrates from geologic reservoirs via depressurization and thermal stimulation. Subsequently, GHRSs have advanced to model more complex production technologies and incorporate geomechanical processes into the existing framework of coupled thermal and hydrologic modeling. This paper contributes to the validation of these recent GHRS developments by providing results from a second GHRS code comparison study, IGHCCS2, also sponsored by DOE/NETL. IGHCCS2 includes participants from an international collection of universities, research institutes, industry, national laboratories, and national geologic surveys. Study participants developed a series of five benchmark problems principally involving gas hydrate processes with geomechanical components. The five problems range from simple geometries with analytical solutions to a representation of the world's first offshore production test of methane hydrates, which was conducted with the depressurization method off the coast of Japan. To identify strengths and limitations in the various GHRSs, study participants submitted solutions for the benchmark problems and discussed differing results via teleconferences. The GHRSs evolved over the course of IGHCCS2 as researchers modified their simulators to reflect new insights, lessons learned, and suggested performance enhancements. The five benchmark problems, final sample solutions, and lessons learned that are presented here document the study outcomes and serve as a reference guide for developing and testing gas hydrate reservoir simulators.

Geologic reservoirs containing gas hydrate occur beneath permafrost environments and within marine continental slope sediments, representing a potentially vast natural gas source. Numerical simulators provide scientists and engineers with tools for understanding how production efficiency depends on the numerous, interdependent (coupled) processes associated with potential production strategies for these gas hydrate reservoirs. Confidence in the modeling and forecasting abilities of these gas hydrate reservoir simulators (GHRSs) grows with successful comparisons against laboratory and field test results, but such results are rare, particularly in natural settings. The hydrate community recognized another approach to building confidence in the GHRS: comparing simulation results between independently developed and executed computer codes on structured problems specifically tailored to the interdependent processes relevant for gas hydrate-bearing systems. The United States Department of Energy, National Energy Technology Laboratory, (DOE/NETL), sponsored the first international gas hydrate code comparison study, IGHCCS1, in the early 2000s. IGHCCS1 focused on coupled thermal and hydrologic processes associated with producing gas hydrates from geologic reservoirs via depressurization and thermal stimulation. Subsequently, GHRSs have advanced to model more complex production technologies and incorporate geomechanical processes into the existing framework of coupled thermal and hydrologic modeling. This paper contributes to the validation of these recent GHRS developments by providing results from a second GHRS code comparison study, IGHCCS2, also sponsored by DOE/NETL. IGHCCS2 includes participants from an international collection of universities, research institutes, industry, national laboratories, and national geologic surveys. Study participants developed a series of five benchmark problems principally involving gas hydrate processes with geomechanical components. The five problems range from simple geometries with analytical solutions to a representation of the world's first offshore production test of methane hydrates, which was conducted with the depressurization method off the coast of Japan. To identify strengths and limitations in the various GHRSs, study participants submitted solutions for the benchmark problems and discussed differing results via teleconferences. The GHRSs evolved over the course of IGHCCS2 as researchers modified their simulators to reflect new insights, lessons learned, and suggested performance enhancements. The five benchmark problems, final sample solutions, and lessons learned that are presented here document the study outcomes and serve as a reference guide for developing and testing gas hydrate reservoir simulators.

Geological carbon storage in deep aquifers involves displacement of resident brine by supercritical CO2 (scCO2), which is an unstable drainage process caused by the invasion of less viscous scCO2. The unstable drainage is greatly complicated by aquifer heterogeneity and anisotropy and regarded as one of the key factors accounting for the uncertainty in storage capacity estimates. The impacts of pore-scale characteristics on the unstable drainage remain poorly understood. In this study, scCO2 drainage experiments were conducted at 40 °C and 9 MPa using a homogeneous elliptical micromodel with low or high anisotropy, a homogeneous/isotropic hexagonal micromodel, and a heterogeneous sandstone-analog micromodel. Each initially water-saturated micromodel was invaded by scCO2 at different rates with logCa (the capillary number)ranging from −7.6 to −4.4, and scCO2/water images were obtained. The measured CO2 saturations in these centimeter-scale micromodels vary considerably from 0.08 to 0.93 depending on the pore-scale characteristics and capillary number. It was also observed that scCO2 drainage follows the classic flow-regime transition from capillary fingering through crossover to viscous fingering for either of the low-anisotropy elliptical and heterogeneous micromodels, but with disparate crossover zones. The crossover zones of scCO2 saturation were then unified with the minimum scCO2 saturation occurring at logCa*=-4.0 using the complete capillary number (Ca*)that considers pore characteristics. For the hexagonal and the high-anisotropy elliptical micromodels, a monotonic increase in scCO2 saturation with increasing Ca* (without crossover)was observed. It appears that the complete capillary number is more appropriate than the classic capillary number when characterizing flow regimes and CO2 saturation in different pore networks.

Describing the pore habit of methane hydrate in sediment matrices is essential for understanding natural distribution of methane hydrate, methane trace (transport and solidification) in the hydrate stability zone, physical properties of hydrate-bearing sediments, and the associated influence on potential gas production. Pore habit visualization in natural media at pore scale even with laboratory synthesized cores has been challenging due to the similar densities of methane hydrate from pore liquid. In this work, we used phase-contrast assisted micro-CT with potassium iodine-doped brine to visualize four phases: sand particles, pore fluid, methane hydrate and methane gas. This study visualizes the pore habit of methane hydrate at various stages including during hydrate formation in excess-gas systems, its evolution after brine injection to replace pore fluid, and hydrate formation in excess-water systems. Hydrate tends to adopt round and smooth surfaces when in contact with water while exhibits relatively angular interfaces when in contact with methane gas. Hydrate formation in excess-gas systems results in a partial cementing and partial mineral-coating pore habit, while hydrate in excess-water systems develops mainly as pore-filling, and locally cementing or mineral-coating where big gas pockets exist at the initial state. Pore liquid replacement from methane gas to brine triggers a shift of hydrate pore habit towards pore-filling. Methane hydrate evolution over time produces bigger hydrate particles but with less contact area with sand particles. The effects of hydrate pore habit become less important as hydrate particle size exceeds the pore size. Additionally, hydrate formation could trap residual methane gas and brine as inclusions.

Describing the pore habit of methane hydrate in sediment matrices is essential for understanding natural distribution of methane hydrate, methane trace (transport and solidification) in the hydrate stability zone, physical properties of hydrate-bearing sediments, and the associated influence on potential gas production. Pore habit visualization in natural media at pore scale even with laboratory synthesized cores has been challenging due to the similar densities of methane hydrate from pore liquid. In this work, we used phase-contrast assisted micro-CT with potassium iodine-doped brine to visualize four phases: sand particles, pore fluid, methane hydrate and methane gas. This study visualizes the pore habit of methane hydrate at various stages including during hydrate formation in excess-gas systems, its evolution after brine injection to replace pore fluid, and hydrate formation in excess-water systems. Hydrate tends to adopt round and smooth surfaces when in contact with water while exhibits relatively angular interfaces when in contact with methane gas. Hydrate formation in excess-gas systems results in a partial cementing and partial mineral-coating pore habit, while hydrate in excess-water systems develops mainly as pore-filling, and locally cementing or mineral-coating where big gas pockets exist at the initial state. Pore liquid replacement from methane gas to brine triggers a shift of hydrate pore habit towards pore-filling. Methane hydrate evolution over time produces bigger hydrate particles but with less contact area with sand particles. The effects of hydrate pore habit become less important as hydrate particle size exceeds the pore size. Additionally, hydrate formation could trap residual methane gas and brine as inclusions.

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

We form methane hydrate by injecting methane gas into a brine-saturated, coarse-grained sample under hydrate-stable thermodynamic conditions. Hydrate forms to a saturation of 11%, which is much lower than that predicted assuming three-phase (gas-hydrate-brine) thermodynamic equilibrium (67%). During hydrate formation, there are temporary flow blockages. We interpret that a hydrate skin forms a physical barrier at the gas-brine interface. The skin fails periodically when the pressure differential exceeds the skin strength. Once the skin is present, further hydrate formation is limited by the rate that methane can diffuse through the solid skin. This process produces distinct thermodynamic states on either side of the skin that allows gas to flow through the sample. This study illuminates how gas can be transported through the hydrate stability zone and thus provides a mechanism for the formation of concentrated hydrate deposits in sand reservoirs. It also illustrates that models that assume local equilibrium at the core-scale and larger may not capture the fundamental behaviors of these gas flow and hydrate formation processes.

We form methane hydrate by injecting methane gas into a brine-saturated, coarse-grained sample under hydrate-stable thermodynamic conditions. Hydrate forms to a saturation of 11%, which is much lower than that predicted assuming three-phase (gas-hydrate-brine) thermodynamic equilibrium (67%). During hydrate formation, there are temporary flow blockages. We interpret that a hydrate skin forms a physical barrier at the gas-brine interface. The skin fails periodically when the pressure differential exceeds the skin strength. Once the skin is present, further hydrate formation is limited by the rate that methane can diffuse through the solid skin. This process produces distinct thermodynamic states on either side of the skin that allows gas to flow through the sample. This study illuminates how gas can be transported through the hydrate stability zone and thus provides a mechanism for the formation of concentrated hydrate deposits in sand reservoirs. It also illustrates that models that assume local equilibrium at the core-scale and larger may not capture the fundamental behaviors of these gas flow and hydrate formation processes.

For better understanding of frequency dependence (dispersion) of seismic wave velocities caused by stress-induced fluid flow, broadband laboratory measurements were performed on a suite of synthetic glass media containing both equant pores and thermal cracks. Complementary forced oscillation, resonant bar, and ultrasonic techniques provided access to millihertz-hertz frequencies, ~1 kHz frequency, and ~1 MHz frequency, respectively. The wave speeds or effective elastic moduli and associated dissipation were measured on samples under dry, argon- or nitrogen-saturated, and water-saturated conditions in sequence. The elastic moduli, in situ permeability, and crack porosity inferred from in situ X-ray computed tomography all attest to strong pressure-induced crack closure for differential (confining-minus-pore) pressures

For better understanding of frequency dependence (dispersion) of seismic wave velocities caused by stress-induced fluid flow, broadband laboratory measurements were performed on a suite of synthetic glass media containing both equant pores and thermal cracks. Complementary forced oscillation, resonant bar, and ultrasonic techniques provided access to millihertz-hertz frequencies, ~1 kHz frequency, and ~1 MHz frequency, respectively. The wave speeds or effective elastic moduli and associated dissipation were measured on samples under dry, argon- or nitrogen-saturated, and water-saturated conditions in sequence. The elastic moduli, in situ permeability, and crack porosity inferred from in situ X-ray computed tomography all attest to strong pressure-induced crack closure for differential (confining-minus-pore) pressures

Depressurization of hydrate-bearing sediments (HBS) can cause the movement of fine particles, and in turn, such fines migration affects fluid flow and mechanical behavior of sediments, ultimately affecting long-term hydrocarbon production and wellbore stability. This study investigated how and to what extent depressurization of HBS causes fines migration using X-ray computed tomography (CT) imaging. Methane hydrate was synthesized in sediments with 10% fines content (FC), composed of sands with silt and/or clay, and the hydrate-bearing samples were stepwisely depressurized while acquiring CT images. The CT images were analyzed to quantify the spatial changes in FC in the host sediment and thus to capture the fines migration during depressurization. It was found that the FC changes began occurring from the hydrate dissociation regions. This confirms that the multiphase flow caused by depressurization accompanies fines migration. Depressurization of HBS with a hydrate saturation of ~20–40% caused FC reduction from ~10% to ~6–9%, and the extent of fines migration differed with the particle sizes of the host sands and the types of fines. It was found that fines migration was more pronounced with coarse sands and with silty fines. Such observed level of FC reduction is estimated to increase sediment permeability by several factors based on the Kozeny-type permeability model. Our results support the notion that the extent of fines migration and its effect on fluid flow behavior need to be assessed in consideration of physical properties of host sediment and fine particles to identify optimum depressurization strategies.

Depressurization of hydrate-bearing sediments (HBS) can cause the movement of fine particles, and in turn, such fines migration affects fluid flow and mechanical behavior of sediments, ultimately affecting long-term hydrocarbon production and wellbore stability. This study investigated how and to what extent depressurization of HBS causes fines migration using X-ray computed tomography (CT) imaging. Methane hydrate was synthesized in sediments with 10% fines content (FC), composed of sands with silt and/or clay, and the hydrate-bearing samples were stepwisely depressurized while acquiring CT images. The CT images were analyzed to quantify the spatial changes in FC in the host sediment and thus to capture the fines migration during depressurization. It was found that the FC changes began occurring from the hydrate dissociation regions. This confirms that the multiphase flow caused by depressurization accompanies fines migration. Depressurization of HBS with a hydrate saturation of ~20–40% caused FC reduction from ~10% to ~6–9%, and the extent of fines migration differed with the particle sizes of the host sands and the types of fines. It was found that fines migration was more pronounced with coarse sands and with silty fines. Such observed level of FC reduction is estimated to increase sediment permeability by several factors based on the Kozeny-type permeability model. Our results support the notion that the extent of fines migration and its effect on fluid flow behavior need to be assessed in consideration of physical properties of host sediment and fine particles to identify optimum depressurization strategies.

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.

Gas flow through the vadose zone from sources at depth involves fractionation effects that can obscure the nature of transport and even the identity of the source. Transport processes are particularly complex in low permeability media but as shown in this study, can be elucidated by measuring the atmospheric noble gases. A series of laboratory column experiments was conducted to evaluate the movement of noble gas from the atmosphere into soil in the presence of a net efflux of CO2, a process that leads to fractionation of the noble gases from their atmospheric abundance ratios. The column packings were designed to simulate natural sedimentary deposition by interlayering low permeability ceramic plates and high permeability beach sand. Gas samples were collected at different depths at CO2 fluxes high enough to cause extreme fractionation of the noble gases (4He/36Ar > 20 times the air ratio). The experimental noble gas fractionation-depth profiles were in good agreement with those predicted by the dusty gas (DG) model, demonstrating the applicability of the DG model across a broad spectrum of environmental conditions. A governing equation based on the dusty gas model was developed to specifically describe noble gas fractionation at each depth that is controlled by the binary diffusion coefficient, Knudsen diffusion coefficient and the ratio of total advection flux to total flux. Finally, the governing equation was used to derive the noble gas fractionation pattern and illustrate how it is influenced by soil CO2 flux, sedimentary sequence, thickness of each sedimentary layer and each layer's physical parameters. Three potential applications of noble gas fractionation are provided: evaluating soil attributes in the path of gas flow from a source at depth to the atmosphere, testing leakage through low permeability barriers used to isolate buried waste, and tracking biological methanogenesis and methane oxidation associated with hydrocarbon degradation.

The rates of mineral dissolution reactions in porous media are difficult to predict, in part because of a lack of understanding of mineral reactive surface area in natural porous media. Common estimates of mineral reactive surface area used in reactive transport models for porous media are typically ad hoc and often based on average grain size, increased to account for surface roughness or decreased by several orders of magnitude to account for reduced surface reactivity of field as opposed to laboratory samples. In this study, accessible mineral surface areas are determined for a sample from the reservoir formation at the Nagaoka pilot CO2 injection site (Japan) using a multi-scale image analysis based on synchrotron X-ray microCT, SEM QEMSCAN, XRD, SANS, and FIB-SEM. This analysis not only accounts for accessibility of mineral surfaces to macro-pores, but also accessibility through connected micro-pores in smectite, the most abundant clay mineral in this sample. While the imaging analysis reveals that most of the micro- and macro-pores are well connected, some pore regions are unconnected and thus inaccessible to fluid flow and diffusion. To evaluate whether mineral accessible surface area accurately reflects reactive surface area a flow-through core experiment is performed and modeled at the continuum scale. The core experiment is performed under conditions replicating the pilot site and the evolution of effluent solutes in the aqueous phase is tracked. Various reactive surface area models are evaluated for their ability to capture the observed effluent chemistry, beginning with parameter values determined as a best fit to a disaggregated sediment experiment (Beckingham et al., 2016) described previously. Simulations that assume that all mineral surfaces are accessible (as in the disaggregated sediment experiment) over-predict the observed mineral reaction rates, suggesting that a reduction of RSA by a factor of 10–20 is required to match the core flood experimental data. While the fit of the effluent chemistry (and inferred mineral dissolution rates) greatly improve when the pore-accessible mineral surface areas are used, it was also necessary to include highly reactive glass phases to match the experimental observations, in agreement with conclusions from the disaggregated sediment experiment. It is hypothesized here that the 10–20 reduction in reactive surface areas based on the limited pore accessibility of reactive phases in core flood experiment may be reasonable for poorly sorted and cemented sediments like those at the Nagaoka site, although this reflects pore rather than larger scale heterogeneity.

Gas flow through the vadose zone from sources at depth involves fractionation effects that can obscure the nature of transport and even the identity of the source. Transport processes are particularly complex in low permeability media but as shown in this study, can be elucidated by measuring the atmospheric noble gases. A series of laboratory column experiments was conducted to evaluate the movement of noble gas from the atmosphere into soil in the presence of a net efflux of CO2, a process that leads to fractionation of the noble gases from their atmospheric abundance ratios. The column packings were designed to simulate natural sedimentary deposition by interlayering low permeability ceramic plates and high permeability beach sand. Gas samples were collected at different depths at CO2 fluxes high enough to cause extreme fractionation of the noble gases (4He/36Ar > 20 times the air ratio). The experimental noble gas fractionation-depth profiles were in good agreement with those predicted by the dusty gas (DG) model, demonstrating the applicability of the DG model across a broad spectrum of environmental conditions. A governing equation based on the dusty gas model was developed to specifically describe noble gas fractionation at each depth that is controlled by the binary diffusion coefficient, Knudsen diffusion coefficient and the ratio of total advection flux to total flux. Finally, the governing equation was used to derive the noble gas fractionation pattern and illustrate how it is influenced by soil CO2 flux, sedimentary sequence, thickness of each sedimentary layer and each layer's physical parameters. Three potential applications of noble gas fractionation are provided: evaluating soil attributes in the path of gas flow from a source at depth to the atmosphere, testing leakage through low permeability barriers used to isolate buried waste, and tracking biological methanogenesis and methane oxidation associated with hydrocarbon degradation.

The rates of mineral dissolution reactions in porous media are difficult to predict, in part because of a lack of understanding of mineral reactive surface area in natural porous media. Common estimates of mineral reactive surface area used in reactive transport models for porous media are typically ad hoc and often based on average grain size, increased to account for surface roughness or decreased by several orders of magnitude to account for reduced surface reactivity of field as opposed to laboratory samples. In this study, accessible mineral surface areas are determined for a sample from the reservoir formation at the Nagaoka pilot CO2 injection site (Japan) using a multi-scale image analysis based on synchrotron X-ray microCT, SEM QEMSCAN, XRD, SANS, and FIB-SEM. This analysis not only accounts for accessibility of mineral surfaces to macro-pores, but also accessibility through connected micro-pores in smectite, the most abundant clay mineral in this sample. While the imaging analysis reveals that most of the micro- and macro-pores are well connected, some pore regions are unconnected and thus inaccessible to fluid flow and diffusion. To evaluate whether mineral accessible surface area accurately reflects reactive surface area a flow-through core experiment is performed and modeled at the continuum scale. The core experiment is performed under conditions replicating the pilot site and the evolution of effluent solutes in the aqueous phase is tracked. Various reactive surface area models are evaluated for their ability to capture the observed effluent chemistry, beginning with parameter values determined as a best fit to a disaggregated sediment experiment (Beckingham et al., 2016) described previously. Simulations that assume that all mineral surfaces are accessible (as in the disaggregated sediment experiment) over-predict the observed mineral reaction rates, suggesting that a reduction of RSA by a factor of 10–20 is required to match the core flood experimental data. While the fit of the effluent chemistry (and inferred mineral dissolution rates) greatly improve when the pore-accessible mineral surface areas are used, it was also necessary to include highly reactive glass phases to match the experimental observations, in agreement with conclusions from the disaggregated sediment experiment. It is hypothesized here that the 10–20 reduction in reactive surface areas based on the limited pore accessibility of reactive phases in core flood experiment may be reasonable for poorly sorted and cemented sediments like those at the Nagaoka site, although this reflects pore rather than larger scale heterogeneity.

Fine-grained sedimentary rocks – namely mudrocks, including their laminated fissile variety — shales – make up about two thirds of all sedimentary rocks in the Earth's crust and a quarter of the continental land mass. Organic-rich shales and mudstones are the source rocks and reservoirs for conventional and unconventional hydrocarbon resources. Mudrocks are relied upon as natural barriers for geological carbon storage and nuclear waste disposal. Consideration of mudrock multi-scale physics and multi-scale spatial and temporal behavior is vital to address emergent phenomena in shale formations perturbed by engineering activities. Unique physical characteristics of shales arise as a result of their layered and highly heterogeneous and anisotropic nature, low permeability fabric, compositional complexity, and nano-scale confined chemical environments. Barriers of lexicon among geoscientists and engineers impede the development and use of conceptual models for the coupled thermal-hydraulic-mechanical-chemical-biological (THMCB) processes in mudrock formations. This manuscript reviews the THMCB process couplings, resulting emergent behavior, and key modeling approaches. We identify future research priorities, in particular fundamental knowledge gaps in understanding the phase behavior under nano-scale confinement, coupled chemo-mechanical effects on fractures, the interplay between physical and chemical processes and their rates, and issues of non-linearity and heterogeneity. We develop recommendations for future research and integrating multi-disciplinary conceptual models for the coupled multi-scale multi-physics behavior of mudrocks. Consistent conceptual models across disciplines are essential for predicting emergent processes in the subsurface, such as self-focusing of flow, time-dependent deformation (creep), fracture network development, and wellbore stability.

Fine-grained sedimentary rocks – namely mudrocks, including their laminated fissile variety — shales – make up about two thirds of all sedimentary rocks in the Earth's crust and a quarter of the continental land mass. Organic-rich shales and mudstones are the source rocks and reservoirs for conventional and unconventional hydrocarbon resources. Mudrocks are relied upon as natural barriers for geological carbon storage and nuclear waste disposal. Consideration of mudrock multi-scale physics and multi-scale spatial and temporal behavior is vital to address emergent phenomena in shale formations perturbed by engineering activities. Unique physical characteristics of shales arise as a result of their layered and highly heterogeneous and anisotropic nature, low permeability fabric, compositional complexity, and nano-scale confined chemical environments. Barriers of lexicon among geoscientists and engineers impede the development and use of conceptual models for the coupled thermal-hydraulic-mechanical-chemical-biological (THMCB) processes in mudrock formations. This manuscript reviews the THMCB process couplings, resulting emergent behavior, and key modeling approaches. We identify future research priorities, in particular fundamental knowledge gaps in understanding the phase behavior under nano-scale confinement, coupled chemo-mechanical effects on fractures, the interplay between physical and chemical processes and their rates, and issues of non-linearity and heterogeneity. We develop recommendations for future research and integrating multi-disciplinary conceptual models for the coupled multi-scale multi-physics behavior of mudrocks. Consistent conceptual models across disciplines are essential for predicting emergent processes in the subsurface, such as self-focusing of flow, time-dependent deformation (creep), fracture network development, and wellbore stability.

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.