Hybrid-dimensional modelling of two-phase flow through fractured porous media with enhanced matrix fracture transmission conditions

In this work, we extend, to two-phase flow, the single-phase Darcy flow model proposed in [26], [12] in which the (d − 1)-dimensional flow in the fractures is coupled with the d-dimensional flow in the matrix. Three types of so called hybrid-dimensional two-phase Darcy flow models are proposed. They all account for fractures acting either as drains or as barriers, since they allow pressure jumps at the matrix-fracture interfaces. The models also permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The three models differ by their transmission conditions at matrix fracture interfaces: while the first model accounts for the nonlinear two-phase Darcy flux conservations, the second and third ones are based on the linear single phase Darcy flux conservations combined with different approximations of the mobilities. We adapt the Vertex Approximate Gradient (VAG) scheme to this problem, in order to account for anisotropy and heterogeneity aspects as well as for applicability on general meshes. Several test cases are presented to compare our hybrid-dimensional models to the generic equi-dimensional model, in which fractures have the same dimension as the matrix, leading to deep insight about the quality of the proposed reduced models.

[1]  Abbas Firoozabadi,et al.  Control-Volume Model for Simulation of Water Injection in Fractured Media: Incorporating Matrix Heterogeneity and Reservoir Wettability Effects , 2007 .

[2]  Alessio Fumagalli,et al.  An Efficient XFEM Approximation of Darcy Flows in Arbitrarily Fractured Porous Media , 2014 .

[3]  Elyes Ahmed,et al.  A reduced fracture model for two-phase flow with different rock types , 2017, Math. Comput. Simul..

[4]  Roland Masson,et al.  Numerical analysis of a two-phase flow discrete fracture model , 2016 .

[5]  Mayya Groza.  Modélisation et discrétisation des écoulements diphasiques en milieux poreux avec réseaux de fractures discrètes , 2016 .

[6]  Mary F. Wheeler,et al.  Decoupling preconditioners in the implicit parallel accurate reservoir simulator (IPARS) , 2001, Numer. Linear Algebra Appl..

[7]  J. Thovert,et al.  Two-phase flow through fractured porous media. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  Raphaèle Herbin,et al.  Small-stencil 3D schemes for diffusive flows in porous media , 2012 .

[9]  Thierry Gallouët,et al.  A model for conductive faults with non-matching grids , 2011, Computational Geosciences.

[10]  L. Durlofsky,et al.  An Efficient Discrete-Fracture Model Applicable for General-Purpose Reservoir Simulators , 2004 .

[11]  Robert Eymard,et al.  Finite volume approximation of a diffusion-dissolution model and application to nuclear waste storage , 2011, Math. Comput. Simul..

[12]  Philippe Angot,et al.  ASYMPTOTIC AND NUMERICAL MODELLING OF FLOWS IN FRACTURED POROUS MEDIA , 2009 .

[13]  Muhammad Sahimi,et al.  Flow and Transport in Porous Media and Fractured Rock: From Classical Methods to Modern Approaches , 1995 .

[14]  Roland Masson,et al.  Gradient schemes for two‐phase flow in heterogeneous porous media and Richards equation , 2014 .

[15]  Anna Scotti,et al.  MIMETIC FINITE DIFFERENCE APPROXIMATION OF FLOWS IN FRACTURED POROUS MEDIA , 2016 .

[16]  Alessio Fumagalli,et al.  A reduced model for Darcy’s problem in networks of fractures , 2014 .

[17]  Hussein Hoteit,et al.  Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures , 2008 .

[18]  Mayur Pal,et al.  Control-volume distributed multi-point flux approximation coupled with a lower-dimensional fracture model , 2015, J. Comput. Phys..

[19]  Mayur Pal,et al.  Control-Volume Distributed Multi-Point Flux Approximation Coupled with a Lower-Dimensional Fracture Model , 2015, ANSS 2015.

[20]  C. D'Angelo,et al.  A mixed finite element method for Darcy flow in fractured porous media with non-matching grids , 2012 .

[21]  Roland Masson,et al.  Vertex Approximate Gradient Scheme for Hybrid Dimensional Two-Phase Darcy Flows in Fractured Porous Media , 2014 .

[22]  Roland Masson,et al.  Domain decomposition for an asymptotic geological fault modeling , 2003 .

[23]  Roland Masson,et al.  Hybrid dimensional two phase flow in fractured porous media , 2016 .

[24]  I. V. Radhakrishna Murthy,et al.  The density difference and generalized programs for two- and three-dimensional gravity modeling , 1990 .

[25]  Roland Masson,et al.  Gradient discretization of hybrid dimensional Darcy flows in fractured porous media , 2014, Numerische Mathematik.

[26]  Roland Masson,et al.  VERTEX CENTRED DISCRETIZATION OF TWO-PHASE DARCY FLOWS ON GENERAL MESHES , 2012 .

[27]  Seong H. Lee,et al.  Efficient Field-Scale Simulation of Black Oil in a Naturally Fractured Reservoir Through Discrete Fracture Networks and Homogenized Media , 2008 .

[28]  Jérôme Jaffré,et al.  Modeling fractures as interfaces for flow and transport in porous media , 2001 .

[29]  R. Eymard,et al.  Discretisation of heterogeneous and anisotropic diffusion problems on general non-conforming meshes. SUSHI: a scheme using stabilisation and hybrid interfaces , 2008, 0801.1430.

[30]  Jérôme Jaffré,et al.  Reduced models for flow in porous media containing faults with discretization using hybrid finite volume schemes , 2015 .

[31]  Jérôme Jaffré,et al.  A discrete fracture model for two-phase flow with matrix-fracture interaction , 2011, ICCS.

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

[33]  Jan M. Nordbotten,et al.  An efficient multi-point flux approximation method for Discrete Fracture-Matrix simulations , 2012, J. Comput. Phys..

[34]  Jeanne Pellerin,et al.  Immiscible two-phase Darcy flow model accounting for vanishing and discontinuous capillary pressures: application to the flow in fractured porous media , 2017, Computational Geosciences.

[35]  Rajat Gupta,et al.  Applied Hydrogeology of Fractured Rocks , 1999 .

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

[37]  Patrick Jenny,et al.  A hierarchical fracture model for the iterative multiscale finite volume method , 2011, J. Comput. Phys..

[38]  AhmedElyes,et al.  A reduced fracture model for two-phase flow with different rock types , 2017 .