skip to main content

This content will become publicly available on August 19, 2024

Title: A resolved CFD–DEM investigation into the onset of suffusion: effect of confining pressure and stress anisotropy

The susceptibility of a granular soil to suffusion is strongly dependent on its grain size distribution (GSD) and the mechanical and hydraulic conditions it is subjected to. This study investigates the onset of suffusion considering the effect of confining pressure and stress anisotropy using a fully resolved computational fluid dynamics and discrete element method (CFD–DEM). Three benchmarks, including the sedimentations of single and two adjacent spheres and the classic one‐dimensional (1D) consolidation are performed to demonstrate the capability of this method for high‐fidelity particle‐fluid simulations. A modified hydraulic criterion for the onset of suffusion considering stress anisotropy is presented. The microstructural changes of soil specimens before and during global suffusion are inspected, with emphasis on the evolutions of particle kinetic energy and displacements, force chain networks, and stress anisotropy. We found that the critical hydraulic gradient is negatively correlated with the confining pressure and the degree of stress anisotropy. Fine particles in the soil matrix are locally detached at small hydraulic gradients before the apparent global suffusion, as manifested by the variation of particle kinetic energy and coordination numbers. The roles of different contact types on force transmission and stress anisotropy in eroded specimens are also examined.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical and Analytical Methods in Geomechanics
Medium: X Size: p. 3018-3043
["p. 3018-3043"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Carbonate sediments play a prominent role on the global geological stage as they store more than $$60\%$$ 60 % of world’s oil and $$40\%$$ 40 % of world’s gas reserves. Prediction of the deformation and failure of porous carbonates is, therefore, essential to minimise reservoir compaction, fault reactivation, or wellbore instability. This relies on our understanding of the mechanisms underlying the observed inelastic response to fluid injection or deviatoric stress perturbations. Understanding the impact of deformation/failure on the hydraulic properties of the rock is also essential as injection/production rates will be affected. In this work, we present new experimental results from triaxial deformation experiments carried out to elucidate the behaviour of a porous limestone reservoir analogue (Savonnières limestone). Drained triaxial and isotropic compression tests were conducted at five different confining pressures in dry and water-saturated conditions. Stress–strain data and X-ray tomography images of the rock indicate two distinct types of deformation and failure regimes: at low confinement (10 MPa) brittle failure in the form of dilatant shear banding was dominant; whereas at higher confinement compaction bands orthogonal to the maximum principal stress formed. In addition to the pore pressure effect, the presence of water in the pore space significantly weakened the rock, thereby shrinking the yield envelope compared to the dry conditions, and shifted the brittle–ductile transition to lower effective confining pressures (from 35 MPa to 29 MPa). Finally, permeability measurements during deformation show a reduction of an order of magnitude in the ductile regime due to the formation of the compaction bands. These results highlight the importance of considering the role of the saturating fluid in the brittle–ductile response of porous rocks and elucidate some of the microstructural processes taking place during this transition. 
    more » « less
  2. 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
  3. Hydraulic fracturing arises as a method to enhance oil and gas production, and also as a way to recover geothermal energy. It is, therefore, essential to understand how injecting a fluid inside a rock reservoir will affect its surroundings. Hydraulic fracturing processes can be strongly affected by the interaction between two mechanisms: the elastic effects caused by the hydraulic pressure applied inside fractures and the poro-mechanical effects caused by the fluid infiltration inside the porous media (i.e. fluid diffusivity); this, in turn, is affected by the injection rate used. The interaction between poro-elastic mechanisms, particularly the effect of the fluid diffusivity, in the hydraulic fracturing processes is not well-understood and is investigated in this paper. This study aims to experimentally and theoretically comprehend the effects of the injection rate on crack propagation and on pore pressures, when flaws pre-fabricated in prismatic gypsum specimens are hydraulically pressurized. In order to accomplish this, laboratory experiments were performed using two injection rates (2 and 20 ml/min), applied by an apparatus consisting of a pressure enclosure with an impermeable membrane in both faces of the specimen, which allowed one to observe the growth of a fluid front from the pre-fabricated flaws to the unsaturated porous media (i.e. rock), before fracturing took place. It was observed that the fracturing pressures and patterns are injection-rate-dependent. This was interpreted to be caused by the different pore pressures that developed in the rock matrix, which resulted from the significantly distinct fluid fronts observed for the two injection rates tested. 
    more » « less
  4. Abstract

    Elevated pore fluid pressure is proposed to contribute to slow earthquakes along shallow subduction plate boundaries. However, the processes that create high fluid pressure, disequilibrium compaction and dehydration reactions, lead to different effective stress paths in fault rocks. These paths are predicted by granular mechanics frameworks to lead to different strengths and deformation modes, yet granular mechanics do not predict their effects on fault stability. To evaluate the role of fluid overpressure on shallow megathrust strength and slip behavior, we conducted triaxial shear experiments on chlorite and celadonite rich gouge layers. Experiments were conducted at constant temperature (130 and 100°C), confining pressure (130 and 140 MPa), and pore fluid pressures (between 10 and 120 MPa). Fluid overpressure due to disequilibrium compaction was simulated by increasing confining and pore fluid pressure synchronously without exceeding the target effective pressure, whereas overpressure due to dehydration reactions was simulated by first loading the sample to a target isotropic effective pressure and then increasing pore fluid pressure to reduce the effective pressure. We find that the effects of fluid pressure and stress path on the mechanical behavior of the chlorite and celadonite gouges can generally be described using the critical state soil mechanics (CSSM) framework. However, path effects are more pronounced and persist to greater displacements in chlorite because its microstructure is more influenced by stress path. Due to its effects on microstructure, the stress path also imparts greater control on the rate‐dependence of chlorite strength, which is not predicted by CSSM.

    more » « less
  5. Natural soil deposits can consist of particles with a wide range of sizes. In current practice, the assessment of shear strength and stress-dilatancy behavior of coarse-grained soils is based on methods developed for poorly graded sands, without explicit consideration for differences in gradation. This paper investigates the influence of the range of particle sizes on the monotonic shear strength and the stress-dilatancy response of poorly- to well-graded soils. Using the 3D discrete element method (DEM), the applicability of commonly used sand-based stress-dilatancy frameworks is assessed for a range of gradations. This DEM investigation employs clumps of spheres to accurately simulate the particle shapes on specimens with coefficients of uniformity (CU) varying between 1.9 and 6.9. These specimens were subjected to isotropically consolidated drained triaxial compression at various relative densities and confining stresses with the objective of isolating the effects of particle size distribution from those of particle shape. The peak and critical state shear strengths and the dilatancy responses of the specimens with different gradations are evaluated. For the same state parameter, the results indicate an increase in the shear strength and rate of dilation as the range of particle sizes increases. However, the critical state line shifts downward, and its slope decreases as CU is increased. The DEM results are compared to Bolton’s stress-dilatancy relationship to highlight the inadequacies of using clean sand-based frameworks in capturing the behavior of well-graded soils. 
    more » « less