杨奇勇 谢运球 罗为群 谷佳慧，2 曾红春，2

(1.中国地质科学院岩溶地质研究所岩溶生态系统与石漠化治理重点实验室， 桂林 541004；2.中国地质大学(北京)地球科学与资源学院， 北京 100083)

引言

土壤是农业生产的基础，也是生态环境的重要组成部分[1]。它是环境中重金属迁移、转化的重要媒介，既是重金属聚集的汇，也是向其他系统迁移的源[2]。土壤对污染物有一定的自净能力，但是这种自净能力是有一定限度的。随着工农业发展、矿山开采等，人类向土壤排放的污染物越来越多，当进入土壤中的污染物超过土壤承受能力时，土壤就会向外界输出污染物，从而危害人类健康[3-4]。因此，研究土壤重金属元素含量的空间分布并对其污染风险进行定量评估有着重要的意义。

地统计学方法在环境污染物的空间分布和风险评估研究方面得到了越来越广泛的应用[5-7]。普通克里格法(Ordinary Kriging, OK)能够对未采样区域实现最佳无偏性估计，是空间预测中最常用的方法。但是，在污染评估与风险管理的实际应用中，人们感兴趣的可能并不是空间某一位置污染物含量的具体数值，而是其超出或低于某一阈值的风险概率分布。指示克里格法(Indicator Kriging，IK)能给出单一变量在一定阈值限制下的风险概率，多变量指示克里格法(Multiple variable indicator Kriging，MVIK)能够将给定阈值的多个变量综合成一个变量，并绘出综合风险概率分布图[8-11]，因而在环境风险评估中有较好的应用前景。如GOOVAERTS等[8]应用IK 法绘制了Swiss Jura地区土壤Cd、Cu和Pb等重金属元……