刘 强， 冯 娜， 张娅娅， 刘雪媛

(天水师范学院 资源与环境工程学院, 甘肃 天水 741000)

基于GIS的天水市秦州区樱桃园土壤重金属元素空间分布与污染评价

刘 强， 冯 娜， 张娅娅， 刘雪媛

(天水师范学院 资源与环境工程学院, 甘肃 天水 741000)

[目的] 获取土壤重金属污染状况,为促进秦州区樱桃园产业的可持续发展提供科学依据，也为其他果树用地的土壤重金属污染防治提供借鉴。[方法] 以天水市秦州区樱桃园主产区的6个镇(太京、玉泉、藉口、关子、牡丹、杨家寺)作为研究点，运用空间分析的方法，研究表层土壤中的重金属污染的空间分布特征，又运用尼梅罗综合污染指数法和污染负荷指数法对土壤重金属环境质量及潜在生态风险进行评价。[结果] 通过应用GIS得到的空间分布图可知，秦州区樱桃园土壤重金属污染程度都较轻，东北部樱桃园土壤重金属污染较西南部污染严重。[结论] 樱桃园土壤重金属污染程度突出元素为Pb，而Zn，Cr，Cu属轻度污染水平，Cd污染最轻，5种重金属元素之间含量差别较明显。土壤重金属污染主要集中在太京镇、玉泉镇、藉口镇，这些区域靠近城市中心，受交通、人为因素影响密切，污染严重。

土壤重金属； GIS； 尼梅罗指数法； 污染负荷指数法； 空间分布

文献参数： 刘强， 冯娜， 张娅娅， 等.基于GIS的天水市秦州区樱桃园土壤重金属元素空间分布与污染评价[J].水土保持通报，2017,37(3)：258-263.DOI:10.13961/j.cnki.stbctb.2017.03.044； Liu Qiang, Feng Na, Zhang Yaya, et al. Spatial distribution and pollution evaluation of heavy metals in cherry orchard soil in Qinzhou District, Tianshui City based on GIS[J]. Bulletin of Soil and Water Conservation, 2017,37(3)：258-263.DOI:10.13961/j.cnki.stbctb.2017.03.044

土壤重金属污染(heavy metal pollution of the soil)是指由于人类活动[1-2]，土壤中的重金属元素在土壤中的含量超过背景值，过量沉积而引起的含量过高的现象[3-4]。近年来随着工业、农业和交通运输业的迅速发展，通过各种途径进入土壤中的有害重金属(Zn, Cd, Pb, Cr, Cu等)不断增加，对农产品造成了日益严重的污染和危害[5-6]。目前，随着全球食品污染问题日益严重和人们环境意识的不断提高[7]，国内外越来越重视对绿色食品和无公害果品土壤环境质量的监测与评价，土壤重金属污染及其危害已成为环境科学领域研究的热点问题[8-10]。

天水市秦州区是典型的半湿润半干旱温带气候，冬无严寒，夏无酷暑，雨量适中，光照充足，肥水中等，土层深厚。得天独厚的自然资源优势，特别适宜樱桃果品生产，是中国最佳的樱桃生产区域之一。通过多年的观察实践表明，秦州区所产樱桃风味独特、色泽艳美，其品质要超过其他樱桃主产区，是专家公认的西北黄土高原地区适宜于露地生产优质樱桃为数不多的地区之一。但是近几年由于秦州区樱桃种植业发展过快，不少果农为了追求最大的经济效益，在缺乏科学理论指导的情况下盲目扩大种植面积、过度使用农药化肥和薄膜，导致了秦州区樱桃园土壤的重金属含量超标，使得大部分樱桃的品质下降，致使果农的经济利益受损，并且过量的重金属元素通过食物链进入人体后将对健康产生极大的危害。因此，本文拟以天水市秦州区樱桃园主产地的6个镇为研究对象，实地采集土壤样本，通过应用地统计学方法，使用半变异函数定量研究樱桃园土壤重金属元素的空间变异特征，然后在此基础上开展污染评价，以期获取土壤重金属污染状况，从而为促进秦州区樱桃产业的可持续发展提供科学依据，也为其他果品用地保护提供借鉴意义。

1 材料与方法

1.1 研究区概况

甘肃省天水市秦州区位于甘肃省东南部，属秦巴山区西秦岭北缘，地处东经105°13′—106°01′，北纬34°05′—34°40′，地跨黄河、长江两大流域，地势西高东低，海拔1 130～2 700 m，土壤主要为冲击、洪积物形成的淤淀土、草甸土，经过开垦耕种熟化而形成以黄绵土、黑垆土。秦州区是典型的半湿润半干旱温带气候，雨量适中，光照充足，肥水中等，土层深厚，冬无严寒，夏无酷暑，气候温和，四季如春，素有“陇上小江南”之称，得天独厚的自然资源优势，使得秦州区成为优质的水果和蔬菜生产基地之一。

1.2 土壤样品的采集

选择秦州区樱桃园主产区的6个镇(太京、玉泉、藉口、关子、牡丹、杨家寺)作为研究点，2016年5月下旬对研究区进行土壤采样，根据每个镇樱桃园面积大小的不同确定采样数目，在樱桃园主产镇太京镇、玉泉镇采样点为20个，其余每个镇采样点为10个，共选取样本点80个。每个样点取0—20 cm的土壤表层，同时使用GPS仪记录采样点地理坐标，以及各采样田块的基本信息。采集后的样品分别装于采集袋中，带回实验室，分出杂物，风干磨碎，过100目的筛，分别装于广口瓶中保存。

1.3 测定方法

樱桃园土壤重金属用盐酸—氢氟酸—高氯酸消解，重金属Cr总量按照火焰原子吸收分光光度法测定，Pb和Cd总量按照石墨炉原子吸收分光光度法测定，Cu和Zn总量按照火焰原子吸收分光光度法测定。

1.4 数据处理

土壤中重金属元素的测定数据，通过应用K-S检验(柯尔莫哥洛夫—斯米诺夫检验)发现土壤Pb和Cr元素含量是不符合正态分布的，这种不符合正态分布的情况在实际经常遇到，这种情况下用传统的半变异函数方法计算出的半变异函数值是有偏差的，这是特异值造成数据有偏的主要原因。在绘制土壤重金属的空间分布图及面积统计时，还采用了地统计学GS+9.0，Surfer软件。

1.5 统计方法

采用半方差函数对采集的样品进行分析，地统计学是以区域化变量理论为基础，以半方差函数为基本工具，研究那些在空间分布上既有随机性又有结构性，或空间相关和依赖性的自然现象，该函数可用下式表示：

(1)

式中：r(h)——半方差函数； h——2个样本间的分离距离； z(xi)，z(xi+h)——随机变量在空间位置xi和xi+h上的取值； N(h)——在分离距离为h时的样本对总数。由r(h)对h作图可作出半方差函数图，该图由一批散点组成，最后根据离差平方和(RSS)最小原则，对试验半方差函数进行拟合，最终得到半方差函数的理论模型。

2 结果与分析

2.1 数据统计与分析

对土壤样本数据进行K-S检验，以算术平均值作为描述总体样本大小特征的统计。表1为秦州区樱桃园土壤重金属元素含量统计分析结果[11]。

表1 甘肃省天水市秦州区土壤重金属元素的描述性分析

从表1可以看出，5种重金属Pb，Cr，Cd，Zn，Cu的含量的最大值分别是40.19，70.45，0.11，103.88，27.50 mg/kg，除Cd以外，其他元素均超过甘肃省土壤背景值，表明樱桃园土壤中Cu，Zn，Pb，Cr有一定程度的积累。其中Zn元素的变异系数最大，离散程度较大。Cu的变异系数最小，离散程度较低。

2.2 土壤中重金属元素的地统计学分析

将所有的样点数据经过数据检验后，采用普通Kriging插值方法，分别用不同类型的模型进行半方差拟和，并选定了恰当的模型后进行Kriging插值。从而得到研究区域的重金属元素分布图[12-14]。

从表2可知，5种土壤重金属元素可分别拟合为高斯模型、高斯模型、指数模型、高斯模型、指数模型。土壤理化性质的块金值均比较大，说明土壤重金属元素的变化主要受结构性随机因素引起的空间异质性的影响。

表2 甘肃省天水市秦州区重金属分布的半方差函数模型

在没有人类活动的纯自然条件下，土壤重金属含量一般具有很强的空间相关性，其含量受气候、母质和地形等结构性因素的影响；而人类干扰如工农业生产则增加了随机因素，使得土壤重金属的空间相关性减弱。块金值与基台值之比称为基底效应，是表示样本间的空间变异程度的重要数据。该值越大，样本间的变异受随机因素影响越大。如果基底效应小于25%，表示变量的空间变异以结构性变异为主，具有较高的空间相关性；介于25%～75%为中等程度空间相关；大于75%时以随机因素为主，为空间弱相关[15-16]。由表2可知，Cu，Zn，Cd，Pb，Cr属于以随机因素为主，为空间弱相关，重金属含量受人类活动影响相对较大。

2.3 土壤中重金属的空间分布格局

利用ArcGIS地统计模块进行 Kriging 插值计算后,土壤重金属含量的插值结果如图1—5所示。

从图1可以看出：Pb污染程度由东北向西南递减；6个镇内部污染程度差别较大。东北部的玉泉镇大部分地区和太京镇北部地区污染最为严重，Pb元素含量38.11～40.19 mg/kg，而西部杨家寺镇的中部地区、中部的牡丹镇的部分地区污染最轻，Pb含量29.79～32.28 mg/kg。

图1 甘肃省天水市秦州区Pb元素空间分布

从图2可以看出，Cr污染程度大致由东向西递减；玉泉镇北部部分地区、太京镇北部部分地区以及藉口镇东部部分地区污染最为严重，Cr含量67.49～70.19 mg/kg。关子镇的部分地区和中部的牡丹镇的一小部分地区污染最轻，Cr含量59.42～62.15 mg/kg。

图2 甘肃省秦州区Cr元素空间分布

从图3可以看出，各地区Cd总体污染程度相当。玉泉镇北部一小部分地区、太京镇北部的一小部分地区、藉口镇中部的极小部分地区以及牡丹镇的极小部分地区污染最为严重，Cd含量0.09～0.11 mg/kg。玉泉镇北部的极小部分地区污染最轻，Cd含量0.04～0.05 mg/kg。

图3 甘肃省秦州区Cd元素空间分布

从图4可以看出：Zn污染程度由东北向西南递减；玉泉镇北部的一小部分地区、太京镇北部的一小部分地区、藉口镇的极小部分地区污染最为严重，Zn含量96.99～103.88 mg/kg。中部牡丹镇的一小部分地区污染最轻，Zn含量40.08～46.98 mg/kg。

从图5可以看出：Cu污染程度总体由东北向西南递减，其中玉泉镇北部部分地区、太京镇北部部分地区、藉口镇的极小部分地区污染最为严重，Cu含量26.88～27.50mg/kg。牡丹镇和关子镇的一小部分地区污染最轻，Cu含量在20.37～21.19 mg/kg，东北部污染较西南部污染严重。

从重金属元素空间分布格局来看，东北部土壤重金属污染较西南部污染严重，主要是因为东北部樱桃园的种植面积比西南部的种植面积大，农药和化肥的使用率高、使用量大，并且这些地区靠近河流，污染的水体通过灌溉进入到樱桃园土壤中，造成土壤重金属的大量积累，对土壤产生了污染。

图4 甘肃省秦州区Zn元素空间分布

图5 甘肃省秦州区Cu元素空间分布

3 樱桃园土壤重金属的尼梅罗指数和污染负荷指数法评价

3.1 尼梅罗指数法

尼梅罗指数法能够较为全面地评价土壤重金属污染水平，可以综合反映重金属对区域的污染程度，不仅反映单个采样点的重金属污染程度，而且可以对某一区域的土壤污染状况进行综合评价[17]。它是一种通过单因子污染指数得出综合污染指数的方法，能够较全面地评判其重金属的污染程度，是人们在评价土壤重金属污染时运用最为广泛的综合指数法。其计算公式为：

(1)

式中：Pi——单向污染指数；Ci——污染物实测值(mg/kg)；Si——根据需要选取的评价标准(mg/kg)。

污染土壤的重金属元素往往不只一种，所以还需要考虑不同重金属元素对于土壤的污染程度，可以用尼罗梅污染指数来衡量，其数学表达式为：

(2)

式中：Pi——单向污染指数的平均值；Pimax——最大单项污染指数。

计算土壤重金属按照尼梅罗指数评价标准划分污染等级。清洁无污染(P≤0.7)；无污染(P≤1)；尚清洁(警界线)(0.73)。

3.2 污染负荷指数法

污染负荷指数是由Tomlinson于1980年提出的，由评价区域所包含的多种重金属成分共同构成。该方法不仅可以反映单个采样点的重金属污染程度，而且可以对某一区域的土壤污染状况进行评价。其参照指标由评价区域的多种重金属共同构成，因此能反映各个重金属对区域污染的贡献程度，还可进一步反映各个重金属污染的时空变化特征[18-19]，其计算公式为：

(3)

式中：Cfi——重金属i的污染系数；Ci——重金属i的测定值；Coi——重金属i的参比值，采用甘肃省土壤背景值；m——参与评价的重金属种类数。

某一区域的污染负荷指数计算公式：

(4)

式中：pa——某污染区域的污染负荷指数；n——该污染区域所含采样点数。

3.3 樱桃园土壤重金属的尼梅罗指数和污染负荷指数法的评价结果

3.3.1 尼梅罗综合污染指数法评价结果 通过应用尼梅罗综合污染指数法计算，秦州区樱桃园土壤重金属得出到以下结果(表3)。

通过表3可以发现秦州区樱桃园土壤重金属测量值与甘肃省土壤背景值进行对比后，可以得出以下结论：

表3 尼梅罗综合污染指数(Pi)评价结果

(1) 除Cr外，Pb，Cr，Zn和Cu的平均含量均高于背景值，分别达到背景值的1.74，1.03，1.10和1.03倍。

(2) Pb的单项污染指数Pi均为1

(3) Pb元素在秦州区樱桃园土壤中污染最为严重，而Cd清洁无污染。

(4) 太京镇、藉口镇、玉泉镇、关子镇综合污染指数均为1

3.3.2 污染负荷指数法评价结果 通过应用污染负荷指数法计算，秦州区樱桃园土壤重金属得出到以下结果(表4)。

(1) 各乡镇土壤重金属污染程度由重到轻依次为：太京镇>玉泉镇>藉口镇>关子镇>牡丹镇>杨家寺。太京镇、藉口镇、玉泉镇土壤重金属的污染负荷指数(PLI)值为1≤P<2，为中等污染水平；关子镇、牡丹镇、杨家寺污染PLI值为P<1，属无污染。

(2) 最小值PLImin=0.84出现在牡丹镇和杨家寺；最大值PLImax=1.15出现在太京镇，两者相差1.37倍。

表4 甘肃省天水市秦州区污染负荷指数(PLI)评价结果

4 结 论

(1) 樱桃园土壤重金属污染受到人为因素的影响较大，重金属元素呈现不同程度的富集，且在空间上呈现出一定的分布规律。总体来讲，重金属元素污染程度为：Pb>Zn>Cr>Cu>Cd；从各乡镇污染程度来看：太京镇>玉泉镇>藉口镇>关子镇>牡丹镇>杨家寺。

(2) 从尼梅罗综合污染评价结果来看：樱桃园土壤重金属污染程度突出元素为Pb，而Zn，Cr，Cu属轻度污染水平，Cd污染最轻。5种重金属元素之间含量差别较明显，其中Pb，Cd的含量差别最大；单就Pb元素而言，其影响范围广泛分布在樱桃园种植的全区，而Cr，Zn，Cu重金属元素重点污染在太京镇、玉泉镇，藉口镇。

(3) 从污染负荷指数法评价结果来看：太京镇土壤重金属污染负荷指数高达1.15，为中等污染水平；玉泉镇、藉口镇、关子镇次之，同为中等污染水平；牡丹镇、杨家寺未受到污染，且污染主要集中在太京镇、玉泉镇、藉口镇，是因为这些区域靠近城市中心，受交通、人为因素影响密切，污染重，而关子镇、牡丹镇、杨家寺离城市中心较远，受影响小，污染轻。

(4) 尼梅罗指数法和污染负荷指数法得出结论基本一致，但在部分结论上稍有差异，2种评价结果均显示秦州区太京镇、玉泉镇污染严重，这是由于这2个镇种植樱桃园历史悠久，农业基础设施完善，适宜樱桃园生产，果农为提高樱桃园产量大量使用化肥农药，导致土壤中重金属元素不断富集。比较2种评价结果可得，关子镇在尼梅罗指数法评价中属于轻污染水平，而在污染负荷指数评价中属于无污染，这是由于2种评价方法的侧重点不同，评价分级不同，进而反映出评价结果的差异。

(5) 通过应用GIS开发技术得到的空间分布图可知，秦州区樱桃园土壤重金属污染程度都较轻。东北部樱桃园土壤重金属污染较西南部污染严重，主要是因为东北部樱桃园的种植面积比西南部的种植面积大，农药和化肥的使用率高、使用量大，并且这些地区靠近河流，污染的水体通过灌溉进入到樱桃园土壤中，造成土壤重金属的大量积累，对土壤产生了严重的污染。

(6) 基于GIS的樱桃园土壤重金属元素空间的空间分布，在采样布设合理的前提下，评价结果可信。

[1] 李冰茹,王纪华,马智宏,等. GIS土壤重金属污染评价中的应用[J].测绘科学,2015,40(2):119-123.

[2] 王幼奇,白一茹,王建宇.基于GIS的银川市不同功能区土壤重金属污染评价及分布特征[J].环境科学,2016,37(2):710-716.

[3] 罗娇赢. 基于GIS与地统计学的哈尔滨市城郊菜地土壤重金属元素空间变异分析[J]. 黑龙江生态工程职业学院学报,2011(3):35-37.

[4] 蔡小冬,孙成胜,蔡立群.基于改进AHP与GIS的白银区耕层土壤重金属污染评价[J].中国农学通报,2014,30(32):194-200.

[5] 郝丽虹,张冬明,吴鹏飞,等.海南岛农用地土壤重金属含量的空间分布[J].中国生态农业学报,2009,17(2):230-234.

[6] 刘庆,杜志勇,史衍玺,等.基于GIS的山东寿光蔬菜产地土壤重金属空间分布特征[J].农业工程学报,2009,25(10):258-263.

[7] 杨东,刘强.基于GIS和地统计学的张掖市甘州区土壤全氮、有机质的空间变异特征分析[J].土壤通报,2011,42(3):593-597.

[8] 陈峰,尹春芹,蒋新,等.基于GIS的南京市典型蔬菜基地土壤重金属污染现状与评价[J].中国环境监测,2008,24(2):40-44.

[9] 李健,郑春江.环境背景数据手册[M].北京:中国环境学出版社,1989.

[10] 罗永清,陈银萍,陶玲,等.兰州市农田土壤重金属污染评价与研究[J].甘肃农业大学学报,2011,46(1):98-104.

[11] 范拴喜,甘卓亭,李美娟,等.土壤重金属污染评价方法进展[J].中国农学通报,2010,26(17):310-315.

[12] 尤冬梅,王纪华,马智宏,等.北京郊区小尺度农产品产地土壤重金属污染性评价[J].上海农业学报,2011,27(3):89-93.

[13] 陈学民,朱阳春,董会平,等.天水苹果园土壤重金属污染调查与评价[J].兰州交通大学学报,2011,30(3):132-135.

[14] 陆安详,王纪华,潘渝春,等.小尺度农田土壤重金属的统计分析与空间分布研究[J].环境科学,2009,28(7):1578-1583.

[15] 戴宇,杨重法,郑袁明.土壤—植物系统中的铬的环境行为及其毒性评价[J].环境科学,2009,30(11):3432-3440.

[16] 徐宁,俞燕芳,毛平生,等.桑树修复土壤重金属污染的研究进展[J].农学学报,2015,5(1):37-40.

[17] 朱爱萍,陈建耀,江涛,等.北江流域横石河—翁江沿岸土壤重金属污染特征分析[J].中国环境科学,2015,35(2):506-515.

[18] 刘勇,岳玲玲,李晋昌,等.太原市土壤重金属污染及其潜在生态风险评价[J].环境科学学报,2011,31(6):1285-1293.

[19] 钟晓兰,周生路,李江涛,等.长江三角洲地区土壤重金属污染的空间变异特征:以江苏省太仓市为例[J].土壤学报,2007,44(1):33-40.

Spatial Distribution and Pollution Evaluation of Heavy Metals in Cherry Orchard Soil in Qinzhou District, Tianshui City Based on GIS

LIU Qiang, FENG Na, ZHANG Yaya, LIU Xueyuan

(CollegeofResourcesandEnvironmentalEngineering,TianshuiNormalUniversity,Tianshui,Gansu741000,China)

[Objective] The status of soil heavy metal pollution was researched to provide scientific basis for promoting the sustainable development of cherry orchard industry in Qinzhou District, and for soil heavy metal pollution prevention and control in other fruit land. [Methods] Six towns of cherry orchard production area in Qinzhou District of Tianshui City(Taijing, Yuquan, Xikou, Guanzi, Mudan and Yangjiasi) were taken as the research sites. Where the spatial distribution characteristics of heavy metal pollution in surface soil was illustrated using method of spatial analysis. And the heavy metal environmental quality and its potential ecological risk were evaluated by Nemerow comprehensive pollution index and pollution load index. [Results] Overall, the soil heavy metal pollution of cherry orchard in Qinzhou District was light, but the soil heavy metal pollution in the northeast region was serious than that in the southwest region of the research sites. [Conclusion] The prominent heavy metal element in the cherry orchard is Pb; and Zn, Cr, Cu pollution also existed but belong to light pollution level; Cd pollution is the most light. It is obvious different among the contents of the five heavy metal elements. The heavy metal polluted towns mainly concentrated in Taijing, Yuquan and Xikou County that are adjacent to downtown, for the reason of traffic and intensive human interference, are seriously polluted.

soil heavy metal; GIS; Nemerow index; pollution load index; spatial distribution

2016-09-27

2016-10-28

天水市科技支撑项目“天水市发展农业循环经济评价研究”阶段性成果(2010476); 天水师范学院中青年教师科研项目“基于GIS的樱桃园土壤重金属元素分析与评价”(TAS1615)

刘强(1983—),男(汉族),甘肃省兰州市人,硕士,讲师,主要从事土壤生态与环境研究。E-mail:guangmingliu1983@163.com。

A

1000-288X(2017)03-0258-06

X825