The appropriate understanding and formulation of rock discontinuities via FEM is still challenging for rock engineering, as continuous algorithms cannot handle the discontinuities in rock mass. Also, different failure modes of rock samples, containing tensile and shear failure, need to be computed separately. In this study, a novel double-phase field damage model was introduced with two independent phase field damage variables. The construction of the proposed model follows the thermodynamics framework from the overall Helmholtz free energy, with elastic, plastic and surface damage components. The proposed model is calibrated via traditional damage variables, based on ultrasonic wave velocity measurement and acoustic emission monitoring, and both show great consistency between simulation results and laboratory observations. Then the double-phase field damage model is applied to COMSOL software to simulate microcrack propagation in a pre-fractured rock sample. Both lateral and wing cracks are observed in this study, manifested as shear- and tensile-dominated cracks. We also observed different microcracking mechanisms in the proposed numerical models, such as tensile and shear cracking, the influence of plastic strain and the percolation between tensile and shear microcracks. Overall, this study provides valuable insights into the mechanics of microcracking in rocks, and the proposed model shows promising results in simulating crack propagation.