徐克全 ， 金立新， 郑勇军， 张 华， 包雨函 ， 时章亮

(1.四川省地质调查院，成都 610081；2.稀有稀土战略资源评价与利用四川省重点实验室，成都 610081)

0 引言

土壤元素含量在空间上具有自相关性，也存在多元素间交叉相关[1-2]，在空间分布上总体受到地质背景的控制，局部可能会受到人为因素的影响。为了获取特定空间范围数据，并对特定空间对象进行评价，需要对有限的空间数据进行插值处理。常用到反距离加权(IDW)、半变异函数的克里金(Kriging)等方法进行插值。但不同方法选取，或半变异函数中的块金、变程、偏基台值的设置，都将对插值结果造成影响，最终影响土壤质量地球化学综合评价面积统计的准确性，不利于土地风险管控和科学利用。笔者通过R语言(运用到sp、gstat、geoR、ggplot2、raster、intensity.analysis等R语言包)，使用麻城镇和摩尼镇1:10 000土地质量地球化学调查评价结果的表层土壤数据[3]，考虑元素空间自相关性可能受到地质背景、土地利用类型、元素的地球化学亲和性、空间距离的影响，评价单元对应赋值方法分别采用地层属性分层克里金法插值、土地利用类型分层克里金法插值、多变量协同克里金法插值、反距离加权平均插值法；考虑“有样”评价单元赋原值以及“无样”评价单元插值方法的更合理性，也是主流方法，增加对“有样”单元均采用局部加权平均(权重为公式3)赋值，而对“无样”单元插值分别采用“地层属性分层克里金插值”、“土地利用类型分层克里金插值”、“多变量协同克里金插值”和“反距离加权平均法(权重公式(2))插值”，共8种赋值方法进行养分单指标和土壤质量地球化学综合评价。……