Propagation of Hydraulic Fractures and Natural Fractures: The Bypassing Behavior in 3D Space

Summary The interaction between natural fractures (NFs) and hydraulic fractures (HFs) in 3D space poses significant challenges to numerical simulation. The interaction behaviors between HFs and NFs in 3D space include crossing, offset, stopping, and bypassing. Many existing numerical simulators are 2D, which inherently limits their ability to account for the vertical growth of fractures. Consequently, they are unable to model the bypassing behavior effectively. In this paper, a mathematical model for the propagation of HFs and NFs in a naturally fractured reservoir is established. The displacement discontinuity method (DDM) is used to solve the rock deformation, while the finite difference method (FDM) is utilized to solve the fluid flow within fractures. The undirected graph structure is used to represent the complex fracture network, and a dynamic adjustment of grid connectivity method is implemented to describe the process by which HFs bypass NFs. By integrating with graphics processing unit (GPU) computing, an efficient 3D simulator for HFs and NFs propagation is developed. The accuracy is verified against analytical solutions and reference solutions. Subsequently, a series of numerical studies on the bypassing behavior are conducted. The primary findings are as follows: (1) The simulator can accurately capture the interactions between HFs and NFs in 3D space, including crossing, arresting, and bypassing behaviors. (2) The bypassing behavior is characterized by three distinct processes—intersection, height growth of HF, and bypassing. (3) During the intersection stage, both the injection pressure and maximum NF width increase. During the height growth stage, the pressure is relatively stable, while the maximum NF width continues to increase. These two stages together result in the “storage” phenomenon. Once the bypassing behavior occurs, both the pressure and the maximum NF width decrease sharply, leading to a “release” phenomenon. Additionally, the “storage” and “release” phenomena may impact proppant transport. (4) Given the presence of bypassing behavior, it is essential to consider the NFs in 3D fracture simulators. The simulator and its findings can provide valuable insights for field design.