Continuum modeling of the effect of surface area growth due to crushing and damage on the permeability of granular rocks

Shiva Esna Ashari,Arghya Das,Giuseppe Buscarnera

DOI: https://doi.org/10.48550/arXiv.1612.03221

2016-12-10

Abstract:This paper discusses a continuum approach to track the evolution of permeability in granular rocks by accounting for the combined effect of porosity changes, grain breakage and cement bond damage. To account for such a broad range of microscopic processes under general loading paths, the Breakage Mechanics theory is used and the computed mechanical response is linked with the Kozeny equation, i.e. a permeability model able to evaluate the reduction of the hydraulic conductivity resulting from the simultaneous loss of porosity and growth of surface area. In particular, the evolution of the internal variables of the model has been linked to idealized geometric schemes at particle scale, with the goal to distinguish the contribution of the fines generated by the disaggregation of the cement matrix from that of the broken fragments resulting from the crushing of the skeleton. Compression/flow experiments available in the literature for different granular rocks are used to validate the proposed methodology. The analyses illustrate that the drop of the permeability of damaged rocks would be severely underestimated without an accurate computation of the growth of surface area, as well as that the distributed fragmentation of skeleton particles tends to have stronger implications than the generation of cement fines. These findings, along with the satisfactory agreement between model predictions and experiments, stress the benefits of adopting microstructure-based constitutive laws for the analysis of coupled hydro-mechanical problems.

Geophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to track the influence of the surface area growth on the permeability of granular rocks during the crushing and damage processes through the continuum model. Specifically, the paper focuses on how to accurately evaluate the influence of these processes on the rock permeability while considering multiple microscopic processes such as porosity change, particle crushing, and cement bond damage.

### Analysis of the Core Problems in the Paper

1. **Problem Background**:
- Fluid flow in granular rocks is usually coupled with multiple deformation mechanisms such as particle crushing, compaction, strain localization, and crack propagation.
- These changes in microstructure directly affect the overall permeability of the rock. Therefore, understanding these relationships is crucial for activities involving fluid injection/extraction, such as oil and gas extraction, underground storage, etc.

2. **Deficiencies in Existing Research**:
- Although there have been a large number of experimental studies and empirical formulas for evaluating the change in rock permeability, most existing continuum models lack coordinated mechanical and hydraulic constitutive laws based on the same microstructure descriptors.
- Many models only rely on the empirical permeability function of porosity change and ignore the important influence of surface area growth.

3. **Objectives of the Paper**:
- Propose a continuum model to quantify the change in the permeability of granular rocks by tracking the changes in microstructure descriptors that occur simultaneously with porosity loss and surface area growth.
- Use the Breakage Mechanics theory to simulate the processes of particle crushing and cement depolymerization, and describe these processes through two internal variables (the breakage index \( B \) and the damage index \( D \)).
- Relate these micromechanical characteristics to the changes in porosity and surface area for the constitutive analysis applicable to hydro - mechanical coupling problems.

### Key Methods and Technologies

1. **Fracture Mechanics Model**:
- Use the fracture mechanics theory to describe the evolution of the Grain Size Distribution (GSD).
- Introduce two internal variables: the breakage index \( B \) and the damage index \( D \), which are associated with the particle size distribution and the cement size distribution respectively.

2. **Permeability Model**:
- Use the Kozeny equation to describe the change in permeability, which can evaluate the influence of porosity loss and surface area growth on hydraulic conductivity.
- Verify the validity of the model through experimental data, especially for different types of granular rocks (such as glass beads, Ottawa sand, and Ottershaw sandstone).

3. **Numerical Simulation**:
- Conduct numerical simulations for different loading paths, including uniaxial compression and radial loading paths.
- Evaluate the accuracy of the model in predicting the change in permeability by comparing the model predictions and experimental results.

### Main Findings

- Model predictions show that ignoring surface area growth will lead to a serious underestimation of the decrease in the permeability of damaged rocks.
- The distributed crushing of skeletal particles has a more significant impact on permeability than the generation of fine cement particles.
- The good agreement between the model predictions and the experimental data demonstrates the importance of using microstructure - based constitutive laws in analyzing coupled hydro - mechanical problems.

### Formula Presentation

1. **Darcy's Law**:
\[
Q = -\frac{\kappa A}{\mu} \frac{\partial P}{\partial x}
\]
where \( Q \) is the flow rate, \( \kappa \) is the intrinsic permeability of the medium, \( A \) is the flow cross - sectional area, \( \mu \) is the viscosity of the fluid, and \( \frac{\partial P}{\partial x} \) is the pressure gradient.

2. **Kozeny Equation**:
\[
\kappa = \frac{c_0 n^3}{S^2}
\]
where \( c_0 \) is the Kozeny constant, \( n \) is the porosity, and \( S \) is the specific surface area (i.e., the total surface area of solid materials per unit volume).

3. **Fracture Energy and Damage Energy**:
\[
E_B = -\fr