Unified pipe network method for simulation of water flow in fractured porous rock

Abstract Rock masses are often conceptualized as dual-permeability media containing fractures or fracture networks with high permeability and porous matrix that is less permeable. In order to overcome the difficulties in simulating fluid flow in a highly discontinuous dual-permeability medium, an effective unified pipe network method is developed, which discretizes the dual-permeability rock mass into a virtual pipe network system. It includes fracture pipe networks and matrix pipe networks. They are constructed separately based on equivalent flow models in a representative area or volume by taking the advantage of the orthogonality of the mesh partition. Numerical examples of fluid flow in 2-D and 3-D domain including porous media and fractured porous media are presented to demonstrate the accuracy, robustness, and effectiveness of the proposed unified pipe network method. Results show that the developed method has good performance even with highly distorted mesh. Water recharge into the fractured rock mass with complex fracture network is studied. It has been found in this case that the effect of aperture change on the water recharge rate is more significant in the early stage compared to the fracture density change.

[1]  Satish Karra,et al.  dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport , 2015, Comput. Geosci..

[2]  R. Arnett,et al.  Modelling fluid flow in fractured‐porous rock masses by finite‐element techniques , 1984 .

[3]  Anil V. Rao,et al.  GPOPS-II , 2014, ACM Trans. Math. Softw..

[4]  Khalid Aziz,et al.  Use of Voronoi Grid in Reservoir Simulation , 1994 .

[5]  Roussos Dimitrakopoulos,et al.  An efficient method for discretizing 3D fractured media for subsurface flow and transport simulations , 2011 .

[6]  A. Firoozabadi,et al.  Control‐volume method for numerical simulation of two‐phase immiscible flow in two‐ and three‐dimensional discrete‐fractured media , 2004 .

[7]  P. A. Witherspoon,et al.  Hydromechanical behavior of a deformable rock fracture subject to normal stress , 1981 .

[8]  Chin-Fu Tsang,et al.  A variable aperture fracture network model for flow and transport in fractured rocks , 1992 .

[9]  Ivars Neretnieks,et al.  Fluid flow and solute transport in a network of channels , 1993 .

[10]  F. Ren,et al.  Pipe network model for unconfined seepage analysis in fractured rock masses , 2016 .

[11]  Martin Vohralík,et al.  Numerical simulation of fracture flow with a mixed-hybrid FEM stochastic discrete fracture network model , 2005 .

[12]  S. Priest.  Discontinuity Analysis for Rock Engineering , 1992 .

[13]  Hussein Hoteit,et al.  Compositional Modeling by the Combined Discontinuous Galerkin and Mixed Methods , 2004 .

[14]  Kassem Mustapha,et al.  A New Approach to Simulating Flow in Discrete Fracture Networks with an Optimized Mesh , 2007, SIAM J. Sci. Comput..

[15]  A. Hoch,et al.  Self‐consistency of a heterogeneous continuum porous medium representation of a fractured medium , 2000 .

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

[17]  B. Rivière,et al.  Part II. Discontinuous Galerkin method applied to a single phase flow in porous media , 2000 .

[18]  Derek Elsworth,et al.  A continuum model for coupled stress and fluid flow in discrete fracture networks , 2016, Geomechanics and Geophysics for Geo-Energy and Geo-Resources.

[19]  Jacques Bodin,et al.  A new method for generating a pipe network to handle channelled flow in fractured rocks , 2007 .

[20]  Olaf Kolditz,et al.  Modelling flow and heat transfer in fractured rocks: conceptual model of a 3-D deterministic fracture network , 1995 .

[21]  Guowei Ma,et al.  A Graph-theoretic Pipe Network Method for water flow simulation in discrete fracture networks: GPNM , 2014 .

[22]  J. Noorishad,et al.  An upstream finite element method for solution of transient transport equation in fractured porous media , 1982 .

[23]  D. Viete,et al.  Applicability of the 'cubic law' for non-Darcian fracture flow , 2011 .

[24]  Jun Yao,et al.  A Two-Phase Flow Simulation of Discrete-Fractured Media using Mimetic Finite Difference Method , 2014 .

[25]  Milind Deo,et al.  Finite element, discrete‐fracture model for multiphase flow in porous media , 2000 .

[26]  Ahmed E. Hassan,et al.  On mapping fracture networks onto continuum , 2008 .

[27]  C. Fidelibus,et al.  Derivation of equivalent pipe network analogues for three‐dimensional discrete fracture networks by the boundary element method , 1999 .

[28]  D. Waldren,et al.  A Multicomponent Isothermal System for Efficient Reservoir Simulation , 1983 .

[29]  Christopher C. Pain,et al.  Non-linear regimes of fluid flow in rock fractures , 2004 .

[30]  H. Jourde,et al.  Flow behavior in a dual fracture network , 2002 .

[31]  Unstructured Finite-volume Meshes for Two-dimensional Flow in Variably Saturated Porous Media , 2011 .

[32]  J. Molenaar,et al.  The effect of capillary forces on immiscible two-phase flow in heterogeneous porous media , 1995 .

[33]  I. Fatt.  The Network Model of Porous Media , 1956 .

[34]  Sebastian Bauer,et al.  Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005 , 2007 .

[35]  Gudmundur S. Bodvarsson,et al.  An active fracture model for unsaturated flow and transport in fractured rocks , 1998 .

[36]  G. Marsily,et al.  Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model , 1990 .

[37]  F. T. Tracy,et al.  Clean two‐ and three‐dimensional analytical solutions of Richards' equation for testing numerical solvers , 2006 .

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

[39]  Donald M. Reeves,et al.  Transport of conservative solutes in simulated fracture networks: 1. Synthetic data generation , 2008 .

[40]  Louis J. Durlofsky,et al.  A general gridding, discretization, and coarsening methodology for modeling flow in porous formations with discrete geological features , 2016 .

[41]  Dae-Sung Cheon,et al.  Breakdown pressure and fracture surface morphology of hydraulic fracturing in shale with H2O, CO2 and N2 , 2016 .

[42]  J. Chilès,et al.  Taking into account the complexity of natural fracture systems in reservoir single-phase flow modelling , 2002 .

[43]  Jean-Raynald de Dreuzy,et al.  Flow Simulation in Three-Dimensional Discrete Fracture Networks , 2009, SIAM J. Sci. Comput..

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

[45]  G. Zhao.  Development of the distinct lattice spring model for large deformation analyses , 2014 .

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

[47]  D. Evans,et al.  Laboratory Studies of Gas Flow Through a Single Natural Fracture , 1986 .

[48]  D. Reeves,et al.  Development and Calibration of Dual‐Permeability Flow Models with Discontinuous Fault Networks , 2014 .

[49]  Jeffrey D. Hyman,et al.  Conforming Delaunay Triangulation of Stochastically Generated Three Dimensional Discrete Fracture Networks: A Feature Rejection Algorithm for Meshing Strategy , 2014, SIAM J. Sci. Comput..

[50]  Pierre M. Adler,et al.  Real Porous Media: Local Geometry and Macroscopic Properties , 1998 .

[51]  Salvatore Torquato,et al.  Morphology and physical properties of Fontainebleau sandstone via a tomographic analysis , 1996 .

[52]  M. Karimi-Fard,et al.  Numerical Simulation of Water Injection in Fractured Media Using the Discrete-Fracture Model and the Galerkin Method , 2003 .

[53]  Sarah Vitel,et al.  Unstructured Upgridding and Transmissibility Upscaling for Preferential Flow Paths in 3D Fractured Reservoirs , 2007 .

[54]  Björn Dverstorp,et al.  Discrete fracture network interpretation of field tracer migration in sparsely fractured rock , 1992 .

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

[56]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[57]  Steven L. Bryant,et al.  Network model evaluation of permeability and spatial correlation in a real random sphere packing , 1993 .

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

[59]  Jiannong Fang,et al.  A MLS-Based Lattice Spring Model for Simulating Elasticity of Materials , 2012 .

[60]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[61]  Hang Si,et al.  TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator , 2015, ACM Trans. Math. Softw..

[62]  Keith W. Jones,et al.  Synchrotron computed microtomography of porous media: Topology and transports. , 1994, Physical review letters.

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

[64]  Shucai Li,et al.  A Graph-theoretic Pipe Network Method for water flow simulation in a porous medium: GPNM , 2014 .

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

[66]  Yu-Shu Wu,et al.  A multiple-porosity method for simulation of naturally fractured petroleum reservoirs , 1988 .

[67]  Urban Svensson,et al.  A continuum representation of fracture networks. Part I: Method and basic test cases , 2001 .

[68]  Jacques Bodin,et al.  Simulation and analysis of solute transport in 2D fracture/pipe networks: the SOLFRAC program. , 2007, Journal of contaminant hydrology.

[69]  H. Si.  Three Dimensional Boundary Conforming Delaunay Mesh Generation , 2008 .

[70]  Yang Li,et al.  An efficient embedded discrete fracture model based on mimetic finite difference method , 2016 .