Numerical Investigation on Mesoscale Evolution of Hydraulic Fractures in Hydrate-Bearing Sediments

Hydraulic fracturing is widely recognized as a potential stimulation technology for the development of challenging natural gas hydrate. However, the fracturing behavior of non-diagenetic hydrate reservoirs has peculiar characteristics that are different from those of conventional oil and gas reservoirs. Herein, a fully coupled fluid-mechanical model for simulating hydraulic fracturing in hydrate-bearing sediments (HBS) was established based on the discrete element method, and the influence of hydrate saturation, in situ stress, and injection rate on the meso-fracture evolution was investigated. The results indicate that with the increase in hydrate saturation, the fracture morphology transitions from bi-wing to multi-branch, thereby enhancing fracture complexity. Both tensile and shear failure modes exist, and the tensile failure between the weakly cemented sediment particles is dominant. The tensile strength of HBS is an exponential function of hydrate saturation, with the breakdown pressure being governed by hydrate saturation and in situ stress, with the form being consistent with the classical Kirsch equation. Additionally, lower in situ stress and higher injection rates are conducive to the generation of microcracks, whereas an excessive injection rate reduces the fracture length. These findings contribute to understanding the meso-evolution mechanism of hydraulic fractures and guide the design of on-site hydraulic fracturing plans of natural gas hydrate reservoirs.