：Long term underground coal mining results in land subsidence and deformation above the goaf area, causing damage to agriculture and degradation of soil, it is crucial to understand the geographical variability of soil nutrients in areas affected by coal mining for effective grain production and targeted reclamation activities. This study utilizes a multi-fractal methodology to clarify the influence of coal mining subsidence on the spatial diversity of soil nutrients in coal and grain composite areas with high groundwater levels. An inquiry was conducted in five selected regions in Zhao Gu Yi Mine, located in Henan Province. In addition to the control check (CK), the selected regions located in the coal mining subsidence area included the perennial subsidence ponding (PSP), seasonal subsidence ponding (SSP), neutral zone (NZ) and horizontal deformation zone (HDZ). Parameters such as generalized dimension parameters( and ) and singularity spectrum parameters( and )were utilized to quantitatively analyze the soil nutrients. The results showed that the heterogeneity of soil organic carbon (SOC) distribution in the horizontal deformation zone increased significantly compared with other regions, and the spatial variability of SOC was dominated in the high value region of the topsoil layer. Total phosphorus (TP) showed great heterogeneity only in the surface layer of seasonal subsidence ponding and horizontal deformation zone, and TP was dominated by low value and high value respectively. Water transport made the two water-soluble nutrients, available nitrogen (AN) and available phosphorus (AP), showed strong heterogeneity in the seasonal subsidence ponding, and showed a more sensitive state in the low value area. Compared with the uncertainty of field management, coal mining subsidence did noy show significant influence on the spatial distribution of soil total nitrogen (TN). Compared with classic geostatistical methods, multi-fractal methods are better suited for evaluating the spatial variability of soil nutrients in coal mining subsidence sites. By acquiring accurate data on the variability of soil nutritional characteristics at various levels, it offers scientific, theoretical, and technical assistance for targeted and precise land reclamation plans.