For underground water-sealed oil storage, the spatial variability of the surrounding rock has a significant impact on the water-sealed effect of a water curtain system. This study presents a methodology for the probabilistic analysis of water curtain performance in underground oil storage, considering the spatial variability of hydraulic conductivity of the surrounding rock based on field data. Anisotropic random fields representing the spatial variability of hydraulic conductivity were established through spatial statistical analysis of field data and introduced into the finite element model of underground oil storage for water-sealed reliability analysis. The water-sealed performance of different water curtain system schemes was studied using Monte Carlo simulation (MCS). The results showed that the difference between the horizontal spatial correlation and the vertical spatial correlation of the surrounding rock has a significant impact on the water-sealed effect of the water curtain system. An excessively large pressure of water curtain boreholes provided a small contribution to improving water curtain performance. The distance between the water curtain holes and the caverns had the less significant affecting the water-sealed reliability of the storage cavern. Finally, the optimal design of the water curtain system is discussed. This study provides valuable insights and a theoretical basis for the optimisation of water curtain system design parameters for underground water-sealed oil storage.