Modeling fluid flow in three-dimensional (3D) Discrete Fracture Networks (DFNs) is of relevance in many engineering applications, such as hydraulic fracturing, oil/gas production, geothermal energy extraction, nuclear waste disposal and CO2 sequestration. A new Boundary Element Method (BEM) technique with discontinuous quadratic elements and a parallel Domain Decomposition Method (DDM) is presented herein for the simulation of the steady-state fluid flow in 3D DFN systems with wellbores, consisting of planar fractures having arbitrary properties and wellbore trajectories. Numerical examples characterized by DFNs of increasing complexity are investigated to evaluate the accuracy and efficiency of the presented technique. The results show that accurate solutions can be obtained with less nodes than with mesh-based methods (e.g. Finite Element Method). In addition, the DDM algorithm used provides a quite fast convergence. The simulation results of the fluid flow around intersections among traces (linear intersections between fractures), intersections between traces and a fracture boundaries, and wellbore intersections is accurate. Source code is available at : https://github.com/BinWang0213/PyDFN3D.
Modeling heat and fluid flow in discrete fracture networks : application to fluid-driven seismicity and engineered geothermal systems