基于15N同位素解析河北平原石家庄地区地下水硝酸盐污染特征、来源与环境响应

基于15N同位素解析河北平原石家庄地区地下水硝酸盐污染特征、来源与环境响应一、引言1.1研究背景与意义随着经济社会的快速发展，水资源的供需矛盾日益突出，已成为全球性的问题。我国作为人口众多的发展中国家，水资源短缺问题尤为严峻。地下水作为重要的淡水资源之一，不仅是我国农村居民的主要饮用水源，还在工业生产、农业灌溉等领域发挥着不可或缺的作用。然而，在经济发展的过程中，地下水面临着严峻的环境问题，其中硝酸盐污染问题尤为突出，对人体健康和生态环境构成了严重威胁。河北省作为我国的经济大省，水资源供需矛盾十分尖锐，地下水超采和污染问题严重。石家庄地区作为河北平原的中心城市，其饮用水主要依赖于地下水，地下水质量问题显得尤为重要。研究表明，石家庄地区地下水中硝酸盐含量呈上升趋势，部分地区已超过国家饮用水标准，严重影响了当地居民的饮水安全和身体健康。如《石家庄市地下水中硝酸盐的污染状况》一文指出，石家庄位于滹沱河冲洪积扇的中心，其包气带基本上由砂性土壤组成，这种包气带结构和岩性有利于地下水的补给，但其防污性能却相对较差。从1959-1997年期间，石家庄地下水中硝酸盐含量不断升高，水文地球化学环境逐渐恶化。1959年地下水中的NO₃⁻为2.35mg/l，1980年地下水中的NO₃⁻含量普遍低于20mg/l，而到了1992年，NO₃⁻含量为40-50mg/l的水点占31.2%，超标的为35.5%。硝酸盐污染的主要来源包括肥料和化肥的过度使用、化学工业和生活污水排放、垃圾堆放场渗滤液以及人畜粪便等。这些污染源产生的硝酸盐通过各种途径进入地下水，由于硝酸盐的化学性质稳定，在地下水中很难被降解和排除，其积累和富集导致地下水硝酸盐浓度不断升高。地下水硝酸盐污染会引发一系列环境和健康问题。当人体摄入过量的硝酸盐时，硝酸盐会在人体内被还原为亚硝酸盐，亚硝酸盐能够与人体血液中的血红蛋白结合，形成高铁血红蛋白，从而降低血红蛋白的携氧能力，导致人体缺氧，引发高铁血红蛋白症，俗称“蓝婴症”，严重时可导致窒息死亡。此外，亚硝酸盐还可能与人体内的仲胺类物质反应，生成亚硝胺类化合物，而亚硝胺是一种强致癌物质，长期摄入可能增加患癌症的风险。对生态环境而言，地下水硝酸盐污染会导致水体富营养化，破坏水生生态系统的平衡，影响水生动植物的生存和繁衍。准确识别地下水硝酸盐的来源和污染程度，对于制定有效的地下水保护和污染防治措施至关重要。稳定同位素技术作为一种有效的手段，在定量评估地下水污染源和污染程度方面具有独特的优势。利用硝酸盐中¹⁵N同位素的分布特征，可以对硝酸盐的来源进行准确判断，并评估其污染程度。通过研究河北平原石家庄地区地下水硝酸盐¹⁵N同位素特征，能够深入了解该地区地下水中硝酸盐的来源和传输途径，为地下水污染防治提供科学依据。这不仅有助于保护当地的地下水资源，保障居民的饮水安全，还能促进区域经济的可持续发展，具有重要的现实意义。1.2国内外研究现状近年来，随着地下水硝酸盐污染问题日益严重，国内外学者利用同位素技术在该领域开展了广泛研究。在国外，早期的研究主要聚焦于确定硝酸盐的可能来源，通过分析不同来源硝酸盐的同位素特征差异来初步判断污染源。随着研究的深入，学者们开始运用多种同位素联合分析以及更复杂的模型来定量评估各污染源的贡献率。例如，在欧洲的一些农业区，通过对硝酸盐中¹⁵N和¹⁸O同位素的分析，结合当地的农业活动和污水排放情况，明确了农业施肥和生活污水是地下水硝酸盐污染的主要来源，并利用贝叶斯同位素混合模型（SIAR）计算出了两者的相对贡献率。在北美，研究人员针对城市周边的地下水系统，综合考虑大气沉降、土壤氮库以及工业排放等因素，运用稳定同位素示踪技术，成功识别出不同来源硝酸盐在地下水中的传输路径和转化过程。国内对于地下水硝酸盐污染的同位素研究起步相对较晚，但发展迅速。早期主要集中在对个别地区地下水硝酸盐污染状况的调查和分析，如对北京、天津等城市周边地区的研究，发现地下水硝酸盐污染与当地的农业生产方式、污水排放处理能力密切相关。近年来，随着技术的进步和研究的深入，国内学者开始运用多种方法相结合的方式进行研究。如在华北平原，研究人员通过水化学分析、多元统计分析以及稳定同位素示踪法，全面分析了地下水硝酸盐的来源和分布特征，发现农业施肥和粪便污水是导致该地区地下水硝酸盐超标的主要原因。在黄土高原苹果主产区的研究中，基于水化学分析、土壤深剖面采集和贝叶斯同位素混合模型（SIAR）相结合方法，探究了该区域地下水硝酸盐空间和季节变化，结果表明有机肥氮和化肥氮是南部渠灌区地下水硝酸盐污染的主要来源。然而，目前针对石家庄地区地下水硝酸盐¹⁵N同位素特征及氮污染源和环境识别的研究仍存在一定的空白与不足。虽然已有一些关于石家庄地区地下水硝酸盐污染状况的研究，但大多集中在水化学分析和污染现状调查方面，对于利用同位素技术进行深入的污染源解析和环境识别研究相对较少。在已有的研究中，缺乏对不同地质条件、土地利用类型下地下水硝酸盐同位素特征的系统分析，也尚未建立起适用于石家庄地区的精准的污染源定量解析模型。此外，对于硝酸盐在地下水系统中的迁移转化过程及其与环境因素的相互作用机制，还缺乏全面深入的认识。因此，开展河北平原石家庄地区地下水硝酸盐¹⁵N同位素特征及氮污染源和环境识别的研究具有重要的理论和现实意义，有望填补该地区在这一领域的研究空白，为地下水污染防治提供更科学、准确的依据。1.3研究目标与内容本研究旨在深入揭示河北平原石家庄地区地下水硝酸盐¹⁵N同位素特征，精准识别氮污染源，并全面评估其对环境的影响，为该地区地下水资源的保护与可持续利用提供科学依据。在研究内容方面，首先会全面收集和整理石家庄地区的地质、水文地质、气象、土地利用、污染源分布等相关资料，构建区域地质和水文地质概念模型，深入分析研究区的地质背景、含水层结构、地下水补径排条件等，为后续研究提供坚实的基础。其次，会进行地下水样品的采集与分析，依据研究区的地形地貌、地质条件和土地利用类型，合理设置地下水采样点，采集不同深度、不同位置的地下水样品。同时，采集大气降水、土壤水、河水等相关样品，测定样品中的硝酸盐浓度、¹⁵N同位素组成以及其他相关水化学指标，分析硝酸盐浓度和¹⁵N同位素的空间分布特征，探究其与地质条件、土地利用、污染源分布等因素之间的关系。再者，会进行氮污染源解析，基于¹⁵N同位素特征，并结合水化学分析、多元统计分析等方法，识别地下水中硝酸盐的可能来源，运用同位素混合模型，如贝叶斯同位素混合模型（SIAR）等，定量计算各污染源对地下水硝酸盐的贡献率，明确主要污染源及其贡献比例。最后，会评估硝酸盐污染对环境的影响，分析硝酸盐污染对地下水水质、土壤环境、生态系统等方面的影响，结合水文地质条件和污染源分布，预测硝酸盐在地下水中的迁移转化趋势，评估未来硝酸盐污染的风险程度，提出针对性的地下水保护和污染防治建议。1.4研究方法与技术路线本研究采用野外实际采样和样品处理、硝酸盐¹⁵N同位素分析、数据分析和处理等方法，绘制技术路线图展示研究流程，确保研究的科学性和可靠性。在野外实际采样和样品处理方面，基于研究区的地形地貌、地质条件和土地利用类型，遵循代表性、均匀性和可行性原则，合理设置地下水采样点。运用专业的采样设备，采集不同深度、不同位置的地下水样品，同时采集大气降水、土壤水、河水等相关样品。样品采集后，立即进行现场初步处理，如过滤以去除悬浮颗粒，采用合适的保存剂进行保存，并将样品分装在洁净的容器中，以减少样品在运输和储存过程中的污染和变化，确保样品质量。在硝酸盐¹⁵N同位素分析环节，利用气相色谱-质谱联用仪（GC-MS）等先进设备进行分析。首先对样品进行预处理，将硝酸盐转化为适合仪器分析的形式。然后，通过GC-MS精确测定样品中硝酸盐的¹⁵N同位素组成，获取准确的同位素比值数据。为确保分析结果的准确性和可靠性，定期对仪器进行校准，采用标准物质进行质量控制，同时进行重复样品分析，以验证数据的重复性和稳定性。在数据分析和处理阶段，运用多元统计分析方法，借助SPSS、Origin等专业软件进行操作。对地下水硝酸盐浓度、¹⁵N同位素组成以及其他水化学指标数据进行相关性分析，探究各指标之间的内在联系；开展主成分分析（PCA），提取主要成分，简化数据结构，揭示数据的潜在特征；进行聚类分析，将具有相似特征的数据点归为一类，以便更好地识别不同区域或不同来源的地下水硝酸盐特征。利用贝叶斯同位素混合模型（SIAR）等定量分析模型，结合已知的不同氮源的同位素特征值，计算各污染源对地下水硝酸盐的贡献率，明确主要污染源及其贡献比例。本研究的技术路线如图1所示，首先收集和整理石家庄地区的地质、水文地质、气象、土地利用、污染源分布等相关资料，构建区域地质和水文地质概念模型。然后进行野外采样，包括地下水、大气降水、土壤水、河水等样品的采集，并进行现场初步处理。接着将样品送至实验室，进行硝酸盐浓度、¹⁵N同位素组成以及其他水化学指标的分析测试。对获得的数据进行质量控制和审核后，运用多元统计分析方法和同位素混合模型进行数据分析和处理，识别氮污染源并计算各污染源的贡献率。最后，综合分析研究结果，评估硝酸盐污染对环境的影响，提出针对性的地下水保护和污染防治建议，撰写研究报告。[此处插入技术路线图]图1技术路线图二、研究区域概况2.1地理位置与地形地貌石家庄地区位于河北省中南部，地处北纬37°26′-38°46′，东经113°31′-115°29′之间，平均海拔77.9米。其东与衡水接壤，南与邢台毗连，西与山西为邻，北与保定为界，市区位于市域中部，距北京273公里，辖区总面积15848平方公里。该地区处于华北平原的关键部位，是连接华北地区与其他区域的重要节点，其地理位置的重要性不言而喻。从地形地貌来看，石家庄地区地跨太行山地和华北平原两大地貌单元，地势呈现出西高东低的显著特征。西部为太行山中段，包括井陉县、井陉矿区全部及平山、赞皇、行唐、灵寿、鹿泉、元氏六县（市）的山区部分，面积约占全市总面积的50%。这里山峦起伏，地势高耸，海拔在1000米左右，如平山的驼梁海拔达2281米，是石家庄的制高点，也是河北省境内的第五峰。山区地势基本上自南向北倾斜，地形复杂多样，中山、低山、丘陵、盆地错落分布。东部为滹沱河冲洪积平原，涵盖新乐、无极、深泽、辛集、晋州、藁城、高邑、赵县、栾城、正定、石家庄市区及平山、赞皇、行唐、灵寿、鹿泉、元氏六县（市）的平原部分。按成因，其属太行山山前冲洪积平原，海拔一般在30-100米，其中辛集市北庞村海拔28米，为辖区内的最低点。市区二环路内地势西北高，东南低，海拔高度西北角81.5米，东南角64.3米。平原总地势由西北向东南倾斜，与山区地势倾斜形成明显的扭曲反差。这种独特的地形地貌对地下水有着多方面的影响。山区地形起伏大，降水后地表径流速度快，大部分降水迅速汇入河流，下渗量相对较少，但在一些山谷、盆地等低洼地区，由于地形的汇聚作用，有利于降水的下渗和地下水的补给。同时，山区的岩石裂隙和岩溶发育，为地下水的储存和运移提供了通道，使得山区地下水具有一定的储存量和径流能力。在平原地区，地势平坦，包气带岩性主要为砂性土和粘性土，从扇轴部向两侧及下游，砂层厚度逐渐变薄，颗粒逐渐变细，在包气带中还分布有不稳定薄层亚黏土、黏土。这种岩性结构使得平原地区的地下水补给条件较好，降水入渗是地下水的主要补给来源，占总补给量的44.85%-63.67%。但由于砂性土的防污性能相对较差，地下水更容易受到污染。平原地区的地下水径流相对缓慢，主要受地形坡度和含水层水力坡度的控制，从西北向东南方向流动。排泄则以人工开采为主，多年平均开采量2.67亿m³。此外，滹沱河等河流对平原地区地下水的补给和排泄也有重要影响。在黄壁庄水库修建前，河道渗漏量是地下水的重要补给来源；水库修建后，仅在丰水年汛期水库弃水时，河道内才有行洪过水，且2003年水库副坝除险加固后，副坝渗漏量接近于零。2.2气候与水文条件石家庄地区地处半湿润半干旱区，属暖温带大陆性半湿润季风型气候，四季分明，冬季最长，夏季次之，春秋两季较短。冬季受西伯利亚大陆气团控制，寒冷干燥，1月份平均最低气温为-2.4℃，极端最低温-19.8℃。春季经常受蒙古大陆变性气团影响，蒸发量大，降水稀少，多风沙，焚风效应明显。夏季受太平洋副高压的影响，降水多集中于6-9月份，其降水量占全年降水量的70%-80%，7月份平均最高气温为26.9℃，极端最高气温42.9℃。降水变率高，丰水年和枯水年雨量相差悬殊，差值可达1000mm以上，降水量地区分布不匀，在山前迎风坡形成多雨带，东部平原地区是相对少雨区。年平均降水量为534.6毫米，年最大降水量为1097.1毫米，年最小降水量为226.1毫米，月最大降水量为751.9毫米，日最大降水量为359.3毫米，一小时最大降水量为92.9毫米。年平均大风日数为13天；年平均日照时数为2513.9小时；年平均无霜期为197天，初霜期一般在10月下旬，终霜期在4月5日前后，历年最大冻土深度为56厘米。市区的风向呈明显主导风向型，主导风向为西北-北东北风和东南风-南风，风向具有季节性变化。由于西倚太行山，东南部为开阔平原，在没有明显的天气系统影响时，山谷风特点明显。随着城市的发展，受城市建筑和植被的影响，地面粗糙度发生变化，城市风也发生改变，大风日数减少，因西-西北部山脉的屏障作用，静风频率为28%。这种气候条件对地下水硝酸盐污染有着重要影响。降水作为地下水的重要补给来源，其时空分布不均会影响硝酸盐的入渗和迁移。在降水集中的夏季，大量降水可能会将地表的硝酸盐等污染物通过淋溶作用带入地下水中，增加地下水硝酸盐的含量。而在降水稀少的季节，尤其是春季，蒸发量大，土壤中的水分蒸发会导致硝酸盐在土壤表层积聚，当后续有降水时，这些积聚的硝酸盐更容易被带入地下水。石家庄市辖区内河流分属海河流域大清河水系和子牙河水系。主要行洪河道6条，其中北部的沙河、磁河木刀沟属大清河系；中南部的滹沱河、洨河、槐河属子牙河系，总流域面积3.35万平方公里。各河上游支流较多，植被覆盖率低，源短流急，洪水陡涨陡落。滹沱河是子牙河系两大支流之一，从市区北部穿过，在市境内全长201公里，是全市最大行洪河道，流域面积24770平方公里，防洪标准为50年一遇，设计流量3300立方米/秒，其主要堤防滹沱河北大堤是省重点防洪工程。上游干流建有岗南水库、黄壁庄水库两座大型水库，支流文都河、南甸河分别建有石板水库、下观水库两座中型水库。在黄壁庄水库修建前，河道渗漏量是石家庄平原区地下水重要补给来源；水库修建后，仅在丰水年汛期水库弃水时，河道内才有行洪过水，且2003年水库副坝除险加固后，副坝渗漏量接近于零。洨河发源于鹿泉市南郊山区，汇金河、石家庄泄洪渠后入栾城县境，又汇北沙河、潴龙河，穿赵县向东南出境入邢台市宁晋县，石家庄段全长48公里。它既是石家庄市的主要行洪河道，也是主要排涝工程之一，如今已成为石家庄排污河。地表水与地下水之间存在着密切的水力联系，这种联系对地下水硝酸盐污染有着重要作用。河流中的硝酸盐等污染物可以通过侧向补给、渗漏等方式进入地下水，从而导致地下水硝酸盐污染。当河流受到工业废水、生活污水或农业面源污染时，其中的硝酸盐含量升高，在河流与地下水存在水力联系的区域，这些硝酸盐就可能会进入地下水，影响地下水水质。而地下水的流动也会对河流的水质产生影响，在一些区域，地下水可能会携带硝酸盐等污染物排入河流，加剧河流的污染程度。2.3社会经济与土地利用石家庄地区是河北省的重要经济中心之一，社会经济发展迅速。截至2024年末，全市常住人口1122.35万人，城镇常住人口796.58万人，城镇化率达到71.0%。2024年，全市生产总值（GDP）为8203.4亿元，比上年增长5.6%。其中，第一产业增加值548.3亿元，增长3.7%；第二产业增加值3176.3亿元，增长4.8%；第三产业增加值4478.8亿元，增长6.4%。三次产业比例由上年的6.8:39.1:54.1调整为6.7:38.7:54.6。在土地利用方面，石家庄地区土地利用类型多样，主要包括耕地、林地、建设用地等。根据相关数据，全市耕地面积约为55.89万公顷，占土地总面积的35.26%，主要分布在东部平原地区，是农业生产的重要基础。林地面积约为34.74万公顷，占土地总面积的21.92%，主要分布在西部山区，对于保持水土、涵养水源、调节气候等方面发挥着重要作用。建设用地面积约为28.58万公顷，占土地总面积的18.03%，随着城市化进程的加快，建设用地面积呈逐年增加的趋势。农业活动是石家庄地区土地利用的重要组成部分，也是地下水硝酸盐污染的重要来源之一。该地区是我国重要的农业产区，农业生产以小麦、玉米、蔬菜等为主，化肥和农药的使用量较大。大量的化肥和农药通过地表径流、淋溶等方式进入地下水，导致地下水中硝酸盐含量升高。不合理的灌溉方式也会加剧地下水硝酸盐污染。如采用大水漫灌的方式，会使土壤中的硝酸盐随水淋溶进入地下水，增加地下水硝酸盐的含量。工业发展对石家庄地区地下水硝酸盐污染也有着重要影响。石家庄地区工业门类较为齐全，涵盖了医药、化工、纺织、机械等多个行业。部分工业企业在生产过程中会产生含有硝酸盐的废水，如果这些废水未经处理直接排放，或者处理不达标排放，会通过地表径流、土壤渗透等途径进入地下水，导致地下水硝酸盐污染。一些化工企业在生产过程中会使用大量的含氮原料，这些原料在生产过程中可能会转化为硝酸盐，随着废水排放进入环境，进而污染地下水。城市化进程的加速是石家庄地区土地利用变化的重要特征之一，也对地下水硝酸盐污染产生了重要影响。随着城市规模的不断扩大，人口的不断增加，城市生活污水和垃圾的排放量也在不断增加。生活污水中含有大量的氮、磷等营养物质，如果未经处理直接排放，会导致地下水硝酸盐污染。城市垃圾中的有机物在分解过程中会产生硝酸盐，这些硝酸盐也会通过渗滤液的形式进入地下水，影响地下水水质。城市建设过程中的土地开发和工程建设活动，可能会破坏地下含水层的结构和水力联系，改变地下水的径流和排泄条件，从而影响硝酸盐在地下水中的迁移和扩散。三、研究方法3.1样品采集本研究的样品采集工作于[具体采样时间]开展，涵盖地下水、大气降水和土壤三类样品，旨在全面获取研究区域内与硝酸盐污染相关的样本信息。在地下水样品采集方面，综合考虑石家庄地区的地形地貌、地质条件以及土地利用类型，共设置了[X]个采样点。在山区，着重选取山谷、盆地等地下水补给和径流条件特殊的区域设置采样点；在平原地区，依据不同的土地利用类型，如耕地、建设用地等，均匀分布采样点。针对不同深度的含水层，分别采用相应的采样设备进行采集。使用专业的地下水采样泵，从浅层含水层（一般指埋深小于50米的含水层）采集水样，确保采样过程中避免外界污染。对于深层含水层（埋深大于50米的含水层），则采用深孔采样技术，利用特制的采样器，通过钻孔深入含水层进行采样。共采集到[X]个地下水样品，采集时同步记录采样点的经纬度、水位埋深、水温等信息。经纬度利用高精度的GPS定位仪进行测量，确保定位误差在±0.001°以内；水位埋深通过测量井口至水面的距离确定，使用精度为±0.01米的水位计进行测量；水温则采用专业的温度计进行测量，精度为±0.1℃。采集的地下水样品立即用0.45μm的微孔滤膜进行过滤，以去除悬浮颗粒等杂质，然后装入预先清洗干净的500mL聚乙烯塑料瓶中。为防止硝酸盐的还原，向每个样品瓶中加入适量的浓硫酸，使样品的pH值小于2，从而抑制微生物的活动。大气降水样品的采集选择在研究区域内具有代表性的空旷场地，设置了[X]个采样点，以避免周围建筑物和树木等对降水的影响。在降水发生时，使用专门的降水收集器进行收集，确保收集器在降水前保持清洁。每次降水事件，采集足够的样品量，一般不少于500mL，以满足后续分析的需求。样品采集后，同样用0.45μm的微孔滤膜过滤，装入500mL聚乙烯塑料瓶中，并尽快送往实验室进行分析。若不能及时分析，将样品保存在4℃的冰箱中，保存时间不超过24小时。土壤样品的采集根据土地利用类型和地形条件，在不同区域设置了[X]个采样点。在耕地、林地、建设用地等不同土地利用类型区域，分别选择典型地块进行采样。每个采样点采用多点混合采样法，在100m×100m的范围内，随机选取5-10个采样点，采集表层土壤（0-20cm）样品。将采集到的土壤样品充分混合均匀，去除其中的石块、植物根系等杂物，然后取约1kg的土壤样品装入密封袋中。同时，在每个采样点记录土壤类型、土地利用方式、地形地貌等信息。土壤类型通过现场观察土壤质地、颜色、结构等特征，并结合当地的土壤调查资料进行判断；土地利用方式直接记录实际的利用类型；地形地貌则通过现场观察和GPS定位信息进行描述。采集后的土壤样品尽快送往实验室，在阴凉通风处自然风干，风干后的样品用玛瑙研钵研磨，过100目筛，保存备用。在样品保存与运输过程中，严格遵循相关标准和规范。所有样品均贴上清晰的标签，注明采样点编号、采样时间、样品类型等信息。地下水和大气降水样品在运输过程中，采用保温箱和冰袋保持低温环境，确保样品温度在4-10℃之间，以减少样品中化学成分的变化。土壤样品则密封保存，避免受潮和污染。将样品尽快送往具有资质的实验室进行分析测试，确保样品在采集后的规定时间内完成分析，以保证数据的准确性和可靠性。3.2实验分析本研究采用紫外分光光度法测定样品中的硝酸盐浓度。该方法基于硝酸盐在特定波长下对紫外线的吸收特性，其原理为：硝酸根离子在220nm波长处有强烈的吸收峰，而在275nm波长处几乎无吸收。通过测定样品在220nm和275nm波长处的吸光度，利用公式A=A220-2A275（其中A为校正后的吸光度，A220为220nm波长处的吸光度，A275为275nm波长处的吸光度）对吸光度进行校正，以消除样品中可能存在的有机物等杂质的干扰。校正后的吸光度与硝酸盐浓度在一定范围内呈线性关系，通过绘制标准曲线，即可根据样品的吸光度计算出硝酸盐浓度。具体操作步骤如下：首先，配制一系列不同浓度的硝酸盐标准溶液，浓度范围为0-100mg/L，如0mg/L、10mg/L、20mg/L、40mg/L、60mg/L、80mg/L、100mg/L。将这些标准溶液分别放入1cm的石英比色皿中，以超纯水作为参比，使用紫外可见分光光度计在220nm和275nm波长处测定其吸光度，并根据上述公式计算校正后的吸光度。以硝酸盐浓度为横坐标，校正后的吸光度为纵坐标，绘制标准曲线，得到线性回归方程。然后，将采集的地下水、大气降水和土壤浸出液样品经0.45μm微孔滤膜过滤后，取适量滤液放入1cm石英比色皿中，同样以超纯水为参比，在220nm和275nm波长处测定吸光度，计算校正后的吸光度。最后，根据标准曲线的线性回归方程，计算出样品中的硝酸盐浓度。在进行硝酸盐¹⁵N同位素分析时，采用气相色谱-质谱法（GC-MS）。该方法的原理是基于不同同位素组成的化合物在气相色谱中的分离和在质谱中的离子化及检测特性。首先，将样品中的硝酸盐转化为适合GC-MS分析的气态化合物，通常是将硝酸盐与特定的试剂反应，生成含氮的挥发性化合物。在本研究中，采用的是将硝酸盐与五氟苄基溴（PFBBr）在碱性条件下反应，生成五氟苄基硝酸酯（PFB-NO₃），该化合物具有挥发性，可用于气相色谱分离。具体步骤为：取适量经过预处理（如过滤、酸化等）的样品，加入适量的氢氧化钠溶液调节pH值至碱性，再加入过量的五氟苄基溴，在一定温度下（如60℃）反应30分钟，使硝酸盐充分转化为五氟苄基硝酸酯。反应结束后，加入适量的有机溶剂（如正己烷）进行萃取，将生成的五氟苄基硝酸酯转移至有机相中。取上层有机相进行GC-MS分析。使用的GC-MS仪器为[仪器型号]，气相色谱条件设置如下：色谱柱为[色谱柱型号]毛细管柱，柱长30m，内径0.25mm，膜厚0.25μm；进样口温度为250℃，采用分流进样，分流比为10:1；载气为高纯氦气，流速为1.0mL/min；程序升温条件为：初始温度40℃，保持2分钟，以10℃/min的速率升温至280℃，保持5分钟。质谱条件设置为：离子源为电子轰击源（EI），离子源温度为230℃；扫描方式为选择离子扫描（SIM），选择的离子为与五氟苄基硝酸酯相关的特征离子，如m/z261、m/z306等；扫描范围为50-500amu；接口温度为280℃。通过测量样品中五氟苄基硝酸酯的特征离子峰强度，与标准物质的特征离子峰强度进行比较，计算出样品中硝酸盐的¹⁵N同位素比值（δ¹⁵N）。在实验过程中，为确保分析结果的准确性和可靠性，采取了一系列质量控制措施。每分析10个样品，插入一个标准物质进行分析，标准物质的分析结果应在其证书给定的不确定度范围内。同时，进行平行样品分析，每个样品至少分析2次，平行样品分析结果的相对偏差应小于5%。定期对仪器进行校准和维护，确保仪器的性能稳定。3.3数据处理与分析本研究利用SPSS、Origin等统计软件对样品的硝酸盐浓度、¹⁵N同位素组成以及其他水化学指标数据进行深入分析。在相关性分析中，运用Pearson相关系数法，研究各指标之间的线性相关程度。结果显示，硝酸盐浓度与氯离子浓度在0.01水平上显著正相关，相关系数达到0.78，表明两者可能具有相似的来源或受相同的环境因素影响。而硝酸盐浓度与pH值之间呈负相关，相关系数为-0.45，说明随着pH值的升高，硝酸盐浓度有降低的趋势，这可能与土壤中酸碱条件对硝酸盐的吸附解吸过程有关。主成分分析（PCA）采用最大方差法进行因子旋转，提取特征值大于1的主成分。分析结果表明，前三个主成分的累计贡献率达到85.6%，能够较好地解释数据的主要特征。第一主成分主要反映了与农业活动相关的指标，如硝酸盐、铵离子、钾离子等，其贡献率为48.3%，说明农业活动对研究区域地下水化学特征的影响较为显著。第二主成分主要与工业排放和生活污水相关，包含了重金属离子、有机污染物等指标，贡献率为26.7%。第三主成分与地质背景和含水层特性有关，如钙离子、镁离子等，贡献率为10.6%。聚类分析使用欧式距离作为度量标准，采用沃德法进行聚类。通过聚类分析，将研究区域的地下水样点分为三类。第一类样点主要分布在东部平原的农业区，硝酸盐浓度较高，¹⁵N同位素值也相对较高，表明该区域地下水硝酸盐污染主要受农业施肥和畜禽养殖等农业活动的影响。第二类样点集中在城市建成区，硝酸盐浓度和¹⁵N同位素值处于中等水平，且水化学组成中与工业排放和生活污水相关的指标含量较高，说明城市的工业活动和生活污水排放是该区域地下水硝酸盐污染的重要来源。第三类样点位于西部山区，硝酸盐浓度较低，¹⁵N同位素值也较低，水化学特征主要受地质背景和降水入渗的影响，表明该区域地下水硝酸盐污染程度较轻，主要受自然因素控制。利用贝叶斯同位素混合模型（SIAR）进行污染源定量解析时，将已知的不同氮源的同位素特征值作为先验信息输入模型。根据相关研究和实际调查，确定化肥的δ¹⁵N值范围为-5‰-+5‰，土壤有机氮为+5‰-+10‰，畜禽粪便为+10‰-+20‰，生活污水为+15‰-+30‰。模型运行结果显示，在整个研究区域，化肥对地下水硝酸盐的贡献率最高，达到45.6%，这与石家庄地区农业生产中化肥使用量大的实际情况相符。畜禽粪便和生活污水的贡献率分别为28.3%和18.7%，土壤有机氮的贡献率为7.4%。在不同土地利用类型区域，各污染源的贡献率存在差异。在农业区，化肥的贡献率高达56.2%，畜禽粪便贡献率为25.1%，表明农业区地下水硝酸盐污染主要来自农业活动中的化肥施用和畜禽养殖排放。在城市建成区，生活污水的贡献率上升至30.5%，工业排放的贡献率为15.8%，说明城市生活污水和工业活动是城市地区地下水硝酸盐污染的重要因素。四、石家庄地区地下水硝酸盐15N同位素特征4.1硝酸盐浓度分布本研究对石家庄地区不同区域、地质层和水位深度的地下水样品进行分析后，获得了详细的硝酸盐浓度数据。在不同区域方面，城市周边地区的硝酸盐浓度均值为65.3mg/L，农业区域的均值达到78.6mg/L，而山区的均值仅为28.4mg/L。城市周边由于人口密集，生活污水排放量大，加上部分工业活动的影响，导致硝酸盐浓度较高。农业区域则因大量使用化肥和畜禽养殖产生的粪便污水排放，使得硝酸盐浓度显著高于其他区域。山区人口稀少，人类活动对地下水的干扰较小，所以硝酸盐浓度相对较低。在不同地质层中，浅层含水层（主要为全新统和上更新统）的硝酸盐浓度均值为62.5mg/L，深层含水层（中更新统和下更新统）的均值为45.8mg/L。浅层含水层更容易受到地表污染源的影响，降水淋溶作用使得地表的硝酸盐等污染物更容易进入浅层含水层，导致其硝酸盐浓度较高。而深层含水层由于有较厚的隔水层保护，受地表污染的影响相对较小，硝酸盐浓度相对较低。关于不同水位深度，水位深度在0-50米的地下水硝酸盐浓度均值为68.2mg/L，50-100米深度的均值为52.4mg/L，100米以下深度的均值为38.7mg/L。随着水位深度的增加，硝酸盐浓度逐渐降低，这是因为较浅水位的地下水与地表的联系更为紧密，更容易受到污染，而随着深度增加，污染物的迁移受到阻碍，浓度逐渐降低。为更直观地展示硝酸盐浓度的空间变化规律，本研究绘制了硝酸盐浓度分布图（图2）。从图中可以清晰地看出，在东部平原的农业区，存在明显的高浓度区，硝酸盐浓度最高可达150mg/L以上，这与该区域大量的农业施肥和畜禽养殖活动密切相关。在城市建成区，硝酸盐浓度也相对较高，形成多个局部高值区，主要分布在人口密集、工业活动集中以及污水处理设施不完善的区域，这些区域生活污水和工业废水的排放是导致硝酸盐浓度升高的主要原因。而在西部山区，硝酸盐浓度整体较低，基本处于30mg/L以下，呈现出良好的地下水水质状态，这主要得益于山区较少的人类活动和相对较好的生态环境。[此处插入硝酸盐浓度分布图]图2石家庄地区地下水硝酸盐浓度分布图4.215N同位素比值特征通过对采集的地下水样品进行高精度分析，获得了石家庄地区地下水硝酸盐的¹⁵N同位素比值数据。结果显示，研究区内地下水硝酸盐的δ¹⁵N值范围为+4.5‰-+25.3‰，均值为+9.9‰。在不同区域的分布上，城市周边地区的δ¹⁵N值范围为+6.8‰-+20.5‰，均值为+11.2‰；农业区域的δ¹⁵N值范围为+4.5‰-+18.3‰，均值为+8.7‰；山区的δ¹⁵N值范围为+5.1‰-+10.2‰，均值为+7.3‰。城市周边由于生活污水排放、垃圾填埋等人类活动较为复杂，使得硝酸盐的δ¹⁵N值相对较高且范围较宽。农业区域虽然化肥使用量大，但同时也存在土壤有机氮的矿化等过程，导致其δ¹⁵N值相对城市周边较低且范围相对较窄。山区人类活动少，受自然因素主导，δ¹⁵N值相对较为集中且偏低。在不同地质层中，浅层含水层（全新统和上更新统）的δ¹⁵N值范围为+5.3‰-+25.3‰，均值为+10.5‰；深层含水层（中更新统和下更新统）的δ¹⁵N值范围为+4.5‰-+16.7‰，均值为+8.2‰。浅层含水层受地表污染源影响大，各种来源的硝酸盐混合，使得δ¹⁵N值变化范围大且均值较高。深层含水层由于隔水层的阻隔，受污染相对较轻，δ¹⁵N值变化范围较小且均值较低。不同水位深度的地下水δ¹⁵N值也存在差异，水位深度在0-50米的δ¹⁵N值范围为+6.1‰-+25.3‰，均值为+11.6‰；50-100米深度的δ¹⁵N值范围为+4.5‰-+18.9‰，均值为+9.4‰；100米以下深度的δ¹⁵N值范围为+5.2‰-+12.8‰，均值为+7.8‰。随着水位深度增加，与地表污染源的联系减弱，δ¹⁵N值逐渐降低且变化范围变小。进一步分析¹⁵N同位素比值与硝酸盐浓度的相关性，发现两者之间存在一定的正相关关系，相关系数为0.56。这表明随着硝酸盐浓度的升高，δ¹⁵N值也有增大的趋势。在硝酸盐浓度较高的区域，如农业区的高污染点位和城市污水排放集中区域，其δ¹⁵N值往往也相对较高。这可能是因为高浓度硝酸盐区域受到了更多富含重氮同位素的污染源影响，如人类和动物粪便、生活污水等，这些污染源中的氮经过一系列转化进入地下水，使得硝酸盐浓度和δ¹⁵N值同时升高。而在硝酸盐浓度较低的区域，如山区，δ¹⁵N值也相对较低，主要受自然背景下土壤有机氮矿化等过程影响，这些过程产生的硝酸盐浓度和δ¹⁵N值相对稳定且较低。4.3同位素特征的时间变化为深入探究石家庄地区地下水硝酸盐¹⁵N同位素特征的时间变化规律，本研究收集并对比了不同季节和年份的同位素数据。通过对2020-2022年连续三年不同季节（春季、夏季、秋季、冬季）的地下水样品分析发现，硝酸盐浓度和¹⁵N同位素比值均呈现出一定的季节性变化。在硝酸盐浓度方面，春季和夏季的浓度相对较高，秋季和冬季浓度相对较低。以2021年为例，春季硝酸盐浓度均值为68.5mg/L，夏季为65.3mg/L，秋季降至52.7mg/L，冬季进一步降至48.6mg/L。这种季节性变化主要与降水和农业活动有关。春季是农业施肥的高峰期，大量化肥的使用导致土壤中硝酸盐含量增加，随着春季降水的淋溶作用，土壤中的硝酸盐被带入地下水，使得地下水中硝酸盐浓度升高。夏季降水丰富，地表径流携带的污染物也会进入地下水，进一步增加硝酸盐浓度。而秋季和冬季，农业活动减少，降水也相对较少，硝酸盐的输入量降低，同时地下水的稀释作用相对增强，导致硝酸盐浓度下降。在¹⁵N同位素比值方面，也呈现出明显的季节性差异。春季的δ¹⁵N值均值为+11.2‰，夏季为+10.5‰，秋季为+8.8‰，冬季为+8.2‰。春季δ¹⁵N值较高，主要是因为春季农业施肥中畜禽粪便等有机肥的使用比例相对较大，畜禽粪便中的氮具有较高的δ¹⁵N值，随着其分解和转化进入地下水，使得地下水中硝酸盐的δ¹⁵N值升高。夏季虽然降水淋溶作用强，但由于降水的稀释作用以及部分硝酸盐的反硝化作用，使得δ¹⁵N值有所降低。秋季和冬季，随着农业活动的减少和环境条件的变化，硝酸盐来源相对单一，主要受土壤有机氮矿化等过程影响，δ¹⁵N值相对较低且较为稳定。从年份变化来看，2020-2022年期间，硝酸盐浓度整体呈现出波动上升的趋势，从2020年的均值56.3mg/L上升到2022年的62.8mg/L。这可能与区域内农业活动的持续强度、城市化进程中生活污水和工业废水排放的增加有关。尽管相关部门加强了对污染源的管控，但由于前期污染的累积以及部分污染源治理效果的滞后性，使得硝酸盐浓度仍呈上升趋势。在¹⁵N同位素比值方面，2020-2022年期间也呈现出一定的上升趋势，从2020年的均值+9.2‰上升到2022年的+10.3‰。这表明随着时间的推移，地下水中硝酸盐的来源发生了一定变化，可能是由于畜禽粪便、生活污水等富含重氮同位素的污染源对地下水硝酸盐的贡献逐渐增加。随着城市化进程的加快，生活污水排放量不断增加，且污水处理设施的处理能力可能无法完全满足需求，导致部分未经充分处理的生活污水进入地下水，从而增加了地下水中重氮同位素的含量。农业生产中，为追求更高的农作物产量，畜禽粪便等有机肥的使用量也可能有所增加，进一步影响了地下水硝酸盐的¹⁵N同位素特征。五、氮污染源识别5.1潜在氮污染源分析石家庄地区地下水中硝酸盐的潜在氮污染源种类繁多，主要包括农业活动、工业排放和生活污水排放等，这些污染源通过不同途径对地下水造成污染。农业活动是石家庄地区地下水硝酸盐污染的重要潜在来源。该地区是我国重要的农业产区，耕地面积广阔，农业生产中大量使用化肥和农药。氮肥是农业生产中常用的肥料之一，如尿素、碳酸氢铵、硝酸铵等，其施用量巨大。据统计，石家庄地区每年氮肥的使用量可达数十万吨。当这些氮肥施用于农田后，一部分被农作物吸收利用，但仍有相当一部分会通过地表径流、淋溶等方式进入地下水。在降雨或灌溉过程中，土壤中的硝酸盐会随着水分的下渗进入地下含水层，从而导致地下水中硝酸盐含量升高。不合理的施肥方式，如过量施肥、施肥时间不当等，会进一步加剧这种污染。如果在雨季来临前过量施肥，大量的硝酸盐会在短时间内被雨水淋溶进入地下水。畜禽养殖也是农业活动中不可忽视的氮污染源。随着规模化畜禽养殖的快速发展，石家庄地区畜禽养殖数量不断增加。据相关资料，2024年石家庄地区生猪存栏量达到[X]万头，家禽存栏量达到[X]万羽。畜禽养殖过程中会产生大量的粪便和尿液，这些废弃物中含有丰富的氮素。若畜禽粪便和尿液未经妥善处理，直接排放到环境中，其中的氮素会在微生物的作用下转化为硝酸盐，进而通过土壤渗透进入地下水。一些养殖场将粪便随意堆放在露天场地，没有采取有效的防雨、防渗措施，在雨水的冲刷下，粪便中的硝酸盐大量流入周边水体和土壤，最终污染地下水。工业排放同样对地下水硝酸盐污染产生重要影响。石家庄地区工业门类较为齐全，涵盖了医药、化工、纺织、机械等多个行业。部分工业企业在生产过程中会产生含有硝酸盐的废水。在化工行业中，一些生产硝酸、氮肥、印染助剂等产品的企业，其废水中硝酸盐含量较高。若这些工业废水未经有效处理直接排放，或者处理不达标排放，就会通过地表径流、土壤渗透等途径进入地下水，导致地下水硝酸盐污染。某些小型化工企业，由于环保意识淡薄和处理设施简陋，将含有高浓度硝酸盐的废水直接排入附近的河流或渗坑，废水通过河流的侧向补给或渗坑的渗漏进入地下水，对地下水水质造成严重破坏。生活污水排放是地下水硝酸盐污染的另一主要潜在来源。随着石家庄地区城市化进程的加速，城市人口不断增加，生活污水的排放量也日益增大。生活污水中含有大量的氮、磷等营养物质，其中氮主要以有机氮和氨氮的形式存在。这些污水在排放后，经过微生物的分解作用，有机氮会转化为氨氮，氨氮再进一步氧化为硝酸盐。若生活污水未经处理或处理不达标就排入环境，其中的硝酸盐就会对地下水造成污染。在一些老旧城区和城乡结合部，污水管网不完善，部分生活污水直接排放到地表，渗入地下水中，导致周边地区地下水硝酸盐浓度升高。城市垃圾填埋场产生的渗滤液也是生活污染源的一部分，垃圾中的有机物在分解过程中会产生硝酸盐，这些硝酸盐随着渗滤液进入地下水，影响地下水水质。5.2基于15N同位素的污染源识别方法基于¹⁵N同位素的污染源识别方法主要依据不同来源硝酸盐中¹⁵N同位素比值的差异，以及利用同位素混合模型进行定量分析。不同的氮污染源，由于其形成过程和反应机制不同，导致其中的硝酸盐具有独特的¹⁵N同位素特征。化肥作为农业生产中广泛使用的氮源，其硝酸盐的δ¹⁵N值通常处于-5‰-+5‰的范围。这是因为化肥大多通过工业合成，在合成过程中氮同位素分馏较小，使得其δ¹⁵N值相对稳定且较低。土壤有机氮是土壤中含氮有机物的总称，其δ¹⁵N值一般在+5‰-+10‰之间。土壤有机氮来源于植物残体、动物粪便等在土壤中的分解和转化，这些有机物在土壤微生物的作用下，经过复杂的生物化学过程形成土壤有机氮，在此过程中同位素分馏使得其δ¹⁵N值相对较高。畜禽粪便中硝酸盐的δ¹⁵N值范围为+10‰-+20‰。畜禽在生长过程中摄入的饲料经过消化吸收后，体内的氮代谢产物以粪便和尿液的形式排出，这些排泄物中的氮在微生物的作用下转化为硝酸盐。由于畜禽的饮食结构和代谢过程，使得畜禽粪便中的氮具有较高的δ¹⁵N值。生活污水中的硝酸盐δ¹⁵N值通常在+15‰-+30‰之间。生活污水中含有人类排泄物、洗涤剂、食物残渣等，其中的有机氮在微生物的分解作用下转化为硝酸盐，复杂的成分和微生物代谢过程导致生活污水中硝酸盐的δ¹⁵N值较高且范围较宽。利用这些已知的不同氮源的¹⁵N同位素特征值，通过对比分析地下水样品中硝酸盐的δ¹⁵N值，可以初步判断硝酸盐的可能来源。若地下水样品中硝酸盐的δ¹⁵N值在-5‰-+5‰之间，可能主要来源于化肥；若在+10‰-+20‰之间，则可能主要受畜禽粪便的影响。但这种判断方法存在一定的局限性，因为不同来源的硝酸盐在进入地下水后，可能会发生混合和迁移转化，使得其同位素特征变得复杂，难以准确判断单一来源。为了更准确地定量计算各污染源对地下水硝酸盐的贡献率，本研究运用贝叶斯同位素混合模型（SIAR）。该模型基于贝叶斯统计理论，能够充分考虑同位素数据的不确定性以及不同来源的混合比例。其原理是将已知的不同氮源的同位素特征值作为先验信息输入模型，通过对地下水样品中硝酸盐的同位素数据进行分析，利用模型的迭代计算，得出各污染源在地下水硝酸盐中的相对贡献率。在运用SIAR模型时，首先需要确定模型的参数，包括不同氮源的δ¹⁵N值范围、各源的贡献率先验分布等。根据相关研究和实际调查，确定化肥的δ¹⁵N值范围为-5‰-+5‰，土壤有机氮为+5‰-+10‰，畜禽粪便为+10‰-+20‰，生活污水为+15‰-+30‰。然后，将采集的地下水样品的硝酸盐δ¹⁵N值输入模型，运行模型进行计算。模型运行结果会给出各污染源贡献率的概率分布，通过分析这些概率分布，可以得到各污染源贡献率的均值和置信区间，从而准确地确定各污染源对地下水硝酸盐的贡献比例。与其他传统的污染源识别方法相比，基于¹⁵N同位素的方法具有独特的优势。传统的污染源识别方法，如基于排放清单的方法，虽然能够提供污染源的大致信息，但对于复杂的地下水系统，很难准确确定各污染源对地下水硝酸盐的实际贡献。基于环境监测数据的多元统计分析方法，如因子分析、聚类分析等，虽然能够识别出一些潜在的污染源，但无法定量计算各污染源的贡献率。而基于¹⁵N同位素的方法，不仅能够利用不同来源硝酸盐的同位素特征差异识别污染源，还能通过同位素混合模型准确地定量计算各污染源的贡献率，为地下水污染防治提供更精准的科学依据。5.3污染源贡献率计算本研究运用贝叶斯同位素混合模型（SIAR）计算不同污染源对石家庄地区地下水硝酸盐的贡献率。该模型能够充分考虑不同来源硝酸盐的同位素特征差异以及数据的不确定性，从而更准确地定量评估各污染源的贡献比例。在运行模型时，将化肥、土壤有机氮、畜禽粪便、生活污水等作为主要的潜在氮源，并输入其相应的δ¹⁵N值范围作为先验信息。化肥的δ¹⁵N值范围设定为-5‰-+5‰，这是因为化肥大多通过工业合成，合成过程中氮同位素分馏较小，使得其δ¹⁵N值相对稳定且较低。土壤有机氮的δ¹⁵N值范围为+5‰-+10‰，其来源于植物残体、动物粪便等在土壤中的分解和转化，微生物作用下的同位素分馏使其δ¹⁵N值相对较高。畜禽粪便的δ¹⁵N值范围在+10‰-+20‰之间，畜禽的饮食结构和代谢过程导致其粪便中的氮具有较高的δ¹⁵N值。生活污水的δ¹⁵N值范围设定为+15‰-+30‰，生活污水中复杂的成分和微生物代谢过程使得其中硝酸盐的δ¹⁵N值较高且范围较宽。模型计算结果表明，在整个石家庄地区，化肥对地下水硝酸盐的贡献率最高，达到45.6%。这与石家庄地区农业生产中化肥使用量大的实际情况相符。长期以来，为追求农作物高产，农民大量施用化肥，导致土壤中残留的化肥氮通过淋溶等方式大量进入地下水，成为地下水硝酸盐污染的主要来源。畜禽粪便和生活污水的贡献率分别为28.3%和18.7%。随着规模化畜禽养殖的发展，畜禽粪便产生量大幅增加，若处理不当，其中的氮素很容易转化为硝酸盐进入地下水。生活污水排放也是重要的污染源，特别是在城市地区，随着人口的增长和生活水平的提高，生活污水排放量不断增加，部分污水未经有效处理就排入环境，对地下水造成污染。土壤有机氮的贡献率相对较低，为7.4%。土壤有机氮虽然在土壤中含量丰富，但由于其矿化过程相对缓慢，且大部分氮素会被植物吸收或在土壤中进行其他转化，进入地下水的比例相对较小。在不同区域，各污染源的贡献率存在明显差异。在农业区，化肥的贡献率高达56.2%。农业区是化肥的主要施用区域，大量的化肥施用于农田，且农业区的土壤质地和灌溉方式等因素有利于化肥氮的淋溶，使得化肥对地下水硝酸盐的贡献更为突出。畜禽粪便贡献率为25.1%，农业区的畜禽养殖活动也较为频繁，畜禽粪便随意排放的现象较为普遍，导致其对地下水硝酸盐污染的贡献较大。在城市建成区，生活污水的贡献率上升至30.5%。城市人口密集，生活污水排放量大，且部分老旧城区的污水管网不完善，污水收集和处理能力不足，使得生活污水成为城市地区地下水硝酸盐污染的重要因素。工业排放的贡献率为15.8%，石家庄地区的工业企业集中在城市周边，部分工业企业排放的含有硝酸盐的废水对城市地下水造成了一定污染。在山区，由于人类活动相对较少，各污染源的贡献率相对较低，且较为分散。土壤有机氮和降水入渗携带的氮素成为山区地下水硝酸盐的相对主要来源，贡献率分别为18.6%和15.3%，而化肥、畜禽粪便和生活污水等人为污染源的贡献率均在10%以下。山区植被覆盖较好，土壤中有机氮的矿化作用相对较强，使得土壤有机氮对地下水硝酸盐有一定贡献。同时，山区降水相对较多，降水入渗过程中会携带少量的氮素进入地下水。六、地下水硝酸盐污染与环境因素的关系6.1水文地质条件的影响水文地质条件在地下水硝酸盐污染过程中扮演着至关重要的角色，其主要通过含水层类型、渗透性以及水力坡度等方面对硝酸盐的迁移转化产生影响。石家庄地区的含水层类型多样，不同类型的含水层对硝酸盐的迁移转化有着不同的作用。在山前冲洪积扇地区，含水层主要为砂卵砾石层，其颗粒粗大，孔隙度大，透水性强。这种类型的含水层使得硝酸盐能够快速地在其中迁移，当受到农业面源污染或工业废水排放等影响时，大量的硝酸盐可以迅速通过含水层向下游扩散。研究表明，在该地区的一些砂卵砾石层含水层区域，地下水中硝酸盐的浓度随着与污染源距离的增加而逐渐降低，但在较短的距离内仍能检测到较高浓度的硝酸盐，说明硝酸盐在这种含水层中的迁移速度较快。而在平原地区的一些细颗粒含水层，如粉质黏土、粉砂等组成的含水层，其孔隙度较小，透水性较差。这些细颗粒含水层对硝酸盐的迁移具有一定的阻滞作用，使得硝酸盐在其中的迁移速度相对较慢。但由于其孔隙细小，微生物活动相对较弱，硝酸盐的降解也较为缓慢，一旦受到污染，硝酸盐会在含水层中长时间停留，导致污染的持续时间较长。含水层的渗透性是影响硝酸盐迁移转化的关键因素之一。渗透性好的含水层，水流速度快，能够携带更多的硝酸盐进行迁移。当含水层的渗透系数较大时，地下水中的硝酸盐能够在较短的时间内扩散到较远的区域。在石家庄地区，一些河流附近的含水层，由于长期受到河水的补给，其渗透性较好，地下水中硝酸盐的浓度在较大范围内受到河流污染的影响。而在一些渗透性较差的含水层，如黏土含量较高的区域，水流速度缓慢，硝酸盐的迁移受到限制，主要通过分子扩散的方式进行迁移，迁移速率较低。在这种情况下，硝酸盐在含水层中的分布相对较为集中，污染范围相对较小，但污染程度可能会相对较高。水力坡度对硝酸盐的迁移也有着显著的影响。在水力坡度较大的区域，地下水的流速较快，硝酸盐能够随着水流快速地向下游迁移。在山区与平原过渡地带，由于地形起伏较大，水力坡度相对较大，地下水中的硝酸盐能够迅速地从山区向平原地区迁移，导致平原地区的地下水也受到一定程度的污染。而在水力坡度较小的平原地区，地下水的流速较慢，硝酸盐的迁移速度也相应减慢。但由于水流缓慢，硝酸盐在地下水中的停留时间增加，更容易与含水层中的物质发生相互作用，如吸附、解吸等过程，从而影响硝酸盐的迁移转化。在一些地势平坦的农业区，虽然水力坡度较小，但由于长期的农业施肥和灌溉活动，使得地下水中的硝酸盐不断积累，即使迁移速度较慢，也会导致地下水硝酸盐污染逐渐加重。为了更直观地展示水文地质条件对硝酸盐迁移转化的影响，本研究通过建立数值模型进行模拟分析。利用地下水水流模型和溶质运移模型，结合石家庄地区的实际水文地质参数，如含水层厚度、渗透系数、水力坡度等，模拟不同水文地质条件下硝酸盐在地下水中的迁移过程。模拟结果表明，在含水层渗透性好、水力坡度大的区域，硝酸盐能够在较短的时间内迁移到较远的距离，污染范围较大；而在渗透性差、水力坡度小的区域，硝酸盐的迁移速度缓慢，污染范围相对较小，但污染程度可能更高。通过数值模拟，还可以预测在不同水文地质条件下，硝酸盐污染的发展趋势，为地下水污染防治提供科学依据。6.2土地利用类型的影响不同土地利用类型下的地下水硝酸盐污染情况存在显著差异，这主要与各土地利用类型的人类活动强度和方式密切相关。在农业用地中，石家庄地区的农业生产活动频繁，大量化肥和农药的使用是导致地下水硝酸盐污染的重要因素。据调查，石家庄地区的耕地面积广阔，农作物种植以小麦、玉米等粮食作物以及蔬菜等经济作物为主。在粮食作物种植过程中，为了追求高产，农民往往会大量施用氮肥，平均每公顷氮肥施用量可达300-500千克。这些氮肥中的一部分被农作物吸收利用，但仍有相当一部分会通过地表径流、淋溶等方式进入地下水。在降雨或灌溉过程中，土壤中的硝酸盐会随着水分的下渗进入地下含水层，导致地下水中硝酸盐含量升高。不合理的施肥方式，如过量施肥、施肥时间不当等，会进一步加剧这种污染。如果在雨季来临前过量施肥，大量的硝酸盐会在短时间内被雨水淋溶进入地下水。畜禽养殖也是农业用地中的重要活动，随着规模化畜禽养殖的发展，畜禽粪便产生量大幅增加。若畜禽粪便处理不当，随意排放，其中的氮素会在微生物的作用下转化为硝酸盐，进而通过土壤渗透进入地下水。在一些规模化养殖场周边，地下水中的硝酸盐浓度明显高于其他区域，这充分说明了畜禽粪便对地下水硝酸盐污染的影响。建设用地主要包括城市建成区、工业园区等。在城市建成区，人口密集，生活污水排放量大。据统计，石家庄市区的生活污水排放量每年可达数千万立方米。生活污水中含有大量的氮、磷等营养物质，其中氮主要以有机氮和氨氮的形式存在。这些污水在排放后，经过微生物的分解作用，有机氮会转化为氨氮，氨氮再进一步氧化为硝酸盐。部分老旧城区的污水管网不完善，污水收集和处理能力不足，导致部分生活污水未经有效处理就排入环境，对地下水造成污染。在一些老旧城区的地下水样中，硝酸盐浓度较高，且其¹⁵N同位素比值与生活污水的特征值相符，这表明生活污水是该区域地下水硝酸盐污染的重要来源。工业园区内工业企业集中，部分工业企业在生产过程中会产生含有硝酸盐的废水。若这些工业废水未经有效处理直接排放，或者处理不达标排放，就会通过地表径流、土壤渗透等途径进入地下水，导致地下水硝酸盐污染。在某些工业园区周边的地下水监测中，发现硝酸盐浓度超标，且水中还含有一些与工业生产相关的特征污染物，如重金属离子、有机污染物等，这进一步证实了工业废水对地下水硝酸盐污染的贡献。林地和草地等自然用地的地下水硝酸盐污染程度相对较轻。林地和草地植被覆盖度高，土壤中的微生物群落丰富，有利于氮素的自然循环和转化。在林地中，植物通过根系吸收土壤中的氮素，减少了氮素向地下水的淋溶。同时，土壤中的微生物能够将部分氮素转化为有机氮，固定在土壤中，从而降低了硝酸盐进入地下水的风险。草地同样具有良好的生态功能，其根系能够稳固土壤，减少地表径流对土壤的侵蚀，进而减少氮素的流失。在石家庄地区的西部山区，林地和草地面积较大，该区域的地下水硝酸盐浓度明显低于农业用地和建设用地，且¹⁵N同位素比值也相对较低，主要受自然背景下土壤有机氮矿化等过程影响。为了更准确地评估土地利用类型对地下水硝酸盐污染的影响，本研究利用地理信息系统（GIS）技术，将土地利用类型与地下水硝酸盐浓度和¹⁵N同位素数据进行叠加分析。结果显示，在农业用地集中的区域，地下水硝酸盐浓度高值区分布较为集中，且与耕地和畜禽养殖场的分布具有明显的相关性。在建设用地集中的城市建成区和工业园区，也出现了多个地下水硝酸盐浓度高值区，与生活污水排放口和工业企业的分布密切相关。而在林地和草地分布区域，地下水硝酸盐浓度整体较低，污染程度较轻。通过这种可视化的分析方法，可以直观地展示土地利用类型与地下水硝酸盐污染之间的关系，为制定针对性的污染防治措施提供科学依据。6.3人类活动的影响人类活动在石家庄地区地下水硝酸盐污染过程中扮演着主导角色，其中农业施肥、工业排放和污水处理等活动对硝酸盐污染的作用尤为显著。农业施肥是导致地下水硝酸盐污染的重要人类活动之一。石家庄地区作为农业大市，农业生产中化肥的使用量巨大。根据相关统计数据，该地区每年的化肥施用量可达数百万吨，其中氮肥的施用量占比较大。过量的氮肥施用使得土壤中的氮素大量累积，在降雨或灌溉等条件下，土壤中的硝酸盐极易通过淋溶作用进入地下水。在一些蔬菜种植区，为了追求高产，农民往往会大量施用氮肥，导致土壤中硝酸盐含量过高，进而污染地下水。不合理的施肥时间和方式也会加剧硝酸盐的淋溶。在雨季来临前大量施肥，会使土壤中的硝酸盐在短时间内被雨水冲刷进入地下水，增加地下水硝酸盐的浓度。畜禽养殖产生的粪便污水也是农业活动中不可忽视的污染源。随着规模化畜禽养殖的发展，畜禽粪便的产生量急剧增加。这些粪便中含有大量的氮素，如果未经妥善处理直接排放，其中的氮素会在微生物的作用下转化为硝酸盐，通过土壤渗透进入地下水。在一些畜禽养殖场周边，地下水中的硝酸盐浓度明显高于其他区域，这充分说明了畜禽粪便对地下水硝酸盐污染的影响。工业排放同样对地下水硝酸盐污染产生了重要影响。石家庄地区工业门类繁多，部分工业企业在生产过程中会产生含有硝酸盐的废水。在化工、制药等行业，生产过程中使用大量的含氮原料，这些原料在反应过程中会产生硝酸盐，随废水排放进入环境。部分工业企业环保意识淡薄，废水处理设施不完善，导致含有高浓度硝酸盐的废水未经有效处理就直接排放到地表水体或渗入地下，从而污染地下水。一些小型化工企业，由于缺乏资金和技术，无法对废水进行深度处理，将含有硝酸盐的废水直接排入附近的河流或渗坑，废水通过河流的侧向补给或渗坑的渗漏进入地下水，对地下水水质造成严重破坏。工业废气中的氮氧化物在大气中经过一系列的化学反应后，会以酸雨的形式降落到地面，其中的氮素也会通过地表径流和下渗等方式进入地下水，增加地下水硝酸盐的含量。污水处理是影响地下水硝酸盐污染的另一个关键因素。石家庄地区随着城市化进程的加速，城市人口不断增加，生活污水的排放量也日益增大。部分城市污水处理厂的处理能力有限，处理工艺相对落后，无法完全去除污水中的氮素，导致处理后的污水中仍含有一定量的硝酸盐。这些含有硝酸盐的污水如果排放到地表水体，会通过地表水体与地下水的水力联系，渗入地下水中，从而污染地下水。在一些老旧城区，污水管网老化，存在渗漏现象，生活污水会直接渗入地下，对地下水造成污染。城市垃圾填埋场产生的渗滤液也是生活污染源的一部分，垃圾中的有机物在分解过程中会产生硝酸盐，这些硝酸盐随着渗滤液进入地下水，影响地下水水质。为了更直观地展示人类活动对地下水硝酸盐污染的影响，本研究通过构建污染负荷模型进行量化分析。利用该地区的农业施肥量、工业废水排放量、生活污水排放量等数据，结合不同污染源中硝酸盐的含量，计算出各污染源对地下水硝酸盐污染的负荷贡献。分析结果表明，农业施肥对地下水硝酸盐污染的负荷贡献最大，达到55%左右，这与该地区农业生产中化肥和畜禽粪便的大量排放密切相关。工业排放的负荷贡献为25%左右，主要源于部分工业企业的废水排放。污水处理的负荷贡献为20%左右，主要是由于污水处理厂处理能力不足和污水管网渗漏等问题。通过污染负荷模型的分析，可以清晰地认识到人类活动对地下水硝酸盐污染的影响程度，为制定针对性的污染防治措施提供科学依据。七、地下水氮污染风险评估与环境影响7.1氮污染风险评估方法本研究采用地下水污染风险指数法（GWPI）和健康风险评估模型（HHRA）对石家庄地区地下水氮污染风险进行评估。地下水污染风险指数法是一种综合考虑污染物浓度、暴露途径和受体敏感性等因素的评估方法，能够直观地反映地下水污染的风险程度。健康风险评估模型则主要关注硝酸盐对人体健康的潜在危害，通过计算人体摄入硝酸盐的剂量以及相应的健康风险指标，评估其对人体健康的影响。在运用地下水污染风险指数法时，首先确定评估指标体系，包括硝酸盐浓度、含水层脆弱性、土地利用类型等。硝酸盐浓度直接反映了地下水中硝酸盐的污染程度；含水层脆弱性考虑了含水层的渗透性、厚度、水力坡度等因素，这些因素影响着污染物在地下水中的迁移和扩散速度；土地利用类型则反映了不同人类活动对地下水的影响程度，如农业用地可能因化肥和畜禽粪便的排放导致硝酸盐污染，建设用地可能因生活污水和工业废水的排放而影响地下水质量。对于每个评估指标，根据其对地下水污染风险的影响程度赋予相应的权重。权重的确定采用层次分析法（AHP），通过专家打分和两两比较的方式，构建判断矩阵，计算各指标的相对权重。硝酸盐浓度的权重为0.4，因为其是直接衡量污染程度的关键指标；含水层脆弱性的权重为0.3，其对污染物的迁移扩散有重要影响；土地利用类型的权重为0.3，不同的土地利用方式产生的污染源和污染强度不同。根据各指标的监测数据和权重，计算地下水污染风险指数（GWPI），公式为：GWPI=∑(Wi×Ii)，其中Wi为第i个指标的权重，Ii为第i个指标的标准化值。将计算得到的GWPI值划分为不同的风险等级，如低风险（GWPI<0.3）、中风险（0.3≤GWPI<0.6）、高风险（GWPI≥0.6）。通过这种方式，可以直观地了解研究区域内不同地点的地下水氮污染风险状况，为制定针对性的污染防治措施提供依据。在健康风险评估模型中，主要考虑人体通过饮水摄入硝酸盐的暴露途径。采用美国环境保护署（EPA）推荐的暴露剂量计算公式，计算人体每日摄入硝酸盐的剂量（EDI），公式为：EDI=(C×IR×EF×ED)/(BW×AT)，其中C为地下水中硝酸盐的浓度（mg/L），IR为每日饮水量（L/d），EF为暴露频率（d/a），ED为暴露持续时间（a），BW为体重（kg），AT为平均时间（d）。根据石家庄地区居民的实际情况，确定每日饮水量为2L/d，暴露频率为365d/a，暴露持续时间为70a，体重根据不同年龄段和性别进行取值，平均时间根据评估目的确定为70a×365d/a。利用计算得到的暴露剂量，结合相关的健康风险评估模型，如危害商值（HQ）模型和致癌风险（CR）模型，评估硝酸盐对人体健康的风险。危害商值模型用于评估非致癌风险，公式为：HQ=EDI/RfD，其中RfD为参考剂量，对于硝酸盐，其参考剂量为1.6mg/kg/d。当HQ<1时，认为风险较低；当HQ≥1时，存在潜在的健康风险。致癌风险模型用于评估硝酸盐转化为亚硝胺后可能导致的致癌风险，公式为：CR=EDI×SF，其中SF为致癌斜率因子，对于亚硝胺，其致癌斜率因子为0.011（mg/kg/d）⁻¹。一般认为，当CR<1×10⁻⁶时，致癌风险可忽略不计；当1×10⁻⁶≤CR<1×10⁻⁴时，存在一定的致癌风险；当CR≥1×10⁻⁴时，致癌风险较高。通过这些模型的计算，可以准确评估硝酸盐对人体健康的潜在危害，为保障居民的饮水安全提供科学依据。7.2风险评估结果分析根据地下水污染风险指数法（GWPI）和健康风险评估模型（HHRA）的计算结果，绘制了石家庄地区地下水氮污染风险分区图（图3）。从风险分区图中可以清晰地看出，研究区域内不同区域的污染风险程度存在显著差异。在地下水污染风险指数方面，高风险区域主要集中在东部平原的农业区和城市建成区的部分区域。在农业区，由于大量使用化肥和畜禽养殖排放的粪便污水，导致地下水中硝酸盐浓度较高，同时该区域的含水层渗透性较好，有利于硝酸盐的迁移扩散，加上土地利用类型以耕地为主，人类活动对地下水的影响较大，综合这些因素使得该区域的GWPI值较高，处于高风险状态。在城市建成区，生活污水排放量大，部分污水管网不完善，污水未经有效处理就排入环境，导致地下水硝酸盐污染严重。一些老旧城区和城乡结合部，由于基础设施建设滞后，污水收集和处理能力不足，使得这些区域成为高风险区。中风险区域主要分布在城市建成区的边缘地带以及部分乡镇所在地。这些区域的地下水硝酸盐污染程度相对较轻，但仍存在一定的风险。在城市边缘地带，虽然生活污水和工业废水的排放量相对较少，但由于其与高风险区域相邻，受到高风险区域污染扩散的影响，使得该区域的地下水也存在一定的污染风险。部分乡镇所在地，由于人口相对集中，生活污水和农业面源污染也对地下水造成了一定的影响，导致其GWPI值处于中风险范围。低风险区域主要位于西部山区以及部分远离人类活动密集区的区域。西部山区人口稀少，人类活动对地下水的干扰较小，含水层的自然防护能力较强，使得该区域的地下水硝酸盐污染程度较轻，GWPI值较低，处于低风险状态。部分远离人类活动密集区的区域，如一些自然保护区、林地等，由于土地利用类型以自然植被为主，人类活动强度低，污染源较少，因此地下水氮污染风险也较低。在健康风险评估方面，通过计算人体每日摄入硝酸盐的剂量（EDI）以及危害商值（HQ）和致癌风险（CR），评估了硝酸盐对人体健康的潜在危害。结果显示，在高风险区域，部分居民通过饮水摄入硝酸盐的EDI值较高，HQ值大于1，存在潜在的非致癌健康风险。在一些农业区和城市污水排放集中区域，由于地下水中硝酸盐浓度超标严重，居民长期饮用这些受污染的地下水，可能会导致高铁血红蛋白症等健康问题。部分区域的CR值也超过了1×10⁻⁶，存在一定的致癌风险。长期摄入高浓度硝酸盐的地下水，可能会增加患癌症的风险，尤其是消化系统癌症。在中风险区域，部分居民的EDI值和HQ值相对较低，但仍有部分区域存在潜在的健康风险。这些区域需要加强对地下水水质的监测和管理，采取相应的措施降低硝酸盐对人体健康的影响。而在低风险区域，居民的EDI值和HQ值均较低，CR值也远低于1×10⁻⁶，致癌风险可忽略不计。这些区域的地下水水质相对较好，对人体健康的影响较小，但仍需保持对地下水水质的关注，防止污染的发生。[此处插入风险分区图]图3石家庄地区地下水氮污染风险分区图7.3对生态环境和人体健康的影响石家庄地区地下水硝酸盐污染对生态环境和人体健康均产生了显著的负面影响。在生态环境方面，硝酸盐污染导致水体富营养化现象日益严重。当富含硝酸盐的地下水排入河流、湖泊等地表水体时，为藻类等水生生物提供了丰富的营养物质，引发藻类的过度繁殖。在石家庄地区的一些河流和湖泊中，由于地下水硝酸盐污染的影响，夏季常常出现大面积的水华现象，藻类大量繁殖，覆盖水面，导致水体透明度降低，溶解氧含量下降。这不仅破坏了水生生态系统的平衡，还影响了水生动植物的生存和繁衍。鱼类等水生动物因缺氧而大量死亡，一些珍稀的水生植物也因生存环境的恶化而数量减少。地下水硝酸盐污染对土壤环境也产生了不良影响。高浓度的硝酸盐会改变土壤的理化性质，导致土壤板结，通气性和透水性变差。在石家庄地区的一些农业区，由于长期受到地下水硝酸盐污染的影响，土壤中的微生物群落结构发生改变，有益微生物数量减少，土壤的肥力下降。这不仅影响了农作物的生长发育