Abstract: Different sample well densities and interpolation methods significantly affect the spatial variability assessment of groundwater quality.The methods of inverse distance weighted (IDW),ordinary Kriging and radial basis function (RBF) were applied to predict the spatial concentration of 3 chemicals (chloride,sulphate,sodium) in the groundwater of Heze City in order to select an appropriate interpolation method producing a reliable map.The mean errors,mean relative errors and root mean square errors were compared to assess the best-fit interpolation method.The results show that in the study area with evenly distributed groundwater sample wells,the ordinary Kriging method gives optimum prediction results for all 3 involved chemicals.The areas of 5 groundwater quality grades were estimated for different sample well densities namely:5 wells,4 wells and 3 wells/10³km².The estimated grade areas with Kriging method have less differences when the prediction results are interpolated with the density of 4 wells per 10³km² or above.The estimated areas of Ⅰ~Ⅲ grade for chloride,sulphate,sodium in the groundwater of Heze City in 2019 are 6 868.0km²,8 125.75km² and 4 770km² while the areas of over the third grade limitations are 5 299.5km²,4 041.75km² and 7 397.50km² respectively.A thorough understanding of the groundwater quality status in the case study area is a significant tool for water resources management and decision-makers.