Abstract：By 8km × 8km grid, the land was divided into 136 sampling points in Liaocheng city. The contents of soil organic matter, total N, total P and total K were monitored. Geostatistics method and Moran’s I of spatial autocorrelation method were used to analyze the spatial features of soil nutrients. Study shows that, the spatial heterogeneity of soil nutrients’ content changes mainly by random factors, for example, human activities such as irrigation, fertilization, cultivation measures, cropping systems, etc. The spatial correlation of the soil nutrients’ content is very weak. On the whole, there is an obvious spatial distribution of the four soil nutrients which gradually declines from the eastern and southern to the western and northern part of Liaocheng city. In the scope of the 40km, Moran's I coefficients of four kinds of soil nutrients are positive which shows a positive spatial correlation. The response of organic matter on the scale is more sensitive than that of other three nutrients. The change of Moran's I of total N and total P is similar, but Moran's I coefficient of total K outside the scope of 64km performs the volatility in the vicinity of 0, which may be related to the local agricultural activities.