Contaminant diffusion in layered soil is very complex. The diffusion process is simplified as a one-dimensional diffusion mode. The Laplace transformation and Laplace numerical inverse transformation are adopted. The proposed solution can be degraded to a single layer case. The obtained results are well consistent with those obtained by the classic theory. The correctness of the presented solution is verified by comparing the calculation results with the results generated by existing numerical software. One case study is conducted and the calculation results have a good agreement with measured data, which prove the reliability of the proposed method. The presented method can be used to solve more complex boundary conditions of contaminant migration.