This paper looked back on the current situation about simulation using PHREEQC of heavy metal transport in soil and the current situation of coupling with other value simulation software. On such a base, the feasibility, advantages and problems of heavy metal reactive solute transport in soil using PHREEQC were discussed.

Through the soil column experiments, the migration behavior of heavy metal copper ions in saturated sandy loam is analyzed, and it is found that the inflow solution with high concentration is conducive to the migration of copper ions. In the meantime, the software package of CXTFIT is utilized to conduct fitting for the outflow dynamics of copper ions with deterministic equilibrium model and non-equilibrium two-zone model, and predict the changing process of copper ion concentration at the different depths of soil column over the time. The results show that in the indoor saturated homogeneous conditions, the simulation accuracy of fitting parameters that applied deterministic equilibrium model is higher than that used the non-equilibrium two-zone model. Moreover, the influence of fixed water on the copper ion migration can be ignored.