The stability problem of a tunnel heading in clay remains a significant challenge in geotechnical engineering. Specifically, when considering the spatial variability of the soil, the stability factor may be influenced by geographically random fields. This study investigates the effect of random fields on a probabilistic analysis of a tunnel heading in undrained clay. The study assumes that the undrained shear strength of the clay increases linearly with depth due to a strength gradient factor. The random adaptive finite element limit analysis is employed to calculate the stability numbers for tunnel headings. Nonstationary random fields with varying vertical correlation lengths are simulated using Monte Carlo simulation technique. The stability analysis takes into account geometry parameters (i.e., cover depth ratio) and nonstationary random field of undrained shear strength parameters. (i.e., strength gradient, coefficient of variation, and vertical correlation length). The results of tunnel face stability using random adaptive finite element limit analysis have also been utilised to assess the probability of design failure over a practical range of deterministic factors of safety. In the context of probabilistic failure analysis, the failure mechanism resulting from varying vertical correlation lengths could influence the probability of design failure. The findings of this study can be of significant interest to tunnel engineering practitioners during the design phase of tunnel heading projects.