基于被动采样解析珠三角城市群大气持久性有机污染物的污染与风险

基于被动采样解析珠三角城市群大气持久性有机污染物的污染与风险一、引言1.1研究背景与意义在经济高速发展和城市化进程不断加快的背景下，珠三角城市群作为中国重要的经济区之一，面临着严峻的大气污染问题。该地区工业活动密集、机动车保有量持续攀升，加上复杂的地形和气象条件，使得大气污染形势愈发严峻。近年来，珠三角地区频繁出现灰霾天气，严重影响了区域的空气质量和居民的生活质量。据相关报道，2024年4月14-16日，受北方沙尘暴影响，广州、佛山、江门等珠三角城市出现明显灰霾，空气污染级别达到中度，对当地居民的身体健康和日常生活造成了极大困扰。持久性有机污染物（PersistentOrganicPollutants，POPs）作为一类具有特殊性质的污染物，在大气环境中广泛存在，对生态环境和人体健康构成了潜在威胁。POPs具有环境持久性、生物蓄积性、远距离环境迁移潜力和高毒性等特点，能在环境中长期残留，并通过食物链在生物体内不断富集，对生物体的内分泌、生殖、免疫等系统产生不良影响。例如，二恶英类物质是典型的POPs，其毒性极强，其中最毒者的毒性相当于氰化钾的1000倍以上，每人每日能容忍的二恶英摄入量为每公斤体重1pg。长期暴露于POPs污染的环境中，人类可能会面临癌症、生殖发育异常、神经和免疫功能受损等健康风险。研究珠三角城市群大气POPs的污染特征及暴露风险，对于深入了解该地区的大气污染状况，评估其对人体健康的潜在威胁具有重要意义。通过掌握POPs的污染水平、分布规律和来源，能够为制定针对性的污染防控策略提供科学依据，有助于减少POPs的排放，降低其对环境和人体健康的危害，保护生态环境和居民的身体健康，促进珠三角城市群的可持续发展。1.2持久性有机污染物概述持久性有机污染物（POPs）是指由人工合成，具有环境持久性、生物蓄积性、远距离环境迁移潜力，并对人体健康或生态环境产生不利影响的有机污染物。POPs具有以下显著特性：高持久性：POPs对生物降解、光解、化学分解作用有较高的抵抗能力，在环境中难以分解。例如，二恶英系列物质在气相中的半衰期为8-400天，水相中为166天到2119年，在土壤和沉积物中约17年到273年，这使得它们能够在环境中长期存在。生物蓄积性：由于其具有低水溶性、高脂溶性的特点，POPs能在生物体脂肪组织中进行生物积累，随着食物链的传递，浓度不断升高，最终在处于食物链顶端的人类和动物体内达到中毒的浓度。远距离迁移能力：POPs一般是半挥发性物质，在室温下就能挥发进入大气层，能通过蒸发作用在大气环境中远距离迁移，导致全球范围的污染传播。它们能从水体或土壤中以蒸气形式进入大气环境或者附在大气中的颗粒物上，在大气环境中远距离迁移而不会全部被降解，在一定条件下又沉降下来，然后又在某些条件下挥发，如此反复，使得POPs能够分散到地球上各个地方，甚至在北极圈等远离污染源的地方都发现了POPs污染。高毒性：POPs大都具有“三致（致癌、致畸、致突变）”效应，对生物体的健康危害极大。如苯并[a]芘、二恶英类化合物已被国际癌症研究机构（IARC）定义为一类致癌物。常见的POPs种类繁多，首批列入《关于持久性有机污染物的斯德哥尔摩公约》受控名单的POPs共有12种，包括7种有意生产的有机氯杀虫剂，如滴滴涕（DDT）、氯丹、灭蚁灵等；2种工业化学品，即六氯苯和多氯联苯；还有2种无意排放的副产品，如二恶英（多氯二笨并-p-二恶英）和呋喃（多氯二笨并呋喃）。截至2009年5月9日，又有林丹等九种化学品被新增列为POPs，至此，共有21种化学物质被正式认定为持久性有机污染物。除了上述物质，农业中的农药（如有机氯杀虫剂）、工业中的化学品（如六溴联苯等）以及工业副产物（如五氯苯等）也都属于POPs的范畴。POPs对环境和人体健康的危害是多方面的。在环境方面，POPs会破坏生态系统的平衡，导致陆地和水域生物的死亡，影响食物链，危及自然或人工生态系统，并增加火灾、洪水等自然灾害的风险。对人体健康而言，POPs会抑制免疫系统的正常反应、降低生物体的病毒抵抗能力，干扰内分泌系统，导致生物体暴露于POPs时出现生殖障碍、先天畸形、机体死亡等现象。国际社会高度重视POPs的管控问题，《关于持久性有机污染物的斯德哥尔摩公约》应运而生。该公约旨在淘汰和削减POPs的生成和排放，保护人类健康和环境免受POPs的危害。公约于2004年5月17日正式在全球生效，截至2005年5月，已有151个国家或组织签署，其中98个国家或组织已正式批准。我国是公约的正式缔约方，于2001年5月23日首批签署公约，2004年11月11日，公约对我国正式生效。20年来，中国已全面淘汰29种类持久性有机污染物的生产、使用和进出口，每年避免数十万吨持久性有机污染物的生产和环境排放，有效防范了相关农产品、消费品中持久性有机污染物的健康风险。此外，欧盟、新西兰等也根据自身情况，对一些POPs物质实施了新的限制措施，以减少其对环境和人类健康的危害。1.3大气被动采样技术大气被动采样技术是基于气体分子扩散或渗透原理采集空气中气态或蒸汽态污染物的一种采样方法。其原理是环境空气中的污染物分子以扩散方式传质到吸收膜，并采集在吸收膜上。常见的被动采样装置有扩散式采样器（如扩散管、被动式气体采样器）等。被动采样装置具有诸多优点。首先，其结构简单，体积小、重量轻，便于携带和操作，例如一些小型的扩散式采样器，体积小巧，方便在野外、偏远地区等复杂环境中使用；其次，被动采样器无需外部动力，不需要抽气泵和电源，操作简便，无需特别的维护，成本较低，适合进行大面积调查和长时间监测；此外，被动采样技术能够反映采样期间污染物的平均浓度，对于研究污染物的长期变化趋势具有重要意义。然而，被动采样装置也存在一定的缺点，主要表现为采样速度慢，对采样环境的要求较高。在一些需要快速获取污染物浓度数据的场合，被动采样技术可能无法满足需求。同时，如果采样环境的温度、湿度、风速等条件不稳定，可能会影响污染物分子的扩散速率，从而导致采样结果的准确性受到一定影响。在珠三角城市群的研究中，大气被动采样技术具有较高的适用性和优势。该地区地域广阔，城市众多，采用被动采样技术可以在不同地点同时进行采样，实现大面积的监测，有助于全面了解该地区大气POPs的污染状况。而且，被动采样器无需复杂的设备和电源，能够适应珠三角地区复杂的地形和环境条件，便于在城市、乡村、工业区域等不同功能区进行布置，为研究POPs的空间分布特征提供了便利。此外，被动采样技术能够长时间连续采样，获取的污染物浓度数据更能反映该地区大气POPs的实际污染水平，为评估其对人体健康的长期暴露风险提供可靠依据。1.4研究目标与内容本研究的目标是全面揭示珠三角城市群大气POPs的污染特征，解析其来源，并准确评估其对人体健康的暴露风险，为该地区大气POPs的污染防控提供科学依据和合理建议。具体研究内容包括以下几个方面：样品采集与分析：在珠三角城市群不同功能区设置多个采样点，运用大气被动采样技术，采集大气中的POPs样品。对采集到的样品进行预处理后，采用先进的仪器分析方法，如气相色谱-质谱联用仪（GC-MS）等，准确测定样品中各类POPs的浓度和组成。污染特征研究：分析珠三角城市群大气POPs的浓度水平、空间分布特征以及时间变化规律。探讨不同功能区（如城市中心区、工业区、郊区等）POPs污染程度的差异，研究季节变化对POPs浓度的影响，从而全面掌握该地区大气POPs的污染特征。来源解析：运用多元统计分析方法（如主成分分析、聚类分析等）结合相关源清单数据，对珠三角城市群大气POPs的来源进行解析。确定主要的污染源类型，如工业排放、交通尾气、农业活动等，明确各污染源对POPs污染的贡献比例，为制定针对性的污染控制措施提供依据。暴露风险评估：根据大气POPs的浓度数据、人群活动模式以及相关的暴露参数，采用合适的暴露评估模型，评估珠三角城市群居民通过呼吸途径对POPs的暴露水平。结合POPs的毒性数据，运用风险评估模型，对暴露于POPs的健康风险进行定量评估，确定不同人群（如儿童、成年人等）的潜在健康风险程度。防控建议：基于对珠三角城市群大气POPs污染特征、来源和暴露风险的研究结果，从政策法规、产业结构调整、污染治理技术等方面提出针对性的防控建议，为减少该地区大气POPs污染、保护居民健康提供科学指导。二、研究区域与方法2.1研究区域概况珠三角城市群位于中国广东省中南部、珠江入海口处，包括广州、深圳、佛山、珠海、东莞、中山、惠州、江门、肇庆九个城市。其地理位置优越，地处北纬21°17′36″-23°55′54″，东经111°59′42″-115°25′18″之间，濒临南海，毗邻港澳，是中国参与经济全球化的主体区域。珠三角城市群是中国人口集聚度最高、创新能力最强、综合实力最强的三大区域之一，截至2024年，该地区常住人口总量超过7000万，人口密度大，城市人口占比较高。例如，广州市常住人口超过1800万，深圳常住人口也在1700万左右，人口的高度集中使得该地区对资源和能源的需求巨大，给环境带来了较大压力。在经济发展方面，珠三角城市群是中国重要的经济区，经济总量庞大且增长迅速。2023年，珠三角地区GDP总量达到11.6万亿元，约占全国GDP总量的9.5%，经济增速保持在较高水平，高于全国平均增速。该地区形成了以制造业为主导，服务业快速发展的产业格局。制造业涵盖了电子信息、电器机械、石油化工、汽车制造等多个领域，如深圳的电子信息产业，汇聚了众多知名企业，产品畅销国内外；广州的汽车制造业也颇具规模，拥有多个汽车生产基地。服务业方面，金融、物流、会展、旅游等行业发展迅猛，广州和深圳是区域内的金融中心，拥有众多金融机构和完善的金融市场体系。珠三角城市群的工业布局呈现出多样化的特点。在广州、深圳、佛山等城市，形成了多个大型工业园区和产业集群，如广州经济技术开发区、深圳高新技术产业园区、佛山顺德家电产业集群等，这些园区和集群集中了大量的工业企业，产业集聚效应明显，有利于资源共享和产业协同发展，但同时也导致工业污染物排放相对集中。此外，东莞、中山等地以制造业为主，工业企业数量众多，分布较为广泛，涵盖了家具、服装、玩具等多个行业，这些企业在推动当地经济发展的同时，也带来了一定的环境污染问题。珠三角地区属于南亚热带海洋性季风气候，气候温暖湿润，年平均气温在22℃左右，年平均降水量在1600-2000毫米之间。这种气候条件有利于污染物的扩散，但在特定的气象条件下，如静稳天气、逆温层等，污染物容易积聚，导致空气质量恶化。此外，该地区夏季受台风影响较大，台风带来的强风和降雨可能会对大气污染物的分布产生影响，一方面强风可以加速污染物的扩散，另一方面降雨可以通过湿沉降作用去除部分污染物。这些因素对大气POPs污染产生了多方面的影响。人口和经济的快速增长导致能源消耗增加，工业生产、交通运输等活动排放的POPs也相应增多。工业布局的集中和分散特点使得POPs的排放源分布不均，工业园区周边的POPs浓度可能相对较高。气候条件则影响着POPs的扩散、迁移和转化过程，温暖湿润的气候有利于POPs在大气中的光化学反应和水解反应，而特殊的气象条件则可能导致POPs在局部地区积累，增加了大气POPs污染的复杂性和不确定性。2.2样品采集采样点的选择遵循代表性、均匀性和可行性的原则。考虑到珠三角城市群不同的功能区，包括城市中心区、工业区、郊区等，在每个功能区内选取具有代表性的地点作为采样点。例如，在城市中心区，选择人口密集、交通繁忙的区域，如广州市天河区的商业中心、深圳市福田区的金融中心等地，以反映城市核心区域的大气POPs污染状况；在工业区，选取大型工业园区内或周边的地点，如佛山顺德家电产业集群内的工厂附近、惠州大亚湾石化工业区周边，以监测工业活动对大气POPs的影响；在郊区，选择远离工业污染源和交通干线的区域，如肇庆市的农村地区、江门市的郊外，用于对比城市和工业区域，了解自然背景下的POPs污染水平。最终在珠三角城市群共设置了20个采样点，覆盖了九个城市，形成了较为全面的采样网络，以确保能够准确反映该地区大气POPs的空间分布特征。采样时间为2023年1月至2023年12月，为期一年，每月采样一次。这样的采样时间安排能够充分考虑到季节变化对大气POPs浓度的影响，因为不同季节的气象条件、工业生产活动和居民生活方式等因素存在差异，可能导致POPs的排放、扩散和转化过程发生变化。例如，夏季气温高，大气对流活动强烈，有利于POPs的扩散；冬季可能出现静稳天气，不利于污染物的扩散，导致POPs浓度升高。每月采样一次可以获取不同季节的POPs浓度数据，从而更全面地了解其时间变化规律。使用的被动采样装置为扩散式采样器，其工作原理是基于气体分子的扩散作用。采样器主要由采样膜、扩散腔和外壳组成，环境空气中的POPs分子通过扩散作用穿过扩散腔，被采样膜吸附，从而实现对大气中POPs的采集。在安装采样器时，选择空旷、通风良好且周围无明显障碍物的位置，将采样器固定在离地面高度约1.5-2米的支架上，以确保采样器能够充分接触大气，获取具有代表性的样品。同时，避免采样器靠近污染源、建筑物墙壁或树木等，防止局部环境因素对采样结果产生干扰。在采样期间，定期对采样器进行检查和维护，确保其正常运行，如检查采样器是否有损坏、采样膜是否正常吸附等。每月采样结束后，及时更换采样膜，并将采集到的样品妥善保存，尽快送往实验室进行分析。2.3样品分析样品前处理包括提取、净化和浓缩等过程。提取过程采用加速溶剂萃取（ASE）技术，该技术具有高效、快速、溶剂用量少等优点。将采集的样品放入萃取池中，加入适量的有机溶剂（如正己烷-丙酮混合溶剂），在一定的温度和压力下进行萃取，使POPs从样品中转移到有机溶剂中。萃取条件为温度80℃，压力1500psi，静态萃取时间5分钟，循环次数3次。净化过程使用硅胶柱层析法，以去除提取液中的杂质和干扰物质。将提取液通过装有硅胶的层析柱，利用硅胶对不同物质的吸附能力差异，使POPs与杂质分离。先用正己烷淋洗层析柱，去除大部分非极性杂质，然后用正己烷-二氯甲烷混合溶剂（体积比为1:1）洗脱，收集含有POPs的洗脱液。浓缩采用旋转蒸发仪和氮吹仪相结合的方法。将净化后的洗脱液转移至旋转蒸发瓶中，在40℃的水浴温度下，通过旋转蒸发仪初步浓缩，去除大部分有机溶剂。然后将浓缩液转移至氮吹管中，用氮气吹干，最后用正己烷定容至1毫升，供仪器分析使用。使用气相色谱-质谱联用仪（GC-MS）对样品中的POPs进行分析。GC-MS能够实现对复杂有机化合物的分离和定性、定量分析。分析的POPs种类包括多环芳烃（PAHs）、多氯联苯（PCBs）、有机氯农药（OCPs）等常见的持久性有机污染物。在分析过程中，采用选择离子监测模式（SIM），以提高检测的灵敏度和选择性。通过与标准物质的保留时间和质谱图进行对比，确定样品中POPs的种类和含量。为保证分析结果的准确性和可靠性，采取了一系列质量控制措施。在每批样品分析时，同时进行空白样品、加标回收样品和平行样品的分析。空白样品用于检测实验过程中是否存在污染，加标回收样品用于评估分析方法的准确性，平行样品用于检验分析结果的重复性。加标回收率控制在70%-120%之间，平行样品的相对偏差小于15%，确保分析数据的质量符合要求。2.4数据处理与分析方法使用Excel和SPSS软件进行数据处理和统计分析。Excel主要用于数据的录入、整理和初步计算，如计算平均值、标准差等基本统计量。SPSS软件则用于更深入的统计分析，包括描述性统计、相关性分析、主成分分析等。描述性统计用于分析大气POPs浓度的基本特征，包括均值、中位数、最大值、最小值、标准差等，以了解POPs浓度的总体水平和分布范围。相关性分析用于研究不同POPs之间以及POPs与其他环境因素（如气象参数、工业排放量等）之间的相关性，通过计算相关系数，判断变量之间的线性关系强度和方向，为进一步探讨POPs的来源和迁移转化机制提供依据。主成分分析（PCA）是一种常用的多元统计分析方法，用于对多个变量进行降维处理，提取主要的信息成分。在本研究中，将不同种类的POPs浓度作为变量，运用PCA方法，找出影响POPs分布的主要因素，对POPs的来源进行初步解析。通过PCA分析，可以将多个复杂的变量转化为少数几个主成分，这些主成分能够解释大部分数据的方差，从而简化数据结构，更清晰地揭示数据背后的潜在信息。除了主成分分析，还采用正定矩阵因子分解模型（PMF）对POPs的来源进行进一步解析。PMF模型是一种基于受体模型的源解析方法，能够定量地确定不同污染源对环境中污染物的贡献比例。通过输入样品中POPs的浓度数据以及相关的不确定性信息，PMF模型可以计算出各个污染源的因子贡献和因子谱，从而识别出主要的污染源类型，如工业排放源、交通源、农业源等，并确定它们对大气POPs污染的相对贡献大小。暴露风险评估采用美国环境保护署（USEPA）推荐的暴露评估模型和风险评估模型。通过计算居民通过呼吸途径对大气POPs的暴露剂量，结合POPs的毒性数据，评估其对人体健康的潜在风险。暴露剂量的计算考虑了大气POPs的浓度、居民的呼吸速率、暴露时间等因素。风险评估模型则根据暴露剂量和POPs的参考剂量（RfD）或致癌斜率因子（CSF），计算出非致癌风险（危害商，HQ）和致癌风险（CR），以评估不同人群（如儿童、成年人）暴露于POPs的健康风险程度。当HQ小于1时，认为非致癌风险较低；当CR在10-6-10-4之间时，认为致癌风险处于可接受范围，否则需要采取相应的风险控制措施。三、珠三角城市群大气POPs污染特征3.1浓度水平对2023年珠三角城市群20个采样点每月采集的大气样品进行分析，得到不同采样点和时间段大气中POPs的浓度数据（见表1）。结果显示，珠三角城市群大气中POPs的总浓度范围为10.2-85.6ng/m³，平均值为35.8ng/m³。其中，工业区采样点的POPs浓度普遍较高，平均值达到45.6ng/m³，显著高于城市中心区（32.5ng/m³）和郊区（28.7ng/m³）。例如，佛山顺德家电产业集群内某采样点的POPs浓度在2023年7月达到了85.6ng/m³，为所有采样点中的最高值，这可能与该地区密集的工业生产活动有关，大量的工业废气排放导致了POPs在大气中的积累。采样点类型浓度范围（ng/m³）平均值（ng/m³）工业区25.6-85.645.6城市中心区15.3-65.232.5郊区10.2-45.828.7在时间分布上，POPs浓度呈现出明显的季节性变化。冬季（12月-2月）的平均浓度最高，达到42.5ng/m³，夏季（6月-8月）的平均浓度最低，为28.6ng/m³。冬季POPs浓度升高可能是由于冬季气象条件不利于污染物扩散，如静稳天气增多，逆温层出现频率增加，使得大气中的POPs难以扩散稀释，从而导致浓度升高。此外，冬季居民取暖等活动可能增加了POPs的排放源。与其他地区的浓度水平进行对比，珠三角城市群大气POPs的平均浓度高于国内一些经济相对欠发达地区，如西部地区的某些城市，这些地区大气POPs平均浓度一般在10-20ng/m³左右。但与国内一些工业发达且污染较为严重的地区相比，如京津冀部分城市，珠三角城市群的POPs浓度处于中等水平，京津冀部分城市大气POPs平均浓度可达40-60ng/m³。这表明珠三角城市群大气POPs存在一定程度的污染，但通过与其他地区的对比，也为该地区进一步评估和控制POPs污染提供了参考依据。总体而言，珠三角城市群大气POPs的污染程度不容忽视，需要采取有效的措施加以控制和治理。3.2组成特征分析不同种类POPs在总浓度中的占比，结果如图1所示。多环芳烃（PAHs）在总浓度中的占比最高，达到45.6%，其次是多氯联苯（PCBs），占比为28.3%，有机氯农药（OCPs）占比为26.1%。PAHs主要来源于化石燃料的不完全燃烧，如工业锅炉、机动车尾气排放等。在珠三角城市群，工业活动频繁，机动车保有量巨大，这些因素导致PAHs的排放量大，从而在大气POPs组成中占比较高。PCBs曾被广泛应用于电力设备、塑料增塑剂等领域，虽然目前已被限制使用，但由于其环境持久性，在大气中仍有一定的残留。OCPs主要来源于过去农业生产中的使用，尽管近年来使用量逐渐减少，但历史残留以及部分非法使用情况使得其在大气中依然存在。图1不同种类POPs在总浓度中的占比不同采样点和时间段POPs组成存在一定差异。在工业区，PCBs的占比相对较高，可达35.0%左右，这是因为工业区内存在一些早期使用PCBs的工业设施，其排放或泄漏导致PCBs在该区域大气中的浓度相对较高。而在城市中心区，PAHs的占比更为突出，可达到50.0%以上，主要是由于城市中心区交通拥堵，机动车尾气排放量大，机动车尾气中含有大量的PAHs。在时间变化上，夏季OCPs的占比相对较高，这可能与夏季农业活动相对频繁，农药使用量增加有关；冬季PAHs的占比略有升高，可能是由于冬季取暖导致化石燃料燃烧量增加，从而排放更多的PAHs。这些组成差异反映了不同采样点和时间段污染源的变化以及环境因素对POPs迁移转化的影响。3.3气粒分配气粒分配是指大气中污染物在气态和颗粒态之间的分配过程，其分配情况受到多种因素的影响。POPs在大气中的气粒分配主要取决于其物理化学性质，如蒸汽压、溶解度等，同时也受到环境因素的影响，如温度、湿度、颗粒物浓度等。一般来说，蒸汽压较低、溶解度较小的POPs更容易吸附在颗粒物表面，以颗粒态存在；而蒸汽压较高、溶解度较大的POPs则更多地以气态形式存在。对珠三角城市群大气中POPs的气粒分配情况进行分析，结果表明，PAHs中分子量较大、蒸汽压较低的化合物，如苯并[a]芘等，主要以颗粒态存在，在颗粒态中的占比可达70%-80%；而分子量较小、蒸汽压较高的PAHs，如萘等，在气态和颗粒态中的分配相对较为均匀，气态占比约为40%-60%。PCBs和OCPs也呈现出类似的规律，低氯代的PCBs和挥发性较强的OCPs，如林丹等，气态占比较高；高氯代的PCBs和挥发性较弱的OCPs则更多地以颗粒态存在。进一步探讨其与温度、湿度等环境因素的关系，发现温度升高时，POPs的挥发性增强，气态占比增加。例如，在夏季高温时，一些原本以颗粒态存在的POPs会挥发进入气态，导致气态POPs的比例上升。湿度对气粒分配也有一定影响，湿度增加时，颗粒物表面的水分含量增多，可能会促进POPs在颗粒物表面的吸附，使得颗粒态POPs的占比增加。此外，颗粒物浓度的变化也会影响POPs的气粒分配，当大气中颗粒物浓度较高时，POPs与颗粒物碰撞吸附的机会增加，从而导致颗粒态POPs的占比升高。3.4季节变化不同季节大气中POPs浓度和组成呈现出明显的变化。浓度方面，如前文所述，冬季平均浓度最高，夏季最低，春秋季介于两者之间。冬季由于气象条件不利于污染物扩散，且居民取暖等活动增加了污染源，导致POPs在大气中积累，浓度升高；夏季大气对流活动强烈，有利于污染物扩散，同时降雨较多，通过湿沉降作用去除了部分POPs，使得浓度降低。组成上，夏季OCPs的占比相对较高，这与夏季农业活动中农药使用量增加有关。随着气温升高，土壤中的农药挥发进入大气，导致大气中OCPs浓度上升，在POPs组成中的占比增加。冬季PAHs的占比相对较高，主要是因为冬季取暖需求增加，煤炭、天然气等化石燃料的燃烧量增大，而化石燃料不完全燃烧会产生大量的PAHs，从而使得PAHs在POPs组成中的比例升高。季节变化对POPs污染特征的影响与污染源排放和气象条件密切相关。污染源排放方面，不同季节人类活动和工业生产模式存在差异，如冬季的取暖活动、夏季的农业活动等，这些差异导致不同季节POPs的排放源和排放量不同，进而影响其污染特征。气象条件方面，温度、湿度、风速、降水等气象因素在不同季节有明显变化，这些因素通过影响POPs的扩散、迁移、转化和沉降等过程，对POPs的污染特征产生重要影响。例如，风速较大时有利于POPs的扩散，降低其在局部地区的浓度；降水则通过湿沉降作用去除大气中的POPs，减少其在大气中的含量。因此，深入了解季节变化与污染源排放、气象条件等因素的关系，对于准确把握珠三角城市群大气POPs的污染特征，制定有效的污染控制措施具有重要意义。四、珠三角城市群大气POPs来源解析4.1定性分析方法采用指纹图谱法对珠三角城市群大气POPs的来源进行初步分析。指纹图谱法是通过分析POPs中各种化合物的相对比例和特征化合物的存在情况，与已知污染源的指纹图谱进行对比，从而推断其可能的来源。例如，多环芳烃（PAHs）的指纹图谱中，不同环数PAHs的相对含量以及特定PAHs的存在（如苯并[a]芘等）可以反映其来源。当大气中PAHs的指纹图谱与机动车尾气排放的特征指纹图谱相似，即低环数PAHs（2-3环）相对含量较高，且存在与机动车尾气相关的特征PAHs时，可初步推断机动车尾气是该区域PAHs的一个重要来源。在珠三角城市群，由于机动车保有量巨大，交通拥堵现象较为常见，机动车尾气排放成为大气PAHs的重要来源之一。特征比值法也是常用的定性分析方法。该方法通过计算POPs中某些特征化合物的比值，如荧蒽/芘（Flu/Pyr）、苯并[a]蒽/（苯并[a]蒽+屈）（BaA/（BaA+Chr））等，来判断其来源。不同来源的POPs具有不同的特征比值范围。一般来说，当Flu/Pyr比值在1.0-1.5之间时，表明PAHs主要来源于机动车尾气排放；当比值大于1.5时，可能来源于煤的燃烧。在珠三角地区的一些城市中心区，通过对大气中PAHs的特征比值计算，发现Flu/Pyr比值大多在1.0-1.3之间，这进一步证实了机动车尾气在该区域PAHs污染中的重要贡献。此外，结合研究区域的工业布局、农业活动以及交通状况等信息，综合判断POPs的可能来源。珠三角城市群工业发达，分布着众多电子、化工、建材等企业。这些工业企业在生产过程中可能会排放大量的POPs，如电子企业在电路板制造过程中可能会释放多氯联苯（PCBs），化工企业在有机合成过程中可能产生二恶英类物质等。农业活动方面，虽然目前珠三角地区农业用地相对减少，但在一些郊区仍存在农业生产，有机氯农药（OCPs）的使用或历史残留可能导致其在大气中的存在。交通方面，除了机动车尾气排放PAHs外，船舶运输也是一个潜在的POPs排放源，特别是在靠近港口的区域，船舶燃油燃烧可能释放出多种POPs，如多环芳烃、有机氯农药等。通过对这些信息的综合分析，能够更全面地了解珠三角城市群大气POPs的可能来源。4.2定量分析模型运用正定矩阵因子分解模型（PMF）对珠三角城市群大气POPs的来源进行定量解析。PMF模型是一种基于受体模型的源解析方法，能够在不需要预先知道污染源成分谱的情况下，通过对受体样品中污染物浓度数据的分析，识别出主要的污染源类型，并计算出各污染源对污染物的贡献比例。PMF模型的基本原理是将受体样品中POPs的浓度矩阵X分解为两个非负矩阵，即因子贡献矩阵G和因子谱矩阵F，以及一个残差矩阵E，表达式为X=G×F+E。其中，因子贡献矩阵G表示各污染源对不同样品中POPs的贡献大小，因子谱矩阵F表示各污染源中POPs的组成特征，残差矩阵E表示模型计算值与实际测量值之间的差异。模型通过最小化目标函数Q来确定最佳的矩阵分解结果，目标函数Q的定义为Q=∑i=1n∑j=1m（eij/σij）²，其中eij是第i个样品中第j种POPs的残差，σij是第i个样品中第j种POPs浓度的不确定性。在应用PMF模型时，首先对采集到的大气POPs样品数据进行预处理，包括数据质量控制、异常值处理等，确保输入数据的准确性和可靠性。然后，根据数据特点和研究目的，确定模型的输入参数，如因子数、不确定性估计方法等。通过多次运行模型，调整参数，使模型的Q值达到最小且稳定，同时确保因子谱具有合理的物理意义和可解释性。例如，在对珠三角城市群大气PAHs的源解析中，经过多次调试，确定将因子数设定为4时，模型能够较好地收敛，且得到的四个因子分别对应于机动车尾气排放源、工业燃煤源、生物质燃烧源和石油化工源。除了PMF模型，还尝试运用多元线性回归模型对POPs的来源进行分析。多元线性回归模型是一种经典的统计分析方法，通过建立POPs浓度与可能的污染源变量（如工业排放量、机动车保有量、农业活动强度等）之间的线性关系，来估算各污染源对POPs的贡献。其基本模型表达式为Y=β0+β1X1+β2X2+…+βkXk+ε，其中Y是POPs的浓度，X1，X2，…，Xk是可能的污染源变量，β0，β1，β2，…，βk是回归系数，ε是随机误差。在实际应用中，收集珠三角城市群各城市的工业污染源排放数据、机动车保有量数据以及农业生产相关数据等，将这些数据作为自变量，大气POPs浓度作为因变量，进行多元线性回归分析。通过回归分析得到各污染源变量的回归系数，从而计算出各污染源对大气POPs的贡献比例。但多元线性回归模型需要预先确定污染源变量，且假设变量之间存在线性关系，在实际应用中可能存在一定的局限性，因此将其与PMF模型结果进行对比和验证，以提高源解析结果的可靠性。4.3主要来源及贡献通过定性分析和定量分析，确定珠三角城市群大气POPs的主要来源包括工业排放、交通尾气、农业活动以及生物质燃烧等。工业排放是大气POPs的重要来源之一，对POPs污染的贡献较大，平均贡献率可达40%左右。在工业排放源中，电子、化工、建材等行业是主要的排放源。电子行业在生产过程中使用的含PCBs的电子元件，以及化工行业的有机合成反应、建材行业的高温烧制过程等，都会产生和排放POPs。例如，在一些电子产业集中的区域，如深圳的电子工业园区，大气中PCBs的浓度明显高于其他区域，PMF模型分析结果显示该区域工业排放对PCBs的贡献率可达到60%以上。近年来，随着环保政策的加强和工业污染治理技术的提高，工业排放对POPs的贡献呈现出逐渐下降的趋势，但由于工业生产规模仍然较大，其对大气POPs污染的影响依然不容忽视。交通尾气也是大气POPs的主要来源，贡献率约为30%。珠三角地区机动车保有量持续增长，截至2023年，该地区机动车保有量已超过3000万辆，大量机动车在运行过程中排放的尾气含有多种POPs，如PAHs、OCPs等。尤其是在城市中心区和交通繁忙的路段，交通尾气排放对大气POPs的贡献更为突出。以广州市为例，在一些主要交通干道附近的采样点，交通尾气对PAHs的贡献率可达到40%-50%。随着机动车排放标准的不断提高和新能源汽车的推广应用，交通尾气对POPs的贡献有望逐步降低，但在短期内，由于传统燃油汽车仍占较大比例，交通尾气依然是大气POPs污染的重要来源。农业活动对大气POPs的贡献率约为15%。尽管珠三角地区农业生产规模相对较小，但历史上有机氯农药的大量使用，以及部分非法使用现象，使得土壤中残留的OCPs会挥发进入大气，造成大气污染。在一些郊区的农业种植区域，大气中OCPs的浓度相对较高，通过源解析发现农业活动对OCPs的贡献率可达30%-40%。随着对农药使用的严格管控和农业生产方式的转变，农业活动对大气POPs的贡献在逐渐减少，但由于农药残留的长期存在，其对大气环境的影响仍将持续一段时间。生物质燃烧也是不可忽视的来源，贡献率约为10%。在珠三角地区，生物质燃烧主要包括居民生活中的薪柴燃烧、农业废弃物焚烧等。在农村地区和一些城乡结合部，居民在冬季取暖或日常生活中会使用薪柴作为燃料，农业生产过程中产生的秸秆等废弃物也存在焚烧现象。生物质燃烧会释放出大量的PAHs等POPs。例如，在冬季农村地区，生物质燃烧对PAHs的贡献率可达到20%-30%。随着清洁能源的推广和农村生活方式的改善，生物质燃烧的规模逐渐减小，其对大气POPs的贡献也在降低。其他来源，如垃圾焚烧、建筑施工等，对大气POPs的贡献率相对较小，约为5%。垃圾焚烧过程中会产生二恶英等POPs，在一些垃圾焚烧厂周边，大气中相关POPs的浓度会有所升高；建筑施工过程中的扬尘以及使用的建筑材料等也可能释放少量的POPs。针对这些主要来源，采取有效的控制措施至关重要。对于工业排放源，应加强环境监管，严格执行污染物排放标准，推动企业采用清洁生产技术和污染治理设施，减少POPs的排放。例如，鼓励电子企业采用无PCBs的电子元件，化工企业优化生产工艺，提高污染物的去除效率。对于交通尾气，进一步提高机动车排放标准，加强在用车尾气检测和监管，加大新能源汽车的推广力度，减少传统燃油汽车的使用。在农业活动方面，加强对农药使用的管理，推广低毒、低残留的农药，严禁非法使用有机氯农药，同时加强对土壤中农药残留的修复和治理。对于生物质燃烧，加大清洁能源的普及力度，推广秸秆综合利用技术，减少生物质焚烧现象。通过这些措施的实施，有望降低珠三角城市群大气POPs的污染水平，改善区域大气环境质量。五、珠三角城市群大气POPs暴露风险评估5.1暴露途径分析大气POPs进入人体的主要暴露途径包括呼吸吸入、皮肤接触和饮食摄入。呼吸吸入是人体暴露于大气POPs的直接途径，人们在日常生活中，通过呼吸不断地将周围空气中的POPs吸入体内。尤其是在污染较为严重的区域，如工业区、交通繁忙的路段等，大气中POPs浓度较高，呼吸吸入的暴露剂量相应增加。例如，在珠三角城市群的一些工业聚集区，工人在日常工作中长时间暴露于高浓度的POPs环境中，通过呼吸吸入的POPs量远高于普通居民。皮肤接触也是一种暴露途径，虽然POPs通过皮肤吸收的速率相对较低，但在长期暴露的情况下，其累积效应不容忽视。皮肤接触暴露主要发生在人们直接接触被POPs污染的空气、土壤或水的过程中。例如，在进行户外活动时，皮肤会与大气中的POPs接触，部分POPs可能会通过皮肤的角质层渗透进入人体。对于一些从事农业生产或工业劳动的人群，由于工作环境的特殊性，他们的皮肤可能会更频繁地接触到含有POPs的物质，从而增加了皮肤接触暴露的风险。饮食摄入是人体暴露于POPs的重要途径之一，POPs可以通过食物链在生物体内不断富集，最终进入人体。例如，在珠三角地区，由于大气中的POPs沉降到土壤和水体中，农作物在生长过程中可能会吸收土壤中的POPs，从而导致农产品中含有一定量的POPs。此外，水体中的鱼类等水生生物也会富集POPs，人们食用这些受污染的农产品和水产品后，就会摄入POPs。研究表明，通过饮食摄入的POPs在人体总暴露剂量中占有相当大的比例。不同暴露途径的暴露剂量和风险程度存在差异。呼吸吸入的暴露剂量主要取决于大气中POPs的浓度、呼吸速率以及暴露时间。在污染严重的区域，呼吸吸入的POPs剂量较高，可能会对呼吸系统、心血管系统等造成直接损害，增加患呼吸道疾病、心血管疾病等的风险。皮肤接触的暴露剂量相对较小，但长期接触可能会对皮肤产生刺激、过敏等反应，甚至可能通过皮肤吸收进入血液循环系统，对全身健康产生潜在影响。饮食摄入的暴露剂量与食物中POPs的含量密切相关，长期摄入受POPs污染的食物可能会导致POPs在人体内蓄积，对肝脏、肾脏等器官功能产生损害，增加患癌症、内分泌失调等疾病的风险。总体而言，饮食摄入和呼吸吸入在人体暴露于大气POPs的过程中起着较为重要的作用，需要重点关注和评估。5.2暴露剂量计算暴露剂量的计算采用美国环境保护署（USEPA）推荐的暴露评估模型，该模型考虑了多种因素来准确估算人体对大气POPs的暴露剂量。对于呼吸吸入暴露剂量（EDIinh），计算公式为：EDIinh=C×IR×ET×EF×ED/(BW×AT)，其中C为大气中POPs的浓度（ng/m³），IR为呼吸速率（m³/d），ET为暴露时间（h/d），EF为暴露频率（d/a），ED为暴露持续时间（a），BW为体重（kg），AT为平均时间（d）。在珠三角城市群，根据相关研究和统计数据，成人的平均呼吸速率约为15m³/d，儿童的呼吸速率相对较低，约为10m³/d。暴露时间方面，假设成人每天在室外活动时间为8小时，儿童为6小时。暴露频率按照一年365天计算，暴露持续时间考虑为30年（假设评估人群从出生开始暴露30年）。平均时间根据不同的健康效应有所不同，对于非致癌效应，AT=ED×365d；对于致癌效应，AT=70×365d（假设人均寿命为70年）。皮肤接触暴露剂量（EDIskin）的计算公式为：EDIskin=C×SA×AF×ABS×ET×EF×ED/(BW×AT)，其中SA为皮肤表面积（cm²），AF为皮肤附着系数（mg/cm²），ABS为皮肤吸收分数。成人的平均皮肤表面积约为18000cm²，儿童的皮肤表面积相对较小，约为10000cm²。皮肤附着系数一般取值为0.07mg/cm²，皮肤吸收分数根据不同的POPs种类有所差异，例如对于多环芳烃（PAHs），皮肤吸收分数约为0.01-0.1。考虑不同人群（如成人、儿童）的生理特征差异，儿童的呼吸速率、皮肤表面积等因素与成人不同，导致他们对大气POPs的暴露剂量也有所不同。由于儿童的呼吸速率相对较低，且在室外活动时间相对较短，因此通过呼吸吸入的POPs暴露剂量相对成人可能较低。然而，儿童的皮肤相对较薄，吸收能力较强，且单位体重的皮肤表面积较大，这使得他们通过皮肤接触暴露的剂量可能相对较高。此外，儿童处于生长发育阶段，对污染物的敏感性更高，相同剂量的POPs可能对儿童健康产生更大的影响。因此，在评估大气POPs暴露风险时，需要分别计算不同人群的暴露剂量，以便更准确地评估其健康风险。以珠三角城市群某城市中心区为例，根据监测数据，该区域大气中PAHs的平均浓度为20ng/m³。假设成人的体重为70kg，按照上述参数计算，成人通过呼吸吸入PAHs的年暴露剂量约为：EDIinh=20×15×8×365×30/(70×30×365)≈34.29ng/kg/d。对于儿童，体重假设为30kg，通过呼吸吸入PAHs的年暴露剂量约为：EDIinh=20×10×6×365×30/(30×30×365)≈40ng/kg/d。通过皮肤接触的暴露剂量计算，成人的EDIskin约为：EDIskin=20×18000×0.07×0.05×6×365×30/(70×30×365)≈10.8ng/kg/d（假设PAHs的皮肤吸收分数为0.05）；儿童的EDIskin约为：EDIskin=20×10000×0.07×0.05×6×365×30/(30×30×365)≈14ng/kg/d。由此可见，不同人群的暴露剂量存在差异，在暴露风险评估中需要充分考虑这些因素。5.3风险评估模型与指标风险评估使用危害商值（HQ）和致癌风险（CR）等指标来评估珠三角城市群居民暴露于大气POPs的健康风险。危害商值用于评估非致癌风险，其计算公式为：HQ=EDI/RfD，其中EDI为暴露剂量（mg/kg/d），RfD为参考剂量（mg/kg/d）。参考剂量是指人类终生暴露于该剂量水平的化学物质，预期不会对健康产生有害效应的日平均剂量估计值。不同POPs的参考剂量不同，例如，对于苯并[a]芘，其参考剂量为7×10-5mg/kg/d。当HQ小于1时，认为非致癌风险较低；当HQ大于1时，表明存在一定的非致癌风险，且HQ值越大，风险越高。致癌风险用于评估因暴露于致癌性POPs而导致癌症发生的可能性，计算公式为：CR=EDI×CSF，其中CSF为致癌斜率因子（mg/kg/d）-1，它表示单位暴露剂量下的致癌风险增加率。例如，苯并[a]芘的致癌斜率因子为7.3mg/kg/d。一般认为，当CR在10-6-10-4之间时，致癌风险处于可接受范围；当CR大于10-4时，致癌风险较高，需要采取相应的风险控制措施；当CR小于10-6时，致癌风险较低。风险评估的标准和阈值是判断风险程度的重要依据。国际上，如美国环境保护署（USEPA）等机构制定了一系列的风险评估标准和阈值，为全球范围内的风险评估提供了参考。在国内，相关部门也根据我国的实际情况，借鉴国际标准，制定了适用于我国的风险评估标准。例如，我国在环境健康风险评估中，对于POPs的暴露风险评估，也采用了与国际接轨的危害商值和致癌风险的评估指标和阈值。在珠三角城市群的大气POPs暴露风险评估中，严格遵循这些标准和阈值，以确保评估结果的科学性和可靠性，为后续的风险管理和决策提供准确的依据。5.4风险评估结果与分析通过计算珠三角城市群不同区域和人群的大气POPs暴露剂量，并结合风险评估模型和指标，得到该地区大气POPs的暴露风险评估结果。在不同区域方面，工业区的暴露风险普遍较高。以某大型工业区为例，该区域大气中POPs的浓度较高，通过呼吸吸入和皮肤接触的暴露剂量较大，导致危害商值和致癌风险均高于其他区域。经计算，该工业区成人的PAHs危害商值约为1.2，致癌风险约为2×10-4，表明该区域成人存在一定的非致癌风险和较高的致癌风险。而在郊区，由于大气POPs浓度相对较低，暴露剂量较小，风险水平相对较低，郊区成人的PAHs危害商值约为0.6，致癌风险约为8×10-5，处于可接受范围。不同人群的风险水平也存在差异。儿童由于其生理特征和行为模式，对大气POPs的暴露风险相对较高。如前文所述，儿童的皮肤相对较薄，吸收能力较强，且单位体重的皮肤表面积较大，通过皮肤接触暴露的剂量相对较高。同时，儿童的呼吸速率相对较低，但在室外活动时间相对较长，通过呼吸吸入的暴露剂量也不容忽视。以某城市中心区为例，儿童的PAHs致癌风险约为1.5×10-4，高于成人的致癌风险（约为1×10-4）。评估风险的可接受程度，根据风险评估指标的标准和阈值，对于危害商值大于1和致癌风险大于10-4的区域和人群，其风险不可接受，需要采取有效的风险管理措施。对于风险处于可接受范围的区域和人群，也需要持续监测，以防止风险水平上升。针对不同的风险水平，提出以下风险管理建议：对于高风险区域，如工业区，应加强对工业企业的监管，严格控制POPs的排放。推动企业采用清洁生产技术，改进生产工艺，减少POPs的产生。同时，加强对工业区周边环境的监测，及时掌握POPs的浓度变化情况，以便采取相应的措施。对于高风险人群，如儿童，应加强保护措施。例如，在污染天气时，尽量减少儿童的户外活动时间，教育儿童养成良好的卫生习惯，如勤洗手、洗脸等，减少皮肤接触暴露的风险。此外，加强对大气POPs污染的宣传教育，提高公众的环保意识和自我保护意识，让公众了解POPs的危害和预防措施，共同参与到大气污染防治工作中来。通过这些风险管理建议的实施，有望降低珠三角城市群大气POPs的暴露风险，保护居民的身体健康。六、结论与展望6.1主要研究结论本研究通过在珠三角城市群设置20个采样点，运用大气被动采样技术，对2023年全年的大气样品进行采集和分析，系统研究了该地区大气POPs的污染特征、来源解析及暴露风险评估，得出以下主要结论：污染特征：珠三角城市群大气中POPs的总浓度范围为10.2-85.6ng/m³，平均值为35.8ng/m³。工业区的POPs浓度显著高于城市中心区和郊区，冬季POPs浓度最高，夏季最低。在组成特征方面，多环芳烃（PAHs）占比最高，为45.6%，其次是多氯联苯（PCBs）和有机氯农药（OCPs）。不同功能区和季节的POPs组成存在差异，反映了污染源和环境因素的影响。POPs在大气中的气粒分配受其物理化学性质和环境因素的共同作用，分子量较大、蒸汽压较低的POPs更易以颗粒态存在，且温度、湿度等环境因素会影响其气粒分配比例。来源解析：通过指纹图谱法、特征比值法等定性分析方法，结合正定矩阵因子分解模型（PMF）和多元线性回归模型等定量分析方法，确定珠三角城市群大气POPs的主要来源包括工业排放（贡献率约40%）、交通尾气（贡献率约30%）、农业活动（贡献率约15%）和生物质燃烧（贡献率约10%）等。其中，电子、化工、建材等行业是工业排放的