Fluid permeability of deformable fracture networks

The authors consider the problem of defining the fracture permeability tensor for each grid lock in a rock mass from maps of natural fractures. For this purpose they implement a statistical model of cracked rock due to M. Oda [1985], where the permeability tensor is related to the crack geometry via a volume average of the contribution from each crack in the population. In this model tectonic stress is implicitly coupled to fluid flow through an assumed relationship between crack aperture and normal stress across the crack. The authors have included the following enhancements to the basic model: (1) a realistic model of crack closure under stress has been added along with the provision to apply tectonic stresses to the fracture system in any orientation, the application of stress results in fracture closure and consequently a reduction in permeability; (2) the fracture permeability can be superimposed onto an arbitrary anisotropic matrix permeability; (3) the fracture surfaces are allowed to slide under the application of shear stress, causing fractures to dilate and result in a permeability increase. Through an example, the authors demonstrate that significant changes in permeability magnitudes and orientations are possible when tectonic stress is applied to a fracture system