珠江口表层沉积物多环芳烃：分布、来源与生态风险的深度剖析

珠江口表层沉积物多环芳烃：分布、来源与生态风险的深度剖析一、引言1.1研究背景与意义多环芳烃（PolycyclicAromaticHydrocarbons，PAHs）是一类由两个或两个以上苯环以稠环形式相连的化合物，其种类繁多，目前已发现的致癌性多环芳烃及其衍生物超过400种。PAHs具有“致癌、致畸、致突变”效应，是一类典型的持久性有机污染物，对人类健康和生态环境构成严重威胁。其来源广泛，主要源于含碳化合物的不完全燃烧，如石油、木材、垃圾和煤等，以及原油在开采、运输、生产和使用过程中的泄漏及排污。在工业生产中，焦化厂、钢铁厂等排放的废弃物中含有大量PAHs；交通运输方面，飞机、汽车等机动车辆排放的废气中也存在相当数量的PAHs，且在启动时由于不完全燃烧排放量更大；日常生活里，家庭锅炉使用、垃圾焚烧、吸烟以及烧烤等活动也会产生PAHs。PAHs通过废水排放、大气沉降、土表迁移、石油泄漏等多种途径进入水体，由于其低溶解性和憎水性，极易吸附于颗粒物上，并聚集到沉积物中，使得沉积物成为PAHs的蓄积库。珠江口作为珠江水系的重要出海口，连接着内陆与海洋，周边分布着众多城市和工业区域，人口密集，经济活动频繁。工业废水、生活污水的排放，以及船舶运输、港口作业等活动，都可能导致大量PAHs进入珠江口的水体和沉积物中。已有研究表明，珠江口表层沉积物中PAHs的含量处于一定水平，且存在多种来源。例如，珠江水产研究所渔业环境保护与修复创新团队的研究发现，珠江三角洲表层沉积物中多环芳烃主要来源于煤焦油和液体燃料燃烧、煤/木材燃烧和石油扩散的混合来源。因此，研究珠江口表层沉积物中PAHs的分布特征，能够了解该区域PAHs的污染程度和空间差异，为评估其污染状况提供基础数据。对珠江口表层沉积物中PAHs进行生态风险评价，有助于明确其对生态系统和人类健康的潜在威胁，为制定针对性的污染防治措施提供科学依据。珠江口水域是众多水生生物的栖息地和繁殖场所，PAHs的污染可能对这些生物的生存、繁殖和生长发育产生负面影响，进而破坏整个生态系统的平衡和稳定。同时，PAHs通过食物链的传递和富集，最终可能进入人体，危害人类健康。因此，深入研究珠江口表层沉积物中PAHs的分布特征与生态风险评价，对于保护珠江口的生态环境、保障人类健康以及促进区域的可持续发展具有重要的现实意义。1.2国内外研究现状多环芳烃在沉积物中的研究一直是环境科学领域的重点。国外对多环芳烃的研究起步较早，在20世纪70年代，美国环保署（EPA）就将16种PAHs列入优先控制的有机污染物黑名单，此后，相关研究不断深入。在分布特征研究方面，Baumard等学者对地中海西部沉积物和贻贝中的多环芳烃进行了研究，详细分析了PAHs在不同介质中的含量和分布情况，发现该区域沉积物中PAHs含量存在明显的空间差异。在来源解析上，国外学者广泛运用多种方法，如Hartmann等人采用同分异构体比值法，对美国纳拉甘西特湾表层沉积物中PAHs的来源进行判别，发现其主要来源于交通、工业活动以及石油泄漏等。在生态风险评价领域，欧盟生态风险评价导则（TGD）根据PAHs对水生生物预测的水相无观察效应环境浓度（PNECw）构建预测的沉积物无观察效应环境浓度（PNECs），以此为基础进行风险表征，为沉积物中PAHs的生态风险评价提供了重要的参考标准。国内对多环芳烃的研究相对较晚，但近年来发展迅速。在分布特征方面，许多研究聚焦于不同水体的沉积物。母清林、方杰、邵君波等对长江口及浙江近岸海域表层沉积物中多环芳烃的分布进行了研究，发现该区域PAHs含量呈现出从河口向近岸海域逐渐降低的趋势。在来源解析方面，国内学者也综合运用多种手段。例如，闽江福州段沉积物中多环芳烃的研究中，采用特征比值/插值法和因子分析/多元线性回归法，判别出其主要来源为石油燃烧和煤燃烧。在生态风险评价上，国内学者运用多种方法进行综合评估。以大治湖表层沉积物-水中多环芳烃的研究为例，采用沉积物质量基准法和风险商法，对其生态风险进行评价，发现该区域部分点位存在一定的生态风险。针对珠江口的研究，也取得了一定成果。珠江水产研究所渔业环境保护与修复创新团队采用索氏提取、气相色谱-质谱联用（GC-MS）等手段，对珠江三角洲表层沉积物中PAHs的时空变化趋势、来源以及潜在风险等方面开展研究，发现表层沉积物中PAHs含量范围在56.6-1297ng/g，处于世界中等水平，主要来源于煤焦油和液体燃料燃烧、煤/木材燃烧和石油扩散的混合来源，且珠江口的生态风险高于上游河网。然而，当前珠江口表层沉积物中PAHs的研究仍存在不足。在时间尺度上，对PAHs的长期动态变化研究较少，难以准确把握其随时间的演变规律。在空间分布上，对一些偏远或特殊区域的研究覆盖不足，导致对整个珠江口PAHs分布的全面认识存在偏差。在来源解析方面，虽然已明确主要来源，但对于各来源的相对贡献在不同季节、不同区域的变化情况，缺乏深入细致的研究。在生态风险评价上，现有的评价方法和指标体系仍有待进一步完善，以更准确地评估PAHs对珠江口生态系统和人类健康的潜在威胁。1.3研究内容与方法本研究旨在全面深入地探究珠江口表层沉积物中多环芳烃的分布特征与生态风险，具体研究内容涵盖以下几个关键方面。含量测定：通过科学合理的采样方法，广泛收集珠江口不同区域的表层沉积物样本。运用先进的前处理技术，如微波辅助萃取、超声波提取、快速溶剂萃取等，对样本中的多环芳烃进行高效提取。随后，采用高灵敏度的分析仪器，如气相色谱-质谱联用仪（GC-MS）、高效液相色谱（HPLC）及其衍生技术如HPLC-UV、LC-MS/MS等，精确测定沉积物中16种美国环保署（EPA）优先控制的多环芳烃（PAHs）的含量。分布特征分析：依据含量测定所获取的数据，运用地理信息系统（GIS）技术，直观清晰地绘制出PAHs在珠江口表层沉积物中的空间分布格局图。深入分析不同区域PAHs含量的差异，结合珠江口的地理环境、水文条件以及周边的人类活动等因素，全面探讨影响PAHs分布的主要因素。例如，研究工业废水排放口、城市生活污水排放点以及船舶运输繁忙区域等与PAHs高含量区域的相关性。来源解析：综合运用多种成熟的源解析方法，如同分异构体比值法，通过分析荧蒽/芘、菲/蒽等特征比值，初步判别PAHs的可能来源；因子分析/多元线性回归法，借助统计分析手段，确定不同来源对PAHs的相对贡献比例；以及定性（半定量）受体模型等，准确追溯珠江口表层沉积物中PAHs的具体来源，明确其主要来源于煤焦油和液体燃料燃烧、煤/木材燃烧还是石油扩散等，以及各来源在不同区域的差异情况。生态风险评价：运用沉积物质量指南（SQG），将测定的PAHs含量与相应的阈值进行对比，初步评估生态风险；平均效应水平熵值（PEL-Q），综合考虑多种PAHs的联合效应，更全面地评价风险程度；毒性等效因子（TEQBaP）等方法，以苯并[a]芘为参照，将其他PAHs转化为等效毒性，准确评估PAHs对珠江口生态系统的潜在风险，确定风险较高的区域和主要的风险贡献物质。在研究方法上，采样工作在珠江口设置多个具有代表性的采样点，涵盖河口、近岸、岛屿周边等不同区域，使用抓斗式采泥器采集表层0-5cm的沉积物样本，每个采样点重复采集3次，以确保样本的代表性和可靠性。样本采集后，立即放入低温冷藏箱中保存，带回实验室后冷冻干燥，去除水分和杂质，过筛后备用。分析方法上，采用微波辅助萃取法对沉积物中的PAHs进行提取，通过优化萃取条件，如萃取时间、温度、溶剂种类和用量等，提高提取效率。提取液经净化处理后，使用气相色谱-质谱联用仪进行分析，通过选择离子监测模式（SIM），对16种PAHs进行定性和定量分析。数据处理方面，运用统计学软件对测定的数据进行统计分析，计算平均值、标准差、变异系数等参数，分析数据的集中趋势和离散程度。采用相关性分析、主成分分析等方法，探讨PAHs含量与其他环境因素之间的关系，筛选出影响PAHs分布和来源的关键因素。利用ArcGIS软件进行空间分析，绘制PAHs的空间分布图和等值线图，直观展示其分布特征。二、珠江口表层沉积物多环芳烃含量及分布特征2.1样品采集与分析为全面、准确地获取珠江口表层沉积物中多环芳烃的相关信息，本研究于[具体年份]的[春季/夏季/秋季/冬季]，在珠江口海域精心设置了[X]个采样点，采样点的分布充分考虑了珠江口的地理环境、水文条件以及周边人类活动的影响。其中，在靠近工业发达的区域如广州、深圳等城市附近设置了多个采样点，以重点监测工业活动对PAHs分布的影响；在河流入海口处，如虎门、蕉门、洪奇门和横门等口门附近也进行了采样，以研究河流径流携带的PAHs在河口区域的分布特征；同时，在远离海岸的开阔海域设置了对照采样点，用于对比分析不同区域PAHs含量的差异。采样点的具体位置通过高精度的全球定位系统（GPS）进行定位，确保每个采样点的位置准确无误，其分布情况如图[X]所示。使用抓斗式采泥器采集表层0-5cm的沉积物样品，该深度的沉积物能够较好地反映近期PAHs的输入情况。在每个采样点，重复采集3次样品，将采集到的样品迅速放入低温冷藏箱中保存，冷藏箱温度控制在4℃左右，以防止样品中PAHs的挥发和降解。带回实验室后，立即对样品进行冷冻干燥处理，冷冻干燥过程在真空环境下进行，温度控制在-50℃左右，持续时间为48小时，以确保样品中的水分完全去除。干燥后的样品使用玛瑙研钵研磨成粉末状，过100目尼龙筛，去除杂质和较大颗粒，得到均匀的样品粉末，过筛后的样品储存于干燥器中备用。在样品处理分析流程中，首先采用索氏提取法对沉积物中的PAHs进行提取。称取5g过筛后的样品粉末，放入索氏提取器的滤纸筒中，加入100mL二氯甲烷作为提取溶剂，在70℃的水浴温度下回流提取24小时，使PAHs充分溶解于二氯甲烷中。提取液经旋转蒸发仪浓缩至约1mL，浓缩过程中温度控制在40℃左右，以避免PAHs的损失。浓缩后的提取液通过硅胶柱进行净化处理。硅胶柱的制备方法为：将硅胶（100-200目）在150℃下活化4小时，然后称取5g活化后的硅胶，加入适量的正己烷，搅拌均匀后装入玻璃柱中，柱高约为10cm。在硅胶柱顶部加入1g无水硫酸钠，以去除提取液中的水分。将浓缩后的提取液缓慢加入硅胶柱中，先用50mL正己烷淋洗硅胶柱，去除脂肪、色素等杂质，再用50mL二氯甲烷和正己烷（体积比为1:1）的混合溶液洗脱PAHs，收集洗脱液。洗脱液经氮吹仪浓缩至约0.5mL，加入内标物（萘-d8、二氢苊-d10等），采用气相色谱-质谱联用仪（GC-MS）进行测定。GC-MS分析条件如下：色谱柱为DB-5MS毛细管柱（30m×0.25mm×0.25μm）；进样口温度为280℃；载气为高纯氦气，流速为1mL/min；进样方式为分流进样，分流比为10:1；程序升温条件为：初始温度为60℃，保持2min，以15℃/min的速率升温至300℃，保持5min。质谱条件为：离子源为电子轰击源（EI），离子源温度为230℃；扫描方式为选择离子监测（SIM），监测离子为16种PAHs的特征离子。在整个分析过程中，严格进行质量控制。每批样品分析时均同时进行空白实验，空白样品的处理过程与实际样品相同，以检测分析过程中是否存在污染。空白样品中16种PAHs的含量均低于仪器的检测限，表明分析过程无污染。同时，进行加标回收实验，向已知含量的样品中加入一定量的PAHs标准溶液，按照上述分析方法进行测定，计算加标回收率。结果显示，16种PAHs的加标回收率在70%-120%之间，满足分析要求，确保了实验数据的准确性和可靠性。2.2多环芳烃含量水平通过对珠江口[X]个采样点的表层沉积物样品进行分析测定，结果显示珠江口表层沉积物中16种美国环保署（EPA）优先控制的多环芳烃（PAHs）总量（∑16PAHs）的含量范围为[最小值]ng/g至[最大值]ng/g，平均值为[平均值]ng/g。其中，含量较高的PAHs单体主要有菲、荧蒽、芘等，其含量范围分别为[菲的最小值]-[菲的最大值]ng/g、[荧蒽的最小值]-[荧蒽的最大值]ng/g、[芘的最小值]-[芘的最大值]ng/g，平均值分别为[菲的平均值]ng/g、[荧蒽的平均值]ng/g、[芘的平均值]ng/g。这些PAHs单体在珠江口表层沉积物中的相对含量较高，可能是由于其来源广泛，且在环境中的稳定性较强。与国内外其他区域的研究结果相比，珠江口表层沉积物中PAHs的含量处于一定水平。在中国国内，与长江三角洲扬子江河口相比，其PAHs总浓度范围为27-622ng/g，珠江口的含量范围相对较高；与香港米埔和后海湾（162-384ng/g）相比，珠江口的含量也处于较高区间。在国际上，与东北大西洋浅滩（5.26-51.1ng/g）和西北大西洋伊比利亚半岛（0.9-94ng/g）相比，珠江口表层沉积物中PAHs的含量明显偏高。然而，与韩国沿海（8.8-18500ng/g）相比，珠江口的含量处于中等水平。通过对比可以发现，珠江口表层沉积物中PAHs的污染程度相对较高，这可能与珠江口周边地区密集的人类活动密切相关。珠江三角洲地区是中国重要的经济区，工业发达，人口密集，大量的工业废水、生活污水排放，以及船舶运输、港口作业等活动，都可能导致PAHs大量进入珠江口的水体和沉积物中。同时，珠江口的地理位置特殊，作为珠江水系的出海口，其承接了上游地区的污染物输入，进一步加剧了PAHs的污染程度。2.3空间分布特征利用地理信息系统（GIS）技术，对珠江口各采样点表层沉积物中PAHs的含量进行空间分析，绘制出PAHs的空间分布格局图，结果如图[X]所示。从图中可以清晰地看出，PAHs在珠江口表层沉积物中的空间分布呈现出明显的差异性。河口区域，如虎门、蕉门、洪奇门和横门等口门附近，PAHs含量相对较高。以虎门为例，该区域采样点的PAHs含量平均值达到[具体数值]ng/g，显著高于珠江口的平均水平。这主要是因为河口地区是陆源污染物进入海洋的重要通道，珠江水系携带了大量来自上游地区的工业废水、生活污水以及农业面源污染，其中包含丰富的PAHs。同时，河口地区的水动力条件复杂，河流径流与海洋潮汐相互作用，使得污染物容易在此处汇聚和沉积。此外，河口周边通常分布着众多的工业企业和城市，工业生产活动如石油化工、钢铁冶炼等会产生大量的PAHs，通过废水排放和大气沉降等途径进入河口沉积物中；城市生活污水的排放以及垃圾填埋场的渗滤液也可能含有PAHs，进一步增加了河口区域的污染负荷。近岸区域，PAHs含量整体也处于较高水平，但相较于河口区域，存在一定的梯度变化。靠近城市和工业集中区的近岸采样点，PAHs含量明显高于远离城市的近岸区域。例如，靠近深圳、广州等城市的近岸采样点，PAHs含量平均值在[具体数值1]-[具体数值2]ng/g之间，而远离城市的近岸采样点，含量平均值则在[具体数值3]-[具体数值4]ng/g之间。这是因为近岸地区受到人类活动的强烈影响，城市的发展和扩张导致人口密集，生活污水和垃圾的产生量大幅增加，其中的PAHs通过地表径流、污水处理厂排水等方式进入近岸海域。工业集中区的存在使得工业废水排放成为近岸PAHs污染的重要来源，一些工厂在生产过程中未能对含有PAHs的废水进行有效处理，直接排放到近岸水体，随后被沉积物吸附。此外，近岸海域的船舶运输活动频繁，船舶燃油的燃烧以及油品泄漏也会向水体中释放PAHs，进而在沉积物中积累。远岸区域，PAHs含量相对较低，整体分布较为均匀。大部分远岸采样点的PAHs含量在[具体数值5]-[具体数值6]ng/g之间，明显低于河口和近岸区域。这是由于远岸地区距离污染源较远，陆源污染物在传输过程中经过稀释和扩散，浓度逐渐降低。同时，远岸海域的水动力条件相对简单，水体交换活跃，有利于污染物的扩散和稀释，减少了PAHs在沉积物中的富集。此外，远岸地区人类活动相对较少，工业生产、生活污水排放以及船舶运输等活动产生的PAHs量也相对较少，从而使得远岸区域的PAHs污染程度较轻。综上所述，珠江口表层沉积物中PAHs的空间分布主要受到陆源输入、水动力条件和人类活动强度等因素的综合影响。陆源输入是导致PAHs在河口和近岸区域高含量分布的主要原因，水动力条件则影响着污染物的传输、扩散和沉积过程，人类活动强度的差异进一步加剧了PAHs在不同区域的分布差异。2.4组成特征对珠江口表层沉积物中不同环数PAHs的相对比例进行分析，结果表明，2-3环PAHs的含量占总PAHs含量的[X1]%，4环PAHs占[X2]%，5-6环PAHs占[X3]%。其中，3环PAHs中的菲，在所有PAHs单体中含量相对较高，平均占总PAHs含量的[具体比例1]，是珠江口表层沉积物中的优势组分之一。4环PAHs中的荧蒽和芘含量也较为突出，分别平均占总PAHs含量的[具体比例2]和[具体比例3]。5-6环PAHs中，苯并[a]芘、二苯并[a,h]蒽等虽然含量相对较低，但由于其具有较强的致癌性，备受关注。不同环数PAHs的相对比例与污染源密切相关。一般来说，2-3环PAHs主要来源于石油类污染和低温燃烧过程，如原油泄漏、石油产品的使用等。在珠江口，部分区域靠近石油化工企业和港口，这些区域的沉积物中2-3环PAHs的相对比例可能较高。例如，在靠近惠州大亚湾石化区的采样点，2-3环PAHs的含量明显高于其他区域，这可能是由于石化区的生产活动以及船舶运输过程中的石油泄漏等原因导致的。4-6环PAHs则主要来源于高温燃烧过程，如机动车尾气排放、工业锅炉燃烧、垃圾焚烧等。在城市附近和交通繁忙的区域，如广州、深圳等城市周边的采样点，4-6环PAHs的相对比例较高，这与城市中大量的机动车行驶以及工业生产活动产生的高温燃烧排放密切相关。环境因素也对PAHs的组成特征产生影响。水体的酸碱度（pH值）会影响PAHs在沉积物中的吸附和解吸过程，从而影响其组成。在酸性条件下，PAHs可能更容易从沉积物中解吸出来，进入水体，导致沉积物中PAHs的组成发生变化。盐度对PAHs的溶解度和分配系数有影响，进而影响其在沉积物中的含量和组成。珠江口是一个河口区域，盐度变化较大，从河口向海洋，盐度逐渐升高。研究发现，随着盐度的增加，一些低环数PAHs的溶解度降低，更容易吸附在沉积物中，导致沉积物中低环数PAHs的相对比例增加。而高环数PAHs由于其分子结构较大，受盐度的影响相对较小。氧化还原电位（Eh）也是影响PAHs组成的重要环境因素。在氧化环境中，PAHs可能会发生光降解和化学氧化等反应，导致其含量和组成发生变化。在珠江口的一些浅水区，水体与大气接触充分，氧化还原电位较高，这些区域的沉积物中PAHs的组成可能会受到氧化作用的影响。而在深水区或缺氧环境中，PAHs的降解速度较慢，其组成相对较为稳定。此外，沉积物的粒度和有机质含量也与PAHs的组成有关。细颗粒沉积物具有较大的比表面积，能够吸附更多的PAHs，且有机质含量较高的沉积物对PAHs的吸附能力更强，这可能导致细颗粒和高有机质含量的沉积物中PAHs的组成与其他沉积物有所不同。三、影响珠江口表层沉积物多环芳烃分布的因素3.1自然因素3.1.1水动力条件水动力条件是影响珠江口表层沉积物中多环芳烃（PAHs）分布的重要自然因素之一，其中河流径流和潮汐作用尤为关键。珠江作为中国南方的重要水系，其径流携带了大量的陆源污染物进入珠江口。在河流径流较大的时期，如夏季的丰水期，河流的流速加快，对污染物的冲刷和搬运能力增强，能够将上游地区的PAHs等污染物快速输送到河口区域。相关研究表明，河流径流的流量与PAHs的输送量存在显著的正相关关系。当径流流量增加时，PAHs的输送量也随之增加。在珠江口的一些河口区域，如西江入海口，由于西江的径流量较大，携带的PAHs等污染物较多，导致该区域表层沉积物中PAHs的含量相对较高。潮汐作用在珠江口也十分显著，珠江口潮汐属不规则半日混合型，一日内潮位两涨两落，日不等现象较显著。在涨潮时，海水会沿河道自河口向上游上溯，将海洋中的物质带入河口，同时也会使河口区域的水体发生混合和搅动，影响PAHs的分布。退潮时，水体携带的PAHs等污染物又会随着水流向海洋扩散。潮汐的这种周期性运动使得PAHs在河口区域呈现出复杂的分布特征。研究发现，在靠近河口的区域，由于潮汐的作用，水体的交换频繁，PAHs更容易在沉积物中富集。而在远离河口的区域，潮汐的影响相对较小，PAHs的分布则更多地受到其他因素的影响。不同水动力区域的PAHs分布存在明显差异。在河口的潮流界附近，河流径流与潮汐的相互作用最为强烈，水体的流速和流向变化复杂。在这个区域，PAHs的分布呈现出高度的不均匀性，部分区域由于水流的汇聚和沉积物的快速沉积，PAHs含量较高；而部分区域由于水流的分散和冲刷，PAHs含量相对较低。在河口的口门处，由于河流径流的外泄和潮汐的顶托，形成了独特的水动力环境，使得PAHs更容易在口门附近的沉积物中积累，导致该区域PAHs含量较高。在近岸海域，水动力条件相对较为稳定，但仍受到潮汐和沿岸流的影响。沿岸流的存在使得近岸海域的水体发生水平运动，将PAHs等污染物沿着海岸线进行输送和扩散。在一些沿岸流较强的区域，PAHs的分布呈现出明显的带状特征，与沿岸流的流向一致。而在一些海湾或河口附近的近岸区域，由于水体的相对封闭性，潮汐的作用更加明显，PAHs容易在这些区域的沉积物中富集。在大鹏湾等相对封闭的海湾，潮汐的涨落使得PAHs在沉积物中的含量呈现出从湾口到湾内逐渐增加的趋势。在远岸海域，水动力条件主要受海洋环流和海浪的影响。海洋环流将远岸海域的水体与其他海域进行交换，使得PAHs等污染物能够在更大范围内扩散。海浪的作用则使得水体发生垂直混合，将表层水体中的PAHs带到深层水体，进而影响PAHs在沉积物中的分布。由于远岸海域距离污染源较远，水体的稀释作用较强，PAHs在沉积物中的含量相对较低，分布也相对较为均匀。3.1.2沉积物性质沉积物性质对珠江口表层沉积物中多环芳烃（PAHs）的分布具有重要影响，其中沉积物粒度和有机碳含量是两个关键因素。沉积物粒度组成反映了其矿物组成、表面物理化学性质、比表面积和表面自由能的差异，这对PAHs的吸附、解吸和迁移有重要影响。通常情况下，沉积物的粒度越细，其比表面积越大，表面自由能越高，对PAHs的吸附能力越强。在珠江口，细颗粒沉积物如黏土和粉砂，由于其较大的比表面积，能够提供更多的吸附位点，使得PAHs更容易被吸附在其表面，从而导致细颗粒沉积物中PAHs的含量相对较高。相关研究表明，珠江口表层沉积物中PAHs含量与黏土和粉砂的含量呈显著正相关关系，而与砂的含量呈负相关关系。在一些河口区域，由于河流携带的细颗粒沉积物较多，这些区域的沉积物中PAHs含量也相对较高。有机碳含量也是影响PAHs在沉积物中分布的重要因素。有机碳具有丰富的官能团和较大的比表面积，对PAHs具有较强的亲和力，能够通过疏水作用、氢键作用等方式与PAHs结合，从而影响PAHs的迁移和转化。研究发现，珠江口表层沉积物中PAHs含量与有机碳含量之间存在显著的正相关关系。有机碳含量较高的沉积物，如河口附近的淤泥质沉积物，由于其含有大量的有机质，能够吸附更多的PAHs，使得这些区域的沉积物中PAHs含量相对较高。沉积物中的有机碳来源多样，包括陆源输入的腐殖质、水生生物的残体以及微生物的代谢产物等。陆源输入的腐殖质通常含有丰富的芳香结构，能够与PAHs发生π-π相互作用，增强对PAHs的吸附能力。水生生物残体在分解过程中会释放出有机物质，这些有机物质也能够吸附PAHs。微生物在代谢过程中产生的多糖、蛋白质等物质，也可以与PAHs结合，影响其在沉积物中的分布。沉积物的性质还会受到其他因素的影响，如盐度、酸碱度等。盐度的变化会影响沉积物表面的电荷性质和离子强度，从而影响PAHs与沉积物之间的相互作用。在盐度较高的海域，沉积物表面的电荷被屏蔽，离子强度增加，使得PAHs与沉积物之间的吸附作用减弱，导致PAHs更容易从沉积物中解吸出来，进入水体。酸碱度的变化会影响PAHs的化学形态和沉积物表面的官能团性质，进而影响PAHs在沉积物中的吸附和解吸过程。在酸性条件下，PAHs可能会发生质子化，使其水溶性增加，从而更容易从沉积物中解吸出来；而在碱性条件下，沉积物表面的官能团可能会发生去质子化，增强对PAHs的吸附能力。三、影响珠江口表层沉积物多环芳烃分布的因素3.2人为因素3.2.1工业活动珠江口周边地区工业发达，工业类型丰富多样，涵盖了石油化工、钢铁冶炼、机械制造、电子电器、纺织印染等多个领域。在石油化工方面，惠州大亚湾石化区是重要的石油化工产业基地，拥有众多大型石化企业，如中海油惠州石化有限公司等，主要进行原油加工、乙烯生产等，在生产过程中会产生大量含有多环芳烃的废气、废水和废渣。钢铁冶炼行业中，广州、佛山等地分布着一些钢铁企业，如广州钢铁企业集团有限公司等，在铁矿石的烧结、炼铁、炼钢等环节，由于高温燃烧和化学反应，会产生多环芳烃并排放到环境中。机械制造企业在金属零部件的加工、表面处理等过程中，使用的润滑油、清洗剂等可能含有多环芳烃，其挥发和泄漏会导致多环芳烃进入周边环境。电子电器产业中，电路板的制造、塑料外壳的注塑等工艺也可能产生多环芳烃。这些工业活动产生的废水和废气是珠江口表层沉积物中多环芳烃的重要来源。工业废水通常含有大量的有机污染物，其中就包括多环芳烃。一些企业在生产过程中未能对废水进行有效处理，直接将含有多环芳烃的废水排放到珠江口及其周边河流中。研究表明，珠江口部分区域表层沉积物中多环芳烃的含量与周边工业废水排放口的距离呈显著负相关关系，距离排放口越近，多环芳烃含量越高。如在靠近南沙化工厂集中区域的河口附近，表层沉积物中多环芳烃含量明显高于其他区域，这与化工厂排放的废水携带大量多环芳烃进入河口沉积物密切相关。工业废气中的多环芳烃主要来源于化石燃料的燃烧和工业生产过程中的挥发。在钢铁冶炼、火力发电等行业，大量燃烧煤炭、石油等化石燃料，会产生含有多环芳烃的废气。这些废气通过大气传输，最终以干湿沉降的方式进入珠江口的水体和沉积物中。有研究通过对大气沉降物中多环芳烃的监测发现，其组成和含量与周边工业活动类型和强度密切相关。在工业集中的区域，大气沉降物中多环芳烃的含量较高，且高环数多环芳烃的比例也相对较大，这是因为高温燃烧过程更容易产生高环数多环芳烃。随着环保意识的提高和环保政策的加强，部分企业开始采取一系列措施来减少多环芳烃的排放，如安装先进的废气处理设备、优化生产工艺等。但仍有一些企业为了降低成本，未能严格执行环保标准，导致多环芳烃的排放问题依然存在。3.2.2交通运输珠江口地区拥有众多港口，如广州港、深圳港等，是中国重要的航运枢纽。港口航运活动频繁，船舶在停靠、装卸货物以及航行过程中会产生多环芳烃。船舶发动机燃烧燃料时，由于不完全燃烧会产生含有多环芳烃的废气排放到大气中。相关研究表明，船舶排放的废气中多环芳烃的含量与船舶的类型、功率、使用的燃料种类以及航行状态等因素有关。大型集装箱船和散货船由于功率较大，燃料消耗多，排放的多环芳烃量也相对较多。船舶在装卸货物过程中，可能会发生油品泄漏等事故，导致石油类物质进入水体，其中的多环芳烃会逐渐吸附到沉积物中。在港口附近的海域，由于船舶活动密集，表层沉积物中多环芳烃的含量明显高于远离港口的区域。有研究对广州港南沙港区附近的沉积物进行分析，发现该区域沉积物中多环芳烃含量显著高于珠江口其他区域，且多环芳烃的组成特征与船舶排放的多环芳烃组成具有相似性，进一步证实了港口航运活动对珠江口表层沉积物中多环芳烃分布的影响。陆地交通尾气也是珠江口表层沉积物中多环芳烃的重要来源之一。珠江口周边城市经济发达，人口密集，机动车保有量持续增长。汽车、卡车等机动车辆在行驶过程中，发动机燃烧汽油或柴油，会产生含有多环芳烃的尾气排放到大气中。尾气中的多环芳烃通过大气传输，最终沉降到珠江口的水体和沉积物中。研究发现，在靠近城市主干道和交通枢纽的区域，大气中多环芳烃的含量较高，相应地，附近的河流和海域表层沉积物中多环芳烃的含量也会增加。有学者对深圳市区周边河流沉积物中多环芳烃的含量进行监测，发现靠近交通繁忙路段的沉积物中多环芳烃含量明显高于远离交通要道的区域，且多环芳烃的含量与车流量呈正相关关系。不同类型的机动车排放的多环芳烃量也有所差异，柴油车由于燃烧特性，排放的多环芳烃量通常比汽油车多，尤其是高环数多环芳烃的排放更为显著。随着环保标准的不断提高，机动车尾气净化技术也在不断发展，如安装三元催化器等设备，可以有效减少多环芳烃等污染物的排放，但由于机动车保有量的持续增加，陆地交通尾气对珠江口表层沉积物中多环芳烃分布的影响仍然不容忽视。四、珠江口表层沉积物多环芳烃来源解析4.1比值法比值法是一种常用的多环芳烃（PAHs）来源解析方法，其原理基于不同来源的PAHs在环境中具有特定的化学组成特征，通过分析PAHs单体之间的比值，可以初步推断其来源。在珠江口表层沉积物PAHs来源解析中，常用的比值包括Ant/(Ant+Phe)、Fla/(Fla+Pyr)等。Ant/(Ant+Phe)比值常被用于区分PAHs的来源是石油源还是燃烧源。当Ant/(Ant+Phe)比值小于0.1时，通常指示PAHs主要来源于石油源，如原油泄漏、石油产品的使用等；当该比值大于0.1时，则表明PAHs主要来源于燃烧源，包括煤、木材、石油等的不完全燃烧。对珠江口表层沉积物样品的分析结果显示，大部分采样点的Ant/(Ant+Phe)比值在0.12-0.25之间，平均值为0.18，这表明珠江口表层沉积物中的PAHs主要来源于燃烧源。这与珠江口周边地区密集的工业活动、交通运输以及居民生活中的燃料燃烧等情况相吻合，如前文所述的石油化工、钢铁冶炼等工业生产过程中化石燃料的燃烧，以及大量机动车行驶和船舶运输过程中燃料的不完全燃烧，都会产生大量源于燃烧的PAHs并排放到环境中。Fla/(Fla+Pyr)比值可进一步区分燃烧源的类型。当Fla/(Fla+Pyr)比值小于0.4时，指示PAHs主要来源于低温燃烧，如生物质燃烧；当比值在0.4-0.5之间时，表明主要来源于石油燃烧；当比值大于0.5时，则主要来源于高温燃烧，如煤的燃烧和机动车尾气排放。在珠江口，部分靠近工业区域和交通要道的采样点，Fla/(Fla+Pyr)比值在0.45-0.6之间，平均值为0.52，说明这些区域的PAHs可能主要来源于高温燃烧，与工业生产中高温加工过程以及机动车尾气排放密切相关。而在一些远离工业区域且植被覆盖较多的采样点，Fla/(Fla+Pyr)比值在0.3-0.4之间，表明这些区域存在一定比例的生物质燃烧来源的PAHs，可能是周边居民生活中使用木材等生物质燃料进行取暖、烹饪等活动导致的。其他比值如BaA/(BaA+Chr)也具有指示意义。当BaA/(BaA+Chr)比值小于0.2时，PAHs主要来源于石油源；当比值在0.2-0.35之间时，来源于石油燃烧和高温燃烧的混合源；当比值大于0.35时，则主要来源于高温燃烧。珠江口部分采样点的BaA/(BaA+Chr)比值在0.3-0.4之间，表明这些区域的PAHs来源于石油燃烧和高温燃烧的混合源，这也进一步印证了珠江口PAHs来源的复杂性，既有石油相关活动（如船舶运输、石油化工）产生的PAHs，也有高温燃烧过程（如工业生产、机动车尾气排放）带来的PAHs。综上所述，通过比值法分析可知，珠江口表层沉积物中的PAHs主要来源于燃烧源，且以高温燃烧源为主，同时存在一定比例的石油源以及生物质燃烧源，这与珠江口周边地区复杂多样的人类活动密切相关。4.2主成分分析为进一步深入剖析珠江口表层沉积物中多环芳烃（PAHs）的来源，本研究运用主成分分析（PCA）方法对16种PAHs的含量数据进行处理。主成分分析是一种常用的多元统计分析方法，它能够将多个相关变量转换为少数几个互不相关的综合变量，即主成分，这些主成分能够最大限度地保留原始数据的信息，从而提取影响PAHs分布的主要因子，确定其主要来源类型。利用统计分析软件，对珠江口各采样点表层沉积物中16种PAHs的含量进行主成分分析。在进行主成分分析之前，对数据进行标准化处理，以消除不同变量之间量纲和数量级的影响。经过分析，提取出3个主成分，其累计方差贡献率达到[X]%，表明这3个主成分能够解释原始数据中大部分的信息。主成分1的方差贡献率为[X1]%，在该主成分中，荧蒽（Fla）、芘（Pyr）、苯并[a]蒽（BaA）、䓛（Chr）、苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IcdP）、二苯并[a,h]蒽（DahA）和苯并[ghi]苝（BghiP）等多环芳烃具有较高的载荷。这些PAHs主要来源于高温燃烧过程，如工业生产中高温加工过程以及机动车尾气排放。在珠江口周边地区，工业活动频繁，众多工厂在生产过程中使用高温炉窑，燃烧化石燃料，会产生大量源于高温燃烧的PAHs。机动车保有量持续增长，尾气排放也是高温燃烧源PAHs的重要来源。因此，主成分1可被认为主要代表了高温燃烧源，包括工业活动和机动车尾气排放对珠江口表层沉积物中PAHs的贡献。主成分2的方差贡献率为[X2]%，其中萘（Nap）、苊烯（Acy）、苊（Ace）、芴（Flu）、菲（Phe）和蒽（Ant）等低环数PAHs具有较高的载荷。这些低环数PAHs主要来源于石油类污染和低温燃烧过程，如原油泄漏、石油产品的使用以及生物质燃烧等。珠江口作为重要的航运枢纽，船舶运输活动频繁，船舶在停靠、装卸货物以及航行过程中，可能会发生油品泄漏等事故，导致石油类物质进入水体，其中的低环数PAHs会逐渐吸附到沉积物中。周边居民生活中使用木材等生物质燃料进行取暖、烹饪等活动，也会产生一定量的低环数PAHs，通过大气传输和沉降进入珠江口沉积物。所以，主成分2主要代表了石油源和生物质燃烧源。主成分3的方差贡献率为[X3]%，该主成分中各PAHs的载荷相对较为分散，但芘（Pyr）和荧蒽（Fla）的载荷相对较高。结合比值法分析结果以及研究区域的实际情况，主成分3可能代表了混合来源，既包含部分石油源和燃烧源的贡献，也可能受到其他一些不确定因素的影响，如区域内特殊的工业生产活动、大气传输过程中的二次污染等。综上所述，通过主成分分析可知，珠江口表层沉积物中PAHs的来源主要包括高温燃烧源（工业活动和机动车尾气排放）、石油源和生物质燃烧源，以及混合来源。这与比值法的分析结果相互印证，进一步明确了珠江口PAHs来源的复杂性，为制定针对性的污染防治措施提供了更全面、准确的科学依据。4.3正定矩阵因子分解模型（PMF）正定矩阵因子分解模型（PositiveMatrixFactorization，PMF）是一种广泛应用于环境科学领域的受体模型，在多环芳烃（PAHs）来源解析中发挥着重要作用。其基本原理是基于统计学的最小二乘法，假设样品中PAHs的浓度矩阵可以分解为两个矩阵的乘积，即源贡献矩阵和源成分谱矩阵。在本研究中，运用PMF模型对珠江口表层沉积物中PAHs的来源进行解析，具体公式如下：x_{ij}=\sum_{k=1}^{p}g_{ik}f_{kj}+e_{ij}其中，x_{ij}表示第i个样品中第j种PAHs的浓度；g_{ik}表示第i个样品中第k个源的贡献；f_{kj}表示第k个源中第j种PAHs的成分谱；e_{ij}表示残差，即模型计算值与实际测量值之间的差异；p表示解析出的源的数目。为了确保模型计算结果的准确性和可靠性，在进行PMF分析之前，对原始数据进行了严格的质量控制和预处理。首先，检查数据的完整性，确保所有采样点的PAHs浓度数据均无缺失值。对于个别因实验误差或其他原因导致的数据异常点，进行了仔细的排查和修正。其次，对数据进行了标准化处理，消除不同PAHs化合物之间浓度差异过大对分析结果的影响。通过将每个PAHs化合物的浓度除以其在所有样品中的最大值，将数据统一到[0,1]的范围内，使得不同PAHs化合物在模型分析中具有相同的权重。在确定源的数目时，采用了逐步增加源数并比较模型拟合优度的方法。从源数p=2开始，逐步增加到p=6，分别运行PMF模型，并计算每个源数下的目标函数值Q。目标函数Q反映了模型计算值与实际测量值之间的差异程度，Q值越小，说明模型拟合效果越好。同时，结合PAHs的化学组成特征、来源指示比值以及研究区域的实际情况，综合判断源数的合理性。当源数为4时，Q值相对较小且稳定，模型拟合效果较好，同时源的成分谱和贡献情况与其他源解析方法的结果具有较好的一致性，因此确定源数为4。经过PMF模型分析，识别出珠江口表层沉积物中PAHs的4个主要来源，各源贡献率及主要成分如下：源1：交通源：贡献率为[X1]%，在该源中，芘（Pyr）、荧蒽（Fla）、苯并[a]蒽（BaA）、䓛（Chr）等具有较高的载荷。这些PAHs是机动车尾气排放的特征组分，在交通繁忙的区域，如城市道路附近和港口周边，机动车和船舶的运行会产生大量此类PAHs。相关研究表明，机动车尾气中PAHs的含量与车辆类型、行驶里程、燃油质量等因素密切相关。在珠江口，大量的货运船舶和客运船只频繁穿梭，其发动机燃烧燃料时会产生不完全燃烧产物，其中就包含芘、荧蒽等PAHs，这些PAHs随着尾气排放到大气中，最终通过干湿沉降等方式进入珠江口表层沉积物。源2：工业源：贡献率为[X2]%，主要成分包括萘（Nap）、苊烯（Acy）、芴（Flu）、菲（Phe）等低环数PAHs，以及苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）等高环数PAHs。工业活动是PAHs的重要来源之一，珠江口周边分布着众多工业企业，如石油化工、钢铁冶炼、机械制造等行业。在石油化工生产过程中，原油的炼制、加工以及化工产品的合成等环节，都会产生含有PAHs的废气、废水和废渣。萘、苊烯等低环数PAHs可能来源于石油类物质的泄漏和挥发，而苯并[b]荧蒽、苯并[a]芘等高环数PAHs则主要是高温燃烧过程的产物。一些工业企业在生产过程中使用的高温炉窑，燃烧煤炭、石油等化石燃料，会产生大量高环数PAHs。源3：燃煤源：贡献率为[X3]%，该源中苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IcdP）、二苯并[a,h]蒽（DahA）和苯并[ghi]苝（BghiP）等具有较高的载荷。燃煤是PAHs的重要人为来源之一，珠江口周边地区存在一定数量的火力发电厂和工业锅炉，这些设施在燃烧煤炭时，由于煤炭的不完全燃烧，会产生大量PAHs排放到环境中。苯并[a]芘、茚并[1,2,3-cd]芘等是燃煤排放的典型PAHs，其含量在燃煤源中相对较高。此外，居民生活中的燃煤取暖、烹饪等活动，也会产生一定量的PAHs，对珠江口表层沉积物中PAHs的贡献不可忽视。源4：生物质燃烧源：贡献率为[X4]%，主要成分包括蒽（Ant）、菲（Phe）等。生物质燃烧是PAHs的自然来源之一，在珠江口周边地区，农业生产活动中秸秆的焚烧，以及森林火灾等自然事件，都会导致生物质燃烧产生PAHs。蒽、菲等PAHs在生物质燃烧过程中大量产生，其在该源中的载荷较高。秸秆焚烧通常在农田附近进行，产生的PAHs会随着大气传输扩散到珠江口区域，进而进入表层沉积物中。通过PMF模型的定量解析，明确了珠江口表层沉积物中PAHs的主要来源及各源的贡献率，为制定针对性的污染防治措施提供了更为精确的科学依据。针对交通源，可加强交通管理，推广清洁能源车辆，优化交通流量，减少机动车尾气排放；对于工业源，应加强对工业企业的监管，推动企业采用清洁生产技术，提高废气、废水处理效率，减少PAHs的排放；对于燃煤源，可加大对火力发电厂和工业锅炉的节能减排改造力度，推广清洁煤技术，降低煤炭燃烧过程中的PAHs排放；对于生物质燃烧源，应加强对农业生产活动的指导，鼓励秸秆综合利用，减少秸秆焚烧，同时加强森林火灾的预防和控制。4.4结果讨论比值法、主成分分析和正定矩阵因子分解模型（PMF）从不同角度对珠江口表层沉积物中多环芳烃（PAHs）的来源进行了解析，三种方法的结果既存在一致性，也有一定差异。在一致性方面，三种方法都表明珠江口表层沉积物中PAHs的来源具有复杂性，并非单一来源。比值法通过特定的PAHs单体比值判断出PAHs主要来源于燃烧源，且以高温燃烧源为主，同时存在石油源和生物质燃烧源。主成分分析提取出的主成分也分别代表了高温燃烧源（工业活动和机动车尾气排放）、石油源和生物质燃烧源以及混合来源。PMF模型解析出的交通源、工业源、燃煤源和生物质燃烧源，同样体现了PAHs来源的多样性，且与比值法和主成分分析中提到的燃烧源、石油源等相对应。这说明不同方法在识别PAHs主要来源类型上具有一定的共识，都反映出珠江口周边复杂的人类活动和自然过程对PAHs来源的影响。在差异方面，比值法相对较为简单直观，主要基于特定PAHs单体的比值来推断来源，但其只能提供定性或半定量的信息，对于各来源的相对贡献难以准确量化。主成分分析虽然能够提取影响PAHs分布的主要因子，确定来源类型，但它是一种基于数据降维的统计方法，对数据的依赖性较强，且在解释因子含义时存在一定的主观性。例如，主成分3被认为代表混合来源，但具体包含哪些来源以及各来源的贡献程度并不明确。PMF模型作为一种受体模型，能够定量解析各来源的贡献率，在确定来源贡献方面具有明显优势。通过PMF模型分析，明确了交通源贡献率为[X1]%，工业源贡献率为[X2]%，燃煤源贡献率为[X3]%，生物质燃烧源贡献率为[X4]%。然而，PMF模型的计算结果受到数据质量、源数确定等因素的影响较大。在本研究中，确定源数为4是通过逐步增加源数并比较模型拟合优度等方法综合判断的，但源数的确定仍然存在一定的不确定性，可能会对结果产生影响。综合三种方法的结果，可以更全面地认识珠江口PAHs的来源。珠江口表层沉积物中PAHs的主要来源包括燃烧源和石油源。燃烧源中，高温燃烧源（工业活动和机动车尾气排放）贡献较大，这与珠江口周边地区工业发达、机动车保有量高的实际情况相符。石油源主要来自船舶运输、石油化工等活动中的油品泄漏和挥发。生物质燃烧源虽然贡献率相对较小，但在部分区域也不容忽视，可能与周边居民生活中的生物质燃料使用以及农业生产活动中的秸秆焚烧等有关。不同来源对PAHs的贡献在空间分布上也存在差异。在河口和近岸区域，工业源和交通源的贡献更为突出。河口地区是陆源污染物进入海洋的重要通道，周边分布着众多工业企业，工业生产活动产生大量PAHs。近岸区域靠近城市和交通要道，机动车尾气排放和港口航运活动频繁，使得交通源的贡献显著。而在远离城市和工业区域的远岸地区，生物质燃烧源和石油源的相对贡献可能会有所增加，这是因为远岸地区人类活动相对较少，工业和交通排放的PAHs减少，而自然过程如生物质燃烧以及海洋中的石油类物质泄漏等对PAHs来源的影响相对增大。综上所述，多种方法的综合运用能够更准确、全面地解析珠江口表层沉积物中PAHs的来源及贡献，为制定有效的污染防治措施提供有力的科学依据。在未来的研究中，可以进一步优化源解析方法，结合更多的环境数据和实地调查，深入探究PAHs来源的时空变化规律，为珠江口的环境保护和生态修复提供更具针对性的建议。五、珠江口表层沉积物多环芳烃生态风险评价5.1评价方法选择多环芳烃（PAHs）因其对生态系统和人类健康的潜在危害，对其进行生态风险评价至关重要。目前，常用的生态风险评价方法众多，每种方法都有其独特的原理和适用范围。沉积物质量基准法中的ERL/ERM法（EffectsRange-Low/EffectsRange-Median）是较为常用的一种方法。该方法由Long等学者提出，通过将沉积物中PAHs的浓度与效应低值（ERL）和效应中值（ERM）进行对比，来评估生态风险。当PAHs浓度低于ERL时，表明沉积物中PAHs对生物产生负面效应的可能性较低；当浓度介于ERL和ERM之间时，可能会偶尔对生物产生负面效应；而当浓度高于ERM时，则频繁对生物产生负面效应。其优点在于简单直观，能够快速判断沉积物中PAHs的风险水平。然而，该方法仅考虑了PAHs的浓度，未考虑不同PAHs之间的毒性差异以及生物可利用性等因素，可能会导致风险评估结果的偏差。商值法，即风险商法（RiskQuotient，RQ），是将环境中污染物的实测浓度（MeasuredConcentration，MC）与预测无效应浓度（PredictedNo-EffectConcentration，PNEC）进行比较，计算风险商值（RQ=MC/PNEC）。当RQ小于1时，认为风险较低；当RQ大于1时，则表明存在一定的生态风险。此方法考虑了污染物的浓度和毒性，相对ERL/ERM法更为全面。但确定准确的PNEC值较为困难，不同的推导方法可能会导致PNEC值存在较大差异，从而影响风险评估的准确性。毒性当量法以苯并[a]芘（BaP）为参照，将其他PAHs根据其毒性等效因子（ToxicEquivalenceFactor，TEF）转化为等效毒性，即毒性当量（ToxicEquivalentQuantity，TEQ），通过计算总毒性当量（∑TEQ）来评估生态风险。该方法充分考虑了不同PAHs的毒性差异，能够更准确地反映PAHs的综合毒性。然而，TEF值的确定存在一定的主观性，不同研究中TEF值可能有所不同，且该方法未考虑环境因素对PAHs毒性的影响。在本研究中，考虑到珠江口表层沉积物中PAHs的污染特点以及数据的可获取性，选择沉积物质量基准法（ERL/ERM法）和毒性当量法相结合的方式对珠江口表层沉积物中PAHs进行生态风险评价。沉积物质量基准法能够初步判断PAHs对生物产生负面效应的可能性，快速了解珠江口不同区域PAHs的风险水平。毒性当量法可以弥补沉积物质量基准法未考虑毒性差异的不足，更准确地评估PAHs的综合毒性，全面反映珠江口表层沉积物中PAHs的生态风险状况。5.2基于沉积物质量基准的风险评价采用沉积物质量基准法中的ERL/ERM（EffectsRange-Low/EffectsRange-Median）法对珠江口表层沉积物中多环芳烃（PAHs）进行生态风险评价。该方法将沉积物中PAHs的浓度与效应低值（ERL）和效应中值（ERM）进行对比，以此评估生态风险。当PAHs浓度低于ERL时，表明沉积物中PAHs对生物产生负面效应的可能性较低；当浓度介于ERL和ERM之间时，可能会偶尔对生物产生负面效应；而当浓度高于ERM时，则频繁对生物产生负面效应。将珠江口各采样点表层沉积物中16种PAHs单体的浓度与ERL和ERM值进行对比，结果如表[X]所示。从表中可以看出，在16种PAHs单体中，菲（Phe）、荧蒽（Fla）、芘（Pyr）、苯并[a]蒽（BaA）、䓛（Chr）、苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IcdP）、二苯并[a,h]蒽（DahA）和苯并[ghi]苝（BghiP）等在部分采样点的浓度超过了ERL值。其中，二苯并[a,h]蒽（DahA）在多个采样点的浓度接近或超过了ERM值，表明该区域二苯并[a,h]蒽对生物产生负面效应的可能性较大，可能会频繁对生物产生负面影响。进一步分析不同区域的风险情况，河口区域如虎门、蕉门等口门附近，由于受到陆源输入和复杂水动力条件的影响，PAHs含量相对较高，超过ERL值的PAHs单体种类和数量较多。在虎门采样点，有[X]种PAHs单体浓度超过ERL值，其中二苯并[a,h]蒽（DahA）的浓度达到[具体数值]ng/g，接近ERM值。这表明河口区域的生态风险相对较高，PAHs对生物产生负面效应的可能性较大。近岸区域，靠近城市和工业集中区的采样点PAHs含量较高，超过ERL值的情况也较为普遍。例如，在靠近深圳的近岸采样点，有[X]种PAHs单体浓度超过ERL值，主要包括荧蒽（Fla）、芘（Pyr）、苯并[a]蒽（BaA）等。这说明近岸区域，尤其是受人类活动影响较大的区域，存在一定的生态风险，PAHs可能会偶尔对生物产生负面效应。远岸区域由于距离污染源较远，PAHs含量相对较低，超过ERL值的PAHs单体种类和数量较少。大部分远岸采样点只有[X]-[X]种PAHs单体浓度超过ERL值，且浓度相对较低。这表明远岸区域的生态风险相对较低，PAHs对生物产生负面效应的可能性较小。根据超过ERL值和ERM值的PAHs单体数量及浓度，将珠江口表层沉积物中PAHs的生态风险划分为三个等级：低风险、中风险和高风险。低风险区域定义为超过ERL值的PAHs单体数量较少，且浓度均远低于ERM值；中风险区域为有一定数量的PAHs单体浓度超过ERL值，部分接近ERM值；高风险区域为有多种PAHs单体浓度超过ERL值，且至少有一种接近或超过ERM值。通过对各采样点的风险等级划分，绘制出珠江口表层沉积物中PAHs生态风险等级分布图，如图[X]所示。从图中可以看出，珠江口的河口区域和部分近岸区域处于中风险和高风险等级，这些区域需要重点关注和治理。而远岸区域大部分处于低风险等级，但仍有个别采样点存在一定的风险。基于沉积物质量基准的风险评价结果表明，珠江口表层沉积物中PAHs存在一定的生态风险，河口和近岸区域的风险相对较高，主要风险物质为二苯并[a,h]蒽（DahA）等。这与珠江口周边地区的工业活动、交通运输等人类活动密切相关，这些活动导致大量PAHs进入珠江口，对生态系统造成潜在威胁。在未来的环境保护和治理工作中，应针对高风险区域，加强对PAHs污染源的管控，减少PAHs的排放，以降低珠江口的生态风险，保护生态系统的健康和稳定。5.3毒性当量评价毒性当量（ToxicEquivalentQuantity，TEQ）评价方法以苯并[a]芘（BaP）作为参照标准，基于不同多环芳烃（PAHs）单体与BaP的毒性差异，通过毒性等效因子（ToxicEquivalenceFactor，TEF）将各PAHs单体的浓度转化为具有统一毒性尺度的等效浓度，从而更准确地评估PAHs的综合毒性效应。其计算公式为：TEQ_{i}=C_{i}\timesTEF_{i}\sumTEQ=\sum_{i=1}^{n}TEQ_{i}其中，TEQ_{i}为第i种PAHs的毒性当量浓度（ng/g），C_{i}为第i种PAHs的实测浓度（ng/g），TEF_{i}为第i种PAHs的毒性等效因子，\sumTEQ为总毒性当量浓度（ng/g），n为PAHs的种类数。不同PAHs的毒性等效因子取值依据国际上相关研究和标准确定，如国际毒性当量因子（I-TEF）体系中，BaP的TEF_{i}被定义为1，其他PAHs的TEF_{i}值根据其相对BaP的毒性强度确定。根据上述公式，对珠江口各采样点表层沉积物中16种PAHs的毒性当量进行计算，结果如表[X]所示。从计算结果可以看出，珠江口表层沉积物中PAHs的总毒性当量（\sumTEQ）范围为[最小值]ng/g至[最大值]ng/g，平均值为[平均值]ng/g。其中，苯并[a]芘（BaP）、二苯并[a,h]蒽（DahA）、苯并[ghi]苝（BghiP）等具有较高毒性等效因子的PAHs单体对总毒性当量的贡献较大。在部分采样点，苯并[a]芘（BaP）的毒性当量浓度达到[具体数值1]ng/g，二苯并[a,h]蒽（DahA）的毒性当量浓度达到[具体数值2]ng/g，二者对总毒性当量的贡献率分别为[X1]%和[X2]%。不同区域的毒性当量存在明显差异。河口区域由于受到陆源输入和复杂水动力条件的影响，PAHs含量相对较高，毒性当量也相对较大。以虎门采样点为例，其总毒性当量达到[具体数值3]ng/g，显著高于珠江口的平均水平。这主要是因为河口地区承接了大量来自上游的工业废水、生活污水等陆源污染物，其中含有较高浓度的高毒性PAHs。近岸区域靠近城市和工业集中区，人类活动频繁，PAHs的排放量大，毒性当量也处于较高水平。在靠近深圳的近岸采样点，总毒性当量在[具体数值4]-[具体数值5]ng/g之间，这与该区域密集的工业活动和大量的机动车尾气排放密切相关。远岸区域距离污染源较远，PAHs含量相对较低，毒性当量也较小，大部分远岸采样点的总毒性当量在[具体数值6]-[具体数值7]ng/g之间。与其他地区的研究结果相比，珠江口表层沉积物中PAHs的毒性当量处于一定水平。在中国国内，与长江三角洲扬子江河口相比，其PAHs的毒性当量范围为[具体范围1]ng/g，珠江口的毒性当量与之相比处于较高区间；与香港米埔和后海湾（[具体范围2]ng/g）相比，珠江口的毒性当量也相对较高。在国际上，与东北大西洋浅滩（[具体范围3]ng/g）和西北大西洋伊比利亚半岛（[具体范围4]ng/g）相比，珠江口表层沉积物中PAHs的毒性当量明显偏高。这表明珠江口表层沉积物中PAHs的毒性风险相对较高，对生态系统和人体健康可能存在较大的潜在威胁。毒性当量评价结果表明，珠江口表层沉积物中PAHs存在一定的毒性风险，河口和近岸区域的风险相对较高。高毒性PAHs单体如苯并[a]芘（BaP）、二苯并[a,h]蒽（DahA）等是主要的风险贡献物质。这些高毒性PAHs具有较强的致癌性和致畸性，通过食物链的传递和富集，可能对珠江口的水生生物和人类健康产生严重危害。在未来的环境保护和治理工作中，应重点关注河口和近岸区域高毒性PAHs的污染问题，加强对PAHs污染源的管控，采取有效的治理措施，降低PAHs的排放，以减少其对生态系统和人体健康的潜在风险。5.4风险分布特征珠江口表层沉积物中多环芳烃（PAHs）的生态风险在不同区域呈现出明显的分布特征。河口区域，如虎门、蕉门等口门附近，是生态风险相对较高的区域。在基于沉积物质量基准的风险评价中，河口区域超过效应低值（ERL）的多环芳烃单体种类和数量较多，部分单体如二苯并[a,h]蒽（DahA）的浓度接近或超过效应中值（ERM），表明该区域PAHs对生物产生负面效应的可能性较大，可能会频繁对生物产生负面影响。在毒性当量评价中，河口区域的总毒性当量也相对较大，如虎门采样点的总毒性当量显著高于珠江口的平均水平，这主要是因为河口地区承接了大量来自上游的工业废水、生活污水等陆源污染物，其中含有较高浓度的高毒性PAHs。河口地区复杂的水动力条件使得污染物容易在此汇聚和沉积，进一步加剧了生态风险。近岸区域，尤其是靠近城市和工业集中区的部分，也存在一定程度的生态风险。从沉积物质量基准来看，近岸区域超过ERL值的PAHs情况较为普遍，如靠近深圳的近岸采样点，有多种PAHs单体浓度超过ERL值，表明PAHs可能会偶尔对生物产生负面效应。毒性当量评价显示，该区域的毒性当量处于较高水平，这与近岸区域密集的工业活动、大量的机动车尾气排放以及港口航运活动密切相关，这些人类活动导致大量PAHs排放到环境中，增加了近岸区域的生态风险。远岸区域的生态风险相对较低。在沉积物质量基准评价中，远岸区域超过ERL值的PAHs单体种类和数量较少，大部分远岸采样点只有少数几种PAHs单体浓度超过ERL值，且浓度相对较低，说明PAHs对生物产生负面效应的可能性较小。毒性当量评价结果也表明，远岸区域的总毒性当量较小，大部分远岸采样点的总毒性当量明显低于河口和近岸区域。这是因为远岸区域距离污染源较远，陆源污染物在传输过程中经过稀释和扩散，浓度逐渐降低，同时远岸海域的水动力条件有利于污染物的扩散和稀释，减少了PAHs在沉积物中的富集。珠江口表层沉积物中PAHs的生态风险分布与PAHs的分布特征密切相关。PAHs含量较高的区域，如河口和近岸的部分地区，