Abstract: In recent years, the land in coal mining subsidence areas has been widely used as construction land.The stability of the building foundation directly determines the feasibility of new buildings in the subsidence area, and the selection of foundation stability analysis methods is an important part of the foundation stability analysis of the building foundation. In order to screen out the optimal method for the stability analysis of the building foundation in the coal mining subsidence area, this paper took the Huaibei Zhonghu coal mining subsidence area as the research object, and discussed the influence of the foundation stability of the newly built buildings in the subsidence area. Based on the main factors of stability, the evaluation criteria of foundation stability were summarized, and probability integral method, mechanical analysis method and empirical formula method were used to evaluate the foundation stability of the study area, and the workload, accuracy and applicable conditions in the comprehensive analysis process were comprehensively analyzed. The methods were compared and analyzed with the model construction factors, and the most suitable method for the stability analysis of building foundation in coal mining subsidence area was pointed out. The research shows that the empirical formula method has small workload, high accuracy, strong applicability, and high model adjust ability in the actual foundation stability analysis. It is more suitable for practical engineering applications of building foundation stability analysis in coal mining subsidence areas. The analytical method was followed by the probability integration method; the evaluation results of the foundation stability of the research area by the probability integration method, the mechanical analytical method and empirical formula method were basically consistent with the engineering reality, and the foundation of the research area was in a stable state; the probability integration method realizes the future of the study area.The ten-year residual surface deformation is expected to be in a slow-moving state in the future, with a maximum residual subsidence of 441 mm, a maximum residual slope of 1.6 mm/m, a maximum residual curvature of 0.02 mm/m2, and a maximum residual horizontal deformation of 1.5 mm/m. The maximum residual horizontal movement is 110 mm.