闽江福州段沉积物多环芳烃：分布、溯源与生态风险解析

闽江福州段沉积物多环芳烃：分布、溯源与生态风险解析一、引言1.1研究背景与意义多环芳烃（PolycyclicAromaticHydrocarbons，PAHs）是一类由两个或两个以上苯环以线性、角状或簇状排列组成的碳氢化合物，是最早被发现且数量最多的一类持久性有机污染物。这类物质性质稳定，难以降解，能在环境中长期存在，并通过大气传输、地表径流等途径广泛分布于各种环境介质中，如土壤、水体、大气以及生物体。PAHs具有亲脂疏水的特性，极易吸附于颗粒物表面，进而在沉积物中大量累积。PAHs对生态环境和人类健康存在严重威胁。其具有较强的急性毒性、遗传毒性、致癌性和光致毒性，会对生物体的免疫系统、呼吸系统、循环系统以及神经系统造成损害，还会对肝脏、肾脏等器官功能产生不良影响。有研究表明，长期暴露于PAHs环境中，人类患癌症的风险显著增加，例如苯并[a]芘被国际癌症研究机构列为1类致癌物，是已知致癌性最高的PAHs之一。在水生生态系统中，PAHs会影响水生生物的生长、发育和繁殖，导致鱼类等生物的畸形、病变甚至死亡，进而破坏整个生态系统的结构和功能。闽江作为福建省的第一大河流，是福州的母亲河，闽江福州段更是福州重要的发展廊道与景观、生态轴线。闽江福州段起于闽清县雄江镇，止于连江县琯头镇长门口，长度约125公里，它不仅承担着饮用水源、工农业用水、航运、渔业养殖等多种重要功能，还在维持区域生态平衡、调节气候等方面发挥着不可替代的作用。然而，随着福州市经济的快速发展和城市化进程的加速，闽江福州段面临着日益严峻的环境污染问题，其中PAHs污染不容忽视。工业废水排放、城市生活污水排放、交通运输业发展以及周边农业活动等，都可能导致PAHs进入闽江福州段水体，并在沉积物中不断积累。研究闽江福州段沉积物中PAHs的分布、来源及其生态风险具有重要的现实意义。通过对PAHs分布特征的研究，可以了解其在闽江福州段的污染现状，明确污染严重的区域，为针对性的污染治理提供依据。解析PAHs的来源，能够准确识别主要污染源，有助于制定有效的污染控制策略，从源头减少PAHs的排放。对PAHs进行生态风险评价，可以评估其对闽江福州段生态系统和人类健康的潜在危害程度，为生态环境保护和风险管理提供科学支撑，保障闽江福州段的生态安全和可持续发展，这对于维护福州地区的生态平衡、保障居民的健康生活以及促进区域经济的绿色发展都具有至关重要的作用。1.2国内外研究现状多环芳烃作为一类备受关注的持久性有机污染物，长期以来都是国内外环境科学领域的研究重点。国外对PAHs的研究起步较早，在20世纪60年代，就有学者开始关注PAHs对环境和生物的潜在危害。早期的研究主要聚焦于PAHs的分析检测技术开发，随着气相色谱-质谱联用（GC-MS）等先进仪器分析技术的不断发展和完善，使得对环境样品中痕量PAHs的准确测定成为可能，为后续的研究奠定了坚实基础。例如，美国环境保护署（USEPA）早在1979年就将16种PAHs列为优先控制污染物，并制定了一系列相关的监测和分析方法标准，推动了全球范围内对PAHs的研究热潮。在PAHs的来源解析方面，国外学者通过大量的研究，建立了多种成熟的源解析方法。如利用特征比值法，通过分析PAHs中不同单体之间的比值，如荧蒽/（荧蒽+芘）、苯并[a]蒽/（苯并[a]蒽+䓛）等，来初步判断PAHs的来源是石油源、燃烧源还是生物源等。美国地质调查局（USGS）的研究人员运用该方法对美国多个地区水体和沉积物中的PAHs进行源解析，发现城市地区PAHs主要来源于交通尾气排放和工业活动，而偏远地区则更多受到自然源和长距离传输的影响。此外，正定矩阵因子分解（PMF）模型等多元统计分析方法也被广泛应用于PAHs的源解析研究，能够更准确地定量解析出不同污染源对环境中PAHs的贡献比例。欧洲的一些研究团队利用PMF模型对多个城市大气中的PAHs进行源解析，识别出机动车尾气排放、生物质燃烧、工业排放等是主要的污染源。对于PAHs的生态风险评价，国外也开展了大量深入的研究。美国国家海洋和大气管理局（NOAA）提出了效应区间低值（ERL）和效应区间中值（ERM）等生态风险评价指标，用于评估沉积物中PAHs对水生生物的潜在危害程度。当沉积物中PAHs含量低于ERL时，被认为对生物产生负面影响的可能性较低；当含量介于ERL和ERM之间时，具有一定的潜在风险；当超过ERM时，则存在较高的生态风险。许多欧洲国家也建立了自己的生态风险评价体系，如荷兰的土壤环境质量标准中对PAHs的限值做出了明确规定，用于指导本国的环境管理和污染治理工作。国内对PAHs的研究虽然起步相对较晚，但近年来发展迅速。在PAHs的分布特征研究方面，国内学者对不同环境介质，如大气、水体、土壤、沉积物等中的PAHs进行了广泛的调查分析。在大气方面，对北京、上海、广州等大城市的研究表明，PAHs浓度呈现明显的季节变化，冬季由于取暖等原因，浓度普遍高于夏季，且在城市中心区域和交通繁忙地段浓度较高。在水体研究中，对长江、黄河、珠江等主要河流的监测发现，不同河段PAHs的含量存在差异，受工业活动和城市污水排放影响较大的河段污染相对严重。在土壤研究方面，对一些工业污染场地和农业用地的调查显示，PAHs含量与土地利用类型、污染源距离等因素密切相关。在PAHs的来源解析方面，国内学者结合我国的实际情况，综合运用多种方法进行研究。例如，通过对不同污染源排放特征的分析，建立了适合我国国情的特征比值数据库，并将其应用于实际样品的源解析中。同时，也积极引入国外先进的源解析模型，如PMF模型、UNMIX模型等，并对模型进行优化和改进，使其更适用于我国复杂的环境状况。在生态风险评价方面，国内在借鉴国外经验的基础上，也在不断探索建立适合我国国情的评价体系。一些研究将PAHs的生态风险与我国的生态保护红线、环境功能区划等相结合，为环境管理提供更具针对性的建议。然而，目前针对闽江福州段沉积物中PAHs的研究仍存在一定的不足。已有研究在采样点的分布上不够全面，无法完整地反映整个闽江福州段PAHs的分布特征，特别是一些偏远支流和生态敏感区域的监测数据相对匮乏。在来源解析方面，虽然已有部分研究采用了特征比值法等方法，但缺乏对多种源解析方法的综合对比和验证，导致对污染源的识别和贡献率的确定存在一定的不确定性。在生态风险评价方面，多采用单一的评价方法，缺乏多种评价方法的综合应用和对比分析，难以全面准确地评估闽江福州段沉积物中PAHs的生态风险。此外，针对闽江福州段独特的地理环境、水文条件以及区域经济发展特点，对PAHs的迁移转化规律及其与生态系统相互作用机制的研究还较为薄弱。本文将在这些方面展开深入研究，弥补现有研究的不足，为闽江福州段的环境保护和污染治理提供更科学、全面的依据。1.3研究内容与方法1.3.1研究内容本研究以闽江福州段沉积物为研究对象，全面深入地开展多环芳烃（PAHs）相关研究，具体内容如下：闽江福州段沉积物中PAHs的分布特征研究：在闽江福州段沿线，综合考虑河流的不同河段、支流汇入情况、周边土地利用类型以及污染源分布等因素，设置多个具有代表性的采样点。使用抓斗式采泥器采集表层沉积物样品，深入分析沉积物中16种美国环保署（USEPA）优先控制的PAHs单体的含量，包括萘（NAP）、苊烯（ANY）、苊（ANA）、芴（FLU）、菲（PHE）、蒽（ANT）、荧蒽（FLT）、芘（PYR）、苯并[a]蒽（BaA）、䓛（CHR）、苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IPY）、二苯并[a,h]蒽（DBA）、苯并[ghi]苝（BPE）。运用地理信息系统（GIS）技术，结合克里金插值法，绘制PAHs的空间分布图，直观展示其在闽江福州段沉积物中的空间分布规律，明确高污染区域和低污染区域。同时，分析不同季节采集的沉积物样品中PAHs含量的变化，探究PAHs含量的季节变化特征及其可能原因，如季节变化导致的河流流量、水温、污染源排放强度的改变等对PAHs分布的影响。闽江福州段沉积物中PAHs的来源解析：综合运用多种源解析方法，对闽江福州段沉积物中PAHs的来源进行准确识别和定量分析。首先采用特征比值法，计算荧蒽/（荧蒽+芘）、苯并[a]蒽/（苯并[a]蒽+䓛）等特征比值，根据不同比值范围对应的PAHs来源指示意义，初步判断PAHs的来源是石油源、燃烧源（包括生物质燃烧、化石燃料燃烧等）还是生物源等。同时，利用主成分分析（PCA）和多元线性回归（MLR）相结合的方法，对PAHs各单体含量数据进行统计分析，提取主要成分，并通过多元线性回归确定各主要成分对应的污染源及其贡献率。此外，还将结合研究区域的实际情况，如工业企业分布、交通运输状况、能源消费结构等信息，对源解析结果进行验证和补充，提高源解析的准确性和可靠性。闽江福州段沉积物中PAHs的生态风险评价：运用多种生态风险评价方法，对闽江福州段沉积物中PAHs的生态风险进行全面、综合的评价。采用美国国家海洋和大气管理局（NOAA）提出的效应区间低值（ERL）和效应区间中值（ERM）作为评价标准，将沉积物中PAHs各单体含量与ERL和ERM值进行对比，判断PAHs对水生生物产生负面影响的可能性及风险程度。当PAHs含量低于ERL时，认为对生物产生负面影响的可能性较低；介于ERL和ERM之间时，具有一定的潜在风险；超过ERM时，则存在较高的生态风险。运用毒性当量因子（TEF）法，将不同PAHs单体的含量换算成具有相同毒性的苯并[a]芘的当量浓度，即毒性当量浓度（TEQ），综合评估PAHs的总体毒性和生态风险。此外，还将采用概率风险评价方法，考虑PAHs含量的不确定性和生物对PAHs毒性响应的不确定性，通过蒙特卡罗模拟等方法，计算不同风险水平下的发生概率，为风险管理提供更全面的信息。闽江福州段沉积物中PAHs分布、来源及生态风险的影响因素分析：深入探讨影响闽江福州段沉积物中PAHs分布、来源及生态风险的多种因素。分析沉积物的理化性质，如粒度组成、有机质含量、阳离子交换容量等对PAHs吸附和解吸过程的影响，进而影响其在沉积物中的分布。研究河流的水文条件，如流量、流速、水位变化等对PAHs的迁移扩散和沉积过程的作用。考虑周边污染源的类型、排放强度和分布情况，以及人类活动，如工业生产、交通运输、农业活动、城市建设等对PAHs输入的影响。通过相关性分析、冗余分析（RDA）等方法，确定各影响因素与PAHs分布、来源及生态风险之间的定量关系，明确主要影响因素。同时，分析不同影响因素之间的相互作用关系，为制定针对性的污染控制和生态保护措施提供科学依据。1.3.2研究方法采样方法：在闽江福州段，按照均匀布点与重点区域加密布点相结合的原则，设置[X]个采样点。使用不锈钢抓斗式采泥器采集表层0-20cm的沉积物样品，每个采样点采集3个子样，将其混合均匀后作为该采样点的代表样品。采样时，详细记录采样点的地理位置（经纬度）、周边环境信息（如是否靠近工厂、居民区、交通要道等）以及采样时间、天气状况、河流流速等现场数据。样品采集后，立即装入经预先清洗和烘干处理的棕色玻璃瓶中，密封保存，并尽快运回实验室，放入冰箱中在4℃下冷藏保存，以防止样品中PAHs的挥发和降解。分析测试方法：将采集的沉积物样品在通风橱中自然风干，去除水分后，用玛瑙研钵研磨至粉末状，过100目筛备用。采用加速溶剂萃取（ASE）技术提取沉积物中的PAHs。称取适量研磨后的样品，加入适量的硅藻土混合均匀，装入萃取池中，以正己烷-丙酮（体积比为1:1）为萃取剂，在10.3-13.8MPa的压力和100-120℃的温度下进行萃取。萃取液经旋转蒸发浓缩后，用硅胶柱进行净化处理。将浓缩后的萃取液转移至硅胶柱顶部，依次用正己烷、正己烷-二氯甲烷（体积比为1:1）等洗脱剂进行洗脱，收集含有PAHs的洗脱液。将洗脱液再次浓缩至近干，用正己烷定容至1mL，转移至进样瓶中，供气相色谱-质谱联用仪（GC-MS）分析。使用GC-MS对PAHs进行定性和定量分析。采用HP-5MS毛细管色谱柱（30m×0.25mm×0.25μm），进样口温度为280℃，分流比为10:1，进样量为1μL。程序升温条件为：初始温度为60℃，保持1min，以15℃/min的速率升温至280℃，保持10min。质谱采用电子轰击离子源（EI），离子源温度为230℃，扫描方式为选择离子扫描（SIM），通过与标准物质的保留时间和质谱图对比进行定性分析，外标法进行定量分析。来源解析方法：特征比值法：根据相关研究文献，确定荧蒽/（荧蒽+芘）、苯并[a]蒽/（苯并[a]蒽+䓛）等特征比值与PAHs来源的对应关系。例如，当荧蒽/（荧蒽+芘）比值小于0.4时，指示PAHs主要来源于石油源；比值在0.4-0.5之间时，为石油源和燃烧源的混合来源；比值大于0.5时，则主要来源于燃烧源。计算闽江福州段沉积物样品中相应的特征比值，初步判断PAHs的来源类型。主成分分析（PCA）和多元线性回归（MLR）：运用统计分析软件（如SPSS）对PAHs各单体含量数据进行主成分分析，提取主成分。通过分析主成分中各PAHs单体的载荷系数，确定主成分所代表的潜在污染源。然后，以各主成分得分作为自变量，以PAHs总量为因变量，进行多元线性回归分析，建立回归方程，计算各潜在污染源对PAHs总量的贡献率。生态风险评价方法：效应区间低值（ERL）和效应区间中值（ERM）法：将沉积物中PAHs各单体含量与NOAA提出的ERL和ERM值进行对比。若PAHs单体含量低于ERL值，判定对水生生物产生负面影响的可能性较低；介于ERL和ERM值之间，具有一定潜在风险；超过ERM值，则存在较高生态风险。同时，计算PAHs总量的ERL商值（PAHs总量/ERL总量）和ERM商值（PAHs总量/ERM总量），进一步评估整体生态风险程度。毒性当量因子（TEF）法：根据国际上公认的TEF值，将沉积物中不同PAHs单体的含量换算成苯并[a]芘的毒性当量浓度（TEQ）。计算公式为：TEQ=\sum_{i=1}^{n}(C_{i}\timesTEF_{i})，其中C_{i}为第i种PAHs单体的浓度，TEF_{i}为第i种PAHs单体的毒性当量因子。通过比较计算得到的TEQ值与相关评价标准，评估PAHs的总体毒性和生态风险水平。概率风险评价方法：利用蒙特卡罗模拟软件（如CrystalBall），考虑PAHs含量的不确定性和生物对PAHs毒性响应的不确定性。设定PAHs含量和毒性响应参数的概率分布函数，进行多次模拟计算，得到不同风险水平下的发生概率。例如，设定PAHs含量服从正态分布，毒性响应参数服从对数正态分布等，通过模拟计算得到沉积物中PAHs对水生生物产生不同程度毒性效应的概率，为风险管理提供更全面的风险信息。二、闽江福州段概况2.1自然地理特征闽江福州段地处福建省东部，位于东经118°08′-120°31′，北纬25°15′-26°29′之间，作为福建省第一大河流闽江的重要组成部分，其地理位置独特，是福州地区重要的生态廊道和经济发展轴线。闽江福州段地势呈现出西北高、东南低的态势，整体地形以山地和丘陵为主，在河流的长期作用下，形成了河谷盆地和冲积平原等多样化的地貌类型。闽江自西北向东南贯穿其中，在流经市区时，于淮安附近分为南、北两支，北支称为北港（又称台江），南支称为南港（又叫乌龙江）。南、北港分别流经南台岛后，在马尾港汇流，之后折向东北，流经闽安镇、亭江，被琅歧岛分隔，最终分别从长门、梅花两水道注入东海。这种复杂的地形地貌对多环芳烃（PAHs）在闽江福州段的分布有着重要影响。一方面，山地和丘陵地区的地表径流会将周边的PAHs携带至河流中，而河谷盆地则可能成为PAHs的沉积区域，使得PAHs在这些区域的沉积物中更容易积累。另一方面，不同地貌类型下的土壤质地和植被覆盖情况不同，也会影响PAHs的来源和迁移转化过程。例如，植被丰富的山地地区，生物质燃烧可能是PAHs的一个重要来源，而在人口密集的平原地区，工业活动和交通运输等人为源则可能占据主导。闽江福州段属于亚热带季风气候，温暖湿润，雨量充沛。年平均气温在16-20℃左右，最冷月（1-2月）平均气温达6-10℃，最热月（7-8月）平均气温为24-29℃。年平均日照数为1700-1980小时，年平均降水量为900-2100毫米，年相对湿度约77%。这种气候条件对PAHs的分布同样有着显著影响。在高温多雨的夏季，河流的流量增大，流速加快，可能会促进PAHs的稀释和扩散，使得沉积物中PAHs的含量相对降低。而在冬季，气温较低，河流流量减少，PAHs更容易在沉积物中富集。此外，降水还可能通过淋溶作用，将大气中的PAHs带入河流，增加河流中PAHs的输入量。闽江福州段的水文特征复杂多样。其多年平均径流量较大，水量丰富，这使得河流对PAHs具有一定的稀释和搬运能力。然而，闽江源短流急，平均约三年就要发生一次超2万立方米每秒的较大洪水。在洪水期，河流的流速和流量急剧增加，可能会导致沉积物的再悬浮和PAHs的重新释放，使得水体中PAHs的浓度升高。同时，洪水还可能将大量的PAHs从上游地区携带至下游，影响PAHs在闽江福州段的空间分布。此外，闽江福州段受潮汐影响明显，在河口地区，潮汐的涨落会改变河流的水流方向和流速，进而影响PAHs的迁移和沉积过程。在涨潮时，海水携带的PAHs可能会进入河流，而在落潮时，河流中的PAHs则可能会随着水流进入海洋。2.2社会经济概况闽江福州段作为福州地区重要的发展轴线，其周边区域的社会经济活动十分活跃，对多环芳烃（PAHs）的排放有着显著影响。在人口分布方面，闽江福州段周边人口密集，福州市作为福建省的省会，2024年末常住人口达850.1万人。市区内众多的居民点沿闽江两岸分布，特别是在北港和南港流经的区域，如台江区、仓山区等，人口密度较大。大量人口的日常生活产生了各类污染物，其中包括PAHs。例如，居民生活中的取暖、烹饪等活动，若使用煤炭、木材等燃料，会产生含有PAHs的烟尘，这些烟尘通过大气沉降等方式进入闽江，进而在沉积物中积累。人口密集区域的生活污水排放中也可能含有一定量的PAHs，如一些塑料制品、个人护理产品等在使用过程中可能会释放出PAHs，随污水进入河流。工业布局对闽江福州段PAHs排放的影响更为直接。闽江福州段沿线分布着多个工业园区和工业集中区，涵盖了化工、机械制造、电子、食品加工等多个行业。化工行业在生产过程中，如石油炼制、有机合成等，会产生大量的PAHs作为副产物或废弃物。机械制造和电子行业在金属加工、表面涂装等环节，可能使用含有PAHs的有机溶剂，这些有机溶剂挥发后会进入大气，最终通过干湿沉降进入闽江。食品加工行业虽然产生的PAHs相对较少，但在一些烟熏、烘焙食品的生产过程中，也会产生一定量的PAHs。这些工业企业的分布和生产规模决定了PAHs的排放强度和空间分布。例如，在工业集中区附近的采样点，沉积物中PAHs的含量往往较高，因为大量的工业排放使得周边环境中的PAHs浓度升高，更容易在沉积物中富集。交通状况也是影响PAHs排放的重要因素。闽江福州段周边交通网络发达，公路、铁路、水路运输繁忙。公路上大量的机动车行驶，其尾气排放是PAHs的重要来源之一。机动车在燃烧汽油或柴油时，由于不完全燃烧，会产生多种PAHs，如苯并[a]芘、荧蒽等。交通流量大的路段，如城市主干道、高速公路出入口附近，PAHs的排放浓度较高。铁路运输中，火车的运行也会产生一定量的PAHs，特别是在柴油机车运行时，其排放的PAHs不容忽视。闽江福州段作为重要的水路运输通道，船舶的航行同样会排放PAHs。船舶发动机燃烧燃料产生的废气中含有PAHs，直接排放到大气中，部分会沉降到江水中，而船舶在维修、保养过程中产生的含油废水等废弃物，若未经妥善处理直接排入江中，也会导致PAHs进入水体和沉积物。农业活动在闽江福州段周边也较为广泛，主要包括种植业和养殖业。在种植业中，农药、化肥的使用可能会引入PAHs。一些有机农药的生产原料或代谢产物中含有PAHs，在使用过程中，这些PAHs会随着雨水冲刷、地表径流等进入闽江。农业废弃物的焚烧也是PAHs的一个来源。农民在田间地头焚烧秸秆、杂草等废弃物时，不完全燃烧会产生大量的PAHs，这些PAHs会通过大气传输进入闽江流域。在养殖业方面，畜禽养殖场产生的粪便若处理不当，其中的有机物在微生物作用下分解，可能会产生PAHs。养殖场附近的水体和土壤中PAHs含量可能会升高，进而通过地表径流等途径影响闽江福州段。闽江福州段周边的人口分布、工业布局、交通状况和农业活动等人类活动，通过不同途径导致了PAHs的排放，这些排放源相互交织，共同影响着闽江福州段沉积物中PAHs的分布、来源和生态风险，是研究闽江福州段PAHs污染问题不可忽视的重要因素。三、多环芳烃分布特征3.1样品采集与分析在闽江福州段的采样工作中，为了全面、准确地获取沉积物中多环芳烃（PAHs）的信息，依据该区域的地理环境、水文条件以及周边人类活动情况，精心规划了采样点的设置。总共在闽江福州段沿线设置了[X]个采样点，涵盖了河流的上、中、下游不同河段。其中，上游区域设置了[X1]个采样点，主要分布在闽清县境内的雄江镇等区域，该区域受人类活动干扰相对较小，设置采样点可获取相对原始状态下沉积物中PAHs的本底值信息。中游区域设置了[X2]个采样点，重点布设在福州市区周边，包括北港和南港流经的台江区、仓山区等区域，这里人口密集、工业活动频繁、交通繁忙，是PAHs的主要排放源所在，通过对这些区域的采样分析，能够了解人类活动对PAHs分布的影响。下游区域设置了[X3]个采样点，集中在马尾区、连江县琯头镇等靠近河口的位置，该区域受潮汐影响明显，通过采样可以研究潮汐作用对PAHs分布的影响。此外，还在一些主要支流与闽江的交汇处设置了[X4]个采样点，以分析支流输入对闽江福州段PAHs分布的贡献。在样品采集方法上，使用专业的不锈钢抓斗式采泥器进行沉积物样品的采集。这种采泥器具有结构坚固、采样效率高、能够有效避免样品受到外界污染等优点。在采集过程中，将采泥器缓慢放入水中，直至到达河底，然后通过机械装置抓取表层0-20cm的沉积物。之所以选择采集表层0-20cm的沉积物，是因为这部分沉积物与水体直接接触，能够更及时、准确地反映水体中PAHs的污染状况，且PAHs在这一深度范围内的分布变化相对较为明显，有利于研究其分布特征。每个采样点采集3个子样，将这3个子样混合均匀后作为该采样点的代表样品。这样做可以减少采样误差，提高样品的代表性，确保分析结果能够真实反映该采样点的PAHs污染情况。在采集过程中，详细记录采样点的地理位置信息，使用高精度的GPS定位仪准确记录采样点的经纬度，精确到小数点后六位。同时，记录周边环境信息，包括是否靠近工厂、居民区、交通要道等，以及采样时间、天气状况、河流流速等现场数据。这些信息对于后续分析PAHs的来源和分布影响因素至关重要。样品采集完成后，立即将其装入经预先清洗和烘干处理的棕色玻璃瓶中，密封保存。选择棕色玻璃瓶是因为其能够有效阻挡光线，减少PAHs在光照条件下可能发生的光化学反应，从而保证样品的稳定性。密封保存则可以防止样品与外界空气接触，避免PAHs的挥发和外界污染物的侵入。尽快将样品运回实验室，并放入冰箱中在4℃下冷藏保存。低温环境可以降低样品中微生物的活性，减缓PAHs的生物降解过程，防止样品中PAHs的含量发生变化，确保后续分析结果的准确性。在实验室中，对沉积物样品进行分析测试。首先将采集的沉积物样品在通风橱中自然风干，去除水分。通风橱能够提供良好的通风条件，加速水分的蒸发，同时避免样品受到室内空气中污染物的污染。自然风干的方式相较于加热烘干等方式，能够更好地保留样品中PAHs的原始形态和含量，减少因高温等因素导致的PAHs损失或变化。风干后的样品用玛瑙研钵研磨至粉末状，过100目筛备用。玛瑙研钵质地坚硬、化学性质稳定，在研磨过程中不会引入杂质，能够保证样品的纯净性。过100目筛可以使样品颗粒更加均匀，有利于后续的提取和分析工作，确保分析结果的准确性和可靠性。采用加速溶剂萃取（ASE）技术提取沉积物中的PAHs。ASE技术具有萃取效率高、速度快、溶剂用量少等优点，能够有效提高PAHs的提取率。具体操作如下：称取适量研磨后的样品，加入适量的硅藻土混合均匀，装入萃取池中。硅藻土具有较大的比表面积和良好的吸附性能，能够帮助分散样品，提高萃取效率。以正己烷-丙酮（体积比为1:1）为萃取剂，在10.3-13.8MPa的压力和100-120℃的温度下进行萃取。这样的压力和温度条件能够使萃取剂更好地渗透到样品内部，与PAHs充分接触，从而提高萃取效果。萃取液经旋转蒸发浓缩后，用硅胶柱进行净化处理。旋转蒸发浓缩可以在较低的温度下快速去除萃取液中的溶剂，减少PAHs的损失。硅胶柱净化能够有效去除萃取液中的杂质，提高PAHs的纯度，为后续的分析测试提供高质量的样品。将浓缩后的萃取液转移至硅胶柱顶部，依次用正己烷、正己烷-二氯甲烷（体积比为1:1）等洗脱剂进行洗脱，收集含有PAHs的洗脱液。通过不同洗脱剂的洗脱，可以将PAHs与其他杂质分离，确保收集到的洗脱液中主要为PAHs。将洗脱液再次浓缩至近干，用正己烷定容至1mL，转移至进样瓶中，供气相色谱-质谱联用仪（GC-MS）分析。使用GC-MS对PAHs进行定性和定量分析。采用HP-5MS毛细管色谱柱（30m×0.25mm×0.25μm），该色谱柱具有良好的分离性能，能够有效分离16种美国环保署（USEPA）优先控制的PAHs单体。进样口温度为280℃，分流比为10:1，进样量为1μL。这样的进样条件能够保证样品在进样口迅速气化，并以合适的比例进入色谱柱进行分离。程序升温条件为：初始温度为60℃，保持1min，以15℃/min的速率升温至280℃，保持10min。通过这样的程序升温，可以使不同沸点的PAHs在色谱柱中得到充分分离，提高分析的准确性。质谱采用电子轰击离子源（EI），离子源温度为230℃，扫描方式为选择离子扫描（SIM）。EI源具有电离效率高、离子碎片丰富等优点，能够提供PAHs的特征质谱信息，便于定性分析。SIM扫描方式则可以提高检测的灵敏度和选择性，确保对痕量PAHs的准确检测。通过与标准物质的保留时间和质谱图对比进行定性分析，外标法进行定量分析。外标法具有操作简单、准确性高的优点，能够准确测定样品中PAHs的含量。在整个分析测试过程中，采取了严格的质量控制措施。每批样品分析时，均同时分析空白样品，空白样品的分析结果应低于方法检出限，以确保实验过程中没有引入外来的PAHs污染。定期对GC-MS进行校准，使用标准PAHs混合溶液进行校准曲线的绘制，校准曲线的相关系数应大于0.995，以保证仪器的准确性和可靠性。对部分样品进行平行样分析，平行样分析结果的相对偏差应小于10%，以评估分析测试结果的精密度。通过这些质量控制措施，确保了分析测试结果的准确性、可靠性和可比性。3.2多环芳烃总量分布对闽江福州段沉积物样品的分析结果显示，该区域沉积物中多环芳烃（PAHs）总量呈现出明显的空间分布特征。PAHs总量的变化范围在[最小值]ng/g至[最大值]ng/g之间，平均值为[平均值]ng/g。从上游到下游来看，总体上呈现出逐渐下降的趋势。在上游区域，PAHs总量相对较高，部分采样点的PAHs总量超过了[上游高值]ng/g。这主要是因为上游地区虽然受人类活动干扰相对较小，但一些工业企业和农业活动仍然会产生一定量的PAHs排放。例如，闽清县境内的一些小型工业企业，在生产过程中可能会排放含有PAHs的废气和废水，这些污染物通过大气沉降和地表径流等方式进入闽江，进而在沉积物中积累。此外，上游地区的农业活动中，农药、化肥的使用以及农业废弃物的焚烧等，也可能导致PAHs的产生和输入。随着河流向下游流动，PAHs总量逐渐降低。在中游区域，PAHs总量平均值为[中游平均值]ng/g，相较于上游有所下降。这一方面是由于河流的稀释作用，随着河水流量的增加，PAHs在水体中的浓度被稀释，进而在沉积物中的积累量也相应减少。另一方面，中游地区的一些自然净化过程，如微生物的降解作用、沉积物的吸附-解吸平衡等，也有助于降低PAHs的含量。然而，在福州市区的北港段，PAHs总量出现了突增的现象。部分采样点的PAHs总量高达[北港高值]ng/g，远高于中游其他区域。这主要是因为福州市区人口密集、工业活动频繁、交通繁忙，是PAHs的主要排放源所在。大量的工业废水排放、城市生活污水排放、机动车尾气排放以及工业废气中的PAHs通过干湿沉降等方式进入闽江，导致北港段沉积物中PAHs的含量显著增加。例如，台江区和仓山区等区域，分布着众多的工业企业和居民点，这些区域的PAHs排放量大，且北港水流相对较缓，有利于PAHs在沉积物中的沉积和积累。到了下游区域，PAHs总量进一步降低，平均值为[下游平均值]ng/g。下游靠近河口，受潮汐影响明显，潮汐的涨落会带来海水的稀释作用，同时也会促进水体的混合和交换，使得PAHs更容易扩散到海洋中，从而减少了在沉积物中的积累。此外，下游地区的生态环境相对较好，植被覆盖度较高，一些植物对PAHs具有一定的吸收和降解能力，也有助于降低PAHs的含量。然而，在一些支流与闽江的交汇处，PAHs总量仍然相对较高。这是因为支流可能携带了来自周边地区的污染物，当支流汇入闽江时，会导致交汇处沉积物中PAHs的含量升高。例如，一些支流周边存在小型加工厂或农业集中区，其排放的污染物会通过支流入江，影响闽江下游沉积物中PAHs的分布。为了更直观地展示PAHs总量的空间分布特征，运用地理信息系统（GIS）技术，结合克里金插值法，绘制了PAHs总量的空间分布图（图1）。从图中可以清晰地看出，PAHs总量高值区域主要集中在上游的部分区域以及福州市区的北港段，而下游大部分区域PAHs总量较低。这种空间分布特征与闽江福州段的自然地理特征、社会经济活动以及水文条件等因素密切相关。通过对PAHs总量分布特征的研究，为进一步解析其来源和评估生态风险提供了重要的基础。3.3单体多环芳烃分布对闽江福州段沉积物中16种美国环保署（USEPA）优先控制的多环芳烃（PAHs）单体含量进行分析，结果显示不同单体的含量存在明显差异，且具有各自独特的分布特征。在这16种PAHs单体中，萘（NAP）、菲（PHE）、荧蒽（FLT）和芘（PYR）的含量相对较高。NAP含量的变化范围在[NAP最小值]ng/g至[NAP最大值]ng/g之间，平均值为[NAP平均值]ng/g。NAP是一种低分子量的PAHs，具有相对较高的挥发性。其在沉积物中的含量相对较高，可能与周边的工业活动和生活污水排放有关。一些工业生产过程，如炼焦、石油化工等，会产生大量的NAP并排放到环境中。生活污水中也可能含有NAP，如一些家庭清洁用品、个人护理产品等在使用过程中可能会释放出NAP。这些来源的NAP通过地表径流、大气沉降等方式进入闽江，进而在沉积物中积累。PHE含量在[PHE最小值]ng/g至[PHE最大值]ng/g之间，平均值为[PHE平均值]ng/g。PHE是一种常见的PAHs，其来源广泛，包括化石燃料的燃烧、生物质燃烧以及石油泄漏等。在闽江福州段，交通尾气排放和工业废气排放可能是PHE的主要来源。大量机动车在道路上行驶，其尾气中含有PHE等PAHs。工业企业在生产过程中，如燃煤发电、钢铁冶炼等，也会排放含有PHE的废气。这些废气中的PHE通过干湿沉降进入闽江，最终在沉积物中富集。FLT含量的变化范围为[FLT最小值]ng/g至[FLT最大值]ng/g，平均值达到[FLT平均值]ng/g。FLT主要来源于高温燃烧过程，如汽车发动机的燃烧、工业锅炉的燃烧等。在闽江福州段，交通繁忙区域和工业集中区附近的沉积物中FLT含量通常较高。例如，福州市区的主要交通干道周边以及工业园区内的采样点，FLT含量明显高于其他区域。这是因为这些区域的机动车尾气排放和工业废气排放量大，使得FLT大量进入环境并在沉积物中积累。PYR含量在[PYR最小值]ng/g至[PYR最大值]ng/g之间，平均值为[PYR平均值]ng/g。PYR与FLT的来源相似，主要也是由燃烧过程产生。在闽江福州段，其分布也受到交通和工业活动的影响。此外，一些垃圾焚烧厂在焚烧垃圾时也会产生PYR，这些PYR通过大气传输进入闽江流域，对沉积物中PYR的含量产生贡献。相比之下，苊烯（ANY）、苊（ANA）、芴（FLU）、蒽（ANT）、苯并[a]蒽（BaA）、䓛（CHR）、苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IPY）、二苯并[a,h]蒽（DBA）和苯并[ghi]苝（BPE）等单体的含量相对较低。其中，BaP作为一种强致癌性的PAHs，其含量在[BaP最小值]ng/g至[BaP最大值]ng/g之间，平均值为[BaP平均值]ng/g。虽然BaP的含量相对较低，但由于其对生态环境和人类健康的严重危害，仍然受到广泛关注。BaP主要来源于高温燃烧过程，如煤炭、石油等化石燃料的不完全燃烧。在闽江福州段，工业锅炉、炉灶等的燃烧以及机动车尾气排放可能是BaP的主要来源。由于其毒性强，即使含量较低，也可能对生态系统和人类健康造成潜在威胁。从空间分布来看，不同PAHs单体在闽江福州段的分布也存在差异。在福州市区的北港段，由于受到工业活动、交通尾气排放和城市生活污水排放等多种污染源的影响，多种PAHs单体的含量都相对较高。例如，NAP、PHE、FLT和PYR等单体在北港段的一些采样点含量明显高于其他区域。而在上游和下游的一些相对偏远区域，PAHs单体的含量相对较低。这是因为这些区域受人类活动干扰较小，污染源相对较少。然而，在一些支流与闽江的交汇处，部分PAHs单体的含量也会出现升高的情况。这是由于支流可能携带了来自周边地区的污染物，当支流汇入闽江时，导致交汇处沉积物中PAHs单体的含量增加。为了进一步分析不同PAHs单体之间的相关性，对各单体含量进行了皮尔逊相关性分析。结果表明，大部分PAHs单体之间呈现显著的正相关关系（P<0.05）。例如，FLT与PYR之间的相关系数达到了[FLT与PYR相关系数]，表明它们可能具有相似的来源和环境行为。这与它们主要都来源于燃烧过程的结论相符合。PHE与ANT之间的相关系数为[PHE与ANT相关系数]，也呈现出较强的相关性。这种相关性分析有助于深入了解PAHs单体之间的内在联系，为源解析和生态风险评价提供更丰富的信息。通过对单体多环芳烃分布的研究，可以更全面地了解闽江福州段沉积物中PAHs的污染特征，为后续的来源解析和生态风险评价提供重要依据。3.4与其他地区对比将闽江福州段沉积物中多环芳烃（PAHs）的含量和分布特征与国内外其他河流、湖泊等水体进行对比，有助于更全面、准确地评估闽江福州段的污染水平，为制定科学合理的污染治理和环境保护措施提供参考依据。与国内一些主要河流相比，闽江福州段沉积物中PAHs的含量处于中等水平。长江作为我国第一大河，其部分河段沉积物中PAHs总量的平均值可达[长江PAHs总量平均值]ng/g，如在一些工业活动密集的区域，PAHs含量甚至超过[长江高值区域PAHs总量]ng/g。这主要是因为长江流域工业发达，众多大型工业企业分布在其沿线，工业废水、废气排放量大，同时交通航运繁忙，船舶尾气排放以及石油类污染物的泄漏等，都导致了PAHs在沉积物中的大量积累。黄河沉积物中PAHs总量平均值约为[黄河PAHs总量平均值]ng/g，其污染程度也较为严重，这与黄河流域的能源结构以煤炭为主，煤炭燃烧过程中产生大量PAHs，以及流域内农业活动中农药、化肥的使用等因素有关。珠江沉积物中PAHs总量平均值在[珠江PAHs总量平均值]ng/g左右，珠江三角洲地区经济快速发展，城市化进程加快，大量的工业活动、交通运输以及城市生活污水排放等，使得珠江沉积物中PAHs含量较高。相比之下，闽江福州段沉积物中PAHs总量平均值为[闽江福州段PAHs总量平均值]ng/g，虽然低于长江、黄河等河流中污染较为严重的区域，但也不容忽视，需要密切关注其污染变化趋势。在与国内湖泊的对比中，太湖沉积物中PAHs总量平均值为[太湖PAHs总量平均值]ng/g，太湖周边工业企业众多，特别是化工、印染等行业，废水排放中含有大量的PAHs。同时，太湖作为重要的旅游胜地，旅游业的发展也带来了一定的污染，如游船尾气排放、游客生活垃圾等，都对太湖沉积物中PAHs的积累产生影响。滇池沉积物中PAHs总量平均值高达[滇池PAHs总量平均值]ng/g，滇池周边城市人口密集，生活污水排放量大，且污水处理设施相对滞后，导致大量含有PAHs的污水直接排入滇池。此外，滇池流域内农业面源污染严重，农药、化肥的不合理使用以及农业废弃物的随意丢弃和焚烧，也增加了PAHs的输入。与这些湖泊相比，闽江福州段沉积物中PAHs的污染程度相对较低，但由于河流的流动性和开放性，PAHs更容易扩散和传播，可能对下游水体和生态系统造成潜在威胁。从国际上看，一些发达国家的河流和湖泊在经过长期的污染治理后，沉积物中PAHs含量相对较低。例如，美国的密西西比河，其沉积物中PAHs总量平均值为[密西西比河PAHs总量平均值]ng/g，美国在环境保护方面投入了大量的资金和技术，制定了严格的环境法规和排放标准，对工业污染源和生活污染源进行了有效的管控，使得密西西比河的水质得到了明显改善，沉积物中PAHs含量处于较低水平。欧洲的莱茵河，在经历了20世纪中叶严重的污染之后，通过实施一系列的污染治理措施，如加强工业废水处理、改善城市污水处理设施、减少农业面源污染等，沉积物中PAHs总量平均值降至[莱茵河PAHs总量平均值]ng/g左右。与之相比，闽江福州段沉积物中PAHs含量相对较高，这反映出我国在环境保护和污染治理方面仍有较大的提升空间，需要借鉴国际先进经验，加强环境管理和污染治理力度。在单体PAHs的对比方面，闽江福州段沉积物中萘（NAP）、菲（PHE）、荧蒽（FLT）和芘（PYR）等单体的含量与其他地区也存在一定差异。在一些工业污染严重的河流中，如我国的海河，NAP含量可高达[海河NAP含量高值]ng/g，远高于闽江福州段的[闽江福州段NAP含量范围]ng/g，这主要是因为海河周边的化工、焦化等行业排放大量含NAP的污染物。在一些城市内河，如上海的苏州河，PHE含量平均值可达[苏州河PHE含量平均值]ng/g，高于闽江福州段的[闽江福州段PHE含量平均值]ng/g，苏州河作为城市内河，受到城市生活污水、工业废水以及交通尾气等多种污染源的影响，PHE的输入量较大。而在一些受石油污染影响较大的水体中，如墨西哥湾因石油泄漏事件后的海域，部分PAHs单体的含量与闽江福州段也有明显不同，其某些与石油源相关的PAHs单体含量显著高于闽江福州段。通过与国内外其他地区的对比可以看出，闽江福州段沉积物中PAHs的污染虽然处于中等水平，但随着福州市社会经济的快速发展，若不加强对PAHs污染源的管控，其污染程度可能会进一步加剧。应充分借鉴国内外先进的污染治理经验，加强环境监测和监管力度，制定严格的环境标准和政策，减少PAHs的排放，保护闽江福州段的生态环境。四、多环芳烃来源解析4.1特征比值法特征比值法是一种常用的多环芳烃（PAHs）来源解析方法，其原理基于不同来源的PAHs在组成和含量上存在差异，导致某些PAHs单体之间的比值具有特定的指示意义。通过分析闽江福州段沉积物中PAHs单体之间的特征比值，可以初步判断PAHs的可能来源。荧蒽（FLT）与芘（PYR）的比值是常用的特征比值之一。一般认为，当FLT/（FLT+PYR）比值小于0.4时，指示PAHs主要来源于石油源。这是因为在石油类物质中，芘的相对含量较高，而荧蒽含量相对较低。当比值在0.4-0.5之间时，表明PAHs为石油源和燃烧源的混合来源。比值大于0.5时，则主要来源于燃烧源。在闽江福州段沉积物样品中，计算得到的FLT/（FLT+PYR）比值范围为[最小值]至[最大值]，平均值为[平均值]。大部分样品的该比值大于0.5，表明闽江福州段沉积物中PAHs主要来源于燃烧源。结合闽江福州段的实际情况，交通尾气排放和工业活动中的高温燃烧过程可能是主要的燃烧源。大量机动车在道路上行驶，其尾气中含有通过燃油不完全燃烧产生的PAHs，这些PAHs通过大气沉降等方式进入闽江。闽江福州段沿线分布着众多工业企业，在工业生产过程中，如锅炉燃烧、金属冶炼等高温燃烧环节，也会产生大量的PAHs并排放到环境中。菲（PHE）与蒽（ANT）的比值也是判断PAHs来源的重要指标。当PHE/ANT比值小于10时，通常认为PAHs主要来源于燃烧源，如生物质燃烧、化石燃料燃烧等。因为在燃烧过程中，蒽的生成量相对较多。当比值大于10时，则可能来源于石油源或石油产品的泄漏等。闽江福州段沉积物样品中PHE/ANT比值在[最小值]至[最大值]之间，平均值为[平均值]。大部分样品的该比值小于10，进一步证实了燃烧源是闽江福州段沉积物中PAHs的主要来源。在福州市区，大量的机动车尾气排放以及周边农村地区的生物质燃烧，如秸秆焚烧等，都可能导致燃烧源PAHs的输入。苯并[a]蒽（BaA）与䓛（CHR）的比值同样具有指示作用。当BaA/（BaA+CHR）比值大于0.35时，暗示PAHs主要来源于交通源，如机动车尾气排放。因为机动车尾气中含有较高比例的苯并[a]蒽。在闽江福州段沉积物中，BaA/（BaA+CHR）比值范围为[最小值]至[最大值]，平均值为[平均值]。部分样品的该比值大于0.35，表明交通源对闽江福州段沉积物中PAHs有一定的贡献。特别是在交通繁忙的区域，如城市主干道、高速公路出入口等附近的采样点，PAHs受交通源的影响更为明显。茚并[1,2,3-cd]芘（IPY）与苯并[ghi]苝（BPE）的比值也可用于源解析。当IPY/（IPY+BPE）比值大于0.5时，指示PAHs主要来源于高温燃烧过程，如工业锅炉的燃烧、垃圾焚烧等。闽江福州段沉积物中IPY/（IPY+BPE）比值在[最小值]至[最大值]之间，平均值为[平均值]。部分样品的该比值大于0.5，说明高温燃烧过程也是闽江福州段沉积物中PAHs的一个重要来源。例如，一些工业企业的锅炉燃烧以及城市垃圾焚烧厂的运行，都会产生通过高温燃烧生成的PAHs。通过对这些特征比值的分析，可以初步判断闽江福州段沉积物中PAHs主要来源于燃烧源，其中交通尾气排放和工业活动中的高温燃烧过程是主要的贡献源。然而，特征比值法存在一定的局限性，它只能提供定性或半定量的信息，且受环境因素的影响较大。因此，为了更准确地解析PAHs的来源，还需要结合其他源解析方法，如主成分分析（PCA）和多元线性回归（MLR）等进行综合分析。4.2因子分析/多元线性回归法为了更准确地定量解析闽江福州段沉积物中多环芳烃（PAHs）的来源，运用因子分析和多元线性回归方法对数据进行深入分析。首先，使用统计分析软件（如SPSS）对闽江福州段沉积物样品中16种PAHs单体的含量数据进行因子分析。因子分析的目的是通过降维的方式，将多个相关变量转化为少数几个不相关的综合变量，即公共因子，这些公共因子能够反映原始变量的主要信息。在进行因子分析前，对数据进行了KMO（Kaiser-Meyer-Olkin）检验和Bartlett球形检验。KMO检验用于衡量变量间的偏相关性，KMO值越接近1，表明变量间的相关性越强，越适合进行因子分析。本研究中KMO值为[KMO值]，大于0.5，说明变量间的相关性较好，适合进行因子分析。Bartlett球形检验用于检验相关矩阵是否为单位矩阵，若检验结果显著（P<0.05），则表明相关矩阵不是单位矩阵，变量间存在相关性，适合进行因子分析。本研究中Bartlett球形检验的P值为[P值]，小于0.05，进一步证实了数据适合进行因子分析。通过因子分析，提取了[因子个数]个公共因子，这[因子个数]个公共因子的累计方差贡献率达到了[累计方差贡献率]%，表明它们能够解释原始数据中大部分的信息。对每个公共因子中各PAHs单体的载荷系数进行分析，确定每个公共因子所代表的潜在污染源。公共因子1中，萘（NAP）、菲（PHE）、蒽（ANT）等低环数PAHs单体具有较高的载荷系数。一般来说，低环数PAHs主要来源于石油类物质的挥发和不完全燃烧。结合闽江福州段的实际情况，周边的工业活动中，如石油化工、机械制造等行业，在生产过程中可能会使用石油类原料，这些原料的挥发和不完全燃烧会产生低环数PAHs。交通运输方面，机动车使用的汽油、柴油等燃料在燃烧过程中也会产生低环数PAHs。因此，公共因子1可能主要代表石油源和交通源的混合来源。公共因子2中，荧蒽（FLT）、芘（PYR）、苯并[a]蒽（BaA）、苯并[b]荧蒽（BbF）等PAHs单体的载荷系数较高。这些PAHs主要来源于高温燃烧过程，如工业锅炉的燃烧、垃圾焚烧以及机动车尾气排放中的高温燃烧部分。在闽江福州段，工业企业的生产活动中，大量的工业锅炉燃烧煤炭、天然气等燃料，会产生这些高温燃烧来源的PAHs。城市中的垃圾焚烧厂在处理垃圾时，也会排放含有这些PAHs的废气。此外，交通繁忙区域的机动车尾气排放也是重要的来源之一。所以，公共因子2可能主要代表高温燃烧源，包括工业燃烧和交通尾气排放中的高温燃烧部分。公共因子3中，茚并[1,2,3-cd]芘（IPY）、二苯并[a,h]蒽（DBA）、苯并[ghi]苝（BPE）等高分子量PAHs单体具有较高的载荷系数。这些高分子量PAHs通常在高温、高压的环境下形成，主要来源于煤炭、生物质等的燃烧。在闽江福州段周边的农村地区，居民生活中使用煤炭取暖、烹饪，以及农业废弃物的焚烧等活动，都可能产生这些高分子量PAHs。因此，公共因子3可能主要代表煤炭和生物质燃烧源。在确定了各公共因子所代表的潜在污染源后，以各公共因子得分作为自变量，以PAHs总量为因变量，进行多元线性回归分析。建立多元线性回归方程：Y=a_1X_1+a_2X_2+a_3X_3+\cdots+a_nX_n+b，其中Y为PAHs总量，X_i为第i个公共因子得分，a_i为第i个公共因子的回归系数，b为常数项。通过回归分析，得到各公共因子的回归系数，并计算出各潜在污染源对PAHs总量的贡献率。计算结果表明，石油源和交通源的混合来源对闽江福州段沉积物中PAHs总量的贡献率为[贡献率1]%。这说明石油类物质的挥发、不完全燃烧以及交通源排放对PAHs的积累有重要影响。高温燃烧源对PAHs总量的贡献率为[贡献率2]%，表明工业燃烧和交通尾气排放中的高温燃烧部分是PAHs的主要来源之一。煤炭和生物质燃烧源对PAHs总量的贡献率为[贡献率3]%，虽然相对前两者贡献率较低，但在部分区域，如农村地区，仍然是不可忽视的PAHs来源。通过因子分析和多元线性回归法，不仅能够识别闽江福州段沉积物中PAHs的主要来源，还能定量计算各来源的贡献率，为制定针对性的污染控制措施提供了更准确的科学依据。与特征比值法相比，该方法能够更全面、深入地解析PAHs的来源，克服了特征比值法只能提供定性或半定量信息的局限性。然而，该方法也存在一定的局限性，如因子分析的结果可能受到数据质量、样本数量等因素的影响，在实际应用中需要结合研究区域的实际情况进行综合分析和验证。4.3来源综合分析综合特征比值法和因子分析/多元线性回归法的结果，可以更全面、准确地剖析闽江福州段沉积物中多环芳烃（PAHs）的主要来源及其分布特征。从特征比值法来看，荧蒽（FLT）与芘（PYR）比值、菲（PHE）与蒽（ANT）比值、苯并[a]蒽（BaA）与䓛（CHR）比值以及茚并[1,2,3-cd]芘（IPY）与苯并[ghi]苝（BPE）比值等的分析结果均表明，燃烧源是闽江福州段沉积物中PAHs的主要来源。其中，交通尾气排放和工业活动中的高温燃烧过程在燃烧源中占据重要地位。而因子分析/多元线性回归法进一步定量地揭示了各污染源的贡献率。石油源和交通源的混合来源对PAHs总量的贡献率为[贡献率1]%，高温燃烧源（包括工业燃烧和交通尾气排放中的高温燃烧部分）的贡献率为[贡献率2]%，煤炭和生物质燃烧源的贡献率为[贡献率3]%。在空间分布上，不同来源的PAHs呈现出各自的特点。在福州市区，尤其是北港段，由于人口密集、工业活动频繁以及交通繁忙，交通源和工业源的PAHs排放显著。这里分布着众多的工业企业，如化工、机械制造等行业，在生产过程中会排放大量含有PAHs的废气、废水。同时，大量机动车在道路上行驶，尾气排放也为该区域带来了丰富的PAHs来源。因此，在北港段的沉积物中，来自交通源和工业源的PAHs含量较高，这与特征比值法中部分样品BaA/（BaA+CHR）比值大于0.35，指示交通源贡献明显，以及因子分析中公共因子2代表的高温燃烧源（包括工业燃烧和交通尾气排放中的高温燃烧部分）贡献率较高的结果相契合。在上游和下游的一些相对偏远区域，受人类活动干扰较小，污染源相对较少。然而，上游地区的一些小型工业企业和农业活动，如农药、化肥的使用以及农业废弃物的焚烧等，仍然会产生一定量的PAHs。这些PAHs主要来源于石油源和煤炭、生物质燃烧源，与因子分析中公共因子1和公共因子3所代表的污染源相对应。下游靠近河口，受潮汐影响明显，虽然总体PAHs含量较低，但在一些支流与闽江的交汇处，由于支流可能携带了来自周边地区的污染物，使得该区域沉积物中PAHs含量升高。这些污染物的来源可能较为复杂，既包括上游地区的工业和农业排放，也可能有支流周边小型加工厂或农业集中区的贡献。闽江福州段沉积物中PAHs的来源主要为燃烧源，其中交通尾气排放、工业活动中的高温燃烧过程以及煤炭和生物质燃烧是主要的贡献源。不同来源的PAHs在空间上呈现出与人类活动和地理环境相关的分布特征。通过两种源解析方法的综合运用，为深入了解闽江福州段PAHs的污染状况提供了更全面的信息，也为制定针对性的污染控制措施提供了坚实的科学依据。在未来的环境保护工作中，应重点针对交通源和工业源进行管控，减少PAHs的排放，同时关注煤炭和生物质燃烧源在部分区域的影响，加强对农业活动的管理，以降低闽江福州段沉积物中PAHs的污染水平，保护闽江的生态环境。五、多环芳烃生态风险评价5.1评价方法选择多环芳烃（PAHs）因其具有毒性、持久性和生物累积性，对生态环境和人类健康存在潜在威胁，因此对闽江福州段沉积物中PAHs进行生态风险评价至关重要。目前，常用的PAHs生态风险评价方法有多种，每种方法都有其特点和适用范围。效应区间低值（ERL）和效应区间中值（ERM）法是较为常用的评价方法之一。该方法由美国国家海洋和大气管理局（NOAA）提出，通过将沉积物中PAHs的浓度与ERL和ERM值进行对比，来判断PAHs对生物产生负面效应的可能性。当PAHs浓度低于ERL时，被认为对生物产生负面影响的可能性较低；当浓度介于ERL和ERM之间时，具有一定的潜在风险；当超过ERM时，则存在较高的生态风险。这种方法简单直观，能够快速地对PAHs的生态风险进行初步评估，且其评价标准是基于大量的生物毒性试验数据得出的，具有较高的可靠性。在对美国多个河口和近岸海域沉积物中PAHs的生态风险评价中，该方法被广泛应用，并为当地的环境管理和污染治理提供了重要依据。毒性当量因子（TEF）法也是常用的评价方法。由于不同PAHs单体的毒性差异较大，TEF法通过将不同PAHs单体的含量换算成具有相同毒性的苯并[a]芘（BaP）的当量浓度，即毒性当量浓度（TEQ），来综合评估PAHs的总体毒性和生态风险。计算公式为：TEQ=\sum_{i=1}^{n}(C_{i}\timesTEF_{i})，其中C_{i}为第i种PAHs单体的浓度，TEF_{i}为第i种PAHs单体的毒性当量因子。这种方法考虑了不同PAHs单体的毒性差异，能够更准确地反映PAHs的综合毒性效应。在对一些工业污染场地和河流沉积物的生态风险评价中，TEF法被用于评估PAHs对生态系统和人类健康的潜在危害，为污染治理和风险管控提供了科学参考。有机碳归一化法是基于PAHs具有亲脂性，易吸附在沉积物中的有机质上这一特性。该方法通过将PAHs含量与沉积物中的有机碳含量进行归一化处理，来消除不同沉积物样品中有机碳含量差异对PAHs含量的影响，从而更准确地评估PAHs的生态风险。研究表明，在一些河流和湖泊的沉积物中，PAHs含量与有机碳含量存在显著的正相关关系，采用有机碳归一化法能够更合理地评价PAHs在不同沉积物样品中的相对污染程度和生态风险。概率风险评价方法则考虑了PAHs含量的不确定性和生物对PAHs毒性响应的不确定性。通过蒙特卡罗模拟等方法，设定PAHs含量和毒性响应参数的概率分布函数，进行多次模拟计算，得到不同风险水平下的发生概率。这种方法能够提供更全面的风险信息，有助于风险管理决策的制定。在对一些复杂环境系统中PAHs的生态风险评价中，概率风险评价方法被用于评估PAHs对生态系统和人类健康的潜在风险，考虑了多种不确定性因素，为环境管理提供了更科学的依据。在本研究中，选择效应区间低值（ERL）和效应区间中值（ERM）法、毒性当量因子（TEF）法以及概率风险评价方法对闽江福州段沉积物中PAHs进行生态风险评价。选择这些方法的依据主要有以下几点：首先，效应区间低值和效应区间中值法简单直观，能够对PAHs的生态风险进行初步的快速评估，且其评价标准具有广泛的认可度和应用案例，便于与其他地区的研究结果进行对比。其次，毒性当量因子法能够考虑不同PAHs单体的毒性差异，更准确地反映PAHs的综合毒性效应，对于全面评估闽江福州段沉积物中PAHs的生态风险具有重要意义。最后，概率风险评价方法考虑了不确定性因素，能够提供更丰富的风险信息，有助于制定更科学合理的风险管理策略。通过多种方法的综合应用，可以从不同角度对闽江福州段沉积物中PAHs的生态风险进行全面、准确的评价，为该区域的环境保护和污染治理提供更可靠的科学依据。5.2生态风险评价结果运用选定的效应区间低值（ERL）和效应区间中值（ERM）法、毒性当量因子（TEF）法以及概率风险评价方法，对闽江福州段沉积物中多环芳烃（PAHs）进行生态风险评价，结果如下：效应区间低值（ERL）和效应区间中值（ERM）法评价结果：将闽江福州段沉积物中16种PAHs单体的含量与美国国家海洋和大气管理局（NOAA）提出的ERL和ERM值进行对比。部分PAHs单体在某些区域呈现出一定的生态风险。其中，荧蒽（FLT）在部分采样点的含量介于ERL值（19.0ng/g）和ERM值（180ng/g）之间，表明这些区域的FLT具有一定的潜在生态风险。例如，在福州市区北港段的[具体采样点名称]，FLT含量达到了[具体含量]ng/g，处于该风险区间内。芘（PYR）在个别采样点的含量也接近ERL值，存在潜在风险。而其他PAHs单体，如萘（NAP）、菲（PHE）、蒽（ANT）等，大部分采样点的含量均低于ERL值，对生物产生负面影响的可能性较低。从PAHs总量来看，计算得到的ERL商值（PAHs总量/ERL总量）范围为[最小值]至[最大值]，平均值为[平均值]，大部分区域的ERL商值小于1，表明总体上闽江福州段沉积物中PAHs总量对生物产生负面影响的可能性较低。但在一些受污染较为严重的区域，如北港段的部分采样点，ERL商值接近或略大于1，存在一定的潜在风险。ERM商值（PAHs总量/ERM总量）范围为[最小值]至[最大值]，平均值为[平均值]，所有采样点的ERM商值均远小于1，说明闽江福州段沉积物中PAHs总量超过ERM值的情况较少，整体不存在较高的生态风险。毒性当量因子（TEF）法评价结果：根据国际上公认的TEF值，将闽江福州段沉积物中不同PAHs单体的含量换算成苯并[a]芘（BaP）的毒性当量浓度（TEQ）。计算得到的TEQ值范围为[最小值]ng/g至[最大值]ng/g，平均值为[平均值]ng/g。与相关评价标准相比，部分采样点的TEQ值超过了阈值，存在一定的生态风险。例如，在[具体采样点名称]，TEQ值达到了[具体含量]ng/g，超过了[阈值]ng/g，表明该区域的PAHs具有一定的毒性风险。从不同区域来看，福州市区北港段的TEQ值相对较高，这与该区域受到工业活动、交通尾气排放等多种污染源的影响有关，使得具有较高毒性的PAHs单体在该区域的沉积物中积累较多。而在上游和下游的一些相对偏远区域，TEQ值相对较低，生态风险相对较小。概率风险评价结果：利用蒙特卡罗模拟软件（如CrystalBall），考虑PAHs含量的不确定性和生物对PAHs毒性响应的不确定性。设定PAHs含量服从正态分布，毒性响应参数服从对数正态分布等，进行多次模拟计算。结果显示，在低风险水平下（如风险概率低于10%），大部分采样点的PAHs对水生生物产生毒性效应的概率较低。但在较高风险水平下（如风险概率高于50%），部分采样点，尤其是福州市区北港段和一些支流与闽江交汇处的采样点，PAHs对水生生物产生毒性效应的概率明显增加。例如，在北港段的[具体采样点名称]，在高风险水平下，PAHs对水生生物产生毒性效应的概率达到了[具体概率]，表明这些区域的PAHs生态风险不容忽视。通过概率风险评价，能够更全面地了解闽江福州段沉积物中PAHs生态风险的不确定性，为风险管理提供更丰富的信息。综合以上三种评价方法的结果，闽江福州段沉积物中PAHs在部分区域存在一定的生态风险，尤其是福州市区北港段以及一些支流与闽江交汇处。这些区域受到人类活动的影响较大，PAHs的来源复杂，含量相对较高，具有一定的潜在生态风险，需要重点关注和加强监测。而在上游和下游的一些相对偏远区域，PAHs的生态风险相对较小。但由于PAHs具有持久性和生物累积性，即使在生态风险相对较小的区域，也不能忽视其长期潜在影响，应加强对整个闽江福州段沉积物中PAHs的管控和治理，保护闽江的生态环境。5.3风险因素分析闽江福州段沉积物中多环芳烃（PAHs）的生态风险受到多种因素的综合影响，深入分析这些因素对于全面理解PAHs的生态风险状况以及制定有效的污染防控策略具有重要意义。PAHs自身的含量和组成是影响生态风险的关键内在因素。从含量上看，随着沉积物中PAHs总量的增加，其对生态系统产生负面影响的可能性增大。当PAHs含量超过一定阈值时，会对水生生物的生理机能产生直接损害，如抑制水生生物的生长、繁殖，影响其免疫和神经系统功能。在闽江福州段，部分采样点较高的PAHs含量，如福州市区北港段，导致该区域生态风险相对较高。在PAHs组成方面，不同环数和结构的PAHs单体毒性差异显著。一般来说，高环数PAHs，如苯并[a]芘（BaP）、二苯并[a,h]蒽（DBA）等，具有更强的致癌、致畸和致突变性，对生态风险的贡献更大。在闽江福州段沉积物中，虽然这些高毒性PAHs单体含量相对较低，但由于其毒性极强，在生态风险评价中不容忽视。它们的存在增加了闽江福州段生态系统的潜在风险，可能对水生生物和人体健康构成长期威胁。沉积物的性质对PAHs生态风险有着重要影响。其中，有机质含量是一个关键因素。PAHs具有亲脂性，易吸附于沉积物中的有机质上。较高的有机质含量能够为PAHs提供更多的吸附位点，从而增加PAHs在沉积物中的富集量。在闽江福州段，沉积物有机质含量与PAHs含量在部分区域呈现显著正相关。这意味着有机质含量高的区域，PAHs更容易积累，进而增加了生态风险。例如，在一些支流与闽江交汇处，由于周边土壤侵蚀等原因，可能导致沉积物中有机质含量增加，同时也伴随着PAHs含量的升高，生态风险增大。沉积物的粒度组成也会影响PAHs的生态风险。细颗粒沉积物比表面积大，对PAHs的吸附能力更强。在闽江福州段，细颗粒沉积物较多的区域，如河湾、缓流区等，PAHs更容易被吸附和积累，使得这些区域的生态风险相对较高。而在水流湍急、粗颗粒沉积物较多的区域，PAHs不易停留和积累，生态风险相对较低。环境条件也是影响PAHs生态风险的重要因素。水文条件，如河流的流量、流速和水位变化等，对PAHs的迁移、扩散和沉积过程产生重要影响。在洪水期，河流流量和流速增大，可能导致沉积物的再悬浮