Grain-scale heterogeneity strongly influences hydraulic fracture initiation and trajectory in crystalline rocks, yet its contributions to non-planar growth and the interaction of multiple nearby cracks remain insufficiently quantified. To address this gap, we perform numerical experiments on a model containing two parallel pre-existing cracks using a hydro-mechanical phase-field framework, systematically quantifying how mineral distribution and axial compression govern non-planar hydraulic fracture growth and inter-fracture competition. The results demonstrate that mineral distribution is the primary driver of fracture complexity. Even within the same Voronoi tessellation, redistributing minerals alone yields markedly different trajectories, deflections, branching patterns, and final morphologies. Furthermore, non-planar growth follows a stepwise, energy-threshold-driven mechanism. When cracks penetrate strong grains or undergo large-angle deflections, propagation is impeded, and injection pressure builds up. Once a critical energy threshold is reached, accumulated energy is rapidly released along the path of minimum incremental energy, manifested as abrupt pressure drops and rapid crack advance. Additionally, the two nearby fractures exhibit strong mechanical competition. Despite negligible hydraulic interference in low-permeability granite, early growth of one fracture redistributes stresses and suppresses the driving force of the other, resulting in asymmetric development. Finally, axial compression primarily governs the overall propagation orientation and influences local failure modes but has a limited effect on peak pressure relative to mineral distribution.
Zhang et al. (Mon,) studied this question.