The Jurassic aquifer in Northwest China is the key aquifer for mine water filling, which is significant due to its loose structure, large porosity, strong rock permeability, and fracture development characteristics. In addition, the water richness in space is extremely uneven, and many coal mine roof water inrush events are closely related to it. A case of evaluation of water-richness of the roof sandstone in the 3-1 coal seam of the Bayangaole minefield was analyzed in depth, and the evaluation index system is established based on lithology and structural characteristics. Specifically, the evaluation indexes are under the influence of the influencing factors of lithology, the density of fault intersection endpoints, and the density of fault scale and the strength of folds as the influencing factors of structure. On this basis, the set pair analysis-variable fuzzy set coupling evaluation method is introduced to form a targeted water-rich evaluation model of a roof sandstone aquifer. By using the coupling method of set pair analysis and variable fuzzy set, a targeted evaluation model is formed to realize the organic integration of indicators. Through the comprehensive analysis of the relative zoning of water abundance and the data from the borehole pumping (drainage) test, the distribution of water abundance grade in the study area is clarified.

The flow of cement‐based slurry within concentric annular geometry is a major problem of study in the field of engineering, especially regarding the prevention of water ingress in an annulus formed during shaft construction utilizing the artificial freezing technique in China. In this study, an analytical governing equation of motion for the axial flow of an incompressible Newtonian fluid in a long concentric annulus is derived under the condition of high groundwater pressure. A stepwise calculation method is proposed to describe the grouting process based on two injection modes, namely, flow rate control and pressure control. The injection time is divided into a series of time segments; correspondingly, the grouted zone is subdivided into infinitesimal elements. Some key parameters, such as the location, dimension, slurry viscosity, and pressure gradient of each element, can be obtained using the developed MATLAB program. On this basis, the distribution of pressure and slurry viscosity in the grouted zone and the variations in injection pressure at the grouting point and grouting flow rate are determined. Two injection mode cases are investigated to reveal the grout propagation in the concentric annulus. Finally, numerical simulations are conducted and employed to validate and calibrate the calculated results. The results obtained by the present stepwise calculation method show good agreement with the numerical results.