Modeling flow and transport in fracture networks using graphs

TL;DR

This paper introduces a graph-based approach to model flow and transport in fracture networks, significantly reducing computational costs while maintaining high accuracy after bias correction, enabling efficient uncertainty quantification.

Contribution

The authors develop and validate a bias-corrected graph algorithm for simulating flow and transport in fracture networks, offering a scalable alternative to high-fidelity DFN models.

Findings

Graph approach is up to 10,000 times faster than DFN models.

Bias correction improves the accuracy of the graph method to closely match DFN results.

Graph method with bias correction effectively supports uncertainty quantification in subsurface flow modeling.

Abstract

Fractures form the main pathways for flow in the subsurface within low-permeability rock. For this reason, accurately predicting flow and transport in fractured systems is vital for improving the performance of subsurface applications. Fracture sizes in these systems can range from millimeters to kilometers. Although, modeling flow and transport using the discrete fracture network (DFN) approach is known to be more accurate due to incorporation of the detailed fracture network structure over continuum-based methods, capturing the flow and transport in such a wide range of scales is still computationally intractable. Furthermore, if one has to quantify uncertainty, hundreds of realizations of these DFN models have to be run. To reduce the computational burden, we solve flow and transport on a graph representation of a DFN. We study the accuracy of the graph approach by comparing breakthrough times and tracer particle statistical data between the graph-based and the high-fidelity DFN approaches, for fracture networks with varying number of fractures and degree of heterogeneity. We show that the graph approach shows a consistent bias with up to an order of magnitude slower breakthrough when compared to the DFN approach. We show that this is due to graph algorithm's under-prediction of the pressure gradients across intersections on a given fracture, leading to slower tracer particle speeds between intersections and longer travel times. We present a bias correction methodology to the graph algorithm that reduces the discrepancy between the DFN and graph predictions. We show that with this bias correction, the graph algorithm predictions significantly improve and the results are very accurate. The good accuracy and the low computational cost, with $O(10^4)$ times lower times than the DFN, makes the graph algorithm, an ideal technique to incorporate in uncertainty quantification methods.