In this study, a simple probabilistic approach using equivalent coefficient of consolidation (c_e) was proposed to consider the spatial variability of coefficient of vertical consolidation (c_v), and the effect of the probability distribution of coefficient of consolidation on degree of consolidation in heterogeneous soil was investigated. The statistical characteristics of consolidation coefficient were estimated from 1,226 field data, and four probability distributions (Normal, Log-normal, Gamma, and Weibull) were applied to consider the effect of probability distribution. The random fields of coefficient of consolidation were generated based on Karhunen-Loeve expansion. Then, the equivalent coefficient of consolidation was calculated from the random field and used as the input value of consolidation analysis. As a result, the probabilistic analysis can be performed effectively by separating random field and numerical analysis, and probabilistic analysis was performed using a Latin hypercube Monte Carlo simulation. The results showed that the statistical properties of c_e were changed by the probability distribution and spatial variability of c_v, and the probability distribution of c_v has considerable effects on the probabilistic results. There was a large difference of failure probability depend on the probability distribution when the autocorrelation distance was small (i.e., highly heterogeneous soil). Therefore, the selection of a suitable probability distribution of c_v is very important for reliable probabilistic analysis of consolidation.