FracKfinder: A MATLAB Toolbox for Computing Three-Dimensional Hydraulic Conductivity Tensors for Fractured Porous Media

Ground Water ◽  
2018 ◽  
Vol 57 (1) ◽  
pp. 75-80 ◽  
Author(s):  
Nathan L. Young ◽  
Jacqueline E. Reber ◽  
William W. Simpkins
2018 ◽  
Vol 97 (1) ◽  
Author(s):  
C. Mezon ◽  
V. V. Mourzenko ◽  
J.-F. Thovert ◽  
R. Antoine ◽  
F. Fontaine ◽  
...  

Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-13 ◽  
Author(s):  
Chi Yao ◽  
Chen He ◽  
Jianhua Yang ◽  
Qinghui Jiang ◽  
Jinsong Huang ◽  
...  

An original 3D numerical approach for fluid flow in fractured porous media is proposed. The whole research domain is discretized by the Delaunay tetrahedron based on the concept of node saturation. Tetrahedral blocks are impermeable, and fluid only flows through the interconnected interfaces between blocks. Fractures and the porous matrix are replaced by the triangular interface network, which is the so-called equivalent matrix-fracture network (EMFN). In this way, the three-dimensional seepage problem becomes a two-dimensional problem. The finite element method is used to solve the steady-state flow problem. The big finding is that the ratio of the macroconductivity of the whole interface network to the local conductivity of an interface is linearly related to the cubic root of the number of nodes used for mesh generation. A formula is presented to describe this relationship. With this formula, we can make sure that the EMFN produces the same macroscopic hydraulic conductivity as the intact rock. The approach is applied in a series of numerical tests to demonstrate its efficiency. Effects of the hydraulic aperture of fracture and connectivity of the fracture network on the effective hydraulic conductivity of fractured rock masses are systematically investigated.


Author(s):  
Alessio Fumagalli ◽  
Eirik Keilegavlen

The accurate description of fluid flow and transport in fractured porous media is of paramount importance to capture the macroscopic behavior of an oil reservoir, a geothermal system, or a CO2 sequestration site, to name few applications. The construction of accurate simulation models for flow in fractures is challenging due to the high ratio between a fracture’s length and width. In this paper, we present a mixed-dimensional Darcy problem which can represent the pressure and Darcy velocity in all the dimensions, i.e. in the rock matrix, in the fractures, and in their intersections. Moreover, we present a mixed-dimensional transport problem which, given the Darcy velocity, describes advection of a passive scalar into the fractured porous media. The approach can handle both conducting and blocking fractures. Our computational grids are created by coarsening of simplex tessellations that conform to the fracture’s surfaces. A suitable choice of the discrete approximation of the previous model, by virtual finite element and finite volume methods, allows us to simulate complex problems with a good balance of accuracy and computational cost. We illustrate the performance of our method by comparing to benchmark studies for two-dimensional fractured porous media, as well as a complex three-dimensional fracture geometry.


Sign in / Sign up

Export Citation Format

Share Document