Abstract: Based on a grid (25 m × 25 m) equidistant sampling, the spatial variability of pH, organic matter, total nitrogen, available phosphorus, CEC and three typical heavy metal elements Cd, As and Pb in soil tillage layer (0-20 cm) were analyzed by using GIS and Geostatistics in the paddy field of 3.56 hm² in Beishan Town, Changsha County, Hunan Province. The results showed that soil pH value and Pb content showed weak variation, and other indexes showed moderate variation. The order of variation was following available phosphorus ＞ Cd ＞ total nitrogen ＞ organic matter ＞ CEC ＞ As ＞ Pb ＞ pH. Results of the semi-variance test showed that the best fitting model of the semi-variance function of organic matter, available phosphorus and As was exponential, and the best semi-variance function of pH, total nitrogen, CEC, Cd, Pb was spherical. All the indicators had a strong spatial correlation except for CEC, which showed moderate spatial correlation. Kriging interpolation analysis showed that pH, total nitrogen, CEC, Pb were plaque distribution, while organic matter, available phosphorus, Cd and As were block and banded distribution. Vegetation, topography and human activities were the main factors driving the variation of soil nutrients and heavy metals in the study area. The correlation between soil nutrients and heavy metals content was significant, among which pH and organic matter, Cd and Pb reached a very significant correlation level.