The talk is concerned with a novel discrete fracture model for single-phase flow in aquifers with highly permeable fractures, which are assumed to be filled with a porous medium different from the porous medium in the surrounding matrix rock. The fractures are treated as a ($n-1$)-dimensional interface in a $n$-dimensional domain. The finite element method couples the flow in the fracture and...
Numerical investigations of subsurface flow in fractured porous media provide information about properties connected to underground matter and heat transport just as characteristics of fluid underground storage capacity. Many diffusion-based models in the literature precisely describe subsurface flow. Nevertheless, pronounced hydro-mechanically-coupled phenomena like inverse water level...
We propose a new formulation based on discontinuous Galerkin (DG) methods in their generalization to polytopic grids for the simulation of flows in fractured porous media. The method that we propose is very flexible from the geometrical point of view, being able to handle meshes made of arbitrarily shaped elements, with edges/faces that may be in arbitrary number (potentially unlimited) and...
Discrete fracture network (DFN) models explicitly use fracture geometry and network topology to simulate flow and transport through fractured systems. Recent advances in high performance computing have opened the door for flow and transport simulations in large explicit three-dimensional DFN. However, this increase in model fidelity and network size comes at a huge computational cost because...
An Algebraic Dynamic Multilevel (ADM) method [1, 2] for fully implicit simulations of multiphase flow in heterogeneous fractured natural porous media is presented. The fine-scale fully-implicit (FIM) system is obtained following the Embedded Discrete Fracture Modelling (EDFM) [3] approach. A set of nested coarse grids at different resolutions (or levels) is constructed independently for each...
Modelling flows within faults is crucial for various applications such as the control of faults on overpressure development or hydrocarbon migration
in sedimentary basins, the recovery of hydrocarbon components, subsurface gas storage, and
the appraisal of the risk of groundwater contamination following an underground nuclear waste disposal.
Faults can be characterized as extended...
We present a mathematical model of flow and solid mechanics in saturated fractured porous media based on the Biot poroelasticity. The fractures are treated as lower-dimensional manifolds on which the system of equations is projected onto the tangent space and coupled to the surrounding through interface conditions. The model can describe porous fractures, fractures filled only by a liquid and...
We consider in this work, the extension to two-phase Darcy flows of the discontinuous pressure models in which the (d-1)-dimensional flow in the fractures is coupled with the d-dimensional flow in the matrix.The model accounts accurately for gravity effects inside the fractures, for discontinuous capillary pressure curves at the matrix fracture (mf) interfaces and for both drains and barriers....
The use of energized fluids in hydraulic fracturing helps minimize formation damage and enhance well productivity especially in water sensitive reservoirs. Many current fracking simulators use single phase and incompressible flow assumptions to model the fracturing fluids. The thermal and compositional effects of using energized fluids on fracture topology are neglected by such simulators.
A...
We are concerned with the application of nonoverlapping domain decomposition methods on reduced fracture models. We provide for this problem a new DD algorithm in which the global problem is reduced to a mortar problem posed on the fractures that is solved by an iterative solver. Classical iterative methods are based on the CG or optimized Schwarz waveform methods which require...
In this work, we present a parallel MPI software, called GENFIELD, to generate stationary Gaussian random fields, based on the circulant embedding method [7, 1, 4].
The advantage of the circulant embedding method over existing methods (Cholesky factorization of the covariance matrix, Karhunen Loève series expansion) is its computational efficiency. It relies on Discrete Fourier Transform that...
The metallic foams are a novel possibility to impact significantly the design of materials taking into account the wide technological applications [1,2]. In this work the procedure to manufacture metallic foams with controlled porosity is presented. A high pressure cell is used to submit the metallic foams to reservoir conditions, which imply high pressures and high temperatures. The injection...
In this talk, experimental results of structural stability of synthetic rocks and metal foams[1] subjected to reservoir conditions are presented. By controlling the injection of water and mineral oils ona confinement cell, the samples are subjected to conditions of high pressure (25KPSI) and high temperature (200°C). By means of micrography techniques the mechanical deformation of foam layers...
Naturally Fractured Reservoirs (NFR) are usually multiscale in nature and exhibit power law length distributions which do not possess any characteristic length scale, rendering the use of continuum methods such as the dual porosity model invalid due to the non-existence of Representative Elementary Volumes (REV) (Berkowitz, 2002). This necessitates the adoption of hybrid models that represent...
