Improved Bayesian model updating of geomaterial parameters for slope reliability assessment considering spatial variability

In engineering practice, Bayesian model updating using field data is often conducted to reduce the substantial inherent epistemic uncertainties in geomaterial properties resulting from complex geological processes. The Bayesian Updating with Subset simulation (BUS) method is commonly employed for this purpose. However, the wealth of field data available for engineers to interpret can lead to challenges associated with the “curse of dimensionality”. Specifically, the value of the likelihood function in the BUS method can become extremely small as the volume of field data increases, potentially falling below the accuracy threshold of computer floating-point operations. This undermines both the computational efficiency and accuracy of Bayesian model updating. To effectively address this technical challenge, this paper proposes an improved BUS method developed based on parallel system reliability analysis. Leveraging the Cholesky decomposition-based midpoint method, the total failure domain in the original BUS method, which involves a low acceptance rate, is subdivided into several sub-failure domains with a high acceptance rate. Facilitated with an improved Metropolis-Hastings algorithm, the improved BUS method enables the consideration of a large volume of field data and spatial variability of geomaterial properties in the probabilistic back analysis. The results of an illustrative soil slope, involving spatially variable undrained shear strength, demonstrate that the improved BUS method is effective in simultaneously incorporating a substantial volume of field measurements and observations in the model updating process. Through a comparison with the original BUS method, the improved BUS method is shown to be useful for Bayesian model updating of high-dimensional spatially variable geomaterial properties and slope reliability assessment.