Modelling the variation in the behaviour of pre-fractured rocks subjected to hydraulic fracturing with permeability of the rock matrix using finite-discrete element method

Abstract

In this study, we present the results of two-dimensional numerical simulations for the effects of rock matrix permeability on the behaviour of hydraulic fractures in intact and pre-fractured rocks. The simulations are performed using the Finite-Discrete Element Method (FDEM). In this method, the deformation and fluid pressure fields within the porous rock blocks, pre-existing fracture network, and hydraulically induced fractures are calculated through a fully coupled hydromechanical scheme. Furthermore, new fractures can initiate in crack elements located between each pair of finite elements and can propagate in any path that the boundary and loading conditions require according to non-linear fracture mechanics criteria. Fluid channels are also defined between pairs of finite elements simulating the inter-connected flow paths through porous media. Four models of the rock mass are created in this study: (i) homogeneous-impermeable, (ii) homogeneous-permeable, (iii) heterogeneous-impermeable matrix, and (iv) heterogeneous-permeable matrix. Heterogeneous rock masses contain a discrete fracture network (natural fractures) in the rock mass structure. Hydraulic fracturing is modelled in domains of 40×40 m2 with the four different structures and mass transport capacities, and the results are compared to each other. The results highlight the significant effect of diffusive fluid flow through rock blocks, in addition to the flow through fracture network, on the global hydromechanical behaviour of the rock mass. These results help to understand the governing hydromechanical processes taking place in fractured rocks with matrix of different permeability, such as granites, shales, carbonate rocks, and sandstones and the extent of complexities required to model their behaviour to achieve reasonable accuracy.