The zoning of copper, nickel and lead heavy metals was investigated by using Kriging method in GIS environment using circular, spherical, exponential and Gaussian variorums. In addition, one-dimensional Hydrus modeling of water flow and heavy metals in the soil environment was simulated up to 50 cm depth for a 210-day period and the concentration of heavy metals to the root .depth was simulated. Distribution of lead element in soil surface with spherical model showed that its variation was in the range of 20 to 70 mg / kg. These values were 50-60 mg / kg for copper and 30 mg / kg for nickel. Investigation of heavy metals in the soil using the Hydrus model showed that the simulated value at the initial 0-15 cm depth has the highest value and at lower depths is decreased. Comparison of the concentrations of these elements with the standard allowed by the WHO showed that the lead element in this region was higher than the permissible level.