2015年我国规模化猪场（年出栏500头以上）生猪出栏占比44%，《全国生猪生产发展规划（2016—2020年）》指出，至2020年，全国目标出栏500头以上生猪的规模化养殖场比重将达到52%。生猪养殖的规模化和集约化发展产生了大量的猪场粪污，而土地利用是实现粪污资源化的有效措施，但不可忽略由此带来的土壤重金属污染风险^[1]。为促进动物生长和防治疾病，重金属通常被高剂量地添加到猪饲料中，以致于猪饲料中存在不同程度的重金属超标现象，超标元素包括Cu、Zn、Cd、Pb、Cr、As等^[2-3]。而猪体对重金属的利用率低，绝大部分残留在粪便中。如潘寻等^[2]发现猪粪中能检测出饲料中所含有的所有元素，其中Cu、Zn的平均含量较高，分别达472.8、1 908.6 mg·kg^-1，均超出《农用污泥中污染物控制标准》的规定值。另外，沼液、沼渣和废水中也有重金属残留。如朱泉雯^[3]检测发现除猪粪中Cu、Zn、As、Cd、Pb超标严重外，沼液中Cu、Zn、As以及沼渣中Cu、Zn、Cd、Pb、Ni也严重超标。章杰等^[4]发现不同养殖模式下的猪场废水中，重金属污染物以Cu、Zn为主，同时含有一定量的Cd、Pb元素，4种重金属的浓度均超出《农田灌溉水质标准》的限值。国外学者研究也证实了养猪沼液和废水中存在一定量的重金属，且以Cu、Zn为主^[5-6]。

综上所述，饲料中的重金属添加剂致使猪场粪污中含有大量的重金属，粪污灌溉可能使这些重金属在环境中迁移、转化后累积于土壤中，直接或间接危害人体健康^[7]。饲料重金属添加剂的种类和土壤类型有地区差异性，为了了解重庆市长期灌溉猪场粪污的土壤是否存在重金属污染状况，有必要基于某一实际猪场对其灌溉区土壤重金属进行分析评价。鉴于重金属的总量是评价土壤污染程度的重要指标，而重金属的形态及比例是决定其在环境中的迁移性和毒性的关键因素^[8]。本文选取粪污中存在量较多的Cu、Zn及毒性较强的Cd、Pb作为目标元素，依托重庆市内在规模上、农业上发展种养模式结合较为典型的一家种猪场，通过采样分析猪场粪污、废水灌溉区的土壤重金属的总量及形态，评价了猪场粪污、废水灌溉区土壤重金属污染程度和污染风险，为确定合理的灌溉体系和防治土壤重金属污染提供科学依据。

1 材料与方法 1.1 研究区域概况

研究区位于重庆市潼南区花岩镇，该种猪场有完善的粪污处理系统，猪粪、猪尿干湿分离后，猪粪、粪渣外运至有机肥厂制成有机肥，尿液、废水经格栅、集水池、固液分离、沉淀池、酸化调节、UASB反应器、二级A/O、除磷沉淀、pH回调池、氧化塘处理后农用，周边配有500 hm²农田来消纳废水。猪场建立前，该地区以养殖散户为主，畜禽粪便几乎全部用作基肥。猪场建立后的近4 a来，每日的排粪量21 t，废水量600 m³，该地区大力发展种养模式结合，以有机肥为主，无机肥为辅，猪场的粪便和废水是周边土壤肥料的主要来源。

1.2 土壤样品采集与测试

土壤样品采样点遍及猪场内及场外，共设置16个，编号为S1~S16。每个采样点视情况采用不同的采样方法，如对污灌水田采用对角线法，对面积小、土壤均匀的地块采用梅花点法，对面积中等、土壤不均匀的地块采用棋盘式法，对坡地采用蛇形法，采集3~ 5个点的表层（0~20 cm）土壤混合为一个代表样。样品剔除石块、沙砾和植物残体等，风干后磨碎，过1 mm尼龙筛用于土壤pH的测定，过0.149 mm筛用于测定土壤重金属总量及分形态含量。

土壤pH采用土液比1:2.5进行配比，pH计法测定。土壤重金属全量的测定采用HF-HNO₃-HClO₄三酸消解法^[9]。土壤重金属形态分析采用BCR三步浸提法^[10]，将土壤分为酸溶态、可还原态、可氧化态。残渣态在BCR的基础上，采用混酸消解方法提取^[11]。

1.3 土壤重金属污染评价方法 1.3.1 单因子污染指数法

该法利用实测数据和标准对比，能评价区域单一重金属的污染程度，计算公式为：

式中：P_i为土壤中重金属i的污染指数；C_i为实测含量；S_i为评价标准，本文以国家土壤环境质量标准（GB 15618—1995）二级标准（以下简称“土壤质量二级标准”）的规定值计算。评判标准：P_i≤1时，未受污染；1 < P_i≤2，轻度污染；2 < P_i≤3，中度污染；P_i>3，重度污染。

1.3.2 内梅罗污染指数法

该法可用来评价重金属的平均污染水平，并突出污染最严重的重金属。计算公式为：

式中：P_N为内梅罗污染指数；P_imax为土壤中重金属的最大污染分指数；P_i为各污染分指数的算术平均值。P_N≤0.7，清洁；0.7 < P_N≤1.0，尚清洁；1.0 < P_N≤2.0，轻度污染；2.0 < P_N≤3.0，中度污染；P_N>3.0，重度污染。

1.3.3 地累积指数法

该法通常称为Muller指数，不仅能反映重金属分布的自然变化特征，而且可以判别人为活动对环境的影响，是区分人为活动影响的重要参数。计算公式为：

式中：I_geo为地累积指数；B_i为重金属i的地球化学背景值，本文以鲍丽然等^[12]调查的重庆渝西经济区土壤元素（指标）地球化学参数的背景值作为重金属的区域背景值；k为考虑岩层差异所引起背景值变化的调整系数，该值取1.5。I_geo分级标准见表 1。

1.3.4 土壤重金属潜在生态污染评价

瑞典学者Hakanson提出的潜在生态污染指数法可反映某一特定环境中各种污染物的影响或多种污染物的综合影响，综合考虑了多元素协同作用、毒性水平、污染浓度以及环境对重金属污染敏感性等因素，消除了区域差异性及异源污染的影响，在土壤重金属风险评价中得到广泛应用^[13]。Hakanson潜在生态风险指数计算公式为：

式中：C_rⁱ为土壤中重金属i相对参比值的污染参数；C_nⁱ为参比值，本文以土壤质量二级标准和区域背景值作参比值；E_rⁱ为重金属i的潜在生态风险参数；T_rⁱ为重金属i的毒性响应系数，参考Hakanson^[14]制定的标准化重金属毒性响应系数，Cu、Zn、Cd、Pb分别为5、1、30、5；RI为土壤重金属综合潜在生态风险指数。重金属潜在生态危害指标与分级标准见表 2^[13]。

1.3.5 土壤重金属生物有效性

重金属的酸溶态（T₁）容易被植物吸收，可还原态（T₂）和可氧化态（T₃）不易被植物吸收，残渣态（T₄）几乎不被植物吸收。根据植物利用的难易程度，计算生物活性系数Ki值，将重金属分为可利用态[K₁=（T₁+T₂）/（T₁+T₂+T₃+T₄）]、中等利用态[K₂=T₃/（T₁+T₂+T₃+T₄）]和难利用态[K₃=T₄/（T₁+T₂+T₃+T₄）]^[15]，从而评价重金属的生物有效性。

1.4 数据处理

实验数据采用Excel 2016和Origin Pro 8.5软件进行分析和处理。

2 结果与讨论 2.1 土壤重金属总量分布特征

根据《土地利用现状分类标准》（GB/T 21010—2017）二级类，将采样点分为不同的土地利用类型。鉴于采样点均是采用猪场粪污灌溉，本文不严格区分采样点来自猪场内还是场外。16个典型土壤样品重金属含量及pH值见表 3。

统计发现，灌溉区Cu、Zn、Cd、Pb的平均含量分别为30.80、91.33、0.55、29.34 mg·kg^-1，土壤平均pH值为8.36，总体呈弱碱性。与土壤质量二级标准相比，各采样点中只有Cd有4个超标点。而Cu、Zn、Cd、Pb的平均含量均超过了区域背景值，超标率分别为26.23%、8.73%、150%、4.79%，超标点分别有12、10、10、10个，其中Cu的超标点最多。

变异系数反映离散程度，该值越大，表示土壤重金属含量受人类活动干扰越强烈。变异系数 < 0.1为弱变异，0.10~0.30为中等变异，>0.30为强变异^[16]。各元素的变异系数Cd（1.06）>Cu（0.31）>Zn（0.23）>Pb（0.14），Cd、Cu属于强变异性水平，Zn、Pb属于中等变异性水平。

调研发现，猪场猪粪中Cu、Zn、Cd、Pb的含量分别为21.94~1 202.36、200~1 512.67、0.45~1.37、2.46~3.98 mg·kg^-1。处理达标的废水中4种金属的含量分别为0.021、0.014 8、0.001 6、0.000 7 mg·L^-1，推测猪粪的施用可能是造成土壤Cu、Zn累积以及离散的主要原因，Cd和Pb的累积和离散可能是由其他人为活动引起的。

2.2 土壤重金属污染评价 2.2.1 单因子污染指数法评价结果

土壤重金属元素单因子污染指数计算结果见表 4。单因子污染指数平均值Cd>Cu≈Zn>Pb。只有Cd属于有污染状态，其中S11、S15属于轻度污染，S4、S7、S12属于中度污染。

2.2.2 内梅罗污染指数法评价结果

采用内梅罗综合污染指数法计算猪场灌区土壤重金属各元素的污染指数结果如下所示：

由计算结果可知，灌溉区土壤Cd属于中度污染，土壤几乎不受Cu、Zn、Pb的污染，土壤环境相对安全，4种金属对土壤的污染风险性为Cd>Cu>Zn>Pb。

2.2.3 地累积指数法评价结果

土壤重金属地累积指数污染等级结果见表 5。Cu有4个点（S3、S14、S8、S12），Zn有2个点（S3、S6），Cd有2个点（S5、S10）可能有污染风险。另外，Cd有2个点（S11、S15）属于中风险状态，3个点（S4、S7、S12）属于中等偏强风险状态。Pb没有污染风险。

综上，采用不同的污染评价方法对土壤重金属总量的评价结果有共性。总体来说，猪场灌区土壤Cd污染风险最大，个别采样点已存在中度污染。鉴于养猪粪污中Cd的含量少，实际农用过程中需要找寻原因以降低Cd风险。土壤虽未受到Cu、Zn污染，但I_geo显示出极个别点介于从无到有的污染等级，鉴于猪场粪污中高含量Cu、Zn的存在，农用时需要密切关注这两种重金属在土壤中的含量变化，防止污染的发生。另外，土壤不存在Pb污染。

2.3 土壤重金属分形态含量分布特征

重金属形态含量百分比是指重金属各形态含量占总量的比例，更准确地指示环境污染对土壤的冲击^[17]。各形态的稳定性有差异，T₁为易变型，T₂、T₃为可变型，T₄为不可变型。在未受污染的自然土壤中，T₁所占比例很低，T₄所占比例较高，在污染的土壤中T₄含量较低，外源重金属会转化为具有活性的形态^[18]。猪场粪污灌区土壤样品重金属各形态质量比数据见表 6，各采样点重金属分形态含量及占比情况见图 1。

灌区内土壤Cu主要以残渣态存在，非残渣态占比约为15%。说明Cu易与土壤中的有机物和矿物质形成难分解的有机络合物和硫化铜等难分解矿物^[19]，但Cu的各形态均处于强变异水平，人为活动对猪场Cu形态分布的影响大。

Zn主要以残渣态存在，但非残渣态的比例较大，约占总量的43%，且处于强变异水平，此部分Zn容易释放出来造成二次污染，因此Zn存在较大的环境风险。

Cd的残渣态占比小，非残渣态比例几近88%。非残渣态中，可变型还原态比例高达53%，接近刘艳萍等^[9]研究养殖废水灌区内土壤重金属Cd形态分布时得出的结论。总的来说，Cd的非残渣态比例大，离散程度也大，表现出了较强的环境风险。

Pb主要以稳定残渣态存在，这可能是因为Pb容易与土壤中的铁锰氧化物结合成螯合物或络合物^[20]。非残渣态约占总量的32%，但变异系数较大，说明Pb虽然在形态上表现出环境风险小，但其形态分布在一定程度上受人类活动的影响。

图 1显示不同重金属元素在同一个采样点中的化学形态分布不同，同一种元素在不同采样点中的形态分布也具有差异性。总体来说，不同元素形态平均变异系数以Cd最大，Cu、Zn次之，Pb最小。Cu、Zn的变异系数由大到小表现为酸溶态、可还原态、可氧化态、残渣态，这与钟晓兰等^[21]的研究结果一致。Cd的变异系数表现为可还原态、酸溶态、残渣态、可氧化态。Pb的变异系数表现为可还原态、可氧化态、酸溶态、残渣态。总体来说，非残渣态较残渣态的变异系数大，说明重金属次生相态含量受外界干扰比较显著，具有较强的空间分异，这种空间分异是耕作情况、管理措施、种植制度、污染源引入等人为活动综合作用的结果。

2.4 土壤重金属生物有效性和潜在生态污染评价 2.4.1 土壤重金属生物有效性

依据表 2数据计算各重金属的生物有效性结果如表 7所示。Cu、Zn和Pb的生物活性系数均表现为K₃>K₁>K₂，表明它们不易被外界生物利用，对生物体造成影响的能力是有限的。但值得注意的是Zn的K₁高达0.4，鉴于Zn在土壤中的总量高，实际中不可忽视Zn的潜在危害。Cd的生物活性系数变化规律为K₁ >K₂>K₃，K₁远大于K₂和K₃，说明土壤中的Cd容易进入生物体，虽然灌区Cd总量小，但形态引起的环境风险高，应该引起足够的重视。4种重金属的K₁值由大到小为Cd>Zn>Pb>Cu，Cu和Pb有效性较低，环境风险小。

2.4.2 土壤重金属潜在生态污染评价

酸溶态、可还原态及可氧化态在环境中的稳定性差，本文以这三者之和为依据做出的生态环境评价结果见表 8。以土壤质量二级标准为参比值，各重金属的E_rⁱ平均值均小于40，单项潜在生态风险均为轻微等级，生态危害程度较低。以区域背景值为参比值时，唯有Cd的平均值介于40~80之间，其单项生态潜在风险达到中等水平。两种参比值下，RI分别为26.74和74.56，均小于80，重金属的综合潜在生态风险为轻微等级，表明采样区整体的重金属污染风险低。参比值的选择影响重金属对综合生态风险指数的贡献程度，参比土壤质量二级标准，贡献率Cd>Cu>Zn>Pb，而参比区域背景值，贡献率Cd>Pb>Cu>Zn，但两种参比值下均以Cd的贡献率最大，说明采样区需要注意Cd的污染防治。

3 结论

（1）养猪场灌溉区土壤呈弱碱性，以重金属的总量来看，Cu、Zn、Cd、Pb的均值均超过了区域背景值，变异系数Cd>Cu>Zn>Pb，Cd和Cu已达强变异水平。

（2）采用单因子污染指数法、内梅罗污染指数法、地累积指数法等污染评价方法对土壤重金属进行污染评价得出，土壤个别样点存在中度的Cd污染，但不存在Cu、Zn、Pb污染。

（3）重金属种类和采样点的不同使土壤重金属形态分布具有差异性，土壤中Cu、Zn、Pb均主要以残渣态存在，Cd主要以可还原态存在。四种重金属的形态平均变异系数以Cd最大，Cu、Zn次之，Pb最小。总体来说，残渣态的变异系数较非残渣态的小。

（4）生物活性系数Cd>Zn>Pb>Cu，Cd最易被外界生物利用。以非残渣态之和为依据，参比土壤质量二级标准，各重金属的潜在生态风险均为轻微级，参比区域背景值，Cd的生态风险为中等水平，其余元素均为轻微级。4种重金属的综合潜在生态风险为轻微等级。