Response of buried flexible pipe-soil system is studied, through numerical analysis, with respect to deflection and buckling in a spatially varying soil media. In numerical modeling procedure, soil parameters are modeled as two-dimensional non-Gaussian homogeneous random field using Cholesky decomposition technique. Numerical analysis is performed using random field theory combined with finite difference numerical code FLAC 5.0 (2D). Monte Carlo simulations are performed to obtain the statistics, i.e., mean and variance of deflection and circumferential (buckling) stresses of buried flexible pipe-soil system in a spatially varying soil media. Results are compared and discussed in the light of available analytical solutions as well as conventional numerical procedures in which soil parameters are considered as uniformly constant. The statistical information obtained from Monte Carlo simulations is further utilized for the reliability analysis of buried flexible pipe-soil system with respect to deflection and buckling. The results of the reliability analysis clearly demonstrate the influence of extent of variation and spatial correlation structure of soil parameters on the performance assessment of buried flexible pipe-soil systems, which is not well captured in conventional procedures.