<p>A THMC coupled phase field framework for modeling fractures in porous rocks is proposed in this study. The framework introduces additionally the damage variable governed synergistically by the phase field and chemical field to account for dual-fracture mechanisms. Through this damage variable, full coupling of the temperature, hydraulic, mechanical, chemical, and phase fields is achieved. Implemented on the COMSOL Multiphysics platform, this multi-field coupling framework is solved by using a staggered iteration algorithm. The proposed framework was verified through fracture propagation induced by various factors. Furthermore, two-dimensional case studies are conducted to investigate the effects of acid concentration, heterogeneity, injection rate, specific surface area, and scale parameters on fracture morphology, fluid pressure distribution, temperature distribution, pressure evolution, and fracture propagation range. Numerical results demonstrate that the predictions of the proposed THMC coupled phase field model for fracture evolution and acid breakthrough consumption align with existing studies, while effectively characterizing the influence of sensitivity parameters.</p>