In order to investigate the structural safety issues of oil and gas pipelines under the coupled hazards of permanent ground deformation (PGD) and corrosion, the mechanical response and failure process of complex corroded pipelines under landslides are explored by effectively mapping the ground deformation motion distribution function onto a 3D nonlinear finite element pipe-soil interaction model. Furthermore, the stress concentration effects under varying corrosion axial spacing, circumferential spacing, and corrosion geometric parameters are revealed in comparison with single-corroded and intact pipelines. The results show that: (1) Corrosion defects have a minor influence on pipe-soil relative displacement but significantly exacerbate stress concentration on the rear sliding surface; (2) The failure process of dual-corroded pipelines is significantly driven by axial spacing and shows a positive correlation with the maximum equivalent stress. At the yield critical point where landslide displacement approaches 0.7 m, the F1 coefficient is most intensely affected by the positive synergy between axial spacing and corrosion depth. Compared with intact pipelines, dual corrosion defects can amplify local stress by up to 2.24 times; (3) The influence of circumferential spacing on maximum stress is relatively limited, but small circumferential spacing enhances the interaction between defects and promotes circumferential penetration of the failure zone; (4) The stress concentration factor prediction formula established based on machine learning genetic programming-symbolic regression (GP-SR) achieves an R2 > 0.95, enabling rapid safety assessment of pipelines under multi-factor coupling. The research findings can provide theoretical support for integrity management and risk assessment of complex hazard-prone pipelines in geologically hazardous environments.