Purpose
Natural gas leak from underground pipelines could lead to serious damage and global warming, whose spreading in soil should be systematically investigated. This paper aims to propose a three-dimensional numerical model to analyze the methane–air transportation in soil. The results could help understand the diffusion process of natural gas in soil, which is essential for locating leak source and reducing damage after leak accident.
Design/methodology/approach
A numerical model using finite element method is proposed to simulate the methane spreading process in porous media after leaking from an underground pipe. Physical models, including fluids transportation in porous media, water evaporation and heat transfer, are taken into account. The numerical results are compared with experimental data to validate the reliability of the simulation model. The effects of methane leaking direction, non-uniform soil porosity, leaking pressure and convective mass transfer coefficient on ground surface are analyzed.
Findings
The methane mole fraction distribution in soil is significantly affected by the leaking direction. Horizontally and vertically non-uniform soil porosity has a stronger effect. Increasing leaking pressure causes increasing methane mole flux and flow rate on the ground surface.
Originality/value
Most existing gas diffusion models in porous media are for one- or two-dimensional simulation, which is not enough for predicting three-dimensional diffusion process after natural gas leak in soil. The heat transfer between gas and soil was also neglected by most researchers, which is very important for predicting the gas-spreading process affected by the soil moisture variation because of water evaporation. In this paper, a three-dimensional numerical model is proposed to further analyze the methane–air transportation in soil using finite element method, with the presence of water evaporation and heat transfer in soil.