A lattice Boltzmann model is proposed to simulate the forming process of artificial frozen soil wall. The enthalpy method is applied to deal with the latent-heat source term, and the adjustable thermal diffusivity is utilized to handle the change of thermophysical parameters. The model is tested by the heat conduction with solid–liquid phase change in semi-infinite space, which shows a good consistence between the numerical and analytical solutions, and the mesh resolution has little effect on the numerical results. Lastly, the development of frozen soil wall is discussed when the freezing pipes are arranged in a square. The results show that the evolution of temperature field with time is closely related to the distance from the freezing pipe. For the soil near freezing pipe, the temperature gradient is larger, the soil temperature drops rapidly and freezes in a short time. The time history curve of temperature is relatively smooth. For the soil far away from freezing pipe, the temperature evolution curve has obvious multistage, which can be divided into four stages: cooling, phase change, partly frozen and completely frozen. The spacing of freezing pipes has a significant influence on the overlapping time of artificial frozen soil wall, and there is a power function relationship between them.