Leakage accidents of buried pipelines have become increasingly common due to the prolonged service of some pipelines which have been in use for more than 150 years. Therefore, there is an urgent need for accurate prediction of pollution scope to aid in the development of emergency remediation strategies. This study investigated the distribution of a light non-aqueous phase liquid (LNAPL) in soils containing gas and water through numerical simulations and laboratory experiments. Firstly, a three-dimensional porous medium model was established using ANSYS FLUENT, and for the first time, the distribution of gas and groundwater in soil environments was simulated in the model. Subsequently, the distribution of the three phases of diesel, gas, and water in soil was studied with different leakage velocities and it was found that the leakage velocity played a significant role in the distribution. The areas of diesel in soils at 60 minutes were 0.112 m2, 0.194 m2, 0.217 m2, and 0.252 m2, with corresponding volumes of 0.028 m3, 0.070 m3, 0.086 m3, and 0.106 m3, respectively, for leakage velocities of 1.3 m/s, 3.4 m/s, 4.6 m/s, and 4.9 m/s. Calculation formulas for distribution areas and volumes were also developed to aid in future prevention and control strategies under different leakage velocities. The study also compared the distribution areas and volumes of diesel in soils with and without groundwater, and it was found that distribution scopes were larger in soils containing groundwater due to capillary force. In order to validate the accuracy of the numerical simulation, laboratory experiments were conducted to study the diffusion of oil, gas, and water under different leakage velocities. The results showed good agreement between the experiments and the simulations. The research findings are of great significance for preventing soil pollution and provide a theoretical basis for developing scientifically sound soil remediation strategies.