For many practical problems in environmental management, information about soil heavy metals, relative to threshold values that may be of practical importance is needed at unsampled sites. The Hangzhou-Jiaxing-Huzhou (HJH) Plain has always been one of the most important rice production areas in Zhejiang province, China, and the soil heavy metal concentration is directly related to the crop quality and ultimately the health of people. Four hundred and fifty soil samples were selected in topsoil in HJH Plain to characterize the spatial variability of Cu, Zn, Pb, Cr and Cd. Ordinary kriging and lognormal kriging were carried out to map the spatial patterns of heavy metals and disjunctive kriging was used to quantify the probability of heavy metal concentrations higher than their guide value. Cokriging method was used to minimize the sampling density for Cu, Zn and Cr. The results of this study could give insight into risk assessment of environmental pollution and decision-making for agriculture.