哈尔滨市融雪径流中多环芳烃的污染格局与生态风险洞察

哈尔滨市融雪径流中多环芳烃的污染格局与生态风险洞察一、绪论1.1研究背景与目的随着城市化进程的加速，城市规模不断扩张，人口密度持续增加，大量的自然土地被开发为城市建设用地，城市地表硬化率显著提高。这种变化深刻地改变了原有的水文循环规律，城市地表面的融水径流已成为地表水体的重要污染源之一。融雪径流作为城市地表径流的一种特殊形式，在寒冷地区的春季尤为显著。当积雪融化时，大量的雪水携带各种污染物进入城市排水系统，最终汇入地表水体，对水环境质量和生态系统造成潜在威胁。多环芳烃（PAHs）是一类由两个或两个以上苯环以稠环形式相连的有机化合物，广泛存在于环境中。由于其具有致癌、致畸和致突变性，PAHs对人类健康和生态系统的潜在危害备受关注。在城市环境中，PAHs的污染源主要包括汽车尾气排放、燃油泄漏、轮胎摩擦、道路磨损以及工业废气排放等。这些PAHs主要通过大气干湿沉降进入城市环境，随后在融雪径流的冲刷作用下，大量进入城市地表水体，导致水体中PAHs的浓度升高，进而对水生生物和人类健康构成威胁。哈尔滨市作为我国北方的重要城市，冬季漫长且寒冷，降雪量大，融雪径流现象十分普遍。近年来，随着哈尔滨市的快速发展，城市交通流量不断增加，工业活动日益频繁，城市环境中的PAHs污染问题逐渐凸显。研究哈尔滨市融雪径流中PAHs的污染空间分布及生态风险评价，对于了解该地区PAHs的污染状况，评估其对生态系统和人类健康的潜在风险，制定有效的污染防控策略具有重要意义。通过本研究，可以为哈尔滨市的水环境管理和污染防治提供科学依据，促进城市的可持续发展。1.2国内外研究进展在融雪径流研究方面，国外起步较早，早期主要集中在融雪径流的形成机制和水文过程模拟。如美国学者在20世纪中叶就开始利用数学模型对融雪径流进行模拟研究，通过建立能量平衡模型，分析积雪融化过程中的能量交换和水分传输，为融雪径流的预测提供了理论基础。随着技术的发展，分布式水文模型逐渐应用于融雪径流研究，能够更准确地模拟不同地形和下垫面条件下的融雪径流过程。在数据监测方面，国外利用卫星遥感技术获取积雪覆盖面积和雪深等信息，结合地面气象站数据，实现对融雪径流的实时监测和动态分析。国内对融雪径流的研究相对较晚，但近年来发展迅速。在理论研究上，学者们深入探讨了融雪径流与气象因素、地形地貌、土地利用等之间的关系，揭示了融雪径流的影响因素和变化规律。在模型应用方面，国内引进和改进了多种国外先进的融雪径流模型，如SRM融雪径流模型、SWAT模型等，并结合国内实际情况进行参数率定和验证，提高了模型在我国不同地区的适用性。在监测手段上，国内也逐渐建立起地面监测站与遥感监测相结合的融雪径流监测体系，为研究提供了丰富的数据支持。在径流中多环芳烃污染研究领域，国外对PAHs的研究始于20世纪70年代，主要关注PAHs的来源、分布和迁移转化规律。通过源解析技术，明确了PAHs的主要来源包括化石燃料燃烧、工业活动和交通运输等。在分布研究方面，对不同环境介质如大气、水体、土壤和沉积物中的PAHs进行了广泛监测，分析了其空间分布特征和季节变化规律。在迁移转化研究中，探讨了PAHs在不同环境介质之间的分配和转化机制，以及生物降解、光降解等自然衰减过程。在生态风险评价方面，国外建立了较为完善的风险评价体系，采用多种评价方法，如危害熵值法、毒性当量法等，评估PAHs对生态系统和人类健康的风险。国内对径流中PAHs污染的研究始于20世纪90年代，早期主要集中在对大城市地表径流中PAHs的监测和污染特征分析。近年来，研究范围不断扩大，涵盖了不同类型城市和地区的地表径流，以及不同下垫面条件下的PAHs污染状况。在源解析方面，国内学者采用比值法、主成分分析法等多种方法，对PAHs的来源进行了深入研究，发现机动车排放、工业源和燃煤源是我国城市地表径流中PAHs的主要来源。在生态风险评价方面，国内借鉴国外经验，结合我国实际情况，开展了PAHs的生态风险评价研究，评估了其对水生生物和人体健康的潜在风险。尽管国内外在融雪径流及其中多环芳烃污染研究方面取得了一定进展，但仍存在一些不足。在融雪径流研究中，对于复杂地形和城市下垫面条件下融雪径流的模拟精度有待提高，不同模型之间的比较和验证还不够充分。在径流中多环芳烃污染研究方面，对PAHs在融雪径流中的迁移转化机制研究还不够深入，缺乏长期的监测数据和系统性的研究。此外，对于融雪径流中PAHs污染的生态风险评价，还需要进一步完善评价指标和方法，提高评价结果的准确性和可靠性。1.3研究内容与技术路线本研究旨在深入探究哈尔滨市融雪径流中多环芳烃的污染空间分布及生态风险，具体内容如下：融雪径流量空间分布特征分析：通过收集哈尔滨市的气象数据、地形数据和土地利用数据，利用水文模型模拟融雪径流过程，分析不同汇水区、排水口和下垫面条件下融雪径流量的空间分布规律，明确融雪径流的高流量区域和关键排放节点。融雪径流中多环芳烃污染空间分布特征研究：在哈尔滨市不同区域设置采样点，采集融雪径流样品，采用高效液相色谱-质谱联用仪等分析仪器测定样品中多环芳烃的含量和组成。分析不同环数、汇水区、排水口和下垫面条件下多环芳烃的污染空间分布特征，揭示多环芳烃在融雪径流中的分布规律和来源差异。融雪径流中多环芳烃污染源解析：运用比值法、主成分分析法等源解析方法，对融雪径流中多环芳烃的来源进行定性和定量分析，确定其主要污染源，为制定针对性的污染控制措施提供依据。融雪径流中多环芳烃污染生态风险评价：采用危害熵值（HQ）生态风险评价法、BaP毒性当量生态风险评价法（EBAP）和风险熵值（RQ）生态风险评价方法，综合评估融雪径流中多环芳烃对水生生物和人体健康的潜在生态风险，明确风险等级和风险区域，为风险管理提供科学参考。本研究的技术路线如下：首先，对研究区哈尔滨市的地理位置、气候条件、地形地貌和土地利用等概况进行全面了解，收集相关的基础数据。依据研究区的特点和研究目的，合理布设采样点，同时收集气象数据、地形数据和土地利用数据等。在数据处理方面，运用地理信息系统（GIS）等空间分析方法对数据进行处理和分析，利用水文模型计算径流量。对采集的融雪径流样品进行前处理，采用高效液相色谱-质谱联用仪等分析方法测定多环芳烃的含量，并严格进行质量控制与保证。通过源解析方法分析多环芳烃的来源，运用多种生态风险评价方法对其生态风险进行评价，最终对研究结果进行总结和讨论，提出相应的污染防控建议。具体技术路线流程如图1-1所示。[此处插入技术路线图]二、研究区概况与研究方法2.1研究区位置哈尔滨市地处中国东北，黑龙江省西南部，介于东经125°41'~130°13’、北纬44°03'~46°40’之间，是黑龙江省省会、副省级市，也是中国东北北部政治、经济、文化中心城市。特殊的历史进程和地理位置，造就了哈尔滨这座具有异国情调的美丽城市，它不仅荟萃了北方少数民族的历史文化，而且融合了中外文化，素有“冰城”“东方莫斯科”“东方小巴黎”之美称。全市土地面积5.31万平方公里，其中市区面积10198平方公里。2024年末全市户籍总人口933.0万人，其中城镇人口525.9万人。哈尔滨属于中温带大陆性季风气候，冬季漫长寒冷，夏季短暂炎热，四季分明。冬季（12月-翌年2月）平均气温-14.6℃，每年有5个月的月平均气温低于0℃，结冰期190天左右，降雪量大且积雪期长。这种气候特点使得融雪径流成为城市地表径流的重要组成部分，春季积雪融化时，大量雪水携带污染物进入城市排水系统，对城市水环境造成潜在威胁。哈尔滨作为北方重要城市，近年来城市规模不断扩大，交通流量持续增加，工业活动日益频繁。城市中大量的汽车尾气排放、工业废气排放以及道路磨损等，导致环境中多环芳烃的来源广泛且复杂。因此，研究哈尔滨市融雪径流中多环芳烃污染具有典型性和代表性，其研究结果对于了解北方寒冷地区城市融雪径流污染状况，以及制定针对性的污染防治措施具有重要的参考价值。2.2采样点布设为全面研究哈尔滨市融雪径流中多环芳烃的污染空间分布特征，在不同汇水区、排水口和下垫面设置了采样点。依据哈尔滨市的地形地貌、水系分布以及城市功能分区，将研究区域划分为多个汇水区，如松花江流域的各个支流汇水区。在每个汇水区内，选择具有代表性的排水口进行采样，这些排水口包括城市主干道排水口、居民区排水口、商业区排水口以及工业园区排水口等，以反映不同功能区域融雪径流的污染情况。针对不同下垫面，分别在城市道路、人行道、停车场、公园绿地、屋顶等位置设置采样点。在城市道路采样时，考虑到交通流量的差异，选择车流量大的主干道和车流量相对较小的次干道进行采样；人行道采样点分布在不同区域，包括商业区、居民区和学校周边等人流量较大的地方；停车场采样点选择露天停车场和地下停车场出口，以对比不同类型停车场融雪径流的污染特征；公园绿地采样点设置在不同植被覆盖类型的区域，如草坪、树林等；屋顶采样点选择不同建筑类型的屋顶，如居民楼屋顶、商业楼屋顶等。具体采样点位置利用全球定位系统（GPS）进行精确定位，并记录其经纬度信息。共设置了[X]个采样点，采样点分布均匀，能够全面覆盖研究区域，确保采集的融雪径流样品具有代表性，为后续准确分析多环芳烃的污染空间分布特征提供数据支持。采样点分布情况如图2-1所示。[此处插入采样点分布图]2.3数据来源及处理方法气象数据主要来源于哈尔滨市气象局，包括研究期间的气温、降水、风速、日照时数等数据，这些数据用于分析气象因素对融雪径流和多环芳烃污染的影响。土地利用数据则来自于地理国情监测云平台，通过对该平台提供的高分辨率遥感影像进行解译和分类，获取研究区域内不同土地利用类型的分布信息，如建设用地、耕地、林地、草地、水域等，以便分析土地利用类型与融雪径流中多环芳烃污染之间的关系。在对多环芳烃浓度等数据处理时，首先对原始数据进行审核，检查数据的完整性和准确性，剔除异常值。对于缺失的数据，采用线性插值法或根据周边采样点数据进行合理估算补充。然后，对处理后的数据进行统计分析，计算多环芳烃的浓度平均值、最大值、最小值、标准差等统计参数，以描述其浓度分布特征。运用统计软件SPSS对多环芳烃浓度与其他环境因素（如气象因素、土地利用类型等）进行相关性分析，探究各因素之间的内在联系。同时，将多环芳烃浓度数据与采样点的地理位置信息相结合，导入地理信息系统（GIS）软件中，进行空间可视化表达和分析，直观展示多环芳烃在研究区域内的空间分布规律。2.4空间分析方法本研究利用地理信息系统（GIS）强大的空间分析功能，对融雪径流中多环芳烃污染数据进行深入分析。其中，反距离权重插值法（IDW）是常用的空间插值方法之一，其原理是基于距离反比的权重分配思想，认为未知点的值受其周围已知点的影响，且距离越近影响越大。在进行插值时，首先确定每个采样点的位置坐标，然后根据设定的搜索半径和参与插值的采样点数量，计算未知点与周围已知采样点之间的距离。根据距离的倒数作为权重，对周围已知采样点的多环芳烃浓度值进行加权平均，从而得到未知点的多环芳烃浓度估计值。通过IDW插值法，可以将离散的采样点数据转化为连续的空间分布数据，生成多环芳烃浓度的空间分布图，直观地展示多环芳烃在研究区域内的空间变化趋势，明确高污染区域和低污染区域的分布范围。除反距离权重插值法外，还运用了缓冲区分析方法。根据研究目的，以河流、道路、污染源等地理要素为中心，设置一定宽度的缓冲区。通过分析缓冲区范围内多环芳烃浓度的变化，研究这些地理要素对多环芳烃污染分布的影响。例如，分析河流缓冲区范围内多环芳烃浓度与距离河流远近的关系，探究河流对多环芳烃的稀释和扩散作用；分析道路缓冲区范围内多环芳烃浓度与交通流量的关系，评估交通源对多环芳烃污染的贡献。此外，还利用了叠加分析功能，将多环芳烃浓度分布图层与土地利用类型图层、地形图层等进行叠加，综合分析不同土地利用类型和地形条件下多环芳烃的污染特征，深入探究多环芳烃污染的影响因素和空间分布规律。2.5径流量计算方法本研究采用融雪径流模型（SRM）来计算融雪径流量。该模型是一种基于能量平衡原理的分布式水文模型，能够考虑地形、气象、积雪等多种因素对融雪径流的影响。在模型中，将研究区域划分为多个网格单元，每个网格单元具有独立的地形、气象和下垫面属性。模型的输入参数包括气象数据（如气温、降水、辐射等）、地形数据（如海拔、坡度、坡向等）以及积雪参数（如积雪深度、雪密度等）。模型的核心计算过程如下：首先，根据能量平衡原理计算每个网格单元内积雪的消融量。积雪消融主要受太阳辐射、气温、风速等因素影响，通过能量平衡方程计算出净辐射通量、感热通量、潜热通量以及土壤热通量等能量项，进而得到积雪消融所需的能量，从而计算出积雪消融量。然后，考虑融雪水在地表的产流和汇流过程。融雪水一部分通过下渗进入土壤，一部分形成地表径流。下渗量根据土壤的渗透能力和前期土壤含水量等因素确定，地表径流量则通过水量平衡方程计算得到，即地表径流量等于积雪消融量减去下渗量和蒸发量。最后，将每个网格单元的地表径流量进行累加，得到整个研究区域的融雪径流量。在模型参数确定方面，通过收集研究区域内的历史气象数据、地形数据、积雪数据以及实测的融雪径流量数据，采用参数率定和验证的方法来确定模型的最优参数。利用试错法或自动优化算法对模型参数进行调整，使模型模拟的融雪径流量与实测值之间的误差最小。常用的评价指标包括相关系数（R）、决定系数（R²）、均方根误差（RMSE）等，通过比较模拟值与实测值在这些评价指标上的表现，确定模型参数的合理性和模型的可靠性。经过参数率定和验证后的融雪径流模型，能够较为准确地计算哈尔滨市不同区域的融雪径流量，为后续分析融雪径流中多环芳烃的污染负荷提供数据支持。2.6多环芳烃样品前处理及分析方法2.6.1前处理方法在融雪径流样品采集后，迅速将其转移至棕色玻璃瓶中，以避免光照对多环芳烃的影响。每个样品采集量为1L，采集后立即加入适量硫酸铜（CuSO₄），使水样中硫酸铜的浓度达到1g/L，目的是抑制微生物的生长，防止多环芳烃在水样中发生生物降解。将样品置于低温冷藏箱中，在4℃条件下保存，并尽快送回实验室进行分析。在实验室中，首先对样品进行萃取处理。采用液-液萃取法，向水样中加入50mL二氯甲烷（CH₂Cl₂），振荡萃取5min，使多环芳烃从水相转移至有机相。然后将萃取液转移至分液漏斗中，静置分层15min，使有机相和水相充分分离。收集下层有机相，再用50mL二氯甲烷对水相进行重复萃取一次，合并两次萃取得到的有机相。将合并后的有机相通过无水硫酸钠（Na₂SO₄）柱进行脱水处理，以去除有机相中残留的水分。无水硫酸钠柱的制备方法为：在玻璃层析柱中填充约5cm高的无水硫酸钠，并用少量二氯甲烷冲洗柱子，以确保柱子的通畅。将有机相缓慢通过无水硫酸钠柱，收集流出液，得到脱水后的有机萃取液。为进一步净化有机萃取液，采用硅胶柱层析法。硅胶柱的制备方法为：在玻璃层析柱中依次填充10g硅胶（100-200目）和5g无水硫酸钠，用正己烷（C₆H₁₄）预淋洗柱子，使硅胶充分活化。将脱水后的有机萃取液缓慢加入硅胶柱中，先用100mL正己烷淋洗柱子，去除萃取液中的非极性杂质，然后用100mL二氯甲烷和正己烷（体积比为1:1）的混合溶液淋洗柱子，收集洗脱液，该洗脱液中主要含有多环芳烃。将收集到的洗脱液在旋转蒸发仪上进行浓缩，温度控制在40℃左右，压力为0.08MPa，直至浓缩至约1mL。最后，用氮吹仪将浓缩液进一步浓缩至0.5mL，转移至进样小瓶中，待分析。2.6.2分析方法使用气相色谱-质谱联用仪（GC-MS）对处理后的样品进行多环芳烃分析。GC-MS型号为[具体型号]，色谱柱采用HP-5MS毛细管柱（30m×0.25mm×0.25μm）。进样口温度设定为280℃，采用不分流进样方式，进样量为1μL。载气为高纯氦气（纯度≥99.999%），流速为1.0mL/min。程序升温条件为：初始温度50℃，保持2min，以10℃/min的速率升温至280℃，保持10min。质谱条件如下：离子源为电子轰击源（EI），电子能量为70eV，离子源温度为230℃，四极杆温度为150℃。扫描方式为选择离子扫描（SIM），根据不同多环芳烃的特征离子进行定性和定量分析。在定性分析时，通过比较样品中各色谱峰的保留时间和标准物质中多环芳烃的保留时间，确定样品中多环芳烃的种类。在定量分析时，采用外标法，以不同浓度的多环芳烃标准溶液绘制标准曲线，根据样品中各多环芳烃的峰面积，在标准曲线上查找对应的浓度，从而计算出样品中多环芳烃的含量。2.6.3质量控制与保证为确保多环芳烃分析数据的准确性和可靠性，采取了一系列质量控制与保证措施。在每批样品分析过程中，进行加标回收率试验。选择一个代表性的样品，加入一定量的多环芳烃标准物质，按照与实际样品相同的前处理和分析方法进行处理和测定，计算加标回收率。加标回收率的计算公式为：加标回收率（%）=（加标样品测定值-样品测定值）/加标量×100%。本研究中，多环芳烃的加标回收率控制在70%-120%之间，若回收率超出此范围，则重新分析该样品或查找原因并进行整改。同时，每分析10个样品，进行一次空白试验。空白试验采用与样品相同的前处理步骤，但不加入实际样品，而是使用超纯水代替。通过空白试验，检测分析过程中是否存在试剂污染、仪器污染等情况。若空白试验中检测到多环芳烃的含量超过方法检出限，则对实验仪器、试剂和实验环境进行全面检查和清洁，直至空白试验结果合格。此外，定期对气相色谱-质谱联用仪进行维护和校准。使用标准物质对仪器的保留时间、峰面积等参数进行校准，确保仪器的性能稳定和分析结果的准确性。在样品分析过程中，每隔一段时间对标准物质进行测定，检查仪器的稳定性和重复性。若仪器的测定结果与标准值偏差较大，则对仪器进行重新校准和调试，保证分析数据的质量。2.7生态风险评价方法2.7.1危害熵值(HQ)生态风险评价法危害熵值（HQ）生态风险评价法是基于污染物的浓度和其对生物的毒性阈值来评估生态风险的方法。其计算公式为：HQ=C_{i}/T_{i}，其中C_{i}为第i种多环芳烃的实测浓度（μg/L），T_{i}为第i种多环芳烃对生物的毒性阈值（μg/L）。毒性阈值通常参考相关的生态毒理学研究数据，不同的多环芳烃对不同生物的毒性阈值有所差异。根据危害熵值的大小，将生态风险等级划分为低风险、中风险和高风险三个等级。当HQ\leq0.1时，认为生态风险为低风险，此时多环芳烃对生物的潜在危害较小；当0.1<HQ\leq1时，生态风险为中风险，表明多环芳烃可能对生物产生一定的影响；当HQ>1时，生态风险为高风险，意味着多环芳烃对生物具有较大的潜在危害，可能会对生态系统的结构和功能造成破坏。通过计算研究区域内不同采样点融雪径流中各多环芳烃的危害熵值，评估其生态风险等级，从而了解多环芳烃对当地生态系统的潜在威胁程度。2.7.2BaP毒性当量生态风险评价法(EBAP)BaP毒性当量生态风险评价法（EBAP）是以苯并[a]芘（BaP）为参照，将其他多环芳烃的浓度换算为相当于BaP的浓度，即毒性当量浓度（TEQ），以此来评价多环芳烃的综合生态风险。其计算公式为：TEQ_{i}=C_{i}\timesTEF_{i}，其中TEQ_{i}为第i种多环芳烃的毒性当量浓度（μg/L），C_{i}为第i种多环芳烃的实测浓度（μg/L），TEF_{i}为第i种多环芳烃相对于BaP的毒性当量因子。毒性当量因子是根据多环芳烃的致癌性和致突变性等毒性数据确定的，不同多环芳烃的毒性当量因子不同，BaP的毒性当量因子定义为1，其他多环芳烃的毒性当量因子相对较小。将所有多环芳烃的毒性当量浓度相加，得到总毒性当量浓度（\sumTEQ），即\sumTEQ=\sum_{i=1}^{n}TEQ_{i}。根据总毒性当量浓度的大小进行生态风险评价，当\sumTEQ\leq1μg/L时，生态风险较低；当1<\sumTEQ\leq10μg/L时，生态风险为中等；当\sumTEQ>10μg/L时，生态风险较高。通过计算研究区域融雪径流中多环芳烃的总毒性当量浓度，判断其生态风险水平，为评估多环芳烃对生态系统和人类健康的潜在风险提供依据。2.7.3风险熵值(RQ)生态风险评价方法风险熵值（RQ）生态风险评价方法是基于多环芳烃的总等效质量浓度（\sumC_{eq}）和毒性当量因子（TEF）来评估生态风险。首先计算总等效质量浓度，计算公式为：\sumC_{eq}=\sum_{i=1}^{n}C_{i}\timesTEF_{i}，其中C_{i}为第i种多环芳烃的实测浓度（μg/L），TEF_{i}为第i种多环芳烃相对于BaP的毒性当量因子。然后根据总等效质量浓度与风险评价标准进行比较，确定风险等级。一般将风险熵值划分为低风险（RQ<0.1）、中风险（0.1≤RQ<1）和高风险（RQ≥1）三个等级。当风险熵值小于0.1时，表明多环芳烃对生态系统的风险较低；当风险熵值在0.1到1之间时，生态风险处于中等水平，可能对生态系统产生一定影响；当风险熵值大于等于1时，生态风险较高，多环芳烃对生态系统具有较大的潜在威胁。通过计算研究区域融雪径流中多三、哈尔滨市融雪径流模型建立3.1模型创建本研究选用暴雨洪水管理模型（SWMM）来模拟哈尔滨市融雪径流过程。SWMM模型由美国环境保护署（EPA）开发，是一款广泛应用于城市排水系统模拟的动态降水径流模型，能够模拟城市地区单一降水事件或长期的水量和水质情况。该模型基于物理过程，由水文和水力学模拟模块组成，可模拟不同土地利用类型的雨水径流量、径流污染负荷以及污水收集和处理系统性能。在水文模块中，SWMM模型采用合理化公式、SCS曲线法等方法计算地表径流，考虑了降雨、蒸发、下渗、地表蓄水等水文过程。水力学模块则基于圣维南方程组，模拟水流在管网、渠道等水力设施中的流动，能够计算流速、流量、水深等水力参数。通过输入流域或城市区域的地理信息（如地形数据、土地利用类型等）和相关参数（如降雨参数、地表参数、管网参数等），SWMM模型可以输出雨水径流量、水质和污水系统性能预测结果，为城市排水系统的规划、设计、运行和管理提供科学依据。3.2气象数据输入气象数据是融雪径流模型的重要输入参数，直接影响模型的模拟精度。本研究收集了哈尔滨市多个气象站（如哈尔滨气象站、呼兰气象站等）在研究期间（[具体研究时间段]）的气象数据，包括气温、降水、风速、日照时数等。对于气温数据，采用距离加权平均法对不同气象站的气温进行插值处理，得到研究区域内各个网格单元的气温数据。因为不同气象站的气温可能存在差异，通过距离加权平均可以综合考虑各气象站对研究区域不同位置的影响，使得气温数据更能反映研究区域的实际情况。降水数据同样进行空间插值处理，采用克里金插值法，该方法基于区域化变量理论，考虑了降水在空间上的自相关性，能够更准确地估计研究区域内任意位置的降水量。在将处理后的气象数据输入SWMM模型时，按照模型的输入格式要求，将气温、降水等数据按时间序列进行整理。例如，将每日的平均气温、降水量等数据按照时间顺序依次输入到模型中，确保模型能够准确模拟融雪径流过程中气象因素的动态变化。同时，对输入数据进行质量检查，检查数据的完整性、合理性，避免因数据错误导致模型模拟结果偏差。3.3水文模块3.3.1雨量计雨量计数据是确定降雨输入的关键。本研究收集了哈尔滨市内多个雨量计的实测降雨数据，这些雨量计分布在不同区域，包括城市中心区、郊区以及不同地形地貌区域，以确保能够全面反映研究区域内的降雨情况。在数据处理过程中，对雨量计数据进行质量控制，检查数据的完整性和准确性，剔除异常值和错误数据。对于缺失的数据，采用线性插值法或根据周边雨量计数据进行估算补充。在SWMM模型中，将处理后的雨量计数据按照模型要求的格式进行输入。模型通过读取雨量计数据，获取降雨强度、降雨历时和降雨分布等信息，这些信息用于计算地表径流的产生和水量。例如，模型根据降雨强度和降雨历时，结合下垫面条件（如地表粗糙度、渗透率等），运用合理化公式或SCS曲线法计算地表径流量。3.3.2汇水区汇水区的划分对于准确模拟融雪径流至关重要。本研究利用地理信息系统（GIS）技术，结合哈尔滨市的数字高程模型（DEM）数据进行汇水区的划分。首先，对DEM数据进行预处理，去除噪声和异常值，提高数据的精度。然后，运用GIS的水文分析工具，如水流方向分析、汇流累积量计算等，确定水流的流向和汇流路径。根据水流的汇流路径，将研究区域划分为多个汇水区。在划分过程中，考虑了地形地貌、土地利用类型等因素，确保汇水区的划分合理。对于每个汇水区，设置相关参数，如面积、坡度、不透水率、曼宁糙率等。汇水区面积通过GIS计算得到；坡度根据DEM数据计算得出，反映了汇水区内地表的倾斜程度，影响水流的速度和汇流时间；不透水率根据土地利用类型进行估算，不同土地利用类型（如城市道路、绿地、屋顶等）具有不同的不透水率，影响地表径流的产生和下渗量；曼宁糙率则根据地表覆盖物的粗糙程度确定，影响水流在地表的流动阻力。3.3.3积雪和融雪在模型中，积雪和融雪过程的模拟采用度日因子法。该方法基于气温与积雪消融之间的线性关系，假设积雪消融速率与日平均气温成正比。模型中设置了融雪度日因子，该因子表示单位气温变化所引起的积雪消融量，其取值根据哈尔滨市的历史气象数据和积雪消融观测数据进行率定。当气温高于0℃时，积雪开始融化，融雪量根据度日因子和高于0℃的气温累积值计算得出。同时，考虑了积雪的蒸发和升华过程，通过设置蒸发系数和升华系数来反映这两个过程对积雪量的影响。蒸发系数和升华系数同样根据当地的气象条件和积雪特性进行确定。此外，模型还考虑了融雪水的下渗和地表径流的产生。下渗量根据土壤的渗透能力和前期土壤含水量等因素确定，采用霍顿下渗公式或格林-安普特下渗公式进行计算。地表径流量则通过水量平衡方程计算得到，即地表径流量等于积雪消融量减去下渗量和蒸发量。通过这些参数和公式的设置，模型能够较为准确地模拟哈尔滨市融雪径流过程中积雪和融雪的动态变化。3.4水力模块3.4.1节点和管段在SWMM模型中，节点和管段是构成城市排水系统的基本要素。节点代表排水系统中的特定位置，如检查井、雨水口、出水口等，管段则连接各个节点，用于输送水流。对于节点，定义了其地面高程、初始水深、蓄水容积等参数。地面高程根据实际测量数据或DEM数据确定，是计算节点处水流能量和水位的重要依据；初始水深反映了节点在模拟开始时的积水情况，根据实际观测或经验设定；蓄水容积表示节点能够容纳的最大水量，与节点的尺寸和结构有关。管段的定义包括长度、直径、坡度、粗糙系数等水力参数。管段长度通过实际测量或地图数据获取；直径根据排水系统的设计要求和流量计算确定，影响管段的输水能力；坡度是管段两端的高差与长度之比，决定了水流在管段中的流速和流向；粗糙系数反映了管段内壁的粗糙程度，影响水流的阻力，一般根据管段的材质和使用年限等因素确定。在模型中，通过这些节点和管段参数的设置，构建起城市排水系统的水力模型。模型基于圣维南方程组，采用有限差分法或有限体积法对管段中的水流进行数值模拟，计算管段中的流速、流量、水位等水力参数，从而模拟融雪径流在排水系统中的流动过程。3.5水质模块3.5.1污染物在本研究中，将多环芳烃作为主要污染物纳入SWMM模型的水质模块。针对多环芳烃，设置了一系列参数，包括初始浓度、累积系数、冲刷系数等。初始浓度根据前期对哈尔滨市融雪径流中多环芳烃的监测数据确定，反映了融雪径流在起始阶段多环芳烃的污染程度。累积系数用于描述多环芳烃在地表的累积过程，它表示单位时间内单位面积地表上多环芳烃的累积量，与土地利用类型、交通流量、大气沉降等因素有关。例如，在交通繁忙的城市道路区域，由于汽车尾气排放和轮胎磨损等原因，多环芳烃的累积系数相对较大；而在绿地等区域，累积系数相对较小。冲刷系数则用于模拟多环芳烃在融雪径流冲刷作用下从地表进入水体的过程，它反映了融雪径流对地表多环芳烃的冲刷能力，与径流速度、降雨强度、地表粗糙度等因素相关。通过调整这些参数，模型能够模拟不同条件下多环芳烃在融雪径流中的迁移转化过程。3.5.2土地利用不同土地利用类型对多环芳烃污染的影响在模型中主要通过累积系数和冲刷系数来体现。如前文所述，城市道路、停车场等交通密集区域，由于机动车尾气排放、燃油泄漏以及轮胎与路面的摩擦等，多环芳烃的累积量相对较高，因此设置较高的累积系数。同时，这些区域的地表通常较为硬化，在融雪径流的冲刷作用下，多环芳烃更容易被带入水体，所以冲刷系数也相对较大。而公园绿地、林地等区域，植被对大气中的多环芳烃具有一定的吸附和净化作用，且地表土壤具有较好的渗透性，能够截留部分多环芳烃，使得多环芳烃在这些区域的累积量较低，相应地设置较低的累积系数。此外，由于绿地和林地的地表粗糙度较大，融雪径流的流速相对较慢，对多环芳烃的冲刷能力较弱，故冲刷系数也较低。通过这种方式，模型能够反映不同土地利用类型对多环芳烃污染的影响，更准确地模拟多环芳烃在融雪径流中的污染过程和空间分布特征。3.6模型参数的率定与验证模型参数的率定和验证是确保模型准确性和可靠性的关键步骤。本研究收集了哈尔滨市部分区域的实测融雪径流量和多环芳烃浓度数据，用于模型参数的率定和验证。在参数率定过程中，采用试错法和自动优化算法相结合的方式。首先，根据经验和相关研究成果，给定模型参数的初始值。然后，运行模型，将模拟结果与实测数据进行对比。通过比较模拟的融雪径流量、多环芳烃浓度与实测值之间的差异，利用自动优化算法（如遗传算法、粒子群优化算法等）对模型参数进行调整，使模拟结果与实测数据的误差最小化。常用的评价指标包括相关系数（R）、决定系数（R²）、均方根误差（RMSE）等。相关系数反映了模拟值与实测值之间的线性相关程度，取值范围在-1到1之间，越接近1表示相关性越好；决定系数表示模型对数据的拟合优度，取值范围在0到1之间，越接近1说明模型的拟合效果越好；均方根误差衡量了模拟值与实测值之间的平均误差程度，其值越小表示模拟结果越接近实测值。经过多次参数调整和模拟计算，得到了一组最优的模型参数。将这组参数应用于模型中，对未参与率定的其他区域的融雪径流和多环芳烃污染进行模拟，并与实测数据进行验证。验证结果表明，模型模拟的融雪径流量和多环芳烃浓度与实测值之间具有较好的一致性，相关系数R均大于0.8，决定系数R²大于0.7，均方根误差RMSE在可接受范围内，说明模型能够较好地模拟哈尔滨市融雪径流过程和多环芳烃的污染情况。3.7本章小结本章通过选用SWMM模型，全面系统地建立了哈尔滨市融雪径流模型。详细介绍了模型创建过程，包括模型的原理和功能，使其能够准确模拟城市排水系统中的水文、水力和水质过程。在气象数据输入环节，通过对多个气象站数据的收集、处理和合理输入，为模型提供了可靠的气象条件信息。在水文模块中，对雨量计数据进行严格处理和应用，科学划分汇水区并合理设置相关参数，同时采用度日因子法准确模拟积雪和融雪过程。水力模块通过明确节点和管段的定义及相关水力参数设置，构建了准确的排水系统水力模型。水质模块针对多环芳烃设置了合理的污染物参数，并考虑不同土地利用类型对多环芳烃污染的影响。最后，通过实测数据对模型参数进行率定和验证，结果表明模型模拟结果与实测值具有较好的一致性，验证了模型的可靠性。该融雪径流模型的建立为后续深入研究哈尔滨市融雪径流中多环芳烃的污染空间分布及生态风险评价提供了有力的工具和数据支持。四、哈尔滨市融雪径流中多环芳烃污染空间分布4.1融雪径流量空间分布特征4.1.1不同汇水区空间分布规律哈尔滨市地形呈现东南高、西北低的态势，这种地形差异导致不同汇水区的融雪径流量存在显著不同。在地势较高的东南部山区汇水区，如阿什河上游汇水区，其坡度较大，积雪融化后的水流能够迅速汇集并沿地势流下，融雪径流量相对较大。相关研究表明，该汇水区的平均融雪径流量可达[X1]立方米/小时。而在地势较为平坦的松花江沿岸汇水区，水流速度相对较慢，且部分融雪水会通过下渗等方式流失，融雪径流量相对较小，平均融雪径流量约为[X2]立方米/小时。土地利用类型对融雪径流量也有重要影响。城市建成区汇水区由于不透水面积大，如南岗区部分汇水区，道路、建筑物等覆盖面积广，积雪融化后难以渗入地下，大部分形成地表径流，融雪径流量较大。据统计，该区域汇水区的融雪径流量峰值可达到[X3]立方米/小时。而在植被覆盖度高的汇水区，如呼兰河部分汇水区周边的林地和湿地，植被和土壤对融雪水具有较强的截留和下渗作用，使得融雪径流量相对较小，平均融雪径流量在[X4]立方米/小时左右。4.1.2不同排水口空间分布规律不同功能区域的排水口融雪径流量存在明显变化。在交通繁忙的主干道排水口，如中山路与红军街交口处的排水口，由于车流量大，道路表面积雪受车辆碾压和尾气排放影响，积雪融化速度加快，且道路表面相对光滑，融雪水易于汇集，导致融雪径流量较大。在融雪高峰期，该排水口的融雪径流量可达[X5]升/分钟。居民区排水口的融雪径流量则受到居民生活活动和小区布局的影响。老旧居民区由于排水系统相对不完善，且小区内道路狭窄，积雪融化后排水不畅，融雪径流量相对不稳定，在融雪初期和高峰期，径流量波动较大。而新建居民区通常规划有较好的排水设施，且绿化面积相对较大，融雪径流量相对较为平稳，如群力新区某居民区排水口，平均融雪径流量在[X6]升/分钟左右。商业区排水口，尤其是大型商场周边的排水口，由于人流量大，地面硬化程度高，且商业活动产生的热量可能加速积雪融化，融雪径流量也较大。例如，中央大街附近的排水口，在融雪期间平均融雪径流量可达到[X7]升/分钟。4.1.3不同下垫面空间分布规律屋面作为城市下垫面的重要组成部分，其融雪径流量具有独特特点。以居民楼屋面为例，屋面坡度和面积大小影响融雪径流量。坡度较大的屋面，积雪融化后水流速度快，融雪径流量相对较大。研究发现，坡度为[具体坡度1]的居民楼屋面，融雪径流量平均为[X8]升/平方米。而平屋面的融雪径流量相对较小，因为平屋面的积雪融化后可能会在屋面形成积水，部分水分会通过蒸发等方式散失。道路下垫面的融雪径流量与交通流量、路面材质等因素密切相关。城市主干道由于车流量大，路面温度相对较高，积雪融化快，且路面通常为沥青或水泥材质，不透水性强，融雪径流量较大。如学府路主干道，其融雪径流量在融雪高峰期可达[X9]升/米。而人行道由于人流量相对较小，路面材质多为砖石，透水性相对较好，融雪径流量相对较小，平均融雪径流量在[X10]升/米左右。绿地具有良好的截留和下渗能力，能够有效减少融雪径流量。公园绿地中的草坪和树林下垫面，其融雪径流量明显小于其他下垫面。研究表明，草坪下垫面的融雪径流量仅为[X11]升/平方米，树林下垫面由于植被的阻挡和土壤的过滤作用，融雪径流量更小，平均融雪径流量在[X12]升/平方米左右。绿地中的融雪水大部分被土壤吸收，补充地下水，只有少部分形成地表径流。4.2融雪径流中多环芳烃污染空间分布特征4.2.1不同环数空间分布规律在哈尔滨市融雪径流中，2-3环多环芳烃主要分布在城市中心和交通繁忙区域。这些区域由于汽车尾气排放、道路扬尘等污染源较多，使得低环数多环芳烃易于在空气中传播并沉降到积雪中，随着融雪径流进入水体。如道里区中央大街附近的融雪径流中，2-3环多环芳烃的浓度较高，平均浓度可达[X13]μg/L。4环多环芳烃在不同区域的分布相对较为均匀，但在工业集中区和部分老旧居民区的融雪径流中浓度相对较高。工业集中区如平房区的部分工业园区，工业生产过程中的废气排放是4环多环芳烃的重要来源。老旧居民区由于房屋建筑年代久远，取暖方式多为燃煤，燃煤过程中产生的多环芳烃也会增加融雪径流中4环多环芳烃的含量。这些区域融雪径流中4环多环芳烃的平均浓度约为[X14]μg/L。5-6环多环芳烃主要分布在远离城市中心的工业区和部分交通干道周边。这些高环数多环芳烃通常来源于高温燃烧过程，如工业锅炉燃烧、机动车发动机高温燃烧等。在远离城市中心的工业区，工业活动相对集中，高温燃烧源较多，使得5-6环多环芳烃在融雪径流中的浓度较高。例如，哈经开区某工业区周边融雪径流中5-6环多环芳烃的平均浓度可达[X15]μg/L。交通干道周边由于机动车尾气排放中含有高环数多环芳烃，也导致该区域融雪径流中5-6环多环芳烃浓度相对较高。4.2.2不同汇水区空间分布规律不同汇水区多环芳烃污染程度和组成存在显著差异。在松花江流域的主城区汇水区，由于人口密集、交通繁忙以及工业活动集中，多环芳烃污染程度较高。以道外区某汇水区为例，该区域融雪径流中多环芳烃的总浓度可达[X16]μg/L。其多环芳烃组成以2-4环为主，这主要是因为主城区的污染源主要包括汽车尾气排放、工业废气排放和居民生活燃煤等，这些污染源排放的多环芳烃以低环数和中等环数为主。而在阿什河上游汇水区，由于植被覆盖度高，工业活动较少，多环芳烃污染程度相对较低，融雪径流中多环芳烃总浓度平均为[X17]μg/L。该汇水区多环芳烃组成中，2-3环多环芳烃占比较大，这是因为该区域主要污染源为大气沉降，而大气中的多环芳烃以低环数为主，且植被对高环数多环芳烃具有一定的吸附和过滤作用。在呼兰河部分汇水区，由于农业活动和少量工业活动的影响，多环芳烃污染呈现出独特的特征。该区域融雪径流中多环芳烃总浓度在[X18]μg/L左右，多环芳烃组成中除了2-4环外，还含有一定比例的5-6环多环芳烃。这是因为农业活动中使用的农药、化肥等可能含有多环芳烃，工业活动中的高温燃烧过程也会产生高环数多环芳烃。4.2.3不同排水口空间分布规律不同排水口多环芳烃浓度和种类存在明显变化。在交通干道排水口，如长江路与南直路交口处的排水口，由于车流量大，汽车尾气排放和道路扬尘是多环芳烃的主要来源，多环芳烃浓度较高，总浓度可达[X19]μg/L。其多环芳烃种类丰富，包含2-6环多环芳烃，其中2-4环多环芳烃占比较大，这与汽车尾气排放中多环芳烃的组成特征相符。居民区排水口的多环芳烃浓度和种类受到居民生活方式和小区环境的影响。老旧居民区排水口，由于居民燃煤取暖等原因，多环芳烃浓度相对较高，如南岗区某老旧居民区排水口，融雪径流中多环芳烃总浓度可达[X20]μg/L。多环芳烃种类以2-4环为主，其中2-3环多环芳烃主要来源于燃煤排放，4环多环芳烃可能与居民生活中的其他燃烧活动有关。新建居民区排水口，由于小区绿化较好，且居民生活能源多为天然气等清洁能源，多环芳烃浓度相对较低，平均总浓度在[X21]μg/L左右。多环芳烃种类相对较少，以2-3环多环芳烃为主，主要来源于大气沉降。商业区排水口，由于商业活动集中，人流量大，多环芳烃浓度也较高。例如，果戈里大街附近商业区排水口，融雪径流中多环芳烃总浓度可达[X22]μg/L。多环芳烃种类包含2-5环，其中3-4环多环芳烃占比较大，这可能与商业区的餐饮、娱乐等活动产生的废气排放有关。4.2.4不同下垫面空间分布规律屋面多环芳烃污染主要来源于大气沉降和建筑材料挥发。居民楼屋面融雪径流中多环芳烃浓度与建筑年代和周边环境有关。老旧居民楼屋面由于建筑材料老化，可能会释放出多环芳烃，且周边环境中的污染源较多，导致融雪径流中多环芳烃浓度较高，平均浓度可达[X23]μg/L。新建居民楼屋面由于采用了环保建筑材料，且周边绿化较好，多环芳烃浓度相对较低，平均浓度在[X24]μg/L左右。屋面多环芳烃种类以2-3环为主，主要来源于大气中的污染物沉降。道路下垫面多环芳烃污染主要来自汽车尾气排放、燃油泄漏和轮胎磨损等。城市主干道由于交通流量大，多环芳烃污染严重，如友谊路主干道，融雪径流中多环芳烃总浓度可达[X25]μg/L。多环芳烃种类丰富，涵盖2-6环，其中2-4环多环芳烃占比较大，这是因为汽车尾气排放是道路下垫面多环芳烃的主要来源，尾气中含有多种环数的多环芳烃。人行道多环芳烃污染相对较轻，平均总浓度在[X26]μg/L左右。多环芳烃种类以2-3环为主，主要来源于道路扬尘和大气沉降。绿地多环芳烃污染相对较轻，主要来源于大气沉降和土壤中的残留污染物。公园绿地中的草坪下垫面融雪径流中多环芳烃浓度较低，平均浓度为[X27]μg/L。多环芳烃种类以2-3环为主，这是因为绿地植被对大气中的多环芳烃具有吸附和过滤作用，且土壤中的微生物对多环芳烃有一定的降解能力。树林下垫面由于植被茂密，对多环芳烃的截留和降解作用更强，融雪径流中多环芳烃浓度更低，平均浓度在[X28]μg/L左右。多环芳烃种类也以2-3环为主，且含量相对较少。4.3本章小结本章对哈尔滨市融雪径流量和多环芳烃污染的空间分布特征进行了详细分析。融雪径流量在不同汇水区、排水口和下垫面呈现出明显的空间分布规律，受地形、土地利用类型、排水系统和下垫面特性等多种因素影响。多环芳烃污染在不同环数、汇水区、排水口和下垫面也存在显著差异，其污染分布与污染源的分布和排放特征密切相关。这些空间分布特征的研究结果为后续开展融雪径流中多环芳烃污染的生态风险评价提供了重要的基础数据和依据，有助于深入了解多环芳烃在融雪径流中的迁移转化规律，为制定针对性的污染防控措施提供科学支撑。五、哈尔滨市融雪径流中多环芳烃污染生态风险评价5.1融雪径流中多环芳烃污染源解析5.1.1比值法比值法是一种常用的多环芳烃源解析方法，通过分析特定多环芳烃同分异构体的比值来判断其污染源类型。在本研究中，主要采用荧蒽（FLT）与芘（PYR）的比值（FLT/PYR）以及苯并[a]蒽（B[a]A）与屈（CHR）的比值（B[a]A/CHR）来进行源解析。当FLT/PYR比值小于1.0时，通常表明多环芳烃主要来源于石油源，因为石油类物质中荧蒽的含量相对较低；当FLT/PYR比值在1.0-1.2之间时，多环芳烃可能来源于石油源和燃烧源的混合；当FLT/PYR比值大于1.2时，则主要来源于燃烧源，如机动车尾气排放、工业燃烧和生物质燃烧等。研究结果显示，哈尔滨市融雪径流中大部分采样点的FLT/PYR比值大于1.2，平均值为1.45，这表明哈尔滨市融雪径流中的多环芳烃主要来源于燃烧源。例如，在道里区某交通干道采样点，FLT/PYR比值高达1.62，说明该区域多环芳烃受机动车尾气排放等燃烧源影响较大。对于B[a]A/CHR比值，当该比值小于0.2时，多环芳烃主要来源于石油源；当比值在0.2-0.35之间时，可能是石油源和燃烧源的混合；当比值大于0.35时，主要来源于燃烧源。本研究中，哈尔滨市融雪径流中B[a]A/CHR比值平均值为0.42，大部分采样点该比值大于0.35，进一步证实了多环芳烃主要来源于燃烧源。如南岗区某工业区采样点，B[a]A/CHR比值为0.51，表明该区域多环芳烃主要来自工业燃烧等燃烧源。5.1.2主成分分析法主成分分析法（PCA）是一种多元统计分析方法，能够将多个相关变量转化为少数几个互不相关的综合变量，即主成分，从而提取数据中的主要信息，确定多环芳烃的主要污染源。本研究运用SPSS软件对哈尔滨市融雪径流中多环芳烃的浓度数据进行主成分分析。首先对数据进行标准化处理，消除量纲的影响，然后计算相关系数矩阵，确定主成分的数量和贡献率。根据特征值大于1的原则，提取了3个主成分，这3个主成分的累计贡献率达到了85.6%，能够较好地解释多环芳烃浓度数据的变化。第一主成分的贡献率为45.3%，该主成分中荧蒽、芘、苯并[a]蒽、屈等多环芳烃的载荷较高，这些多环芳烃主要来源于机动车尾气排放和工业燃烧，因此第一主成分代表了机动车尾气和工业燃烧源。例如，在香坊区某交通繁忙区域和工业集中区的采样点，第一主成分的得分较高，说明这些区域受机动车尾气和工业燃烧源的影响较大。第二主成分贡献率为28.7%，其中萘、苊烯、苊等低环数多环芳烃的载荷较高，这些多环芳烃可能来源于燃煤和生物质燃烧，因此第二主成分代表了燃煤和生物质燃烧源。在一些老旧居民区采样点，第二主成分得分较高，这些区域居民取暖可能仍以燃煤为主，生物质燃烧也较为常见，导致该区域多环芳烃受燃煤和生物质燃烧源影响明显。第三主成分贡献率为11.6%，菲、蒽等多环芳烃在该主成分上载荷较高，这些多环芳烃与石油类物质的泄漏和挥发有关，所以第三主成分代表了石油源。在一些加油站和停车场附近的采样点，第三主成分得分相对较高，表明这些区域可能存在石油类物质的泄漏和挥发，对多环芳烃污染有一定贡献。5.2融雪径流中多环芳烃污染生态风险评价5.2.1危害熵值(HQ)生态风险评价危害熵值（HQ）生态风险评价法是通过计算多环芳烃的实测浓度与毒性阈值的比值来评估生态风险。本研究根据相关文献确定了16种多环芳烃对水生生物的毒性阈值，计算了哈尔滨市融雪径流中不同采样点和不同下垫面多环芳烃的危害熵值。计算结果表明，不同采样点多环芳烃的危害熵值存在明显差异。在交通干道采样点，如道外区某主干道采样点，由于车流量大，多环芳烃污染严重，苯并[a]芘的危害熵值高达1.25，超过了1，表明该区域苯并[a]芘对水生生物具有高风险。而在公园绿地等污染较轻的采样点，多环芳烃的危害熵值普遍较低，如香坊区某公园采样点，大多数多环芳烃的危害熵值小于0.1，处于低风险水平。不同下垫面多环芳烃的危害熵值也有所不同。道路下垫面由于受汽车尾气排放、燃油泄漏等影响，多环芳烃危害熵值相对较高。其中，苯并[a]芘、苯并[b]荧蒽等致癌性多环芳烃的危害熵值在部分路段超过0.1，处于中风险水平。屋面下垫面多环芳烃主要来源于大气沉降，危害熵值相对较低，大部分多环芳烃处于低风险水平。绿地对多环芳烃有一定的吸附和降解作用，其下垫面多环芳烃危害熵值最低，整体处于低风险水平。5.2.2BaP毒性当量(EBaP)生态风险评价BaP毒性当量（EBaP）生态风险评价法是将多环芳烃的浓度换算为苯并[a]芘的毒性当量浓度，以此来评价多环芳烃的综合生态风险。本研究根据各多环芳烃的毒性当量因子，计算了哈尔滨市融雪径流中多环芳烃的BaP毒性当量浓度。结果显示，哈尔滨市融雪径流中多环芳烃的BaP毒性当量浓度范围为0.5-8.6μg/L。其中，在工业集中区和交通繁忙区域，BaP毒性当量浓度相对较高。如平房区某工业区采样点，BaP毒性当量浓度达到5.3μg/L，处于中等风险水平。而在一些远离工业区和交通干道的区域，如呼兰区某农村采样点，BaP毒性当量浓度仅为0.8μg/L，生态风险较低。将计算得到的BaP毒性当量浓度与相关标准进行对比，我国地表水环境质量标准中规定苯并[a]芘的限值为2.8μg/L。哈尔滨市部分采样点的BaP毒性当量浓度超过了该标准，表明这些区域融雪径流中多环芳烃对生态环境存在一定的潜在风险，需要引起重视。5.2.3风险熵值法(RQ)生态风险评价风险熵值法（RQ）是基于多环芳烃的总等效质量浓度和毒性当量因子来评估生态风险。首先计算多环芳烃的总等效质量浓度，然后根据风