At present, in the vast majority of coal mine production processes in China, the degree of hydrogeological exploration often lags behind geological exploration. The main difficulty in evaluating the water richness of coal seam top and bottom water-bearing beds is that the existing evaluation methods often rely on less hydrogeological investigation data. How to utilize the abundant geological exploration data in the mining area to appraise the water-rich distribution of sandstone aquifers is a feasible and challenging methodology. At present, some experts and scholars have tried to use multivariate factor analysis to solve the problem of water-richness evaluation, and they have achieved certain results, but there are some shortcomings: (1) The prediction results are mostly qualitative estimations of the water-richness grade, and there is a lack of quantitative analysis of the units-inflow; and (2) at present, the more advanced prediction methods, such as the back propagation (BP) neural network model, have the disadvantages of low accuracy, requiring many iterations, and slow convergence speed. Therefore, with geological exploration data of the 1503E working face of the Yili No.1 coal mine as the basis., this paper uses grey correlation analysis to screen out the factor indexes suitable for the evaluation of the water richness of a weakly cemented sandstone aquifer, and it combines principal component analysis (PCA) with a BP neural network. Based on the selected factor indexes, a prediction model of the water richness of a weakly cemented sandstone aquifer is established. The results show that compared with the existing methods, the prediction accuracy is higher and has a certain application value.