松花江表层沉积物中16种多环芳烃空间分布特征及生态风险评价

李泽文， 王海燕， 孔秀琴， 杨艳艳， 钟 震， 李 莉， 吴 桐

1.兰州理工大学， 甘肃 兰州 730050

2.中国环境科学研究院环境污染控制工程技术研究中心， 北京 100012

3.中国环境科学研究院， 环境基准与风险评估国家重点实验室， 北京 100012

沉积物是河流生态系统的重要组成部分，作为污染物载体时刻与河水进行着物质和能量交换[1]. PAHs(多环芳烃)主要来源于环境中有机物的不完全燃烧[2]，具有低溶解性和疏水性，很容易与水体中悬浮物结合，沉降于河流沉积物中，因此沉积物是PAHs的主要聚集地，具有吸附和释放PAHs的功效[3-4]. PAHs污染具有来源广、残毒时间长、致癌性强的特点，对河流水生生物和人类健康造成了潜在威胁[5-6]，尤其是Nap(萘)、Ace(苊烯)、Acp(苊)等16种PAHs(∑16PAHs)为美国[7]和我国的优先控制污染物[8]，更应引起关注和研究.

国内外学者对沉积物PAHs污染分布特征及生态风险评价展开了大量研究，涉及河流、湖泊和近海水域等[9-14]. 沈小明等[9]研究长江河口沉积物中w(∑14PAHs)为83.4～5 206.9 ngg，部分点位具有潜在生态风险；陈明华等[10]利用比值法分析太湖竺山湾湖滨带w(∑16PAHs)主要来源于化石燃料的高温燃烧；冉涛等[11]发现，渤海湾中部海域w(∑16PAHs)呈现由近岸区向深海区递减的趋势；Singare等[12]等发现，印度Mithi河沉积物中w(∑16PAHs)为 1 206.0～3 705.0 ngg，具有一定的生态风险；Souza等[13]等发现，巴西Poxim河w(∑16PAHs)为2.2～28.4 ngg，主要来源于燃烧源；Guinan等[14]等发现，北爱尔兰Larne和Strangford海湾沉积物中w(∑13PAHs)为83～2 300 ngg，主要来源于燃烧源，部分点位具有潜在生态风险；Echols等[15]曾报道，美国Missouri河下游沉积物中w(∑27PAHs)〔除2个城市点位分别为 1 100 和 4 000 ngg外〕为250～700 ngg，主要来源于燃烧源；Duodu等[16]研究了澳大利亚Brisbane河表层沉积物15种PAHs的污染状况，认为其具有一定的致癌风险. 综上，不同地区不同水体中PAHs污染分布特征及生态风险均有所不同，但多数来源于燃烧源.

松花江是我国七大水系之一，属国家重点流域，也是东北地区居民生产、生活的主要水源[17]. 松花江来自北源嫩江和南源第二松花江，在吉林省松原市的三岔河附近汇合后，称松花江干流(该研究称松花江主段干流). 近年来，随着松花江流域社会经济的快速发展，松花江沉积物的污染物逐渐增多. 2015年黑龙江和吉林两省废水排放量占松花江全流域的95%左右[18-19]，但目前有关松花江流域沉积物中PAHs空间分布特征的研究主要集中在第二松花江的部分河段[20]和松花江主段干流[21-22]. CUI等[21]报道2014年松花江主段干流及第二松花江干流表层沉积物中w(∑16PAHs)为34～4 456 ngg，主要来源于煤炭燃烧；ZHAO等[22]研究松花江主干流(包括嫩江和第二松花江干流各1个点位、松花江主段干流13个点位)2007年表层沉积物中w(∑16PAHs)为68.25～654.15 ngg，总体处于低生态风险水平. 但有关松花江全流域(黑龙江省和吉林省)干流和支流表层沉积物中16种PAHs分布特征及其生态风险评价的研究尚未见报道，缺乏针对第二松花江、松花江主段干流及其典型支流沉积物中16种PAHs的污染特征揭示和风险状况评估.

因此，该研究以松花江全流域(黑龙江省和吉林省)干流和典型支流为对象，揭示了16种PAHs在松花江全流域干支流表层沉积物的空间分布特征，采用沉积物质量基准法和沉积物质量标准法对其生态风险进行了评价，并运用比值法和主成分分析法对其来源进行解析，以期为松花江水生态恢复研究提供科学基础和理论依据.

1 材料与方法

1.1 样品采集

2017年9月在松花江全流域(具体为第二松花江干流、松花江主段干流及其典型支流)进行表层沉积物样品的采集，共布设26个采样点，分布在流域上、中、下游区域(见图1和表1). 使用抓斗式采样器采集沉积物样品，每个采样点随机采集3份沉积物平行样. 采集的样品用锡箔纸包好，置于不锈钢盒中，在实验室冷冻保存，备用.

1.2 仪器与试剂

仪器：冷冻干燥机(美国Labconco公司)；SPE-DEX 4790全自动固相萃取仪(美国Horizon公司)；ASE-350 加速溶剂萃取仪(美国Dionex公司)；氮吹浓缩仪( 美国Caliper公司)；Agilent 7890A型气相色谱-质谱仪(GC-MS)(美国Agilent 公司).

主要试剂包括16种PAHs混标、代标和内标. 无水硫酸钠、硅藻土(干燥剂)和石英砂，使用前马弗炉450 ℃烘烤4 h后冷却至室温. 铜粉使用前先用1 molL硝酸氧化，去除铜粉表面的氧化物，再用甲醇冲洗3遍，氮气吹干. 氦气，纯度≥99.99%. 正己烷、环己烷、丙酮，均为分析纯.

1.3 样品预处理与分析

样品预处理包括样品冷冻干燥、研磨过筛、称量、加替代标准物、加速溶剂萃取、氮吹、PSA+硅胶小柱净化、氮吹、加内标、定容、GC-MS测定等步骤[23]：①将 冷冻干燥后的样品研磨成粉，过0.149 mm筛后称取5 g，加2 g硅藻土后置于加速溶剂萃取仪中进行萃取，萃取溶剂为体积比为1∶1的正己烷与丙酮溶液，萃取仪载气压力为0.8 MPa，加热温度为100 ℃，萃取池压力为8.3～13.8 MPa，预加热平衡5 min，静态萃取5 min，溶剂淋洗体积为60%池体积；氮气吹扫时间为60 s；静态萃取2次. ②将萃取后的萃取液氮吹浓缩至2 mL，进行铜粉脱硫，然后采用Aglient硅胶小柱串联PSA小柱，对萃取液进行纯化，以去除腐殖酸、色素和极性杂质等. 连接好固相萃取装置，在串联好的小柱上加入2 g无水硫酸钠脱水；然后加入20 mL正己烷进行柱子的活化，活化后进行上样，用40 mL体积比为1∶1的丙酮与正己烷溶液洗脱PAHs，收集洗脱液，然后再次氮吹缓慢吹至近干，用正己烷溶液定容到1 mL，用于GC-MS分析.

图1 松花江流域采样点分布Fig.1 Distribution of sampling sites in Songhua River Basin

表1 松花江流域采样点信息

注：松花江主段干流，简称主段干流；第二松花江简称二松.

GC-MS分析条件：进样口温度为280 ℃，不分流；进样量为1.0 μL；柱流量为1 mLmin；柱温为80 ℃保留2 min，以20 ℃min升至180 ℃，保持5 min，再以10 ℃min升至290 ℃，保持5 min. 离子源温度为230 ℃；离子化能量为70 eV；接口温度为280 ℃；四级杆温度为150 ℃；用选择离子扫描进行定量检测.

1.4 生态风险评价方法

1.4.1沉积物质量基准法

沉积物质量基准法常用的两个指标是ERL(effects range low, 效应范围低值)和ERM(effects range media, 效应范围中值)[24]. ERL和ERM将PAHs的生态风险分为3个等级：当有机物浓度高于ERM时，生态风险发生的概率较高(大于50%)，对周边生物及生态环境会产生经常性生态风险；当污染物浓度大于ERL小于ERM时，对周边生物体及生态环境偶尔产生生态风险；当研究区域污染物浓度低于ERL时，发生生态风险的概率较低(小于10%)，几乎不会对周边人类及生态环境产生危害[25].

1.4.2沉积物质量标准法

沉积物质量标准法也是一种广泛用于评价沉积物中PAHs生态风险的方法[26]. 目前我国尚未颁布，一般采用2006年加拿大魁北克省颁布的沉积物质量标准法，该标准设定REL(罕见效应浓度值)、TEL(临界效应浓度值)、OEL(偶然效应浓度值)、PEL(可能效应浓度值)和FEL(频繁效应浓度值)5个阈值来评价沉积物中PAHs的污染程度和生态风险(见表2).

表2 加拿大魁北克省淡水沉积物中PAHs的质量评价标准

1.5 来源分析

1.5.1特征分子比值法

PAHs的污染来源常运用同分异构体特征比值法来判断[27]，常用比值有An(Phe+An)、Flu(Flu+Pyr)和BaA(BaA+Chr)，其不同特征比值范围代表沉积物中PAHs的主要来源(见表3).

表3 同分异构比值法鉴定PAHs污染来源

1.5.2主成分分析法

主成分分析是多元统计分析中一种非常重要的分析方法，对原始数据进行降维处理分析，将变量之间的相关性系数组合成较少的因子，通过因子之间的相关性，来确定沉积物中PAHs的污染来源[28].

2 结果与讨论

2.1 沉积物PAHs含量分布

松花江全流域26个采样点表层沉积物中16种PAHs均有检出，其质量分数(以干质量计)范围为169.76～3 769.19 ngg，平均值为 1 458.43 ngg(见表4). 其中w(BaA)为19.80～950.06 ngg，平均值为327.81 ngg，明显高于其他化合物.w(Fl)、w(Acp)、w(Ace)和w(BkF)相对较低，变化范围分别为3.77～60.44、22.67～113.89、11.19～95.32、11.45～130.15 ngg，平均值分别为12.21、55.22、38.57和62.83 ngg. 其中干流采样点w(∑16PAHs)为459.92～2 092.58 ngg，平均值为 1 173.67 ngg，支流采样点w(∑16PAHs)为169.76～3 769.19 ngg，平均值为 1 598.41 ngg，支流采样点w(∑16PAHs)明显高于干流，并且支流采样点每种PAHs的平均值均高于干流.

松花江上、中、下游采样点w(∑16PAHs)分别为478.06～3 769.19、459.92～2 522.94 和169.76～1 269.82 ngg，平均值分别为 1 818.13、1 364.72 和616.62 ngg，总体呈从上游到下游逐渐降低的趋势. 这可能是因为松花江上游曾集中了大量的石油化工等工业企业，而且有石油开采行业[18-19]，因此上游w(∑16PAHs)相对偏高；松花江下游工业相对没有上游发达，PAHs污染主要受居民生活污染等因素影响，w(∑16PAHs)相对较低.w(∑16PAHs)的最大值出现在S8号采样点(松花江上游第二松花江支流辉发河楞场村)，这可能与该采样点位于白山市(人口密集、交通发达和工业企业较多)下游有关；最小值出现在S24号采样点(松花江下游一级支流倭肯河上的安兴村)，这与该采样点周围自然环境较好、远离城市和受人类活动影响较小有关.

表4 松花江表层沉积物中PAHs质量分数

从各江段w(∑16PAHs)来看，松花江主段支流(1 622.47 ngg)>松花江主段干流(1 509.86 ngg)，第二松花江支流(2 143.01 ngg)>第二松花江干流(928.72 ngg)，针对最大支流牡丹江，牡丹江支流(1 302.54 ngg)>牡丹江干流(1 069.95 ngg)，总体而言，支流w(∑16PAHs)高于干流，这主要是因为支流多为城市内河，接纳的污染物相对集中，同时由于支流流量较小，PAHs等有机物更容易被积累于沉积物中，而干流流量相对较大、含沙量高、有机物易稀释扩散；并且第二松花江干支流的w(∑16PAHs)(1 651.89 ngg)高于松花江主段干支流(1 588.68 ngg)，可能与第二松花江石化企业集中有关；松花江主段干流和第二松花江干流表层沉积物中w(∑16PAHs)与CUI等[21]研究结果接近.

研究区域沉积物中PAHs的质量分数和污染程度与当地的经济发展和人类活动密切相关，工业越发达、人口越密集，PAHs质量分数越高. 与国内外其他河流相比(见表5)，松花江表层沉积物污染水平处于中低水平，其PAHs质量分数低于我国的珠江、大辽河、墨水河和印度的Mithi河，与南黄海和第二松花江接近，但高于太湖、莱州湾、渤海北部和滦河口.

表5 松花江表层沉积物中PAHs质量分数与国内外其他水体对比

2.2 沉积物PAHs组成特征

PAHs根据苯环数可以分为2环、3环、4环、5环和6环PAHs. 一般而言，小分子的2环和3环PAHs主要来源于石油化工污染源，而高分子的PAHs(4环以上)主要来源于燃烧过程[39]. 由于PAHs性质的差异，其环境行为也有所不同，表现出的毒性和致癌性也存在差异[40-41]，因此研究PAHs的组成特征，对于保障河流生态系统健康具有重要意义.

松花江流域干支流26个表层沉积物中2、3、4、5、6环的PAHs分别占w(∑16PAHs)的3.77%、22.98%、45.16%、14.08%和14.01%，占比从大到小依次为4环>3环>5环>6环>2环. 除了S13、S21和S24号采样点以3环为主，占比范围为26%～39%，其他采样点都是以4环为主，其中S2、S12、S14和S17号采样点4环占比均高于50%. 总体而言，松花江表层沉积物4环以上PAHs占比较高，为45.16%(见图2)，可以初步判断PAHs主要来源于各种物质的燃烧过程[39]. 与国内其他河湖相比，研究区沉积物各环PAHs组成分布特征与长江口[42]及千岛湖[43]完全一致(4环>3环>5环>6环>2环).

2.3 PAHs的来源分析

2.3.1特征分子比值法

2.3.2主成分分析

对松花江流域26个表层沉积物进行主成分分析，提取3个主成分(PC1、PC2、PC3)，计算的累计方差贡献率如表6所示. 由表6可见，主成分PC1方差贡献率为62.22%，主要由2环～6环的PAHs构成，Nap、Ace、Acp、Fl、Phe、An、Flu、Pyr、Chr、BbF、BkF、BaP、DbA和BgP具有较大的负荷值，其中Nap、Ace、Acp、Fl和Flu主要是木材、秸秆等生物质燃烧的指示物[46]，而高环Phe、Pyr、Chr、BbF与BkF主要来源于煤炭燃烧和交通污染[47]，因此推断PC1是燃烧混合源. PC2方差贡献率为9.73%，InP具有较大的负荷值. Larsen等[47]发现，InP主要来自柴油燃烧源，因此PC2为柴油燃烧源. PC3方差贡献率为7.99%，BaA和InP具有较大的负荷值，BaA是天然气排放的重要指标[29]，因此推断PC3为煤炭和天然气燃烧污染源. 主成分分析结果表明，松花江流域沉积物中的PAHs主要受木材、煤炭、柴油和天然气等燃烧源为主的混合源影响，石油泄漏污染较少，与2.3.1节解析结果相一致.

图3 松花江表层沉积物中PAHs来源诊断Fig.3 Diagnosis of PAHs sources in surface sediments of Songhua River

表6 松花江表层沉积物中PAHs主成分分析

2.4 沉积物PAHs的生态风险评价

2.4.1沉积物质量基准法

沉积物质量基准法相对比较成熟[48-49]，该研究采用沉积物质量基准法对松花江表层沉积物PAHs的生态风险进行评价，结果如表7所示.

表7 松花江表层沉积物中PAHs质量分数与ERL和ERM的比较

Table 7 Comparison of PAHs concentration with ERL and ERM values in surface sediments of Songhua River ngg

表7 松花江表层沉积物中PAHs质量分数与ERL和ERM的比较

PAHsERLERM最大值最小值平均值Nap1602 100173.674.1443.51Ace4464095.3211.1738.57Acp16500113.8922.6755.22Fl1954060.443.7712.22Phe2401 500633.2818.11160.19An85.31 100151.992.8268.90Flu6005 100287.892.6988.93Pyr6652 600227.4424.4275.23BaA2611 600950.5819.8339.34Chr3842 800318.689.8497.19BaP4301 600367.524.25136.13DbA63.4260204.955.6963.97

由表7可见，只有w(Acp)、w(BaA)和w(DbA)超过了ERL，位于ERL和ERM之间，表明其对周边生物体及生态环境会偶尔产生生态风险[24]，但BaA和DbA是化学强致癌物，未来应重点研究该污染物在松花江流域的来源和途径，并提出相应的管控措施；其余PAHs组分的平均浓度均在ERL之下，表明这些PAHs组分具有较低的生态风险水平，发生生态风险的概率较低(小于10%).

Nap、Ace、Acp、An、Fl、Phe、BaA和DbA这8种PAHs组分质量分数的最大值均超过了ERL，位于ERL和ERM之间，说明研究区内有个别采样点对周边生物体及生态环境会偶尔产生生态风险. 由图4可见，S1、S2、S3、S7、S8、S14、S15和S21号采样点的8种PAHs组分质量分数较高，其中S1、S2、S3、S7和S8号采样点位于松花江上游，表明松花江上游PAHs质量分数较高. 这些地区曾聚集大量的石化企业，可能是PAHs污染的主要来源之一，建议未来应控制流域周边的石化企业的数量，同时加强对上述企业PAHs等有毒有机物的管控，从源头上减少PAHs的排放. 倭肯河上的S24号采样点的PAHs质量分数最低，为138.53 ngg，表明受到人类活动影响较小. Flu、Pyr、Chr和Bap组分质量分数的最大值均低于ERL，表明研究区内26个采样点的Flu、Pyr、Chr和Bap组分具有较低的生态风险水平，发生生态风险的概率较低. BbF、BkF、InP和BgP这4种PAHs组分没有规定最低安全值，需加以重视对其来源的管控.

2.4.2沉积物质量标准法

根据沉积物质量标准法，对松花江流域26个表层沉积物的12种PAHs生态风险状况进行评价，把26个采样点中至少有一种PAHs组分质量分数大于FEL的采样点列于表8第一行，至少有一种PAHs组分质量分数介于PEL和FEL之间的采样点列于表8第二行，依次类推，将各采样点的污染程度区分在各阈值之间(见表8).

由表8可见，松花江流域26个表层沉积物中，有3个采样点(S2、S14和S15)存在高生态风险，其分别位于松花江流域中上游的拉林河和蚂蚁河支流上，其中S2号采样点位于哈尔滨市双城区周边，S14和S15号采样点位于哈尔滨市下游的通河县附近，均处于人类活动密集、交通发达的城市段，PAHs的污染来源广，因此沉积物中PAHs污染程度较高，相关部门应该予以重视，采取一定的措施，减少S2、S14和S15号采样点PAHs污染物的排放. 有6个采样点(S1、S3、S5、S8、S17和S21)表层沉积物(23.07%)至少有一种PAHs的平均浓度介于PEL和FEL之间，潜在风险几率较高；有10个采样点(S4、S7、S9、S10、S13、S18、S19、S20、S23和S25)表层沉积物(38.46%)至少有一种PAHs的平均浓度介于OEL和PEL之间，潜在风险几率中等. 剩余7个采样点介于TEL和OEL之间，潜在风险几率较低. 研究区沉积物PAHs生态风险污染程度与长江口[50]基本持平，略高于太湖地区[51]，潜在风险水平中等. 总体而言，松花江流域有部分采样点存在一定的生态风险，相关部门应该采取相应的措施减少上述采样点PAHs的排放和治理.

表8 松花江表层沉积物中PAHs的污染程度分析

3 结论

a) 松花江干支流表层沉积物中w(∑16PAHs)范围为169.76～3 769.19 ngg，平均值为 1 458.43 ngg，以3～6环高环为主；并且支流w(∑16PAHs)(范围为169.76～3 769.19 ngg，平均值为 1 598.41 ngg)高于干流(范围为459.92～2 092.58 ngg，平均值为 1 173.67 ngg)，与国内外其他河流相比，处于中低污染水平. 沉积物中w(∑16PAHs)总体呈现从上游到下游逐渐降低趋势，并且支流的含量大于干流.

b) 松花江流域表层沉积物PAHs主要来源于生物质燃烧和石油燃烧，这与研究区域内产业布局和能源结构一致.

c) 松花江流域表层沉积物PAHs整体风险水平为低生态风险，位于中上游支流的3个采样点至少一种PAHs组分质量分数存在高生态风险，可能与石油化工企业等的排放及燃烧源相关，应采取相应的管控措施.