skip to main content

Explore Research Products in the NSF-PAR

Title: Viscoplastic modeling of elastic and creep deformation of fractured Berea sandstone.
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.  more » « less
Award ID(s):
1829597 1727661
NSF-PAR ID:
10252459
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
55th US Rock Mechanics/Geomechanics Symposium
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , 56th U.S. Rock Mechanics/Geomechanics Symposium, Santa Fe, New Mexico, USA, June 2022)
    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
  2. ; ; ( , 56th U.S. Rock Mechanics/Geomechanics Symposium)
    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 of 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).

     
    more » « less
  3. ; ; ( , 55th US Rock Mechanics/Geomechanics Symposium)
    null (Ed.)
    Fracturing in brittle rocks exhibits a significant nonlinear region surrounding the crack tip called the fracture process zone (FPZ). In this study, the evolution of the FPZ under pure mode II loading using notched deep beam under three-point loading was investigated. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique to characterize various crack characteristics such as its type and FPZ evolution in Barre granite specimens. Both displacement and strain approaches of the two-dimensional digital image correlation were used to identify the mode of fracture under pure mode II loading. Both approaches showed that the crack initiation occur under mode I despite the pure mode II loading at the notch tip. The displacement approach was used for characterizing the evolution of the FPZ which analyzed the crack tip opening displacement and crack tip sliding displacement to identify the transition between the three stages of FPZ evolution, namely, (a) elastic stage, (b) formation of the FPZ, and (c) the macro-crack initiation. The results showed that the evolution of the FPZ of mode I fracture under pure mode II loading is similar to cases of pure mode I loading of the same rock. 
    more » « less
  4. ; ( , Journal of Geophysical Research: Solid Earth)
    Abstract

    Brittle fracture propagation in rocks is a complex process due to significant grain‐scale heterogeneity and evolving stress states under dynamic loading conditions. In this work, we use digital image correlation and linear elastic fracture mechanics to make instantaneous measurements of the opening (mode I) and in plane shear (mode II) components of the stress intensity field during dynamic mixed mode crack initiation and propagation in crystalline and granular rocks. Both rock types display some similar fracture behaviors as observed in engineered materials, including rate dependent fracture initiation toughness and a direct relationship between propagation toughness and crack velocity; however, measured propagation toughness is higher than quasi‐static values at crack velocities well below the branching velocity in both rocks. Additionally, due to grain scale controls on the fracture process, mixed mode crack propagation is fundamentally different between these two rock types. Mixed mode propagation is energetically more favorable than pure opening mode propagation in sandstone, while the opposite is true in granite. Furthermore, following initiation, propagation in granite occurs so as to minimize the mode II contribution, irrespective of the initiation conditions, while fractures in sandstone maintain a non‐negligible mode II contribution during propagation across the sample.

     
    more » « less
  5. ; ( , Journal of Geophysical Research: Solid Earth)
    Abstract

    Critical processes including seismic faulting, reservoir compartmentalization, and borehole failure involve high‐pressure mechanical behavior and strain localization of sedimentary rocks such as sandstone. Sand is often used as a model material to study the mechanical behavior of poorly lithified sandstone. Recent studies exploring the multi‐scale mechanics of sand have characterized the brittle, low‐pressure regime of behavior; however, limited work has provided insights into the ductile, high‐pressure regime of behavior viain‐situmeasurements. Critical features of the ductile regime, including grain breakage, grain micromechanics, and volumetric strain behavior therefore remain under‐explored. Here, we use a new high‐pressure triaxial apparatus within‐situx‐ray tomography to provide new insights into deformation banding, grain breakage, and grain micromechanics in Ottawa sand subjected to triaxial compression under confining pressures between 10 and 45 MPa. We observed strain‐hardening at pressures above 15 MPa and strain‐neutral responses at pressures below 15 MPa. Compacting shear bands and grain breakage were observed at all pressures with no significant variation due to grain size, except for minor increases in breakage in less‐rounded sands. Grain breakage emerged at stress levels lower than the assumed yield threshold and more intense breakage was associated with thinner deformation bands. Contact sliding at inter‐grain contacts demonstrated a bifurcation into a bimodal distribution, with intense sliding within deformation bands and reduced but non‐negligible sliding outside of deformation bands, suggesting that off‐band zones remain mechanically active during strain hardening.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email