Abstract: Aiming at the migration and diffusion of heavy metal pollutants in soil， the finite volume method is used to discretize the basic equations that heavy metal pollutants transport in two-dimensional steady water flow， and then the calculation model for migration in saturated soil is obtained. The dynamic process of heavy metal pollutants transporting in different conditions is simulated by Matlab to study the migration laws of heavy metal pollutants. The research shows that the concentration contours of heavy metal pollutants gradually changes from central parabolic into peripheral elliptical over time， and concentration decreases continuously in both horizontal and vertical directions， when the source of pollution always exists; however， the concentration contours of heavy metal pollutants mainly presents elliptical over time， and pollution range gradually becomes larger with concentration decreasing continuously and the maximum concentration at the ellipse center， when the source of pollution lasts for only 10 days.

Key words: saturated soil, heavy metal pollutants, migration, finite volume method, numerical simulation