Computational geometry: an introduction
Computational geometry: an introduction
Derivation of the double porosity model of single phase flow via homogenization theory
SIAM Journal on Mathematical Analysis
Microstructure models of porous media
Homogenization and porous media
Multiscale finite element methods for high-contrast problems using local spectral basis functions
Journal of Computational Physics
Synthetic benchmark for modeling flow in 3D fractured media
Computers & Geosciences
An optimization approach for large scale simulations of discrete fracture network flows
Journal of Computational Physics
| Hi-index | 31.45 |
Determining flows and transfers in natural Discrete Fracture Networks (DFNs) receives more and more attention in applied geosciences, from oil and gas industry, geothermal recovery, to CO"2 sequestration issues. Geological characterization of these DFN is becoming more and more accurate. So, accounting for the details of the flow inside the fracture network, and later modelling the transfers between the matrix and the fractures, with the same level of accuracy is an important issue in simulation. In this paper, we focus on the modelling of transient pressure diffusion in the fracture network only. The matrix is assumed to be impervious, and the fluid is slightly compressible. The boundary condition to be applied at the fracture intersections in the limit of vanishing thickness is analysed first. We build a systematic series of approximation schemes that are relevant for characteristic times greater than the typical diffusion time over one fracture. The lowest order approximation appears under the form of solving a transient problem in a resistor/capacitor network. The topology of this network is the same as the network of geometrical intersections between fractures. Formal expressions are given for determining the capacity (or mass) matrix, and the transmissitivity matrix. In practice, these coefficients can be determined by solving steady state local boundary value problems over uncoupled fractures. A mass condensation (or mass lumping) strategy is proposed to simplify the large scale resolution, yielding existing heuristic methods in 2D. An application to the determination of the up-scaled permeability of a DFN is also presented.