An Integrally Embedded Discrete Fracture Model for Flow Simulation in Anisotropic Formations

The embedded discrete fracture model (EDFM), among different flow simulation models, achieves a good balance between efficiency and accuracy. In the EDFM, micro-scale fractures that cannot be characterized individually need to be homogenized into the matrix, which may bring anisotropy into the matrix. However, the simplified matrix–fracture fluid exchange assumption makes it difficult for EDFM to address the anisotropic flow. In this paper, an integrally embedded discrete fracture model (iEDFM) suitable for anisotropic formations is proposed. Structured mesh is employed for the anisotropic matrix, and the fracture element, which consists of a group of connected fractures, is integrally embedded in the matrix grid. An analytic pressure distribution is derived for the point source in anisotropic formation expressed by permeability tensor, and applied to the matrix–fracture transmissibility calculation. Two case studies were conducted and compared with the analytic solution or fine grid result to demonstrate the advantage and applicability of iEDFM to address anisotropic formation. In addition, a two-phase flow example with a reported dataset was studied to analyze the effect of the matrix anisotropy on the simulation result, which also showed the feasibility of iEDFM to address anisotropic formation with complex fracture networks.