泾惠渠灌区浅层地下水重金属污染特征与健康风险的深度剖析

泾惠渠灌区浅层地下水重金属污染特征与健康风险的深度剖析一、引言1.1研究背景与意义泾惠渠灌区位于陕西省关中平原中部，是陕西重要的粮农基地，承担着泾阳、三原、高陵、临潼、阎良、富平6县（区）145万亩农田的灌溉任务。该灌区以全省2.5%的耕地，生产出全省5.8%的粮食，年提供商品粮1.9亿公斤，素有“关中白菜心”之称，在陕西省乃至全国的农业生产和经济建设中都发挥着重要的支柱作用。引泾灌溉历史最早可追溯到公元前246年的秦郑国渠，距今已有2200多年，历经汉白渠、唐郑白渠等历代发展，1932年重修后改称泾惠渠并沿用至今，堪称“水利天然博物馆”。然而，近年来随着灌区经济的迅猛发展，城市化进程加快和人类活动增加，对水资源的需求不断增长，地下水开采量逐年递增，一系列环境问题随之而来。一方面，工业废水、生活污水的排放以及农业面源污染，导致该地区地下水环境受到严重污染，灌区浅层地下水质量整体较差，且地下水盐分呈逐年上升趋势。据调查统计，2023年灌区微咸水、咸水面积占比96.57%，地下水盐分较2009年增加1.5-3倍，土壤盐分较2009年增加1-2倍。《陕西统计年鉴》（2022年）数据显示，灌区所在部分地区化学药剂、化肥、农用塑料薄膜使用量均超全省平均水平，2018年硝酸盐超出地下水水质Ⅲ类标准的面积高达79.58%，超出生活饮用水卫生标准的面积达99.4%，灌区土壤等均存在不同程度的重金属污染。另一方面，地下水的超采引起水位下降，降落漏斗区不断扩大，导致地下水位大幅下降，灌区内出现降落漏斗等一系列环境地质问题，对地下水水质产生了负面影响。重金属具有毒性大、难降解、易在生物体内富集等特点，一旦进入地下水，会在水体、土壤和生物体内不断迁移转化，通过食物链进入人体，对人体健康造成潜在威胁。长期暴露于被重金属污染的地下水环境中，可能会导致人体出现各种疾病，如神经系统损伤、癌症、肾脏疾病等。例如，As是一种明确的致癌物质，长期摄入会增加患皮肤癌、肺癌等癌症的风险；Cr（Ⅵ）具有强氧化性和毒性，对人体的皮肤、呼吸道、胃肠道等都有损害作用。因此，开展泾惠渠灌区浅层地下水中重金属的健康风险评价具有重要的现实意义。准确评估泾惠渠灌区浅层地下水中重金属对人体健康的潜在风险，有助于全面了解该地区的水环境质量状况和健康风险水平。这不仅能够为当地政府和相关部门制定科学合理的水资源保护政策、水污染防治措施提供重要的科学依据，促进水资源的合理利用和管理，还能提高公众对地下水污染问题的认识和重视程度，增强公众的环保意识和健康意识，从而有效保护当地居民的身体健康和生态环境安全，保障灌区的可持续发展。1.2国内外研究现状在地下水重金属污染研究方面，国内外学者已取得了丰富的成果。国外研究起步较早，在重金属污染的来源解析、迁移转化规律等方面进行了深入探究。例如，美国地质调查局（USGS）对全国多个地区的地下水进行长期监测，分析重金属的来源和分布特征，发现工业排放、农业活动和矿业开采是地下水重金属污染的主要来源。在欧洲，研究人员利用数值模型模拟重金属在地下水中的迁移转化过程，为污染治理提供理论支持。国内学者在地下水重金属污染研究方面也做了大量工作。众多学者对不同地区的地下水进行了详细的调查分析，揭示了重金属污染的现状和特征。比如，通过对华北平原、珠江三角洲等地区的研究发现，部分区域地下水存在不同程度的重金属超标现象，且与当地的经济发展模式和人类活动密切相关。在健康风险评价领域，国外发展较为成熟，美国环境保护署（USEPA）早在20世纪80年代就提出了一套完整的健康风险评价体系，包括危害识别、暴露评估、剂量-反应评估和风险表征四个步骤，被广泛应用于各种环境污染物的健康风险评价。欧洲一些国家也建立了适合本国国情的健康风险评价模型和方法，注重多污染物联合作用和人群易感性差异对健康风险的影响。国内健康风险评价研究虽然起步相对较晚，但近年来发展迅速。学者们在借鉴国外先进经验的基础上，结合我国实际情况，对健康风险评价模型进行了改进和完善。如针对不同地区的环境特点和人群暴露特征，建立了本地化的健康风险评价模型，使其更符合我国的实际情况。同时，在评价方法上不断创新，将多种分析技术和手段应用于健康风险评价中，提高了评价的准确性和可靠性。然而，针对泾惠渠灌区浅层地下水中重金属的健康风险评价研究相对较少。已有的研究主要集中在灌区地下水的水质评价、水位变化及灌溉效益等方面，对于地下水中重金属的污染特征及对人体健康的潜在风险研究不够深入。虽然张艳、徐斌、刘秀花等人对陕西省泾惠渠灌区地下水污染与人体健康风险进行了评价，分析了As、Cr（Ⅵ）和NO₃⁻-N等污染物的风险状况，但对于其他重金属元素的研究仍有欠缺，且在评价过程中可能未充分考虑灌区复杂的环境因素和人群暴露途径的多样性。因此，有必要开展系统的研究，全面评估泾惠渠灌区浅层地下水中重金属的健康风险，为该地区的水资源保护和居民健康保障提供科学依据。1.3研究内容与方法1.3.1研究内容泾惠渠灌区浅层地下水重金属含量及分布特征：系统采集泾惠渠灌区浅层地下水样本，运用先进的分析仪器和方法，准确测定其中As、Cr、Cd、Pb、Hg等重金属元素的含量。通过数据分析，深入探究重金属含量在空间上的分布规律，包括不同区域、不同深度的含量差异，绘制详细的重金属含量空间分布图，直观展示其分布特征。泾惠渠灌区浅层地下水重金属来源解析：综合运用多元统计分析方法，如主成分分析、因子分析等，对重金属含量数据及相关环境因素数据进行处理，初步识别重金属的潜在来源。结合同位素示踪技术，利用重金属元素的同位素组成特征，追溯其具体来源，明确各来源对地下水重金属污染的贡献比例，为针对性治理提供依据。泾惠渠灌区浅层地下水重金属健康风险评价：基于美国环境保护署（USEPA）推荐的健康风险评价模型，充分考虑灌区居民的生活习惯、用水方式等因素，准确确定暴露参数，对灌区浅层地下水中重金属对人体健康的潜在风险进行全面评估。分别计算致癌风险和非致癌风险，明确不同重金属元素、不同暴露途径（饮水、皮肤接触等）对不同人群（成人、儿童等）的健康风险程度，筛选出高风险区域和关键风险因子。泾惠渠灌区浅层地下水重金属污染防治建议：根据重金属污染特征和健康风险评价结果，从源头控制、过程阻断、末端治理等环节出发，提出科学合理的污染防治建议。制定严格的工业排放标准和农业面源污染控制措施，减少重金属的排放；加强对重点污染源的监管，建立完善的监测体系，实时掌握污染动态；推广绿色环保的修复技术，对已污染的地下水进行有效治理，逐步改善地下水环境质量。同时，结合当地实际情况，制定切实可行的水资源保护规划，保障居民用水安全。1.3.2研究方法采样与分析方法：在泾惠渠灌区范围内，综合考虑地形地貌、土地利用类型、人口分布以及现有监测井的分布情况，科学合理地设置地下水采样点，确保样本具有代表性。使用专业的采样设备，严格按照相关标准和规范采集浅层地下水样本，并及时将样本送往实验室进行分析。采用电感耦合等离子体质谱仪（ICP-MS）等先进仪器测定地下水中重金属元素的含量，确保分析结果的准确性和可靠性。同时，对水样的pH值、溶解氧、电导率等常规水质参数进行同步测定，为后续分析提供基础数据。污染评价方法：运用单因子指数法，分别计算各重金属元素的污染指数，直观判断每种重金属的污染程度。采用内梅罗综合污染指数法，综合考虑多种重金属的协同作用，全面评价地下水的污染状况，确定污染等级。结合地质统计学方法，利用克里金插值等技术，对重金属含量数据进行空间插值处理，绘制重金属含量的空间分布等值线图，清晰展示其空间变化趋势，分析污染的空间分布特征。健康风险评价方法：依据USEPA推荐的健康风险评价模型，结合灌区实际情况，确定暴露剂量、暴露频率、暴露时间等关键参数。对于致癌风险，采用致癌风险模型计算终生致癌风险；对于非致癌风险，利用危害商模型计算危害商值。通过蒙特卡罗模拟等不确定性分析方法，评估风险评价结果的不确定性，使评价结果更加科学合理，准确反映实际健康风险水平。二、泾惠渠灌区概况2.1自然地理条件泾惠渠灌区位于陕西省关中平原中部，地处东经108°4′-109°2′，北纬34°2′-34°5′之间，是一个从泾河自流引水的大（Ⅱ）型灌区。灌区北依仲山和黄土台原，西、南、东三面分别被泾河、渭河、石川河环绕，清河则自西向东穿流而过。灌区东西长度约为70km，南北最宽处达20km，总面积约1180km²。其地势呈现出西北高、东南低的态势，海拔高程处于350-450m之间，地面坡降在1/300-1/600之间，属于典型的北方平原灌区。这种地势特征使得水流能够在重力作用下较为顺畅地自流灌溉，为农业生产提供了良好的水利条件，但也可能导致地下水在流动过程中受到地形影响，出现局部的水位变化和水力坡度差异。该地区属于大陆性半干旱季风气候区，具有显著的气候特点。多年平均降水量为538.9mm，降水主要集中在夏季，且多以暴雨形式出现，降水的时空分布不均，容易引发季节性干旱和洪涝灾害，对地下水的补给产生较大影响。例如，在降水集中的时段，大量雨水可能迅速形成地表径流，来不及充分下渗补充地下水；而在干旱季节，地下水的蒸发和开采量大于补给量，导致水位下降。年蒸发量高达1212mm，远大于降水量，加剧了水资源的短缺状况，使得地下水处于持续的消耗状态。总日照时数为2200h，充足的光照有利于农作物的光合作用和生长发育，但也会增加土壤水分的蒸发，间接影响地下水与土壤水之间的水分交换。多年平均气温为13.4℃，气温的变化会影响土壤水分的物理性质和地下水的水温，进而对地下水的化学性质和微生物活动产生一定作用。在水文方面，泾河是灌区的主要水源，发源于宁夏六盘山东麓，全长455.1km，流域面积45421km²。泾河水量丰富，为灌区提供了稳定的灌溉水源，但河水的含沙量较大，在引水灌溉过程中容易造成渠道淤积，影响灌溉效率和水质，同时也可能携带一些泥沙中的污染物进入地下水系统。灌区除了泾河外，还有石川河、清河等河流，这些河流相互连通，构成了复杂的地表水系网络。地表水与地下水之间存在着密切的水力联系，在丰水期，地表水会通过渗漏等方式补给地下水；而在枯水期，地下水则会向地表水排泄，维持河流水量。此外，灌区还分布着一些水库和池塘，如中型水库2座（张家山水库、西郊水库），小型水库1座（贺兰水库），总库容达4790.7万m³。这些水利设施在调节地表水资源、补充地下水以及保障农业灌溉用水等方面发挥着重要作用，通过合理的蓄放水操作，可以有效调控地表水与地下水的水位和水量平衡。2.2社会经济状况泾惠渠灌区涉及西安、咸阳、渭南三市的临潼、阎良、高陵、泾阳、三原、富平6县（区），人口众多，据不完全统计，灌区总人口达数百万。众多的人口对水资源的需求巨大，在日常生活、农业灌溉和工业生产等方面都离不开水，这使得地下水的开采量不断增加，从而对地下水环境产生了较大压力。例如，随着人口的增长，生活污水的排放量也相应增加，如果处理不当，这些污水中的污染物可能会渗入地下，对浅层地下水造成污染，其中就可能包含重金属污染物。在产业结构方面，灌区以农业为主导产业，是陕西省重要的粮棉油生产基地。农作物以小麦、玉米、棉花、蔬菜等为主，复种指数在1.85以上，农业生产高度依赖灌溉水源，每年用于灌溉的水量占总用水量的较大比例。长期的农业灌溉活动，使得地下水与土壤水之间频繁进行物质交换，土壤中的一些重金属元素可能会随着灌溉水的下渗进入地下水，导致地下水中重金属含量升高。同时，为了提高农作物产量，农民大量使用化肥、农药和农用塑料薄膜，这些农业投入品中的重金属成分，如磷肥中的镉、农药中的汞等，在长期的农业生产过程中，通过地表径流和淋溶作用进入地下水系统，成为地下水重金属污染的重要来源。灌区的工业发展也较为迅速，形成了以食品加工、建材、机械制造、化工等为主的工业体系。这些工业企业在生产过程中会产生大量的工业废水，其中含有多种重金属污染物，如化工企业排放的废水中可能含有汞、镉、铅等重金属，机械制造企业的废水可能含有铬等重金属。如果这些工业废水未经有效处理就直接排放，会对地表水和土壤造成污染，进而通过渗漏等方式影响浅层地下水水质。部分工业企业的选址靠近河流或地下水补给区，其排放的污染物更容易进入地下水系统，加剧了地下水重金属污染的风险。2.3水文地质条件泾惠渠灌区地处关中盆地中部，在地质构造上位于鄂尔多斯台向斜南缘渭河断陷的一部分，是一个新生代以来的沉降盆地。其含水层结构较为复杂，主要由第四系松散堆积物组成，自下而上可分为下更新统（Q1）、中更新统（Q2）、上更新统（Q3）和全新统（Q4）。下更新统（Q1）主要为一套湖相沉积的黏土、粉质黏土夹砂层，厚度较大，一般在100-300m之间，是深层地下水的主要含水层。中更新统（Q2）为冲洪积成因的砂卵砾石层与黏土、粉质黏土互层，厚度在50-150m之间，是灌区浅层地下水的重要含水层之一，其透水性较好，富水性较强，能够储存和传输大量的地下水。上更新统（Q3）主要由冲洪积的粉细砂、中粗砂及粉质黏土组成，厚度相对较薄，一般在20-80m之间，与中更新统含水层存在水力联系，共同构成了浅层地下水的含水层系统。全新统（Q4）主要分布在地表，为近代河流冲积物，以砂、砾石和粉质黏土为主，厚度较薄，一般在5-20m之间，对浅层地下水的补给和径流有一定影响。灌区地下水类型主要为孔隙水，根据埋藏条件和水力特征，可进一步分为浅层潜水和深层承压水。浅层潜水主要赋存于第四系上更新统和全新统地层中，含水层岩性以粉细砂、中粗砂为主，水位埋深受地形、降水、灌溉等因素影响较大。在灌区的不同区域，浅层潜水水位埋深存在明显差异。一般来说，在地势较高的区域，如北部的黄土台原地区，水位埋深较大，可达10-30m；而在地势较低的区域，如靠近河流和灌区中心的平原地带，水位埋深较浅，多在3-10m之间。深层承压水主要赋存于第四系下更新统和中更新统的砂卵砾石层中，由于其上覆有较厚的黏性土层作为隔水层，承压水具有一定的水头压力。深层承压水水位埋深相对较稳定，但近年来随着开采量的增加，水位也呈下降趋势。地下水的补给来源主要包括大气降水入渗、地表水渗漏、灌溉水回渗以及侧向径流补给等。大气降水是灌区地下水的重要补给来源之一，虽然该地区年降水量相对较少，但降水集中在夏季，且多以暴雨形式出现，在短时间内能够形成较大的地表径流，部分雨水通过地表入渗补给地下水。根据相关研究，大气降水入渗补给系数在0.1-0.3之间。地表水渗漏也是重要的补给方式，泾河及其支流作为灌区的主要地表水源，河水在流动过程中通过河床渗漏补给地下水。灌区的渠道渗漏也较为明显，由于部分渠道为土渠，防渗性能较差，灌溉水在输送过程中大量渗漏，成为地下水的重要补给源。据估算，渠道渗漏补给量占地下水总补给量的30%-50%。灌溉水回渗同样不可忽视，灌区农业灌溉用水量较大，灌溉水在田间入渗后，一部分形成土壤水，另一部分下渗补给地下水。随着灌溉技术的改进，灌溉水回渗量有所减少，但仍是地下水补给的重要组成部分。侧向径流补给主要来自灌区周边山区的地下水径流，由于山区地势较高，地下水在重力作用下向灌区侧向流动，补给灌区地下水，但侧向径流补给量相对较小。在径流方面，地下水的流动方向总体上与地形坡度一致，由西北向东南流动。在灌区的不同区域，由于含水层岩性、厚度以及水力坡度的差异，地下水的径流速度有所不同。在含水层透水性较好、水力坡度较大的区域，如靠近河流的砂质平原地区，地下水径流速度较快，一般在0.5-2m/d之间；而在含水层透水性较差、水力坡度较小的区域，如黄土台原地区，地下水径流速度较慢，多在0.1-0.5m/d之间。地下水的排泄方式主要有蒸发、人工开采和侧向径流排泄。蒸发排泄主要发生在浅层潜水，在干旱的气候条件下，浅层潜水通过土壤孔隙蒸发进入大气，尤其是在夏季高温时段，蒸发量较大。人工开采是灌区地下水排泄的主要方式之一，随着工农业生产和生活用水需求的增加，地下水开采量不断增大。目前，灌区机井数量众多，开采深度涵盖浅层潜水和深层承压水，过度开采导致部分区域地下水位下降，形成降落漏斗。侧向径流排泄则是指灌区地下水向周边地势较低的区域排泄，与周边地区的地下水进行水量交换，但排泄量相对较小。这种水文地质条件对重金属的迁移转化有着重要作用。首先，含水层的岩性和结构影响着重金属的吸附和解吸过程。例如，黏土和粉质黏土具有较大的比表面积和阳离子交换容量，能够吸附大量的重金属离子，从而减缓其在地下水中的迁移速度。而砂质含水层透水性好，重金属离子在其中的迁移相对较快。其次，地下水的水位变化和水力坡度决定了重金属的迁移方向和速度。在水位上升期，地下水的流速加快，可能携带更多的重金属离子迁移；而在水位下降期，部分重金属可能会因吸附作用而在含水层中沉淀下来。此外，地下水的补给和排泄过程也会影响重金属的分布。如灌溉水回渗可能会将土壤中的重金属带入地下水，而蒸发排泄则会使地下水中的重金属浓度相对升高。三、研究方法3.1样品采集与分析为全面、准确地掌握泾惠渠灌区浅层地下水中重金属的含量及分布特征，本研究在采样点设置上遵循科学合理的原则。综合考虑灌区的地形地貌、土地利用类型、人口分布以及现有监测井的分布情况，在整个灌区范围内共设置了[X]个采样点。其中，在地势较高的黄土台原区域设置了[X1]个采样点，这些区域地下水的补给和径流条件相对复杂，且受人类活动影响相对较小，有助于分析自然因素对重金属含量的影响；在靠近河流的平原区域设置了[X2]个采样点，河流附近地表水与地下水的水力联系密切，且人类活动频繁，工业、农业和生活污染排放集中，能较好地反映这些因素对地下水重金属含量的综合影响；在灌区中心的农业灌溉集中区域设置了[X3]个采样点，农业灌溉活动是灌区的主要用水方式，大量的灌溉水回渗可能携带土壤中的重金属进入地下水，该区域的采样点可重点研究农业活动对地下水中重金属的影响。在采样时间的选择上，考虑到泾惠渠灌区的气候特点和地下水的动态变化规律，分别在丰水期（7-9月）和枯水期（1-3月）进行采样。丰水期降水较多，地表径流和灌溉水量大，可能会对地下水中重金属的含量和分布产生稀释或冲刷作用；枯水期降水稀少，地下水主要靠侧向径流和少量的灌溉水回渗补给，水位相对稳定，此时采样能反映出地下水在相对稳定状态下的重金属含量情况。通过在不同时期采样，可以更全面地了解地下水中重金属含量随时间的变化特征。使用专业的地下水采样设备进行样品采集。在每个采样点，首先利用便携式水位计测量地下水水位，记录水位埋深数据。然后，采用贝勒管进行水样采集，贝勒管材质为高密度聚乙烯，具有化学稳定性好、不易吸附重金属等优点。将贝勒管缓慢放入井中，到达预定深度后，通过提拉绳索使贝勒管开口关闭，确保采集到的水样具有代表性。每个采样点采集[X4]L水样，分别装入不同的采样瓶中。用于测定重金属含量的水样，采集后立即加入适量的优级纯硝酸，使水样的pH值小于2，以防止重金属离子沉淀和吸附在容器壁上。样品采集后，及时送往实验室进行分析测定。在运输过程中，将水样放置在装有冰袋的保温箱中，保持低温状态，以减少水样中化学成分的变化。同时，对每个水样进行编号和详细记录，包括采样点位置、采样时间、采样深度、水位等信息，确保样品的可追溯性。在实验室中，采用电感耦合等离子体质谱仪（ICP-MS）测定地下水中As、Cr、Cd、Pb、Hg等重金属元素的含量。ICP-MS具有灵敏度高、分析速度快、可同时测定多种元素等优点，能够准确地测定地下水中痕量重金属的含量。在测定之前，对仪器进行严格的调试和校准，使用标准溶液绘制校准曲线，确保测定结果的准确性和可靠性。同时，对每批样品进行空白试验和加标回收试验，空白试验用于检测实验室环境和试剂对测定结果的影响，加标回收试验用于评估分析方法的准确性和可靠性。加标回收率控制在80%-120%之间，以保证测定数据的质量。除了测定重金属含量外，还对水样的pH值、溶解氧、电导率等常规水质参数进行同步测定。使用pH计测定水样的pH值，pH计经过标准缓冲溶液校准，确保测量精度在±0.01范围内；采用溶解氧仪测定溶解氧含量，溶解氧仪在使用前进行零点校准和满度校准；利用电导率仪测定电导率，电导率仪根据标准溶液进行校准。这些常规水质参数的测定，有助于分析地下水的化学性质和环境条件，为深入研究重金属在地下水中的迁移转化规律提供基础数据。3.2评价方法3.2.1单因子污染指数法单因子污染指数法是一种简单直观的污染评价方法，用于判断单一重金属在地下水中的污染程度。其原理是通过将地下水中某种重金属的实测浓度与相应的评价标准进行对比，从而得出该重金属的污染状况。该方法的计算公式为：P_i=\frac{C_i}{S_i}其中，P_i为第i种重金属的污染指数；C_i为第i种重金属的实测浓度（mg/L）；S_i为第i种重金属的评价标准（mg/L），本研究采用《地下水质量标准》（GB/T14848-2017）中的Ⅲ类标准作为评价标准，该标准是地下水适用于集中式生活饮用水水源及工农业用水的水质限定标准，能够有效反映地下水中重金属对人体健康和一般工业用水的影响程度。当P_i\leq1时，表示该重金属未超标，地下水未受到该重金属的污染；当P_i>1时，则表明该重金属超标，地下水受到了该重金属的污染，且P_i值越大，污染程度越严重。例如，若某采样点地下水中Cd的实测浓度为0.001mg/L，而《地下水质量标准》（GB/T14848-2017）中Ⅲ类标准规定Cd的限值为0.005mg/L，则Cd的污染指数P_{Cd}=\frac{0.001}{0.005}=0.2\leq1，说明该采样点地下水中Cd未超标，未受到Cd的污染。单因子污染指数法能够清晰地反映出每种重金属的污染状况，为后续深入分析地下水重金属污染问题提供了基础数据和初步判断。3.2.2地累积指数法地累积指数法由德国科学家Muller于20世纪60年代晚期提出，是一种用于研究沉积物及其它物质中重金属污染程度的定量指标，该方法充分考虑了人为污染因素、环境地球化学背景值以及自然成岩作用引起背景值变动的因素，弥补了其它评价方法的不足。在本研究中，运用该方法来判断泾惠渠灌区浅层地下水中重金属的污染等级。其计算公式为：I_{geo}=\log_2\left[\frac{C_i}{K\timesB_i}\right]其中，I_{geo}为地累积指数；C_i为元素i在地下水中的实测含量（mg/L）；B_i为元素i的地球化学背景值（mg/L），本研究采用陕西省土壤中重金属元素的背景值作为参考；K为考虑各地岩石差异可能会引起背景值变动而取的系数，取值为1.5。地累积指数的计算结果按照以下评价标准划分污染等级：当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时，为极严重污染。例如，若某采样点地下水中Pb的实测含量为0.05mg/L，陕西省土壤中Pb的背景值为26mg/kg（换算为以L为单位的浓度），则地累积指数I_{geo(Pb)}=\log_2\left[\frac{0.05}{1.5\timesB_{Pb}}\right]，通过计算得到的I_{geo(Pb)}值，对照上述污染等级划分标准，即可判断该采样点地下水中Pb的污染等级。地累积指数法能够综合考虑多种因素，较为全面地反映地下水中重金属的污染程度和污染等级，有助于准确评估灌区浅层地下水的污染状况。3.2.3聚类分析和主成分分析聚类分析和主成分分析是多元统计分析中的重要方法，在本研究中用于解析泾惠渠灌区浅层地下水中重金属的污染来源。聚类分析是根据事物本身的特性研究个体分类的方法，它能够将具有相似性质或特征的数据对象归为一类。在本研究中，通过对地下水中不同重金属元素含量之间的相似性进行分析，将重金属元素划分为不同的类别。具体步骤为：首先，计算各重金属元素之间的距离或相似性度量，常用的距离度量方法有欧氏距离、曼哈顿距离等，相似性度量方法有相关系数等；然后，根据距离或相似性度量结果，采用系统聚类法、K-均值聚类法等聚类算法进行聚类分析。例如，若通过计算发现As、Cd两种重金属元素的含量在不同采样点之间的变化趋势较为相似，它们之间的距离或相似性度量值满足一定的聚类条件，那么在聚类分析结果中，As和Cd可能会被归为同一类，这表明它们可能具有相似的来源或受到相似的环境因素影响。聚类分析能够初步将具有相似来源或行为的重金属元素聚集在一起，为进一步分析污染来源提供线索。主成分分析是把原来多个变量化为少数几个综合指标的一种统计分析方法，从数学角度来看，这是一种降维处理技术。在本研究中，主成分分析用于提取地下水中重金属数据的主要特征，找出影响重金属含量变化的主要因素，从而解析污染来源。其原理是通过线性变换，将原来的多个相关变量转换为少数几个互不相关的综合变量，即主成分，这些主成分能够尽可能多地保留原来变量的信息。具体计算步骤如下：计算相关系数矩阵：首先对地下水中各重金属元素的含量数据进行标准化处理，消除量纲和数量级的影响。然后计算各重金属元素之间的相关系数，得到相关系数矩阵R，其中r_{ij}（i,j=1,2,\cdots,p）为第i种重金属与第j种重金属的相关系数，计算公式为r_{ij}=\frac{\sum_{k=1}^{n}(x_{ki}-\overline{x_i})(x_{kj}-\overline{x_j})}{\sqrt{\sum_{k=1}^{n}(x_{ki}-\overline{x_i})^2\sum_{k=1}^{n}(x_{kj}-\overline{x_j})^2}}，x_{ki}为第k个采样点中第i种重金属的含量，\overline{x_i}为第i种重金属含量的平均值，n为采样点个数，p为重金属元素个数。计算特征值与特征向量：解特征方程|\lambdaI-R|=0，求出特征值\lambda_i（i=1,2,\cdots,p），并使其按大小顺序排列，即\lambda_1\geq\lambda_2\geq\cdots\geq\lambda_p\geq0；然后分别求出对应于特征值\lambda_i的特征向量e_i（i=1,2,\cdots,p）。计算主成分贡献率及累计贡献率：主成分贡献率\alpha_i=\frac{\lambda_i}{\sum_{j=1}^{p}\lambda_j}，累计贡献率\beta_m=\sum_{i=1}^{m}\alpha_i（m\leqp）。一般取累计贡献率达85%-95%的特征值\lambda_1,\lambda_2,\cdots,\lambda_m所对应的第一，第二，……，第m（m\leqp）个主成分。计算主成分载荷：主成分载荷l_{ij}=\sqrt{\lambda_i}e_{ij}，其中l_{ij}为第i个主成分在第j个变量上的载荷，e_{ij}为第i个特征向量的第j个分量。通过分析主成分载荷，可以确定每个主成分与原始重金属变量之间的关系，从而识别出影响重金属含量的主要因素。例如，若经过主成分分析得到第一主成分，其贡献率较高，且在该主成分中，Cr、Pb等重金属元素的载荷较大，这说明Cr、Pb等重金属元素在第一主成分中起主要作用，它们可能受到同一来源的影响，如工业排放；而第二主成分中，As、Cd等重金属元素的载荷较大，表明As、Cd等重金属元素可能受到另一种来源的影响，如农业活动。主成分分析能够有效地从复杂的重金属数据中提取主要信息，识别出潜在的污染来源，为针对性地制定污染防治措施提供科学依据。3.2.4健康风险评价模型本研究采用美国环保署（USEPA）推荐的健康风险评价模型，对泾惠渠灌区浅层地下水中重金属对人体健康的潜在风险进行评估。该模型将健康风险分为致癌风险和非致癌风险，充分考虑了不同重金属元素的毒性、人体暴露途径以及暴露参数等因素，能够较为全面地评估重金属对人体健康的影响。对于致癌风险，通常认为人体在低剂量暴露条件下，暴露剂量率和人体致癌风险之间呈线性关系。致癌风险值（CR）的计算公式如下：CR=\sum_{i=1}^{n}\left(\frac{C_i\timesIR_i\timesEF\timesED}{BW\timesAT}\timesSF_i\right)其中，C_i为第i种致癌重金属在地下水中的浓度（mg/L）；IR_i为第i种重金属的摄入速率（L/d），对于饮水暴露途径，成人的饮水摄入速率一般取值为2L/d，儿童的饮水摄入速率一般取值为0.7L/d；EF为暴露频率（d/年），考虑到灌区居民全年都可能接触地下水，取值为365d/年；ED为暴露持续时间（年），成人的暴露持续时间一般取值为70年，儿童的暴露持续时间一般取值为6年；BW为体重（kg），成人的平均体重一般取值为70kg，儿童的平均体重一般取值为15kg；AT为平均暴露时间（d），对于致癌风险，AT=ED\times365；SF_i为第i种致癌重金属的致癌斜率因子（mg/(kg・d)），不同重金属的致癌斜率因子可从USEPA的相关数据库中获取，如As的致癌斜率因子为1.5（mg/(kg・d)）。对于非致癌风险，以参考剂量（RfD）值为衡量标准，暴露剂量率和参考剂量的关系用危害商（HQ）来表示。危害商的计算公式如下：HQ=\sum_{i=1}^{n}\frac{C_i\timesIR_i\timesEF\timesED}{BW\timesAT\timesRfD_i}其中，RfD_i为第i种非致癌重金属的参考剂量（mg/(kg・d)），同样可从USEPA的相关数据库中获取，如Cr的参考剂量为0.003mg/(kg・d)，其余参数含义与致癌风险计算公式中相同。当HQ\leq1时，表示非致癌风险处于可接受水平；当HQ>1时，则表明存在非致癌风险，且HQ值越大，风险越高。对于致癌风险，通常设定10^{-6}为可接受致癌风险水平下限，10^{-4}为可接受致癌风险水平上限。若计算得到的致癌风险值在10^{-6}-10^{-4}之间，表明致癌风险处于可接受范围；若小于10^{-6}，则致癌风险较低；若大于10^{-4}，则致癌风险较高。通过该健康风险评价模型，能够准确计算出泾惠渠灌区浅层地下水中重金属对不同人群（成人、儿童）通过饮水暴露途径的致癌风险和非致癌风险，为评估当地居民的健康风险水平提供科学依据。四、泾惠渠灌区浅层地下水重金属含量与分布特征4.1重金属含量分析通过对泾惠渠灌区浅层地下水样品的分析测定，得到了各采样点地下水中Cr、Cd、Pb、Ni、Cu、Fe、Mn、As等重金属的含量数据，具体结果如表1所示。表1泾惠渠灌区浅层地下水重金属含量（mg/L）采样点编号CrCdPbNiCuFeMnAs1[具体数值1][具体数值2][具体数值3][具体数值4][具体数值5][具体数值6][具体数值7][具体数值8]2[具体数值9][具体数值10][具体数值11][具体数值12][具体数值13][具体数值14][具体数值15][具体数值16]3[具体数值17][具体数值18][具体数值19][具体数值20][具体数值21][具体数值22][具体数值23][具体数值24]...........................[X][具体数值X1][具体数值X2][具体数值X3][具体数值X4][具体数值X5][具体数值X6][具体数值X7][具体数值X8]将这些重金属含量数据与《地下水质量标准》（GB/T14848-2017）中的Ⅲ类标准进行对比，判断其超标情况。该标准中，Cr的限值为0.05mg/L，Cd的限值为0.005mg/L，Pb的限值为0.01mg/L，Ni的限值为0.02mg/L，Cu的限值为1.0mg/L，Fe的限值为0.3mg/L，Mn的限值为0.1mg/L，As的限值为0.01mg/L。统计结果显示，在所有采样点中，Cr含量的平均值为[具体平均数值1]mg/L，最大值为[具体最大数值1]mg/L，最小值为[具体最小数值1]mg/L，有[X5]个采样点的Cr含量超过了Ⅲ类标准，超标率为[X5]/[X]×100%=[具体百分比1]%。Cd含量的平均值为[具体平均数值2]mg/L，最大值为[具体最大数值2]mg/L，最小值为[具体最小数值2]mg/L，有[X6]个采样点的Cd含量超标，超标率为[X6]/[X]×100%=[具体百分比2]%。Pb含量的平均值为[具体平均数值3]mg/L，最大值为[具体最大数值3]mg/L，最小值为[具体最小数值3]mg/L，超标采样点有[X7]个，超标率为[X7]/[X]×100%=[具体百分比3]%。Ni含量的平均值为[具体平均数值4]mg/L，最大值为[具体最大数值4]mg/L，最小值为[具体最小数值4]mg/L，超标率为[X8]/[X]×100%=[具体百分比4]%，超标采样点数为[X8]。Cu含量的平均值为[具体平均数值5]mg/L，最大值为[具体最大数值5]mg/L，最小值为[具体最小数值5]mg/L，所有采样点的Cu含量均未超过标准限值。Fe含量的平均值为[具体平均数值6]mg/L，最大值为[具体最大数值6]mg/L，最小值为[具体最小数值6]mg/L，超标采样点数量为[X9]，超标率为[X9]/[X]×100%=[具体百分比5]%。Mn含量的平均值为[具体平均数值7]mg/L，最大值为[具体最大数值7]mg/L，最小值为[具体最小数值7]mg/L，有[X10]个采样点的Mn含量超标，超标率为[X10]/[X]×100%=[具体百分比6]%。As含量的平均值为[具体平均数值8]mg/L，最大值为[具体最大数值8]mg/L，最小值为[具体最小数值8]mg/L，超标采样点数为[X11]，超标率为[X11]/[X]×100%=[具体百分比7]%。从上述数据可以看出，泾惠渠灌区浅层地下水中Cr、Cd、Pb、Ni、Fe、Mn、As等重金属均存在不同程度的超标现象，其中[超标率最高的重金属元素]的超标情况较为突出，这表明灌区浅层地下水已受到一定程度的重金属污染，需要引起高度重视。4.2空间分布特征运用地统计学和GIS技术，对泾惠渠灌区浅层地下水中重金属含量的空间分布特征进行深入分析。地统计学中的变异函数能够定量描述区域化变量的空间变异性和相关性，通过计算变异函数，可以了解重金属含量在空间上的变化规律。在本研究中，首先对各重金属含量数据进行半变异函数分析，以确定其空间自相关范围和结构特征。以Cr元素为例，计算得到的半变异函数云图及拟合模型参数如下（表2）：表2Cr元素半变异函数拟合参数模型块金值（C0）基台值（C0+C）变程（A0）决定系数（R²）[具体模型名称][具体数值1][具体数值2][具体数值3][具体数值4]从半变异函数拟合结果可以看出，Cr元素的块金值与基台值之比（C0/（C0+C））为[具体比例值]，该比值反映了随机因素对空间变异性的影响程度。当C0/（C0+C）<25%时，表示空间相关性强，即Cr含量在空间上的变化主要受结构性因素控制，如地质背景、土壤类型等；当25%≤C0/（C0+C）≤75%时，空间相关性中等，说明随机因素和结构性因素对Cr含量的空间分布都有一定影响；当C0/（C0+C）>75%时，空间相关性弱，随机因素起主导作用。在本研究中，若Cr元素的C0/（C0+C）处于中等范围，这表明Cr含量的空间分布既受到地质等结构性因素的影响，也受到人类活动等随机因素的干扰。变程（A0）表示在该距离范围内，Cr含量具有空间自相关性，超出此距离，空间自相关性消失。通过半变异函数分析得到的这些参数，为后续的空间插值和分布特征分析提供了重要依据。利用克里金插值法，基于半变异函数模型对各重金属含量进行空间插值，得到泾惠渠灌区浅层地下水中重金属含量的空间分布图（图1）。[此处插入Cr、Cd、Pb等重金属含量空间分布图]从Cr含量空间分布图（图1a）可以看出，在灌区的[具体高值区域位置]，如靠近某工业集中区附近，Cr含量较高，呈现出明显的高值聚集区。这可能是由于该区域工业活动频繁，工业废水排放中含有大量的Cr元素，通过地表径流、土壤淋溶等途径进入地下水，导致该区域地下水中Cr含量升高。而在灌区的[具体低值区域位置]，如远离工业源的农业种植区，Cr含量相对较低。这说明工业活动是影响该区域地下水中Cr含量空间分布的重要因素之一。Cd含量空间分布图（图1b）显示，高值区主要集中在[具体高值区域]，这些区域可能受到农业面源污染和部分小型乡镇企业污染的影响。农业生产中大量使用的磷肥、农药以及农用塑料薄膜中可能含有Cd元素，随着时间的推移，这些Cd元素逐渐进入土壤，并通过淋溶作用进入地下水。同时，部分小型乡镇企业在生产过程中缺乏有效的污染治理措施，其排放的废水和废渣中也可能含有Cd，进一步加重了周边地下水的污染。低值区则主要分布在[具体低值区域]，这些区域生态环境相对较好，人类活动强度较小，对地下水的污染相对较轻。Pb含量空间分布图（图1c）呈现出在[具体高值区域位置]，如交通干线附近和城镇周边，Pb含量较高的特征。交通干线附近的高Pb含量可能与汽车尾气排放、道路扬尘等有关，汽车尾气中含有一定量的Pb，随着大气沉降进入土壤和水体，进而影响地下水。城镇周边人口密集，生活污水排放、垃圾填埋等人类活动也可能导致Pb进入地下水。而在[具体低值区域位置]，如山区或生态保护区，Pb含量较低，这些区域受人类活动干扰较小，自然生态系统相对完整，对地下水起到了较好的保护作用。通过对各重金属含量空间分布图的分析，可以发现泾惠渠灌区浅层地下水中重金属含量的空间分布呈现出明显的区域性差异，这种差异与地形地貌、土地利用类型、工业布局、农业活动以及人口分布等因素密切相关。在地势低洼、排水不畅的区域，由于污染物容易聚集，地下水中重金属含量相对较高；在工业集中区和城镇周边，人类活动强度大，污染物排放多，重金属污染较为严重；而在农业种植区和生态保护区，人类活动相对较少，地下水中重金属含量相对较低。这些空间分布特征的揭示，有助于准确识别重金属污染的高风险区域，为针对性地制定污染防治措施提供科学依据。五、泾惠渠灌区浅层地下水重金属污染评价5.1单因子污染指数评价根据单因子污染指数法的计算公式，对泾惠渠灌区浅层地下水中的Cr、Cd、Pb、Ni、Fe、Mn、As等重金属进行单因子污染指数计算，计算结果如表3所示。表3泾惠渠灌区浅层地下水重金属单因子污染指数采样点编号Cr污染指数Cd污染指数Pb污染指数Ni污染指数Fe污染指数Mn污染指数As污染指数1[具体数值11][具体数值12][具体数值13][具体数值14][具体数值15][具体数值16][具体数值17]2[具体数值21][具体数值22][具体数值23][具体数值24][具体数值25][具体数值26][具体数值27]3[具体数值31][具体数值32][具体数值33][具体数值34][具体数值35][具体数值36][具体数值37]........................[X][具体数值X11][具体数值X12][具体数值X13][具体数值X14][具体数值X15][具体数值X16][具体数值X17]通过对各采样点重金属单因子污染指数的分析，可以判断出不同重金属在不同区域的污染程度。当污染指数P_i\leq1时，表明该重金属未超标，地下水未受到该重金属的污染；当P_i>1时，则说明该重金属超标，地下水受到了该重金属的污染，且P_i值越大，污染程度越严重。统计各重金属污染指数的平均值和最大值，结果显示，Cr污染指数的平均值为[具体平均数值11]，最大值为[具体最大数值11]，在[具体采样点编号]采样点处出现最大值，该点的Cr污染指数远大于1，表明该区域地下水中Cr污染较为严重，可能受到了周边工业排放或其他污染源的强烈影响。Cd污染指数的平均值为[具体平均数值12]，最大值为[具体最大数值12]，在[具体采样点编号]采样点达到最大值，该点的Cd污染指数超出1较多，说明该区域存在较严重的Cd污染，可能与农业面源污染或小型工业污染有关。Pb污染指数的平均值为[具体平均数值13]，最大值为[具体最大数值13]，在[具体采样点编号]采样点处最大值超出1一定范围，表明该区域的Pb污染不容忽视，可能与交通污染、工业废弃物排放等因素有关。Ni污染指数的平均值为[具体平均数值14]，最大值为[具体最大数值14]，在[具体采样点编号]采样点出现较大值，虽然整体超标情况相对较轻，但个别区域仍存在一定程度的Ni污染，可能受到土壤母质和工业活动的共同影响。Fe污染指数的平均值为[具体平均数值15]，最大值为[具体最大数值15]，在[具体采样点编号]采样点最大值远超1，说明该区域Fe污染严重，可能是由于地下水与富含铁的岩石或土壤发生相互作用，或者受到工业废水排放的影响。Mn污染指数的平均值为[具体平均数值16]，最大值为[具体最大数值16]，在[具体采样点编号]采样点处出现高值，该区域的Mn污染可能与地质条件和农业活动有关。As污染指数的平均值为[具体平均数值17]，最大值为[具体最大数值17]，在[具体采样点编号]采样点最大值显著超出1，表明该区域As污染严重，As作为一种致癌重金属，其污染对人体健康具有较大威胁，可能与矿业活动、含砷农药使用等因素有关。综合各重金属的单因子污染指数情况，[污染指数最高的重金属元素]的污染程度相对较高，是泾惠渠灌区浅层地下水的主要污染因子之一。通过对各采样点污染指数的分析，可以明确[具体污染区域位置]，如[具体污染严重的工业集中区、农业种植区等]，这些区域是重金属污染的高风险区域，需要重点关注和治理。同时，不同区域的污染因子存在差异，在制定污染防治措施时，应根据各区域的污染特征，有针对性地采取措施，以有效降低地下水重金属污染程度，保障地下水水质安全。5.2地累积指数评价根据第三章中给出的地累积指数计算公式，计算泾惠渠灌区浅层地下水中各重金属的地累积指数，结果如表4所示。表4泾惠渠灌区浅层地下水重金属地累积指数采样点编号Cr地累积指数Cd地累积指数Pb地累积指数Ni地累积指数Fe地累积指数Mn地累积指数As地累积指数1[具体数值21][具体数值22][具体数值23][具体数值24][具体数值25][具体数值26][具体数值27]2[具体数值31][具体数值32][具体数值33][具体数值34][具体数值35][具体数值36][具体数值37]3[具体数值41][具体数值42][具体数值43][具体数值44][具体数值45][具体数值46][具体数值47]........................[X][具体数值X21][具体数值X22][具体数值X23][具体数值X24][具体数值X25][具体数值X26][具体数值X27]按照地累积指数的评价标准，对各重金属的污染等级进行划分。地累积指数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时，为极严重污染。统计各重金属不同污染等级的采样点数量及比例，结果如表5所示。表5泾惠渠灌区浅层地下水重金属污染等级统计重金属无污染（%）轻度-中等污染（%）中等污染（%）中等-强污染（%）强污染（%）强-极严重污染（%）极严重污染（%）Cr[具体比例1][具体比例2][具体比例3][具体比例4][具体比例5][具体比例6][具体比例7]Cd[具体比例8][具体比例9][具体比例10][具体比例11][具体比例12][具体比例13][具体比例14]Pb[具体比例15][具体比例16][具体比例17][具体比例18][具体比例19][具体比例20][具体比例21]Ni[具体比例22][具体比例23][具体比例24][具体比例25][具体比例26][具体比例27][具体比例28]Fe[具体比例29][具体比例30][具体比例31][具体比例32][具体比例33][具体比例34][具体比例35]Mn[具体比例36][具体比例37][具体比例38][具体比例39][具体比例40][具体比例41][具体比例42]As[具体比例43][具体比例44][具体比例45][具体比例46][具体比例47][具体比例48][具体比例49]从表5可以看出，泾惠渠灌区浅层地下水中不同重金属的污染等级存在差异。其中，Cr的污染等级以轻度-中等污染和中等污染为主，在[具体采样点编号]等采样点出现了中等-强污染的情况，表明这些区域的Cr污染相对较为严重，可能与周边的工业活动或其他污染源有关。Cd的污染等级分布较为分散，无污染、轻度-中等污染和中等污染的采样点均占有一定比例，但也存在个别采样点达到中等-强污染水平，这可能与农业面源污染和部分小型工业排放有关。Pb的污染主要集中在无污染和轻度-中等污染等级，但仍有部分采样点存在中等污染情况，可能受到交通污染、工业废弃物排放等因素的影响。Ni的污染程度相对较轻，大部分采样点处于无污染状态，仅有少数采样点出现轻度-中等污染，这说明Ni在该区域地下水中的污染情况相对不突出，可能与土壤母质和工业活动的影响相对较小有关。Fe的污染等级以中等污染和中等-强污染为主，在[具体采样点编号]等采样点污染较为严重，可能是由于地下水与富含铁的岩石或土壤发生相互作用，或者受到工业废水排放的影响。Mn的污染主要表现为轻度-中等污染和中等污染，在[具体采样点编号]等采样点存在一定程度的污染，可能与地质条件和农业活动有关。As的污染较为严重，部分采样点达到强污染和强-极严重污染等级，As作为一种致癌重金属，其高污染等级对人体健康具有较大威胁，可能与矿业活动、含砷农药使用等因素有关。综合各重金属的地累积指数评价结果，[污染等级最高的重金属元素]的污染等级相对较高，是泾惠渠灌区浅层地下水的主要污染因子之一。通过地累积指数评价，能够更全面地了解各重金属在地下水中的污染程度和潜在生态风险，为制定合理的污染防治措施提供科学依据。在[具体高污染区域位置]，如[具体污染严重的工业集中区、农业种植区等]，应加强对重金属污染的监测和治理，采取有效的措施减少重金属的排放和迁移，降低其对地下水环境和人体健康的危害。5.3污染来源解析为深入探究泾惠渠灌区浅层地下水中重金属的污染来源，本研究运用聚类分析和主成分分析方法，对地下水中重金属含量数据及相关环境因素数据进行处理和分析。聚类分析结果（表6）显示，泾惠渠灌区浅层地下水中的重金属可分为3类。第1类包括Cr和Pb，这两种重金属在工业生产中广泛应用，例如在电镀、金属加工、化工等行业，Cr常用于镀铬工艺，Pb则用于电池制造、涂料生产等。在这些工业活动中，若废水未经有效处理直接排放，会导致Cr和Pb进入地表水和土壤，进而通过渗漏等方式污染地下水。第2类包含Cd和As，农业活动是这两种重金属的主要来源之一。Cd常存在于磷肥中，长期大量使用磷肥会使土壤中Cd含量增加，通过淋溶作用进入地下水。As则可能来源于含砷农药的使用，以及一些含砷的农业废弃物的不合理处置。此外，部分矿业活动也可能导致As和Cd的释放，污染周边的地下水。第3类为Fe和Mn，它们主要来源于自然地质过程。该地区的岩石和土壤中富含Fe和Mn元素，在地下水的长期溶滤作用下，这些元素会逐渐溶解进入地下水中。同时，一些人为活动，如建筑施工、道路建设等，可能破坏地表植被和土壤结构，加速岩石和土壤的风化，从而增加Fe和Mn向地下水中的释放量。表6聚类分析结果类别重金属元素可能来源1Cr、Pb工业排放（电镀、金属加工、化工、电池制造、涂料生产等行业废水排放）2Cd、As农业活动（磷肥使用、含砷农药使用、农业废弃物处置）、矿业活动3Fe、Mn自然地质过程（岩石和土壤溶滤）、人为活动（建筑施工、道路建设等导致的岩石和土壤风化加速）主成分分析结果表明，前3个主成分的累计贡献率达到[具体累计贡献率数值]%，能够较好地解释原始数据的大部分信息。第1主成分贡献率为[具体贡献率数值1]%，在该主成分中，Cr、Pb、Zn等重金属元素的载荷较大，与工业污染源密切相关。如前所述，工业生产过程中产生的废水、废气和废渣中含有大量的重金属，这些重金属通过各种途径进入环境，最终污染地下水。在泾惠渠灌区，部分工业企业集中分布在特定区域，其排放的污染物可能在周边地下水中积累，导致这些区域地下水中Cr、Pb等重金属含量升高。通过对工业企业的调查发现，一些电镀厂、金属加工厂的废水处理设施不完善，存在废水超标排放的情况，这可能是导致第1主成分中相关重金属载荷较高的原因之一。第2主成分贡献率为[具体贡献率数值2]%，As、Cd等重金属元素在该主成分上有较高载荷，主要与农业活动和矿业活动有关。农业活动中，除了上述提到的磷肥和农药使用外，畜禽养殖废弃物的排放也可能含有一定量的As和Cd。畜禽饲料中常添加一些含重金属的添加剂，畜禽摄入后，部分重金属会通过粪便排出，若这些粪便未经妥善处理直接用于农田施肥，会增加土壤和地下水中As和Cd的含量。矿业活动方面，灌区周边可能存在一些小型矿山，在开采、选矿等过程中，会产生大量的尾矿和废渣，其中含有丰富的As、Cd等重金属。这些尾矿和废渣若随意堆放，在雨水淋溶作用下，重金属会随淋滤液进入地下水，造成污染。第3主成分贡献率为[具体贡献率数值3]%，Fe、Mn等重金属元素在该主成分上的载荷较大，主要受自然地质因素的影响。从地质构造上看，泾惠渠灌区位于特定的地质区域，其地下岩石和土壤中Fe、Mn的含量相对较高。地下水在与这些岩石和土壤相互作用过程中，会溶解其中的Fe、Mn元素，使其进入地下水中。此外，地下水的酸碱度、氧化还原条件等也会影响Fe、Mn的溶解和迁移。在一些偏酸性的地下水中，Fe、Mn的溶解度会增加，从而导致地下水中Fe、Mn含量升高。通过计算各主成分得分与重金属含量之间的相关性系数，进一步验证了上述结论。各主成分得分与对应重金属含量之间的相关性系数较高，表明主成分分析结果能够较好地反映重金属的污染来源。根据主成分分析结果，计算各污染源对地下水中重金属污染的贡献比例。结果显示，工业污染源对地下水中重金属污染的贡献比例为[具体贡献比例1]%，农业和矿业污染源的贡献比例为[具体贡献比例2]%，自然地质因素的贡献比例为[具体贡献比例3]%。综上所述，泾惠渠灌区浅层地下水中重金属的污染来源主要包括工业排放、农业活动、矿业活动以及自然地质过程，其中工业污染源和农业及矿业污染源对重金属污染的贡献相对较大。针对不同的污染来源，应采取相应的污染防治措施，如加强工业企业的环境监管，提高废水处理水平；推广绿色农业生产方式，合理使用化肥和农药；加强矿业活动的管理，规范尾矿和废渣的处置等，以有效减少地下水中重金属的污染。六、泾惠渠灌区浅层地下水重金属健康风险评价6.1健康风险评价模型本研究采用美国环保署（USEPA）推荐的健康风险评价模型，该模型是目前国际上广泛应用且认可度较高的评估工具，其原理基于对人体暴露于污染物的剂量与健康效应之间关系的定量分析。通过对污染物在环境中的浓度、人体暴露途径、暴露时间以及污染物的毒性等多方面因素的综合考量，准确评估重金属对人体健康产生的潜在风险。在该模型中，人体暴露途径主要考虑饮水摄入和皮肤接触两种方式。饮水摄入是人体接触地下水中重金属的重要途径之一，人们日常饮用的地下水若受到重金属污染，重金属会随着饮水直接进入人体消化系统，进而被人体吸收。皮肤接触则是指人体在日常生活中，如洗漱、沐浴等过程中，皮肤与含有重金属的地下水直接接触，重金属通过皮肤渗透进入人体。这两种暴露途径在灌区居民的日常生活中较为常见，且对健康风险的贡献较大，因此在本研究中重点考虑这两种途径。模型中涉及多个关键参数，每个参数都具有特定的含义和取值依据。如日均暴露剂量（ADD），它表示人体每天通过各种暴露途径接触到的重金属的平均剂量，计算公式为：ADD=\frac{C\timesIR\timesEF\timesED}{BW\timesAT}其中，C为地下水中重金属的浓度（mg/L），这是通过对泾惠渠灌区浅层地下水样品的实际检测分析得到的，能真实反映灌区地下水中重金属的含量情况；IR为摄入速率（L/d），对于饮水暴露途径，成人的饮水摄入速率一般取值为2L/d，儿童的饮水摄入速率一般取值为0.7L/d，这是基于相关的人体生理研究和实际生活调查得出的，不同年龄段人群的饮水量存在差异，儿童由于身体发育尚未完全，饮水量相对较少；EF为暴露频率（d/年），考虑到灌区居民全年都可能接触地下水，取值为365d/年，这是符合实际生活情况的设定；ED为暴露持续时间（年），成人的暴露持续时间一般取值为70年，儿童的暴露持续时间一般取值为6年，这是根据不同年龄段人群的预期寿命和实际生活阶段确定的，成人的寿命较长，暴露时间相应较长，而儿童在成长过程中，暴露于地下水的时间相对较短；BW为体重（kg），成人的平均体重一般取值为70kg，儿童的平均体重一般取值为15kg，这是基于大量的人体生理数据统计得出的，不同年龄段人群的体重差异明显，体重因素会影响人体对重金属的吸收和代谢；AT为平均暴露时间（d），对于致癌风险，AT=ED\times365，对于非致癌风险，AT取值为暴露持续时间内的总天数，这样的取值能够准确反映不同风险类型下人体的平均暴露时间。参考剂量（RfD）是指人类长期暴露于某污染物而不产生可检测到的有害健康效应的日平均暴露剂量估计值，单位为mg/(kg・d)。不同重金属的参考剂量不同，例如Cr的参考剂量为0.003mg/(kg・d)，这是通过大量的毒理学实验和研究，综合考虑重金属的毒性、生物可利用性以及人体对其的耐受性等多方面因素确定的。致癌斜率因子（SF）是指在低剂量范围内，单位剂量的化学物质暴露引起的终身超额致癌风险的统计估计值，单位为mg/(kg・d)。如As的致癌斜率因子为1.5（mg/(kg・d)），它是基于大量的流行病学研究和动物实验数据，通过科学的统计分析和评估得出的，用于衡量As对人体致癌风险的潜在影响程度。通过这些参数的准确取值和模型的合理运用，能够全面、准确地评估泾惠渠灌区浅层地下水中重金属对人体健康的潜在风险，为后续制定科学合理的污染防治措施和健康风险管理策略提供坚实的科学依据。6.2暴露参数确定暴露参数是健康风险评价中的关键要素，其准确性直接影响到风险评价结果的可靠性。在本研究中，针对不同暴露途径下成人和儿童的暴露参数进行了详细确定，具体取值方法和依据如下：饮水摄入途径：摄入速率（IR）：成人的饮水摄入速率取值为2L/d，这是基于大量的人体生理研究和实际生活调查得出的。相关研究表明，成人在正常生活状态下，每天通过饮水摄入的水量平均约为2L。儿童的饮水摄入速率一般取值为0.7L/d，儿童由于身体发育尚未完全，新陈代谢相对较快，但身体体积和需水量相对较小，因此饮水量明显低于成人。暴露频率（EF）：考虑到泾惠渠灌区居民全年都可能接触地下水用于日常生活饮用，所以暴露频率取值为365d/年，以准确反映居民全年通过饮水接触地下水中重金属的实际情况。暴露持续时间（ED）：成人的暴露持续时间一般取值为70年，这是根据人类平均预期寿命以及居民长期生活在该地区的实际情况确定的，成人在其一生中较长时间都可能暴露于当地的地下水环境中。儿童的暴露持续时间一般取值为6年，这是考虑到儿童在成长过程中，在幼儿和童年阶段主要接触当地的地下水，随着年龄增长，可能会因生活环境变化等因素，接触当地地下水的时间相对减少。体重（BW）：成人的平均体重一般取值为70kg，这是基于大量的人体生理数据统计得出的，不同地区成人的平均体重虽有一定差异，但70kg是一个被广泛接受的代表值。儿童的平均体重一般取值为15kg，儿童在成长过程中体重不断变化，15kg是根据儿童在相应年龄段的平均体重统计结果确定的，能够较好地反映儿童在该阶段的身体特征。皮肤接触途径：皮肤接触面积（SA）：成人的皮肤接触面积取值为1.73m²，这是根据人体表面积的计算公式和大量人体测量数据统计得出的。人体表面积与身高、体重等因素相关，通过相关公式计算并结合实际统计，1.73m²是成人皮肤接触面积的常见取值。儿童的皮肤接触面积取值为0.34m²，儿童身体较小，皮肤表面积也相应较小，根据儿童身体发育特点和相关测量数据，确定此值以准确反映儿童皮肤接触的实际情况。皮肤渗透系数（Kp）：不同重金属的皮肤渗透系数不同，如As的皮肤渗透系数为0.000014cm/h，Cr的皮肤渗透系数为0.000002cm/h，这些数值是通过大量的实验室研究和毒理学实验获得的，反映了不同重金属通过皮肤渗透进入人体的难易程度。暴露时间（ET）：对于成人和儿童，在计算皮肤接触暴露参数时，假设每次接触地下水的时间为0.5h，这是考虑到居民在日常生活中，如洗漱、沐浴等与地下水接触的活动，每次持续时间一般在0.5h左右。暴露频率同样取值为365d/年，以反映全年的皮肤接触情况。通过对不同暴露途径下成人和儿童暴露参数的准确确定，能够更真实地反映泾惠渠灌区居民通过饮水和皮肤接触途径暴露于浅层地下水中重金属的实际情况，为后续准确计算健康风险提供了可靠的基础数据。6.3健康风险计算结果基于前文确定的健康风险评价模型和暴露参数，分别计算泾惠渠灌区浅层地下水中重金属通过饮水和皮肤接触途径对成人和儿童产生的致癌风险和非致癌风险值，具体计算结果如表7所示。表7泾惠渠灌区浅层地下水重金属健康风险计算结果重金属暴露途径成人儿童致癌风险非致癌风险致癌风险非致癌风险As饮水[具体数值1][具体数值2][具体数值3][具体数值4]皮肤接触[具体数值5][具体数值6][具体数值7][具体数值8]Cr饮水[具体数值9][具体数值10][具体数值11][具体数值12]皮肤接触[具体数值13][具体数值14][具体数值15][具体数值16]Cd饮水[具体数值17][具体数值18][具体数值19][具体数值20]皮肤接触[具体数值21][具体数值22][具体数值23][具体数值24]..................从计算结果可以看出，对于致癌风险，As通过饮水途径对成人和儿童产生的致癌风险值相对较高，分别为[具体数值1]和[具体数值3]，均超过了可接受致癌风险水平上限10^{-4}，表明As通过饮水途径对人体致癌风险较高，需要引起高度关注。这可能是由于泾惠渠灌区部分区域地下水中As含量超标严重，且As是一种明确的致癌物质，长期通过饮水摄入会在人体内逐渐积累，增加患癌风险。Cr和Cd通过饮水途径对成人和儿童的致癌风险值相对较低，均在可接受范围内，但也不容忽视，仍需持续关注其对人体健康的潜在影响。通过皮肤接触途径，As、Cr和Cd对成人和儿童的致癌风险值均较低，远低于可接受致癌风险水平下限10^{-6}，说明皮肤接触途径对人体的致癌风险相对较小。在非致癌风险方面，As通过饮水途径对成人和儿童的非致癌风险值分别为[具体数值2]和[具体数值4]，均大于1，表明存在非致癌风险，且风险相对较高。这可能是由于As的毒性较大，即使在较低浓度下长期摄入也可能对人体健康产生不良影响。Cr通过饮水途径对成人和儿童的非致癌风险值分别为[具体数值10]和[具体数值12]，虽小于1，但已接近1，说明存在一定的潜在风险。Cd通过饮水途径对成人和儿童的非致癌风险值相对较小，均远小于1，表明非致癌风险较低。通过皮肤接触途径，As、Cr和Cd对成人和儿童的非致癌风险值均较小，远小于1，说明皮肤接触途径对人体的非致癌风险相对较低。总体而言，泾惠渠灌区浅层地下水中重金属对人体健康存在一定风险，其中As通过饮水途径产生的致癌风险和非致癌风险较为突出，是需要重点关注和防控的关键重金属元素。在制定污染防治措施和保障居民健康的策略时，应优先考虑降低As的污染风险，加强对地下水中As含量的监测和治理，同时关注其他重金属元素的潜在风险，采取综合措施减少重金属对人体健康的危害。6.4风险评价结果分析将计算得到的泾惠渠灌区浅层地下水中重金属的健康风险值与国际上普遍认可的可接受风险水平进行对比，以便更直观地评估风险程度。国际上通常设定10^{-6}为可接受致癌风险水平下限，10^{-4}为可接受致癌风险水平上限；对于非致癌风险，当危害商HQ\leq1时，表示非致癌风险处于可接受水平。通过对比发现，不同重金属对人体健康的风险存在显著差异。As的致癌风险和非致癌风险在所有重金属中表现最为突出。在致癌风险方面，As通过饮水途径对成人和儿童产生的致癌风险值均超过了可接受致癌风险水平上限10^{-4}，分别达到[具体数值1]和[具体数值3]，这表明长期饮用受As污染的地下水，居民患癌的风险较高。从非致癌风险来看，As通过饮水途径对成人和儿童的非致癌风险值分别为[具体数值2]