To address the synergistic stability evaluation of coal pillars and caving deposits in room-and-pillar mining of nearly vertical coal seams, this study takes the 101 Coal Mine (104th Regiment, Xishan Area, Urumqi, Xinjiang) as the engineering background. It combines physical similarity simulation and theoretical analysis to explore the synergistic bearing mechanism of coal pillars and caved deposits. Based on limit equilibrium theory, a combined instability criterion considering roof mudstone’s bending-toppling and shear-sliding is established; the Rankine earth pressure theory is modified, and a stability coefficient K_s (reflecting synergistic bearing effect) is proposed to realize quantitative evaluation of goaf stability. A model experiment simulates the mining of a 73° nearly vertical coal seam. Results show the roof instability mode (under coal pillars and caved deposits) is equivalent to anti-dip slope’s toppling-sliding composite failure. Experimental and theoretical results agree well, verifying the model’s rationality and applicability. This study provides a theoretical basis and analytical method for calculating the synergistic stability of coal pillars and caving deposits in nearly vertical coal seams.

Particle Flow Code (PFC) numerical simulations were adopted to simulate the mining process and the process of goaf collapse and to predict the residual subsidence of abandoned goafs in steeply inclined multi-seam coal mines, taking the No. 101 Coal Mine in the Xishan Mining Area of Urumqi, China, as an example. Scaled physical simulations were also employed to simulate the evolution of voids in the coal–rock mixture in the goaf. The results show that after mining, the roof of shallow coal seams becomes unstable and collapses in the anti-dip direction, causing the materials within the unconsolidated layer to fall and backfill the goaf, which further leads to ground subsidence. The mining of deep coal seams is also accompanied by the overall movement of overlying strata along the dip direction of the coal seams and surface subsidence. The content of voids within the broken coal–rock mass in the goaf tends to decrease with increasing pressure, showing a negative exponential correlation. Based on the observed relationship between displacement and void content obtained from the simulation experiments, it is inferred that the residual displacement under the current conditions of the study area accounts for approximately 10.5% of the total displacement. Combining the results of the PFC simulation and the evolution law of void content, the residual subsidence of the goaf in the study area since mine closure is predicted to range from 0 to 1 m, with a high-value zone distributed in the northeastern part of the study area. Deep goafs within the B₇–B_11–12 and B₁₄–B₁₈ coal seam groups mainly contribute to the residual subsidence. The distribution of goaf collapse pits, as revealed by field investigation, also verifies the reliability of the prediction results.