Graph-based flow modeling approach adapted to multiscale discrete-fracture-network models.

Fractured rocks are often modeled as multiscale populations of interconnected discrete fractures (discrete fracture network, DFN). Graph representations of DFNs reduce their complexity to their connectivity structure by forming an assembly of nodes connected by links (edges) to which physical properties, like a conductance, can be assigned. In this study, we address the issue of using graphs to predict flow as a fast and relevant substitute to classical DFNs. We consider two types of graphs, whether the nodes represent the fractures (fracture graph) or the intersections between fractures (intersection graph). We introduce an edge conductance expression that accounts for both the portion of the fracture surface that carries flow and fracture transmissivity. We find that including the fracture size in the expression improves the prediction of flow compared to expressions used in previous studies that did not. The two graph types yield very different results. The fracture graph systematically underestimates local flow values. In contrast, the intersection graph overestimates the flow in each fracture because of the connectivity redundancy in fractures with multiple intersections. We address the latter issue by introducing a correction factor into the conductances based on the number of intersections on each fracture. We test the robustness of the proposed conductance model by comparing flow properties of the graph with high-fidelity DFN simulations over a wide range of network types. The good agreement found between the intersection graph and test suite indicates that this representation could be useful for predictive purposes.

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

[2]  Jeffrey D. Hyman,et al.  Emergence of Stable Laws for First Passage Times in Three-Dimensional Random Fracture Networks. , 2019, Physical review letters.

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

[4]  Olivier Bour,et al.  On the connectivity of three‐dimensional fault networks , 1998 .

[5]  J. Hyman,et al.  Multilevel Monte Carlo Predictions of First Passage Times in Three‐Dimensional Discrete Fracture Networks: A Graph‐Based Approach , 2020, Water Resources Research.

[6]  J. Latham,et al.  The use of discrete fracture networks for modelling coupled geomechanical and hydrological behaviour of fractured rocks , 2017 .

[7]  Gowri Srinivasan,et al.  Predictions of first passage times in sparse discrete fracture networks using graph-based reductions. , 2017, Physical review. E.

[8]  J. Selroos,et al.  Overview of hydrogeological safety assessment modeling conducted for the proposed high-level nuclear waste repository site at Forsmark, Sweden , 2014, Hydrogeology Journal.

[9]  Xudong Han,et al.  A stepwise approach for 3D fracture intersection analysis and application to a hydropower station in Southwest China , 2016 .

[10]  Satish Karra,et al.  Advancing Graph‐Based Algorithms for Predicting Flow and Transport in Fractured Rock , 2018, Water Resources Research.

[11]  Alex Hansen,et al.  Topological impact of constrained fracture growth , 2015, Front. Phys..

[12]  P. Witherspoon,et al.  Porous media equivalents for networks of discontinuous fractures , 1982 .

[13]  M E Newman,et al.  Scientific collaboration networks. I. Network construction and fundamental results. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  A. Barabasi,et al.  Lethality and centrality in protein networks , 2001, Nature.

[15]  R. Freeze A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media , 1975 .

[16]  Caroline Darcel,et al.  Connectivity, permeability, and channeling in randomly distributed and kinematically defined discrete fracture network models , 2016 .

[17]  Satish Karra,et al.  Model reduction for fractured porous media: a machine learning approach for identifying main flow pathways , 2019, Computational Geosciences.

[18]  S. Geiger,et al.  Partitioning Thresholds in Hybrid Implicit‐Explicit Representations of Naturally Fractured Reservoirs , 2020, Water Resources Research.

[19]  Jeffrey D. Hyman,et al.  Identifying Backbones in Three-Dimensional Discrete Fracture Networks: A Bipartite Graph-Based Approach , 2018, Multiscale Model. Simul..

[20]  Albert-László Barabási,et al.  Statistical mechanics of complex networks , 2001, ArXiv.

[21]  P. Davy,et al.  Percolation parameter and percolation-threshold estimates for three-dimensional random ellipses with widely scattered distributions of eccentricity and size , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[22]  C. Darcel,et al.  Flow in multiscale fractal fracture networks , 2006, Geological Society, London, Special Publications.

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

[24]  A. Barabasi,et al.  Scale-free characteristics of random networks: the topology of the world-wide web , 2000 .

[25]  Satish Karra,et al.  Effect of advective flow in fractures and matrix diffusion on natural gas production , 2015 .

[26]  H. S. Viswanathan,et al.  Understanding hydraulic fracturing: a multi-scale problem , 2016, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[27]  P. Davy,et al.  A model of fracture nucleation, growth and arrest, and consequences for fracture density and scaling , 2013 .

[28]  S Karra,et al.  Modeling flow and transport in fracture networks using graphs. , 2018, Physical review. E.

[29]  Satish Karra,et al.  Efficient Monte Carlo With Graph-Based Subsurface Flow and Transport Models , 2018 .

[30]  M. Dentz,et al.  Linking Structural and Transport Properties in Three‐Dimensional Fracture Networks , 2019, Journal of Geophysical Research: Solid Earth.

[31]  P. Taylor,et al.  Initial analysis of air temperature and related data from the Phoenix MET station and their use in estimating turbulent heat fluxes , 2010 .

[32]  Herbert H. Einstein,et al.  Characterizing rock joint geometry with joint system models , 1988 .

[33]  Vito Adrian Cantu,et al.  Machine learning for graph-based representations of three-dimensional discrete fracture networks , 2017, Computational Geosciences.

[34]  Luca Valentini,et al.  The “small-world” topology of rock fracture networks , 2007 .

[35]  Sigmund Mongstad Hope,et al.  Topology of fracture networks , 2012, Front. Physics.

[36]  Paul A. Witherspoon,et al.  A Model for Steady Fluid Flow in Random Three‐Dimensional Networks of Disc‐Shaped Fractures , 1985 .

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

[38]  N. Odling,et al.  Scaling of fracture systems in geological media , 2001 .

[39]  Emmanuel Ledoux,et al.  Modeling fracture flow with a stochastic discrete fracture network: Calibration and validation: 2. The transport model , 1990 .

[40]  Olivier Bour,et al.  Connectivity of random fault networks following a power law fault length distribution , 1997 .