Fluid flow modeling of naturally fractured reservoirs remains a challenge because of the complex nature of fracture systems controlled by various chemical and physical phenomena. A discrete fracture network (DFN) model represents an approach to capturing the relationship of fractures in a fracture system. Topology represents the connectivity aspect of the fracture planes, which have a fundamental role in flow simulation in geomaterials involving fractures and the rock matrix. Therefore, one of the most-used methods to treat fractured reservoirs is the double porosity-double permeability model. This approach requires the shape factor calculation, a key parameter used to determine the effects of coupled fracture-matrix fluid flow on the mass transfer between different domains. This paper presents a numerical investigation that aimed to evaluate the impact of fracture topology on the shape factor and equivalent permeability through hydraulic connectivity (f). This study was based on numerical simulations of flow performed in discrete fracture network (DFN) models embedded in finite element meshes (FEM). Modeled cases represent four hypothetical examples of fractured media and three real scenarios extracted from a Brazilian pre-salt carbonate reservoir model. We have compared the results of the numerical simulations with data obtained using Oda’s analytical model and Oda’s correction approach, considering the hydraulic connectivity f. The simulations showed that the equivalent permeability and the shape factor are strongly influenced by the hydraulic connectivity (f) in synthetic scenarios for X and Y-node topological patterns, which showed the higher value for f (0.81) and more expressive values for upscaled permeability (kx-node = 0.1151 and ky-node = 0.1153) and shape factor (25.6 and 14.5), respectively. We have shown that the analytical methods are not efficient for estimating the equivalent permeability of the fractured medium, including when these methods were corrected using topological aspects.

Oil production in offshore regions involves the transportation of oil and gas in submarine pipelines, which are vulnerable to geological processes triggered by subsurface oil production like fault reactivation. The fault reactivation process can lead to phenomena that impact the seabed, like subsidence and fluid exudation, and can trigger instability of submarine slopes, which can result in environmental and economic damage. The present work addresses a coupled hydromechanical numerical modeling of a hypothetical case involving fault reactivation caused by oil reservoir production and its impact on an overlying submarine slope. The hypothetical case was simulated using a finite element model. The case involves a reservoir which is cut by a fault zone that reaches the seabed. The slope instability studied was induced by the injection and production of fluids in the reservoir. The fault zone is assumed to be a sealing region and a geomechanical and pressure field discontinuity within the reservoir. Int this work was used the Mohr-Coulomb elastoplastic model with Perzyna viscoplastic regularization to represent the behavior of the fault zone and the overlying submarine slope. Results showed that the fault reactivation, caused by the reservoir production, developed shear stress and shear plastic strain along the fault and through the submarine slope, causing horizontal and vertical displacements in the slope mass and acting as a trigger factor for slope stability. Pore pressure increase at the bottom of the slope structure correlated with the injection pressure artificially increased into the reservoir.