Search

Understanding deformations and fluid flow in fractured rocks is of central importance for many subsurface flow applications. Thus, numerical frameworks are needed that capture the coupled mechanical and hydraulic behaviour, including scenarios with complex fracture networks. This paper employs the extended finite volume method to represent fracture manifolds in a poroelastic matrix domain and to compute shear and tensile displacements of fracture segments. However, using embedded fractures with non‐conforming grids can lead to severe convergence issues while computing shear slip and tensile opening of intersecting and parallel fractures, particularly if they have similar slopes. Our proposed solution of this problem is to slightly deform the fracture geometry by merging critical segments. We discuss and show examples of which attributes the merged segment has to adopt to ensure the correct slip behaviour at an intersection. Thereby it is crucial that the flow topology remains unaltered. Analysing failures of kinked fractures and fracture intersections show the flexibility of the devised algorithm. Moreover, it is demonstrated that deformation and opening patterns of a natural network with more than 200 fractures are qualitatively predicted very well. This method is a simple solution to treat multiple fracture segments in a grid cell, hence it can be used for any model embedding fractures in non‐conforming grids. [ABSTRACT FROM AUTHOR]

Fractured-reservoir relative permeability, water breakthrough, and recovery cannot be extrapolated from core samples, but computer simulations allow their quantification through the use of discrete fracture models at an intermediate scale. For this purpose, we represent intersecting naturally and stochastically generated fractures in massive or layered porous rock with an unstructured hybrid finite-element (FE) grid. We compute two-phase flow with an implicit FE/finite volume (FV) method (FE/FVM) to identify the emergent properties of this complex system. The results offer many important insights: Flow velocity varies by three to seven orders of magnitude and velocity spectra are multimodal, with significant overlaps between fracture- and matrix- flow domains. Residual saturations greatly exceed those that were initially assigned to the rock matrix. Total mobility is low over a wide saturation range and is very sensitive to small saturation changes. When fractures dominate the flow, but fracture porosity is low (10-3 to 1%), gridblock average relative permeabilities. kr,avg>, cross over during saturation changes of less than 1%. Such upscaled kr,avg> yield a convex, highly dispersive fractional-flow function without a shock. Its shape cannot be matched with any conventional model, and a new formalism based on the fracture/ matrix flux ratio is proposed. Spontaneous imbibition during waterflooding occurs only over a small fraction of the total fracture/matrix-interface area because water imbibes only a limited number of fractures. Yet in some of these, flow will be sufficiently fast for this process to enhance recovery significantly. We also observe that a rate dependence of recovery and water breakthrough occurs earlier in transient-state flow than in steady-state flow. [ABSTRACT FROM AUTHOR]