赵县土壤重金属污染：基于多方法评价与源解析的研究

赵县土壤重金属污染：基于多方法评价与源解析的研究一、引言1.1研究背景与意义随着我国工业化、城镇化和农业现代化进程的加快，土壤重金属污染问题日益凸显。土壤作为生态系统的重要组成部分，不仅是人类赖以生存的物质基础，还在维持生态平衡、保障食品安全等方面发挥着关键作用。然而，工业生产活动以及农业上化肥农药等的大量、超量使用，使得大量重金属积聚于土壤中，这严重影响到土壤生态环境安全，直接或间接威胁人类健康。目前全国有超过8％的耕地存在重金属污染现象，且这一比例仍然在不断变大。在诸多重金属中，人们普遍关注并加以研究的重金属元素是对生物有较大毒性的Cr、Cd、Hg、Pb、Ni、Zn、Cu以及类金属As。赵县位于河北省中南部，是传统的农业大县，也是全国粮食生产先进县和中国优质梨果生产基地重点县。近年来，赵县的生物制药、化工、纺织印染、造纸包装等产业发展迅速，在推动经济增长的同时，也给当地的生态环境带来了一定压力。农业生产中，不合理的化肥、农药使用以及污水灌溉等现象普遍存在，这些因素都可能导致土壤重金属污染，进而影响农产品质量和生态环境安全。对赵县土壤重金属污染进行评价及源解析具有重要的现实意义。准确评估赵县土壤重金属污染状况，能够为当地土壤环境保护和污染治理提供科学依据，有助于制定针对性的污染防治措施，有效控制污染的进一步扩散，保障土壤生态环境安全。解析土壤重金属的来源，可以明确污染的主要成因，为从源头上减少重金属排放、实现精准治理提供指导，促进农业可持续发展。本研究还能为赵县国土空间的优化、生态红线的划定等提供数据支持和理论参考，推动区域经济与环境的协调发展。1.2国内外研究现状土壤重金属污染问题一直是国内外环境科学领域的研究热点。国外在土壤重金属污染研究方面起步较早，在污染特征、迁移转化规律以及治理修复技术等方面取得了丰硕成果。在污染特征研究方面，国外学者通过大量的实地监测和数据分析，对不同区域、不同土地利用类型下土壤重金属的含量、分布及变化趋势进行了系统研究。例如，美国地质调查局（USGS）长期开展土壤重金属监测工作，积累了丰富的数据资源，为深入了解土壤重金属污染特征提供了有力支持。欧洲一些国家也对本国土壤重金属污染状况进行了全面调查，明确了不同地区土壤重金属的污染程度和空间分布差异。在重金属迁移转化规律研究方面，国外学者借助先进的分析技术和模型，深入探讨了重金属在土壤-植物系统中的迁移途径、转化机制以及影响因素。如利用同位素示踪技术研究重金属在土壤中的迁移过程，通过建立数学模型预测重金属的迁移转化趋势。在治理修复技术方面，国外已研发出多种成熟的方法。物理修复技术如电动修复法，利用电场作用使土壤中的重金属离子向电极方向迁移，从而达到去除重金属的目的；化学修复技术如化学淋洗法，通过向土壤中添加淋洗剂，将重金属从土壤中溶解出来，再进行分离和处理；生物修复技术则利用植物、微生物或其联合体的作用，降低土壤中重金属的含量或毒性。例如，利用超富集植物吸收土壤中的重金属，通过微生物的代谢活动改变重金属的形态，降低其生物有效性。国内对土壤重金属污染的研究始于20世纪70年代，近年来随着环境问题的日益突出，相关研究得到了快速发展。在污染评价方面，国内学者结合我国土壤环境特点，建立了多种适合国情的评价方法和标准体系。如单因子污染指数法、内梅罗综合污染指数法、潜在生态危害指数法等，这些方法从不同角度对土壤重金属污染程度进行评价，为污染治理提供了科学依据。在源解析方面，国内学者综合运用多元统计分析、同位素示踪、受体模型等技术手段，对土壤重金属的来源进行了深入研究。例如，通过主成分分析和聚类分析，将土壤重金属污染源分为自然源和人为源，并进一步解析人为源的具体类型，如工业排放、交通污染、农业活动等。针对赵县的土壤重金属污染研究，目前已有一些成果。赵红安等人采集赵县表层土壤样品，检测Cr、Cd、Hg、Pb、Ni、Zn、Cu和As含量，采用单因子污染指数法和内梅罗指数法对赵县土壤的重金属污染状况进行评价，引用受体模型UNMIX模型对土壤中重金属来源进行解析。结果表明，Cu元素超标率为24.10%，污染程度最严重；赵县东北部谢庄乡和范庄镇是污染主要区域；土壤中重金属污染源包括“土壤源”“农药源”“交通-污灌源”，其中“交通-污灌源”是主要污染源。然而，现有研究仍存在一定不足。在污染评价方面，评价指标和方法相对单一，缺乏对土壤重金属污染的全面、综合评价；在源解析方面，虽然确定了主要污染源，但对各污染源的贡献率定量分析不够准确，且缺乏对污染源时空变化规律的研究。此外，针对赵县土壤重金属污染的治理修复研究较少，尚未形成一套完整的治理方案。因此，有必要进一步深入研究赵县土壤重金属污染状况，完善污染评价体系，精确解析污染源，为污染治理提供更科学、更有效的理论支持和技术指导。1.3研究内容与方法本研究旨在全面、深入地了解赵县土壤重金属污染状况，精确解析污染源，为赵县土壤环境保护和污染治理提供科学依据和技术支持。具体研究内容包括：土壤重金属污染评价：通过实地采样，获取赵县不同区域的土壤样品，运用电感耦合等离子体质谱（ICP-MS）等先进分析技术，准确测定土壤中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等重金属元素的含量。采用单因子污染指数法，对各重金属元素的污染程度进行单独评价，明确每种重金属在土壤中的污染状况。运用内梅罗综合污染指数法，综合考虑多种重金属元素的影响，全面评估赵县土壤的整体污染程度，确定污染严重的区域。土壤重金属源解析：运用多元统计分析中的主成分分析和聚类分析方法，对土壤重金属含量数据进行处理，初步识别土壤重金属的潜在污染源类型。基于正定矩阵因子分解（PMF）模型，结合赵县的工业布局、农业生产活动以及交通状况等实际情况，对土壤重金属的来源进行定量解析，确定各污染源对土壤重金属污染的贡献率。土壤重金属空间分布特征研究：借助地理信息系统（GIS）技术，将土壤重金属含量数据与地理空间信息相结合，绘制土壤重金属含量的空间分布图，直观展示赵县土壤中重金属的空间分布规律。运用地统计学方法，分析土壤重金属含量的空间变异特征，探讨影响土壤重金属空间分布的因素，如地形、土壤类型、土地利用方式等。本研究采用的主要方法如下：采样与分析方法：根据赵县的地形地貌、土地利用类型以及工业分布等因素，采用网格布点法和随机抽样法相结合的方式，在赵县境内共设置[X]个采样点，采集表层土壤样品（0-20cm）。每个采样点采集的土壤样品经自然风干、研磨、过筛后，采用盐酸-硝酸-氢氟酸-高氯酸消解体系进行消解，然后使用电感耦合等离子体质谱（ICP-MS）测定土壤中重金属元素的含量，并通过国家标准物质（GBW07405）进行质量控制，确保分析结果的准确性和可靠性。污染评价方法：单因子污染指数法是通过计算土壤中某一重金属元素的实测含量与该元素的评价标准值之比，来评价该元素的污染程度，公式为：P_i=\frac{C_i}{S_i}，其中P_i为第i种重金属的单因子污染指数，C_i为第i种重金属的实测含量，S_i为第i种重金属的评价标准值。内梅罗综合污染指数法综合考虑了单因子污染指数的最大值和平均值，能够更全面地反映土壤的污染状况，公式为：P_{综}=\sqrt{\frac{(P_{i\max}^2+\overline{P_i}^2)}{2}}，其中P_{综}为内梅罗综合污染指数，P_{i\max}为单因子污染指数的最大值，\overline{P_i}为单因子污染指数的平均值。源解析方法：主成分分析（PCA）是一种多元统计分析方法，通过对原始数据进行降维处理，将多个相关变量转化为少数几个不相关的主成分，从而提取数据的主要特征。聚类分析则是根据数据的相似性，将样品或变量分为不同的类别，以便对数据进行分类和解释。正定矩阵因子分解（PMF）模型是一种受体模型，它通过对观测数据进行因子分解，将观测数据分解为不同的源贡献和源成分，从而实现对污染源的解析。在本研究中，利用PMF模型对土壤重金属含量数据进行分析，确定土壤重金属的主要污染源及其贡献率。空间分析方法：地理信息系统（GIS）技术具有强大的空间数据处理和分析能力，能够将土壤重金属含量数据与地理空间信息进行整合，实现数据的可视化和空间分析。在本研究中，利用GIS软件绘制土壤重金属含量的空间分布图、污染等级图等，直观展示土壤重金属的空间分布特征。地统计学方法则是研究空间数据变异特征的有效工具，通过半方差函数分析土壤重金属含量的空间自相关性和变异性，确定其空间分布的结构特征和随机特征，为土壤重金属污染的空间预测和制图提供理论依据。二、研究区域概况与研究方法2.1赵县自然与社会经济概况赵县位于河北省中南部，地处东经114°36′-115°02′，北纬37°37′-37°53′之间，距省会石家庄主城区30公里，是石家庄都市圈的重要组成部分。县域总面积675平方公里，地势西北高东南低，地形开阔平坦，平均海拔41.8米，属于太行山东麓中段冲洪积形成的平原区，具体分为滹沱河冲洪积扇、洨河冲洪积扇、沙河冲洪积扇三大地貌单元。这种地形地貌特征使得赵县的土壤分布具有一定的规律性，为土壤重金属的迁移和积累提供了基础条件。赵县属东部季风气候区暖温带半湿润地区，四季分明，春季干燥多风，夏季炎热多雨，秋季天高气爽，冬季寒冷少雪。年平均气温13.9℃，年平均降水量498.9毫米，降水主要集中在夏季。光热充足，地下水丰富，有利于农业生产的发展。然而，这种气候条件也可能影响土壤中重金属的活性和迁移转化，例如降水可能导致重金属的淋溶和扩散，气温变化可能影响土壤微生物的活性，进而影响重金属的形态和生物有效性。赵县的土壤类型主要有潮土、褐土和风沙土。潮土是全县面积最大的土壤类型，分布广泛，主要集中在河流冲积平原地区，其土壤质地适中，保水保肥能力较强，是农业生产的主要土壤类型；褐土主要分布在县域的西北部和东南部，土壤肥力较高，适合多种农作物生长；风沙土主要分布在故河道及河滩地区，土壤质地较粗，保水保肥能力差，植被覆盖度较低。不同土壤类型的理化性质差异较大，对重金属的吸附、解吸和固定能力也不同，这对土壤重金属的污染状况和生态风险有着重要影响。赵县境内河流水系较为发达，主要河流有洨河、沙河、滹沱河故道等。这些河流不仅为农业灌溉和居民生活提供了水源，也在土壤重金属的迁移扩散过程中起到了重要作用。河流的水流搬运作用可能将上游地区的重金属污染物带到下游，导致下游地区土壤重金属含量升高；同时，河流的灌溉用水也可能将水中的重金属带入农田土壤，造成土壤污染。赵县是传统的农业大县，也是全国粮食生产先进县和中国优质梨果生产基地重点县。全县耕地面积58.45万亩，园地20.7万亩，主要农作物有小麦、玉米、棉花、花生等，特产雪花梨、鸭梨等。近年来，随着农业现代化进程的加快，赵县的农业生产方式发生了较大变化，化肥、农药的使用量不断增加，这可能导致土壤中重金属的积累。同时，农业灌溉用水的质量也可能影响土壤重金属含量，不合理的灌溉方式可能加剧土壤重金属污染。在工业方面，赵县初步形成了以医药产业为主导，现代食品加工、纺织服装、农机装备制造等特色产业协同发展的经济体系。现有规上企业115家、科技型中小企业662家、高新技术企业68家。赵县经济开发区是省政府批准设立的首批省级工业聚集区，实行“一区两园”模式，西区以医药制造业为主导产业，巩固提升现代食品、纺织服装、装备制造等传统产业；东区以医药制造为主，辅之发展相关的精细化工产业。工业生产过程中产生的废气、废水和废渣可能含有大量重金属，若未经有效处理直接排放，会对周边土壤环境造成严重污染。例如，化工企业排放的废水中可能含有汞、镉、铅等重金属，这些重金属随废水进入土壤后，会在土壤中积累，对土壤生态系统和农作物生长产生危害。赵县的服务业主要由文化旅游业为主导，以及住宿餐饮业、电子商务等来带动。赵县历史悠久，文物众多，有赵州桥、永通桥、柏林寺塔、陀罗尼经幢等国保文物6处、省保文物6处，每年吸引大量游客前来观光旅游。文化旅游业的发展虽然对土壤重金属污染的直接影响较小，但游客的大量涌入可能会带来一些间接影响，如旅游设施建设可能改变土地利用方式，从而影响土壤重金属的分布；游客的活动可能导致土壤扰动，增加重金属的迁移扩散风险。2.2土壤样品采集与分析在本研究中，为全面、准确地了解赵县土壤重金属污染状况，依据赵县的地形地貌、土地利用类型以及工业分布等因素，采用网格布点法和随机抽样法相结合的方式确定采样点。网格布点法能保证采样的均匀性，全面覆盖研究区域；随机抽样法则增加了采样的随机性，减少人为因素的干扰，使采样结果更具代表性。在赵县境内共设置了[X]个采样点，这些采样点广泛分布于赵县的各个乡镇，涵盖了农田、果园、林地以及靠近工业区域的土壤。采样过程严格按照相关标准和规范进行，确保采集的土壤样品能够真实反映当地的土壤状况。采样深度设定为0-20cm，这是因为表层土壤是与人类活动和生态系统交互最为频繁的部分，重金属在这一层的积累和分布情况对环境和人类健康的影响最为直接。每个采样点均按照梅花法进行采样，在10cm×10cm正方形的四个顶点和中心共五处各采集0-20cm耕层土壤1kg，将这五处采集的土壤充分混合，组成混合土样。这种采样方法能够综合考虑采样点周围土壤的情况，减少局部差异对结果的影响。采集的混合土样充分混合后，采用四分法反复取舍，最后保留1kg左右土样作为该点的混合样品。四分法是一种经典的样品缩减方法，通过将样品分成四等份，取对角的两份合并，再重复操作，直至得到所需的样品量，能有效保证样品的代表性。共采集了[X]个这样的混合样品，将其送往实验室进行分析。在实验室中，首先对土样进行自然风干处理，这一步骤在通风良好、无阳光直射的室内进行，目的是让土壤中的水分自然挥发，避免因烘干等高温处理导致土壤中重金属形态发生变化。风干后的土壤去除残根等杂物，这些杂物可能会影响土壤重金属含量的测定结果，必须予以清除。然后用木棍碾压，将较大的土块破碎，再用玛瑙研钵研细，过100目尼龙网筛，使土壤样品达到均匀、细腻的状态，便于后续的分析测试。为了准确测定土壤中重金属元素的含量，采用盐酸-硝酸-氢氟酸-高氯酸消解体系对土壤样品进行消解。该消解体系能够有效地破坏土壤的矿物结构，使其中的重金属元素释放出来，转化为可测定的离子态。具体操作过程为：称取一定量的过筛土壤样品于聚四氟乙烯坩埚中，加入适量的盐酸、硝酸、氢氟酸和高氯酸，在电热板上逐渐升温进行消解。消解过程中，密切观察溶液的颜色和状态变化，确保消解完全。消解完成后，将溶液冷却至室温，转移至容量瓶中，用超纯水定容至刻度线，待测。重金属含量的测定采用电感耦合等离子体质谱（ICP-MS）技术。ICP-MS具有灵敏度高、分析速度快、可同时测定多种元素等优点，能够准确测定土壤中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等重金属元素的含量。在测定过程中，为确保分析结果的准确性和可靠性，采取了严格的质量控制措施。每批样品均同时分析国家标准物质（GBW07405），国家标准物质的测定结果应在其标准值的不确定度范围内，以此来验证分析方法的准确性。同时，每10个样品插入一个空白样品，空白样品的测定结果应低于方法检出限，以监控分析过程中是否存在污染。对部分样品进行平行测定，平行样的相对偏差应小于10%，以保证测定结果的精密度。通过这些质量控制措施，有效确保了土壤重金属含量测定结果的可靠性，为后续的污染评价和源解析提供了准确的数据支持。2.3土壤重金属污染评价方法2.3.1单因子污染指数法单因子污染指数法是土壤重金属污染评价中最为基础和常用的方法之一，它通过将土壤中某一重金属元素的实测含量与该元素的评价标准值进行对比，从而对该元素的污染程度进行单独评价。这种方法的原理简单直观，能够清晰地反映出每种重金属在土壤中的污染状况，其计算公式为：P_i=\frac{C_i}{S_i}其中，P_i为第i种重金属的单因子污染指数，C_i为第i种重金属的实测含量（mg/kg），S_i为第i种重金属的评价标准值（mg/kg）。在实际应用中，评价标准值的选择至关重要，通常会根据研究目的、土壤类型以及相关的环境质量标准来确定。例如，在农业土壤评价中，常采用《土壤环境质量农用地土壤污染风险管控标准（试行）》（GB15618-2018）中的筛选值作为评价标准。当P_i\leq1时，表明土壤中该重金属元素的含量处于安全范围，未受到污染；当1<P_i\leq2时，说明土壤受到轻度污染；当2<P_i\leq3时，意味着土壤受到中度污染；当P_i>3时，则表示土壤受到重度污染。单因子污染指数法的优点是计算简便，能够快速准确地判断出单一重金属元素的污染程度，缺点是无法综合考虑多种重金属元素之间的相互作用和协同效应，不能全面反映土壤的整体污染状况。2.3.2内梅罗指数法内梅罗综合污染指数法综合考虑了单因子污染指数的最大值和平均值，能够更全面地反映土壤的污染状况。它克服了单因子污染指数法的局限性，从整体上对土壤的污染程度进行评价，其计算公式为：P_{综}=\sqrt{\frac{(P_{i\max}^2+\overline{P_i}^2)}{2}}其中，P_{综}为内梅罗综合污染指数，P_{i\max}为单因子污染指数的最大值，\overline{P_i}为单因子污染指数的平均值。内梅罗综合污染指数法通过对单因子污染指数的综合处理，既突出了污染最严重的重金属元素的影响，又考虑了其他重金属元素的平均污染水平，使得评价结果更加客观、全面。根据内梅罗综合污染指数的大小，可以将土壤污染程度划分为不同等级。当P_{综}\leq0.7时，土壤处于安全清洁状态；当0.7<P_{综}\leq1.0时，土壤处于警戒限，尚清洁；当1.0<P_{综}\leq2.0时，土壤受到轻度污染；当2.0<P_{综}\leq3.0时，土壤受到中度污染；当P_{综}>3.0时，土壤受到重度污染。内梅罗指数法在土壤重金属污染评价中应用广泛，能够为土壤污染的综合防治提供科学依据，但是该方法对高污染值较为敏感，可能会夸大某些污染程度较高区域的污染状况。2.3.3地积累指数法地积累指数法由德国科学家Müller提出，它不仅考虑了土壤中重金属元素的实测含量，还引入了背景值的概念，同时也考虑了自然成土过程中可能存在的背景值波动，能够更准确地反映土壤中重金属的污染程度以及人为活动对土壤重金属积累的影响。地积累指数法的计算公式为：I_{geo}=\log_2\frac{C_i}{1.5B_i}其中，I_{geo}为地积累指数，C_i为第i种重金属的实测含量（mg/kg），B_i为第i种重金属的地球化学背景值（mg/kg），1.5是考虑到自然成土过程中背景值波动而引入的修正系数。地积累指数法将土壤重金属污染程度分为7个等级，I_{geo}\leq0时，为无污染；0<I_{geo}\leq1时，为轻度污染；1<I_{geo}\leq2时，为偏中度污染；2<I_{geo}\leq3时，为中度污染；3<I_{geo}\leq4时，为偏重污染；4<I_{geo}\leq5时，为重度污染；I_{geo}>5时，为严重污染。地积累指数法在评价土壤重金属污染时，能够较好地反映出土壤中重金属的自然来源和人为来源的差异，对于研究土壤污染的历史演变和人为活动对土壤环境的影响具有重要意义。该方法依赖于准确的背景值数据，背景值的选择可能会对评价结果产生较大影响，在不同地区应用时需要根据实际情况合理确定背景值。2.3.4潜在生态风险指数法潜在生态风险指数法由瑞典科学家Hakanson提出，它综合考虑了重金属的毒性、含量以及环境对重金属污染的敏感性等因素，能够全面评估土壤中多种重金属对生态环境造成的潜在风险。潜在生态风险指数法的计算公式为：RI=\sum_{i=1}^{n}E_r^i=\sum_{i=1}^{n}T_r^i\times\frac{C_i}{C_n^i}其中，RI为潜在生态风险指数，E_r^i为第i种重金属的潜在生态风险系数，T_r^i为第i种重金属的毒性响应系数，C_i为第i种重金属的实测含量（mg/kg），C_n^i为第i种重金属的参比含量（mg/kg）。在计算过程中，不同重金属的毒性响应系数不同，Hg的毒性响应系数最高，为40；Cd的毒性响应系数为30；As的毒性响应系数为10；Pb、Cu、Ni的毒性响应系数均为5；Cr的毒性响应系数为2；Zn的毒性响应系数为1。参比含量通常选择当地土壤背景值或全球土壤背景值。根据潜在生态风险指数的大小，将潜在生态风险程度划分为5个等级，RI<150时，为低风险；150\leqRI<300时，为中等风险；300\leqRI<600时，为较高风险；600\leqRI<1200时，为高风险；RI\geq1200时，为很高风险。潜在生态风险指数法能够直观地反映出土壤中多种重金属对生态环境的潜在危害程度，为土壤污染的生态风险评估和治理提供了重要参考，该方法在确定毒性响应系数和参比含量时存在一定的主观性，不同的取值可能会导致评价结果有所差异。2.4土壤重金属源解析方法2.4.1多元统计分析多元统计分析是土壤重金属源解析中常用的方法之一，它能够对多个变量之间的复杂关系进行综合分析，从而揭示土壤重金属的潜在来源。主成分分析（PCA）和聚类分析（CA）是其中最具代表性的两种方法。主成分分析是一种通过降维技术将多个相关变量转化为少数几个不相关的主成分的多元统计方法。在土壤重金属源解析中，它基于土壤中各重金属元素含量之间的相关性，将众多重金属变量转换为少数几个综合指标，即主成分。每个主成分都是原始变量的线性组合，且各主成分之间相互独立。这些主成分能够最大限度地保留原始数据的信息，其中第一主成分通常包含了数据中最主要的变异信息，后续主成分依次包含较少的变异信息。通过对主成分的分析，可以提取出数据中的主要特征，进而识别出土壤重金属的潜在污染源类型。例如，如果某一主成分中某些重金属元素的载荷较高，说明这些重金属元素之间具有较强的相关性，可能来自同一污染源。聚类分析则是根据样品或变量之间的相似性，将其划分为不同类别的统计方法。在土壤重金属源解析中，它可以对土壤样品或重金属元素进行聚类。对土壤样品进行聚类时，相似的土壤样品会被归为一类，这些样品可能受到相同污染源的影响；对重金属元素进行聚类时，具有相似含量分布特征的重金属元素会被聚在一起，表明它们可能具有相同的来源。聚类分析能够直观地展示土壤重金属的分类情况，为污染源的识别提供重要线索。在实际应用中，主成分分析和聚类分析常常结合使用。首先通过主成分分析提取土壤重金属数据的主要特征，初步确定潜在的污染源类型；然后利用聚类分析对土壤样品或重金属元素进行分类，进一步验证和细化主成分分析的结果。这种结合使用的方法能够充分发挥两种方法的优势，提高土壤重金属源解析的准确性和可靠性。例如，在对某地区土壤重金属污染研究中，通过主成分分析发现土壤中重金属主要与工业排放、农业活动和自然源有关；再通过聚类分析，将土壤样品分为不同类别，分别对应不同的污染源，从而更准确地确定了各污染源对土壤重金属污染的贡献。2.4.2正定矩阵因子分解（PMF）模型正定矩阵因子分解（PMF）模型是一种广泛应用于土壤重金属源解析的受体模型，它基于矩阵分解的原理，能够有效解析土壤重金属的来源及各污染源的贡献率。PMF模型的基本原理是将观测到的土壤重金属浓度数据矩阵X分解为两个非负矩阵G和F，以及一个残差矩阵E，即X_{ij}=\sum_{k=1}^{p}G_{ik}F_{kj}+E_{ij}。其中，X_{ij}表示第i个样品中第j种重金属的浓度；G_{ik}表示第k个因子对第i个样品的贡献；F_{kj}表示第k个因子中第j种重金属的含量；E_{ij}表示残差。模型通过最小化目标函数Q=\sum_{i=1}^{n}\sum_{j=1}^{m}(\frac{E_{ij}}{u_{ij}})^2来确定G和F矩阵，其中u_{ij}是测量误差。在赵县土壤重金属源解析中应用PMF模型，首先需要对原始数据进行预处理。由于土壤中重金属含量数据可能存在异常值和缺失值，这些数据会影响模型的准确性和稳定性，因此需要对其进行处理。异常值可通过统计方法如3\sigma准则进行识别和剔除，缺失值则可采用均值插补、回归插补等方法进行填补。同时，为了消除不同重金属元素浓度量纲和数量级的差异，需要对数据进行标准化处理，使各元素的数据具有可比性。确定因子数是PMF模型应用中的关键步骤。因子数过多会导致模型过拟合，使结果过于复杂且难以解释；因子数过少则可能无法充分反映土壤重金属的来源信息。通常可采用平行运行PMF模型、观察目标函数Q值的变化、结合特征值分析以及参考相关研究经验等方法来确定合适的因子数。在赵县的研究中，经过多次试验和分析，最终确定了[X]个因子，分别代表不同的污染源。模型运行过程中，还需对结果进行质量控制和验证。通过检查残差分布是否符合正态分布、计算不确定性参数如G和F矩阵元素的不确定性范围等方式，来评估模型结果的可靠性。若残差分布不符合正态分布，可能表示模型存在缺陷，需要进一步调整；不确定性参数过大则说明结果的可靠性较低，需要重新分析数据或调整模型参数。只有经过严格质量控制和验证的结果，才能用于后续的污染源解析和分析。2.4.3UNMIX模型UNMIX模型是一种基于非负约束最小二乘法的受体模型，在土壤重金属源解析中具有独特的优势，能够有效解析土壤重金属的污染源。该模型假设土壤中重金属的浓度是由多个污染源的贡献叠加而成，通过对土壤样品中重金属浓度的分析，结合源成分谱信息，运用非负约束最小二乘法求解，从而确定各污染源对土壤重金属污染的贡献率。在赵县土壤重金属源解析研究中，应用UNMIX模型的步骤如下：首先，对土壤样品中重金属浓度数据进行全面、细致的收集和整理，确保数据的准确性和完整性。同时，广泛收集赵县地区可能的污染源信息，如工业企业的排放数据、农业生产中化肥农药的使用情况以及交通流量等，以此确定源成分谱。这些源成分谱数据是模型解析的重要依据，其准确性直接影响解析结果的可靠性。其次，根据数据特点和研究需求，合理设置UNMIX模型的参数，如迭代次数、收敛标准等。迭代次数决定了模型计算的精度，收敛标准则确保模型在合理的范围内收敛，避免出现过度拟合或不收敛的情况。然后，将整理好的数据和设置好的参数输入UNMIX模型进行运算。在运算过程中，模型会不断调整计算结果，直到满足设定的收敛标准。最后，对模型输出的结果进行深入分析和验证。通过与实际调查情况对比、结合其他源解析方法的结果进行综合判断等方式，评估解析结果的合理性和可靠性。若发现结果存在异常或不合理之处，需重新检查数据和模型设置，进行必要的调整和优化，以确保解析结果能够真实反映赵县土壤重金属的污染源情况。三、赵县土壤重金属污染评价结果与分析3.1土壤重金属含量统计特征对赵县[X]个土壤样品中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等8种重金属元素的含量进行测定后，得到了详细的统计数据，其结果如表1所示。从平均值来看，赵县土壤中Cr含量平均值为[Cr均值]mg/kg，Cd含量平均值为[Cd均值]mg/kg，Hg含量平均值为[Hg均值]mg/kg，Pb含量平均值为[Pb均值]mg/kg，Ni含量平均值为[Ni均值]mg/kg，Zn含量平均值为[Zn均值]mg/kg，Cu含量平均值为[Cu均值]mg/kg，As含量平均值为[As均值]mg/kg。通过与河北省土壤背景值以及《土壤环境质量农用地土壤污染风险管控标准（试行）》（GB15618-2018）中的筛选值进行对比分析，可以更清晰地了解赵县土壤重金属的污染状况。与河北省土壤背景值相比，赵县土壤中[具体重金属1]的平均含量略高于背景值，这可能暗示着该重金属受到了一定程度的人为活动影响，如工业排放、农业活动等；而[具体重金属2]的平均含量则与背景值相近，表明其含量较为稳定，受外界干扰较小。与GB15618-2018筛选值相比，[具体重金属3]的平均值远低于筛选值，处于安全范围之内；然而，[具体重金属4]的部分样品含量超过了筛选值，存在一定的污染风险，需要引起高度关注。从最大值和最小值来看，各重金属元素的含量范围存在较大差异。其中，[具体重金属5]的最大值与最小值相差最为悬殊，这表明在赵县不同区域，该重金属的含量分布极不均匀，可能与局部的污染源分布有关。例如，在靠近工业区域或交通要道的采样点，[具体重金属5]的含量可能较高；而在远离污染源的农业区域，其含量则相对较低。相比之下，[具体重金属6]的最大值与最小值差距较小，说明该重金属在赵县土壤中的分布相对较为均匀。标准差反映了数据的离散程度，标准差越大，说明数据的离散程度越大，即土壤中重金属含量的空间变异性越大。从统计结果来看，[具体重金属7]的标准差较大，这意味着该重金属在赵县土壤中的含量在空间上存在较大的变异性，可能受到多种因素的影响，如地形、土壤类型、土地利用方式以及污染源的分布等。而[具体重金属8]的标准差较小，表明其在土壤中的含量相对稳定，空间变异性较小。变异系数是标准差与平均值的比值，它消除了数据量纲的影响，能够更准确地反映数据的离散程度。一般来说，变异系数大于15%时，表明数据具有较强的空间变异性。在赵县土壤中，[具体重金属9]的变异系数大于15%，属于强变异性，这进一步证实了该重金属在空间分布上的不均匀性，可能受到强烈的人为活动或特殊的地质条件影响。而[具体重金属10]的变异系数小于15%，属于弱变异性，说明其含量相对稳定，受外界因素影响较小。通过对这些统计参数的深入分析，可以初步了解赵县土壤中重金属的含量特征、分布规律以及可能受到的影响因素，为后续的污染评价和源解析提供重要的数据支持。表1赵县土壤重金属含量统计特征（mg/kg）重金属元素平均值最大值最小值标准差变异系数（%）河北省土壤背景值GB15618-2018筛选值Cr[Cr均值][Cr最大值][Cr最小值][Cr标准差][Cr变异系数][河北省Cr背景值][GB15618Cr筛选值]Cd[Cd均值][Cd最大值][Cd最小值][Cd标准差][Cd变异系数][河北省Cd背景值][GB15618Cd筛选值]Hg[Hg均值][Hg最大值][Hg最小值][Hg标准差][Hg变异系数][河北省Hg背景值][GB15618Hg筛选值]Pb[Pb均值][Pb最大值][Pb最小值][Pb标准差][Pb变异系数][河北省Pb背景值][GB15618Pb筛选值]Ni[Ni均值][Ni最大值][Ni最小值][Ni标准差][Ni变异系数][河北省Ni背景值][GB15618Ni筛选值]Zn[Zn均值][Zn最大值][Zn最小值][Zn标准差][Zn变异系数][河北省Zn背景值][GB15618Zn筛选值]Cu[Cu均值][Cu最大值][Cu最小值][Cu标准差][Cu变异系数][河北省Cu背景值][GB15618Cu筛选值]As[As均值][As最大值][As最小值][As标准差][As变异系数][河北省As背景值][GB15618As筛选值]3.2单因子污染指数评价结果根据单因子污染指数法的计算公式P_i=\frac{C_i}{S_i}，以《土壤环境质量农用地土壤污染风险管控标准（试行）》（GB15618-2018）中的筛选值作为评价标准，计算赵县土壤中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等8种重金属元素的单因子污染指数，计算结果如表2所示。从表2可以看出，赵县土壤中不同重金属元素的单因子污染指数存在明显差异。其中，Cu元素的单因子污染指数最大值达到了[Cu最大污染指数]，超标率为[Cu超标率]，且有[重度污染点数]个重度污染位点（P_i>3），表明Cu元素在赵县部分区域的土壤中污染较为严重。这些重度污染位点可能与当地的工业活动、农业生产方式以及特殊的地质条件有关。例如，某些工业企业在生产过程中可能会排放含有铜的废水、废气和废渣，这些污染物进入土壤后，导致土壤中铜含量升高。农业生产中，不合理地使用含铜的农药、化肥，也可能造成土壤中铜的积累。Cr元素的单因子污染指数范围为[Cr污染指数范围]，平均值为[Cr平均污染指数]，所有样品的污染指数均小于1，说明Cr元素在赵县土壤中整体处于安全水平，未受到明显污染。这可能是由于赵县地区的地质背景中Cr元素含量较低，且工业活动对Cr元素的排放控制较好，没有对土壤造成显著影响。Cd元素的单因子污染指数最大值为[Cd最大污染指数]，平均值为[Cd平均污染指数]，超标率为[Cd超标率]，存在轻度污染的情况。虽然超标率相对较低，但镉是一种毒性较强的重金属，即使是轻度污染也可能对土壤生态系统和农作物生长产生潜在危害。其污染可能与工业排放、污水灌溉以及含镉农药的使用有关。例如，一些化工企业排放的废水中可能含有镉，通过污水灌溉进入农田土壤，导致土壤中镉含量升高。Hg元素的单因子污染指数最大值为[Hg最大污染指数]，平均值为[Hg平均污染指数]，超标率为[Hg超标率]，部分区域存在轻度污染。汞具有挥发性和生物累积性，对人体健康和生态环境危害较大。其污染来源可能包括工业废气排放、含汞农药和化肥的使用以及废旧电池等电子垃圾的不当处理。在赵县，一些小型工业企业可能缺乏有效的废气处理设施，导致汞排放到大气中，随后通过大气沉降进入土壤。Pb元素的单因子污染指数范围较窄，最大值为[Pb最大污染指数]，平均值为[Pb平均污染指数]，所有样品均未超标，表明Pb元素在赵县土壤中的污染程度较轻，处于安全状态。这可能得益于当地对铅污染源的有效管控，如加强对工业企业的监管，减少铅的排放；同时，随着环保意识的提高，含铅汽油的使用逐渐减少，也降低了交通源对土壤铅污染的贡献。Ni元素的单因子污染指数最大值为[Ni最大污染指数]，平均值为[Ni平均污染指数]，未出现超标情况，说明Ni元素在赵县土壤中未受到明显污染，土壤质量良好。这可能与赵县的自然地质条件以及工业布局有关，当地的地质背景中镍含量较低，且工业生产中涉及镍排放的企业较少。Zn元素的单因子污染指数最大值为[Zn最大污染指数]，平均值为[Zn平均污染指数]，所有样品均未超标，表明Zn元素在赵县土壤中处于安全水平，未受到污染。锌是植物生长所必需的微量元素，但过量的锌也会对植物产生毒害作用。赵县土壤中锌含量未超标，可能是由于农业生产中对含锌化肥的使用较为合理，且工业排放对锌的影响较小。As元素的单因子污染指数最大值为[As最大污染指数]，平均值为[As平均污染指数]，未出现超标情况，说明As元素在赵县土壤中未受到污染，土壤环境质量较好。砷是一种类金属元素，具有毒性，其污染来源主要包括工业废水、废气排放以及含砷农药的使用。赵县在工业发展过程中，对砷污染源的控制较为严格，同时农业生产中含砷农药的使用量较少，从而使得土壤中砷含量处于安全范围内。综合以上分析，在赵县土壤中，Cu元素的污染程度最为严重，其次是Cd和Hg元素，存在一定的污染风险；而Cr、Pb、Ni、Zn和As元素整体处于安全水平，未受到明显污染。针对污染较为严重的重金属元素，需要进一步深入研究其污染来源和迁移转化规律，以便采取有效的治理措施，保护土壤生态环境。表2赵县土壤重金属单因子污染指数统计结果重金属元素最大值最小值平均值超标率（%）Cr[Cr最大污染指数][Cr最小污染指数][Cr平均污染指数]0Cd[Cd最大污染指数][Cd最小污染指数][Cd平均污染指数][Cd超标率]Hg[Hg最大污染指数][Hg最小污染指数][Hg平均污染指数][Hg超标率]Pb[Pb最大污染指数][Pb最小污染指数][Pb平均污染指数]0Ni[Ni最大污染指数][Ni最小污染指数][Ni平均污染指数]0Zn[Zn最大污染指数][Zn最小污染指数][Zn平均污染指数]0Cu[Cu最大污染指数][Cu最小污染指数][Cu平均污染指数][Cu超标率]As[As最大污染指数][As最小污染指数][As平均污染指数]03.3内梅罗指数评价结果在内梅罗综合污染指数评价中，首先依据单因子污染指数的计算结果，获取每个采样点处单因子污染指数的最大值P_{i\max}以及单因子污染指数的平均值\overline{P_i}，再代入内梅罗综合污染指数公式P_{综}=\sqrt{\frac{(P_{i\max}^2+\overline{P_i}^2)}{2}}进行计算。通过该公式计算得出赵县各个采样点的内梅罗综合污染指数，其计算结果如表3所示。表3赵县土壤内梅罗综合污染指数统计结果统计参数内梅罗综合污染指数最大值[P综最大值]最小值[P综最小值]平均值[P综平均值]从统计结果来看，赵县土壤内梅罗综合污染指数的最大值为[P综最大值]，出现在[具体点位名称]，表明该点位的土壤污染较为严重；最小值为[P综最小值]，位于[具体点位名称]，说明此点位土壤污染程度较轻。平均值为[P综平均值]，整体处于[整体污染程度描述，如警戒限与轻度污染之间]。依据内梅罗综合污染指数的分级标准，将赵县土壤污染等级进行划分。当P_{综}\leq0.7时，土壤处于安全清洁状态；当0.7<P_{综}\leq1.0时，土壤处于警戒限，尚清洁；当1.0<P_{综}\leq2.0时，土壤受到轻度污染；当2.0<P_{综}\leq3.0时，土壤受到中度污染；当P_{综}>3.0时，土壤受到重度污染。经过统计，处于安全清洁状态的采样点有[安全清洁点数]个，占总采样点数的[安全清洁占比]%；处于警戒限的采样点有[警戒限点数]个，占比为[警戒限占比]%；受到轻度污染的采样点有[轻度污染点数]个，占[轻度污染占比]%；受到中度污染的采样点有[中度污染点数]个，占[中度污染占比]%；受到重度污染的采样点有[重度污染点数]个，占[重度污染占比]%。由此可见，赵县土壤整体以轻度污染和警戒限状态为主，部分区域存在中度和重度污染情况。为更直观地展示赵县土壤重金属污染的空间分布特征，借助地理信息系统（GIS）技术，将内梅罗综合污染指数与地理空间信息相结合，绘制赵县土壤重金属污染空间分布图，如图1所示。从图中可以清晰地看出，赵县东北部的谢庄乡和范庄镇是污染较为集中的区域，这些区域的内梅罗综合污染指数较高，部分地区达到中度甚至重度污染水平。这可能与该地区的工业活动、农业生产方式以及交通状况密切相关。赵县东北部有较多的工业企业，工业生产过程中排放的废气、废水和废渣可能含有大量重金属，这些污染物进入土壤后，导致土壤重金属含量升高。该地区是赵县的主要农业产区，农业生产中大量使用化肥、农药，以及不合理的污水灌溉等，也可能加重土壤重金属污染。交通干道附近的土壤污染程度也相对较高，可能是由于汽车尾气排放、道路扬尘等原因，使得重金属在土壤中积累。相比之下，赵县北部的前大章乡和韩村镇等地的土壤污染程度较轻，大部分区域处于安全清洁或警戒限状态。这些地区工业活动相对较少，农业生产方式较为传统，对土壤环境的影响较小。地形地貌、土壤类型等自然因素也可能对土壤重金属的分布产生影响。例如，土壤的质地、酸碱度、有机质含量等都会影响重金属在土壤中的吸附、解吸和迁移转化，从而导致土壤重金属含量在空间上的差异。综上所述，内梅罗综合污染指数评价结果全面地反映了赵县土壤重金属的污染状况和空间分布特征，为赵县土壤环境保护和污染治理提供了重要的科学依据。通过对污染区域的分析，可以明确重点治理区域，针对不同污染程度的区域采取相应的治理措施，如对重度污染区域进行土壤修复，对轻度污染和警戒限区域加强监测和管理，减少污染物的排放，防止污染进一步加重。[此处插入赵县土壤重金属污染空间分布图]图1赵县土壤重金属污染空间分布图3.4地积累指数评价结果依据地积累指数法的计算公式I_{geo}=\log_2\frac{C_i}{1.5B_i}，以河北省土壤背景值作为地球化学背景值B_i，计算赵县土壤中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等8种重金属元素的地积累指数，计算结果如表4所示。表4赵县土壤重金属地积累指数统计结果重金属元素最大值最小值平均值污染等级占比（无污染，%）污染等级占比（轻度污染，%）污染等级占比（偏中度污染，%）污染等级占比（中度污染，%）污染等级占比（偏重污染，%）污染等级占比（重度污染，%）污染等级占比（严重污染，%）Cr[Cr最大地积累指数][Cr最小地积累指数][Cr平均地积累指数][Cr无污染占比][Cr轻度污染占比][Cr偏中度污染占比][Cr中度污染占比][Cr偏重污染占比][Cr重度污染占比][Cr严重污染占比]Cd[Cd最大地积累指数][Cd最小地积累指数][Cd平均地积累指数][Cd无污染占比][Cd轻度污染占比][Cd偏中度污染占比][Cd中度污染占比][Cd偏重污染占比][Cd重度污染占比][Cd严重污染占比]Hg[Hg最大地积累指数][Hg最小地积累指数][Hg平均地积累指数][Hg无污染占比][Hg轻度污染占比][Hg偏中度污染占比][Hg中度污染占比][Hg偏重污染占比][Hg重度污染占比][Hg严重污染占比]Pb[Pb最大地积累指数][Pb最小地积累指数][Pb平均地积累指数][Pb无污染占比][Pb轻度污染占比][Pb偏中度污染占比][Pb中度污染占比][Pb偏重污染占比][Pb重度污染占比][Pb严重污染占比]Ni[Ni最大地积累指数][Ni最小地积累指数][Ni平均地积累指数][Ni无污染占比][Ni轻度污染占比][Ni偏中度污染占比][Ni中度污染占比][Ni偏重污染占比][Ni重度污染占比][Ni严重污染占比]Zn[Zn最大地积累指数][Zn最小地积累指数][Zn平均地积累指数][Zn无污染占比][Zn轻度污染占比][Zn偏中度污染占比][Zn中度污染占比][Zn偏重污染占比][Zn重度污染占比][Zn严重污染占比]Cu[Cu最大地积累指数][Cu最小地积累指数][Cu平均地积累指数][Cu无污染占比][Cu轻度污染占比][Cu偏中度污染占比][Cu中度污染占比][Cu偏重污染占比][Cu重度污染占比][Cu严重污染占比]As[As最大地积累指数][As最小地积累指数][As平均地积累指数][As无污染占比][As轻度污染占比][As偏中度污染占比][As中度污染占比][As偏重污染占比][As重度污染占比][As严重污染占比]从表4中可以看出，赵县土壤中不同重金属元素的地积累指数存在显著差异，表明它们的污染程度和来源有所不同。其中，Cu元素的地积累指数平均值为[Cu平均地积累指数]，处于中度污染水平。最大值达到了[Cu最大地积累指数]，属于重度污染，这与单因子污染指数和内梅罗指数的评价结果一致，进一步说明Cu元素在赵县部分区域的土壤中污染较为严重。其污染可能主要源于工业活动中含铜废水、废气和废渣的排放，以及农业生产中含铜农药、化肥的不合理使用。在赵县的一些工业园区，金属加工、电子制造等企业在生产过程中会产生大量含铜污染物，如果这些污染物未经有效处理直接排放，很容易导致周边土壤铜含量超标。农业生产中，为了防治病虫害，部分农民可能会过量使用含铜农药，长期积累下来，也会造成土壤铜污染。Cd元素的地积累指数平均值为[Cd平均地积累指数]，处于轻度污染水平。最大值为[Cd最大地积累指数]，达到偏中度污染，表明部分区域存在一定程度的镉污染。镉的毒性较强，即使是轻度污染也可能对土壤生态系统和农作物生长产生潜在危害。其污染来源可能与工业排放、污水灌溉以及含镉农药的使用有关。一些化工企业排放的废水中含有镉，通过污水灌溉进入农田土壤，导致土壤中镉含量升高。含镉农药在农业生产中的使用，也会增加土壤中镉的输入。Hg元素的地积累指数平均值为[Hg平均地积累指数]，处于轻度污染水平。最大值为[Hg最大地积累指数]，达到偏中度污染，说明部分区域存在汞污染。汞具有挥发性和生物累积性，对人体健康和生态环境危害较大。其污染可能与工业废气排放、含汞农药和化肥的使用以及废旧电池等电子垃圾的不当处理有关。在赵县，一些小型工业企业可能缺乏有效的废气处理设施，导致汞排放到大气中，随后通过大气沉降进入土壤。含汞农药和化肥的使用，以及废旧电池等电子垃圾随意丢弃在土壤中，也会导致土壤汞污染。Cr、Pb、Ni、Zn和As元素的地积累指数平均值均小于0，处于无污染水平，表明这些元素在赵县土壤中整体未受到明显污染。Cr元素的地积累指数范围为[Cr地积累指数范围]，所有样品均无污染，说明赵县土壤中铬含量相对稳定，受外界干扰较小。这可能是由于赵县地区的地质背景中Cr元素含量较低，且工业活动对Cr元素的排放控制较好。Pb元素的地积累指数范围较窄，最大值为[Pb最大地积累指数]，所有样品均无污染，表明赵县土壤中铅污染程度较轻，处于安全状态。这可能得益于当地对铅污染源的有效管控，如加强对工业企业的监管，减少铅的排放；同时，随着环保意识的提高，含铅汽油的使用逐渐减少，也降低了交通源对土壤铅污染的贡献。Ni元素的地积累指数最大值为[Ni最大地积累指数]，所有样品均无污染，说明赵县土壤中镍未受到明显污染，土壤质量良好。这可能与赵县的自然地质条件以及工业布局有关，当地的地质背景中镍含量较低，且工业生产中涉及镍排放的企业较少。Zn元素的地积累指数最大值为[Zn最大地积累指数]，所有样品均无污染，表明赵县土壤中锌处于安全水平，未受到污染。锌是植物生长所必需的微量元素，但过量的锌也会对植物产生毒害作用。赵县土壤中锌含量未超标，可能是由于农业生产中对含锌化肥的使用较为合理，且工业排放对锌的影响较小。As元素的地积累指数最大值为[As最大地积累指数]，所有样品均无污染，说明赵县土壤中砷未受到污染，土壤环境质量较好。砷是一种类金属元素，具有毒性，其污染来源主要包括工业废水、废气排放以及含砷农药的使用。赵县在工业发展过程中，对砷污染源的控制较为严格，同时农业生产中含砷农药的使用量较少，从而使得土壤中砷含量处于安全范围内。为了更直观地展示赵县土壤重金属污染的空间分布特征，利用地理信息系统（GIS）技术，将地积累指数与地理空间信息相结合，绘制赵县土壤重金属地积累指数空间分布图，如图2所示。从图中可以清晰地看出，赵县东北部的谢庄乡和范庄镇是Cu、Cd和Hg等重金属污染相对较重的区域。这些区域的地积累指数较高，部分地区达到中度甚至重度污染水平。这可能与该地区的工业活动、农业生产方式以及交通状况密切相关。赵县东北部有较多的工业企业，工业生产过程中排放的废气、废水和废渣可能含有大量重金属，这些污染物进入土壤后，导致土壤重金属含量升高。该地区是赵县的主要农业产区，农业生产中大量使用化肥、农药，以及不合理的污水灌溉等，也可能加重土壤重金属污染。交通干道附近的土壤污染程度也相对较高，可能是由于汽车尾气排放、道路扬尘等原因，使得重金属在土壤中积累。相比之下，赵县北部的前大章乡和韩村镇等地的土壤污染程度较轻，大部分区域处于无污染或轻度污染状态。这些地区工业活动相对较少，农业生产方式较为传统，对土壤环境的影响较小。地形地貌、土壤类型等自然因素也可能对土壤重金属的分布产生影响。例如，土壤的质地、酸碱度、有机质含量等都会影响重金属在土壤中的吸附、解吸和迁移转化，从而导致土壤重金属含量在空间上的差异。综上所述，地积累指数评价结果全面地反映了赵县土壤重金属的污染状况和空间分布特征，为赵县土壤环境保护和污染治理提供了重要的科学依据。通过对污染区域的分析，可以明确重点治理区域，针对不同污染程度的区域采取相应的治理措施，如对重度污染区域进行土壤修复，对轻度污染和无污染区域加强监测和管理，减少污染物的排放，防止污染进一步加重。[此处插入赵县土壤重金属地积累指数空间分布图]图2赵县土壤重金属地积累指数空间分布图3.5潜在生态风险指数评价结果依据潜在生态风险指数法的计算公式RI=\sum_{i=1}^{n}E_r^i=\sum_{i=1}^{n}T_r^i\times\frac{C_i}{C_n^i}，以河北省土壤背景值作为参比含量C_n^i，计算赵县土壤中Cr、Cd、Hg、Pb、Ni、Zn、Cu和As等8种重金属元素的潜在生态风险系数E_r^i以及潜在生态风险指数RI，计算结果如表5所示。表5赵县土壤重金属潜在生态风险系数及风险指数统计结果重金属元素潜在生态风险系数最大值潜在生态风险系数最小值潜在生态风险系数平均值潜在生态风险等级（平均值）潜在生态风险指数最大值潜在生态风险指数最小值潜在生态风险指数平均值潜在生态风险等级（平均值）Cr[Cr最大风险系数][Cr最小风险系数][Cr平均风险系数][Cr平均风险等级][RI最大含Cr][RI最小含Cr][RI平均含Cr][RI平均等级含Cr]Cd[Cd最大风险系数][Cd最小风险系数][Cd平均风险系数][Cd平均风险等级][RI最大含Cd][RI最小含Cd][RI平均含Cd][RI平均等级含Cd]Hg[Hg最大风险系数][Hg最小风险系数][Hg平均风险系数][Hg平均风险等级][RI最大含Hg][RI最小含Hg][RI平均含Hg][RI平均等级含Hg]Pb[Pb最大风险系数][Pb最小风险系数][Pb平均风险系数][Pb平均风险等级][RI最大含Pb][RI最小含Pb][RI平均含Pb][RI平均等级含Pb]Ni[Ni最大风险系数][Ni最小风险系数][Ni平均风险系数][Ni平均风险等级][RI最大含Ni][RI最小含Ni][RI平均含Ni][RI平均等级含Ni]Zn[Zn最大风险系数][Zn最小风险系数][Zn平均风险系数][Zn平均风险等级][RI最大含Zn][RI最小含Zn][RI平均含Zn][RI平均等级含Zn]Cu[Cu最大风险系数][Cu最小风险系数][Cu平均风险系数][Cu平均风险等级][RI最大含Cu][RI最小含Cu][RI平均含Cu][RI平均等级含Cu]As[As最大风险系数][As最小风险系数][As平均风险系数][As平均风险等级][RI最大含As][RI最小含As][RI平均含As][RI平均等级含As]综合[综合最大风险系数][综合最小风险系数][综合平均风险系数]-[RI最大值][RI最小值][RI平均值][RI平均等级]从表5可以看出，赵县土壤中不同重金属元素的潜在生态风险系数和风险指数存在明显差异。其中，Hg元素的潜在生态风险系数平均值为[Hg平均风险系数]，处于中等风险水平，最大值达到了[Hg最大风险系数]，属于较高风险，表明部分区域存在较高的汞污染生态风险。汞具有挥发性和生物累积性，对人体健康和生态环境危害较大，其污染可能与工业废气排放、含汞农药和化肥的使用以及废旧电池等电子垃圾的不当处理有关。在赵县，一些小型工业企业可能缺乏有效的废气处理设施，导致汞排放到大气中，随后通过大气沉降进入土壤。含汞农药和化肥的使用，以及废旧电池等电子垃圾随意丢弃在土壤中，也会导致土壤汞污染。Cd元素的潜在生态风险系数平均值为[Cd平均风险系数]，处于中等风险水平，最大值为[Cd最大风险系数]，达到较高风险，说明部分区域存在一定程度的镉污染生态风险。镉的毒性较强，即使是中等风险也可能对土壤生态系统和农作物生长产生潜在危害。其污染来源可能与工业排放、污水灌溉以及含镉农药的使用有关。一些化工企业排放的废水中含有镉，通过污水灌溉进入农田土壤，导致土壤中镉含量升高。含镉农药在农业生产中的使用，也会增加土壤中镉的输入。Cu元素的潜在生态风险系数平均值为[Cu平均风险系数]，处于低风险水平，最大值为[Cu最大风险系数]，达到中等风险，表明部分区域存在一定的铜污染生态风险。虽然整体风险较低，但在污染较重的区域仍需关注。其污染可能主要源于工业活动中含铜废水、废气和废渣的排放，以及农业生产中含铜农药、化肥的不合理使用。在赵县的一些工业园区，金属加工、电子制造等企业在生产过程中会产生大量含铜污染物，如果这些污染物未经有效处理直接排放，很容易导致周边土壤铜含量超标。农业生产中，为了防治病虫害，部分农民可能会过量使用含铜农药，长期积累下来，也会造成土壤铜污染。Cr、Pb、Ni、Zn和As元素的潜在生态风险系数平均值均较低，处于低风险水平，表明这些元素在赵县土壤中整体对生态环境的潜在危害较小。Cr元素的潜在生态风险系数范围为[Cr风险系数范围]，所有样品均处于低风险，说明赵县土壤中铬对生态环境的影响较小。这可能是由于赵县地区的地质背景中Cr元素含量较低，且工业活动对Cr元素的排放控制较好。Pb元素的潜在生态风险系数范围较窄，最大值为[Pb最大风险系数]，所有样品均处于低风险，表明赵县土壤中铅对生态环境的潜在危害较轻。这可能得益于当地对铅污染源的有效管控，如加强对工业企业的监管，减少铅的排放；同时，随着环保意识的提高，含铅汽油的使用逐渐减少，也降低了交通源对土壤铅污染的贡献。Ni元素的潜在生态风险系数最大值为[Ni最大风险系数]，所有样品均处于低风险，说明赵县土壤中镍对生态环境未产生明显危害，土壤质量良好。这可能与赵县的自然地质条件以及工业布局有关，当地的地质背景中镍含量较低，且工业生产中涉及镍排放的企业较少。Zn元素的潜在生态风险系数最大值为[Zn最大风险系数]，所有样品均处于低风险，表明赵县土壤中锌对生态环境处于安全水平，未受到污染。锌是植物生长所必需的微量元素，但过量的锌也会对植物产生毒害作用。赵县土壤中锌含量未超标，可能是由于农业生产中对含锌化肥的使用较为合理，且工业排放对锌的影响较小。As元素的潜在生态风险系数最大值为[As最大风险系数]，所有样品均处于低风险，说明赵县土壤中砷对生态环境未造成污染，土壤环境质量较好。砷是一种类金属元素，具有毒性，其污染来源主要包括工业废水、废气排放以及含砷农药的使用。赵县在工业发展过程中，对砷污染源的控制较为严格，同时农业生产中含砷农药的使用量较少，从而使得土壤中砷含量处于安全范围内。综合考虑8种重金属元素，赵县土壤潜在生态风险指数的平均值为[RI平均值]，处于中等风险水平。最大值为[RI最大值]，达到较高风险水平，最小值为[RI最小值]，处于低风险水平。这表明赵县部分区域的土壤存在较高的潜在生态风险，需要引起重视。为了更直观地展示赵县土壤潜在生态风险的空间分布特征，利用地理信息系统（GIS）技术，将潜在生态风险指数与地理空间信息相结合，绘制赵县土壤潜在生态风险指数空间分布图，如图3所示。从图中可以清晰地看出，赵县东北部的谢庄乡和范庄镇是潜在生态风险相对较高的区域，这些区域的潜在生态风险指数较高，部分地区达到较高风险水平。这可能与该地区的工业活动、农业生产方式以及交通状况密切相关。赵县东北部有较多的工业企业，工业生产过程中排放的废气、废水和废渣可能含有大量重金属，这些污染物进入土壤后，导致土壤重金属含量升高，从而增加了潜在生态风险。该地区是赵县的主要农业产区，农业生产中大量使用化肥、农药，以及不合理的污水灌溉等，也可能加重土壤重金属污染，进而提高潜在生态风险。交通干道附近的土壤潜在生态风险程度也相对较高，可能是由于汽车尾气排放、道路扬尘等原因，使得重金属在土壤中积累，增加了潜在生态风险。相比之下，赵县北部的前大章乡和韩村镇等地的土壤潜在生态风险程度较轻，大部分区域处于低风险或中等风险状态。这些地区工业活动相对较少，农业生产方式较为传统，对土壤环境的影响较小。地形地貌、土壤类型等自然因素也可能对土壤潜在生态风险的分布产生影响。例如，土壤的质地、酸碱度、有机质含量等都会影响重金属在土壤中的吸附、解吸和迁移转化，从而导致土壤潜在生态风险在空间上的差异。综上所述，潜在生态风险指数评价结果全面地反映了赵县土壤重金属对生态环境的潜在危害状况和空间分布特征，为赵县土壤环境保护和污染治理提供了重要的科学依据。通过对潜在生态风险较高区域的分析，可以明确重点防控区域，针对不同风险程度的区域采取相应的治理措施，如对较高风险区域进行土壤修复，对中等风险和低风险区域加强监测和管理，减少污染物的排放，防止潜在生态风险进一步加剧。[此处插入赵县土壤潜在生态风险指数空间分布图]图3赵县土壤潜在生态风险指数空间分布图四、赵县土壤重金属源解析结果与分析4.1UNMIX模型解析结果将赵县土壤样品中重金属浓度数据以及收集到的源成分谱信息输入UNMIX模型进行运算，经过多次调试和优化参数，确保模型收敛且结果合理可靠。模型运行结果表明，赵县土壤重金属主要来源于[X]个不同的污染源，各污染源对土壤重金属的贡献率如表6所示。表6UNMIX模型解析出的赵县土壤重金属污染源及贡献率污染源贡献率（%）主要贡献的重金属元素污染源1[贡献率1][主要贡献重金属1、主要贡献重金属2等]污染源2[贡献率2][主要贡献重金属3、主要贡献重金属4等]污染源3[贡献率3][主要贡献重金属5、主要贡献重金属6等]………………污染源1的贡献率为[贡献率1]%，是赵县土壤重金属污染的主要来源之一。该污染源主要贡献的重金属元素为[主要贡献重金属1、主要贡献重金属2等]。通过对赵县的实际情况分析，结合工业布局和生产活动特点，推测污染源1可能与赵县的工业生产活动密切相关。赵县有多个工业企业，在生产过程中，如金属冶炼、化工制造等行业会产生含有重金属的废气、废水和废渣。金属冶炼企业在矿石熔炼过程中，会释放出大量含有铜、铅、锌等重金属的废气，这些废气在大气中扩散后，通过大气沉降进入土壤，导致土壤中相应重金属含量升高。化工企业排放的废水中可能含有汞、镉等重金属，若未经有效处理直接排放，会通过地表径流或地下水渗透进入土壤，造成土壤污染。废渣的随意堆放也会使其中的重金属逐渐释放到周围土壤中，加重土壤污染程度。污染源2的贡献率为[贡献率2]%，是另一重要污染源。其主要贡献的重金属元素为[主要贡献重金属3、主要贡献重金属4等]。考虑到赵县作为农业大县的实际情况，该污染源可能主要来自农业活动。在农业生产中，大量使用化肥、农药和农膜。化肥中可能含有一定量的重金属，如磷肥中常含有镉，长期大量施用磷肥会导致土壤中镉含量增加。农药中也可能含有重金属，如含铜、锌的杀菌剂，在防治病虫害的同时，会使土壤中铜、锌含量升高。农膜的使用虽然对农作物生长有一定的促进作用，但一些农膜中含有铅、镉等重金属，随着农膜的老化和破碎，这些重金属会逐渐释放到土壤中，对土壤环境造成污染。污染源3的贡献率为[贡献率3]%，主要贡献的重金属元素为[主要贡献重金属5、主要贡献重金属6等]。根据赵县的交通状况和相关研究，该污染源可能与交通-污灌有关。赵县交通网络较为发达，汽车尾气排放和道路扬尘是土壤重金属污染的重要来源之一。汽车尾气中含有铅、锌、镉等重金属，随着汽车的行驶，尾气排放到大气中，部分重金属会通过大气沉降进入土壤。道路扬尘中也可能含有重金属，在风力作用下，扬尘中的重金属会扩散到周围土壤中。污水灌溉也是导致土壤重金属污染的重要因素。赵县部分地区存在利用未经处理或处理不达标污水进行灌溉的情况，这些污水中含有大量重金属，如汞、镉、铅等，长期灌溉会使土壤中重金属含量不断积累，造成土壤污染。4.2各污染源对重金属的贡献分析通过UNMIX模型解析出的不同污染源对土壤中重金属的贡献存在显著差异，这种差异对于深入理解赵县土壤重金属污染的形成机制和制定针对性的防治措施具有重要意义。“土壤源”对土壤重金属的贡献率为[贡献率数值]%，在自然因素主导下，其对Cr、Ni等重金属元素的贡献较为突出。从地质背景来看，赵县位于太行山东麓中段冲洪积形成的平原区，成土母质的特性在很大程度上决定了土壤中重金属的本底含量。研究表明，不同的成土母质其矿物组成和化学元素含量存在差异，赵县的成土母质中可能富含Cr、Ni等元素，从而使得“土壤源”对这两种重金属元素的贡献较为明显。Cr和Ni在土壤中的化学性质相对稳定，其迁移性较弱，在长期的成土过程中，主要受地质因素的影响，因此“土壤源”成为它们的主要来源。“农药源”对土壤重金属的贡献率为[贡献率数值]%，在农业活动的影响下，对Cu、Zn等重金属元素的贡献较为显著。在赵县的农业生产中，农药的使用较为普遍。据相关调查数据显示，部分农药中含有一定量的Cu、Zn等重金属元素。例如，一些含铜的杀菌剂在防治农作物病害方面应用广泛，长期大量使用此类农药会导致土壤中Cu含量增加。有机锌类农药也被广泛用于农业生产，其在土壤中的累积使得“农药源”对Zn元素的贡献较为突出。从空间分布上看，“农药源”对赵县东北部地区的影响较大。这是因为该地区是赵县的主要农业产区，农作物种植面积广，农药使用量相对较大。长期的农药施用使得土壤中Cu、Zn等重金属逐渐积累，从而加重了该地区的土壤重金属污染程度。“交通-污灌源”对土壤重金属的贡献率最高，达到[贡献率数值]%，是赵县土壤重金属污染的主要污染源。在交通活动和污水灌溉的双重作用下，对Cd、Hg、Pb等重金属元素的贡献尤为明显。赵县交通网络较为发达，汽车尾气排放和道路扬尘是土壤重金属污染的重要来源之一。汽车尾气中含有大量的Cd、Hg、Pb等重金属，随着汽车的行驶，尾气排放到大气中，部分重金属会通过大气沉降进入土壤。相关研究表明，在交通干道附近，土壤中Cd、Hg、Pb等重金属的含量明显高于其他区域，这充分说明了交通活动对土壤重金属污染的影响。污水灌溉也是导致土壤重金属污染的重要因素。赵县部分地区存在利用未经处理或处理不达标污水进行灌溉的情况，这些污水中含有大量的Cd、Hg、Pb等重金属。长期的污水灌溉使得土壤中重金属含量不断积累，造成了严重的土壤污染。从空间分布来看，“交通-污灌源”对赵县的整体土壤环境都有较大影响，尤其是在交通干道沿线和污水灌溉区域，土壤中Cd、Hg、Pb等重金属含量较高，污染程度较为严重。4.3污染源成因探讨赵县土壤重金属污染源的形成是多种因素共同作用的结果，与当地的工农业生产、交通状况以及灌溉方式等密切相关。在工业生产方面，赵县初步形成了以医药产业为主导，现代食品加工、纺织服装、农机装备制造等特色产业协同发展的经济体系。赵县经济开发区实行“一区两园”模式，西区以医药制造业为主导产业，巩固提升现代食品、纺织服装、装备制造等传统产业；东区以医药制造为主，辅之发展相关的精细化工产业。这些工业企业在生产过程中会产生大量含有重金属的污染物。金属冶炼企业在矿石熔炼过程中，会释放出含有铜、铅、锌等重金属的废气，这些废气在大气中扩散后，通过大气沉降进入土壤，导致土壤中相应重金属含量升高。化工企业排放的废水中可能含有汞、镉等重金属，若未经有效处理直接排放，会通过地表径流或地下水渗透进入土壤，造成土壤污染。工业废渣的随意堆放也会使其中的重金属逐渐释放到周围土壤中，加重土壤污染程度。农业活动也是赵县土壤重金属污染的重要来源。赵县是传统的农业大县，全县耕地面积58.45万亩，园地20.7万亩，主要农作物有小麦、玉米、棉花、花生等，特产雪花梨、鸭梨等。在农业生产中，大量使用化肥、农药和农膜。化肥中可能含有一定量的重金属，如磷肥中常含有镉，长期大量施用磷肥会导致土壤中镉含量增加。据统计，赵县每年磷肥的施用量较大，这无疑增加了土壤中镉污染的风险。