Finite element simulations of interactions between multiple hydraulic fractures in a poroelastic rock

Abstract

A fully coupled three-dimensional finite-element model for hydraulic fracturing in permeable rocks is utilised to investigate the interaction between multiple simultaneous and sequential hydraulic fractures. Fractures are modelled as surface discontinuities within a three-dimensional matrix. This model simultaneously accounts for laminar flow within the fracture, Darcy flow within the rock matrix, poroelastic deformation of the rock, and the propagation of fractures using a linear elastic fracture mechanics framework. The leakoff of fracturing fluid into the surrounding rocks is defined as a function of the pressure gradient at the fracture surface, the fluid viscosity, and the matrix permeability. The coupled equations are solved numerically using the finite element method. Quadratic tetrahedral and triangle elements are used for spatial discretisation of volumes and surfaces, respectively. The model is validated against various analytical solutions for plane-strain and penny-shaped hydraulic fractures. Several cases of simultaneous fracturing of multiple hydraulic fractures are simulated in which the effects of the various parameters (the in situ stresses, the distance between fractures, the permeability of the matrix, the Biot poroelastic coefficient, and the number of the fractures in a group) are investigated. The results show that the stress induced by the opening of the fractures, and the stress induced by the fluid leakoff, each have the effect of locally altering the magnitudes and orientations of the principal stresses, hence altering the propagation direction of the fractures. Opening of a fracture induces excessive compression (also known as the “stress shadow”) that causes adjacent fractures to curve away from each other. This excessive compression competes against the differential in situ stresses, which tend to cause the fracture to grow in the plane normal to the minimum in situ stress. The stress shadow effect is reduced by increasing the distance between fractures, and is increased by increasing the leakoff, which may be due to increased permeability of the rock, or an increase in the Biot coefficient.