A multi-layer reduced model for flow in porous media with a fault and surrounding damage zones

Alessio Fumagalli,Anna Scotti

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

2019-03-04

Abstract:In this work we present a new conceptual model to describe fluid flow in a porous media system in presence of a large fault. Geological faults are often modeled simply as interfaces in the rock matrix, but they are complex structure where the high strain core is surrounded by the so called damage zones, characterized by the presence of smaller fractures which enhance the permeability of the medium. To obtain reliable simulation outcomes these damage zone, as well as the fault, have to be accurately described. The new model proposed in this work considers both these two regions as lower dimensional and embedded in the rock matrix. The model is presented, analyzed, and tested in several configurations to prove its robustness and ability to capture many important features, such as hight contrast and heterogeneity of permeability.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the fluid flow problem in the presence of large faults and their surrounding damaged zones in porous media. Specifically, geological faults are usually simply modeled as interfaces in the rock matrix, but they are actually complex structures, with a high - strain core (i.e., the fault core) surrounded by so - called damaged zones, and these zones enhance the permeability of the medium due to the presence of smaller fractures. In order to obtain reliable simulation results, these damaged zones as well as the faults must be accurately described.

### Main problems and challenges

1. **Description of complex structures**: A fault is not just a simple interface; it includes a high - strain core area (fault core), and this core is surrounded by highly fragmented areas (damaged zones). These areas contain many fractures that are smaller in scale than the fault itself, and these fractures will change the nature of the local flow paths.

2. **Differences in thickness**: The thickness of a fault is usually several orders of magnitude smaller than other characteristic sizes in the porous medium, but its physical properties may be very different from those of the porous medium. Due to filling processes and chemical reactions, these objects may be partially or completely blocked, so they can serve as preferential channels for fluid flow or geological barriers.

3. **Limitations of existing models**: Traditional methods characterize this area by assigning different permeabilities to the faults, but because of its small thickness, it is impractical to incorporate it into the grid representing the rock matrix. A common method is to consider the fault as a low - dimensional object and derive new conceptual models to describe the flow and pressure behavior inside and across these objects.

### Goals of the paper

This paper proposes a new multi - layer simplified model, aiming to:
- Consider both the fault and the damaged zone as low - dimensional objects embedded in the rock matrix simultaneously.
- Verify the robustness of this model and its ability to capture important features, such as high - contrast and heterogeneous permeability, through numerical analysis and testing.
- Provide a method that can simulate different scenarios without extreme grid refinement, where the rock matrix, the damaged zone, and the fault may have different permeability values.

By introducing this multi - layer interface law, the model can handle different openings and asymmetries on both sides of the fault, thereby better simulating the fluid flow in actual geological situations.