In this paper, the supercritical carbonation process of cement-based materials is modelled by introducing a random porosity field to simulate the heterogeneous geometry of the carbonation profile. The suitability of two different random fields of porosity, based on the probability density function (PDF) and the ellipsoidal autocorrelation function (EAF) methods, are investigated, respectively, in simulating the distribution of porosity in cement mortar. After incorporating the above random fields into an established supercritical carbonation model, it is found that with some modifications, the EAF method with consideration of spatial correlation produces better simulation of the irregularities of the carbonation zones that have been observed from experimental results. It is also found that for given average porosity and coefficient of variation, the predicted average and maximum carbonation depths have much smaller coefficients of variation. The validated EAF supercritical carbonation model is then used in parametric studies that are conducted to assess the effect of various factors on the carbonation depth of the chemical process.