A traditional grid model for soil sampling may suffer from poor efficiency and low accuracy. With a nonferrous metal processing plant as the study area, a three-dimensional kriging interpolation model was built based on this plant's preliminary investigation data for arsenic (As), and a detailed survey sampling programme was proposed. The sampling density at the pollution interval of the surface soil was estimated by the coefficient of variation method, and the sampling depth was determined by the pollution interval of the vertical prediction results. The results showed that the encrypted soil sampling distribution optimisation method obtains greater pointing accuracy with fewer points. The sampling accuracy was 87.62% after optimising the depth of pointing. Moreover, this approach could save 66.13% of the sampling costs and 56.93% of the testing costs compared to a full deployment programme. This study provides a new and cost-effective method for predicting the extent of contamination exceedance at a site and provides valuable information to guide post-remediation strategies for contaminated sites.