Mineral-Scale Mechanical Properties of Carbonate Rocks Based on Nanoindentation

Carbonate reservoirs in the Shunbei area develop pronounced fracture networks after acidized hydraulic fracturing and thus have the potential to be repurposed as underground gas storage (UGS) after hydrocarbon depletion. Characterizing their mechanical behavior is essential for safe UGS operation; however, deep to ultra-deep natural cores are difficult to obtain, and conventional macroscopic tests often cannot provide parameters that meet engineering requirements. To address this issue, nanoindentation combined with QEMSCAN (Quantitative Evaluation of Minerals by Scanning Electron Microscopy) was employed to quantify microscale mineral distributions and the mechanical properties of the major constituents. The investigated rock is calcite-dominated (89.62%), with minor quartz (9.89%) and trace feldspar-group minerals (1.89%). Minerals are randomly embedded, and soft–hard phase boundaries are widely distributed. A finite–discrete element method (FDEM) model was then constructed and calibrated in ABAQUS. The discrepancies in uniaxial compressive strength and elastic modulus relative to laboratory results were 6.51% and 9.91%, respectively, indicating good agreement in both mechanical response and failure mode. Parametric analyses using three additional models with different mineral proportions show that damage preferentially initiates at mineral phase boundaries and stress concentration zones induced by end constraints. Microcracks then propagate and coalesce into a dominant compressive–shear band, and final failure is mainly governed by slip along the shear band with localized tensile cracking. With increasing quartz and feldspar contents, enhanced heterogeneity and a higher density of phase boundaries lead to a higher density of crack nucleation sites and increased crack branching, and the failure pattern transitions from a single shear-band–controlled mode to a more network-like fracture system. Moreover, macroscopic strength is not determined solely by the intrinsic strength of individual minerals; heterogeneity and phase-boundary characteristics strongly govern microcrack behavior, such that higher hard-phase contents may result in a lower peak strength.