Deep soil mixing is a widely used ground improvement technique for the enhancement of bearing capacity and settlement behaviour of soft soil deposits. Literature suggests that the strength of deep cement mixed column (DCM) improved ground is not homogenous based on tests conducted on core samples collected from the field. However, due to the challenges associated with the strength and settlement determination, heterogeneity is disregarded in the approaches used in most of the present research. The current investigation focuses on how spatial variability affects the strength characteristics of DCM column-improved ground by considering a unit cell. To account for the spatial variability, the material properties such as cohesion (c) and angle of internal friction (φ) of soft soil and DCM column are varied with respect to coefficient of variation (COV) and correlation length (ϴ) with the help of a user-defined soil model (UDSM) using PLAXIS 2D. Stochastic analysis is conducted by varying the COV values from 0 to 0.4 to analyse the effect of spatial variability on the load-bearing behaviour of the DCM column-improved soil in the unit cell model. The mean bearing capacity for a unit cell model is determined considering the variability of properties, and the values are compared with deterministic analysis.