Abstract： Based on the random field theory, a random finite element model of coupled unsaturated consolidation in spatial random soils was established using the covariance matrix decomposition method and Monte Carlo simulation method. The effects of spatial variability of Young’s modulus and saturated permeability were investigated. The results show that the mean and standard deviation of the settlement are significantly affected by the coefficient of variation of Young’s modulus, and COVE. When COVE increases, the mean and standard deviation of the settlement increase, which means the uncertainty of the settlement are increased due to the spatial variability of Young’s modulus. When the correlation length of Young’s modulus δln E decreases, the standard deviation of excess pore water pressure and the surface settlement decrease due to the increased spatial averaging effect. The spatial averaging effect is more significant when the thickness of soil layer is smaller than δln E. In general, the effect of spatial variability of ks is much less significant than that of elastic modulus, which mainly affects the unsaturated consolidation in process.