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

 
more » « less
Award ID(s):
2045242
NSF-PAR ID:
10395568
Author(s) / Creator(s):
;
Date Published:
Journal Name:
56th U.S. Rock Mechanics/Geomechanics Symposium, Santa Fe, New Mexico, USA, June 2022
Format(s):
Medium: X
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, while the composition of the fluids plays only a secondary role under the studied conditions. 
    more » « 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 both 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. 
    more » « less
  3. Abstract

    To study the mechanical behavior of polymineralic rocks, we performed deformation experiments on two‐phase aggregates of olivine (Ol) + ferropericlase (Per) with periclase fractions (fPer) between 0.1 and 0.8. Each sample was deformed in torsion atT = 1523 K,P = 300 MPa at a constant strain rate to a final shear strain ofγ = 6 to 7. The stress‐strain data and calculated values of the stress exponent,n, indicate that Ol in our samples deformed by dislocation‐accommodated sliding along grain interfaces while Per deformed via dislocation creep. At shear strains ofγ < 1, the strengths of samples withfPer > 0.5 match model predictions for both phases deforming at the same stress, the lower‐strength bound for two‐phase materials, while the strengths of samples withfPer < 0.5 are greater than predicted by models for both phases deforming at the same strain rate, the upper‐strength bound. These observations suggest a transition from a weak‐phase supported to a strong‐phase supported regime with decreasingfPer. Aboveγ = 4, however, the strength of all two‐phase samples is greater than those predicted by either the uniform‐stress or the uniform‐strain rate bound. We hypothesize that the high strengths in the Ol + Per system are due to the presence of phase boundaries in two‐phase samples, for which deformation is rate limited by dislocation motion along interfacial boundaries. This observation contrasts with the mechanical behavior of samples consisting of Ol + pyroxene, which are weaker, possibly due to impurities at phase boundaries.

     
    more » « less
  4. Abstract

    The constitutive behavior of faults is central to many interconnected aspects of earthquake science, from fault dynamics to induced seismicity, to seismic hazards characterization. Yet, a friction law applicable to the range of temperatures found in the brittle crust and upper mantle is still missing. In particular, rocks often exhibit a transition from steady‐state velocity‐strengthening at room temperature to velocity‐weakening in warmer conditions that is poorly understood. Here, we investigate the effect of competing healing mechanisms on the evolution of frictional resistance in a physical model of rate‐, state‐, and temperature‐dependent friction. The yield strength for fault slip depends on the real area of contact, which is modulated by the competition between the growth and erosion of interfacial micro‐asperities. Incorporating multiple healing mechanisms and rock‐forming minerals with different thermodynamic properties allows a transition of the velocity‐ and temperature‐dependence of friction at steady‐state with varying temperatures. We explain the mechanical data for granite, pyroxene, amphibole, shale, and natural fault gouges with activation energies and stress power exponent for weakening of 10–50 kJ/mol and 55–150, respectively, compatible with subcritical crack growth and inter‐granular flow in the active slip zone. Activation energies for the time‐dependent healing process in the range 90–130 kJ/mol in dry conditions and 20–65 kJ/mol in wet conditions indicate the prominence of viscoelastic collapse of microasperities in the absence of water and of pressure‐solution creep, crack healing, and cementation when assisted by pore fluids.

     
    more » « less
  5. Cold sintering densification and coarsening mechanisms are considered from the perspective of the nonequilibrium chemo-mechanical process known in Earth Sciences as pressure solution creep (or dissolutionprecipitation creep). This is an important mechanism of densification and deformation in many geological rock formations in the Earth’s upper crust, and although very slow in nature, it is of direct relevance to the cold sintering process. In cold sintering, we select particulate materials and identify experimental processing parameters to significantly accelerate the kinetics of dissolution-precipitation phenomena, with appropriate consideration of chemistry, applied stress, particle size and temperatures. In the theory of pressure solution, pressure-driven densification is considered to involve the consecutive stages of dissolution at grain contact points, then diffusive transport along the grain boundaries towards open pore surfaces, and then precipitation, all driven by chemical potential gradients. In this study, it is shown that cold sintering of BaTiO3, ZnO and KH2PO4 (KDP) ceramic materials proceeds by the same type of serial process, with the pressure solution creep rate being controlled by the slowest kinetic step. This is demonstrated by the values of activation energy (Ea) for densification, which are in good agreement with the existing literature on dissolution, precipitation, or coarsening. The influence of pressure on the morphology of ZnO grains also supports the pressure solution mechanism. Other characteristics that can be understood qualitatively in terms of pressure solution are observed in the in systems such as BaTiO3 and KDP. We further consider activation energies for grain growth with respect to the precipitation process, as well as evidence for coalescence and Ostwald ripening during cold sintering. For completeness we also consider materials that show significant plastic deformation under compression. Our findings point the way for new advances in densification, microstructural control, and reductions in cold sintering pressure, via the use of appropriate transient solvents in metals and hybrid organic-inorganic systems, such as the Methylammonium lead bromide (MAPBr) perovskite. 
    more » « less