Abstract: The pore structure and permeability of surrounding rock in water-rich mine often undergo significant changes under mining disturbances, which can induce the formation of water-conducting fracture zones and even lead to water inrush disasters. The formation of water-conducting fracture zones is not only highly unpredictable, but also involves the coupling effect of stress and seepage in the fractured surrounding rock during the mining process, which makes the model development challenging. To study the formation law of the fracture zones under the coupling effect of stress and seepage during the mining process, the surrounding rock containing a large number of pores and fractures is simplified as the dual-porosity medium, which is formed by the overlapping of porous matrix and fractures. Then, based on the Mixture Coupling Theory, a thermodynamically consistent constitutive model, namely the HMD (hydro-mechanical-damage) constitutive model, is established with non-equilibrium thermodynamic methods that considers the coupling of stress, seepage and damage. Meanwhile, the dynamic evolution equations of stress, porosity and fracture volume fraction under HMD coupling conditions are derived. Finally, the derived coupling equations are applied to practical engineering scenarios, and the evolution laws of pore water pressure, stress, vertical displacement and damage of the surrounding rock are numerically simulated and analyzed with finite element method software. The possible locations of water-conducting fracture zones during the mining process of the water-rich mine are also identified.