Mining activities play a significant role in environmental pollution by producing large amounts of tailings which comprise heavy metals. The impressive increase in mining activities in recent decades, due to their high influence on the industry of developing countries, duplicates the need for a substantial effort to develop and apply efficient measures of pollution control, mitigation, and abatement. In this study, our objective was to investigate the effect of simulation of the leachate, pH and inflow intensity of transport of $Pb^{2+}$, $Zn^{2+}$, and $Cd^{2+}$ through Lakan lead and zinc plant tailings, in Iran, and to evaluate the modeling efficiency by comparing the modeling results and the results obtained from previous column studies. We used the HP1 model and the PHREEQC database to simulate metals transport through a saturated soil column during a 15 day time period. The simulations assumed local equilibrium. As expected, a lower pH and inflow intensity increased metal transport. The retardation of heavy metals followed the order $Zn^{2+}$ > $Pb^{2+}$ > $Cd^{2+}$ and the removal concentrations of Cd, Pb, and Zn at the inflow intensity critical scenario, and Cd and Pb at inflow acidity critical scenario exceeded the allowable EPA and Iranian's 1053 standard thresholds. However, although the simulation results generally agreed well with the results of the column study, improvements are expected by using multi-dimensional models and a kinetic modeling approach for the reactions involved. The results of such investigations will be highly useful for designing preventative strategies to control reactive transport of hazardous metals and minimize their environmental effects.