Under the influence of external factors such as high temperature, chemical corrosion, and loading-unloading cycles, micro-cracks and pores may develop within rock structures. Continuum damage mechanics (CDM) characterizes these micro-defects arising within rocks, investigating the physico-mechanical behaviors of rocks in damaged states. This study, grounded in the theoretical framework of CDM and utilizing the Mohr-Coulomb failure criterion, integrates different probability distributions and damage coupling methods to develop three statistical damage models. The investigation primarily explores the effects of strength criteria, probability functions, and damage coupling mechanisms on the statistical damage models. In addressing the complex interplay between rock micro-defects and external stressors such as high temperatures and chemical corrosion, this study leverages the principles of CDM within the Mohr-Coulomb criterion framework. It advances the field by developing three innovative statistical damage models, enriched by diverse probability distributions and damage coupling methods. The research delineates the substantial impact of strength criteria, probability functions, and damage coupling mechanisms on the behavior of these models. Notably, it achieves a synthesis of theoretical constructs and experimental validation, demonstrating the models' capacity to reflect the nuanced behaviors of damaged rock accurately. Furthermore, this investigation contributes a methodological blueprint for statistical damage model construction, underscoring the critical selection of strength criteria and probability distributions. These efforts fortify the theoretical underpinnings necessary for the enhanced engineering application of statistics-driven damage analysis.