Abstract: In this paper，the Wucaiwan，desert zone and Jiangjun temple in Zhundong region was used as the research target and used multivariate statistical analysis and spatial analysis method for correlation analysis and spatial distribution of soil heavy metals between Wucaiwan，desert zone and Jiangjun temple in Zhundong.The results showed that the amplitude of Hg and As were low in these three region，soil heavy metal content was more concentrated，on the contrast，the amplitude of Cu，Cr，Pb and Zn were high.The three regions were extracted four principal components by PCA；Hg，As，Cu，Cr，Pb and Zn elements were the main factors that affected the contents of heavy metals in the soil area of the colorful Bay and the desert.As，Cu，Cr and Pb elements were the main factors that affected contents of heavy metals in Jiang Jun Temple.By analyzing the correlation of soil heavy metals in different regions，Zn elements in Wucaiwan were negatively correlated with Zn elements in Jiangjun temple，and the correlation coefficient was -0.753，Cr elements in the desert zone and Hg，Cr elements were significantly positive correlation in Jiangjun temple，the correlation coefficient were 0.776，0.751，significant level of correlation p< 0.01，Cu elements in the desert zone and Pb elements had a significant positive correlation in Jiangjun temple，the correlation coefficient was 0.792.then the significant correlation elements were fitted，three fitting relation between Hg elements in Jiangjun temple and Cr elements in desert zone were the best，and the correlation coefficient R²=0.9491.The spatial interpolation of each element was obvious，and the high value area was mainly distributed in Wucaiwan and Jiangjun temple.