渤海及邻近海域表层沉积物中多环芳烃的溯源与风险评估：基于环境与健康双重视角

渤海及邻近海域表层沉积物中多环芳烃的溯源与风险评估：基于环境与健康双重视角一、引言1.1研究背景与意义多环芳烃（PolycyclicAromaticHydrocarbons，PAHs）是一类由两个或两个以上苯环以稠环形式相连的有机化合物，广泛存在于自然环境和人类社会中。PAHs的来源较为复杂，可分为自然源和人为源。自然源涵盖火山喷发、森林火灾、海洋生物分解以及石油和天然气的自然释放等。而人为源则主要包括工业生产过程中的废气、废水排放，农业生产中农药、化肥的使用，城市生活污水和垃圾的排放，以及交通运输中汽车尾气、船舶废气等。在各类污染源中，人类活动所产生的PAHs在环境中的富集程度相对较高，已然成为PAHs污染的主要诱因。PAHs具有相对稳定性和难降解性，这使得它们易于在环境中富集，并通过食物链进行生物积累。PAHs对人类健康和生态系统安全会产生诸多不良影响。在生态系统方面，PAHs作为典型的海洋有机污染物，会对海洋生态系统造成广泛影响。比如，PAHs能够干扰海洋生物的内分泌系统，影响其生长、发育和繁殖，进而威胁海洋生态系统的稳定性，高浓度的PAHs甚至可能导致海洋生态系统的崩溃。对人类健康而言，PAHs具有致癌、致畸、致突变的“三致”效应。由于PAHs易富集在水产食品中，人类通过食物链的传递间接摄入PAHs，进而对人体健康产生危害。例如，长期食用富含PAHs的鱼类，PAHs会在人体内蓄积，可能引发肝癌、结肠癌、胃癌等多种疾病，对于孕妇、儿童等特殊人群，其危害更为显著。渤海及邻近海域作为中国北方重要的经济区和海洋生态系统，具有极为重要的战略地位和生态价值。该海域周边分布着众多的工业城市和港口，工业活动频繁，交通航运发达，这些人类活动使得渤海及邻近海域面临着严峻的PAHs污染问题。相关研究表明，环渤海经济圈的PAHs排放占中国总排放量的20%。同时，渤海是一个浅的、半封闭的边缘海，水交换能力较低，这限制了污染物通过与开放海洋的混合而稀释，导致PAHs在该海域的沉积物中不断积累。长期以来的粗放式发展模式下的污染排放，使得PAHs在渤海周围的土壤和沉积物中大量累积，通过地表径流等方式，这些残留污染物成为渤海重要的PAHs二次来源。特别是在夏季（6-8月），渤海地区的降水占全年降水总量的70-80%，强降雨冲刷会使环渤海周边地表环境中的颗粒态PAHs大量进入渤海，进一步加剧了该海域的PAHs污染。鉴于PAHs对生态系统和人类健康的潜在威胁，以及渤海及邻近海域PAHs污染的严峻现状，开展对该海域表层沉积物中PAHs的研究显得尤为重要。通过研究PAHs的来源，可以明确污染的主要成因，从而为制定针对性的污染控制措施提供科学依据；评估PAHs的生态风险和健康风险，能够全面了解其对海洋生态系统和人类健康的危害程度，为环境保护和人类健康防护提供决策支持。这不仅有助于保护渤海及邻近海域的生态环境，维护海洋生态系统的平衡和稳定，还对保障当地居民的身体健康、促进区域经济的可持续发展具有深远的意义。1.2国内外研究现状在国际上，针对海洋沉积物中PAHs的研究开展得较早且较为广泛。早期研究主要集中在PAHs的检测分析方法上，随着仪器分析技术的不断进步，气相色谱-质谱联用（GC-MS）、高效液相色谱（HPLC）等技术逐渐成为分析PAHs的主流方法，这些技术的应用使得对PAHs的定性和定量分析更加准确和灵敏。在PAHs来源解析方面，国外学者采用了多种方法。例如，通过研究PAHs的同分异构体比值，如菲/蒽（Phe/Ant）、荧蒽/芘（Flu/Pyr）等，来判断PAHs的来源是石油源还是燃烧源。Yunker等对加拿大弗雷泽河流域的研究发现，PAH比值在判断PAH来源和组成方面具有重要作用，当Phe/Ant比值大于15时，表明PAHs主要来源于石油源，而当该比值小于10时，则主要来源于燃烧源。此外，主成分分析（PCA）、正定矩阵因子分解（PMF）等多元统计分析方法也被广泛应用于PAHs的源解析，这些方法能够综合考虑多种PAHs的浓度数据，更全面地解析PAHs的来源。关于PAHs的生态风险评估，国外学者建立了一系列评估方法和指标。效应区间低/中值法（ERL/ERM）是常用的评估方法之一，该方法通过将沉积物中PAHs的浓度与ERL和ERM值进行比较，来判断PAHs对生物产生负面效应的可能性。Long等通过对大量海洋和河口沉积物的研究，确定了不同PAHs的ERL和ERM值，为生态风险评估提供了重要参考。此外，物种敏感度分布（SSD）法也被用于评估PAHs对不同生物物种的毒性风险，该方法考虑了不同生物对PAHs的敏感性差异，能够更准确地评估PAHs对生态系统的潜在影响。在健康风险评估方面，国外学者主要关注PAHs通过食物链传递对人类健康的影响。研究发现，PAHs在水生生物体内具有生物富集作用，人类食用受污染的水产品可能会摄入大量PAHs，从而增加患癌症等疾病的风险。例如，对波罗的海地区的研究表明，当地居民通过食用受PAHs污染的鱼类，体内PAHs的含量明显升高，健康风险增加。在国内，渤海及邻近海域作为重要的经济和生态区域，其PAHs污染问题也受到了众多学者的关注。李加付等利用GC/MS测定了渤海及其邻近海域表层沉积物中的16种优控多环芳烃，采用多种数据分析技术解析了PAHs的来源，结果表明该海域PAHs主要来自煤炭燃烧源、交通源（石油燃烧）、焦化源和石油源，其贡献分别为54.3%、28.6%、13.4%和3.7%。林田教授团队针对雨季前后渤海水体中多环芳烃的分布、生态风险进行了分析和讨论，研究发现雨季期间，低分子量PAHs从颗粒相转移到溶解相，成为水体中PAHs的次生污染源。国内学者在PAHs的生态风险评估方面也开展了大量研究。如采用ERL/ERM法对渤海部分海域沉积物中PAHs的生态风险进行评估，发现部分站位的PAHs浓度超过了ERL值，存在一定的生态风险。在健康风险评估方面，研究主要集中在通过膳食摄入途径对人体健康的影响，通过建立暴露评估模型，评估了居民通过食用海产品摄入PAHs的健康风险。尽管国内外在渤海及邻近海域PAHs的研究方面取得了一定的成果，但仍存在一些不足之处。一方面，目前的研究多集中在特定区域或特定季节，缺乏对整个渤海及邻近海域长时间序列、全方位的监测和研究，难以全面了解PAHs的时空分布特征和变化规律。另一方面，在PAHs的源解析方面，虽然采用了多种方法，但不同方法之间的结果存在一定差异，且对于一些复杂的混合污染源，解析结果的准确性还有待提高。在生态风险和健康风险评估方面，现有的评估方法大多基于实验室数据和模型模拟，缺乏实际环境中生物和人体暴露的监测数据，导致评估结果与实际情况可能存在偏差。因此，开展更全面、深入的研究，对于准确评估渤海及邻近海域PAHs的污染状况和风险具有重要意义，这也正是本研究的出发点和必要性所在。1.3研究目的与内容本研究旨在全面、深入地了解渤海及邻近海域表层沉积物中多环芳烃（PAHs）的污染状况，准确解析其来源，并科学评估其生态风险和健康风险，为该海域的环境保护和污染治理提供坚实的科学依据和有效的决策支持。具体研究内容如下：PAHs的含量与分布特征研究：系统分析渤海及邻近海域表层沉积物中16种美国环保署（EPA）优先控制的PAHs的含量，利用地理信息系统（GIS）技术绘制其空间分布格局图，深入探究PAHs在不同区域的浓度差异，全面了解PAHs在该海域表层沉积物中的整体污染水平和分布规律。PAHs的来源解析：综合运用分子诊断比值法和主成分分析-多元线性回归（PCA-MLR）模型等多种方法，对PAHs的来源进行定性和定量分析。通过分析PAHs同分异构体的比值，初步判断PAHs的来源类型，再结合PCA-MLR模型，综合考虑多种PAHs的浓度数据，更精确地解析出不同来源对PAHs污染的贡献比例，明确主要污染源，为制定针对性的污染控制措施提供科学依据。PAHs的生态风险评估：采用效应区间低/中值法（ERL/ERM）和平均效应区间中值商法（M-ERM-Q）等方法，对渤海及邻近海域表层沉积物中PAHs的生态风险进行全面评估。通过将沉积物中PAHs的浓度与ERL和ERM值进行对比，判断PAHs对生物产生负面效应的可能性；利用M-ERM-Q法计算综合生态风险指数，更全面地评估PAHs对整个生态系统的潜在影响程度，确定该海域生态风险的高值区域和低值区域，为生态保护提供决策支持。PAHs的健康风险评估：构建基于食物链传递的健康风险评估模型，充分考虑当地居民的膳食结构和海产品摄入量等因素，全面评估PAHs通过食物链对人体健康产生的潜在风险。计算居民通过食用海产品摄入PAHs的日均暴露剂量，结合PAHs的致癌毒性和非致癌毒性数据，评估其对人体健康的危害程度，包括致癌风险和非致癌风险，为保障当地居民的身体健康提供科学依据。1.4研究方法与技术路线本研究采用了多种科学、系统的研究方法，以确保研究结果的准确性和可靠性，具体如下：样品采集：在渤海及邻近海域设置多个采样站位，依据海洋调查规范，利用抓斗式采泥器采集表层沉积物样品。为确保样品的代表性，每个站位采集3-5个子样，并进行混合。同时，详细记录采样时间、地点、经纬度、水深等信息，以便后续分析。样品分析：将采集的沉积物样品冷冻干燥后，研磨过筛。采用索氏提取法，使用正己烷-二氯甲烷混合溶剂对样品中的PAHs进行提取，提取液经硅胶柱和弗罗里硅土柱净化后，使用气相色谱-质谱联用仪（GC-MS）进行分析测定。通过选择离子监测模式（SIM）对16种美国环保署（EPA）优先控制的PAHs进行定性和定量分析，外标法定量。为保证分析结果的准确性，每批样品分析时均设置空白样品、平行样品和加标回收样品，确保回收率在70%-120%之间，相对标准偏差小于10%。来源解析方法：采用分子诊断比值法，通过分析PAHs同分异构体的比值，如菲/蒽（Phe/Ant）、荧蒽/芘（Flu/Pyr）、苯并[a]蒽/（苯并[a]蒽+䓛）（BaA/（BaA+CHR））等，初步判断PAHs的来源类型是石油源、燃烧源还是其他源。同时，运用主成分分析-多元线性回归（PCA-MLR）模型，对多种PAHs的浓度数据进行处理。首先，通过主成分分析提取主要成分，然后将主要成分作为自变量，总PAHs浓度作为因变量，进行多元线性回归分析，从而定量解析不同来源对PAHs污染的贡献比例。生态风险评估方法：运用效应区间低/中值法（ERL/ERM），将沉积物中PAHs的浓度与ERL和ERM值进行对比。当PAHs浓度低于ERL值时，表明PAHs对生物产生负面效应的可能性较低；当PAHs浓度介于ERL和ERM值之间时，可能偶尔会对生物产生负面效应；当PAHs浓度高于ERM值时，则很可能经常对生物产生负面效应。同时，采用平均效应区间中值商法（M-ERM-Q）计算综合生态风险指数，该指数考虑了多种PAHs的联合作用，能够更全面地评估PAHs对整个生态系统的潜在影响程度。M-ERM-Q的计算公式为：M-ERM-Q=\frac{1}{n}\sum_{i=1}^{n}\frac{C_{i}}{ERM_{i}}，其中，n为PAHs的种类数，C_{i}为第i种PAHs的浓度，ERM_{i}为第i种PAHs的效应区间中值。健康风险评估方法：构建基于食物链传递的健康风险评估模型，考虑当地居民的膳食结构和海产品摄入量等因素。通过问卷调查等方式获取当地居民食用海产品的种类、频率和摄入量等信息。计算居民通过食用海产品摄入PAHs的日均暴露剂量（ADD），公式为：ADD=\frac{C\timesIR\timesEF\timesED}{BW\timesAT}，其中，C为海产品中PAHs的浓度，IR为海产品的日均摄入量，EF为暴露频率，ED为暴露持续时间，BW为平均体重，AT为平均暴露时间。结合PAHs的致癌毒性和非致癌毒性数据，评估其对人体健康的危害程度，包括致癌风险（CR）和非致癌风险（HQ）。致癌风险的计算公式为：CR=\sum_{i=1}^{n}ADD_{i}\timesSF_{i}，其中，ADD_{i}为第i种PAHs的日均暴露剂量，SF_{i}为第i种PAHs的致癌斜率因子。非致癌风险的计算公式为：HQ=\sum_{i=1}^{n}\frac{ADD_{i}}{RfD_{i}}，其中，RfD_{i}为第i种PAHs的参考剂量。当CR大于1\times10^{-6}时，表明存在致癌风险；当HQ大于1时，表明存在非致癌风险。本研究的技术路线如图1-1所示，首先进行样品采集与分析，获取渤海及邻近海域表层沉积物中PAHs的浓度数据；然后，运用分子诊断比值法和PCA-MLR模型对PAHs的来源进行解析；接着，采用ERL/ERM法和M-ERM-Q法对PAHs的生态风险进行评估；最后，构建健康风险评估模型，评估PAHs通过食物链对人体健康的潜在风险。通过这一系列步骤，全面、系统地研究渤海及邻近海域表层沉积物中PAHs的来源、生态风险和健康风险。\begin{figure}[htbp]\centering\includegraphics[width=12cm]{技术路线图.jpg}\caption{ç

”究技术路线图}\label{fig:技术路线图}\end{figure}二、研究区域与研究方法2.1研究区域概况渤海地处中国大陆东部北端，是一个近封闭的内海，宛如一颗镶嵌在华夏大地东北边缘的蓝色明珠，其经纬度范围为北纬37°07′-41°0′，东经117°35′-121°10′。它北邻辽宁省，西邻河北省和天津市，南邻山东省，宛如被三省一市温柔环抱。渤海东南面通过渤海海峡与黄海相通，恰似一条纽带连接着两片海域，其海岸线长约3800公里，东西宽约346公里，南北长约550公里，总面积约7.7万平方千米，由辽东湾、渤海湾、莱州湾、中央盆地和渤海海峡5部分组成，每一部分都有着独特的地理风貌和生态特征，共同构成了渤海丰富多彩的海洋生态画卷。渤海海底地形较为平坦，中部微微下凹，仿佛一个巨大的浅碟。其地貌类型多样，主要包括潮滩、平原、沙脊、洼地等。潮滩是海陆交互作用的前沿地带，随着潮汐的涨落，时而被海水淹没，时而露出水面，孕育着独特的潮间带生物群落；平原地势开阔，为海洋生物提供了广阔的栖息和觅食空间；沙脊则像是海底的脊梁，对海水的流动和沉积物的分布产生着重要影响；洼地则成为了一些特殊生物的避风港和栖息地。渤海平均水深18米，最大水深84米，相对较浅的水深使得渤海的水动力条件和生态环境与其他深海区域有着明显的差异。渤海属温带季风气候，四季分明，仿佛一位性格鲜明的艺术家，用不同的笔触描绘着每个季节的独特景致。冬季干寒，凛冽的北风呼啸而过，给海面带来了寒冷的气息，部分海域甚至会出现结冰现象，冰层宛如一面面巨大的镜子，反射着冬日的阳光；夏季湿暖，东南季风带来了丰富的水汽，使得渤海的气温升高，海水温度适宜，为海洋生物的繁殖和生长提供了良好的条件；春秋为过渡季节，气候宜人，仿佛是大自然在两个极端之间的温柔过渡，此时的渤海，海天一色，风景如画。渤海的水文特征独特。沿海多为不正规半日潮，海水的涨落如同大地的呼吸，有着自己独特的节奏，盐度为30‰，这一盐度条件适宜多种海洋生物的生存和繁衍。海况以0-4级为主，相对较为平稳的海况使得渤海成为了重要的航运通道和渔业产区。然而，渤海的水交换能力较低，这限制了污染物通过与开放海洋的混合而稀释，导致污染物容易在该海域的沉积物中不断积累，就像一个缺乏新陈代谢能力的身体，逐渐被污染物所侵蚀。渤海及邻近海域周边分布着众多的工业城市和港口，如大连、天津、唐山、青岛等，这些城市宛如一颗颗璀璨的明珠，镶嵌在渤海的周边。它们是经济发展的重要引擎，同时也是人类活动的密集区域。发达的工业生产、频繁的交通航运以及大量的城市生活污水和垃圾排放，使得渤海及邻近海域面临着严峻的PAHs污染问题。例如，一些工业企业在生产过程中会排放含有PAHs的废气和废水，这些污染物通过大气沉降和地表径流等方式进入渤海；船舶在航行过程中也会排放含有PAHs的废气和废水，直接污染了渤海的水体；城市生活污水和垃圾中含有的PAHs，经过处理后仍有部分残留，最终也会流入渤海。渤海的生态环境丰富多样，宛如一个巨大的生态宝库。海域生物种类繁多，共有国家一、二级保护鸟类60余种，如白鹤、黑鹳、海鸬鹚、黄嘴白鹭、白琵鹭等，它们在渤海的天空中翱翔，为这片海域增添了生机与活力；海域还包含游泳动物114种、浮游植物50种等，这些生物相互依存，构成了复杂的食物链和生态系统。渤海沿岸分布有大连斑海豹国家级自然保护区、双台河口国家级自然保护区等多个生态保护区，这些保护区宛如一个个绿色的卫士，守护着渤海的生态平衡和生物多样性。然而，PAHs的污染对渤海的生态环境产生了严重的威胁，它会干扰海洋生物的内分泌系统，影响其生长、发育和繁殖，导致海洋生物的数量减少和物种多样性降低，甚至可能引发海洋生态系统的崩溃。渤海复杂的地理位置、独特的水文特征以及丰富的生态环境，都对PAHs的分布和迁移产生了深远的影响。相对封闭的海域和较低的水交换能力，使得PAHs容易在沉积物中积累；周边的工业城市和港口等人类活动密集区域，为PAHs的输入提供了丰富的来源；而渤海的气候、水文条件以及生物活动等，又会影响PAHs在水体和沉积物中的迁移、转化和降解过程。因此，深入研究渤海及邻近海域的这些特征，对于准确理解PAHs的污染状况和风险具有重要意义。2.2样品采集与分析本研究的样品采集工作于[具体年份]的[具体月份]进行，此时渤海及邻近海域的水文和气象条件相对稳定，能够较好地反映该海域的常规状况，且避开了如台风、暴雨等极端天气以及主要的渔业捕捞季节，从而保证采集的样品具有代表性。采样范围覆盖了渤海的辽东湾、渤海湾、莱州湾、中央盆地以及邻近海域的主要区域，共设置了[X]个采样站位，站位的分布充分考虑了海域的不同地形地貌、水文条件以及周边人类活动的影响。在辽东湾，站位主要分布在靠近河口和沿岸工业城市的区域，以监测河流输入和陆源污染对该区域的影响；在渤海湾，着重在天津、唐山等工业发达城市的近岸海域设置站位，以评估工业活动产生的污染物排放对海域的影响；莱州湾的站位则兼顾了黄河入海口以及周边渔业养殖区，旨在研究河流输入和渔业活动对该海域的影响；中央盆地的站位相对均匀分布，以获取该区域的背景值；邻近海域的站位则主要设置在与黄海的交界处以及主要航运通道附近，以了解不同海域之间的物质交换以及航运活动对海域的影响。样品采集使用抓斗式采泥器，该采泥器具有操作简便、采集效率高、能够获取较完整的表层沉积物样品等优点。每个站位采集3-5个子样，然后将这些子样充分混合，以减少样品的空间异质性，提高样品的代表性。在采集过程中，使用全球定位系统（GPS）精确定位采样点的经纬度，确保采样位置的准确性，并详细记录采样时间、地点、水深、水温、盐度等现场环境参数。水深通过船载测深仪测量，水温、盐度则使用多参数水质监测仪现场测定，这些参数的记录对于后续分析PAHs在沉积物中的分布与环境因素的关系具有重要意义。样品采集后，立即用聚乙烯密封袋密封，并放置在低温冷藏箱中保存，以防止样品受到外界污染和PAHs的挥发损失。回到实验室后，将样品转移至冰箱中，在4℃条件下冷藏保存，待进一步处理。在样品保存和运输过程中，严格遵守相关的样品保存和运输规范，确保样品的完整性和稳定性。样品处理和分析过程如下：首先，将采集的沉积物样品置于冷冻干燥机中进行冷冻干燥，冷冻干燥能够在低温下除去样品中的水分，避免PAHs在高温干燥过程中发生挥发或分解。干燥后的样品使用玛瑙研钵研磨，使其充分粉碎，并通过100目筛网过筛，以保证样品的均匀性，便于后续的提取和分析。接着，采用索氏提取法对样品中的PAHs进行提取。索氏提取法是一种经典的提取方法，具有提取效率高、提取完全等优点。将研磨过筛后的样品放入滤纸筒中，置于索氏提取器中，加入正己烷-二氯甲烷（体积比为1:1）混合溶剂作为提取剂，在恒温水浴锅中回流提取24小时。在提取过程中，提取剂不断循环，能够充分与样品接触，将样品中的PAHs溶解并提取出来。提取液经过硅胶柱和弗罗里硅土柱净化，以去除其中的杂质和干扰物质，提高分析结果的准确性。硅胶柱和弗罗里硅土柱具有良好的吸附性能，能够有效地吸附提取液中的脂肪、色素、蜡质等杂质。净化后的提取液使用旋转蒸发仪进行浓缩，将提取液的体积减小至1-2mL，以便后续的仪器分析。最后，使用气相色谱-质谱联用仪（GC-MS）对浓缩后的提取液进行分析测定。GC-MS结合了气相色谱的高分离效率和质谱的高灵敏度、高选择性，能够对PAHs进行准确的定性和定量分析。在分析过程中，采用选择离子监测模式（SIM），对16种美国环保署（EPA）优先控制的PAHs进行检测。通过与标准物质的保留时间和质谱图进行对比，确定样品中PAHs的种类，并采用外标法定量，计算出样品中PAHs的含量。为保证分析结果的准确性和可靠性，每批样品分析时均设置空白样品、平行样品和加标回收样品。空白样品用于检测分析过程中是否存在污染，平行样品用于评估分析方法的重复性，加标回收样品用于检验分析方法的准确性和可靠性。实验结果表明，空白样品中未检测到目标PAHs，说明分析过程中无污染；平行样品的相对标准偏差小于10%，表明分析方法具有良好的重复性；加标回收样品的回收率在70%-120%之间，符合分析要求，说明分析方法准确可靠，能够满足本研究的需求。2.3数据分析方法本研究运用了多种数据分析方法，对渤海及邻近海域表层沉积物中PAHs的来源、生态风险和健康风险进行全面且深入的剖析。在PAHs来源解析方面，分子诊断比值法发挥着重要作用。该方法基于不同来源的PAHs在同分异构体比值上存在差异的原理。例如，菲/蒽（Phe/Ant）比值常用于判断PAHs的来源，当Phe/Ant比值大于15时，PAHs主要来源于石油源；当该比值小于10时，则主要来源于燃烧源。荧蒽/芘（Flu/Pyr）比值也具有重要指示意义，当Flu/Pyr比值小于0.5时，表明PAHs主要来源于石油源；当比值在0.5-1.0之间时，来源于石油燃烧源；当比值大于1.0时，则来源于煤和生物质燃烧源。苯并[a]蒽/（苯并[a]蒽+䓛）（BaA/（BaA+CHR））比值同样可用于源解析，当BaA/（BaA+CHR）比值小于0.2时，PAHs主要来源于石油源；当比值在0.2-0.35之间时，来源于石油燃烧和煤燃烧的混合源；当比值大于0.35时，主要来源于煤和生物质燃烧源。通过分析这些特征比值，可以初步判断PAHs的来源类型，为深入解析PAHs的来源提供重要线索。主成分分析-多元线性回归（PCA-MLR）模型则是一种更为复杂和全面的源解析方法。主成分分析（PCA）能够将多个相关变量转化为少数几个不相关的综合变量，即主成分。这些主成分能够最大限度地保留原始数据的信息，同时降低数据的维度，便于后续分析。在本研究中，PCA通过对多种PAHs的浓度数据进行处理，提取出主要成分，这些主要成分代表了不同的PAHs来源类型。例如，某一主成分可能主要反映了石油源PAHs的浓度变化，另一主成分可能主要反映了燃烧源PAHs的浓度变化。多元线性回归（MLR）则是在主成分分析的基础上，将提取出的主成分作为自变量，总PAHs浓度作为因变量，建立回归模型。通过回归分析，可以确定每个主成分对总PAHs浓度的贡献系数，从而定量解析不同来源对PAHs污染的贡献比例。例如，如果某一主成分对应的贡献系数较大，说明该主成分所代表的PAHs来源对总PAHs浓度的贡献较大，即该来源是主要污染源之一。PCA-MLR模型综合了主成分分析和多元线性回归的优势，能够更准确、全面地解析PAHs的来源，为制定针对性的污染控制措施提供科学依据。在生态风险评估方面，效应区间低/中值法（ERL/ERM）是常用的评估方法之一。该方法通过将沉积物中PAHs的浓度与ERL和ERM值进行比较，来判断PAHs对生物产生负面效应的可能性。ERL值是指沉积物中污染物浓度低于该值时，对生物产生负面效应的可能性较低；ERM值是指沉积物中污染物浓度高于该值时，很可能经常对生物产生负面效应；当PAHs浓度介于ERL和ERM值之间时，可能偶尔会对生物产生负面效应。例如，对于萘（Nap），其ERL值为160μg/kg，ERM值为2100μg/kg。如果某一采样站位沉积物中萘的浓度为100μg/kg，低于ERL值，说明该站位萘对生物产生负面效应的可能性较低；如果浓度为1000μg/kg，介于ERL和ERM值之间，则可能偶尔会对生物产生负面效应；如果浓度为3000μg/kg，高于ERM值，则很可能经常对生物产生负面效应。平均效应区间中值商法（M-ERM-Q）则是一种综合考虑多种PAHs联合作用的生态风险评估方法。该方法通过计算M-ERM-Q值，来评估PAHs对整个生态系统的潜在影响程度。M-ERM-Q的计算公式为：M-ERM-Q=\frac{1}{n}\sum_{i=1}^{n}\frac{C_{i}}{ERM_{i}}，其中，n为PAHs的种类数，C_{i}为第i种PAHs的浓度，ERM_{i}为第i种PAHs的效应区间中值。M-ERM-Q值越大，表明PAHs对生态系统的潜在影响程度越高。例如，当M-ERM-Q值小于0.1时，生态风险较低；当值在0.1-0.5之间时，生态风险中等；当值大于0.5时，生态风险较高。通过计算M-ERM-Q值，可以更全面地了解PAHs对生态系统的综合影响，为生态保护提供更有针对性的决策支持。在健康风险评估方面，构建基于食物链传递的健康风险评估模型是关键。该模型充分考虑当地居民的膳食结构和海产品摄入量等因素，通过计算居民通过食用海产品摄入PAHs的日均暴露剂量（ADD），结合PAHs的致癌毒性和非致癌毒性数据，评估其对人体健康的危害程度，包括致癌风险（CR）和非致癌风险（HQ）。日均暴露剂量（ADD）的计算公式为：ADD=\frac{C\timesIR\timesEF\timesED}{BW\timesAT}，其中，C为海产品中PAHs的浓度，IR为海产品的日均摄入量，EF为暴露频率，ED为暴露持续时间，BW为平均体重，AT为平均暴露时间。通过问卷调查等方式获取当地居民食用海产品的种类、频率和摄入量等信息，代入公式即可计算出ADD值。例如，某居民每日食用海产品的量为100g，该海产品中PAHs的浓度为10μg/kg，暴露频率为365天/年，暴露持续时间为70年，平均体重为60kg，平均暴露时间为70年×365天/年，则该居民通过食用该海产品摄入PAHs的日均暴露剂量为：ADD=\frac{10\times100\times365\times70}{60\times70\times365}=16.67μg/（kg·d）。致癌风险（CR）的计算公式为：CR=\sum_{i=1}^{n}ADD_{i}\timesSF_{i}，其中，ADD_{i}为第i种PAHs的日均暴露剂量，SF_{i}为第i种PAHs的致癌斜率因子。致癌斜率因子是指单位剂量的化学物质引起癌症的概率，不同PAHs的致癌斜率因子不同。例如，苯并[a]芘（BaP）的致癌斜率因子为7.3×10⁻²（mg/kg・d）⁻¹。如果某居民通过食用海产品摄入苯并[a]芘的日均暴露剂量为1μg/（kg・d），则其致癌风险为：CR=1\times10^{-3}\times7.3\times10^{-2}=7.3\times10^{-5}。当CR大于1\times10^{-6}时，表明存在致癌风险。非致癌风险（HQ）的计算公式为：HQ=\sum_{i=1}^{n}\frac{ADD_{i}}{RfD_{i}}，其中，RfD_{i}为第i种PAHs的参考剂量。参考剂量是指人类长期暴露于某种化学物质而不产生明显健康危害的每日最大摄入量。例如，萘（Nap）的参考剂量为2×10⁻¹mg/（kg・d）。如果某居民通过食用海产品摄入萘的日均暴露剂量为10μg/（kg・d），则其非致癌风险为：HQ=\frac{10\times10^{-3}}{2\times10^{-1}}=0.05。当HQ大于1时，表明存在非致癌风险。通过这些数据分析方法，本研究能够全面、深入地了解渤海及邻近海域表层沉积物中PAHs的来源、生态风险和健康风险，为该海域的环境保护和污染治理提供科学依据和决策支持。三、渤海及邻近海域表层沉积物中多环芳烃的来源解析3.1多环芳烃的组成与分布特征对渤海及邻近海域[X]个采样站位的表层沉积物样品进行分析后，检测出16种美国环保署（EPA）优先控制的PAHs，分别为萘（Nap）、苊烯（Acy）、苊（Ace）、芴（Flu）、菲（Phe）、蒽（Ant）、荧蒽（Fla）、芘（Pyr）、苯并[a]蒽（BaA）、䓛（Chr）、苯并[b]荧蒽（BbF）、苯并[k]荧蒽（BkF）、苯并[a]芘（BaP）、茚并[1,2,3-cd]芘（IcdP）、二苯并[a,h]蒽（DahA）和苯并[g,h,i]苝（BghiP）。其中，Nap、Phe、Fla、Pyr等低环PAHs在样品中的检出率较高，且浓度相对较大；而Acy、Ace、Ant等在部分样品中的浓度较低，甚至未被检出；高环PAHs如IcdP、DahA、BghiP等的浓度相对较低，但它们具有较强的致癌性和生物累积性，对生态环境和人类健康的潜在危害不容忽视。从PAHs的环数分布来看，3-4环PAHs在总PAHs中所占比例较高，平均值达到[X]%，这表明该海域表层沉积物中的PAHs主要来源于不完全燃烧过程。低环PAHs（2-3环）相对分子质量较小，挥发性较强，更容易通过大气传输和地表径流等方式进入海洋环境；而高环PAHs（5-6环）相对分子质量较大，挥发性较低，更倾向于吸附在沉积物颗粒表面，在海洋环境中具有较高的稳定性。将本研究结果与国内外其他海域的研究数据进行对比，渤海及邻近海域表层沉积物中PAHs的总浓度处于中等污染水平。例如，与珠江口海域相比，渤海及邻近海域PAHs的总浓度相对较低，珠江口海域由于受到工业废水排放、船舶运输等多种因素的影响，PAHs污染较为严重，总浓度较高；而与一些偏远的大洋海域相比，渤海及邻近海域PAHs的总浓度则明显较高，这主要是由于渤海周边人类活动密集，PAHs输入量大。通过地理信息系统（GIS）技术绘制渤海及邻近海域表层沉积物中PAHs的空间分布格局图（图3-1），可以清晰地看出PAHs的浓度在不同区域存在明显差异。总体上呈现出近岸高、远岸低，河口及港湾地区高、开阔海域低的分布特征。在辽东湾，靠近辽河河口和大连等城市的近岸海域PAHs浓度较高，这是因为辽河作为辽东湾的主要入海河流，携带了大量来自流域内工业废水、生活污水以及农业面源污染中的PAHs；大连是重要的工业城市和港口，工业生产和船舶运输等活动排放的PAHs也会对周边海域造成污染。在渤海湾，天津、唐山等城市的近岸海域PAHs浓度显著高于其他区域，天津是中国重要的工业基地和港口城市，工业活动频繁，交通航运发达，PAHs的排放量较大；唐山的钢铁、化工等产业也较为发达，其排放的污染物通过地表径流和大气沉降等方式进入渤海湾，导致该区域PAHs污染较为严重。莱州湾的黄河入海口附近PAHs浓度较高，黄河是中国的第二长河，其流域面积广，流经多个工业城市和农业产区，携带了大量的污染物，其中包括PAHs；同时，莱州湾周边的渔业养殖和盐业生产等活动也可能对海域造成一定的污染。中央盆地和开阔海域的PAHs浓度相对较低，这是因为这些区域距离污染源较远，受到人类活动的影响较小，且水动力条件相对较强，有利于污染物的扩散和稀释。\begin{figure}[htbp]\centering\includegraphics[width=12cm]{PAHs空间分布æ

¼å±€å›¾.jpg}\caption{渤海及邻近海域表层沉积物中PAHs的空间分布æ

¼å±€å›¾}\label{fig:PAHs空间分布æ

¼å±€å›¾}\end{figure}影响PAHs在渤海及邻近海域表层沉积物中分布的因素是多方面的，主要包括污染源分布、水动力条件和沉积物性质等。污染源分布是影响PAHs浓度空间分布的关键因素，如前文所述，近岸和河口地区由于受到工业废水排放、生活污水排放、船舶运输、大气沉降等多种污染源的影响，PAHs输入量大，导致这些区域的PAHs浓度较高。水动力条件对PAHs的分布也有着重要影响，渤海的海流主要由渤海海峡流入的黄海暖流余脉和沿岸流组成，海流的运动可以携带PAHs在海域内扩散。在水动力较强的区域，如开阔海域和海峡附近，PAHs能够更快地扩散和稀释，浓度相对较低；而在水动力较弱的区域，如港湾和河口地区，PAHs容易积聚，浓度相对较高。此外，潮汐的涨落也会影响PAHs在沉积物中的分布，潮汐可以促进水体与沉积物之间的物质交换，使PAHs在沉积物中的含量发生变化。沉积物性质也是影响PAHs分布的重要因素之一，沉积物中的有机碳含量、粒度等会影响PAHs的吸附和解析行为。一般来说，有机碳含量较高的沉积物对PAHs具有较强的吸附能力，能够固定更多的PAHs，从而使沉积物中的PAHs浓度升高；而粒度较细的沉积物比表面积较大，也有利于PAHs的吸附。研究表明，渤海及邻近海域表层沉积物中PAHs的浓度与有机碳含量呈显著正相关，与沉积物粒度呈负相关。3.2多环芳烃的来源识别方法为准确剖析渤海及邻近海域表层沉积物中PAHs的来源，本研究综合运用了多种方法，包括分子诊断比值法和主成分分析-多元线性回归（PCA-MLR）模型等，这些方法相互补充，从不同角度揭示了PAHs的来源信息。分子诊断比值法是基于不同来源的PAHs在同分异构体比值上存在差异这一特性来判断PAHs的来源类型。例如，菲/蒽（Phe/Ant）比值常用于判断PAHs的来源，当Phe/Ant比值大于15时，PAHs主要来源于石油源；当该比值小于10时，则主要来源于燃烧源。这是因为在石油的形成和开采过程中，菲和蒽的生成机制相对稳定，使得石油源PAHs中Phe/Ant比值较高；而在燃烧过程中，由于高温和复杂的化学反应，菲和蒽的生成比例发生变化，导致燃烧源PAHs中Phe/Ant比值较低。荧蒽/芘（Flu/Pyr）比值也具有重要指示意义，当Flu/Pyr比值小于0.5时，表明PAHs主要来源于石油源；当比值在0.5-1.0之间时，来源于石油燃烧源；当比值大于1.0时，则来源于煤和生物质燃烧源。这是因为不同的燃烧过程会产生不同比例的荧蒽和芘，石油燃烧过程中产生的荧蒽相对较少，而煤和生物质燃烧过程中产生的荧蒽相对较多。苯并[a]蒽/（苯并[a]蒽+䓛）（BaA/（BaA+CHR））比值同样可用于源解析，当BaA/（BaA+CHR）比值小于0.2时，PAHs主要来源于石油源；当比值在0.2-0.35之间时，来源于石油燃烧和煤燃烧的混合源；当比值大于0.35时，主要来源于煤和生物质燃烧源。这是因为不同来源的PAHs在苯并[a]蒽和䓛的生成比例上存在差异，石油源PAHs中苯并[a]蒽的含量相对较低，而煤和生物质燃烧源PAHs中苯并[a]蒽的含量相对较高。主成分分析-多元线性回归（PCA-MLR）模型则是一种更为复杂和全面的源解析方法。主成分分析（PCA）能够将多个相关变量转化为少数几个不相关的综合变量，即主成分。这些主成分能够最大限度地保留原始数据的信息，同时降低数据的维度，便于后续分析。在本研究中，PCA通过对多种PAHs的浓度数据进行处理，提取出主要成分，这些主要成分代表了不同的PAHs来源类型。例如，某一主成分可能主要反映了石油源PAHs的浓度变化，另一主成分可能主要反映了燃烧源PAHs的浓度变化。多元线性回归（MLR）则是在主成分分析的基础上，将提取出的主成分作为自变量，总PAHs浓度作为因变量，建立回归模型。通过回归分析，可以确定每个主成分对总PAHs浓度的贡献系数，从而定量解析不同来源对PAHs污染的贡献比例。例如，如果某一主成分对应的贡献系数较大，说明该主成分所代表的PAHs来源对总PAHs浓度的贡献较大，即该来源是主要污染源之一。PCA-MLR模型综合了主成分分析和多元线性回归的优势，能够更准确、全面地解析PAHs的来源，为制定针对性的污染控制措施提供科学依据。3.3自然来源分析自然来源对渤海及邻近海域表层沉积物中PAHs的贡献不可忽视，其中火山喷发和森林火灾是较为重要的自然因素。火山喷发是一种强烈的地质活动，在喷发过程中，大量的岩浆和火山灰被释放到大气中，这些物质中可能含有PAHs。例如，冰岛的埃亚菲亚德拉火山喷发时，火山灰中检测出了多种PAHs。虽然渤海及邻近海域目前没有近期的火山喷发记录，但在地质历史时期，该区域可能受到过火山活动的影响，这些历史时期火山喷发产生的PAHs可能会残留在表层沉积物中。森林火灾也是PAHs的一个重要自然来源。当森林发生火灾时，树木和植被中的有机物质在高温下不完全燃烧，会产生大量的PAHs。这些PAHs会随着大气传输，最终通过干湿沉降等方式进入渤海及邻近海域。以澳大利亚的森林大火为例，2019-2020年澳大利亚发生了大规模的森林火灾，火灾产生的大量烟雾中含有丰富的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.4人为来源分析3.4.1工业污染工业活动是渤海及邻近海域表层沉积物中PAHs的重要人为来源之一，其中石化、钢铁等行业的影响尤为显著。石化行业在原油开采、炼制、运输以及石油产品的生产过程中，会不可避免地产生PAHs。例如，在原油开采过程中，油井的泄漏和溢油事故会导致石油类物质进入海洋环境，其中就包含大量的PAHs。据统计，全球每年因石油开采和运输过程中的泄漏而进入海洋的石油量可达数百万吨，这些石油中的PAHs会在海洋中扩散和迁移，最终吸附在沉积物颗粒表面，进入表层沉积物。在石油炼制过程中，高温裂解、催化重整等工艺会使石油中的有机物质发生复杂的化学反应，产生多种PAHs。某大型石化企业的研究表明，其炼油过程中排放的废气和废水中含有大量的PAHs，其中萘、菲、芘等低环PAHs的浓度较高。这些PAHs通过大气沉降和地表径流等方式进入渤海及邻近海域，对该海域的生态环境造成了严重威胁。钢铁行业也是PAHs的重要排放源。在钢铁生产过程中，烧结、焦化、炼铁、炼钢等工序都会产生PAHs。以烧结工序为例，铁矿石、焦炭等原料在高温烧结过程中，其中的有机物质会发生不完全燃烧，产生大量的PAHs。有研究表明，钢铁企业烧结机排放的废气中，PAHs的浓度可高达数百微克每立方米。这些废气中的PAHs通过大气传输，最终沉降到渤海及邻近海域，成为该海域PAHs的重要来源之一。焦化过程中，煤炭在高温干馏条件下分解，会产生大量的煤焦油和煤气，其中煤焦油中含有丰富的PAHs。据估算，每生产1吨焦炭，大约会产生10-30千克的煤焦油，而煤焦油中PAHs的含量可高达30%-50%。这些PAHs随着焦化废水和废气的排放进入环境，对渤海及邻近海域的生态环境产生了潜在的危害。此外，化工、电力等行业也会排放一定量的PAHs。化工行业在生产有机化学品、塑料、橡胶等产品时，会使用大量的有机原料，这些原料在反应过程中可能会产生PAHs。电力行业中的火力发电，尤其是燃煤发电，煤炭在燃烧过程中会产生大量的PAHs。某燃煤发电厂的研究发现，其烟囱排放的飞灰中含有多种PAHs，其中苯并[a]芘等强致癌性PAHs的含量较高。这些飞灰通过大气沉降进入渤海及邻近海域，增加了该海域PAHs的污染负荷。为了评估工业污染对渤海及邻近海域PAHs污染的贡献，本研究运用主成分分析-多元线性回归（PCA-MLR）模型对数据进行分析。结果显示，工业源对该海域表层沉积物中PAHs的贡献率达到[X]%。在辽东湾靠近石化企业和钢铁厂的海域，PAHs的浓度明显高于其他区域，这进一步证实了工业污染是该海域PAHs的重要来源。在某石化企业附近的采样站位，PAHs的总浓度达到了[具体浓度]，远高于渤海及邻近海域的平均水平。这表明，工业活动排放的PAHs在局部海域造成了严重的污染，对海洋生态环境和人类健康构成了潜在威胁。3.4.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的比例相对较高，而汽油车排放的PAHs中，低环PAHs的比例相对较高。为了评估交通运输对渤海及邻近海域PAHs污染的贡献，本研究运用主成分分析-多元线性回归（PCA-MLR）模型对数据进行分析。结果显示，交通运输源对该海域表层沉积物中PAHs的贡献率达到[X]%。在渤海湾的天津港和唐山港等港口附近，PAHs的浓度明显高于其他区域，这主要是由于港口船舶运输活动频繁，船舶排放的PAHs较多。在天津港附近的采样站位，PAHs的总浓度达到了[具体浓度]，远高于渤海及邻近海域的平均水平。这表明，船舶运输对该区域的PAHs污染产生了重要影响。在靠近城市道路的海域，PAHs的浓度也相对较高，这与汽车尾气排放密切相关。在某城市沿海公路附近的采样站位，PAHs的总浓度比远离公路的站位高出[X]%。这进一步证实了汽车尾气排放是渤海及邻近海域PAHs的重要来源之一。3.4.3其他人为来源农业污染也是渤海及邻近海域表层沉积物中PAHs的一个不可忽视的人为来源。农业生产中广泛使用的农药和化肥，其中部分成分含有PAHs。例如，一些有机氯农药和有机磷农药在生产过程中可能会产生PAHs杂质。这些农药在使用过程中，会通过大气漂移、地表径流等方式进入海洋环境。有研究表明，在使用农药较多的农田附近的河流中，PAHs的浓度明显高于未使用农药的区域。这些河流最终汇入渤海，导致渤海及邻近海域的PAHs污染加重。此外，农业废弃物的焚烧也是PAHs的一个来源。在农村地区，农作物秸秆、杂草等农业废弃物常常被露天焚烧，在焚烧过程中，有机物质不完全燃烧会产生大量的PAHs。这些PAHs会随着大气传输，最终沉降到渤海及邻近海域。据估算，每年因农业废弃物焚烧排放到大气中的PAHs量可达数万吨。城市生活污染同样对渤海及邻近海域的PAHs污染有重要影响。城市居民日常生活中产生的污水和垃圾是PAHs的重要来源。城市污水中含有大量的有机物质，这些有机物质在微生物的作用下分解，可能会产生PAHs。此外，污水中还可能含有来自家庭清洁用品、化妆品等的PAHs。城市垃圾中也含有各种有机物质，在垃圾填埋或焚烧过程中，会产生PAHs。研究表明，垃圾焚烧产生的飞灰中含有大量的PAHs，其中一些高环PAHs具有较强的致癌性。这些飞灰如果未经妥善处理，会通过大气沉降等方式进入渤海及邻近海域。为了评估这些其他人为来源对渤海及邻近海域PAHs污染的贡献，本研究运用主成分分析-多元线性回归（PCA-MLR）模型对数据进行分析。结果显示，农业源和生活源对该海域表层沉积物中PAHs的贡献率分别达到[X]%和[X]%。在莱州湾靠近黄河入海口的区域，由于黄河流域农业活动频繁，该区域的PAHs浓度相对较高。在某靠近黄河入海口的采样站位，PAHs的总浓度比远离入海口的站位高出[X]%。这表明，农业污染对该区域的PAHs污染产生了重要影响。在城市周边的海域，由于城市生活污水和垃圾排放的影响，PAHs的浓度也相对较高。在某城市附近的采样站位，PAHs的总浓度比远离城市的站位高出[X]%。这进一步证实了城市生活污染是渤海及邻近海域PAHs的重要来源之一。四、渤海及邻近海域表层沉积物中多环芳烃的生态风险评估4.1生态风险评估方法效应区间低值/效应区间中值法（ERL/ERM）是一种被广泛应用于评估沉积物中污染物生态风险的方法。该方法由Long等学者于1995年提出，其核心原理是通过将沉积物中污染物的浓度与预先设定的效应区间低值（ERL）和效应区间中值（ERM）进行对比，以此来判断污染物对生物产生负面效应的可能性。当沉积物中PAHs的浓度低于ERL值时，表明PAHs对生物产生负面效应的可能性较低；当浓度介于ERL和ERM值之间时，可能偶尔会对生物产生负面效应；而当浓度高于ERM值时，则很可能经常对生物产生负面效应。这种方法的优势在于其简单直观，能够快速地对生态风险进行初步评估。例如，在对某海域沉积物中PAHs的生态风险评估中，通过将各采样点PAHs的浓度与ERL和ERM值进行比较，能够清晰地判断出哪些区域的生态风险较高，哪些区域相对较低。然而，该方法也存在一定的局限性。它只是基于污染物的浓度进行评估，没有充分考虑到生物的敏感性差异、污染物之间的协同作用以及环境因素的影响。不同生物对PAHs的敏感性可能存在很大差异，有些生物可能对较低浓度的PAHs就非常敏感，而有些生物则可能具有较强的耐受性。此外，环境因素如温度、盐度、pH值等也会影响PAHs的毒性和生物可利用性。在高温环境下，PAHs的挥发性可能增强，从而降低其在沉积物中的浓度，但同时也可能增加其对生物的暴露风险。平均效应区间中值商法（M-ERM-Q）则是一种综合考虑多种PAHs联合作用的生态风险评估方法。该方法的计算公式为M-ERM-Q=\frac{1}{n}\sum_{i=1}^{n}\frac{C_{i}}{ERM_{i}}，其中，n为PAHs的种类数，C_{i}为第i种PAHs的浓度，ERM_{i}为第i种PAHs的效应区间中值。M-ERM-Q值越大，表明PAHs对生态系统的潜在影响程度越高。当M-ERM-Q值小于0.1时，生态风险较低；当值在0.1-0.5之间时，生态风险中等；当值大于0.5时，生态风险较高。与ERL/ERM法相比，M-ERM-Q法的优势在于它考虑了多种PAHs的联合毒性效应。在实际环境中，PAHs往往是以混合物的形式存在，不同PAHs之间可能存在协同作用，使得它们的联合毒性大于单个PAHs的毒性之和。M-ERM-Q法能够更全面地评估PAHs对生态系统的综合影响。然而，该方法也有其不足之处。它需要准确获取每种PAHs的ERM值，而这些值的确定往往需要大量的实验研究和数据积累，不同研究得到的ERM值可能存在差异，这会影响评估结果的准确性。此外，M-ERM-Q法同样没有充分考虑生物敏感性和环境因素的影响。除了上述两种常用方法外，物种敏感度分布（SSD）法也是一种重要的生态风险评估方法。该方法通过构建不同生物物种对PAHs的敏感度分布曲线，来评估PAHs对生态系统中不同生物物种的毒性风险。它考虑了不同生物对PAHs的敏感性差异，能够更准确地评估PAHs对生态系统的潜在影响。例如，通过对多种海洋生物进行毒性测试，获取它们对PAHs的半数致死浓度（LC50）或半数抑制浓度（IC50）等数据，然后利用这些数据构建SSD曲线。根据SSD曲线，可以确定在一定风险水平下，PAHs对多少比例的生物物种会产生毒性效应。SSD法的优势在于它能够为生态保护提供更具针对性的信息，有助于制定更合理的环境保护政策。然而，该方法也面临一些挑战。获取大量不同生物物种的毒性数据需要耗费大量的时间、人力和物力，而且不同生物之间的实验条件难以完全统一，这可能会导致数据的可比性受到影响。此外，SSD法假设生物物种之间是相互独立的，但在实际生态系统中，生物之间存在复杂的相互关系，这种假设可能会使评估结果与实际情况存在一定偏差。在本研究中，综合运用ERL/ERM法和M-ERM-Q法对渤海及邻近海域表层沉积物中PAHs的生态风险进行评估。首先，使用ERL/ERM法对每个采样站位的PAHs浓度进行单独评估，判断每种PAHs对生物产生负面效应的可能性。然后，运用M-ERM-Q法计算每个采样站位的综合生态风险指数，全面评估多种PAHs的联合作用对生态系统的潜在影响。通过这两种方法的结合，可以更全面、准确地了解渤海及邻近海域表层沉积物中PAHs的生态风险状况。4.2对海洋生物的影响PAHs对渤海及邻近海域的海洋生物产生了多方面的负面影响，严重威胁着海洋生态系统的稳定和生物多样性。在生存方面，PAHs的毒性会导致海洋生物死亡率上升。研究表明，当海洋生物暴露于含有PAHs的环境中时，其生理机能会受到损害，从而影响其生存能力。以渤海常见的鱼类为例，在PAHs污染较为严重的海域，某些鱼类的死亡率明显高于无污染海域。有研究对辽东湾某PAHs污染区域的鱼类进行调查，发现该区域鱼类的死亡率比对照区域高出[X]%。这是因为PAHs进入鱼类体内后，会干扰其正常的生理代谢过程，影响其呼吸、循环等系统的功能，导致鱼类无法正常生存。在发育方面，PAHs会干扰海洋生物的正常发育过程，导致发育异常。对海洋鱼类胚胎发育的研究发现，PAHs会引发胚胎死亡、畸形变化、发育延迟等问题。在对渤海某海域的研究中，发现受到PAHs污染的鱼类胚胎，出现了脊柱弯曲、围心腔水肿等畸形现象，胚胎死亡率也显著增加。当鱼类胚胎暴露于含有苯并[a]芘（BaP）的环境中时，随着BaP浓度的升高，胚胎死亡率明显上升，畸形率也显著增加。这是因为PAHs可以被鱼类体内的肝脏氧化酶氧化成更加有毒的代谢产物，如羟基化产物等，这些代谢产物会损害胚胎细胞的结构和功能，干扰胚胎的正常发育。PAHs对海洋生物的繁殖也产生了显著的影响。研究表明，PAHs会降低海洋生物的繁殖能力，影响其生殖系统的发育和功能。在对渤海海域的贝类研究中发现，长期暴露于PAHs污染环境中的贝类，其生殖腺发育受到抑制，产卵量明显减少。对某PAHs污染海域的贝类进行监测，发现其产卵量比无污染海域的贝类减少了[X]%。此外，PAHs还会影响海洋生物的精子质量和卵子受精率，导致生殖失败。以日本虎斑猛水蚤为受试生物的研究表明，菲（Phe）暴露会显著降低其挂卵率，在高浓度Phe暴露组甚至未出现挂卵现象，导致生殖失败。这是因为PAHs会干扰海洋生物体内的内分泌系统，影响生殖激素的合成和分泌，从而影响生殖过程。PAHs对海洋生物的影响还具有食物链传递效应。低营养级的海洋生物如浮游生物、藻类等，会吸收环境中的PAHs。这些生物被高营养级的生物捕食后，PAHs会在食物链中逐渐积累和放大，对高营养级生物产生更大的危害。在渤海的食物链中，浮游生物是初级生产者，它们会吸收海水中的PAHs。小鱼以浮游生物为食，PAHs会在小鱼体内积累。大鱼又以小鱼为食，PAHs在大鱼体内的浓度会进一步升高。最终，处于食物链顶端的海洋哺乳动物和人类，通过食用受污染的海产品，会摄入大量的PAHs，从而对健康产生潜在威胁。有研究表明，渤海海域的一些海洋哺乳动物体内检测出了较高浓度的PAHs，这与它们长期食用受污染的鱼类密切相关。4.3生态风险的空间分布特征利用地理信息系统（GIS）技术，绘制渤海及邻近海域表层沉积物中PAHs生态风险的空间分布图（图4-1），以便直观地展示生态风险的空间分布特征。从图中可以清晰地看出，渤海及邻近海域表层沉积物中PAHs的生态风险呈现出明显的空间差异。在辽东湾，靠近辽河河口和大连等城市的近岸海域生态风险相对较高。以辽河河口附近的站位为例，该站位的M-ERM-Q值达到了0.45，处于中等风险水平。这主要是因为辽河作为辽东湾的主要入海河流，携带了大量来自流域内工业废水、生活污水以及农业面源污染中的PAHs。这些PAHs在河口附近大量沉积，导致该区域的PAHs浓度升高，从而增加了生态风险。此外，大连作为重要的工业城市和港口，工业生产和船舶运输等活动排放的PAHs也会对周边海域造成污染，进一步加剧了该区域的生态风险。渤海湾的天津、唐山等城市的近岸海域是生态风险的高值区域。在天津港附近的站位，M-ERM-Q值高达0.62，属于高风险水平。天津是中国重要的工业基地和港口城市，工业活动频繁，交通航运发达，PAHs的排放量较大。这些PAHs通过大气沉降、地表径流和船舶排放等方式进入渤海湾，使得该区域的PAHs浓度显著高于其他区域，生态风险也随之升高。唐山的钢铁、化工等产业也较为发达，其排放的污染物通过地表径流和大气沉降等方式进入渤海湾，对该区域的生态环境造成了严重威胁，进一步提高了该区域的生态风险。莱州湾的黄河入海口附近生态风险也相对较高。在黄河入海口附近的站位，M-ERM-Q值为0.38，处于中等风险水平。黄河是中国的第二长河，其流域面积广，流经多个工业城市和农业产区，携带了大量的污染物，其中包括PAHs。这些PAHs随着黄河水进入莱州湾，在入海口附近大量沉积，导致该区域的PAHs浓度升高，生态风险增加。此外，莱州湾周边的渔业养殖和盐业生产等活动也可能对海域造成一定的污染，对该区域的生态风险产生一定的影响。中央盆地和开阔海域的生态风险相对较低。在中央盆地的一些站位，M-ERM-Q值仅为0.08，处于低风险水平。这是因为这些区域距离污染源较远，受到人类活动的影响较小，且水动力条件相对较强，有利于污染物的扩散和稀释。海流的运动可以将PAHs携带到其他区域，降低了该区域的PAHs浓度，从而降低了生态风险。此外，开阔海域的水体交换能力较强，能够更快地将污染物稀释，减少了PAHs在沉积物中的积累，进一步降低了生态风险。总体而言，渤海及邻近海域表层沉积物中PAHs的生态风险呈现出近岸高、远岸低，河口及港湾地区高、开阔海域低的分布特征。这种空间分布特征与PAHs的来源和迁移转化规律密切相关。近岸和河口地区由于受到工业废水排放、生活污水排放、船舶运输、大气沉降等多种污染源的影响，PAHs输入量大，导致这些区域的PAHs浓度较高，生态风险也相应较高。而远岸和开阔海域距离污染源较远，受到人类活动的影响较小，且水动力条件相对较强，有利于污染物的扩散和稀释，因此生态风险相对较低。\begin{figure}[htbp]\centering\includegraphics[width=12cm]{生态风险空间分布图.jpg}\caption{渤海及邻近海域表层沉积物中PAHs生态风险的空间分布图}\label{fig:生态风险空间分布图}\end{figure}4.4不确定性分析在本研究中，评估过程中存在多个不确定性因素，这些因素对结果产生了不同程度的影响。从数据误差方面来看，样品采集的代表性是一个关键问题。尽管在采样时设置了多个站位，并对每个站位采集多个子样进行混合，但渤海及邻近海域面积广阔，地形和环境复杂多样，仍然难以完全确保采集的样品能够准确代表整个海域的情况。在某些偏远或难以到达的区域，可能无法进行采样，这就导致这些区域的数据缺失，从而影响评估结果的全面性和准确性。分析检测过程也可能引入误差。虽然采用了先进的气相色谱-质谱联用仪（GC-MS）进行分析测定，并严格控制实验条件，但仪器本身的精度限制以及实验操作过程中的微小差异，都可能导致检测结果的偏差。在样品提取过程中，提取效率可能无法达到100%，导致部分PAHs未被完全提取出来，从而使检测结果偏低；在仪器分析过程中，噪音、基线漂移等因素也可能影响峰的识别和定量，导致分析结果出现误差。在模型假设方面，无论是来源解析模型还是风险评估模型，都存在一定的假设条件。在主成分分析-多元线性回归（PCA-MLR）模型中，假设不同来源的PAHs之间相互独立，不存在交叉影响。然而，在实际环境中，不同来源的PAHs可能会发生相互作用，例如在大气传输过程中，石油源和燃烧源的PAHs可能会混合在一起，然后通过大气沉降进入海洋，这就违背了模型的假设条件，从而可能导致源解析结果的偏差。在生态风险评估中，效应区间低/中值法（ERL/ERM）和平均效应区间中值商法（M-ERM-Q）等方法假设PAHs的毒性效应是线性的，且不考虑生物的适应性和修复能力。但实际上，生物在长期暴露于PAHs的环境中可能会产生适应性，其对PAHs的耐受性可能会增强；同时，生物自身也具有一定的修复能力，能够在一定程度上减轻PAHs对其造成的损害。这些因素都没有在模型中得到充分考虑，可能导致生态风险评估结果的不准确。在健康风险评估中，基于食物链传递的健康风险评估模型假设当地居民的膳食结构和海产品摄入量在一段时间内保持稳定。然而，随着经济的发展和人们生活水平的提高，居民的膳食结构可能会发生变化，海产品的摄入量也可能会有所波动。此外，模型中使用的PAHs的致癌斜率因子和参考剂量等参数，是基于大量的实验研究和统计分析得到的，但这些参数可能存在一定的不确定性，不同的研究结果可能会有所差异，这也会影响健康风险评估结果的准确性。为了减少不确定性因素对结果的影响，本研究采取了一系列措施。在样品采集方面，尽可能增加采样站位的数量和分布范围，提高样品的代表性；在分析检测过程中，严格按照标准操作规程进行实验，增加平行样品的数量，对分析结果进行多次验证，以降低误差。对于模型假设，在分析结果时充分考虑模型的局限性，结合实际情况进行综合判断，并