Towards effective flow simulations in realistic Discrete Fracture Networks

We focus on the simulation of underground flow in fractured media, modeled by means of Discrete Fracture Networks. Focusing on a new recent numerical approach proposed by the authors for tackling the problem avoiding mesh generation problems, we further improve the new family of methods making a step further towards effective simulations of large, multi-scale, heterogeneous networks. Namely, we tackle the imposition of Dirichlet boundary conditions in weak form, in such a way that geometrical complexity of the DFN is not an issue; we effectively solve DFN problems with fracture transmissivities spanning many orders of magnitude and approaching zero; furthermore, we address several numerical issues for improving the numerical solution also in quite challenging networks

Simulation of subseismic joint and fault networks using a heuristic mechanical model

Flow simulations of fractured and faulted reservoirs require representation of subseismic structures about which subsurface data are limited. We describe a method for simulating fracture growth that is mechanically based but heuristic, allowing for realistic modelling of fracture networks with reasonable run times. The method takes a triangulated meshed surface as input, together with an initial stress field. Fractures initiate and grow based on the stress field, and the growing fractures relieve the stress in the mesh. We show that a wide range of bedding-plane joint networks can be modelled simply by varying the distribution and anisotropy of the initial stress field. The results are in good qualitative agreement with natural joint patterns. We then apply the method to a set of parallel veins and demonstrate how the variations in thickness of the veins can be represented. Lastly, we apply the method to the simulation of normal fault patterns on salt domes. We derive the stress field on the bedding surface using the horizon curvature. The modelled fault network shows both radial and concentric faults. The new method provides an effective means of modelling joint and fault networks that can be imported to the flow simulator

Dual virtual element method for discrete fractures networks

Discrete fracture networks is a key ingredient in the simulation of physical processes which involve fluid flow in the underground, when the surrounding rock matrix is considered impervious. In this paper we present two different models to compute the pressure field and Darcy velocity in the system. The first allows a normal flow out of a fracture at the intersections, while the second grants also a tangential flow along the intersections. For the numerical discretization, we use the mixed virtual finite element method as it is known to handle grid elements of, almost, any arbitrary shape. The flexibility of the discretization allows us to loosen the requirements on grid construction, and thus significantly simplify the flow discretization compared to traditional discrete fracture network models. A coarsening algorithm, from the algebraic multigrid literature, is also considered to further speed up the computation. The performance of the method is validated by numerical experiments

Fine Scale Simulation of Fractured Reservoirs: Applications and Comparison

Imperial Users onl

From invasion percolation to flow in rock fracture networks

The main purpose of this work is to simulate two-phase flow in the form of immiscible displacement through anisotropic, three-dimensional (3D) discrete fracture networks (DFN). The considered DFNs are artificially generated, based on a general distribution function or are conditioned on measured data from deep geological investigations. We introduce several modifications to the invasion percolation (MIP) to incorporate fracture inclinations, intersection lines, as well as the hydraulic path length inside the fractures. Additionally a trapping algorithm is implemented that forbids any advance of the invading fluid into a region, where the defending fluid is completely encircled by the invader and has no escape route. We study invasion, saturation, and flow through artificial fracture networks, with varying anisotropy and size and finally compare our findings to well studied, conditioned fracture networks.Comment: 18 pages, 10 figure

Dynamic development of hydrofracture

Many natural examples of complex joint and vein networks in layered sedimentary rocks are hydrofractures that form by a combination of pore fluid overpressure and tectonic stresses. In this paper, a two-dimensional hybrid hydro-mechanical formulation is proposed to model the dynamic development of natural hydrofractures. The numerical scheme combines a discrete element model (DEM) framework that represents a porous solid medium with a supplementary Darcy based pore-pressure diffusion as continuum description for the fluid. This combination yields a porosity controlled coupling between an evolving fracture network and the associated hydraulic field. The model is tested on some basic cases of hydro-driven fracturing commonly found in nature, e.g., fracturing due to local fluid overpressure in rocks subjected to hydrostatic and nonhydrostatic tectonic loadings. In our models we find that seepage forces created by hydraulic pressure gradients together with poroelastic feedback upon discrete fracturing play a significant role in subsurface rock deformation. These forces manipulate the growth and geometry of hydrofractures in addition to tectonic stresses and the mechanical properties of the porous rocks. Our results show characteristic failure patterns that reflect different tectonic and lithological conditions and are qualitatively consistent with existing analogue and numerical studies as well as field observations. The applied scheme is numerically efficient, can be applied at various scales and is computational cost effective with the least involvement of sophisticated mathematical computation of hydrodynamic flow between the solid grains

Upscaling of a dual-permeability Monte Carlo simulation model for contaminant transport in fractured networks by genetic algorithm parameter identification

International audienceThe transport of radionuclides in fractured media plays a fundamental role in determining the level of risk offered by a radioactive waste repository in terms of expected doses. Discrete Fracture Networks (DFN) methods can provide detailed solutions to the problem of modeling the contaminant transport in fractured media. However, within the framework of the performance assessment (PA) of radioactive waste repositories, the computational efforts required are not compatible with the repeated calculations that need to be performed for the probabilistic uncertainty and sensitivity analyses of PA. In this paper, we present a novel upscaling approach, which consists in computing the detailed numerical fractured flow and transport solutions on a small scale and use the results to derive the equivalent continuum parameters of a lean, one-dimensional Dual-Permeability, Monte Carlo Simulation (DPMCS) model by means of a Genetic Algorithm search. The proposed upscaling procedure is illustrated with reference to a realistic case study of migration taken from literature