Determining the equivalent permeability tensor for fractured rock masses using a stochastic REV approach: Method and application to the field data from Sellafield, UK

A numerical procedure to determine the equivalent permeability tensor of a fractured rock is presented, using a stochastic REV (Representative Elementary Volume) concept that uses multiple realizations of stochastic DFN (Discrete Fracture Network) models. Ten square DFN models are generated using the Monte Carlo simulations of the fracture system based on the data obtained from a site characterization program at Sellafield, Cumbria, UK. Smaller models with varying sizes of from 0.25 mx0.25 m to 10 mx10 m are extracted from the generated DFN models and are used as two-dimensional geometrical models for calculation of equivalent permeability tensor. The DFN models are also rotated in 30degrees intervals to evaluate the tensor characteristics of calculated directional permeability. Results show that the variance of the calculated permeability values decreases significantly as the side lengths of the DFN models increase, which justifies the existence of a REV. The REV side length found in this analysis is about 5 m and 8 m with 20% and 10% acceptable variations, respectively. The calculated directional permeability values at the REV size have tensor characteristic that is confirmed by a close approximation of an ellipse in a polar plot of the reciprocal of square roots of the directional permeability.