本研究采用地统计学与GIS相结合的方法,研究了四川省双流县农田土壤铜、锌、铬、铅含量的空间变异特征及其影响因子,并对其采样数量、污染状况、风险评估及超标可能性预测等方面进行了深入分析,取得了以下主要结论：(1)研究区内铜、锌、铅含量受到外源因素深刻影响,积累非常明显,值域分布范围较为广泛。铜、锌、铬均具有中等空间相关性,其空间变异受到结构性因子和随机性因子共同影响；而铅含量空间相关性很弱,其空间变异主要受随机性因子影响。各向异性下,四种重金属含量空间相关性均有一定程度提高。土壤铜、锌、铬含量分布趋势较为相似,均为西北高,东南低,且分别具有极显著的正相关性。铅含量在区域内的分布趋势并不明显,不同地貌类型下均有高值区分布。交互检验表明,铜、锌、铬含量插值精度较高而铅含量降低,主要受到其空间相关性的影响。(2)相对于常规随机抽样方法,协同克里格法确保了样点缩减后变量的空间相关性。在相同取样数量下,协同克里格法的预测精度要高于普通克里格法。在初始样点数量缩减10%的情况下,铜、铬含量协同克里格插值的预测精度仍高于初始623个样点的普通克里格插值,且两种方法得到空间分布图具有极高的相似性。(3)根据《食用农产品产地环境质量评价标准》中土壤环境质量划分标准限值,样点铜、锌含量单项污染指数多处于清洁等级,而铅含量单项污染指数和综合污染指数多处于尚清洁等级。面积统计表明,区域内土壤环境污染并不严重,多处于尚清洁水平。其中,铅含量单项污染指数处于超标等级的面积比例较大。潜在生态风险评价表明,铜、锌、铅潜在生态危害较小,均处于轻微生态风险等级。三种重金属含量超过临界值的概率均较低,其中,铅含量超标概率的值域范围最为集中。铜,锌含量超标概率分布的规律性相对较强,而铅含量超标概率分布较为离散。(4)土壤铜、锌、铬含量在三种地貌类型和三种成土母质中的变化顺序分别表现出一致性但影响程度的差异有所不同,其中,地貌类型是影响铜含量的重要因子,而成土母质对铜和锌含量的影响较深；有机肥使用量对农田土壤铜、锌、铬含量空间变异有重要影响；农业综合条件中农业基础条件对土壤铜、锌、铬含量空间变异的影响较大,而铅含量的空间变异主要受农药使用量和运输业规模的影响；缓冲区分析表明,土壤铅含量随着与道路距离的增加表现出先增加后降低的趋势。