A new equi-dimensional fracture model using polyhedral cells for microseismic data sets

Abstract We present a method for modeling flow in porous media in the presence of complex fracture networks. The approach utilizes the Mimetic Finite Difference (MFD) method. We employ a novel equi-dimensional approach for meshing fractures. By using polyhedral cells we avoid the common challenge in equi-dimensional fracture modeling of creating small cells at the intersection point. We also demonstrate how polyhedra can mesh complex fractures without introducing a large number of cells. We use polyhedra and the MFD method a second time for embedding fracture boundaries in the matrix domain using a “cut-cell” paradigm. The embedding approach has the advantage of being simple and localizes irregular cells to the area around the fractures. It also circumvents the need for conventional mesh generation, which can be challenging when applied to complex fracture geometries. We present numerical results confirming the validity of our approach for complex fracture networks and for different flow models. In our first example, we compare our method to the popular dual-porosity technique. Our second example compares our method with directly meshed fractures (single-porosity) for two-phase flow. The third example demonstrates two-phase flow for the case of intersecting ellipsoid fractures in three-dimensions, which are typical in microseismic analysis of fractures. Finally, we demonstrate our method on a two-dimensional fracture network produced from microseismic field data.

[1]  Christopher A. Wright,et al.  Integrating Fracture Mapping Technologies to Optimize Stimulations in the Barnett Shale , 2002 .

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

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

[4]  S. Srinivasan,et al.  Stochastic simulation of fracture strikes using seismic anisotropy induced velocity anomalies , 2009 .

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

[6]  K. Aziz,et al.  Matrix-fracture transfer shape factors for dual-porosity simulators , 1995 .

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

[8]  Stefano Berrone,et al.  A PDE-Constrained Optimization Formulation for Discrete Fracture Network Flows , 2013, SIAM J. Sci. Comput..

[9]  Mary F. Wheeler,et al.  A lubrication fracture model in a poro-elastic medium , 2015 .

[10]  K. H. Coats,et al.  Implicit Compositional Simulation of Single-Porosity and Dual-Porosity Reservoirs , 1989 .

[11]  Mary F. Wheeler,et al.  A generalized Mimetic Finite Difference Method and Two-Point Flux schemes over Voronoi diagrams , 2017 .

[12]  Stefano Berrone,et al.  A Parallel Solver for Large Scale DFN Flow Simulations , 2015, SIAM J. Sci. Comput..

[13]  Brian Mirtich,et al.  Fast and Accurate Computation of Polyhedral Mass Properties , 1996, J. Graphics, GPU, & Game Tools.

[14]  Maria José Jurado,et al.  Two-dimensional geostatistical modeling and prediction of the fracture system in the Albala Granitic Pluton, SW Iberian Massif, Spain , 2001 .

[15]  H. Kazemi,et al.  NUMERICAL SIMULATION OF WATER-OIL FLOW IN NATURALLY FRACTURED RESERVOIRS , 1976 .

[16]  R. Eymard,et al.  A UNIFIED APPROACH TO MIMETIC FINITE DIFFERENCE, HYBRID FINITE VOLUME AND MIXED FINITE VOLUME METHODS , 2008, 0812.2097.

[17]  Anna Scotti,et al.  Optimal Techniques to Simulate Flow in Fractured Reservoir , 2014 .

[18]  Thomas Viard,et al.  Unstructured Cut-Cell Grids for Modeling Complex Reservoirs , 2014 .

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

[20]  Georg Kosakowski,et al.  Flow pattern variability in natural fracture intersections , 1999 .

[21]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

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

[23]  Konstantin Lipnikov,et al.  Convergence of the Mimetic Finite Difference Method for Diffusion Problems on Polyhedral Meshes , 2005, SIAM J. Numer. Anal..

[24]  Amy Henderson,et al.  The ParaView Guide: A Parallel Visualization Application , 2004 .

[25]  Rainer Helmig,et al.  Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system , 2009 .

[26]  Brian Berkowitz,et al.  Effects of junction transfer characteristics on transport in fracture networks , 2001 .

[27]  Ivan Yotov,et al.  Local flux mimetic finite difference methods , 2009, Numerische Mathematik.

[28]  Stein Krogstad,et al.  Open-source MATLAB implementation of consistent discretisations on complex grids , 2012, Computational Geosciences.

[29]  R. Hinkelmann,et al.  Adaptive methods for the equidimensional modelling of flow and transport processes in fractured aquifers , 2002 .

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

[31]  Gergina Pencheva,et al.  Modeling Multiphase Flow With Nonplanar Fractures , 2013, ANSS 2013.

[32]  Stefano Berrone,et al.  Towards effective flow simulations in realistic discrete fracture networks , 2016, J. Comput. Phys..