[Objective] This study was to explore the difference of kriging interpolation and sequential Gaussian simulation on analyzing soil heavy metal pollution with a view to provide references for analyzing the heavy metal pollution of soil. [Method] The sampling data of soil copper from a county of Liaocheng, Shandong Province was set as the study objective. Kriging interpolation and sequential Gaussian simu- lation were used to simulate the spatial distribution of soil copper. And 30 sampling points were selected as the cross-validation data set to compare the two interpola- tion methods. [Result] Kriging method and Gaussian sequential simulation have their own advantages on simulating mean segment and extreme segment, therefore, re- searchers should choose the proper method based on the characteristics of test data and application purposes. [Conclusion] Analysis of soil heavy metal pollution is the prerequisite for soil management and ecological restoration. The result of this study is of important significance for choosing different interpolating and simulating methods to analyze soil heavy metal pollution based on different purposes.