杨杰, 董静, 宋洲,2, 杨成梅, 刘田, 周顺超, 胡核, 黄聪
(1.湖北省地质实验测试中心, 湖北 武汉 430034;2.资源与生态环境地质湖北省重点实验室(湖北省地质局), 湖北 武汉 430034;3.湖北华祥地质环境检测科技有限公司, 湖北 武汉 430000)
中国铅锌矿产资源较为丰富[1],但铅锌矿床主要为中小型,且矿石组分复杂,通常伴生Cu、Cd、As等重金属元素[2],矿业活动过程中这些重金属元素通过地表径流、风力传送和雨水淋滤等方式扩散到周边河流、土壤,继而下渗到地下水,造成周边土壤重金属污染[3-5],甚至引发粮食安全问题[6-8],因此客观评价尾矿库周边农田土壤和农作物污染状况对后期尾矿库监管、土壤污染防治和安全利用有着重要的指导意义。
目前,土壤重金属污染评价方法有很多,按照评价对象不同,分为土壤单独评价法和土壤农-作物系统复合评价法两大类,两者区别在于土壤-农作物复合评价法考虑土壤作为农作物赖以生存的环境要素,将土壤和农作物有机联系起来综合评价,而前者仅对土壤自身评价。土壤单独评价法发展成熟,主要有单因子污染指数法[9-10]、内梅罗污染指数法[11]、潜在生态风险指数法[12-13]、地累积指数法[14-15]、富集因子法[16]、模糊数学法[17-18]、灰色聚类法[19-20]、层次分析法[21-22]等,其中,应用最广泛的是单因子污染指数法、地累积指数法以及潜在生态风险指数法。土壤-农作物复合评价法发展相对较晚,但随着《土壤环境质量 农用地土壤污染风险管控标准(试行)》(GB 15618—2018)正式实施,该标准以确保农作物质量安全为主要目标,农用地土壤污染评价由单一的土壤监测转变成土壤和农产品协同监测,复合评价法很有必要。目前土壤-农作物复合评价法应用最广泛的有王玉军等[23]提出的综合质量影响指数法[24-26]和美国USEPA开发的健康风险评价法[27-30]。
位于湖北西部村庄的某铜铅锌矿主要矿物成分有闪锌矿、方铅矿、黄铜矿、斑铜矿、黄铁矿及白铁矿,矿区开采时间长,且由于历史原因污染防治设施落后,铜铅锌尾矿成分及镉、砷等伴生元素可能会排放到环境中,对周围居民生活、农田种植带来风险,迄今为止该尾矿库周边农田土壤重金属污染现状和对人体的健康风险尚不明晰。基于此,本文选取尾矿库周边农田土壤和水稻为研究对象,分析土壤和水稻中8项重金属的污染状况,并采用潜在生态风险指数法对周边农田土壤重金属污染程度进行评价,通过适合中国公民的人体健康风险评价模型评估土壤和水稻重金属的人体健康风险,综合评价了尾矿库对周围环境造成的影响,为该尾矿库周边农田生态环境保护、生态农业建设和可持续发展提供科学依据。
该铜铅锌尾矿库位于湖北省西部村庄,地处沮漳河中下游,大巴山脉东麓,荆山山脉以南,是鄂西山地向江汉平原的过渡地带。所在地区属丘陵地貌,以南北向低缓丘陵及冲沟为主。区内地势北高南低,海拔2600~4500m,高差悬殊,地形坡度较陡,多悬崖绝壁。区内无常流性地表水体,地表水系不发育,多为小型冲沟、堰塘及小水库等,主要为大气降水补给,通过地表径流由高到低向堰塘、水库及冲沟排泄。区内属亚热带季风性湿润气候,四季分明,雨热同季,风向季节性变化较大,冬季多偏北风,夏季多偏南风,春秋两季则偏东风。研究区所在乡镇常驻人口约6万人,区内村民居住分散。区内地表植被发育多为松树及灌木。土地利用以林地、草地、园地和耕地等农业用地为主。主要农作物为水稻、玉米、小麦和蔬菜等。研究区铜铅锌矿尾矿库建成于1986年,2015年关闭后成为无主尾矿库,库区堆积总库容约9.5万方,宽50~140m,长约400m,属于山谷型尾矿库,库区渗滤液收集后汇入南面下游的渗滤液收集池,后进入地表排水系统。
根据研究区农田分布情况,总体按100m×100m网格布设监测点,尾矿库下游南面和紧邻的东面500m范围内农田加密布点,按50m×50m网格设置,共计布设50个监测点,同步采集土壤样品和水稻稻穗样品,采集表层土壤(深度为0~20cm)和水稻稻穗样品各55组(各含10%的平行样)。采样点位分布如图1。
土壤样品采样按照《土壤环境监测技术规范》(HJ/T 166—2004)进行,每个土壤监测点土壤样品为该网格范围内随机5个点位的混合样,表层土壤采样深度为0~20cm。采样时去除土壤样品中的砾石、动植物残渣以及肥料团块等杂物,每件样品大于2kg。同时现场用GPS记录采样点位置信息,填写采样记录。水稻稻穗采样参照《农、畜、水产品污染监测技术规范》(NY/T 398—2000)执行,于收获期采集水稻籽实,采样点位与对应的土壤点位保持一致,在同一田块不同位置采集的稻穗等量混匀组成一个混合样品(1kg)。
图1 研究区及采样点位分布图Fig.1 Map of sampling points distribution and the location of the study area
土壤样品经室温风干后研磨,过100目标准筛,再消解后上机测试。电感耦合等离子体质谱仪(X2型,美国ThermoFisher公司)测定Cd、Pb、Cu、Ni、Cr、Zn,检测依据《区域地球化学样品分析方法 第5部分:镉含量测定 电感耦合等离子体质谱法》(DZ/T 0279.5—2016)、《区域地球化学样品分析方法 第3部分:钡、铍、铋等15个元素量测定 电感耦合等离子体质谱法》(DZ/T 0279.3—2016)、《区域地球化学样品分析方法 第2部分:氧化钙等27个成分量测定 电感耦合等离子体原子发射光谱法》(DZ/T 0279.2—2016),检出限分别为0.021、0.5、0.6、0.6、0.2、0.03mg/kg;原子荧光光谱法(AFS-8500型,北京海光)测定As和Hg,检测依据《土壤和沉积物汞、砷、硒、铋、锑的测定 微波消解/原子荧光法》(HJ 680—2013)、《土壤质量 总汞、总砷、总铅的测定 原子荧光法 第1部分:土壤中总汞的测定》(GB/T 22105.1—2008),检出限分别为0.01mg/kg、0.002mg/kg;酸度计(PHS-3C型,上海雷磁)测定土壤样品pH值,检测依据《土壤 pH值的测定 电位法》(HJ 962—2018)。
水稻稻穗经脱粒后于70℃恒温烘箱烘干,脱壳粉碎后消解,再上机测试。电感耦合等离子体质谱仪(X2型,美国热电)测定Cd、Pb、Cu、Ni、Cr、Zn,检测依据《食品安全国家标准 食品中多元素的测定》(GB 5009.268—2016),检出限分别为0.002、0.02、0.05、0.2、0.05、0.5mg/kg;原子荧光光谱法(AFS-3100,北京海光)测定As和Hg,检测依据《食品安全国家标准 食品中总砷及无机砷的测定 第一篇 总砷的测定 第二法 氢化物发生原子荧光光谱法》(GB 5009.11—2014)、《食品安全国家标准 食品中总汞及有机汞的测定 第一篇 食品中总汞的测定 第一法 原子荧光光谱分析法》(GB 5009.17—2014),检出限分别为0.003、0.0005mg/kg。
所有样品检测均按照规范要求进行全过程质量控制,每批次样品插入10%比例的空白样、平行样(水稻样品比例为100%)以及国家一级有证标准物质试验,平行试验相对误差小于20%,有证标准物质均在推荐值范围内,数据质量可靠。全部样品测试工作在湖北省地质实验测试中心化学分析室完成,所有检测项目均通过国家CMA计量认证。
通过单项污染指数法和潜在生态风险指数法评价该铜铅锌尾矿库周边农田土壤污染状况,采用美国环境保护署(USEPA)开发的健康风险评估模型评估土壤和水稻中重金属人体健康风险。
2.3.1潜在生态风险指数法
潜在生态风险指数法是1980年瑞典科学家Hakanson从沉积学角度提出的评价重金属污染程度和潜在生态危害的一种方法[31]。该方法除了考虑重金属的含量之外,还将重金属的生态效应、环境效应和毒性联系,突出污染较严重、毒性较强的重金属的作用,采用具有可比的、等价指数分级法进行评价。其表达式为:
2.3.2重金属健康风险评价模型
本研究采用美国USEPA开发的健康风险评价模型对铜铅锌尾矿库周边土壤-水稻系统进行健康风险评价,模型计算公式具体如表1所示。
CRi为致癌重金属i的单项健康风险指数;CDIij为重金属i第j种暴露途径日均暴露量;SFij为致癌重金属i第j种暴露途径的斜率系数;TCR为致癌重金属通过3种途径的总致癌风险。当CR或TCR≤1×10-6时,认为不存在致癌风险或风险较小,风险可忽略不计;当CR或TCR>1×10-6时,认为存在较高的致癌风险。HQi为非致癌重金属的单项健康风险指数,RfDij为非致癌重金属i第j种暴露途径的参考剂量,HI为8种重金属通过3种途径的总非致癌风险。HI值低于1可认为对敏感人群不存在风险,HI值大于1可认为超过人体健康可接受阈值,不可接受。各参数含义见表2,数值见表3。各参数数值主要来源于《建设用地土壤污染风险评估技术导则》(HJ 25.3—2019)和前人研究成果[34]。
表1 健康风险评价模型计算公式
表2 健康风险评价暴露参数
表3 重金属不同暴露途径参考计量(RfD)和致癌斜率因子(SF)
2.3.3数据整理和分析
数据整理和分析采用Excel 2016和Origin 2018(OriginLab公司,美国)完成,研究区采样点分布图采用Arcgis 10.2(ERSI公司,美国)软件绘制完成。
尾矿库周边农田土壤pH值范围为4.97~8.15,其中pH≤5.5点位仅1个,5.5
表4 研究区土壤重金属含量和pH值结果统计
表5 研究区水稻重金属含量特征
对研究区土壤重金属超标点位空间分布情况分析发现,超标区主要集中分布在尾矿库下游300m范围内农田,且随着农田与尾矿库的距离增加,土壤重金属含量也逐渐降低,表明土壤重金属污染程度在逐渐减轻。地表径流、风力传送和雨水淋滤是主要扩散方式,考虑到当地风向季节性变化较大,冬季多偏北风,夏季多偏南风,而超标区集中在尾矿库南面渗滤池下游,尾矿库北面上游农田土壤重金属含量均不超标,这在一定程度上说明土壤重金属的污染可能与尾矿库渗滤液的泄漏有关。尾矿库渗滤液中含有大量重金属,在自然降雨的作用下极易通过地表径流向下游扩散,从而造成周边土壤重金属含量超标[17]。尾矿库周边农田土壤Cd、Cu、Pb、Zn、As重金属均存在环境风险,Cd元素最突出,需重点关注农产品的Cd污染风险。
研究区尾矿库周边农田水稻重金属含量结果见表5,鉴于《食品安全国家标准 食品中污染物限量》(GB 2762—2017)未给出Cu、Zn和Ni元素的标准限量,本研究仅对水稻中As、Cd、Pb、Hg、Cr元素的污染状况进行评价,50个水稻点位中有7个重金属含量超过标准限值,超标率14%,且超标元素均为Cd,最大超标倍数为9.39倍,其他元素均不超标,这可能与水稻对Cd具有较强的富集能力有关。研究区种植的水稻Cd易超标,存在粮食安全风险,在中国南方其他地区也发现水稻Cd超标较高的现象[35-36]。有研究表明与其他谷类作物相比,水稻具有更强的Cd积累能力,而作为主要粮食作物之一的水稻也是经膳食摄入Cd的主要来源,因此水稻Cd超标对人体具有较高的健康风险[37]。
在空间分布上,5个水稻Cd超标点位集中分布在尾矿库下游300m范围内,与土壤重金属高含量区分布基本一致;随着与尾矿库距离增加,水稻Cd含量降低。本研究区为丘陵地貌,以南北向低缓丘陵及冲沟为主,地势相对平缓处有稻田分布,中间间隔有地势较高的林地,由此推测降雨和地表径流对尾矿库的侵蚀和冲刷是造成下游农田土壤重金属污染的主要原因。
研究区土壤污染评价结果与农作物污染评价结果并不完全一致,可能由于土壤中重金属具有隐蔽性和滞后性等特点,如果仅以土壤中重金属的风险等级划分污染风险区和污染管控区,可能无法与农作物的污染风险区对应,说明仅依靠土壤污染风险管控标准(GB 15618—2018)开展农田土壤-农作物系统的污染评价可能存在偏差,应当结合农作物生物效应开展复合评价。
元素含量的相关性分析可以判定不同元素是否来自相同污染源,有研究表明来自相同污染源的重金属之间的相关性较大[38]。研究区50个土壤采样点重金属含量相关性分析结果表明(图2),土壤中As、Cd、Cu、Pb、Zn和Ni六种元素含量之间均存在极显著正相关关系(p<0.01),说明土壤中的As、Cd、Cu、Pb、Zn和Ni六种元素具有明显的同源性[39],由于该尾矿库为铜铅锌尾矿库,As和Cd为伴生元素,故推测土壤中As、Cd、Cu、Pb、Zn和Ni六种重金属元素均来源于尾矿库。由于当地风向季节性变化较大,冬季多偏北风,夏季多偏南风,而重金属超标点位几乎全部分布在尾矿库下游,尾矿库上游农田土壤重金属含量均不超标,可初步推测尾矿库扬尘不是主要污染途径,是尾矿库渗滤液泄漏,通过地表径流的途径扩散。而土壤中的Hg和Cr元素与其他几种重金属元素之间相关性较弱,且均与湖北省土壤元素背景值相近[32],推测这两种元素来源于土壤母质。
对水稻与土壤中重金属含量进行相关性分析的结果表明(图2),水稻与土壤中对应的As、Cd、Cu、Pb、Zn、Ni和Cr七种元素含量均具有一定的正相关关系,且水稻中的Cd、Cu、Zn三种元素具有显著正相关关系,说明水稻富集的重金属元素主要来源于土壤,并且与其他元素相比,水稻对土壤Cd、Cu、Zn三种元素具有较高的富集能力,导致水稻中这三种元素的含量相对较高。
利用潜在生态风险指数法评价研究区农田土壤的生态风险,从评价结果(表6)可知,8项重金属平均富集系数由大到小依次为:Cd(4.41)>Pb(2.55)>Zn(2.17)>Cu(1.50)>As=1.00>Ni(0.91)>Cr(0.87)>Hg(0.72),表明Cd、Pb、Zn、Cu四个元素显著富集,这与尾矿库主要污染元素相符,进一步说明尾矿库周边农田土壤重金属主要来源于尾矿库。尾矿库在降雨、地表径流以及风化等自然因素作用下会导致重金属等污染物向周围环境中不断扩散,从而增加尾矿库周边土壤中重金属的潜在生态风险[40]。
表6 农田土壤不同重金属富集系数和潜在生态风险(Ei)评价结果
土壤中8项重金属的单项潜在生态风险平均值由大到小依次为:Cd(132)>Hg(28.7)>Pb(12.8)>As(10.0)>Cu(7.52)>Ni(4.54)>Zn(2.17)>Cr(1.74),其中Cd单项潜在生态风险指数(RI)均值最大,达强风险等级(80 综合潜在生态风险RI平均值为200,处于中度风险。整体上研究区农田土壤78%采样点重金属潜在生态风险为轻微,14%点位处于中度生态风险,2%点位处于强生态风险和很强生态风险,4%点位处于极强生态风险(图3)。潜在生态风险主要集中在尾矿库下游农田。因此,尾矿库生态环境监测需重点关注Cd元素的生态风险,以防止Cd污染进一步扩散。 图3 农田土壤重金属潜在生态风险等级分布Fig.3 Class distribution of ecological risk of heavy metals in farmland soil 利用健康风险评价模型对尾矿库周边农田土壤-水稻系统进行评价,结果见表7。单一重金属不同暴露途径的致癌风险CR平均值(表7)由大到小依次均为:经口摄入>皮肤接触>呼吸吸入,且皮肤接触和呼吸吸入途径比经口摄入的风险小约10倍,即经口摄入是致癌风险的主要途径。土壤重金属总致癌风险TCR平均值为3.50×10-5大于可接受水平(10-6),全部点位土壤总致癌风险大于可接受水平,说明该区域整体存在土壤重金属致癌风险。图4a显示四种致癌重金属元素的单项TCR平均值由大到小为:As>Cd>Pb>Ni,As和Cd致癌风险贡献较大,贡献平均值分别为85.7%和11.7%,两者贡献之和占比超过95%。 对单一重金属不同暴露途径的非致癌风险,HQ值大小依次均为:经口摄入>呼吸吸入>皮肤接触,即经口摄入也是非致癌风险的主要途径。研究区土壤重金属总非致癌风险HI平均值为1.74,最大值为10.1,50个点位中有49个点位HI值大于1,说明研究区普遍存在重金属非致癌风险。由图4b可知8种重金属元素HQ平均值由大到小为:As>Pb>Cr>Ni>Cd>Cu>Zn>Hg。比较各元素非致癌风险贡献值(图5),发现As为主要贡献元素,其次是Pb。值得注意的是,As和Pb在部分点位HI值异常高,且大于1,说明位于尾矿库下游距离渗滤池较近的农田非致癌风险较高,应引起足够的重视。 表7 农田土壤和水稻中不同暴露途径重金属致癌风险和非致癌风险 图4 研究区农田土壤重金属健康风险评价结果 图5 研究区农田土壤中重金属非致癌风险贡献值Fig.5 Non-carcinogenic risk contributions of varied metals in farmland soil from the study area 根据研究结果证实,不同重金属经口摄入普遍是皮肤接触和呼吸吸入的10倍,经口摄入是主要的暴露途径,这与前人研究结果基本一致[41],因此要关注与当地土壤重金属直接相关的农作物中的重金属含量,避免同一暴露途径下土壤和农作物中重金属累计摄入,降低对人体健康的影响。本研究在评估农田土壤健康风险的基础上进一步评估了农作物摄入途径对人体的健康风险。当地居民以水稻为主要食物,本文食物暴露途径主要考虑水稻的摄入。研究区水稻重金属总致癌风险TCR平均值为3.92×10-6,略大于可接受水平(10-6),几乎全部点位TCR值大于可接受水平,说明该区域水稻重金属整体存在致癌风险。三种致癌重金属元素的单项致癌风险CR平均值由大到小为:Cd>As>Pb,Cd和As元素CR值大于可接受水平,Pb元素小于可接受水平。与土壤重金属三类摄入途径相比,水稻摄入途径的TCR值与呼吸吸入和皮肤接触途径相当,占经口摄入途径的1/10。非致癌风险方面,水稻重金属HI平均值为0.028,所有点位HI值全部小于1。与土壤重金属三类摄入途径相比,水稻摄入途径的HQ值与呼吸吸入途径相当。 本文以鄂西某铜铅锌尾矿库为研究对象,调查其周边农田土壤和水稻中8项重金属分布,采用单项指数法和潜在生态风险指数法评价其污染状况,并进一步评估了土壤-水稻系统的健康风险。结果表明:①研究区土壤As、Cd、Cu、Pb、Zn存在超标,其中Cd超标率20%最大;水稻中重金属仅Cd超过食品安全国家标准限值,超标率14%。②土壤中重金属有相同的来源,尾矿库渗滤液泄漏是可能的污染源;水稻与土壤重金属含量具有一定的正相关性,其中Cd元素相关性最好,可能由于水稻对土壤Cd吸收能力强。③研究区农田土壤Cd、Pb、Zn、Cu显著富集,Cd富集系数达4.41;研究区整体处在中度生态风险,6%点位具有极强潜在生态风险。④研究区全部点位土壤总致癌风险大于可接受水平,存在重金属致癌风险,As和Cd风险较大;几乎全部点位总非致癌风险大于可接受水平,8种重金属中有且仅有As元素的单项非致癌风险大于可接受水平,As为主要贡献元素,其次是Pb和Cr。研究区几乎全部点位水稻总致癌风险大于可接受水平,水稻中重金属整体存在致癌风险,最大贡献来自Cd元素;所有点位水稻总非致癌风险全部在可接受水平内。 综合认为,该尾矿库周边农田土壤和水稻已受到重金属污染,存在一定的生态风险,对当地居民健康造成的风险值得足够重视,尤其是对Cd元素,应优先控制。3.5 农田土壤和水稻健康风险评价
4 结论