Heterogeneity at the grain scale strongly influences hydraulic fracturing in crystalline rock; however, systematic studies quantifying its impacts on the evolution of injection pressure and crack propagation remain limited. To address this gap, we employ a hydro-mechanical phase-field model incorporating Voronoi-based microstructures to systematically quantify the effects of grain-scale heterogeneity on hydraulic fracturing. Two numerical experimental programs are designed to examine the effects of (i) mean grain size and (ii) mineral distribution under different axial stresses. The simulations reveal a close coupling between injection pressure and crack-length evolution, and both responses are strongly governed by grain-scale heterogeneity. When the fracture enters weak minerals, it advances rapidly and pressure drops; when it encounters on strong minerals, growth slows or arrests and pressure builds until a threshold triggers the next advance. Moreover, peak pressure statistics further indicate that mineral distribution dominates the response scatter, while axial stress plays a secondary role. Specifically, the mean peak pressures at 0 and 10 MPa are similar (about 14.31 and 14.21 MPa), whereas rearranging minerals within the same Voronoi tessellation changes peak pressure by more than 4 MPa. Higher peaks occur when strong minerals lie ahead of the initial crack tip, increasing resistance to initiation and early growth. Finally, the stress state modulates fracture trajectories: under low axial stress, fractures preferentially follow mineral boundaries, whereas higher axial stress strengthens macroscopic stress guidance and shifts the path toward a direction closer to being perpendicular to the maximum principal stress. This trend is consistent with energy minimization, since interface detouring under high axial stress incurs a larger elastic free energy penalty.

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.