土壤作为一种重要的自然体，是关系到人类生存发展的重要资源，土壤中元素含量、空间结构以及空间分布特征可以作为区分自然环境受人类活动影响强度的重要依据，已经引起相关领域学者的广泛关注^[1-3]。研究表明重金属可通过施肥、污水灌溉、大气沉降、动物粪便和生物体利用等途径进入土壤，导致土壤中重金属含量超标，进而通过食物链和生物富集作用危害人体健康^[4-6]。

近些年来，国内外学者从不同角度分析了土壤中元素特征，并着重研究危害较大的重金属元素，通过BCR法^[7-9]分形态分析，多元统计分析法^[10]分析其来源，地统计分析法^[11-12]分析其空间特征，并评价其潜在生态风险。大多数研究集中于经济发展水平较高、受人类活动影响强度大的地区。我国西北干旱区地处内陆，随着西部大开发的进程不断深入，工农业生产活动不断加强，人类活动对土壤的干扰日益突出。已有对西北干旱区的绿洲表层土壤的重金属进行研究^[13-15]，海米提等^[16]对焉耆盆地表土的6种重金属的来源、分布及风险进行了研究，本文在此基础上从全元素角度分析，整体把握流域内绿洲土壤元素空间特征及重金属污染特征。

本文选取天山南麓博斯腾湖流域绿洲区为研究区域，通过数理统计学法和地统计学方法分析重金属的含量特征、空间变异特征、空间分布特征，并用潜在风险指数法、内梅罗污染指数法结合新疆土壤背景值评价在西部大开发不断深化过程下，干旱区绿洲土壤重金属含量与污染状况，为该区域的生态环境保护和资源的可持续开发利用提供科学依据。

博斯腾湖流域(40°25′~43°21′ N，82°57′~90°39′ E)位于新疆维吾尔自治区巴音郭楞蒙古自治州(简称巴州)境内，地处塔里木盆地和塔克拉玛干沙漠东北缘。流域总面积为7.7×10⁴ km²，开都河是流域的主要水源地，发源于天山中段依连哈比尔尕山南坡，高山区河流终年积雪，流经巴音布鲁克，最后注入焉耆盆地的博斯腾湖。博斯腾湖是流域内最大的天然集水区，水域辽阔，东西长约55 km，南北宽约20 km，水位变化较大，水位1 048.75 m，水面面积1 002.4 km²，平均水深8.8 m，最大水深17 m^[17]。孔雀河为博斯腾湖的出湖河流，流经库尔勒市。流域内还有黄水沟、清水河、乌什塔拉河等，开都河中下游流经焉耆回族自治县、博湖县，黄水沟流经和静县，清水河流经和硕县，承担着流域内居民生产生活用水任务。本研究区位于博斯腾湖流域焉耆盆地绿洲区(图 1)，辖和静县、和硕县、焉耆县、博湖县，区域内有G314、G218、S206、S325、S306纵横交错，是主要的人类活动区域，工农业发展水平较高。

图 1 研究区及采样点分布示意图 Figure 1 Distribution of sampling sites in Bostan Lake Basin

图选项

2016年6月在博斯腾湖流域内西边绿洲区等设置10 km×10 km的网格，在北部按照20 km的间距沿对角采样，总共77个采样点，主要涉及农田、林地、草地等绿洲土壤。利用环刀采取表层0~20 cm土壤，标记编号并装入塑料袋带回实验室，具体采样点位置见图 1。样品经自然风干后去除植物残体和碎石，研磨过100目筛，经105 ℃烘箱烘干后，取0.1~0.2 g样品于消解罐中，用去离子水润湿后，加入2 mL HCl、4 mL HNO₃、2 mL HF-H₂O，在德国Berghof MWS-3微波消解装置中进行消解。升温至120 ℃，保持1 min，然后再升温至160 ℃，保持5 min，最后再升温至195 ℃，维持25 min。消解后冷却至80 ℃，将消解液转移到聚四氟乙烯坩埚内，置于电热板上，逐渐升温加热至2 mL左右，冷却后，加入1%HNO₃定容至25 mL，摇匀待测。用电感耦合等离子体原子发射光谱仪(美国Leeman Labs Profile，ICP-AES)测得Al、Ca、Fe、Mn、V、Zn、As、Mo、Sb、Tl共10种元素含量；用电感耦合等离子体质谱仪(美国Agilent Technologies，ICP-MS)测得Cr、Co、Ni、Cu、Cd、Pb共6种元素含量；Hg采用AFS-920型双道原子荧光仪测定其含量；采用美国SPEX CertiprePTM Custom Assurance Standard多元素标准溶液，中国土壤标准物质GSS-12作为标准参考物质，测试完后取20%进行重复测定，使得误差小于5%。

元素含量数据的描述性统计分析、正态分布检验(k-s检验)、等数理统计分析均在SPSS 22.0软件中完成。半方差函数模型在Gs⁺9.0中完成，反距离加权插值法(IDW)在ArcGIS 10.3地统计模块中进行。

内梅罗污染指数是在单因子污染指数法的基础上发展而来的，首先计算单个重金属的污染指数，其计算公式如下：

(1)

式中：P_i为元素i的污染指数值，C_i为土壤样品中元素i含量的实测值，S_i为元素i的背景值，本文中所有重金属的背景值均选取新疆土壤环境背景值^[18]。在单因子污染指数法的基础上计算内梅罗综合污染指数，是目前比较常用的一种评价方法^[19-20]。计算公式为：

(2)

式中：P为采样点的内梅罗综合污染指数，P_{i max}为因子i的最大污染指数，P_{i ave}为因子i的平均污染指数。分成安全(P≤0.7)、警戒(0.7＜P≤1)、轻微污染(1＜P≤2)、中等污染(2＜P≤3)、重污染(P＞3)5种水平。

潜在生态风险指数法(RI)是瑞典著名化学家Hakanson于1980年提出的^[21]。该方法因综合考虑了区域背景值、重金属的生态毒理学特征等而被广泛运用于水体、水环境沉积物和土壤中重金属的潜在生态风险评价^[22-24]。其计算公式如下：

(3)

式中：C_fⁱ为重金属i的污染指数；Cⁱ为重金属i的实测含量；C_nⁱ为重金属i的土壤背景值；E_rⁱ为第i种重金属的单因子潜在生态危害指数；T_rⁱ为第i种重金属的毒性系数。Hakanson潜在生态风险评价标准是基于As、Hg、Cd、Cr、Cu、Pb、Zn和PCB 8种元素，与本文所研究的重金属元素种类不完全相同，采用其原有的分级标准会产生较大的偏差^[23]，因此本文参考Fernandez等^[25]的方法(表 1)。

表 1 潜在生态风险指数法分级标准 Table 1 Classification criteria of the potential ecological risk index

表选项

博斯腾湖流域77个表层土壤样品中17种元素含量统计结果见表 2。结果表明，17种元素平均含量顺序为Ca＞Al＞Fe＞Mn＞V＞Zn＞Cr＞Ni＞Cu＞Pb＞As＞Co＞Sb＞Mo＞Tl＞Cd＞Hg。与新疆土壤元素背景值^[18]相比，Al、Ca、Fe含量远远超过背景值，Cr、Cd、Hg平均含量略高于背景值，其他11种元素平均含量均低于背景值。其中，Co、Mo元素含量未超标，12种元素超标率排序为Cd＞Hg＞Cr＞V＞Sb＞As＞Zn=Ni＞Tl＞Pb＞Cu＞Mn，Cd、Hg超标率达77%、62%，Tl、Pb、Cu、Mn的超标率低于10%。

表 2 流域表土元素含量统计结果 Table 2 Statistics results of element in topsoil of the Bostan Lake Basin

表选项

变异系数(CV)是各元素标准偏差与平均值的百分比，反映了不同采样点元素含量的离散程度，受人类活动影响程度不同，元素的空间分布差异较大，CV越大离散程度越高，反之则越低^[26-27]。一般可将样本的变异程度分为3级：CV＜10%为弱变异性，10%≤CV≤100%为中等变异性，CV＞100%为强变异性。从表 2可看出，除了Al为弱变异外，其他16种元素均为中等变异程度。变异系数最高的为Mo，达46%。表明博斯腾湖流域绿洲土壤样品中各元素在空间分布上存在一定程度的差异但是并不显著。与新疆博尔塔拉河流域^[14]相比，土壤元素含量相近，空间分布特征类似。

利用SPSS 22.0对元素含量数据进行正态分布检验(k-s检验)，对于符合正态分布的含量数据[P(k-s)＞0.05]直接进行地统计学分析，服从正态分布的元素有Al、Ca、Fe、V、Zn、Cr、Cu、As、Cd、Sb、Tl、Pb、Hg，将不符合正态分布的Mn、Co、Mo经过Minitab 17软件进行Johnson变换后P值均大于0.05，服从正态分布^[28]。Ni经变换后仍然不符合正态分布，地统计分析的结果仅供参考。

通过对17种元素进行半方差函数的曲线拟合，确定了最优的理论模型和拟合曲线，半方差模型及其参数值如表 3所示，部分元素变异函数模型如图 2所示。元素Al、Ca、Fe、V、Zn、Hg的拟合结果均符合高斯模型。球状模型表明了元素的空间聚集分布程度，其空间结构是当采样间距小于变程时样点间的空间相关性随间距增大而降低^[29]。符合该分布的元素有Mn、Ni、Sb、Tl。如果半方差函数模型是非水平直线型，基台值是渐近线的则为指数模型，符合该模型的元素有Cr、Co、Cu、As、Cd、Mo、Pb。

表 3 表土元素半方差函数最优模型及拟合参数(n=77) Table 3 Semivariance model and fitting parameters of elements in topsoil(n=77)

表选项

图 2 土壤元素Al、Ca、Zn、As、Sb、Pb变异函数模型 Figure 2 Semivariograms for element Al, Ca, Zn, As, Sb and Pb of topsoil

图选项

当距离(h)为0时，半方差函数值γ(0)称为块金值(C₀)，主要由测定误差和小于最小采样尺度的非连续性变异引起，属于随机性的变异因素；由结构性因素引起的变异函数值变化称为结构方差，用C表示。基台值(C₀+C)是半方差函数从初始的C₀随h增大而达到一个相对恒定的值，表示系统内的总变异^[30]。从表 3可看出，17种元素的块金值都小于1，说明本研究的土壤元素含量采样方案和采样间距较好地反映了研究区元素的变异程度。

土壤元素的空间分布是由区域的结构性因素(如成土母质、地形、气候等)和随机性因素(如工农业生产、交通等)共同作用的结果，结构性因素会使土壤元素的空间相关性增强，而随机性因素则会使其空间相关性减弱^[31]。块金值和基台值之比是反映区域化变量空间异质性程度的重要指标，该比值用以反映空间变异影响因素中区域结构性因素(成土母质、地形、气候等)和随机性因素(人为活动等)的作用，称为块金效应^[32]。其比值越小，空间相关性越强。李哈滨等^[30]研究认为，当比值小于0.25，以结构性变异为主，具有强烈的空间自相关性；0.25~0.75为中等强度的空间相关性；大于0.75为弱的空间相关性，变量的空间变异以随机性因素为主；如果接近于1，变量在研究尺度上具有恒定的变异。研究区16种元素的块金效应均小于0.75，其中Fe、V、Zn、Co、Cu、As、Mo、Cd、Tl的块金效应小于0.25，具有较强的空间相关性；Al、Ca、Sb、Pb、Hg的块金效应大于0.3，具有中等强度的空间相关性。变程是表示变量空间相关性范围的变量，变程以内的元素具有空间自相关性，变程以外的则不存在空间自相关。Fe、Mn、V、Cr、Co、Zn、As、Tl、Hg的变程均大于50 km，在较大范围内存在相关性。说明Fe、Mn、V、Zn、Cr、Co、As、Tl、Hg受结构性因素影响较大，并且其空间自相关性良好，空间分布未受到人为因素影响。Cu、Mo、Cd、Sb、Pb变程较小，仅在小范围内存在空间相关性，说明整体结构还在，受区域随机因素(农业灌溉、交通、工业生产)影响较大，影响因素复杂。Fe的变程为165.58 km，变程较大，是地壳中稳定元素，是自然风化产物，受区域地质作用影响^[33]。Al的变程为18.01 km，变程较小，除了与成土母质有关，还受随机因素的影响。Ca的变程为59.76 km，是活动性元素，易在土壤中发生迁移，主要受成土母质的影响^{[29, 34]}。

为分析重金属元素含量的分布特征，对Mn、Zn、V、Cr、Co、Ni、Cu、Tl、Mo、Sb、As、Cd、Pb、Hg进行空间分析，结果表明(图 3)，Mn、Mo、Tl的空间分布较为均匀，区域之间的差异并不明显，仅在流域北部和硕县境内含量稍高于其他地区，自然来源为主。Zn、V、Cr、Co、Cu的空间分布特征类似，在东部和北部区域县城周边、交通干线附近含量较高。Sb、Hg、Cd、Pb的空间分布特征类似，呈斑块状分布，集中在和硕县、焉耆县和博湖县周边局部地区，和局部工农业活动有关。这说明交通尾气排放、城市和工厂固体垃圾堆放以及农药化肥过量使用等造成区域重金属含量较高^{[5, 15]}。As的高值区分布在西部开都河出山口和东部县城、街道周边，说明As既有自然来源，也有人为来源。研究表明^[35-36]，天山山脉富含砷矿，在大气搬运、降水冲刷、淋滤等作用下尘埃中与溶解出的As和其他重金属元素被河流搬运沉积于博斯腾湖流域，可能是造成区域As等重金属含量较高的重要原因。

图 3 博斯腾湖流域绿洲表土重金属元素空间分布图 Figure 3 Spatial distribution of heavy metals in oasis topsoil of the Bosten Lake Basin

图选项

综上所述，重金属含量在县城周边、省国道交通干线附近以及局部的工厂附近具有较高值，东部地区高于西部地区，东部的工农业、城镇等较西部密集，人类活动的影响程度较高。

为分析不同环境介质中重金属的危害，众多科学家提出了不同的评价方法，本文运用单因子和内梅罗污染指数法以及潜在生态风险指数法评价博斯腾湖流域绿洲表层土壤中重金属的污染程度，综合两种评价结果，以准确把握流域绿洲土壤的质量状况，为生态环境保护提供参考。

单因子指数结果显示(表 4)，Co的单因子指数小于1，没有污染，Mn、V、Zn、Cr、Ni、Cu、As、Cd、Sb、Tl、Pb、Hg的单因子污染指数均值在0~2之间，个别样点存在污染，其中被污染的样点数排序为Cd(58)＞Hg(46)＞Cr(36)＞V(31)＞Sb(28)＞As(21)＞Zn(17)＞Ni(16)＞Tl(6)＞Pb(5)＞Mn(3)＞Cu(2)。根据单因子污染指数结果计算内梅罗综合污染指数值，并评价其污染程度。结果表明Co属于安全级别，Mn、Cu、Tl、Pb的污染程度在警戒限内。V、Zn、Cr、Ni、As、Cd、Sb、Hg的内梅罗污染指数在1~2之间，属于轻微污染程度，其中Cd、Hg的污染指数值接近1.5，轻微污染程度稍高于其余重金属。

表 4 博斯腾湖流域绿洲表土重金属内梅罗污染指数 Table 4 Nemerow pollution index of heavy metal in oasis topsoil of the Bosten Lake Basin

表选项

Cd、Hg的内梅罗污染指数呈斑块状分布(图 4)，高值区均分布在流域东南部，集中在焉耆县和博湖县县城周边、交通线附近，说明其多分布在受人类活动影响强度高的区域。

图 4 Cd、Hg的内梅罗污染指数空间分布 Figure 4 Spatial distribution of Nemerow index of Cd and Hg element

图选项

利用潜在生态风险指数法计算的结果如表 5所示，博斯腾湖流域绿洲表土11种重金属平均E_rⁱ的大小排序为Hg(45.39)＞Cd(36.97)＞As(8.60)＞Pb(4.36)＞Ni(4.35)＞Cu(3.38)＞Co(2.82)＞Cr(2.02)＞V(1.95)＞Zn(0.83)＞Mn(0.81)，Cd的平均E_rⁱ值为36.97，约37.7%的样点E_rⁱ值大于40，存在轻微生态风险，Hg的平均E_rⁱ值为45.39，最大值达到74.7，约62.3%的样点风险值大于40，存在轻微生态风险，其余重金属元素的E_rⁱ值均小于40，无生态风险。综合潜在生态风险指数RI值范围为55.1~161.85，按照综合风险程度标准，27%的样点RI值小于100，为轻微风险；73%样点的RI值介于100~250之间，为中等风险。总体上看，博斯腾湖流域表层土壤重金属Cd、Hg的风险程度较高，应该引起相关部门的重视，其他重金属并无较大程度的生态风险，由于该流域工农业生产等人类活动的影响程度小，经济欠发达，因此并未对流域土壤等生态环境造成较大风险。

表 5 博斯腾湖流域绿洲表土重金属元素潜在生态风险评价 Table 5 The potential ecological risk assement of heavy metal in oasis topsoil of Bosten Lake Basin

表选项

从单因子Hg(Cd不满足插值要求)和综合潜在生态风险空间分布特征来看(图 5)，东南部明显高于西北部，风险高值区域分布在焉耆县和博湖县周边、农田分布密集区域、工厂等周围。

图 5 Hg的潜在生态危害指数和综合生态危害指数的空间分布 Figure 5 Spatial distribution of potential ecological risk index of Hg element and integration

图选项

(1) 利用经典统计学方法分析博斯腾湖流域表土中元素含量特征，结果显示17种元素平均含量顺序为Ca＞Al＞Fe＞Mn＞V＞Zn＞Cr＞Ni＞Cu＞Pb＞As＞Co＞Sb＞Mo＞Tl＞Cd＞Hg。其中常量元素Al、Ca、Fe含量远超新疆土壤背景值，除重金属元素Co、Mo含量低于背景值外，其余12种重金属元素存在不同程度的超标，Cd、Hg超标率高达77%、62%。

(2) 利用地统计方法分析了土壤元素的空间变异特征、空间分布特征。结果表明Fe、Mn、V、Cr、Co、Zn、As、Tl、Hg受结构性因素影响较大，其空间自相关性良好；Cu、Mo、Cd、Sb、Pb变程较小，仅在小范围内存在空间相关性，说明整体结构还在，受区域随机因素(农业灌溉、交通、工业生产等)影响较大，影响因素复杂。Al、Fe是地壳中稳定元素，受区域地质的影响，Al变程较小，说明还受随机因素影响。Ca是活动性元素，受土壤化学淋溶、风化影响为主。空间分布结果为大多数重金属元素在县城周边、工厂以及交通干线附近具有较高值，As在西部农业稀疏区和东南部含量较高，为自然和人为双重来源。

(3) 分别运用2种风险评价方法分析重金属的潜在生态风险；内梅罗污染指数法结果显示，Co属于安全级别，Mn、Cu、Tl、Pb的污染程度在警戒限内。V、Zn、Cr、Ni、As、Cd、Sb、Hg的内梅罗污染指数在1~2之间，属于轻微污染程度，其中Cd、Hg的污染指数值接近1.5，轻微污染程度稍高于其余重金属；单因子生态风险结果表明Cd、Hg具有轻微生态风险，综合生态风险结果为27%样点存在轻微风险，73%样点存在中等风险。风险值的空间分布表现为东南部区域较高，集中于焉耆县、博湖县县城周边城镇、交通干线附近等。总体看来，博斯腾湖流域的Cd、Hg风险程度较大，需要引起重视，其他重金属风险较低，人类工农业生产活动并未对流域土壤产生较大危害。