Diverse geostatistical tools such as kriging have been used to estimate the volume and spatial coverage of contaminated soil needed for remediation. However, many approaches frequently yield estimation errors, due to inherent geostatistical uncertainties. Such errors may yield over- or under-estimation of the amounts of polluted soils, which cause an over-estimation of remediation cost as well as an incomplete clean-up of a contaminated land. Therefore, it is very important to use a better estimation tool considering uncertainties arising from incomplete field investigation (i.e., contamination survey) and mathematical spatial estimation. In the current work, as better estimation tools we propose stochastic simulation approaches which allow the remediation volume to be assessed more accurately along with uncertainty estimation. To test the efficiency of proposed methods, heavy metals (esp., Pb) contaminated soil of a shooting range area was selected. In addition, we suggest a quantitative method to delineate the confident interval of estimated volume (and spatial extent) of polluted soil based on the spatial aspect of uncertainty. The methods proposed in this work can improve a better decision making on soil remediation.