Firstly, using the MATLAB software to map out the contour map of eight heavy metals concentration, we conclude the spatial distribution of heavy metals in soil, then adopt the Nemerow comprehensive pollution index method to evaluate soil heavy metals pollution degree of each factor and know the five kind of areas and pollution level of eight heavy metals. We analysis the propagation characteristic of heavy metal pollutants, build a set of differential equations diffusion model, and finally determine the position of the pollution sources.