Thermal effects on geologic carbon storage

The final publication is available at Springer via http://dx.doi.org/10.1016/j.earscirev.2016.12.011One of the most promising ways to significantly reduce greenhouse gases emissions, while carbon-free energy sources are developed, is Carbon Capture and Storage (CCS). Non-isothermal effects play a major role in all stages of CCS. In this paper, we review the literature on thermal effects related to CCS, which is receiving an increasing interest as a result of the awareness that the comprehension of non-isothermal processes is crucial for a successful deployment of CCS projects. We start by reviewing CO2 transport, which connects the regions where CO2 is captured with suitable geostorage sites. The optimal conditions for CO2 transport, both onshore (through pipelines) and offshore (through pipelines or ships), are such that CO2 stays in liquid state. To minimize costs, CO2 should ideally be injected at the wellhead in similar pressure and temperature conditions as it is delivered by transport. To optimize the injection conditions, coupled wellbore and reservoir simulators that solve the strongly non-linear problem of CO2 pressure, temperature and density within the wellbore and non-isothermal two-phase flow within the storage formation have been developed. CO2 in its way down the injection well heats up due to compression and friction at a lower rate than the geothermal gradient, and thus, reaches the storage formation at a lower temperature than that of the rock. Inside the storage formation, CO2 injection induces temperature changes due to the advection of the cool injected CO2, the Joule-Thomson cooling effect, endothermic water vaporization and exothermic CO2 dissolution. These thermal effects lead to thermo-hydro-mechanical-chemical coupled processes with non-trivial interpretations. These coupled processes also play a relevant role in “Utilization” options that may provide an added value to the injected CO2, such as Enhanced Oil Recovery (EOR), Enhanced Coal Bed Methane (ECBM) and geothermal energy extraction combined with CO2 storage. If the injected CO2 leaks through faults, the caprock or wellbores, strong cooling will occur due to the expansion of CO2 as pressure decreases with depth. Finally, we conclude by identifying research gaps and challenges of thermal effects related to CCS.Peer ReviewedPostprint (author's final draft

Development of a Numerical Simulator for Analyzing the Geomechanical Performance of Hydrate-Bearing Sediments

In this paper, we describe the development and application of a numerical simulator that analyzes the geomechanical performance of hydrate-bearing sediments, which may become an important future energy supply. The simulator is developed by coupling a robust numerical simulator of coupled fluid flow, hydrate thermodynamics, and phase behavior in geologic media (TOUGH+HYDRATE) with an established geomechanical code (FLAC3D). We demonstrate the current simulator capabilities and applicability for two examples of geomechanical responses of hydrate bearing sediments during production-induced hydrate dissociation. In these applications, the coupled geomechanical behavior within hydrate-bearing seducements are considered through a Mohr-Coulomb constitutive model, corrected for changes in pore-filling hydrate and ice content, based on laboratory data. The results demonstrate how depressurization-based gas production from oceanic hydrate deposits may lead to severe geomechanical problems unless care is taken in designing the production scheme. We conclude that the coupled simulator can be used to design production strategies for optimizing production, while avoiding damaging geomechanical problems

Analysis of Injection-Induced Micro-Earthquakes in a Geothermal Steam Reservoir, The Geysers Geothermal Field, California

In this study we analyze relative contributions to the cause and mechanism of injection-induced micro-earthquakes (MEQs) at The Geysers geothermal field, California. We estimated the potential for inducing seismicity by coupled thermal-hydrological-mechanical analysis of the geothermal steam production and cold water injection to calculate changes in stress (in time and space) and investigated if those changes could induce a rock mechanical failure and associated MEQs. An important aspect of the analysis is the concept of a rock mass that is critically stressed for shear failure. This means that shear stress in the region is near the rock-mass frictional strength, and therefore very small perturbations of the stress field can trigger an MEQ. Our analysis shows that the most important cause for injection-induced MEQs at The Geysers is cooling and associated thermal-elastic shrinkage of the rock around the injected fluid that changes the stress state in such a way that mechanical failure and seismicity can be induced. Specifically, the cooling shrinkage results in unloading and associated loss of shear strength in critically shear-stressed fractures, which are then reactivated. Thus, our analysis shows that cooling-induced shear slip along fractures is the dominant mechanism of injection-induced MEQs at The Geysers

A New Parameter to Assess Hydromechanical Effect in Single-hole Hydraulic Testing and Grouting

Grouting or filling of the open voids in fractured rock is done by introducing a fluid, a grout, through boreholes under pressure. The grout may be either a Newtonian fluid or a Bingham fluid. The penetration of the grout and the resulting pressure profile may give rise to hydromechanical effects, which depends on factors such as the fracture aperture, pressure at the borehole and the rheological properties of the grout. In this paper, we postulate that a new parameter, {angstrom}, which is the integral of the fluid pressure change in the fracture plane, is an appropriate measure to describe the change in fracture aperture volume due to a change in effective stress. In many cases, analytic expressions are available to calculate pressure profiles for relevant input data and the {angstrom} parameter. The approach is verified against a fully coupled hydromechanical simulator for the case of a Newtonian fluid. Results of the verification exercise show that the new approach is reasonable and that the {angstrom}-parameter is a good measure for the fracture volume change: i.e., the larger the {angstrom}-parameter, the larger the fracture volume change, in an almost linear fashion. To demonstrate the application of the approach, short duration hydraulic tests and constant pressure grouting are studied. Concluded is that using analytic expressions for penetration lengths and pressure profiles to calculate the {angstrom} parameter provides a possibility to describe a complex situation and compare, discuss and weigh the impact of hydromechanical couplings for different alternatives. Further, the analyses identify an effect of high-pressure grouting, where uncontrolled grouting of larger fractures and insufficient (or less-than-expected) sealing of finer fractures is a potential result

Uncertainty in the Maximum Principal Stress Estimated from Hydraulic Fracturing Measurements Due to the Presence of the Induced Fracture

Abstract The classical theory for hydraulic fracturing stress measurements assumes an ideal case with a linear elastic, homogenous, and isotropic medium; and a fracture that reopens distinctly when the minimum tangential borehole stress is exceeded. The induced fracture disturbs this ideal picture in several aspects, which are important for the evaluation of the maximum horizontal principal stress using the fracture reopening pressure. This disturbance can be attributed to the fracture normal stiffness and the initial hydraulic fracture permeability. In this paper, the hydraulic fracturing reopening test is studied by coupled hydromechanical modeling that takes into account an induced fracture that is incompletely closed. The result shows that with realistic equipment compliance, the apparent fracture reopening evaluated from the well-pressure is close to the magnitude of the minimum horizontal principal stress with little or no correlation to the maximum horizontal principal stress. This observation suggests that determination of maximum principal stress by hydraulic fracturing using the reopening pressure is very uncertain

On the Relationship between Stress and Elastic Strain for Porous and Fractured Rock

Modeling the mechanical deformations of porous and fractured rocks requires a stress-strain relationship. Experience with inherently heterogeneous earth materials suggests that different varieties of Hook's law should be applied within regions of the rock having significantly different stress-strain behavior, e.g., such as solid phase and various void geometries. We apply this idea by dividing a rock body conceptually into two distinct parts. The natural strain (volume change divided by rock volume at the current stress state), rather than the engineering strain (volume change divided by the unstressed rock volume), should be used in Hooke's law for accurate modeling of the elastic deformation of that part of the pore volume subject to a relatively large degree of relative deformation (i.e., cracks or fractures). This approach permits the derivation of constitutive relations between stress and a variety of mechanical and/or hydraulic rock properties. We show that the theoretical predictions of this method are generally consistent with empirical expressions (from field data) and also laboratory rock experimental data

The Influence of an Interlayer on Dual Hydraulic Fractures Propagation

Multi-cluster hydraulic fracturing of long-range horizontal wells is an approach for enhancing the productivity of low-permeability shale reservoirs. In this study, RFPA-Petrol (rock failure process analysis on petroleum problems) is applied for modeling hydraulic fracture propagation in multilayered formations. RFPA-Petrol based on coupled hydraulic-mechanical-damage (HMD) modeling was first tested by modeling a laboratory scale experiment on a physical (cement) model with a single completion. The modeling demonstrated the capability of RFPA-Petrol for simulating hydraulic fracture propagation. Then, we used RFPA-Petrol to investigate how the difference in material properties between oil-bearing layers and interlayers and the fracturing fluid properties influence the propagation of dual fractures in multilayered laboratory-scale models. In this case, the models with geological discontinuities in the vertical direction are strongly heterogeneous and RFPA-Petrol simulations successfully modeled the fracture configurations

Implementation of the Barcelona Basic Model into TOUGH–FLAC for simulations of the geomechanical behavior of unsaturated soils

This paper presents the implementation of the Barcelona Basic Model (BBM) into the TOUGH-FLAC simulator analyzing the geomechanical behavior of unsaturated soils. We implemented the BBM into TOUGH-FLAC by (1) extending an existing FLAC{sup 3D} module for the Modified Cam-Clay (MCC) model in FLAC{sup 3D} and (2) adding computational routines for suction-dependent strain and net stress (i.e., total stress minus gas pressure) for unsaturated soils. We implemented a thermo-elasto-plastic version of the BBM, wherein the soil strength depends on both suction and temperature. The implementation of the BBM into TOUGH-FLAC was verified and tested against several published numerical model simulations and laboratory experiments involving the coupled thermal-hydrological-mechanical (THM) behavior of unsaturated soils. The simulation tests included modeling the mechanical behavior of bentonite-sand mixtures, which are being considered as back-fill and buffer materials for geological disposal of spent nuclear fuel. We also tested and demonstrated the use of the BBM and TOUGH-FLAC for a problem involving the coupled THM processes within a bentonite-backfilled nuclear waste emplacement tunnel. The simulation results indicated complex geomechanical behavior of the bentonite backfill, including a nonuniform distribution of buffer porosity and density that could not be captured in an alternative, simplified, linear-elastic swelling model. As a result of the work presented in this paper, TOUGH-FLAC with BBM is now fully operational and ready to be applied to problems associated with nuclear waste disposal in bentonite-backfilled tunnels, as well as other scientific and engineering problems related to the mechanical behavior of unsaturated soils

Geomechanical Performance of Hydrate-Bearing Sediment in Offshore Environments

The objective of this multi-year, multi-institutional research project was to develop the knowledge base and quantitative predictive capability for the description of geomechanical performance of hydrate-bearing sediments (hereafter referred to as HBS) in oceanic environments. The focus was on the determination of the envelope of hydrate stability under conditions typical of those related to the construction and operation of offshore platforms. We have developed a robust numerical simulator of hydrate behavior in geologic media by coupling a reservoir model with a commercial geomechanical code. We also investigated the geomechanical behavior of oceanic HBS using pore-scale models (conceptual and mathematical) of fluid flow, stress analysis, and damage propagation. The objective of the UC Berkeley work was to develop a grain-scale model of hydrate-bearing sediments. Hydrate dissociation alters the strength of HBS. In particular, transformation of hydrate clusters into gas and liquid water weakens the skeleton and, simultaneously, reduces the effective stress by increasing the pore pressure. The large-scale objective of the study is evaluation of geomechanical stability of offshore oil and gas production infrastructure. At Lawrence Berkeley National Laboratory (LBNL), we have developed the numerical model TOUGH + Hydrate + FLAC3D to evaluate how the formation and disassociation of hydrates in seafloor sediments affects seafloor stability. Several technical papers were published using results from this model. LBNL also developed laboratory equipment and methods to produce realistic laboratory samples of sediments containing gas hydrates so that mechanical properties could be measured in the laboratory. These properties are required to run TOUGH + Hydrate + FLAC3D to evaluate seafloor stability issues. At Texas A&M University we performed a detailed literature review to determine what gas hydrate formation properties had been measured and reported in the literature. We then used TOUGH + Hydrate to simulate the observed gas production and reservoir pressure field data at Messoyakha. We simulated various scenarios that help to explain the field behavior. We have evaluated the effect of reservoir parameters on gas recovery from hydrates. Our work should be beneficial to others who are investigating how to produce gas from a hydrate capped gas reservoir. The results also can be used to better evaluate the process of producing gas from offshore hydrates. The Schlumberger PETREL model is used in industry to the description of geologic horizons and the special distribution of properties. An interface between FLAC3D and Petrel was built by Schlumberger to allow for efficient data entry into TOUGH + Hydrate + FLAC3D

Shear-slip analysis in multiphase fluid-flow reservoir engineeringap plications using TOUGH-FLAC

This paper describes and demonstrates the use of the coupledTOUGH-FLAC simulator for geomechanical shear-slip (failure) analysis inmultiphase fluid-flow reservoir-engineering applications. Two approachesfor analyzing shear-slip are described, one using continuum stress-strainanalysis and another using discrete fault analysis. The use of shear-slipanalysis in TOUGH-FLAC is demonstrated on application examples related toCO2 sequestration and geothermal energy extraction. In the case of CO2sequestration, the shear-slip analysis is used to evaluate maximumsustainable CO2-injection pressure under increasing reservoir pressure,whereas in the case of geothermal energy extraction, the shear-slipanalysis is used to study induced seismicity during steam productionunder decreasing reservoir pressure and temperature