Source code and examples for the paper
"A multi-layer reduced model for flow in porous media with a fault and surrounding damage zones" by Alessio Fumagalli, Anna Scotti.
See arXiv pre-print.
Runscripts for all test cases of the work available example1 and example2.
Note that you may have to revert to an older version of PorePy to run the examples.
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.
If you use this work in your research, we ask you to cite the following publication arXiv pre-print.
If you want to run the code you need to install PorePy and you might revert it.
Newer versions of PorePy may not be compatible with this repository.
See license.