skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Osmotic Ion Concentration Control of Steady-State Subcritical Fracture Growth in Shale
ABSTRACT:The mechanism of formation of natural cracks in sedimentary rocks in the geologic past is an important problem in hydraulic fracturing. Why are the natural cracks roughly parallel and equidistant, and why is the spacing in the order of 10 cm rather than 1 cm or 100 cm? Fracture mechanics alone cannot answer these questions. Here it is proposed that fracture mechanics must be coupled with the diffusion of solute ions (Na+ and Cl− are considered here), driven by an osmotic pressure gradient. Parallel equidistant cracks are considered to be subcritical and governed by the Charles-Evans law. The evolution in solute concentration also affects the solvent pressure in the pores and cracks, altering the resistance to frictional sliding. Only steady-state propagation and periodic cracks are studied. An analytical solution of the crack spacing as a function of the properties of the rock as well as the solvent and solute, and the imposed far-field deformation is obtained. Finally, the stability of the growth of parallel cracks is proven by examining the second variation of free energy. Stability of the periodic growth state is also considered. 1 INTRODUCTIONThe deep layers of sedimentary rocks such as shale and sandstone are usually intersected by systems of nearly parallel natural cracks either filled by mineral deposits or closed by creep over a million year life span. Their spacing is roughly uniform and is on the order of 0.1 m (rather than 1 m or 0.01m). These cracks likely play an important role in hydraulic fracturing for gas or oil recovery (aka fracking, fraccing or frac) (Rahimi-Aghdam et al., 2019, e.g.). Therefore, understanding the mechanism of their formation in the distant geologic past is of interest.What controls the spacing of the nearly parallel cracks in shale? According to the fracture mechanics alone, the crack spacing is arbitrary. If propagating parallel equidistant cracks are in a critical state, stability analysis shows that many cracks would have to stop growing, causing a great increase of their average spacing, which was obviously not the case (Bažant et al., 2014).  more » « less
Award ID(s):
2029641
PAR ID:
10626899
Author(s) / Creator(s):
;
Publisher / Repository:
ARMA
Date Published:
Format(s):
Medium: X
Location:
Santa Fe, New Mexico, USA
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, ARMA)
    ABSTRACT:Creation of a fracture network in a hydraulic fracturing process is essential for subsurface energy extraction and CO2 sequestration. It is facilitated by reactivation of pre-existing intersecting weak layers and cemented cracks in the rock. In this study, a poromechanical model is developed for the hydraulic fracturing process in rocks containing such pre-existing weak layers. Based on the mixture theory, the crack band model is used to simulate the growth of a crack system. The governing equations with the parameters for hydromechanical coupling are derived, to describe the evolution of the opening and branching of cracks caused by water injection. Microplane model M7 is adopted to characterize the deformation and fracturing of the solid skeleton of the rock, and the Poiseuille law is used to characterize fluid flow through the hydraulic fractures. Numerical simulations are performed to reproduce and interpret recently published laboratory-scale hydraulic fracturing experiments conducted at Los Alamos National Laboratory (LANL). In these experiments, the rock was represented by confined plaster slabs containing orthogonal intersecting weak layers of higher porosity. Numerical simulations reveal how poromechanical characteristics such as the Biot coefficient and the fluid injection rate lead to various typical fracture modes observed in the experiments. These modes include formation of one dominant planar crack or various orthogonal fracture networks. 
    more » « less
  2. ; (, 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. INTRODUCTIONDetermination 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
  3. ; ; (, 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. INTRODUCTIONThe 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
  4. ; ; ; ; (, International Journal for Numerical and Analytical Methods in Geomechanics)
    Abstract Many geo‐engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase‐field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase‐field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy‐based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase‐field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase‐field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. The comparison of numerical results for two‐dimensional test cases with existing analytical solutions demonstrates that the proposed phase‐field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three‐dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios. 
    more » « less
  5. (, International journal of fracture)
    The interaction between cracks and inclusions plays an important role in the fracture behavior of particulate composites. It is commonly recognized that an inclusion stiffer than the matrix tends to deflect an approaching crack away while a softer inclusion attracts the crack. Here, we demonstrate by analytical modeling and numerical simulations that the crack-inclusion interaction can be tuned by an applied T-stress. Under a sufficiently large compressive applied T-stress, cracks can be attracted to stiffer inclusions while repelled by softer ones, thus reversing the conventional trend. Potential applications of this work include composite electrodes in lithium-ion batteries and hydraulic fracturing. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email