Abstract
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.