乌兹别克斯坦阿姆河流域水体中多环芳烃的分布、来源及风险评估*

金 苗，吴敬禄，占水娥，Shakhimardan Shaniyazov

(1：中国科学院南京地理与湖泊研究所,湖泊与环境国家重点实验室，南京 210008) (2：中国科学院大学，北京 100049) (3：Berdakh Karakalpak State University, Nukus 230100，Uzbekistan)

多环芳烃(polycyclic aromatic hydrocarbons, PAHs)是具有致癌、致突变和致畸变的一类持久性有机污染物. 其高毒性、持久性以及生物富集效应对人类健康及生态环境造成极大危害, 其中16种单体多环芳烃早在1980s已被美国环境保护署(USEPA)确定为环境中优先控制污染物[1]. 近年来, 国内外学者分析了不同区域多环芳烃的分布、来源，并对当地环境、人类健康等进行了风险评价, 指出环境中多环芳烃的污染来源主要为煤炭、石油、木材等不完全燃烧、工业排放以及石油类物质泄露等. 同时, 水体中多环芳烃也可以通过地表径流、大气干湿沉降、生活及工业污水排放等方式进入环境介质并通过大气-水-土等途径在不同环境介质中迁移[2-4].

阿姆河是中亚水量最大的内陆河, 是咸海的两大水源之一, 同时也是乌兹别克斯坦境内的主要水源, 为当地提供工农业及生活用水. 阿姆河发源于帕米尔高原, 最终汇入乌兹别克斯坦境内的咸海, 在河口形成大规模的湿地三角洲. 苏联时期咸海流域地区曾大力建设运河、灌渠及水库, 水资源被大量节流用于发展工农业[5], 导致阿姆河径流量骤减, 咸海水体面积快速萎缩[6]. 另一方面, 工业和农业领域排放的污染物扩散到三角洲地区, 尤其是单一的棉花种植业, 为保持稳定的高产量逐年增加农药和化肥的使用, 致使咸海流域生态环境恶化, 严重威胁该地区居民的生活乃至生命健康[7]. 已有研究显示, 在咸海流域周边居民的血液、头发及尿液等组织中发现超标的砷及有机氯农药含量[8-11], 同时咸海流域的地表和地下水中氮及其他元素如铜、铅、铬和铀浓度也均高于世界卫生组织的限定值[12-14]. PAHs作为与人类活动最为相关的持久性有机污染物已被广泛关注, 但目前对咸海流域的PAHs污染水平, 尤其是水体污染的研究较少. 本文分析了乌兹别克斯坦境内阿姆河地区不同水体多环芳烃浓度、分布及来源, 并对其潜在风险进行了评估. 研究区域覆盖了乌兹别克斯坦境内阿姆河下游河段以及阿姆河三角洲地区, 其中下游河段区域主要位于花拉子模州境内, 该地种植业发达, 棉产量为全国最高, 也是水稻的主要产区, 而阿姆河三角洲区域则处于行政区卡拉卡尔帕克斯坦共和国内, 该区域首府努库斯及周边城市人口密集、交通条件优越, 同时也集中了化工厂、发电站等大型企业. 了解不同区域环境中PAHs的组成特征以及受人类活动及生态环境影响的因素, 可为该地区水资源及人类生态环境综合评估提供科学依据.

1 材料与方法

1.1 试剂与仪器

主要仪器：Agilent1200高效液相色谱-二极管阵列检测器串联荧光检测器(HPLC-DAD-FLD, 美国安捷伦科技有限公司)、固相萃取真空装置(美国色谱科公司)、纯水仪(MILLIPORE Elix Essential-5, 美国密里博公司)和旋转蒸发仪(BUCHI R-300, 瑞士步崎有限公司).

主要试剂及药品：德国默克色谱纯试剂：正己烷(n-Hexane)、乙腈(Acetonitrile)、甲醇(Methanol)、二氯甲烷(Dichloromethane). 16种美国环保部(USEPA)优先控制的多环芳烃(PAHs)混合标样，包括萘(naphthalene, Nap)、苊(acenaphthene, Ace)、苊烯(acenaphthylene，Acy)、芴(fluorine, Flu)、菲(phenanthrene, Phe)、蒽(anthracene, Ant)、荧蒽(fluoranthene, Fla)、芘(pyrene, Pyr)、苯并[a]蒽(benzo(a)anthracene, BaA)、(chrysene, Chr)、苯并[b]荧蒽(benzo(b)fluranthene, BbF)、苯并芘(benzo(a) pyrene, BaP)、茚并[1,2,3-c,d]芘(indeno(1,2,3-cd)pyrene, InP)、二苯并[a,h]蒽(dibenz(ah)anthracene, DahA)和苯并[g,h,i]苝(benzo(g,h,i)perylene，BghiP)，均购于美国Accustandard公司.

1.2 样品采集

在乌兹别克斯坦境内的阿姆河下游河段及阿姆河三角洲范围内共设置采样点50个，其中阿姆河三角洲地区采样点45个；阿姆河下游河段采样点5个. 阿姆河三角洲位于阿姆河流域下游河段，在城市努库斯(42°27′11″N, 59°36′37″W)以下，河流分数支入咸海，形成面积达1万km2的广阔河口三角洲，是发达的农业灌溉区[15]. 阿姆河主要由冰雪融水和雨水供给，本次水样采集于2019年8月，具体采样位置及空间分布见图1. 依据《水质采样技术指导》(HJ 494-2009) 相关内容，利用表层采样器采集水样，控制每个样点采样深度为0.50 m左右，体积大于1 L，装入密封玻璃瓶中低温冷藏备用.

图1 阿姆河流域采样点分布Fig.1 Distribution of sampling sites in Amu Darya Basin

1.3 前处理方法与仪器分析

1 L水样经0.45 μm混合纤维素滤膜过滤去除杂质后, 利用固相萃取装置进行富集，选用富集柱为C18柱(Supelco, 500 mg/6 mL). 富集前用甲醇和纯水各5 mL对小柱浸泡淋洗进行活化，水样富集后抽干30 min，用二氯甲烷与正己烷(V∶V=7∶3)混合液20 mL进行洗脱，浓缩后用甲醇定容1 mL待测.

采用高效液相色谱二极管阵列检测器串联荧光检测器进行PAHs含量分析，荧光检测器除对Acy没有响应外，对其余15种PAHs响应灵敏，因此利用二极管阵列检测器测定Acy浓度，同时串联二极管阵列检测器测定的16种PAHs可以与荧光检测器依次对应，避免分析测试中的假阳性. 色谱柱为PAHs(4.6 μm ⅹ 250 mm，waters)专用柱；流动相为水和乙腈，流速为1 mL/min; 梯度洗脱程序：20 min内乙腈由初始60%上升至100%，维持10 min；DAD检测波长为238 nm，荧光检测波长为梯度变化[16].

1.4 质量控制与质量保证

按照1.3节分别利用16种多环芳烃标准样品获得各物质对应的保留时间，同时采用多点校正外标曲线法对实际样品中多环芳烃物质进行定性定量. 方法在1～200 μg/L间线性良好，R2均大于0.99，各物质检出限在0.12～0.41 μg/L之间. 实验过程通过方法空白、平行样品、标样插标及加标回收率进行质量控制与质量保证. 16种物质相对标准偏差RSD值在0.5%～2.3%之间，方法空白以无目标物检出为准，加标回收率在80%～104%之间，符合质控要求. 加标回收率是在空白样品中加入固定浓度的标准物质再进行完整前处理操作，通过仪器定量分析空白样品中标样浓度可以指示整个前处理操作流程，确定方法适用性与稳定性; 另外在实际测试过程中随机抽取已完成定量的样品或者对测试结果存疑的样品加入标准物质进行复测，确保样品定性分析准确性，详细参数见附表Ⅰ.

2 结果与讨论

2.1 含量与分布特征

由表1可知，16种PAHs在大多数水体样品中均有检出，研究区水体中∑PAHs浓度范围为3.19～779 ng/L，平均值和中位值分别为98.4和40.1 ng/L，单体浓度范围0～333 ng/L，检出浓度最高的单体为Acy. 从PAHs各单体组分检测结果分析可以看出，有5种单体物质检出率为100%，分别是Flu、Ant、Fla、Pyr和Chr，单体BbF的检出率为98%，但检出总量最高，水样中总浓度为786 ng/L，平均值为15.7 ng/L，中值为2.79 ng/L. 另外，DahA、BghiP、InP 3种单体物质的检出率较低，分别为76%、73%和61%，这主要由于低分子量多环芳烃具有较高的蒸汽压和水溶性，更易存在于水体中，而高分子量多环芳烃具有较低的水溶性和疏水性，致使水体中检出率较低[17]. 表1中高环PAHs (5～6环)总浓度与中低环PAHs (2～4环)总浓度相近，但各单体物质变异系数较高，高环PAHs浓度范围区间较大，表明研究区水体PAHs浓度空间差异大.

表1 阿姆河流域水体中PAHs检出情况统计

“ND”表示未检出；“-”表示无对应值.

从各采样点PAHs的总浓度和16种单体浓度(图2)可以看出，PAHs浓度较高的采样点主要集中在三角洲区域，平均值为104 ng/L，中值为41.5 ng/L. 下游河段区域PAHs浓度较低，平均值为47.4 ng/L，中值为33.4 ng/L，与文献中下游河口浓度高于中上游地区的特点相似[18-19]，其中三角洲地区PAHs浓度范围为12.6～779 ng/L，浓度较高的采样点主要集中在城市周边、农业灌溉区及阿姆河下游靠近咸海区域，最高浓度位于31#采样点城市钦博伊(Chimboy)附近及10#咸海滩涂，浓度分别为779和573 ng/L; 三角洲区域浓度较低的采样点主要位于阿姆河支流及三角洲边缘非城镇区域. 下游河段采样点PAHs浓度有明显的从上至下递增的趋势，浓度峰值出现在46#采样点，浓度为76.5 ng/L. 已有研究表明[20-21]，水中PAHs污染程度通常可以按PAHs总浓度分为4个水平：轻微污染(10～50 ng/L)、轻度污染(50～250 ng/L)、中度污染(250～1000 ng/L)和重度污染水平(>1000 ng/L). 研究区域50个采样点中达到轻度污染的共有18个，另有5个采样点处于中度PAHs污染水平.

图2 阿姆河流域不同水体PAHs及单体浓度(ADD：三角洲区域；ADL：下游河段区域)Fig.2 Total and mean PAHs concentrations in Amu Darya Basin (ADD: Amu Darya Delta; ADL: Amu Darya Lower stream)

表2列出了研究区水体与世界范围不同区域湖泊、河流流域水体中PAHs的浓度范围及平均浓度情况[18,22-34]. 通常人口稠密的城市或工业化地区水体中PAHs浓度普遍较高，如阿姆河流域PAHs平均浓度水平低于黄河三角洲、伊朗侯尔河(Kor River)地区，并远低于建有水利渡槽的印度戈麦蒂河(Gomti River)及肯尼亚维多利亚湖维南湾(Winam Gulf，Lake Victoria)港口区域中PAHs浓度，但高于美国东部沿海地区萨斯奎哈纳河(Susquehanna River)和德国易北河和威瑟河(Elbe and Weser Rivers)等非城市地区；同时阿姆河流域的PAHs总浓度和平均浓度均高于同处于中亚地区的哈萨克斯坦巴尔喀什湖流域(Ili-Balkhash Basin)，与墨西哥托斯湾(Todos Santos Bay)、中国长江下游支流以及意大利台伯河及河口 (Tiber Rivers and estuary)等区域浓度相近，总体而言，与其他研究区域相比，阿姆河流域PAHs浓度处于中等水平.

表2 近年来不同区域水体多环芳烃总浓度及平均浓度的比较

2.2 不同水体PAHs组成特征及来源解析

根据苯环个数将PAHs进行组成分类，利用不同苯环数PAHs相对丰度法可初步判断PAHs来源[35]，下游河段和三角洲两个区域水体中PAHs的组成特征如图3所示，所有采样点水体中PAHs的组成以中低环PAHs (2～4环)占比略高，其中下游河段水体PAHs组成中3环PAHs (Acy、Ace、Flu、Phe、Ant)占比为42%，检出率最高的化合物为Phe和Flu. Phe通常与焦油生产有关，而Flu的产生则与生物质燃烧相关，表明该区域可能存在焦油燃烧和生物质燃烧污染；三角洲区域水体PAHs组成中5环物质(BbF、BkF、BaP、DahA)占比最高，为32%，其次为4环(Fla、Pyr、BaA、Chr)和3环PAHs，所占比例分别为23%和18%，表明该区域水体中PAHs多来源于化石燃料和木材等生物质在高温下的不完全燃烧. 通常水体中4环以下PAHs的浓度较高，而5环和6环PAHs占优的水体，如西班牙内陆浅水湖、埃及尼罗河水体及黄河河口流域等，其污染源主要是化石燃料的高温燃烧污染[36-39]. 在三角洲区域中有12个采样点的水体中高环PAHs浓度高于低环PAHs，并且PAHs的污染程度也处于较高的水平，这些采样点位主要位于城镇、集中农业灌溉区、废弃码头及近咸海水面区域，可能受到交通、工业及煤炭能源的污染[40]. 阿姆河三角洲和下游河段两区域水体PAHs各成分间没有明显相关性(P>0.05)，表明三角洲区域水体PAHs的污染除下游河段汇入外，还存在其他污染源的输入，同时三角洲地区也因为常年径流量减少，总蒸发量大，大部分地区没有形成有效循环[41-43]，工业用水需求不断加大及农业灌溉用水回流等原因影响整体区域水环境中PAHs浓度.

图3 阿姆河下游河段区域(a)以及三角洲区域(b)PAHs物质组成Fig.3 Composition of PAHs in the Amu Darya lower stream area (a) and Amu Darya Delta (b)

目前，常用的源解析方法还包括同分异构体比值法[44]，研究区不同水样中同分异构体Ant/(Phe+Ant)和Flu/(Flu+Pyr)的比值见图4. 一般认为[45-46]，当Flu/(Flu+Pyr)比值小于0.4时，PAHs主要来源于石油及石油物质泄漏，当Flu/(Flu+Pyr)比值大于0.5时，PAHs主要来源于煤燃烧和包括草木在内的生物质燃烧；当Flu/(Flu+Pyr)比值介于0.4～0.5之间时，表明PAHs主要来源于高温燃烧，包括汽油、柴油和原油的燃烧以及机动车尾气排放，同时Ant/(Ant+Phe)的比值小于0.1和大于0.1也分别代表PAHs的石油来源和燃烧来源. 由图4可知，在阿姆河下游河段区域PAHs的Ant/(Phe+Ant)和Flu/(Flu+Pyr)的比值分别在0.1～0.3和0.4～0.6之间，主要集中在燃烧区域，可见该地区水体中PAHs主要来源于燃烧，且主要为煤炭和草木等生物质燃烧；而三角洲区域水体中PAHs的Flu/(Flu+Pyr)比值中有部分样点小于0.4，说明该区域存在着石油源污染，结合Ant/(Phe+Ant)比值可以判别阿姆河三角洲流域中部分地区存在着PAHs的石油与燃烧混合源污染，约占采样点位的15%左右，其他大部分地区的PAHs仍来源于燃烧，约40%的地区水体样品中PAHs来源与交通污染有关，比值法较好地补充验证了相对丰度法对两个区域PAHs来源的判断.

图4 乌兹别克斯坦阿姆河地区不同水体中PAHs来源同分异构体比值Fig.4 Cross plots for PAHs isomeric ratios in sample sites of the Amu Darya Basin，Uzbekistan

阿姆河三角洲地区PMF运行结果见图5，模型运行结果采用4因子，迭代数为17时定量解析PAHs来源，其Q(Robust)值136与Q(True)值140最为接近. 主因子1中Nap具有最高的因子载荷，其通常作为化工石油类物质泄漏指示物. 另外，具有较高因子载荷的Flu和Ace也同属低环PAHs，一般认为其来自石油泄漏及石油化工产品污染[50]；BbF和Inp等高环PAHs则为汽柴油高温燃烧产物[51]. 三角洲地区人口众多，努库斯、莫伊纳克、钦博伊等大中型城市每天产生大量生活污水和工业废水，同时沿咸海地带发展渔业，渔业船舶航行及沿岸排放废水可以引入大量石油类物质，因此化工石油类污染源为三角洲地区PAHs主要来源. 主因子2中主要载荷Acy是木柴燃烧的指示物，Ant的产生与燃煤有关; 该地区农业生产较为发达，其所属的卡拉卡尔帕克斯坦共和国的耕地面积是乌兹别克斯坦全国耕地面积第2大行政区[52]，木材秸秆等生物质也是其主要燃料来源. 主因子3中主要载荷BaA和Chr来自天然气燃烧，根据《BP世界能源统计年鉴(2019)》数据: 2019年乌兹别克斯坦天然气产量为566亿m3，占世界总产量的1.5%，消费量为426亿m3，是该区域的主要初级能源，Fla、Pyr和Phe则主要来源于焦炭燃煤燃烧. 主因子4中具有较高载荷因子主要为高环PAHs，与柴油汽油燃烧相关，通常为交通运输过程中机动车柴油汽油的燃烧、泄漏及尾气排放[53-55].

图5 阿姆河三角洲地区PMF源成分谱Fig.5 Source profiles obtained from the PMF model in the Amu Darya Delta area

阿姆河下游河段地区PMF运行结果如图6，模型运行结果采用2因子，Q(Robust)值与Q(True)值均为21.6. 主因子1中具有较高载荷的Acy为木柴等生物质燃烧源指示物，主因子2中高环PAHs单体InP、Bghip、DahA具有较高载荷，这些单体均来自于柴油高温燃烧的交通污染. 下游河段地区位于乌兹别克斯坦花拉子模州，主要依靠农业经济，棉植业是该区域农业的基础，同时畜牧业也在区域内占有举足轻重的经济地位，该行政区农村人口远高于城市人口[52]，煤炭、秸秆等为该区域居民的主要燃料来源，同时行政区内有通向首都塔什干的主要路线，交通运输活动繁忙，可以解释下游河段区域PAHs主要来源.

图6 阿姆河下游河段PMF源成分谱Fig.6 Source profiles obtained from the PMF model in the Amu Darya lower stream area

土壤沉积物中PAHs通常是长时间积累的结果，而水体环境中PAHs则可以在一定程度上反映近期环境行为. 近年来乌兹别克斯坦工业发展迅速，其工业地位在中亚地区举足轻重，天然气、有色金属等产业都较为发达; 同时作为农业大国，农业产值占国内生产总值的25%～30%，种植业尤其是棉花种植，仍是本国的经济支柱之一. 综合研究区域水体中PAHs分布特征，结合相对丰度法、同分异构体比值法和正定矩阵分解法3种源解析方法得出，三角洲地区工业生产过程中煤和石油天然气能源的燃烧、石油化学品泄漏及汽车尾气的排放等均导致了区域水体环境中PAHs浓度升高，阿姆河三角洲及下游河段地区水体中PAHs来源均有木材等生物质低温燃烧的贡献，但农业活动产生的生物质燃烧以及作为主要生活能源的木柴燃烧产生的PAHs占下游河段地区水体中PAHs的比重更大.

一般来说，枯水期水量下降、循环不畅等原因可能导致环境水体中PAHs浓度增加[56]，同时丰水期由于降水集中，雨量及地表径流量增加，促使包括工农业污染、城市交通污染在内的污染物随地表径流以及灌渠进入水体环境[57]，因此不同季节水体中PAHs浓度和组成变化各不相同. 乌兹别克斯坦境内阿姆河流域属大陆性气候，降雨主要集中在冬、春季节，本次样品采集在8月完成，处于该区域全年炎热干燥时期，水体蒸发作用增强同时PAHs从沉积物和悬浮颗粒物上解析作用加大[58]，可能导致PAHs浓度和组成变化，因此有必要对其他季节阿姆河流域水体的PAHs情况进一步研究.

2.3 水体污染水平及风险评价

USEPA水质标准规定了地表水中13种PAHs浓度限值[59]，用于评估水体中PAHs污染程度. 本研究中水体样品BaP的浓度均超过USEPA标准，同时分别有86%、75%、69%、45%和18%的水样中单体BbF、DahA、InP、BaA和BkF浓度超过标准限制；另外11个样点的∑PAHs总浓度超过欧盟规定的允许水生生物暴露安全限制的最大浓度[60]，占总采样点的21%，表明目前该区域水体中PAHs可能通过生物富集作用对人类健康产生影响.

本研究同时采用Kalf风险熵值法结合曹志国等和Li等[61-62]的研究，以忽略浓度(negligible concentrations，NCs)和最大允许浓度(maximum permissible concentrations，MPCs)作为16 种PAHs单体的参考浓度，通过计算风险熵值对阿姆河流域水体多环芳烃的生态风险进行评价，风险熵值(risk quotient，RQ)计算公式为:

RQ=cPAHs/cQV

(1)

RQNCs=cPAHs/cQV(NCs)

(2)

RQMPCs=cPAHs/cQV(MPCs)

(3)

式中，cPAHs为水体中各单体PAHs的质量浓度；cQV为各单体PAHs所对应的风险标准值，其中cQV(NCs)和cQV(MPCs)分别为最低风险标准值和最高风险标准值；RQNCs和RQMPCs则分别表示最低风险和最高风险浓度风险熵值. 对于单体PAHs,当RQNCs<1时，一般认为风险较低可以忽略，当RQNCs≥1且RQMPCs<1时，处于中等风险，应采取措施防止进一步污染；而当RQMPCs≥1时，可能已经造成了污染，需引起重视并积极干预控制，其中PAHs单体及∑PAHs风险分级情况见表3.

表3 单体PAHs与∑PAHs风险等级分级

阿姆河流域表层水体中单体PAHs主要处于3个等级(图7)：单体Nap、Ace、Phe和Ant的RQMPCs和RQNCs值均小于1，为低风险等级；单体BbF的RQNCs值和RQMPCs值均大于1，提示该类PAHs单体存在的风险较高；其余PAHs单体的RQMPCs值小于1而RQNCs大于1，提示已经具有一定的生态风险，应引起重视. 另外，根据研究区域内ΣPAHs生态风险熵值，采样点位中有12和8个点位分别处于中等风险2和高风险等级，高风险等级主要集中在阿姆河三角洲地区，由于阿姆河最终汇聚蓄积咸海区域，可能会加重咸海地区污染水平，因此需制定合理的防控计划.

图7 采样点风险熵值：单体PAHs生态风险熵值(a)；不同采样点ΣPAHs生态风险熵值(b)Fig.7 Ecological risks of PAHs in water samples collected in the study area: RQNCs and RQMPCs for single pollutants (a)，and ∑RQNCs and ∑RQMPCs for total pollutants (b)

3 结论

1)乌兹别克境内阿姆河流域采集的水体样品中，∑PAHs浓度在3.19～779 ng/L之间，均值为98.4 ng/L，中值40.1 ng/L，与世界不同区域的湖泊流域水体相比，PAHs浓度处于中等水平.

2)50个采样点主要集中于阿姆河三角洲和阿姆河下游河段两个区域，其中三角洲区域既有下游河段汇入蓄积又存在其他污染源的输入. 高环PAHs占比略高于中低环PAHs，但各采样点间浓度差异较大.

3)源解析分析PAHs来源显示，阿姆河下游河段地区PAHs主要来源于煤炭和草木等生物质燃烧；三角洲区域水体受生物质燃烧、石油燃烧的混合来源的污染，部分地区水体污染也与交通污染有关.

4)生态风险评估结果显示，PAHs总浓度超过欧盟最大允许浓度的规定限值，表明水体中PAHs可能会通过水生生物富集作用对人类健康构成一定威胁. 经风险熵值法分析，单体PAHs中Nap、Ace、Phe和Ant处于较低风险，而BbF的RQNCs值和RQMPCs值均大于1，污染程度较为严重，其余单体均列为中等风险等级;∑PAHs 风险熵值显示研究区域内近半数点位风险等级相对较低，中高风险等级点位主要集中在阿姆河三角洲地区，需引起重视.

4 附录

附表Ⅰ见电子版(DOI：10.18307/2022.0312).