Numerical Investigation of Hydraulic Fracture Propagation Based on Cohesive Zone Model in Naturally Fractured Formations

Complex propagation patterns of hydraulic fractures often play important roles in naturally fractured formations due to complex mechanisms. Therefore, understanding propagation patterns and the geometry of fractures is essential for hydraulic fracturing design. In this work, a seepage-stress-damage coupled model based on the finite pore pressure cohesive zone (PPCZ) method was developed to investigate hydraulic fracture propagation behavior in a naturally fractured reservoir. Compared with the traditional finite element method, the coupled model with global insertion cohesive elements realizes arbitrary propagation of fluid-driven fractures. Numerical simulations of multiple-cluster hydraulic fracturing were carried out to investigate the sensitivities of a multitude of parameters. The results reveal that stress interference from multiple-clusters is responsible for serious suppression and diversion of the fracture network. A lower stress difference benefits the fracture network and helps open natural fractures. By comparing the mechanism of fluid injection, the maximal fracture network can be achieved with various injection rates and viscosities at different fracturing stages. Cluster parameters, including the number of clusters and their spacing, were optimal, satisfying the requirement of creating a large fracture network. These results offer new insights into the propagation pattern of fluid driven fractures and should act as a guide for multiple-cluster hydraulic fracturing, which can help increase the hydraulic fracture volume in naturally fractured reservoirs.