ASYMPTOTIC AND NUMERICAL MODELLING OF FLOWS IN FRACTURED POROUS MEDIA

This study concerns some asymptotic models used to compute the flow outside and inside fractures in a bidimensional porous medium. The flow is governed by the Darcy law both in the fractures and in the porous matrix with large discontinuities in the permeability tensor. These fractures are supposed to have a small thickness with respect to the macroscopic length scale, so that we can asymptotically reduce them to immersed polygonal fault interfaces and the model finally consists in a coupling between a 2D elliptic problem and a 1D equation on the sharp interfaces modelling the fractures. A cell-centered finite volume scheme on general polygonal meshes fitting the interfaces is derived to solve the set of equations with the additional differential transmission conditions linking both pressure and normal velocity jumps through the interfaces. We prove the convergence of the FV scheme for any set of data and parameters of the models and derive existence and uniqueness of the solution to the asymptotic models proposed. The models are then numerically experimented for highly or partially immersed fractures. Some numerical results are reported showing different kinds of flows in the case of impermeable or partially/highly permeable fractures. The influence of the variation of the aperture of the fractures is also investigated. The numerical solutions of the asymptotic models are validated by comparing them to the solutions of the global Darcy model or to some analytic solutions.

[1]  P. Angot FINITE VOLUME METHODS FOR NON SMOOTH SOLUTION OF DIFFUSION MODELS;: APPLICATION TO IMPERFECT CONTACT PROBLEMS , 1999 .

[2]  F. Hermeline,et al.  Approximation of diffusion operators with discontinuous tensor coefficients on distorted meshes , 2003 .

[3]  F. Boyer,et al.  Discrete duality finite volume schemes for Leray−Lions−type elliptic problems on general 2D meshes , 2007 .

[4]  Thierry Gallouët,et al.  H-Convergence and Numerical Schemes for Elliptic Problems , 2003, SIAM J. Numer. Anal..

[5]  P. Adler,et al.  Darcy flow around a two-dimensional permeable lens , 2006 .

[6]  Vincent Martin,et al.  Modeling Fractures and Barriers as Interfaces for Flow in Porous Media , 2005, SIAM J. Sci. Comput..

[7]  Philippe Angot A model of fracture for elliptic problems with flux and solution jumps , 2003 .

[8]  Monique Dauge,et al.  Polynomials in the Sobolev World , 2007 .

[9]  Christine Bernardi,et al.  Compatibilité de traces aux arêtes et coins d'un polyèdre , 2000 .

[10]  P. Fabrie,et al.  Computation of large-scale parameters for dispersion in fissured porous medium using finite-volume method , 2001 .

[11]  Pierre M. Adler,et al.  Effective permeability of fractured porous media in steady state flow , 2003 .

[12]  J. Bear,et al.  1 – Modeling Flow and Contaminant Transport in Fractured Rocks , 1993 .

[13]  R. Eymard,et al.  Finite Volume Methods , 2019, Computational Methods for Fluid Dynamics.

[14]  R. Helmig,et al.  A mixed-dimensional finite volume method for two-phase flow in fractured porous media , 2006 .

[15]  Pierre M. Adler,et al.  Fractures and Fracture Networks , 1999 .

[16]  Franck Boyer,et al.  Finite Volume Method for 2D Linear and Nonlinear Elliptic Problems with Discontinuities , 2008, SIAM J. Numer. Anal..

[17]  B. Berkowitz Characterizing flow and transport in fractured geological media: A review , 2002 .

[18]  Sally M. Benson,et al.  Dynamics of fluids in fractured rock , 2000 .

[19]  Michel Quintard,et al.  Implementation of a Finite-Volume Method for the Determination of Effective Parameters in Fissured Porous Media , 2000 .

[20]  Pascal Omnes,et al.  A FINITE VOLUME METHOD FOR THE LAPLACE EQUATION ON ALMOST ARBITRARY TWO-DIMENSIONAL GRIDS , 2005 .

[21]  Robert Eymard,et al.  H-convergence and numerical schemes for elliptic equations SIAM Journal on Numerical Analysis , 2000 .

[22]  Chin-Fu Tsang,et al.  Flow and Contaminant Transport in Fractured Rock , 1993 .

[23]  P. Grisvard Elliptic Problems in Nonsmooth Domains , 1985 .