运用普通克里格、泛克里格、协同克里格和回归克里格4种方法,结合由DEM获取的高程因子以及土壤全氮和阳离子交换量(CEC),预测了黑龙江省海伦市耕地有机质含量的空间分布。不同样点数量下海伦市土壤有机质含量的空间变异结构分析表明,样点数量多并不一定能够识别土壤有机质含量的结构性连续组分,最优化的布置采样点位置可能比单纯增加采样点的数量更重要。不同样点数量下的空间预测精度分析表明,以全氮为协同变量的协同克里格预测精度最高,能反映出海伦市土壤有机质含量空间分布的局部变异细节,并且能够解释超过50%的有机质方差;而普通克里格、泛克里格和回归克里格3种方法预测误差均较高,回归克里格并不一定能提高预测精度。获取与土壤有机质含量相关性高且完整覆盖研究区的辅助数据是提高有机质含量空间预测精度,降低采样点数量的重要途径之一。