A Model for Steady Fluid Flow in Random Three‐Dimensional Networks of Disc‐Shaped Fractures

A model for steady fluid flow in three-dimensional, random networks of fractures has been developed. In this model the fractures are disc shaped discontinuities in an impermeable matrix. The fracture discs can be arbitrarily located within the rock volume and can have any desired distribution of aperture, radius orientation, and density. Thus where the disc model is appropriate it is possible to calculate flow through fracture networks which are statistically similar to those that occur in nature. After the boundary conditions and the desired fractures are specified, the intersections (nodes) between these discs (elements) are identified. Then steady flow through the network is calculated using a mixed analytical-numerical technique. In each fracture, analytic equations for flow into or out of each node as a function of the average head at each node are developed. The equations are based on image theory and the assumption that each node is a source (or sink) of uniform strength. A set of mass balance equations is constructed which equate flow into a node from one of its associated fractures to flow out of the node into the other associated fracture. These equations are solved for the average head at each node, and flux between fractures can then be calculated by substituting the average head values back into the analytical equations. The model has been successfully checked against analytical results for several cases of two and three intersecting fractures. We plan to use these techniques to measure the permeability of fracture networks.