This paper reports on analyses of fluid exchange between fractures in the simulation of flow through fractured porous media. Flow in the porous media is modeled with traditional Darcy’s equations and the coupling between flow in the porous media and fractures is based on the conceptual Discrete-Fracture-Matrix (DFM) representation, where the fractures are idealized as lower-dimensional elements at the interface of matrix elements. The DFNMesh algorithm is adopted to generate the Finite Element meshes and a novel methodology to handle overlapping fractures in the context of Mixed Finite Element methods is explored. Flux approximation with H(div)-conforming spaces are adopted which are particularly appealing for this analysis because of its inherent properties such as local mass conservation and strong divergence-free enforcement for incompressible flows. The analyses are carried out using simple two-fracture domains and a benchmark problem with eight fractures.
Month: January 2026
Few-shot and continuous online learning for forecasting in the energy industry
Experimental Investigation of Polymer Degradation and Its Effects on Electrical Submersible Pump Operation
Challenges and suggestions for defining electrofacies models: The problem with well-log resolution and the resulting “shoulder bed effects” applied to a carbonate reservoir
Unravelling the origin of reworked deposits in Aptian lacustrine carbonate reservoirs of the Santos Basin, SE Brazil
Chapter Five – Extending h adaptivity with refinement patterns
This contribution introduces the idea of refinement patterns for the generation of optimal meshes in the context of the Finite Element Method. The main idea is to generate a library of possible patterns on which elements can be refined and use this library to inform an h adaptive code on how to handle complex refinements in regions of interest. There are no restrictions on the type of elements that can be refined, and the patterns can be generated for any element type. The main advantage of this approach is that it allows for the generation of optimal meshes in a systematic way where, even if a certain pattern is not available, it can easily be included through a simple text file with nodes and sub-elements. The contribution presents a detailed methodology for incorporating refinement patterns into h adaptive Finite Element Method codes and demonstrates the effectiveness of the approach through mesh refinement of problems with complex geometries.
DFNMesh: Finite element meshing for discrete fracture matrix models
This paper presents a novel methodology for robust Finite Element (FE) mesh generation of Discrete Fracture Matrix (DFM) models. The method can handle the complete multi-dimensional domain from the tridimensional porous rock matrix, through fracture surfaces, down to open curves for fracture–fracture intersections. The accompanying open-source code is written in modern C++ with a JSON interface and largely relies on two state-of-the-art FE libraries: NeoPZ and Gmsh. Starting from a user-defined coarse mesh, fractures are sequentially read as convex polygons and introduced without distortions to the coarse geometry. The main steps involve: intersect edges by checking for nodes on opposite sides of the fracture plane, extend intersections from edges to faces, coalesce intersections to closest existing nodes (given a tolerance), refine interface elements to conform to the fracture, identify subsets of fracture surface, mesh the surface, and locate boundaries and intersections where they arise. Finally, the space around fractures is filled with the fine-scale unstructured mesh, which is kept conformal. The robustness of the implementation is derived from the consistent background of well-defined and simple premises like convexity and side-specific element neighborhood. Results show that the proposed technique can construct adequate 3D DFM grids, while still giving users freedom to adjust between geometrical fidelity and mesh quality through more aggressive feature rejection.
Core scale investigation of fluid flow in the heterogeneous porous media based on X-ray computed tomography images: Upscaling and history matching approaches
In this paper, experiments and simulations were performed on outcrop samples from Lagoa Salgada in Rio de Janeiro, Brazil, as a possible analog to one of the most typical Brazilian Pre-salt carbonate reservoirs rocks. The rocks were microbial carbonates where plugs comprising two main facies were sampled, simplified as fine-grained and vugular facies. The plugs were utilized to study the impact of pore geometry with both experimental and simulation approaches on recovery factor, saturation profile, and relative permeability estimations. To provide direct visualization of the geometry, description of pore structure, and calculation of concentration profiles, computed tomography (CT) imaging was integrated with experimental measurements of petrography and core flooding. The injection of two pore volumes of formation water resulted in a recovery factor between 28 and 34 percent for the plug samples. Furthermore, based on porosity generated by dry and wet CT, as well as saturation profiles resulting from CT data collected during drainage and imbibition processes along the length of the plugs, it is revealed that the distribution of these properties was diverse and heterogeneous. An algorithm was used to process the 2D tomography images of the samples to remove the region related to the exterior parts. The images were then stacked to create a 3D fine-scale grid to simulate the porous media and the fluid flow by applying rules for the segmentation of rock types, porosity, and permeability estimations of each grid block. Course-scale grids were created by applying upscaling techniques to reduce computation time. Simulated produced fluid cuts for different upscaled models were compared with the experimental results from core flooding. A history-matching technique was then applied to match experimental and simulation results, calculating the relative permeability of two main defined facies and creating an updated model capable of assessing past and present performance and future forecasting. Since relative permeability is essential for accurate simulation, estimating these curves in the heterogeneous pre-salt reservoir considering different facies, greatly influences reasonable future prediction performance and the ability to make informed operational decisions.
- « Previous
- 1
- 2
- 3
- 4