Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Abstract

AbstractFlow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low‐permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three‐dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady‐state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance‐ranking/hierarchy of different flow paths in the network using graph‐based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network‐scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high‐density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.