Applies the Kriging interpolation method to draw the spatial concentration distribution of heavy metal elements, and employs geo accumulation index method to calculate the average value and grade frequency of heavy metal accumulation index of each functional district; Then, taking the example of the living area, analyzes the correlation of heavy metals in soil; Finally discusses the prediction model of heavy metal pollutants accumulation in soil and the feasibility of various heavy metal removal methods.