We present a method for porous media flow in the presence of complex fracture networks. The approach uses the Mimetic Finite Difference method (MFD) and takes advantage of MFD's ability to solve over a general set of polyhedral cells. This flexibility is used to mesh fracture intersections in two and three-dimensional settings without creating small cells at the intersection point. We also demonstrate how to use general polyhedra for embedding fracture boundaries in the reservoir domain. The target application is representing fracture networks inferred from microseismic analysis.
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The code used to generate the results was written in Python and used the SciPy and NumPy libraries[Jones et al. (2001)]. The Python package Shapely was used for two-dimensional polygon calculations.Visualization was done using Paraview [Henderson (2007)]. For cell volume and centroid computations,the code uses the algorithm defined in Mirtich (1996). The Authors would like to acknowledge help fromDr. Ivan Yotov, Dr. Mark Mear, and Dr. Xin Yang. A portion of this research was supported by the KingAbdullah University of Science and Technology Academic Excellence Alliance and Saudi Aramco.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.