A quasi steady state method for solving transient Darcy flow in complex 3D fractured networks

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.

[1]  Yu-Shu Wu.  Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs , 2000 .

[2]  Wilkinson,et al.  Hydrodynamic dispersion in network models of porous media. , 1986, Physical review letters.

[3]  Brian Berkowitz,et al.  PERCOLATION THEORY AND ITS APPLICATION TO GROUNDWATER HYDROLOGY , 1993 .

[4]  J. Erhel,et al.  A mixed hybrid Mortar method for solving flow in discrete fracture networks , 2010 .

[5]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

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

[7]  M. Panfilov.  Averaged model-type transition in flows through multiple heterogeneous porous media , 1994 .

[8]  Y. Yortsos,et al.  Application of Fractal Geometry to the Study of Networks of Fractures and Their Pressure Transient , 1995 .

[9]  A. M. Saidi,et al.  Simulation of Naturally Fractured Reservoirs , 1983 .

[10]  J. F. Mccarthy,et al.  Effective permeability of sandstone-shale reservoirs by a random walk method , 1990 .

[11]  Sylvain Sarda,et al.  Hydraulic Characterization of Fractured Reservoirs: Simulation on Discrete Fracture Models , 2002 .

[12]  P. Lemonnier,et al.  Simulation of Naturally Fractured Reservoirs. State of the Art - Part 2 – Matrix-Fracture Transfers and Typical Features of Numerical Studies , 2010 .

[13]  Gudmundur S. Bodvarsson,et al.  A numerical dual‐porosity model with semianalytical treatment of fracture/matrix flow , 1993 .

[14]  Bernard Bourbiaux,et al.  Fractured Reservoir Simulation: a Challenging and Rewarding Issue , 2010 .

[15]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[16]  T. Estebenet,et al.  Up-scaling Flow in Fractured Media: Equivalence Between the Large Scale Averaging Theory and the Continuous Time Random Walk Method , 2001 .

[17]  Jean-Raynald de Dreuzy,et al.  Use of power averaging for quantifying the influence of structure organization on permeability upscaling in on‐lattice networks under mean parallel flow , 2010 .

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

[19]  Michel Quintard,et al.  Transport in chemically and mechanically heterogeneous porous media. I: Theoretical development of region-averaged equations for slightly compressible single-phase flow , 1996 .

[20]  U. Hornung.  Homogenization and porous media , 1996 .

[21]  Bernard Bourbiaux,et al.  Hydraulic Characterization of Faults and Fractures Using a Dual Medium Discrete Fracture Network Simulator , 2004 .

[22]  Adler,et al.  Geometry of simulated fractures. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[23]  Benoît Nœtinger,et al.  Up-Scaling of Double Porosity Fractured Media Using Continuous-Time Random Walks Methods , 2000 .

[24]  Michel Quintard,et al.  Quasi-steady two-equation models for diffusive transport in fractured porous media: large-scale properties for densely fractured systems , 2001 .

[25]  P. Renard,et al.  Calculating equivalent permeability: a review , 1997 .

[26]  André Fourno,et al.  Smeared fractures: a promising approach to model transfers in fractured media , 2004 .

[27]  Todd Arbogast,et al.  Derivation of the double porosity model of single phase flow via homogenization theory , 1990 .

[28]  Sylvain Sarda,et al.  A New Approach of Fractured Reservoirs , 1998 .

[29]  Michael Ian Shamos,et al.  Computational geometry: an introduction , 1985 .

[30]  Pierre M. Adler,et al.  Pressure drawdown well tests in fractured porous media , 2003 .

[31]  T. N. Narasimhan,et al.  MINC: AN APPROACH FOR ANALYZING TRANSPORT IN STRONGLY HETEROGENEOUS SYSTEMS , 1988 .

[32]  Michel Quintard,et al.  A two-phase flow simulation of a fractured reservoir using a new fissure element method , 2001 .

[33]  Claude Boutin,et al.  Macroscopic modeling of double-porosity reservoirs , 1996 .

[34]  Yalchin Efendiev,et al.  Multiscale finite element methods for high-contrast problems using local spectral basis functions , 2011, J. Comput. Phys..

[35]  Brian Berkowitz,et al.  Mass transfer at fracture intersections: An evaluation of mixing models , 1994 .

[36]  J. F. Mccarthy,et al.  Continuous-time random walks on random media , 1993 .

[37]  Michel Quintard,et al.  Transport in chemically and mechanically heterogeneous porous media. II: Comparison with numerical experiments for slightly compressible single-phase flow , 1996 .

[38]  Felipe Pereira,et al.  A parallelizable method for two‐phase flows in naturally‐fractured reservoirs , 1997 .

[39]  Wilkinson,et al.  Transport and dispersion in random networks with percolation disorder. , 1988, Physical review. A, General physics.

[40]  Martin S. Fridson,et al.  Trends , 1948, Bankmagazin.

[41]  Bernard Bourbiaux,et al.  Simulation of Naturally Fractured Reservoirs. State of the Art - Part 2 – Matrix-Fracture Transfers and Typical Features of Numerical Studies , 2010 .

[42]  Benoit Noetinger,et al.  Scaling Up Matrix-Fracture Transfers in Dual-Porosity Models: Theory and Application , 1999 .

[43]  S. P. Neuman,et al.  Trends, prospects and challenges in quantifying flow and transport through fractured rocks , 2005 .

[44]  B. Nœtinger,et al.  A Direct Determination of the Transient Exchange Term of Fractured Media Using a Continuous Time Random Walk Method , 2001 .

[45]  Y. Yortsos,et al.  Pressure Transient Analysis of Fractal Reservoirs , 1990 .