<p>This study performed a probabilistic analysis of the undrained lateral bearing capacity for a system of two circular piles at great depths. The adopted methodology combined finite element limit analysis with random field theory, considering the undrained shear strength as a random variable. The investigation systematically examined the influence of probabilistic parameters, notably the coefficient of variation of shear strength (CovSᵤ) and the spatial correlation length, on the mean bearing capacity factors of the pile group (µNₚ). The analysis also incorporated the effect of inter-pile spacing and adhesion factors at the soil-pile interface. The results revealed a significant inverse correlation between CovSᵤ and µNₚ for different spacing configurations. This research established a robust methodological framework for the probabilistic analysis of deep foundations, offering promising perspectives for their optimized design.</p>