A numerical method for two-phase flow in fractured porous media with non-matching grids

Abstract We propose a novel computational method for the efficient simulation of two-phase flow in fractured porous media. Instead of refining the grid to capture the flow along the faults or fractures, we represent the latter as immersed interfaces, using a reduced model for the flow and suitable coupling conditions. We allow for non matching grids between the porous matrix and the fractures to increase the flexibility of the method in realistic cases. We employ the extended finite element method for the Darcy problem and a finite volume method that is able to handle cut cells and matrix-fracture interactions for the saturation equation. Moreover, we address through numerical experiments the problem of the choice of a suitable numerical flux in the case of a discontinuous flux function at the interface between the fracture and the porous matrix. A wrong approximate solution of the Riemann problem can yield unphysical solutions even in simple cases.

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

[2]  Alessio Fumagalli,et al.  A Reduced Model for Flow and Transport in Fractured Porous Media with Non-matching Grids , 2013 .

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

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

[5]  Hussein Hoteit,et al.  An efficient numerical model for incompressible two-phase flow in fractured media , 2008 .

[6]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[7]  John D. Towers,et al.  L¹ STABILITY FOR ENTROPY SOLUTIONS OF NONLINEAR DEGENERATE PARABOLIC CONVECTION-DIFFUSION EQUATIONS WITH DISCONTINUOUS COEFFICIENTS , 2003 .

[8]  Ruben Juanes,et al.  CO2 migration in saline aquifers. Part 1. Capillary trapping under slope and groundwater flow , 2010, Journal of Fluid Mechanics.

[9]  G. Chavent.  Mathematical models and finite elements for reservoir simulation , 1986 .

[10]  Jérôme Jaffré,et al.  Domain Decomposition for Some Transmission Problems in Flow in Porous Media , 2000 .

[11]  Laurent Trenty,et al.  A benchmark study on problems related to CO2 storage in geologic formations , 2009 .

[12]  Mark A. Herkommer,et al.  Data-volume reduction of data gathered along lines using the correlation coefficient to determine breakpoints , 1985 .

[13]  G. Chavent,et al.  A finite element simulator for incompressible two-phase flow , 1987 .

[14]  Siam Staff,et al.  Godunov-Type Methods for Conservation Laws with a Flux Function Discontinuous in Space , 2004 .

[15]  E. F. Kaasschieter.  Solving the Buckley–Leverett equation with gravity in a heterogeneous porous medium , 1999 .

[16]  Jérôme Jaffré,et al.  Upstream differencing for multiphase flow in reservoir simulation , 1991 .

[17]  Emil M. Constantinescu,et al.  Multirate Timestepping Methods for Hyperbolic Conservation Laws , 2007, J. Sci. Comput..

[18]  N. W. Lanfredi,et al.  HP 67/97 calculator waves application programs , 1987 .

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

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

[21]  Rainer Helmig,et al.  CO2 leakage through an abandoned well: problem-oriented benchmarks , 2007 .

[22]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[23]  P. Hansbo,et al.  An unfitted finite element method, based on Nitsche's method, for elliptic interface problems , 2002 .

[24]  Jérôme Jaffré,et al.  On the upstream mobility scheme for two-phase flow in porous media , 2009, ArXiv.

[25]  Ivar Aavatsmark,et al.  Errors in the upstream mobility scheme for countercurrent two-phase flow in heterogeneous porous media , 2012, Computational Geosciences.

[26]  Alessio Fumagalli,et al.  Numerical modelling of multiphase subsurface flow in the presence of fractures , 2011 .