Leakage and Diffusion Law and Risk Assessment of Buried Natural Gas Pipelines Considering Soil Stratification and Permeability Difference

This study investigates methane leakage and diffusion from a buried high-pressure natural gas pipeline (8 MPa, 1000 mm diameter) using CFD simulations with the DES turbulence model. Based on homogeneous and layered soil models, the influences of soil porosity (0.46 to 0.54), particle size (10 μm to 100 μm), and soil stratification on the spatial and temporal characteristics of methane diffusion are systematically explored. The simulation results show that (1) methane diffuses from the leak hole to the surrounding soil in an ellipsoidal pattern, with the fastest diffusion speed along the pipeline’s axial direction. (2) In homogeneous soil, within the range of soil parameter values considered in this study, the absolute changes in risk assessment indices (FDR, GDR) caused by soil particle size were more significant; whereas the relative percentage changes in risk assessment indicators caused by soil porosity were more pronounced. (3) In layered soil, the permeability contrast between adjacent layers creates the permeability discontinuity interface effect. When a fine-grained or low-porosity layer overlies a coarse-grained layer, the upper layer acts as a hydraulic barrier, prolonging FDT from 130 s to 354 s while promoting significant horizontal spread at the interface. Conversely, a coarse-grained or high-porosity upper layer accelerates vertical breakthrough. These findings provide a scientific basis for risk assessment, monitoring site optimization, and emergency response planning, particularly in regions with heterogeneous stratified soils.