The invention discloses a soil organic carbon density spatial heterogeneity-considered Kriging interpolation model. According to the Kriging interpolation model, influences, on common Kriging methods, of low-spatial heterogeneity and spatial abnormal values are reduced by combining land utilization information; and according to different trend term functions, the Kriging interpolation model comprises a virtual variable regression Kriging method, a mean centralized Kriging method and a median centralized Kriging method. By utilizing the model, high-precision SOCD prediction and mapping are realized, and finally, an optimum prediction method and a prediction result thereof are selected according to model evaluation indexes. The reliability, prediction rationality, prediction precision and decision coefficient R2 of the model are remarkably improved; and meanwhile, the model realizes high-precision SOCD spatial mapping, so that significance is provided for realizing precision agriculture and improving the agricultural production efficiency.