Probabilistic estimation of variogram parameters of geotechnical properties with a trend based on Bayesian inference using Markov chain Monte Carlo simulation

When geotechnical properties show a trend with spatial coordinates, estimation of a variogram model is a challenging task. In the previous studies, the trend-removal method based on ordinary least-squares approach has been commonly used. However, the obtained variogram is biased because the residuals are assumed to be statistically independent. In this study, the ability of Bayesian inference using Markov chain Monte Carlo (MCMC) simulation to estimate the variogram of geotechnical properties with a trend is explored using cone penetration resistance (q(c)) data of piezocone penetration tests (CPTU). The results show that the Bayesian inference method can estimate variogram parameters and the coefficients of trend function accurately. Based on the posterior variogram models, the predictive uncertainty and total uncertainty bounds are presented using kriging and sequential Gaussian simulation (SGS) methods. 96% validation points lie within the 95% confidence intervals of the total uncertainty based on 200 measurements of q(c) at NS31. The median and mean SSD of the prediction are 0.34 and 0.87, which is closer to the SSD criterion than the trend-removal method. The model uncertainty of the variograms, the predictive and total uncertainty of prediction all decrease as the sampling density increases at NS31 and NS12.