蓟运河河岸带土壤重金属污染特征、来源解析及环境健康风险评价

蓟运河河岸带土壤重金属污染特征、来源解析及环境健康风险评价一、引言1.1研究背景与意义土壤作为人类赖以生存的重要自然资源，是生态系统的关键组成部分，其质量状况直接关系到生态平衡、食品安全以及人类健康。然而，随着全球工业化、城市化进程的加速推进，大量含有重金属的工业废水、废气和废渣未经有效处理便被排放到环境中，加之农业生产中农药、化肥以及污水灌溉的不合理使用，土壤重金属污染问题日益严峻，已成为全球关注的环境焦点之一。重金属是指密度大于4.5g/cm³的金属元素，如汞（Hg）、镉（Cd）、铅（Pb）、铬（Cr）、砷（As）、锌（Zn）、铜（Cu）、镍（Ni）等。这些重金属在土壤中具有高度的稳定性，难以被自然降解，且容易在土壤中不断累积。一旦土壤遭受重金属污染，不仅会导致土壤物理化学性质恶化，土壤肥力下降，影响植物的正常生长和发育，降低农作物的产量与质量，还可能通过食物链的富集作用进入人体，对人体健康造成潜在的严重威胁。例如，长期食用被镉污染的大米，可能引发“痛痛病”，导致人体骨骼疼痛、骨质疏松甚至骨折；铅污染则可能损害人体神经系统，影响儿童的智力发育；汞污染会对人体的肝脏、肾脏和神经系统产生毒性作用，严重时可导致中毒死亡。河岸带作为陆地生态系统与水生生态系统之间的过渡地带，具有独特的生态功能。它不仅能够调节河流水文过程、净化水质、维护生物多样性，还在物质循环和能量流动中发挥着关键作用。然而，河岸带也极易受到人类活动和河流污染的双重影响，成为重金属等污染物的重要汇聚区域。蓟运河作为海河流域的重要支流，是天津市重要的水源地之一，其河岸带生态环境对于保障区域生态安全和水资源可持续利用至关重要。近年来，随着蓟运河流域内工业生产、农业活动以及城市化进程的不断发展，蓟运河河岸带土壤面临着越来越严重的重金属污染风险。工业排放的废水、废渣中含有的重金属，通过地表径流、大气沉降等途径进入河岸带土壤；农业生产中使用的农药、化肥以及污水灌溉，也可能导致土壤中重金属含量超标。这些重金属在土壤中的累积，不仅会破坏河岸带的生态平衡，影响植被的生长和物种多样性，还可能通过食物链对人类健康构成潜在威胁。因此，开展蓟运河河岸带土壤重金属污染与环境健康风险影响评价研究具有极其重要的现实意义。通过对蓟运河河岸带土壤重金属污染状况的全面调查和分析，能够准确掌握重金属的含量、分布特征以及污染来源，为制定针对性的污染防治措施提供科学依据。同时，对土壤重金属污染可能带来的环境健康风险进行评估，有助于提前预警潜在的健康危害，保障当地居民的身体健康和生态环境安全，促进区域的可持续发展。1.2国内外研究现状土壤重金属污染问题长期以来一直是环境科学领域的研究重点，国内外众多学者围绕这一课题展开了广泛而深入的研究，取得了丰硕的成果。国外对土壤重金属污染的研究起步较早，在基础理论和分析方法等方面积累了深厚的经验。早在20世纪60年代，随着工业化进程的加速，土壤重金属污染问题逐渐凸显，欧美等发达国家率先开展了相关研究。早期的研究主要集中在重金属在土壤中的含量测定和分布特征分析，通过对不同地区土壤样品的采集和分析，初步掌握了重金属的污染状况。随着研究的深入，学者们开始关注重金属在土壤中的化学形态、迁移转化规律以及生物有效性等方面。例如，Tessier等[1]提出的连续提取法，将土壤中的重金属分为可交换态、碳酸盐结合态、铁锰氧化物结合态、有机结合态和残渣态，为研究重金属的化学形态提供了重要的方法，使得人们能够更深入地了解重金属在土壤中的赋存状态和潜在风险。在重金属污染来源解析方面，多元统计分析、同位素示踪等技术被广泛应用，通过对土壤中重金属含量与环境因素的相关性分析以及同位素组成的测定，准确识别出重金属的污染源，为污染治理提供了科学依据。在土壤重金属污染的风险评价方面，国外也建立了较为完善的体系。风险评价模型不断发展和创新，从最初简单的单因子评价法，逐渐发展到综合考虑多种因素的多因子评价法，如地累积指数法、潜在生态危害指数法、健康风险评价模型等。这些模型能够更全面、准确地评估土壤重金属污染对生态环境和人体健康的潜在风险。美国环保局（EPA）提出的健康风险评价模型，考虑了重金属通过不同途径进入人体的暴露剂量以及人体对重金属的吸收、代谢等因素，为评估土壤重金属污染对人体健康的风险提供了科学的方法和标准。国内对土壤重金属污染的研究始于20世纪80年代，虽然起步相对较晚，但发展迅速。在土壤重金属污染调查方面，国家先后开展了一系列大规模的土壤污染状况调查，如全国土壤污染状况调查、多目标地球化学调查等，基本摸清了我国土壤重金属污染的总体状况和区域分布特征。研究表明，我国部分地区土壤重金属污染较为严重，尤其是在一些工矿区、城市周边和农业产区。在污染来源方面，国内研究发现，工业排放、矿业开采、农业活动以及交通污染等是导致土壤重金属污染的主要原因。例如，矿业开采过程中产生的废渣、废水未经有效处理直接排放，会导致周边土壤中重金属含量大幅升高；农业生产中过量使用含重金属的农药、化肥以及污水灌溉，也会造成土壤重金属污染。在污染治理和修复技术方面，国内学者也进行了大量的研究和实践。物理修复、化学修复、生物修复等多种修复技术不断涌现，其中生物修复技术由于具有环境友好、成本低等优点，受到了广泛关注。植物修复技术通过利用超富集植物对重金属的吸收、转运和积累能力，将土壤中的重金属去除或降低到可接受水平；微生物修复技术则利用微生物对重金属的吸附、转化等作用，降低重金属的生物有效性和毒性。此外，联合修复技术，如植物-微生物联合修复、化学-生物联合修复等，也成为研究的热点，通过多种修复技术的协同作用，提高修复效果和效率。针对蓟运河河岸带土壤重金属污染，已有一些研究开展。杨金翠等通过对蓟运河汉沽段河岸带土壤样品的采集和分析，研究了铜、锌、铅、镉、铬、镍、汞、砷等8种重金属元素的含量及分布特征，发现该区域土壤中部分重金属存在不同程度的超标现象，且锌与铬、镍、镉和铜呈显著正相关，初步判断这些元素可能具有同源性或存在复合污染。王斌等采用地积累指数法和潜在生态危害系数法对天津河流入海口沉积物重金属污染状况进行评价，结果表明蓟运河入海口沉积物重金属污染较重，处于强风险水平，其中Hg污染不容忽视，多数入海口沉积物存在不同程度的Hg和Cd污染生态风险。然而，目前关于蓟运河河岸带土壤重金属污染的研究仍存在一些不足之处。现有研究多集中在个别重金属元素或局部区域，缺乏对整个蓟运河河岸带土壤重金属污染的系统全面调查和分析；在污染来源解析方面，虽然已有一些研究运用相关性分析等方法进行初步判断，但对于复杂的污染来源，仍需进一步结合更先进的技术手段，如正定矩阵因子分解模型（PMF）、稳定同位素示踪技术等，进行准确识别；在环境健康风险评价方面，现有的研究多采用传统的评价模型，对一些新型风险因素，如重金属的复合污染效应、低剂量长期暴露风险等考虑不足，评价结果的准确性和可靠性有待进一步提高；此外，针对蓟运河河岸带土壤重金属污染的治理和修复技术研究相对较少，缺乏适合该区域特点的有效治理措施和技术方案。综上所述，国内外在土壤重金属污染研究方面已取得了众多成果，但蓟运河河岸带土壤重金属污染研究仍存在一定的局限性。进一步深入开展蓟运河河岸带土壤重金属污染的全面调查、准确解析污染来源、完善环境健康风险评价体系以及研发针对性的治理修复技术，具有重要的理论和现实意义。1.3研究内容与方法1.3.1研究内容蓟运河河岸带土壤样品采集与重金属含量分析：沿着蓟运河河岸带，综合考虑河流的不同河段、土地利用类型以及人类活动强度等因素，采用科学合理的布点方法，如系统布点法与随机布点法相结合，设置多个采样点。在每个采样点，按照一定的深度间隔，采集表层土壤（0-20cm）和深层土壤（20-60cm）样品。将采集到的土壤样品进行自然风干、去除杂质、研磨过筛等预处理后，运用先进的分析测试技术，如电感耦合等离子体质谱仪（ICP-MS）、原子吸收光谱仪（AAS）等，精确测定土壤中汞（Hg）、镉（Cd）、铅（Pb）、铬（Cr）、砷（As）、锌（Zn）、铜（Cu）、镍（Ni）等重金属元素的含量。土壤重金属污染评价：运用多种污染评价方法，如单因子污染指数法、内梅罗综合污染指数法、地累积指数法、潜在生态危害指数法等，对蓟运河河岸带土壤重金属污染程度和潜在生态风险进行全面、系统的评价。单因子污染指数法可直观反映单一重金属元素的污染状况；内梅罗综合污染指数法能综合考虑多种重金属元素的污染程度；地累积指数法可判断重金属元素的污染程度与自然背景值的关系；潜在生态危害指数法可评估重金属污染对生态环境的潜在危害程度。通过这些方法的综合运用，准确揭示蓟运河河岸带土壤重金属污染的空间分布特征和潜在生态风险等级。土壤重金属来源解析：运用多元统计分析方法，如主成分分析（PCA）、因子分析（FA）等，对土壤重金属含量数据进行处理，初步识别重金属的主要来源。通过分析不同重金属元素之间的相关性，提取主成分和因子，确定可能的污染源类型。结合正定矩阵因子分解模型（PMF），进一步解析土壤重金属的具体来源，量化各污染源对土壤重金属污染的贡献率。利用铅、硫等稳定同位素示踪技术，对重金属污染源进行更精确的溯源，通过测定土壤中重金属的同位素组成，并与潜在污染源的同位素特征进行对比，确定其来源，为制定针对性的污染防治措施提供科学依据。土壤重金属环境健康风险评价：基于美国环保局（EPA）推荐的健康风险评价模型，结合蓟运河河岸带当地的实际情况，确定相关参数，如土壤重金属含量、土壤摄取率、暴露时间、暴露频率、人体体重等，对土壤重金属通过不同途径（如经口摄入、皮肤接触、呼吸吸入）进入人体所导致的健康风险进行定量评估。分别计算重金属的非致癌风险和致癌风险，通过比较风险值与相应的风险阈值，判断土壤重金属污染对当地居民健康的潜在威胁程度。同时，考虑重金属的复合污染效应，采用联合毒性模型，评估多种重金属共同作用下的健康风险，使评价结果更符合实际情况。基于评价结果的污染防治建议：根据蓟运河河岸带土壤重金属污染评价、来源解析和环境健康风险评价的结果，针对性地提出切实可行的污染防治建议和措施。对于不同污染程度和来源的区域，制定差异化的治理策略。在源头控制方面，加强对工业企业的监管，严格控制重金属污染物的排放；优化农业生产方式，减少农药、化肥的不合理使用，推广绿色农业技术。在污染治理方面，针对轻度污染区域，采用植物修复、微生物修复等原位修复技术；对于中度和重度污染区域，可考虑采用化学淋洗、固化稳定化等异位修复技术或联合修复技术。同时，加强河岸带生态系统的保护和恢复，提高土壤的自净能力，降低重金属污染对生态环境和人体健康的风险。1.3.2研究方法野外采样与室内分析方法：野外采样过程严格遵循相关标准和规范，确保样品的代表性和可靠性。在采样前，对采样区域进行详细的勘察和记录，包括地理位置、土地利用类型、周边环境等信息。使用GPS定位仪准确记录采样点的坐标。采集的土壤样品装入洁净的样品袋中，及时运回实验室进行处理。室内分析过程中，对土壤样品的预处理、消解和测试等环节进行严格质量控制，采用国家标准方法或国际公认的分析方法进行操作。每批样品分析时，均设置空白样品、平行样品和标准参考物质，以确保分析结果的准确性和精密度。污染评价方法：单因子污染指数法计算公式为P_{i}=C_{i}/S_{i}，其中P_{i}为第i种重金属的单因子污染指数，C_{i}为第i种重金属的实测含量，S_{i}为第i种重金属的评价标准。当P_{i}\leq1时，表明土壤未受到该重金属污染；当P_{i}>1时，表明土壤受到该重金属污染，且P_{i}值越大，污染程度越严重。内梅罗综合污染指数法计算公式为P_{综}=\sqrt{\frac{(P_{imax})^2+(P_{i均})^2}{2}}，其中P_{综}为内梅罗综合污染指数，P_{imax}为单因子污染指数中的最大值，P_{i均}为单因子污染指数的平均值。该方法综合考虑了单因子污染指数的最大值和平均值，能更全面地反映土壤的综合污染状况。地累积指数法计算公式为I_{geo}=\log_{2}(\frac{C_{i}}{1.5B_{i}})，其中I_{geo}为地累积指数，C_{i}为第i种重金属的实测含量，B_{i}为第i种重金属的地球化学背景值，1.5为考虑到自然成岩作用可能引起背景值变动而取的系数。根据I_{geo}值的大小，将土壤重金属污染程度划分为7个等级，从无污染到极强污染。潜在生态危害指数法计算公式为RI=\sum_{i=1}^{n}E_{r}^{i}，其中RI为潜在生态危害指数，E_{r}^{i}为第i种重金属的潜在生态危害系数，E_{r}^{i}=T_{r}^{i}\times\frac{C_{i}}{S_{i}}，T_{r}^{i}为第i种重金属的毒性响应系数，反映重金属的毒性水平和生物可利用性。通过RI值的大小，将潜在生态危害程度划分为轻微、中等、强、很强、极强5个等级。来源解析方法：主成分分析（PCA）是一种将多个变量转化为少数几个综合变量（主成分）的多元统计分析方法。通过对土壤重金属含量数据进行PCA分析，可提取出主成分，并计算各主成分的贡献率和因子载荷。根据因子载荷的大小，确定各主成分所代表的重金属元素组合，从而初步判断重金属的来源。例如，若某主成分中铅、锌、镉等元素的因子载荷较大，可能表明该主成分代表的污染源与工业活动或交通污染有关。因子分析（FA）是在主成分分析的基础上，进一步对主成分进行旋转，使因子的含义更加明确。通过因子分析，可得到各因子的得分和因子得分系数矩阵，从而确定各因子对土壤重金属污染的影响程度。正定矩阵因子分解模型（PMF）是一种基于受体模型的源解析方法，它通过对监测数据的拟合，将观测到的污染物浓度分解为不同的源贡献和源成分谱。在PMF模型中，假设观测到的污染物浓度矩阵X可以表示为源成分谱矩阵G和源贡献矩阵F的乘积加上残差矩阵E，即X=GF+E。通过优化算法，求解出G和F，从而确定土壤重金属的污染源类型和各污染源的贡献率。稳定同位素示踪技术是利用稳定同位素的自然丰度或人为标记的同位素作为示踪剂，追踪物质在环境中的来源、迁移和转化过程。例如，铅的稳定同位素^{206}Pb、^{207}Pb、^{208}Pb具有不同的原子质量和相对丰度，通过测定土壤中铅的同位素组成，并与潜在污染源（如工业废气、废水、废渣，以及汽车尾气等）的铅同位素特征进行对比，可确定土壤中铅的来源。健康风险评价方法：经口摄入途径的暴露剂量计算公式为ADD_{ing}=\frac{C\timesIR\timesEF\timesED}{BW\timesAT}\times10^{-6}，其中ADD_{ing}为经口摄入途径的日均暴露剂量（mg/kg・d），C为土壤中重金属的含量（mg/kg），IR为土壤摄取率（mg/d），EF为暴露频率（d/a），ED为暴露时间（a），BW为人体体重（kg），AT为平均暴露时间（d）。皮肤接触途径的暴露剂量计算公式为ADD_{derm}=\frac{C\timesSA\timesAF\timesABS\timesEF\timesED}{BW\timesAT}\times10^{-6}，其中ADD_{derm}为皮肤接触途径的日均暴露剂量（mg/kg・d），SA为皮肤暴露面积（cm²），AF为皮肤表面土壤附着系数（mg/cm²），ABS为皮肤对重金属的吸收效率系数。呼吸吸入途径的暴露剂量计算公式为ADD_{inh}=\frac{C\timesInhR\timesEF\timesED}{PEF\timesBW\timesAT}，其中ADD_{inh}为呼吸吸入途径的日均暴露剂量（mg/kg・d），InhR为呼吸速率（m³/d），PEF为颗粒物排放因子（m³/kg）。非致癌风险采用危害商（HQ）来评价，计算公式为HQ=\frac{ADD}{RfD}，其中HQ为危害商，ADD为日均暴露剂量，RfD为参考剂量。当HQ<1时，表明非致癌风险较低；当HQ\geq1时，表明存在一定的非致癌风险。致癌风险采用致癌风险指数（CR）来评价，计算公式为CR=ADD\timesSF，其中CR为致癌风险指数，SF为致癌斜率因子。一般认为，当CR在10^{-6}-10^{-4}之间时，致癌风险处于可接受范围；当CR>10^{-4}时，致癌风险较高。对于重金属复合污染效应的评估，采用联合毒性模型，如浓度加和模型（CA）和独立作用模型（IA）。浓度加和模型假设混合物中各成分的毒性作用具有加和性，独立作用模型假设各成分的毒性作用相互独立。通过比较两种模型的计算结果，综合评估重金属复合污染对人体健康的风险。二、研究区域概况与研究方法2.1研究区域概况蓟运河是中国海河流域北系的主要河流之一，地理位置处于滦河以西，潮白河以东，其水系发源于燕山南麓，跨越河北省、北京市和天津市3个省级行政区，流域面积达10228平方千米。蓟运河干流全长157千米，河道宽度在300-500米之间，自北向南流经天津市蓟州区、宝坻区、宁河区、滨海新区，最终在滨海新区蓟运河防潮闸处注入渤海。蓟运河的主要支流众多，包括泃河、州河、还乡河、还乡新河、黎河、兰泉河、箭杆河、双成河等。其中，泃河作为上游西支，发源于河北省兴隆县大青山南侧将军关外青灰岭，全长180公里，流域面积2687平方公里，其上游建有海子水库；州河为上游东支，发源于迁西县的鹿耳岭和遵化县的秋水岭，全长112公里，上有黎河、沙河、果河、淋河等分支，于蓟州区东北始入天津蓟州区境，汇合后进入于桥水库库区。蓟运河的水文特征较为复杂。其上游山区属于华北的多雨区，河水径流量较大；而下游为平原，地势平坦，坡降锐减，河道宽浅，且入海口受海潮顶托，河水宣泄困难，历史上常常决堤泛滥为患。为了调节水量、防洪减灾，上游建有邱庄、海子和于桥3座大型水库，下游建有蓟运河防潮闸，并在防潮闸上游另辟新河道入海。在干旱少雨的年份，上游河流除汛期外基本无下泄流量进入蓟运河。研究区域属暖温带半湿润大陆型季风气候，四季分明。年均降水量在610-780毫米之间，降水主要集中在夏季，且多以暴雨形式出现。夏季气温较高，平均气温可达25℃-28℃；冬季较为寒冷，平均气温在-5℃-0℃之间。这种气候条件对蓟运河的水文过程和土壤重金属的迁移转化有着重要影响。夏季的强降水可能导致地表径流增大，加速土壤中重金属的冲刷和迁移；而冬季的低温则可能影响土壤微生物的活性，进而影响重金属的生物有效性。蓟运河周边的土地利用类型丰富多样。在山区，主要以林地为主，植被覆盖率较高，能够起到保持水土、涵养水源的作用；在平原地区，耕地是主要的土地利用类型，种植着小麦、玉米、水稻等农作物，农业生产活动频繁，农药、化肥的使用以及污水灌溉等可能导致土壤重金属污染；此外，还有部分区域为建设用地，分布着城镇、村庄以及工业企业，工业排放、交通污染等也是土壤重金属污染的重要来源。在河岸带附近，还存在一定面积的湿地，湿地生态系统具有独特的生态功能，能够对河流水质进行净化，但其也容易受到重金属污染的影响。2.2样品采集与分析为全面、准确地掌握蓟运河河岸带土壤重金属污染状况，本研究依据相关标准和规范，结合研究区域的实际特点，精心开展了样品采集与分析工作。在采样点布设方面，采用带状布点法沿着蓟运河河岸带进行设置。从蓟运河的上游至下游，根据河流的不同河段特点、土地利用类型以及人类活动强度的差异，每隔一定距离选取一个采样区域。在每个采样区域内，从水体至陆地方向0-10米的范围内，均匀设置多个采样点，确保涵盖了不同生态环境和受污染程度可能不同的区域。共设置了[X]个采样点，使用GPS定位仪精确记录每个采样点的经纬度坐标，保证采样点位置的准确性和可追溯性。土壤样品采集时，对于每个采样点，分别采集表层土壤（0-20cm）和深层土壤（20-60cm）样品。使用不锈钢土钻垂直插入土壤，按照规定深度采集土样，以确保采集的样品能够代表相应深度土层的真实情况。每个采样点的取土深度及采样量严格保持均匀一致，土样上层与下层的比例相同，以避免因采样不均导致的数据偏差。在每个采样点，按照“随机”“等量”和“多点混合”的原则，采集5-10个子样，将这些子样充分混合，组成一个混合样品，以增强样品的代表性。每个混合样品的重量约为1kg，若样品数量过多，采用四分法将多余的土壤弃去。采集后的土壤样品及时运回实验室进行处理。首先，将土壤样品放置在通风良好、阴凉干燥的地方自然风干，避免阳光直射和高温环境，防止土壤中重金属形态发生变化。风干过程中，定期翻动土壤样品，加速风干进程，并使其均匀干燥。待土壤样品完全风干后，用镊子仔细去除其中的石块、动植物残余物等杂质，确保样品的纯净度。随后，使用玛瑙棒将土壤样品研磨，使其通过2毫米尼龙筛，进一步混匀，以保证样品的均匀性。为了精确测定土壤中铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等重金属元素的含量，依据国家标准方法（如GB/T17138-1997《土壤质量铜、锌的测定火焰原子吸收分光光度法》、GB/T17141-1997《土壤质量铅、镉的测定石墨炉原子吸收分光光度法》等）对土壤样品进行前处理。对于铜、锌、铅、镉、铬、镍等重金属元素，采用硝酸-氢氟酸-高氯酸消解体系，将土壤样品中的有机物完全消解，使重金属元素以离子态存在于溶液中，以便后续的仪器分析。具体操作步骤为：称取适量经过研磨过筛的土壤样品于聚四氟乙烯坩埚中，加入一定量的硝酸、氢氟酸和高氯酸，在电热板上缓慢升温消解，直至样品完全溶解，溶液澄清透明。消解完成后，将溶液转移至容量瓶中，用超纯水定容至刻度线，摇匀备用。对于汞和砷元素，由于其具有易挥发的特性，采用王水水浴消解的方法。称取适量土壤样品于具塞比色管中，加入王水，在沸水浴中加热消解一定时间，使汞和砷元素充分溶解在溶液中。消解结束后，冷却至室温，用超纯水定容，过滤后取上清液进行分析。采用电感耦合等离子体质谱仪（ICP-MS）、原子吸收光谱仪（AAS）等先进的分析仪器对消解后的土壤样品溶液进行测定。在分析过程中，严格按照仪器操作规程进行操作，确保仪器的稳定性和准确性。每批样品分析时，均设置空白样品、平行样品和标准参考物质。空白样品用于检测实验过程中是否存在污染；平行样品用于检验分析结果的精密度；标准参考物质用于验证分析方法的准确性和可靠性。通过对空白样品、平行样品和标准参考物质的分析，有效保证了土壤重金属含量测定结果的准确性和可靠性。2.3数据处理与分析方法本研究运用了多种科学严谨的数据处理与分析方法，以确保对蓟运河河岸带土壤重金属污染状况及环境健康风险的评估准确可靠。在土壤重金属含量的统计分析方面，采用了一系列常用的统计指标，以全面了解重金属含量的分布特征。通过计算最大值和最小值，能够明确土壤中各重金属元素含量的取值范围，直观反映出含量的极值情况，有助于发现可能存在的异常高值或低值区域，为后续深入分析提供线索。平均值作为数据的集中趋势代表，能够反映出整体的含量水平；中值则不受极端值的影响，更能体现数据的中间位置情况，在数据存在异常值时，中值能提供更稳健的参考。标准偏差用于衡量数据的离散程度，其值越大，说明数据的分散性越大，即不同采样点间重金属含量的差异越大，这有助于评估土壤重金属污染的空间异质性。变异系数是标准偏差与平均值的比值，它消除了量纲的影响，能够更准确地比较不同重金属元素含量的变异程度，对于判断哪些重金属元素在土壤中的分布更不均匀具有重要意义。相关性分析是探究不同重金属元素之间关系的重要手段。本研究采用Pearson相关系数法进行相关性分析，该方法通过计算两个变量之间的线性相关程度，来判断它们之间是否存在显著的相关性。若相关系数的绝对值越接近1，则表明两个变量之间的线性相关性越强；若相关系数接近0，则表示两者之间几乎不存在线性相关关系。通过对不同重金属元素含量数据进行相关性分析，能够初步判断它们的同源性或异源性。例如，若某几种重金属元素之间呈现显著的正相关关系，可能暗示它们具有相同的污染源，或者在环境中存在共同的迁移转化机制，导致它们在土壤中的含量变化具有一致性；反之，若相关性不明显，则说明它们的来源或迁移过程可能相互独立。在重金属污染评价方面，选用了多种评价方法，从不同角度对土壤重金属污染程度进行评估。地积累指数法是一种广泛应用的评价方法，它通过将土壤中重金属的实测含量与地球化学背景值进行比较，并考虑自然成岩作用对背景值的影响，能够直观地判断土壤中重金属的污染程度与自然背景值的偏离情况。根据计算得到的地积累指数值，可将污染程度划分为7个等级，从无污染到极强污染，为污染程度的量化提供了明确的标准。污染负荷指数法则综合考虑了多种重金属元素的污染情况，通过计算单一金属最高污染系数、某点污染负荷指数以及某区域（流域）的污染负荷指数，能够直观地反映各个重金属对污染的贡献程度，以及重金属在时间和空间上的变化趋势。该方法不仅能评估整体的污染水平，还能分析不同重金属元素在污染过程中的相对重要性，有助于针对性地制定污染防治策略。环境健康风险评价对于保障公众健康至关重要。潜在生态风险指数法基于重金属的毒性响应系数和污染程度，综合评估重金属污染对生态环境的潜在危害程度。该方法考虑了重金属的生物毒性和在环境中的累积效应，通过计算潜在生态危害指数，将潜在生态危害程度划分为轻微、中等、强、很强、极强5个等级，为生态风险的评估提供了量化依据。健康风险评估模型则依据美国环保局（EPA）推荐的方法，结合蓟运河河岸带当地的实际情况，考虑土壤重金属通过经口摄入、皮肤接触、呼吸吸入等不同途径进入人体的暴露剂量，以及人体对重金属的吸收、代谢等因素，分别计算非致癌风险和致癌风险。通过将计算得到的风险值与相应的风险阈值进行比较，能够准确判断土壤重金属污染对当地居民健康的潜在威胁程度，为采取有效的健康防护措施提供科学依据。在评估过程中，还充分考虑了重金属的复合污染效应，采用联合毒性模型，如浓度加和模型（CA）和独立作用模型（IA），来综合评估多种重金属共同作用下的健康风险，使评价结果更符合实际情况。三、蓟运河河岸带土壤重金属污染现状分析3.1土壤重金属含量及分布特征对采集的蓟运河河岸带[X]个土壤样品进行分析，获取了铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等8种重金属元素的含量数据，具体统计结果如表1所示。重金属元素最大值（mg/kg）最小值（mg/kg）平均值（mg/kg）中值（mg/kg）标准偏差变异系数（%）Cu[Cu最大值][Cu最小值][Cu平均值][Cu中值][Cu标准偏差][Cu变异系数]Zn[Zn最大值][Zn最小值][Zn平均值][Zn中值][Zn标准偏差][Zn变异系数]Pb[Pb最大值][Pb最小值][Pb平均值][Pb中值][Pb标准偏差][Pb变异系数]Cd[Cd最大值][Cd最小值][Cd平均值][Cd中值][Cd标准偏差][Cd变异系数]Cr[Cr最大值][Cr最小值][Cr平均值][Cr中值][Cr标准偏差][Cr变异系数]Ni[Ni最大值][Ni最小值][Ni平均值][Ni中值][Ni标准偏差][Ni变异系数]Hg[Hg最大值][Hg最小值][Hg平均值][Hg中值][Hg标准偏差][Hg变异系数]As[As最大值][As最小值][As平均值][As中值][As标准偏差][As变异系数]从表1中可以看出，蓟运河河岸带土壤中各重金属元素含量存在一定差异。其中，Zn的含量相对较高，平均值达到[Zn平均值]mg/kg，最大值为[Zn最大值]mg/kg；而Cd的含量相对较低，平均值为[Cd平均值]mg/kg，最小值仅为[Cd最小值]mg/kg。各重金属元素含量的变异系数也有所不同，变异系数反映了数据的离散程度，变异系数越大，说明该元素在不同采样点的含量差异越大。例如，Hg的变异系数达到[Hg变异系数]%，表明Hg在蓟运河河岸带土壤中的分布极不均匀，不同采样点间的含量差异较大；而Cr的变异系数相对较小，为[Cr变异系数]%，说明Cr在土壤中的分布相对较为均匀。为更直观地展示不同采样点土壤重金属含量的空间分布情况，绘制了重金属含量的空间分布图（图1）。在图1中，横坐标表示采样点的位置，从蓟运河的上游至下游依次排列；纵坐标表示重金属的含量。通过该图可以清晰地看到，不同重金属元素在不同采样点的含量呈现出不同的变化趋势。在蓟运河的上游部分采样点，Cu、Zn等重金属元素含量相对较低，随着河流向下游流动，在一些靠近城镇和工业区域的采样点，这些重金属元素的含量明显升高，可能是由于城镇和工业活动的影响，如工业废水排放、固体废弃物堆放等，导致重金属在土壤中累积。而对于Hg元素，在个别采样点出现了含量极高的峰值，这些采样点可能受到了特殊污染源的影响，如附近存在汞矿开采或相关化工企业。此外，还可以观察到，在河流的入海口附近，部分重金属元素（如Cd、Pb等）的含量也相对较高。这可能是因为入海口地区受到海洋潮汐和河流携带污染物的共同作用，使得重金属在此处汇聚和沉积。同时，入海口周边的土地利用类型可能较为复杂，包括农业用地、养殖用地等，农业生产中的农药、化肥使用以及养殖废水排放等也可能对土壤重金属含量产生影响。通过对土壤重金属含量及分布特征的分析，初步揭示了蓟运河河岸带土壤重金属污染的一些基本情况。不同重金属元素的含量和分布存在差异，且受到多种因素的影响，为后续进一步深入分析土壤重金属污染程度、来源以及环境健康风险奠定了基础。3.2河岸带土壤重金属含量的相关性分析运用Pearson相关系数法对蓟运河河岸带土壤中8种重金属元素（铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As））含量进行相关性分析，结果如表2所示。重金属元素CuZnPbCdCrNiHgAsCu1[Zn与Cu的相关系数][Pb与Cu的相关系数][Cd与Cu的相关系数][Cr与Cu的相关系数][Ni与Cu的相关系数][Hg与Cu的相关系数][As与Cu的相关系数]Zn[Zn与Cu的相关系数]1[Pb与Zn的相关系数][Cd与Zn的相关系数][Cr与Zn的相关系数][Ni与Zn的相关系数][Hg与Zn的相关系数][As与Zn的相关系数]Pb[Pb与Cu的相关系数][Pb与Zn的相关系数]1[Cd与Pb的相关系数][Cr与Pb的相关系数][Ni与Pb的相关系数][Hg与Pb的相关系数][As与Pb的相关系数]Cd[Cd与Cu的相关系数][Cd与Zn的相关系数][Cd与Pb的相关系数]1[Cr与Cd的相关系数][Ni与Cd的相关系数][Hg与Cd的相关系数][As与Cd的相关系数]Cr[Cr与Cu的相关系数][Cr与Zn的相关系数][Cr与Pb的相关系数][Cr与Cd的相关系数]1[Ni与Cr的相关系数][Hg与Cr的相关系数][As与Cr的相关系数]Ni[Ni与Cu的相关系数][Ni与Zn的相关系数][Ni与Pb的相关系数][Ni与Cd的相关系数][Ni与Cr的相关系数]1[Hg与Ni的相关系数][As与Ni的相关系数]Hg[Hg与Cu的相关系数][Hg与Zn的相关系数][Hg与Pb的相关系数][Hg与Cd的相关系数][Hg与Cr的相关系数][Hg与Ni的相关系数]1[As与Hg的相关系数]As[As与Cu的相关系数][As与Zn的相关系数][As与Pb的相关系数][As与Cd的相关系数][As与Cr的相关系数][As与Ni的相关系数][As与Hg的相关系数]1从表2可以看出，Zn与Cr、Ni、Cd和Cu之间呈现出显著的正相关关系，相关系数分别为[Cr与Zn的相关系数]、[Ni与Zn的相关系数]、[Cd与Zn的相关系数]、[Cu与Zn的相关系数]，均通过了显著性水平为0.01的检验。这表明这些元素在蓟运河河岸带土壤中的含量变化具有较强的一致性，可能存在共同的污染源或相似的迁移转化机制。例如，工业生产中的金属冶炼、电镀等行业，在生产过程中会排放含有多种重金属的废水、废气和废渣，这些污染物进入土壤后，可能导致Zn、Cr、Ni、Cd和Cu等元素在土壤中同时累积。此外，农业活动中使用的含重金属的农药、化肥以及污水灌溉，也可能是这些元素共同污染的来源之一。Ni与Cu、Zn、Cr和Cd也具有显著的正相关性，相关系数依次为[Ni与Cu的相关系数]、[Ni与Zn的相关系数]、[Ni与Cr的相关系数]、[Ni与Cd的相关系数]，同样在0.01水平上显著。这进一步印证了这些元素之间存在紧密的联系，它们可能在环境中相互作用，共同参与了土壤的污染过程。在自然环境中，土壤的成土母质可能富含这些重金属元素，在长期的风化、淋溶等地质作用下，这些元素逐渐释放到土壤中，使得它们在土壤中的含量具有相关性。Cr与Zn、Ni和Cd之间的显著正相关关系也十分明显，相关系数分别为[Cr与Zn的相关系数]、[Cr与Ni的相关系数]、[Cr与Cd的相关系数]，在0.01水平上显著。这说明Cr与这些元素在土壤中的行为具有相似性，它们可能受到相同的环境因素影响，或者在土壤中存在共同的载体，使得它们能够同时迁移和积累。例如，土壤中的黏土矿物、有机质等对重金属具有吸附作用，某些黏土矿物或有机质可能对Cr、Zn、Ni和Cd具有较强的亲和力，从而使它们在土壤中共同富集。Cd和Zn、Cr、Ni之间存在非常显著的正相关关系，相关系数较高，分别为[Cd与Zn的相关系数]、[Cd与Cr的相关系数]、[Cd与Ni的相关系数]，在0.01水平上显著。这表明Cd与这些元素的同源性可能性较大，它们在土壤中的来源和迁移过程可能紧密相关。研究表明，交通污染也是土壤重金属污染的一个重要来源，汽车尾气排放、轮胎磨损以及道路扬尘等，都可能导致Cd、Zn、Cr、Ni等重金属在土壤中积累，使得它们之间呈现出显著的正相关关系。然而，其他重金属之间的相关性并不明显。例如，Pb与其他重金属元素之间的相关系数大多较小，且未通过显著性检验，说明Pb在蓟运河河岸带土壤中的来源和迁移过程可能与其他元素相互独立。这可能是因为Pb的污染源较为特殊，主要来源于工业生产中的铅冶炼、蓄电池制造等行业，以及交通污染中的含铅汽油的使用。随着无铅汽油的推广使用，交通源对土壤中Pb的贡献逐渐减少，但工业排放仍然是重要的污染源，其排放的Pb可能通过不同的途径进入土壤，与其他重金属元素的来源和迁移过程不同。Hg与大多数重金属元素之间的相关性也不显著，相关系数均较小。Hg在土壤中的行为较为复杂，它具有较强的挥发性和生物累积性，其污染源可能包括工业废气排放、汞矿开采以及农业中含汞农药的使用等。这些污染源的排放方式和迁移途径与其他重金属元素存在差异，导致Hg与其他元素在土壤中的含量变化没有明显的相关性。As与其他重金属元素之间同样未呈现出显著的相关性。As的来源广泛，包括自然来源（如岩石风化）和人为来源（如采矿、冶金、农药使用等）。不同来源的As在土壤中的化学形态和迁移转化规律不同，可能导致其与其他重金属元素在土壤中的分布和含量变化不一致。根据重金属含量的相关性分析可初步判断重金属的同源性及异源性。本研究中，极为显著正相关的元素（如Cr与Cd、Zn与Cr等）具有一定的同源性或存在复合污染。这些元素可能来自相同的污染源，或者在环境中通过共同的迁移转化过程在土壤中累积。而其他相关性不太显著的元素（如Pb、Hg、As）的来源与Zn、Cr、Ni、Cd、Cu等可能不同，它们各自受到不同污染源和环境因素的影响。这一结果为后续进一步运用多元统计分析、正定矩阵因子分解模型（PMF）以及稳定同位素示踪技术等方法，准确解析蓟运河河岸带土壤重金属的来源提供了重要的线索和依据。通过深入探究重金属的来源，能够为制定针对性的污染防治措施提供科学指导，有效减少土壤重金属污染，保护蓟运河河岸带的生态环境。四、蓟运河河岸带土壤重金属污染评价4.1地积累指数法评价地积累指数法（IndexofGeoaccumulation，简称I_{geo}）由德国学者Muller于1979年提出，是一种用于研究沉积物及其它物质中重金属污染程度的定量指标。该方法不仅考虑了人为污染因素，还充分考虑了自然成岩作用可能引起的背景值变动，弥补了其他评价方法的不足，在土壤重金属污染评价中得到了广泛应用。其计算公式为：I_{geo}=\log_{2}(\frac{C_{i}}{1.5B_{i}})式中，I_{geo}为地积累指数；C_{i}为第i种重金属的实测含量（mg/kg）；B_{i}为第i种重金属的地球化学背景值（mg/kg），本研究采用天津市土壤重金属元素背景值作为参考；1.5为考虑到自然成岩作用可能引起背景值变动而取的系数，用于表征沉积特征、岩石地质及其它相关联的影响。根据I_{geo}值的大小，将土壤重金属污染程度划分为7个等级，具体评价标准如表3所示：地积累指数I_{geo}范围污染等级污染程度描述I_{geo}\leq00级无污染0\ltI_{geo}\leq11级轻度-中等污染1\ltI_{geo}\leq22级中等污染2\ltI_{geo}\leq33级中等-强污染3\ltI_{geo}\leq44级强污染4\ltI_{geo}\leq55级强-极严重污染I_{geo}\gt56级极严重污染运用上述公式，对蓟运河河岸带土壤中铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等8种重金属元素的地积累指数进行计算，结果如表4所示：重金属元素地积累指数平均值I_{geo}最小值I_{geo}最大值污染等级Cu[Cu地积累指数平均值][Cu地积累指数最小值][Cu地积累指数最大值][Cu污染等级]Zn[Zn地积累指数平均值][Zn地积累指数最小值][Zn地积累指数最大值][Zn污染等级]Pb[Pb地积累指数平均值][Pb地积累指数最小值][Pb地积累指数最大值][Pb污染等级]Cd[Cd地积累指数平均值][Cd地积累指数最小值][Cd地积累指数最大值][Cd污染等级]Cr[Cr地积累指数平均值][Cr地积累指数最小值][Cr地积累指数最大值][Cr污染等级]Ni[Ni地积累指数平均值][Ni地积累指数最小值][Ni地积累指数最大值][Ni污染等级]Hg[Hg地积累指数平均值][Hg地积累指数最小值][Hg地积累指数最大值][Hg污染等级]As[As地积累指数平均值][As地积累指数最小值][As地积累指数最大值][As污染等级]从计算结果来看，蓟运河河岸带土壤中不同重金属元素的地积累指数存在明显差异，反映出其污染程度各不相同。其中，Cd的地积累指数平均值相对较高，达到[Cd地积累指数平均值]，处于[Cd污染等级]，表明Cd在部分采样点的污染程度较为严重，可能存在较强的人为污染源影响。从最大值来看，Cd的I_{geo}最大值为[Cd地积累指数最大值]，处于[Cd最大污染等级]，说明在个别采样点，Cd的污染已达到[Cd最大污染等级描述]水平，需重点关注这些区域的污染来源和治理措施。Hg的地积累指数平均值为[Hg地积累指数平均值]，处于[Hg污染等级]，且其最大值为[Hg地积累指数最大值]，处于[Hg最大污染等级]，说明Hg在蓟运河河岸带土壤中也存在一定程度的污染，在部分区域污染较为突出。Hg具有较强的挥发性和生物累积性，其污染来源可能包括工业废气排放、汞矿开采以及农业中含汞农药的使用等。相比之下，Cu、Zn、Pb、Cr、Ni、As等重金属元素的地积累指数平均值相对较低，大多处于无污染或轻度-中等污染水平。例如，Cu的地积累指数平均值为[Cu地积累指数平均值]，处于[Cu污染等级]，表明Cu在整体上污染程度较轻，但在个别采样点，其I_{geo}最大值为[Cu地积累指数最大值]，达到了[Cu最大污染等级]，说明在这些局部区域，Cu的污染情况也不容忽视。通过地积累指数法的评价，明确了蓟运河河岸带土壤中主要污染元素为Cd和Hg，它们在部分区域的污染程度已达到中等及以上水平。而其他重金属元素虽然整体污染程度相对较轻，但在个别采样点也存在一定的污染风险。这一评价结果为进一步分析土壤重金属污染来源和制定针对性的污染防治措施提供了重要依据。后续可结合其他评价方法和分析技术，深入探究污染原因，为保护蓟运河河岸带土壤生态环境提供科学支撑。4.2污染负荷指数法评价污染负荷指数（PollutionLoadIndex，PLI）由Tomlinson在1980年提出，该方法能够综合反映土壤中多种重金属的污染程度，同时还能直观地展现各个重金属对污染的贡献程度，以及重金属在时间和空间上的变化趋势。其计算过程如下：单一金属最高污染系数：单一金属最高污染系数C_{f}^{i}的计算公式为C_{f}^{i}=C_{s}^{i}/C_{n}^{i}。其中，C_{s}^{i}表示第i种重金属的实测含量（mg/kg）；C_{n}^{i}表示第i种重金属的参比值，本研究采用天津市土壤重金属元素背景值作为参比值。C_{f}^{i}值越大，表明该种重金属的污染程度越高。某点污染负荷指数：某点污染负荷指数PLI_{site}的计算公式为PLI_{site}=\sqrt[n]{\prod_{i=1}^{n}C_{f}^{i}}。这里，n代表参与评价的重金属种类数量。PLI_{site}综合考虑了多种重金属在某一点的污染情况，其值越大，说明该点土壤受到多种重金属的综合污染程度越严重。某区域（流域）的污染负荷指数：某区域（流域）的污染负荷指数PLI_{zone}的计算公式为PLI_{zone}=\sqrt[m]{\prod_{j=1}^{m}PLI_{sitej}}。其中，m表示某区域（流域）内采样点的数量。PLI_{zone}能够反映整个区域（流域）的土壤重金属污染综合水平，通过比较不同区域的PLI_{zone}值，可以评估不同区域土壤重金属污染程度的差异。根据上述公式，对蓟运河河岸带土壤中铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等8种重金属元素进行污染负荷指数计算，结果如表5所示：采样点编号C_{f}^{Cu}C_{f}^{Zn}C_{f}^{Pb}C_{f}^{Cd}C_{f}^{Cr}C_{f}^{Ni}C_{f}^{Hg}C_{f}^{As}PLI_{site}1[C_{f}^{Cu1}][C_{f}^{Zn1}][C_{f}^{Pb1}][C_{f}^{Cd1}][C_{f}^{Cr1}][C_{f}^{Ni1}][C_{f}^{Hg1}][C_{f}^{As1}][PLI_{site1}]2[C_{f}^{Cu2}][C_{f}^{Zn2}][C_{f}^{Pb2}][C_{f}^{Cd2}][C_{f}^{Cr2}][C_{f}^{Ni2}][C_{f}^{Hg2}][C_{f}^{As2}][PLI_{site2}]..............................[X][C_{f}^{CuX}][C_{f}^{ZnX}][C_{f}^{PbX}][C_{f}^{CdX}][C_{f}^{CrX}][C_{f}^{NiX}][C_{f}^{HgX}][C_{f}^{AsX}][PLI_{siteX}]经计算，蓟运河河岸带土壤的PLI_{zone}值为[PLI_{zone}]。一般认为，当PLI_{zone}\leq1时，表明研究区域土壤整体未受到明显的重金属污染；当PLI_{zone}>1时，则表明土壤受到了一定程度的重金属污染，且PLI_{zone}值越大，污染程度越严重。本研究中，[PLI_{zone}]大于1，说明蓟运河河岸带土壤整体存在一定程度的重金属污染。从不同区域来看，[具体区域1]的PLI_{site}平均值为[PLI_{site1}平均值]，[具体区域2]的PLI_{site}平均值为[PLI_{site2}平均值]。[具体区域1]的PLI_{site}平均值相对较高，表明该区域土壤受到重金属的综合污染程度相对较重，可能与该区域附近存在较多的工业企业、交通流量较大或农业活动中农药化肥使用量较大等因素有关。而[具体区域2]的PLI_{site}平均值相对较低，说明该区域土壤污染程度相对较轻，可能是由于该区域的生态环境较好，人类活动干扰较少。进一步分析各重金属的C_{f}^{i}值，发现Cd和Hg的C_{f}^{i}值在多个采样点均相对较高，对PLI_{site}的贡献较大。这与地积累指数法的评价结果一致，再次表明Cd和Hg是蓟运河河岸带土壤中的主要污染元素。其中，Cd的C_{f}^{i}最大值达到[C_{f}^{Cd}最大值]，说明在部分采样点，Cd的污染程度较为突出，可能存在特定的污染源，如附近存在镉矿开采或相关工业生产活动。Hg的C_{f}^{i}最大值为[C_{f}^{Hg}最大值]，也显示出在个别区域，Hg的污染较为严重，可能与汞矿开采、含汞工业废水排放或含汞农药使用等因素有关。通过污染负荷指数法的评价，全面了解了蓟运河河岸带土壤重金属的综合污染程度以及各重金属的污染贡献情况，明确了不同区域的污染负荷差异。这为后续制定针对性的污染防治措施提供了重要依据，对于污染程度较重的区域，应加强污染源排查和治理；对于主要污染元素Cd和Hg，需重点关注其来源和迁移转化规律，采取有效的措施降低其在土壤中的含量，以保护蓟运河河岸带的生态环境。4.3潜在生态风险指数法评价潜在生态风险指数法由瑞典科学家Hakanson于1980年提出，该方法从沉积学角度出发，综合考虑了土壤重金属含量、多元素协同作用、毒性水平、污染浓度以及环境对重金属污染敏感性等因素，能够全面、客观地评估土壤重金属污染对生态环境的潜在危害程度，在环境风险评价中得到了广泛应用。其计算公式如下：C_{f}^{i}=\frac{C_{s}^{i}}{C_{n}^{i}}E_{r}^{i}=T_{r}^{i}\timesC_{f}^{i}RI=\sum_{i=1}^{n}E_{r}^{i}=\sum_{i=1}^{n}T_{r}^{i}\timesC_{f}^{i}=\sum_{i=1}^{n}T_{r}^{i}\times\frac{C_{s}^{i}}{C_{n}^{i}}式中：RI为多元素环境风险综合指数，即潜在生态危害指数；E_{r}^{i}为第i种重金属的潜在生态危害系数；C_{f}^{i}为重金属i相对参比值的污染系数；C_{s}^{i}为重金属i的实测浓度（mg/kg）；C_{n}^{i}为重金属i的评价参比值，本研究采用天津市土壤重金属元素背景值作为评价参比值；T_{r}^{i}为重金属i的毒性响应系数，它主要反映重金属毒性水平和环境对重金属污染的敏感程度。常见重金属的毒性响应系数取值为：Hg=40，Cd=30，As=10，Pb=5，Cr=2，Cu=5，Ni=5，Zn=1。根据RI值的大小，将潜在生态危害程度划分为5个等级，具体评价标准如表6所示：潜在生态危害指数RI范围潜在生态危害程度等级RI\lt150轻微生态危害150\leqRI\lt300中等生态危害300\leqRI\lt600强生态危害600\leqRI\lt1200很强生态危害RI\geq1200极强生态危害运用上述公式，对蓟运河河岸带土壤中铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等8种重金属元素的潜在生态危害系数E_{r}^{i}和潜在生态危害指数RI进行计算，结果如表7所示：采样点编号E_{r}^{Cu}E_{r}^{Zn}E_{r}^{Pb}E_{r}^{Cd}E_{r}^{Cr}E_{r}^{Ni}E_{r}^{Hg}E_{r}^{As}RI1[E_{r}^{Cu1}][E_{r}^{Zn1}][E_{r}^{Pb1}][E_{r}^{Cd1}][E_{r}^{Cr1}][E_{r}^{Ni1}][E_{r}^{Hg1}][E_{r}^{As1}][RI1]2[E_{r}^{Cu2}][E_{r}^{Zn2}][E_{r}^{Pb2}][E_{r}^{Cd2}][E_{r}^{Cr2}][E_{r}^{Ni2}][E_{r}^{Hg2}][E_{r}^{As2}][RI2]..............................[X][E_{r}^{CuX}][E_{r}^{ZnX}][E_{r}^{PbX}][E_{r}^{CdX}][E_{r}^{CrX}][E_{r}^{NiX}][E_{r}^{HgX}][E_{r}^{AsX}][RIX]从计算结果可以看出，蓟运河河岸带土壤中不同重金属元素的潜在生态危害系数存在较大差异。其中，Hg和Cd的潜在生态危害系数相对较高，对潜在生态危害指数RI的贡献较大。在部分采样点，Hg的E_{r}^{Hg}值达到[E_{r}^{Hg}最大值]，Cd的E_{r}^{Cd}值达到[E_{r}^{Cd}最大值]，表明这两种重金属在这些区域具有较强的潜在生态危害。这主要是因为Hg和Cd具有较高的毒性响应系数，且在土壤中的含量超过背景值的倍数相对较大。例如，在一些工业活动频繁的区域，汞矿开采、含汞化工产品生产以及镉的冶炼等活动，会导致大量的Hg和Cd排放到环境中，通过大气沉降、地表径流等途径进入蓟运河河岸带土壤，从而增加了土壤中这两种重金属的含量，加大了潜在生态危害。而Zn、Cu、Pb、Cr、Ni、As等重金属元素的潜在生态危害系数相对较低。其中，Zn的E_{r}^{Zn}平均值为[E_{r}^{Zn}平均值]，Cu的E_{r}^{Cu}平均值为[E_{r}^{Cu}平均值]，Pb的E_{r}^{Pb}平均值为[E_{r}^{Pb}平均值]，Cr的E_{r}^{Cr}平均值为[E_{r}^{Cr}平均值]，Ni的E_{r}^{Ni}平均值为[E_{r}^{Ni}平均值]，As的E_{r}^{As}平均值为[E_{r}^{As}平均值]。这些重金属虽然在土壤中也有一定的含量，但由于其毒性响应系数相对较低，且大部分采样点的含量超过背景值的倍数较小，因此对潜在生态危害指数的贡献相对较小。从潜在生态危害指数RI来看，蓟运河河岸带土壤的潜在生态危害程度整体呈现出一定的空间差异。在[具体区域1]，RI平均值为[RI平均值1]，处于[潜在生态危害程度等级1]，表明该区域土壤的潜在生态危害程度相对较高，可能与该区域附近存在较多的工业污染源、交通繁忙以及农业生产中农药化肥的不合理使用等因素有关。而在[具体区域2]，RI平均值为[RI平均值2]，处于[潜在生态危害程度等级2]，潜在生态危害程度相对较轻，可能是因为该区域生态环境较好，人类活动干扰较少。综合分析，蓟运河河岸带土壤中主要的潜在生态风险贡献元素为Hg和Cd，这与地积累指数法和污染负荷指数法的评价结果一致。这些重金属的污染来源主要包括工业排放、矿业开采以及农业活动等。工业生产中的废气、废水和废渣排放，矿业开采过程中产生的尾矿和废石，以及农业生产中使用的含重金属农药、化肥和污水灌溉等，都可能导致Hg和Cd在土壤中积累，增加潜在生态风险。通过潜在生态风险指数法的评价，全面了解了蓟运河河岸带土壤重金属的潜在生态风险状况，明确了主要的风险贡献元素和高风险区域。这为制定针对性的生态保护和污染防治措施提供了重要依据，对于潜在生态危害程度较高的区域，应加强对污染源的管控，减少重金属的排放；对于主要风险贡献元素，需采取有效的治理措施，降低其在土壤中的含量和潜在生态危害，以保护蓟运河河岸带的生态环境和生态系统的健康稳定。五、蓟运河河岸带土壤重金属污染来源解析5.1主成分分析主成分分析（PrincipalComponentAnalysis，PCA）是一种多元统计分析方法，能够将多个相关变量转换为少数几个互不相关的综合变量，即主成分。这些主成分能够最大程度地保留原始数据的信息，从而简化数据结构，揭示数据之间的内在关系。在土壤重金属污染来源解析中，主成分分析通过对不同重金属元素含量数据的分析，提取出主成分，并根据主成分与重金属元素之间的相关性，判断可能的污染来源。运用主成分分析方法对蓟运河河岸带土壤中铜（Cu）、锌（Zn）、铅（Pb）、镉（Cd）、铬（Cr）、镍（Ni）、汞（Hg）、砷（As）等8种重金属元素含量数据进行处理。首先对数据进行标准化处理，消除量纲和数量级的影响，以确保分析结果的准确性和可靠性。然后计算相关系数矩阵，通过特征值分解，提取主成分。根据主成分的特征值大于1且累计贡献率大于80%的原则，确定了3个主成分，其贡献率分别为[主成分1贡献率]、[主成分2贡献率]、[主成分3贡献率]，累计贡献率达到[累计贡献率]，这表明这3个主成分能够解释原始数据中大部分的信息。各主成分的因子载荷矩阵如表8所示：重金属元素主成分1主成分2主成分3Cu[Cu在主成分1的因子载荷][Cu在主成分2的因子载荷][Cu在主成分3的因子载荷]Zn[Zn在主成分1的因子载荷][Zn在主成分2的因子载荷][Zn在主成分3的因子载荷]Pb[Pb在主成分1的因子载荷][Pb在主成分2的因子载荷][Pb在主成分3的因子载荷]Cd[Cd在主成分1的因子载荷][Cd在主成分2的因子载荷][Cd在主成分3的因子载荷]Cr[Cr在主成分1的因子载荷][Cr在主成分2的因子载荷][Cr在主成分3的因子载荷]Ni[Ni在主成分1的因子载荷][Ni在主成分2的因子载荷][Ni在主成分3的因子载荷]Hg[Hg在主成分1的因子载荷][Hg在主成分2的因子载荷][Hg在主成分3的因子载荷]As[As在主成分1的因子载荷][As在主成分2的因子载荷][As在主成分3的因子载荷]主成分1的贡献率最高，为[主成分1贡献率]。在主成分1中，Zn、Cr、Ni、Cd和Cu具有较高的正因子载荷，分别为[Zn在主成分1的因子载荷]、[Cr在主成分1的因子载荷]、[Ni在主成分1的因子载荷]、[Cd在主成分1的因子载荷]、[Cu在主成分1的因子载荷]。结合相关性分析结果，这些元素之间呈现显著的正相关关系，表明它们可能具有共同的污染源。工业活动是土壤重金属污染的重要来源之一，金属冶炼、电镀、化工等行业在生产过程中会产生大量含有多种重金属的废水、废气和废渣。这些污染物未经有效处理直接排放到环境中，通过大气沉降、地表径流等途径进入蓟运河河岸带土壤，导致Zn、Cr、Ni、Cd和Cu等元素在土壤中积累。例如，金属冶炼过程中，矿石中的重金属在高温熔炼等工艺下被释放出来，随废气排放到大气中，最终沉降到土壤表面；电镀行业产生的含重金属废水若未经达标处理直接排放，会通过地表径流污染河岸带土壤。因此，主成分1可能主要代表工业污染源。主成分2的贡献率为[主成分2贡献率]，其中Pb和Hg具有较高的因子载荷，分别为[Pb在主成分2的因子载荷]、[Hg在主成分2的因子载荷]。Pb的污染源较为复杂，工业生产中的铅冶炼、蓄电池制造等行业是其主要来源之一。随着无铅汽油的推广使用，交通源对土壤中Pb的贡献逐渐减少，但在一些交通繁忙的区域，汽车尾气排放、轮胎磨损以及道路扬尘等仍可能导致Pb在土壤中积累。Hg主要来源于工业废气排放、汞矿开采以及农业中含汞农药的使用等。在蓟运河河岸带，部分区域可能存在小型工业企业，其生产过程中排放的含Hg废气以及含Pb废渣等，可能是导致该区域土壤中Pb和Hg含量升高的原因。此外，农业活动中使用含汞农药也可能对土壤Hg含量产生一定影响。因此，主成分2可能代表工业和农业混合污染源。主成分3的贡献率为[主成分3贡献率]，As在该主成分中具有较高的因子载荷，为[As在主成分3的因子载荷]。As的来源广泛，自然来源包括岩石风化，人为来源包括采矿、冶金、农药使用等。在蓟运河河岸带，土壤中As可能主要来源于自然源和农业源。该区域的土壤母质中可能含有一定量的As，在长期的风化作用下，As逐渐释放到土壤中。同时，农业生产中使用的含砷农药和化肥，也可能增加土壤中As的含量。虽然As与其他重金属元素之间的相关性不明显，但通过主成分分析可以看出，As在主成分3中具有独特的贡献，表明其来源与其他元素有所不同。因此，主成分3可能主要代表自然源和农业污染源。通过主成分分析，初步确定了蓟运河河岸带土壤重金属的主要污染来源。主成分1代表工业污染源，主成分2代表工业和农业混合污染源，主成分3代表自然源和农业污染源。这一结果为进一步深入研究土壤重金属污染来源提供了重要线索，后续可结合正定矩阵因子分解模型（PMF）、稳定同位素示踪技术等方法，更准确地解析各污染源对土壤重金属污染的贡献率，为制定针对性的污染防治措施提供科学依据。5.2相关性分析辅助来源判断在对蓟运河河岸带土壤重金属进行来源解析的过程中，相关性分析是一项关键的初步判断手段，它与主成分分析等方法相互配合，能更准确地识别重金属的污染来源。通过前文的相关性分析可知，土壤中部分重金属元素之间呈现出显著的正相关关系，这为判断它们的同源性提供了重要线索。其中，Zn与Cr、Ni、Cd和Cu之间的显著正相关关系十分突出。这种相关性表明，这些元素在蓟运河河岸带土壤中的积累过程可能紧密相连，极有可能来源于相同的污染源。工业活动作为土壤重金属污染的重要源头之一，在蓟运河周边地区，金属冶炼、电镀、化工等行业分布较为广泛。以金属冶炼厂为例，在其生产过程中，矿石中的多种重金属元素如Zn、Cr、Ni、Cd和Cu会在高温熔炼、精炼等工艺环节中被释放出来。这些重金属一部分会随着工业废气排放到大气中，随后通过大气沉降的方式进入河岸带土壤；另一部分则可能随着工业废水排放到周边水体，再经由地表径流或灌溉用水进入土壤。电镀行业在生产过程中会使用大量含重金属的电镀液，若废水处理不当，直接排放到环境中，也会导致土壤中这些重金属元素的含量显著增加。因此，基于它们之间的强相关性，可以合理推测，这些元素可能主要源于工业生产活动。Ni与Cu、Zn、Cr和Cd之间同样存在显著的正相关性，这进一步支持了它们具有共同污染源的推断。在自然环境中，虽然土壤的成土母质可能含有这些重金属元素，但仅靠自然风化等地质作用，难以解释它们在土壤中如此紧密的相关性。而人为活动，尤其是工业活动，能够将这些重金属大量释放到环境中，并通过多种途径使其在土壤中共同积累。例如，一些小型的金属加工厂，由于生产设备简陋、环保措施不到位，在生产过程中会同时排放出含有Ni、Cu、Zn、Cr和Cd的污染物，这些污染物进入土壤后，使得这些元素在土壤中的含量呈现出同步变化的趋势。Cr与Zn、Ni和Cd之间的显著正相关关系，也从侧面反映了它们在环境行为上的相似性。土壤中的黏土矿物和有机质对重金属具有吸附作用，某些特定的黏土矿物或有机质可能对Cr、Zn、Ni和Cd具有较强的亲和力，从而使它们在土壤中共同迁移和富集。这也进一步说明，这些元素可能受到相同的环境因素影响，或者在土壤中存在共同的载体，使得它们能够同时迁移和积累。而工业排放的污染物中往往含有多种重金属，这些污染物进入土壤后，通过与土壤中的黏土矿物、有机质等相互作用，导致Cr、Zn、Ni和Cd在土壤中呈现出显著的正相关关系。Cd和Zn、Cr、Ni之间非常显著的正相关关系，更是表明它们的同源性可能性极大。除了工业活动，交通污染也是土壤重金属污染的一个重要来源。在蓟运河河岸带周边，交通道路密集，汽车尾气排放、轮胎磨损以及道路扬尘等，都会导致Cd、Zn、Cr、Ni等重金属在土壤中积累。汽车尾气中含有多种重金属，如轮胎磨损过程中会释放出Zn，刹车片磨损会产生Cd等。这些重金属随着交通活动的进行，不断在道路周边的土壤中积累，使得Cd与Zn、Cr、Ni之间呈现出显著的正相关关系。相比之下，Pb与其他重金属元素之间的相关性大多不明显，说明其来源具有独特性。Pb的污染源较为复杂，工业生产中的铅冶炼、蓄电池制造等行业是其主要来源之一。尽管随着无铅汽油的广泛使用，交通源对土壤中Pb的贡献有所减少，但在一些交通繁忙的区域，汽车尾气排放、轮胎磨损以及道路扬尘等仍可能导致Pb在土壤中积累。在蓟运河河岸带，部分区域可能存在小型的铅冶炼厂或蓄电池生产企业，这些企业排放的含Pb废气、废水和废渣，是导致该区域土壤中Pb含量升高的重要原因。由于其污染源与其他重金属元素的来源存在差异，所以在土壤中的含量变化与其他元素没有明显的相关性。Hg与大多数重金属元素之间的相关性也不显著，这与Hg的特殊性质和污染来源密切相关。Hg具有较强的挥发性和生物累积性，其污染源主要包括工业废气排放、汞矿开采以及农业中含汞农药的使用等。在蓟运河河岸带，若附近存在汞矿开采活动，开采过程中会释放大量的Hg，这些Hg会通过大气传播，沉降到土壤中。含汞化工产品生产企业排放的废气中也含有Hg，同样会对周边土壤造成污染。农业活动中使用含汞农药，虽然使用量相对较少，但长期积累也可能对土壤Hg含量产生一定影响。由于这些污染源的排放方式和迁移途径与其他重金属元素不同，导致Hg与其他元素在土壤中的含量变化没有明显的相关性。As与其他重金属元素之间未呈现出显著的相关性，这表明As的来源与其他元素有所不同。As的来源广泛，自然来源包括岩石风化，人为来源包括采矿、冶金、农药使用等。在蓟运河河岸带，土壤中As可能主要来源于自然源和农业源。该区域的土壤母质中可能含有一定量的As，在长期的风化作用下，As逐渐释放到土壤中。同时，农业生产中使用的含砷农药和化肥，也可能增加土壤中As的含量。由于As的来源相对独立，其在土壤中的化学形态和迁移转化规律与其他重金属元素不同，所以与其他元素在土壤中的分布和含量变化不一致。相关性分析结果与主成分分析相互印证，进一步明确了重金属的污染来源。在主成分分析中，主成分1中Zn、Cr、Ni、Cd和Cu具有较高的正因子载荷，与相关性分析中这些元素之间的显著正相关关系一致，表明主成分1代表的工业污染源是合理的推断。主成分2中Pb和Hg