skip to main content

Title: Computational Modeling of Creep Behavior in Shales Induced by Fluid-Rock Interaction
ABSTRACT:

The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones of more » reacted shale.

1. INTRODUCTION

Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.

« less
Authors:
;
Award ID(s):
2045242
Publication Date:
NSF-PAR ID:
10395568
Journal Name:
56th U.S. Rock Mechanics/Geomechanics Symposium, Santa Fe, New Mexico, USA, June 2022
Sponsoring Org:
National Science Foundation
More Like this
  1. Geological carbon sequestration provides permanent CO2 storage to mitigate the current high concentration of CO2 in the atmosphere. CO2 mineralization in basalts has been proven to be one of the most secure storage options. For successful implementation and future improvements of this technology, the time-dependent deformation behavior of reservoir rocks in presence of reactive fluids needs to be studied in detail. We conducted load stepping creep experiments on basalts from the CarbFix site (Iceland) under several pore fluid conditions (dry, H2O-saturated and H2O+CO2-saturated) at temperature, T≈80 °C and effective pressure, Peff = 50 MPa, during which we collected mechanical, acoustic and pore fluid chemistry data. We observed transient creep at stresses as low as 11% of the ultimate failure strength, well below the stress level at the onset of bulk dilatancy. Acoustic emissions (AEs) correlated strongly with strain accumulation, indicating that the creep deformation was a brittle process in agreement with microstructural observations. The rate and magnitude of AEs were higher in fluid-saturated experiments than in dry conditions. We infer that the predominant mechanism governing creep deformation is time- and stress-dependent sub-critical dilatant cracking. Our results suggest that the presence of aqueous fluids exerts first order control on creep deformation of basaltic rocks, whilemore »the composition of the fluids plays only a secondary role under the studied conditions.« less
  2. Abstract The vagina is a viscoelastic fibromuscular organ that provides support to the pelvic organs. The viscoelastic properties of the vagina are understudied but may be critical for pelvic stability. Most studies evaluate vaginal viscoelasticity under a single uniaxial load; however, the vagina is subjected to dynamic multiaxial loading in the body. It is unknown how varied multiaxial loading conditions affect vaginal viscoelastic behavior and which microstructural processes dictate the viscoelastic response. Therefore, the objective was to develop methods using extension-inflation protocols to quantify vaginal viscoelastic creep under various circumferential and axial loads. Then, the protocol was applied to quantify vaginal creep and collagen microstructure in the fibulin-5 wildtype and haploinsufficient vaginas. To evaluate pressure-dependent creep, the fibulin-5 wildtype and haploinsufficient vaginas (n = 7/genotype) were subjected to various constant pressures at the physiologic length for 100 s. For axial length-dependent creep, the vaginas (n = 7/genotype) were extended to various fixed axial lengths then subjected to the mean in vivo pressure for 100 s. Second-harmonic generation imaging was performed to quantify collagen fiber organization and undulation (n = 3/genotype). Increased pressure significantly increased creep strain in the wildtype, but not the haploinsufficient vagina. The axial length did not significantly affect the creep rate or strain in bothmore »genotypes. Collagen undulation varied through the depth of the subepithelium but not between genotypes. These findings suggest that the creep response to loading may vary with biological processes and pathologies, therefore, evaluating vaginal creep under various circumferential loads may be important to understand vaginal function.« less
  3. ABSTRACT:

    Due to rock mass being commonly subjected to compressive or shear loading, the mode II fracture toughness is an important material parameter for rocks. Fracturing in rocks is governed by the behavior of a nonlinear region surrounding the crack tip called the fracture process zone (FPZ). However, the characteristics of mode II fracture are still determined based on the linear elastic fracture mechanics (LEFM), which assumes that a pure mode II loading results in a pure mode II fracture. In this study, the FPZ development in Barre granite specimens under mode II loading was investigated using the short beam compression (SBC) test. Additionally, the influence of lateral confinement on various characteristics of mode II fracture was studied. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique (2D-DIC) to identify fracture mode and characterize the FPZ evolution in Barre granite specimens. The 2D-DIC analysis showed a dominant mixed-mode I/II fracture in the ligament between two notches, irrespective of confinement level on the SBC specimens. The influence of confinement on the SBC specimens was assessed by analyzing the evolution of crack displacement and changes in value of mode II fracture toughness. Larger levels ofmore »damage in confined specimens were observed prior to the failure than the unconfined specimens, indicating an increase in the fracture resistance and therefore mode II fracture toughness with the confining stress.

    1. INTRODUCTION

    The fracturing in laboratory-scale rock specimens is often characterized by the deformation of the inelastic region surrounding the crack tips, also known as the fracture process zone (FPZ) (Backers et al., 2005; Ghamgosar and Erarslan, 2016). While the influence of the FPZ on mode I fracture in rocks has been extensively investigated, there are limited studies on FPZ development in rocks under pure mode II loading (Ji et al., 2016; Lin et al., 2020; Garg et al., 2021; Li et al., 2021).

    « less
  4. We observed and modeled the elastic, inelastic and time-dependent viscous properties of damaged Berea Sandstone samples to investigate the impact of damage on the rheological properties of rocks. Cylindrical Berea Sandstone plugs were prepared both parallel and perpendicular to bedding. We impacted the samples with Split Hopkinson Pressure Bar to pervasively fracture the specimens at different strain rates. Longitudinal mode-I fractures are dominant in specimens impacted at relatively low strain rates (about 130 /s), whereas shear fractures also form in specimens deformed at high strain rates (up to 250 /s). The damaged rocks were subjected to multiple steps of differential stress loading and hold stages under 15 MPa confining pressure. A key observation is that higher damaged specimens showed greater axial and volumetric creep strain deformation during loading and hold stages. Poisson ratio also increase with increasing damage. We modeled the volumetric strain of the sandstone specimens using a Perzyna viscoplasticity law that employs the Modified Cam Clay model as the yield criterion (Haghighat et al. 2020). We deduced that fractured rocks undergo substantial bulk time-dependent deformation due to volumetric compaction and fracture closure. Damage increase results in decrease of the effective viscosity of the material.
  5. Abstract. Geological carbon sequestration provides permanentCO2 storage to mitigate the current high concentration of CO2 inthe atmosphere. CO2 mineralization in basalts has been proven to be oneof the most secure storage options. For successful implementation and futureimprovements of this technology, the time-dependent deformation behavior ofreservoir rocks in the presence of reactive fluids needs to be studied indetail. We conducted load-stepping creep experiments on basalts from theCarbFix site (Iceland) under several pore fluid conditions (dry,H2O saturated and H2O + CO2 saturated) at temperature,T≈80 ∘C and effective pressure, Peff=50 MPa,during which we collected mechanical, acoustic and pore fluid chemistrydata. We observed transient creep at stresses as low as 11 % of thefailure strength. Acoustic emissions (AEs) correlated strongly with strainaccumulation, indicating that the creep deformation was a brittle process inagreement with microstructural observations. The rate and magnitude of AEswere higher in fluid-saturated experiments than in dry conditions. We inferthat the predominant mechanism governing creep deformation is time- andstress-dependent subcritical dilatant cracking. Our results suggest thatthe presence of aqueous fluids exerts first-order control on creepdeformation of basaltic rocks, while the composition of the fluids playsonly a secondary role under the studied conditions.