In order to solve the related problems of roof collapse in deep thick loose seam, taking the actual geological conditions of 1302 working face of Zhaolou Coal Mine as the geological background. A mechanical calculation and analysis model for a deep thick loose seam is established. A mining dynamic simulation was carried out using the numerical simulation software FLAC3D. Fracture development height was observed using a multi-stage plugging and leak detection system. Sublevel solid compaction filling technology was used to fill the goaf. Research shows that the caving zone, fissure development zone, and bending subsidence zone are formed during thick loose seam mining. The caving zone area is larger when the ground stress on the roof is larger, and fracture zones are better developed because the density of the loose seam area is smaller. Numerical simulation results are slightly larger than the measured results, so the roof overburden activity can be effectively predicted using the numerical simulation. The problem of serious roof collapse and surface subsidence in goafs can be effectively controlled using sectional compact filling technology.

A coupled thermal–hydraulic–chemical (THC) model was carried out in this paper to study the influence of rock heterogeneity and the coupling effect of temperature, groundwater, and hydrochemistry on rock damage. Firstly, the hydrochemical and hydraulic erosion equations were established. The equations of the coupled THC model were established by combining the hydrochemical and hydraulic erosion equations, the flow equations, and the heat transfer equations. Weibull distribution was adopted to govern the heterogeneity of initial rock porosity distribution. Secondly, the influence of the hydrochemistry, the temperature and the initial porosity heterogeneity on porosity and fluid velocity change was studied. Then the rock damage rule changed with time at different pH values and temperature was studied. Finally, an actual deep coal mine model was established to apply the THC model to predict water inrush. Results indicate that: (1) The average porosity and average fluid velocity approximately show linear growth and exponential growth with time, respectively, and their growth rates increase with decreasing pH value and increasing temperature in a certain acidity and temperature range. (2) The increase of initial porosity heterogeneity has little influence on porosity change, but it can increase the fluid velocity growth rate. The porosity heterogeneity and fluid velocity heterogeneity approximately show exponential growth with increasing time, and the rock heterogeneity growth contributes to form cracks. The increase of temperature and decrease of pH value have little influence on the porosity heterogeneity, but they can increase the growth rate of the fluid velocity heterogeneity. (3) The rock damage shows linear growth with time, and its growth rate increases with decreasing pH value and increasing temperature in a certain acidity range and temperature range. (4) The increase of rock heterogeneity can increase the possibility of water inrush.