Control‐volume method for numerical simulation of two‐phase immiscible flow in two‐ and three‐dimensional discrete‐fractured media

Abstract

[1] We provide a numerical procedure for the simulation of two-phase immiscible and incompressible flow in two-and three-dimensional discrete-fractured media. The concept of cross-flow equilibrium is used to reduce the fracture dimension from n to (n-1) in the calculation of flow in the fractures. This concept, which is often referred to as the discrete-fracture model, has a significant effect on the reduction of computational time. The spatial discretization is performed with the control-volume method. This method is locally conservative and allows the use of unstructured grids to represent complex geometries, such as discrete-fracture configurations. The relative permeability is upwinded with a criterion based on the evaluation of the flux direction at the boundaries of the control volumes, which is consistent with the physics of fluid flow. The system of partial differential equations is decoupled and solved using the implicit-pressure, explicitsaturation (IMPES) approach. The algorithm has been successfully tested in two-and three-dimensional numerical simulations of wetting phase fluid injection (such as water) in discrete-fractured media saturated by a nonwetting phase (such as nonaqueous phase liquid or oil) with mild to high nonlinearity in relative permeability and capillary pressure. To the best of our knowledge, results for simulations of two-phase immiscible and incompressible flow in three-dimensional discrete-fractured media, including capillary and gravity effects, are the first to appear in the literature.

Generation of coarse‐scale continuum flow models from detailed fracture characterizations

“…[25] For simulations of the fully resolved (fine scale) model and the local upscaling calculations (i.e., solutions of equations (7) and (8)), any discrete fracture simulation procedure [e.g., Bogdanov et al, 2003aBogdanov et al, , 2003bMonteagudo and Firoozabadi, 2004;Matthäi et al, 2005;Hoteit and Firoozabadi, 2005] could be applied. In this work we use a recently developed finite volume based discrete fracture model, presented by Karimi-Fard et al [2004] and also described (for a different application involving flow in systems characterized by thin but extensive lowpermeability compaction bands, which act as ''antifractures'') by Sternlof et al [2006].…”

Development of an Efficient Embedded Discrete Fracture Model for 3D Compositional Reservoir Simulation in Fractured Reservoirs

“…We used the above matrix capillary pressure relation since this was the relation used by Monteagudo and Firoozabadi (2004). The capillary pressure equation presented by Brooks and Corey (1964) would have also been an alternative equation to be used in this example.…”

Discrete fracture network modelling: A computational algorithms

“…Для построения конечномерной задачи в уравнениях (30), (31) производится замена функций и на их непрерывные, кусочно-линейные (10) и кусочно-постоянные ана-логи (12) [10]- [12]:…”