The excavation of a super-large foundation pit is a greatly risky project in coastal water-rich silt strata, and it is of great significance to establish a numerical model to analyze and predict the stability of the foundation pit retaining its structure for safe construction. Based on a geological investigation report and key node monitoring data, the model parameters have been modified. Subsequently, the random field theory and numerical analysis were employed to proceed with deterministic analysis and uncertainty analysis, respectively. Using the uniform distribution method, Gaussian distribution method and the covariance matrix method to generate a random field model, finite difference software was applied to analyze the impact of spatial variability of cohesion and the internal friction angle on structural deformation. The study shows that the overall distribution of axial force is small on both sides and large in the center, and the axial force is larger near the shaped region. Due to the principle of “lever”, there is a tendency for the horizontal displacement of surrounding piles to partially rebound from the pit when the bearing platform pit is excavated. The spatial variability of the internal friction angle and cohesion has an important influence on the numerical value of the enclosure structure and surface deformation, and the variation pattern is basically unchanged.