INVESTIGATIONS ON REPRESENTATIVE ELEMENTARY VOLUME AND DIRECTIONAL PERMEABILITY OF FRACTAL-BASED FRACTURE NETWORKS USING POLYGON SUB-MODELS

Abstract

Since the directional permeability of fractured rock masses is significantly dependent on the geometric properties of fractures, in this work, a numerical study was performed to analyze the relationships between them, in which fracture length follows a fractal distribution. A method to estimate the representative elementary volume (REV) size and directional permeability ([Formula: see text] by extracting regular polygon sub-models with different orientation angles ([Formula: see text] and side lengths ([Formula: see text] from an original discrete fracture network (DFN) model was developed. The results show that the fracture number has a power law relationship with the fracture length and the evolution of the exponent agrees well with that reported in previous studies, which confirms the reliability of the proposed fractal length distribution and stochastically generated DFN models. The [Formula: see text] varies significantly due to the influence of random numbers utilized to generate fracture location, orientation and length when [Formula: see text] is small. When [Formula: see text] exceeds some certain values, [Formula: see text] holds a constant value despite of [Formula: see text], in which the model scale is regarded as the REV size and the corresponding area of DFN model is represented by [Formula: see text] (in 2D). The directional permeability contours for DFN models plotted in the polar coordinate system approximate to circles when the model size is greater than the REV size. The [Formula: see text] decreases with the increment of fractal dimension of fracture length distribution ([Formula: see text]. However, the decreasing rate of [Formula: see text] (79.5%) when [Formula: see text] increases from 1.4 to 1.5 changes more significantly than that (34.8%) when [Formula: see text] increases from 1.5 to 1.6 for regular hexagon sub-models. This indicates that the small non-persistent fractures dominate the preferential flow paths; thereafter, the flow rate distribution becomes more homogeneous when [Formula: see text] exceeds a certain value (i.e. 1.5). A larger [Formula: see text] results in a denser fracture network and a stronger conductivity.