Hydraulic Fracture Propagation in Layered Heterogeneous Rocks with Spatially Non-Gaussian Random Hydromechanical Features

Abstract

This study proposes a stochastic method to analyse the propagation of hydraulic fractures affected by layered heterogeneity in rocks in a toughness-dominated regime. The study utilises the phase-field method in the context of two-dimensional finite element analysis to model the hydraulic fracture (HF) propagation in rock materials in laboratory scale. Field data on hydrogeologic properties of some rocks reveal that material heterogeneity may appear in the form of leptokurtic marginal distributions. Generalised sub-Gaussian (GSG) model is capable of capturing physical characteristics of such rocks, and it is employed to stochastically model rocks with layered lithologic heterogeneity by generating a large number of auto- and cross-correlated random fields for hydro-geomechanical properties. To investigate the sensitivity of the cracking response to the inherent characteristics of material heterogeneity, various GSG distribution forms are considered in Monte Carlo (MC) analyses. The HF’s deviation from the theoretically predicted direction, which is perpendicular to the direction of the minimum in situ stress, is correlated with the distribution of hydro-geomechanical properties, showing a Gaussian-type distribution. This study concludes that the differential stress and the bedding orientation are the main factors affecting the HF deviation and the required breakdown pressure for initiating the HF propagation from a borehole. In the application of directional hydraulic fracturing (DHF), the effect of bedding layers becomes dominant when the bedding orientation is aligned with the direction of perforations in the boreholes.