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
This content will become publicly available on August 19, 2024
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):
- 2113474
- NSF-PAR ID:
- 10444085
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical and Analytical Methods in Geomechanics
- Volume:
- 47
- Issue:
- 16
- ISSN:
- 0363-9061
- Format(s):
- Medium: X Size: p. 3018-3043
- Size(s):
- ["p. 3018-3043"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
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.
-
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
-
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
-
more » « less
Abstract The spatial scale and intensity of Earth’s magnetotail current sheet determine the magnetotail configuration, which is critical to one of the most energetically powerful phenomena in the Earth’s magnetosphere, substorms. In the absence of statistical information about plasma currents, theories of the magnetotail current sheets were mostly based on the isotropic stress balance. Such models suggest that thin current sheets cannot be long and should have strong plasma pressure gradients along the magnetotail. Using Magnetospheric Multiscale and THEMIS observations and global simulations, we explore realistic configuration of the magnetotail current sheet. We find that the magnetotail current sheet is thinner than expected from theories that assume isotropic stress balance. Observed plasma pressure gradients in thin current sheets are insufficiently strong (i.e., current sheets are too long) to balance the magnetic field line tension force. Therefore, pressure anisotropy is essential in the configuration of thin current sheets where instability precedes substorm onset.
