Numerical simulation of crack propagation in layered formations

In the present study, a boundary element method based on the higher order displacement discontinuity formulation is presented to solve the general problem of hydraulic fracture propagation in layered formations. Displacement collocation technique is employed to model the higher order displacement variation along the crack and the special crack tip element near its ends. The hydraulic fracture propagation and its interaction with the layer interface in non-homogenous rock materials are studied by the proposed semi-analytical (hybridized boundary element-boundary collocation) method. The maximum tangential stress criterion (or sigma-criterion) of fracture mechanics considering different elastic constants (Young modulus and Poisson's ratio) is used to obtain the fracture path. The fracture propagation from stiff to soft and soft to stiff media for cracks having different inclination angles is modeled, and the effects of elastic constants on the hydraulic fracture propagation is studied. The results show that if the hydraulic fracture originates in the stiffer layer, its capability to cross the layer increases and is vice versa for the softer material. The comparison of the results gained from the numerical method with those in the literature show a good performance of the method in the case of propagation of hydraulic fracture in layered formations.