skip to main content


Title: Micromechanical modeling for rate‐dependent behavior of salt rock under cyclic loading
Summary

The dependence of rock behavior on the deformation rate is still not well understood. In salt rock, the fundamental mechanisms that drive the accumulation of irreversible deformation, the reduction of stiffness, and the development of hysteresis during cyclic loading are usually attributed to intracrystalline plasticity and diffusion. We hypothesize that at low pressure and low temperature, the rate‐dependent behavior of salt rock is governed by water‐assisted diffusion along grain boundaries. Accordingly, a chemo‐mechanical homogenization framework is proposed in which the representative elementary volume (REV) is viewed as a homogeneous polycrystalline matrix that contains sliding grain‐boundary cracks. The slip is related to the mass of salt ions that diffuse along the crack surface. The relationship between fluid inclusion‐scale and REV‐scale stresses and strains is established by using the Mori–Tanaka homogenization scheme. It is noted from the model that a lower strain rate and a larger number of sliding cracks enhance stiffness reduction and hysteresis. Thinner sliding cracks (i.e., thinner brine films) promote stiffness reduction and accelerate stress redistributions. The larger the volume fraction of the crack inclusions, the larger the REV deformation and the larger the hysteresis. Results presented in this study shed light on the mechanical behavior of salt rock that is pertinent to the design of geological storage facilities that undergo cyclic unloading, which could help optimize the energy production cycle with low carbon emissions.

 
more » « less
Award ID(s):
1552368 1361996
NSF-PAR ID:
10455427
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical and Analytical Methods in Geomechanics
Volume:
45
Issue:
1
ISSN:
0363-9061
Page Range / eLocation ID:
p. 28-44
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The operation of fracture, diffusion, and intracrystalline‐plastic micromechanisms during semibrittle deformation of rock is directly relevant to understanding mechanical behavior across the brittle‐plastic transition in the crust. An outstanding question is whether (1) the micromechanisms of semibrittle flow can be considered to operate independently, as represented in typical crustal strength profiles across the brittle to plastic transition, or (2) the micromechanisms are coupled such that the transition is represented by a distinct rheology with dependency on effective pressure, temperature, and strain rate. We employ triaxial stress‐cycling experiments to investigate elastic‐plastic and viscoelastic behaviors during semibrittle flow in two distinctly different monomineralic, polycrystalline, synthetic salt‐rocks. During semibrittle flow at high differential stress, granular, low‐porosity, work‐hardened salt‐rocks deform predominantly by grain‐boundary sliding and wing‐crack opening accompanied by minor intragranular dislocation glide. In contrast, fully annealed, near‐zero porosity salt‐rocks flow at lower differential stress by intragranular dislocation glide accompanied by grain‐boundary sliding and opening. Grain‐boundary sliding is frictional during semibrittle flow at higher strain rates, but the associated dispersal of water from fluid inclusions along boundaries can activate fluid‐assisted diffusional sliding at lower strain rates. Changes in elastic properties with semibrittle flow largely reflect activation of sliding along closed grain boundaries. Observed microstructures, pronounced hysteresis and anelasticity during cyclic stressing after semibrittle flow, and stress relaxation behaviors indicate coupled operation of micromechanisms leading to a distinct rheology (hypothesis 2 above).

     
    more » « less
  2. Abstract

    Crack interactions leading to shear localization were quantified using microstructural analysis for brittle faults and high‐temperature ductile faults formed during experiments on quartz sandstone. In both faulting regimes, the nucleation of macroscopic faults results from the interactions of microfractures at two length scalesin ensemble. Brittle faults nucleate when the longest mesoscale shear fractures and transgranular tensile cracks critically interact. In contrast, ductile faults nucleate when the longest mesoscale shear fractures and multi‐grain scale intergranular shear cracks critically interact. For both faulting regimes, we conclude the interaction and coalescence of the longest mesoscale shear fractures is the fundamental process responsible for fault nucleation. Hence, mesoscale shear fractures, which accommodate the majority of axial strain prior to shear localization in both faulting regimes, also serve as the nucleus of macroscopic faults. Locally, the growth of the mesoscale shear fractures is promoted by the interaction and coalescence of the multi‐grain scale cracks in both faulting regimes. We hypothesize that attainment of a critical microstructure for shear localization (i.e., local clustering of the longest microfractures) requires a characteristic amount of plastic axial strain, which depends on deformation conditions. In brittle faulting, distributed microfracturing is confined within limited regions of the rock volume, which expedites crack clustering and fault nucleation at low characteristic strains. In ductile faulting, distributed microfracturing occurs more uniformly throughout the rock volume, delaying shear localization to high characteristic strains. Accurate prediction of shear localization requires models that describe crack interactions of thelargestflaws that account for crackclustering.

     
    more » « less
  3. 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
  4. Abstract

    Subduction interfaces exhibit a variety of slip behaviors, including megathrust and slow earthquakes. Field observations are consistent with crack‐seal deformation, in which tensile cracks are sealed by fluid‐transported solute. However, there are few constraints on the mass fluxes and length and time scales of such deformation and attendant increases in cohesion within the seismogenic zone. Here, we present a systematic geochemical investigation of mass transport associated with development of crack‐seal veins in the Shimanto Belt—an accretionary complex that preserves a record of plate boundary slip behavior at temperatures relevant to the seismogenic zone: 150–350°C. These mélanges show evidence for shear across decameter‐scale zones of deformation dominated by anastomosing scaly fabrics and pervasive veins. We use meso‐ and microstructural observations with geochemical observations of scaly fabrics and crack‐seal veins to evaluate the role of silica redistribution in healing fracture porosity along the plate interface and modulating slip behavior. Crack‐seal veins contain primarily quartz with albite and calcite, and vein textures provide evidence for partial sealing. Bulk rock analyses determined that the amount of phyllosilicates, specifically illite and chlorite, increases with temperature. A simple mass‐balance model based on the immobile chemical component TiO2shows a systematic trend in mobility for all three mélanges and increased element mobility as a function of temperature. Scaly fabrics and veins show compositional evidence for locally sourced mineral redistribution. This study supports a model where development of scaly slip surfaces and fracture healing through temperature‐dependent mineral redistribution can impact slip behavior and fluid flow along the subduction plate interface.

     
    more » « less
  5. Abstract

    For brittle friction and rock deformation, the coefficientαin the general effective stress relationσe = σ − αPpcan be approximated as unity with sufficient accuracy. However, it is uncertain ifαdeviates from unity for semibrittle flow when both brittle and intracrystalline‐plastic deformation is involved. We conducted triaxial and isostatic compression experiments on synthetic salt‐rocks (∼300 ppm water) at room temperature to test the effective stress relation in the semibrittle regime using silicone oil and argon gas as pore fluids. Confining and pore pressures were cycled while their difference (differential pressure) was kept constant, such that changes in the mechanical behavior would indicate deviation ofαfrom unity. Microstructural observations were used to determine the dependence ofαon true area of grain contact from asperity yielding. In triaxial compression experiments, semibrittle flow involves grain boundary cracking and sliding, and intragranular dislocation glide and cracking. Flow strength remains constant for changes in pore fluid pressure of more than two orders of magnitude. In isostatic compression experiments, samples show combined processes of microcracking, grain boundary sliding, dislocation glide, and fluid‐assisted grain boundary migration recrystallization. Volumetric strain depends directly on the differential pressures (i.e.,αequals one). Analysis of grain‐contact area in both experiments indicates thatαis independent of the true area of contact defined by plastic yielding at grain boundaries. The observation ofαeffectively equals one may be explained by operation of pressure‐independent intracrystalline‐plastic mechanisms and transmission of pore pressure at grain boundaries through thin fluid films.

     
    more » « less