旌阳区天元镇水源区水文地球化学特征与氨氮污染模拟预测研究

旌阳区天元镇水源区水文地球化学特征与氨氮污染模拟预测研究一、引言1.1研究背景与意义水是生命之源，是人类社会赖以生存和发展的重要自然资源。随着经济的快速发展和人口的持续增长，水资源的合理开发、利用与保护已成为全球关注的焦点问题。旌阳区天元镇水源区作为当地重要的供水水源地，其水资源的质量和可持续性直接关系到区域内居民的生活用水安全以及工农业生产的稳定发展。水文地球化学是研究地下水化学成分的形成、分布、迁移及其与环境相互作用的学科，对于深入理解地下水系统的演化规律、水资源的形成机制以及合理开发利用具有重要意义。通过对水文地球化学特征的研究，可以揭示地下水与岩石、土壤之间的物质交换过程，了解地下水的补给、径流和排泄条件，为水资源评价和管理提供科学依据。同时，水文地球化学指标还可以作为地下水污染的重要指示标志，帮助识别潜在的污染源和污染途径，为地下水污染防治提供技术支持。在众多的地下水污染物中，氨氮是一种常见且危害较大的污染物。氨氮主要来源于生活污水、工业废水、农业面源污染以及垃圾填埋场渗滤液等。当氨氮进入水体后，会消耗水中的溶解氧，导致水体富营养化，破坏水生生态系统的平衡。此外，氨氮在一定条件下还可能转化为亚硝酸盐和硝酸盐，这些物质对人体健康具有潜在危害，如引发高铁血红蛋白症和产生致癌的亚硝胺等。因此，对氨氮污染的研究和控制对于保障水资源质量和生态环境安全至关重要。对于旌阳区天元镇水源区而言，开展水文地球化学特征成因分析及氨氮污染模拟预测研究具有多重重要意义。从保障居民用水安全角度来看，通过深入了解水源区水文地球化学特征，可以准确评估地下水的水质状况，及时发现潜在的水质问题，为居民提供安全可靠的饮用水源。同时，对氨氮污染的模拟预测能够提前预警氨氮污染的发展趋势，以便采取有效的防控措施，防止氨氮污染对居民健康造成威胁。从水资源保护与可持续利用方面考虑，水文地球化学特征的研究有助于揭示水源区地下水的形成和演化规律，合理确定地下水的开采量和开采方案，实现水资源的科学管理和可持续利用。而对氨氮污染的研究则可以为制定针对性的污染防治措施提供依据，减少氨氮污染对水资源的破坏，保护水源区的生态环境。从区域经济发展角度分析，稳定可靠的水资源供应是区域经济健康发展的基础。保障水源区的水质安全，能够为当地工农业生产提供优质的水资源，促进区域经济的可持续发展。反之，如果水源区受到氨氮污染等问题的困扰，将会增加水处理成本，影响工农业生产的正常进行，进而制约区域经济的发展。1.2国内外研究现状1.2.1水文地球化学特征研究进展水文地球化学特征的研究历史久远，随着科学技术的不断进步，其研究方法和分析技术也在持续革新与完善。早期，研究者主要借助化学分析手段对地下水的化学成分展开测定，以此来了解地下水的基本化学特征。随着离子色谱、电感耦合等离子体质谱（ICP-MS）等先进分析仪器的问世，对地下水中各种离子、微量元素以及同位素的分析精度和效率得到了大幅提升，能够更准确地揭示地下水的化学组成和演化过程。在研究成果方面，国内外学者已在不同地区、不同类型的地下水系统中取得了丰硕的成果。在岩溶地区，学者们通过对岩溶水的水文地球化学特征研究，揭示了岩溶作用对地下水化学成分的控制机制，以及岩溶水的循环特征与演化规律。例如，在我国西南岩溶地区，研究发现岩溶水的水化学类型主要为Ca-HCO₃型，其化学成分主要受碳酸盐岩的溶解和岩溶作用的影响。在深部地下水研究领域，相关研究成果为深入理解地球深部物质循环和地下水的形成演化提供了重要依据。有研究通过对深部热液的水文地球化学分析，探讨了深部热液与围岩之间的相互作用过程，以及深部热液的来源和运移机制。在滨海地区，针对海水入侵对地下水水文地球化学特征的影响也有大量研究，明确了海水入侵过程中地下水中盐分增加、化学类型改变等特征。然而，针对旌阳区天元镇水源区的水文地球化学特征研究仍存在明显的空白。该区域独特的地质构造、地层岩性以及人类活动影响，使得其水文地球化学特征可能具有独特性，但目前尚未有专门针对该区域的系统研究。这就需要开展深入的实地调查和分析测试工作，以揭示该区域水文地球化学特征的形成机制和分布规律。1.2.2氨氮污染研究现状氨氮污染作为水环境领域的重要研究内容，其来源广泛。生活污水中含氮有机物的分解是氨氮的重要来源之一，随着城市人口的增长和生活水平的提高，生活污水排放量不断增加，其中的氨氮含量也相应上升。工业废水如焦化、合成氨、食品加工等行业排放的废水中含有大量氨氮，这些工业废水若未经有效处理直接排放，将对水环境造成严重污染。农业面源污染方面，大量使用的氮肥以及畜禽养殖产生的粪便等，通过地表径流和淋溶作用进入水体，也会导致水体中氨氮含量升高。氨氮污染对生态环境和人类健康具有诸多危害。在生态环境方面，氨氮会消耗水中的溶解氧，导致水体富营养化，引发藻类等浮游生物大量繁殖，形成水华现象，破坏水生生态系统的平衡，使水体出现黑臭现象，影响水体景观和生态功能。对人类健康而言，氨氮在一定条件下会转化为亚硝酸盐和硝酸盐，亚硝酸盐可与人体血红蛋白结合形成高铁血红蛋白，导致人体缺氧中毒；硝酸盐则可能在人体内被还原为亚硝酸盐，进而与食物中的仲胺类物质反应生成具有致癌性的亚硝胺，严重威胁人体健康。在氨氮污染治理方面，国内外已发展出多种技术，如物理法中的吹脱法，通过调节废水pH值并向其中通入气体，使氨分子挥发出来，达到去除氨氮的目的，该方法操作相对简单，但受温度、气液比等因素影响较大；化学法中的化学沉淀法，向废水中投加Mg²⁺与PO₄³⁻药剂，与NH₄⁺形成难溶的磷酸氨镁沉淀，从而实现脱氮，此方法适用于高浓度氨氮废水处理，且产生的沉淀可作为肥料再利用；生物法中的硝化反硝化技术应用广泛，利用微生物的作用将氨氮转化为氮气，包括好氧硝化和缺氧反硝化两个阶段，新型脱氮技术如短程硝化反硝化、同时硝化反硝化等，能够提高处理效率，减少处理成本。在氨氮污染模拟预测方面，国内外学者也开展了大量研究。通过建立数学模型，如基于质量守恒定律和化学反应动力学原理的水质模型，能够模拟氨氮在水体中的迁移、转化和扩散过程。然而，目前的模拟预测研究仍存在一些不足，对于复杂的水文地质条件和多种污染源相互作用的情况，模型的准确性和可靠性有待进一步提高。不同地区的水文地质条件和污染源特征差异较大，如何根据具体情况准确选取模型参数，使模型能够更真实地反映实际情况，仍是当前研究的难点之一。1.2.3数值模拟技术应用数值模拟技术在水文地质领域和氨氮污染模拟预测中发挥着重要作用。在水文地质领域，数值模拟技术可用于模拟地下水流场的变化，分析地下水的补给、径流和排泄过程。通过建立地下水流数值模型，能够预测不同开采方案下地下水位的变化趋势，为水资源的合理开发利用提供科学依据。例如，在某矿区的研究中，利用数值模拟技术构建了三维非稳定地下水流模型，分析了矿井疏干过程中地下水位的下降情况，以及对周边水环境的影响。在氨氮污染模拟预测方面，数值模拟技术能够模拟氨氮在地下水中的迁移转化过程，预测氨氮污染的扩散范围和程度。通过将水文地质模型与氨氮迁移转化模型相结合，可以考虑地下水流动、化学反应、吸附解吸等多种因素对氨氮迁移的影响。在研究某垃圾填埋场周边地下水氨氮污染时，运用数值模拟技术建立了地下水水流和氨氮迁移耦合模型，模拟了氨氮从垃圾填埋场渗漏进入地下水后的迁移过程，预测了污染的发展趋势，为制定污染防控措施提供了有力支持。在本研究中，数值模拟技术将用于构建旌阳区天元镇水源区的水文地质模型和氨氮迁移模型。通过对水文地质条件的准确刻画和模型参数的合理选取，模拟地下水流场和氨氮在地下水中的迁移转化过程，预测氨氮污染的发展趋势，为水源区的水资源保护和污染防治提供科学的决策依据。数值模拟技术还可以通过情景分析，评估不同防控措施对氨氮污染的控制效果，为选择最优的防控方案提供参考。1.3研究内容与方法1.3.1研究内容本研究聚焦旌阳区天元镇水源区，围绕水文地球化学特征、氨氮污染现状与迁移模拟预测展开深入研究，具体内容如下：水源区水文地球化学特征研究：对水源区的地下水及地表水进行系统采样，精确测定包括主要离子（如Ca²⁺、Mg²⁺、Na⁺、K⁺、HCO₃⁻、SO₄²⁻、Cl⁻等）、微量元素（如Fe、Mn、Cu、Zn等）以及同位素（如氢氧同位素、碳同位素等）在内的各项化学成分。运用数理统计方法，如描述性统计、相关性分析等，对数据进行处理与分析，以揭示水源区地下水和地表水的化学组成特征及分布规律。绘制Piper三线图、Schollerer图等多种水文地球化学图件，直观展示水化学类型的分布情况，深入分析控制水文地球化学特征的主要因素，包括岩石风化、阳离子交换、蒸发浓缩、人类活动等。氨氮污染现状调查：通过实地调查与监测，全面了解水源区内氨氮污染的来源，包括生活污水排放点的分布、排放量及污水中氨氮含量；工业企业的生产类型、废水排放情况及氨氮浓度；农业面源污染，如农田施肥种类、施肥量、畜禽养殖规模及粪便处理方式等对氨氮污染的贡献。对水源区内不同监测点的地下水和地表水进行氨氮含量的测定，分析氨氮在空间上的分布特征，确定氨氮污染的高值区和低值区，并结合地形地貌、地质条件以及污染源分布，探讨氨氮分布差异的原因。氨氮迁移模拟预测：基于对水源区水文地质条件的详细勘察和分析，构建准确的水文地质概念模型，明确含水层的分布、厚度、渗透系数等参数，以及边界条件和源汇项。选用合适的数值模拟软件，如VisualMODFLOW等，建立氨氮在地下水中的迁移转化模型，充分考虑对流、弥散、吸附解吸、生物化学作用等过程对氨氮迁移的影响。通过对模型进行参数率定和验证，确保模型能够准确反映实际情况。利用验证后的模型，对不同情景下氨氮污染的发展趋势进行模拟预测，如在现有污染源不变的情况下，预测未来一定时期内氨氮污染的扩散范围和浓度变化；分析在采取不同污染防控措施后，氨氮污染的改善情况，为制定合理的污染防治方案提供科学依据。1.3.2研究方法为实现上述研究内容，本研究将综合运用多种研究方法，具体如下：现场调查法：对旌阳区天元镇水源区进行全面的实地调查，详细记录地形地貌特征，包括山地、平原、河流、湖泊等的分布情况，以及它们对水流和污染物迁移的影响。收集区域内的气象数据，如降水量、蒸发量、气温、风速等，分析气象条件对水资源补给和污染物扩散的作用。深入调查地质构造，了解断层、褶皱等地质构造的分布及其对地下水流动的控制作用；详细勘查地层岩性，确定不同地层的岩石类型、渗透性等特征，为后续研究提供基础地质信息。全面排查水源区内的污染源，记录生活污水排放口、工业废水排放源、农业面源污染的分布位置、排放方式及污染强度等信息。样品分析法：在水源区内按照科学的布点原则，采集地下水和地表水样品。运用离子色谱仪精确测定样品中的主要离子浓度，确保分析结果的准确性和可靠性；使用电感耦合等离子体质谱仪（ICP-MS）测定微量元素含量，能够检测出极低浓度的微量元素；采用同位素质谱仪分析氢氧同位素、碳同位素等，通过同位素分析揭示地下水的补给来源、径流路径以及与其他水体的相互关系。通过这些先进的分析仪器和方法，获取全面、准确的样品化学组成信息。数据分析方法：运用描述性统计分析，计算各项化学指标的平均值、中位数、最大值、最小值、标准差等统计参数，对数据的集中趋势、离散程度等进行初步分析，了解数据的基本特征。开展相关性分析，研究不同化学指标之间的相关性，判断它们之间是否存在内在联系，以及这种联系的紧密程度，从而揭示水文地球化学过程和氨氮污染的影响因素。进行主成分分析（PCA）和因子分析（FA），将多个相关变量转化为少数几个综合因子，提取数据中的主要信息，简化数据结构，更清晰地识别影响水文地球化学特征和氨氮污染的主要因素。数值模拟方法：借助专业的数值模拟软件VisualMODFLOW，它是一款广泛应用于地下水模拟的软件，具有强大的功能和良好的用户界面。根据构建的水文地质概念模型，输入含水层参数、边界条件、源汇项等数据，建立三维地下水流模型，模拟地下水的流动状态，包括流速、流向、水位分布等。将氨氮迁移转化模型与地下水流模型进行耦合，考虑氨氮在地下水中的对流、弥散、吸附解吸、生物化学降解等过程，模拟氨氮在地下水中的迁移转化过程。通过不断调整模型参数，使模拟结果与实际监测数据相匹配，实现模型的率定和验证，确保模型的准确性和可靠性，为氨氮污染的预测和防治提供有力的技术支持。1.4技术路线本研究技术路线如图1-1所示：资料收集与现场调研：广泛收集旌阳区天元镇水源区的地形地貌、气象水文、地质构造、地层岩性等相关资料，全面了解区域的自然地理和地质背景。开展现场调查，详细记录污染源的分布、排放情况以及周边环境状况，为后续研究提供基础数据。样品采集与测试分析：在水源区内科学布置采样点，系统采集地下水和地表水样品。运用先进的分析仪器，准确测定样品中的主要离子、微量元素、同位素以及氨氮含量等指标，获取全面的样品化学组成信息。水文地球化学特征分析：对测试分析得到的数据进行数理统计分析，运用描述性统计、相关性分析等方法，初步了解数据特征和变量之间的关系。绘制Piper三线图、Schollerer图等水文地球化学图件，直观展示水化学类型和特征的分布情况，深入分析控制水文地球化学特征的主要因素。氨氮污染现状分析：根据氨氮含量的监测数据，分析氨氮在水源区内的空间分布特征，确定氨氮污染的高值区和低值区。结合地形地貌、地质条件以及污染源分布，探讨氨氮分布差异的原因，全面掌握氨氮污染现状。模型构建与模拟预测：基于对水文地质条件的深入理解，构建准确的水文地质概念模型。选用合适的数值模拟软件，建立三维地下水流模型和氨氮迁移转化模型。通过对模型进行参数率定和验证，确保模型的准确性和可靠性。利用验证后的模型，对不同情景下氨氮污染的发展趋势进行模拟预测。结果分析与对策建议：对模拟预测结果进行深入分析，评估氨氮污染对水源区的影响程度和发展趋势。结合分析结果，提出针对性的污染防治措施和水资源保护建议，为水源区的可持续发展提供科学依据。[此处插入技术路线图1-1，图中应清晰展示各环节之间的逻辑关系和数据流向，包括资料收集、样品采集与分析、数据分析、模型构建与模拟、结果分析与对策建议等步骤，各步骤之间用箭头连接表示先后顺序][此处插入技术路线图1-1，图中应清晰展示各环节之间的逻辑关系和数据流向，包括资料收集、样品采集与分析、数据分析、模型构建与模拟、结果分析与对策建议等步骤，各步骤之间用箭头连接表示先后顺序]二、研究区域概况2.1自然地理条件2.1.1地理位置与交通旌阳区天元镇水源区位于四川省德阳市旌阳区，地处东经104.3372°，北纬31.1281°，距离市区中心4公里，是进出德阳市市西大门的重要门户。其东与德阳市区以宝成铁路为界，南与旌阳区八角井镇为邻，西与广汉市金轮镇隔河相望，北同孝泉镇、扬嘉镇相连，是全省经济强县旌阳区、什邡市、绵竹市、广汉市的交汇地，区域优势显著。在交通方面，宝成铁路复线、德绵路、德茂路纵贯天元镇全境。宝成铁路作为连接中国西北与西南地区的重要交通干线，为天元镇与外界的物资运输和人员往来提供了便利。德绵路和德茂路则加强了天元镇与周边城镇的联系，促进了区域间的经济交流与合作。此外，市区长江西路延伸段、亭江路的复工建设工程，进一步完善了天元镇的交通网络。市一环路、二环路横穿而过，将整个天元镇纳入“10分钟经济圈”，极大地缩短了天元镇与德阳市区及其他区域的时空距离。便捷的交通条件不仅有利于人员和物资的流动，还为水源区的水资源开发、利用与保护提供了便利，同时也便于对水源区进行监测和管理，及时应对可能出现的水源污染等问题。2021年1月26日，德阳裕兴公共交通有限责任公司新开37路公交线路，从天元镇凯元检测站出发，经天鹤南路、东海路等路段，最终到达德中青衣江校区(德阳八中），解决了天元镇武庙社区、天鹤南路及周边区域市民公交出行需求，加强了天元镇与市区其他区域的公共交通联系，方便了居民的日常出行。这些交通线路的布局，使得天元镇在区域水资源调配和保护中具有重要的枢纽地位，能够更好地整合区域水资源，实现水资源的合理配置。2.1.2地形地貌天元镇水源区位于成都平原北部边缘，区内地貌大致以石亭江为界，整体呈现出平原地貌特征，海拔在460-520m之间。主要由绵远河、石亭江等河流冲洪积所形成的漫滩I级阶地和上更新统冰水-流水堆积的II级阶地组成，地势较为平坦开阔。平原地貌对水源区水文地质条件产生了多方面的影响。在地下水补给方面，平坦的地形有利于大气降水的入渗，使得降水能够较为均匀地渗透到地下，补充地下水。地表水体如河流、渠道等，也因地势平坦，水流速度相对缓慢，增加了地表水与地下水的水力联系，促进了地表水对地下水的补给。在地下水径流方面，平原地形使得地下水水力坡度较小，径流速度相对较慢，地下水在含水层中流动较为平稳。这也导致了地下水的排泄方式主要以侧向径流排泄和蒸发排泄为主，在地势低洼处，地下水可能会以泉的形式出露地表，形成排泄点。在河流阶地分布上，不同阶地的岩性和结构差异对地下水的赋存和运移有着重要影响。漫滩I级阶地多由全新统砂砾卵石组成，透水性强，富水性好，是地下水的主要储存和运移空间。上更新统冰水-流水堆积的II级阶地，其岩性相对较细，可能含有较多的粉质粘土、亚粘土等，透水性相对较弱，对地下水的径流产生一定的阻滞作用，使得地下水在不同阶地之间的流动存在差异，影响着地下水的分布和动态变化。这种地形地貌特征，决定了水源区地下水的形成、分布和运移规律，对后续的水文地球化学特征和氨氮污染研究具有重要的基础作用。2.1.3气象水文天元镇水源区属四川盆地亚热带湿润季风气候区，气候温和，降水充沛，四季分明。多年平均气温16.1℃，年平均无霜期271天，日照小时数适中。这种温和的气候条件，使得水源区的水体不易受到极端温度的影响，有利于水资源的稳定储存和循环。主要灾害性天气有低温寒潮、干旱、洪涝、冰雹和秋绵雨等，这些灾害性天气会对水资源的形成和分布产生显著影响。例如，干旱会导致降水量减少，地表水体干涸，地下水补给不足，水位下降；洪涝则会使地表径流量增大，可能引发河水倒灌，影响地下水水质，还可能造成水土流失，改变地表和地下的水文地质条件。在水文方面，旌阳区内河流主要为沱江上游绵远河，由北向南流经城区，属宽浅型复式河床，游荡性河流。河床宽200-400m，比降2‰左右，最大流量2235m³（1943年），常年洪水流量600m³左右，一般流量40-45m³，枯季流量仅5m³。西南边界为沱江支流石亭江，区内河段长，最大流量450m³，枯季流量较小。河流的流量变化和水位波动，直接影响着地表水与地下水的相互转化关系。在洪水期，河流水位高于地下水位，河水补给地下水；在枯水期，地下水位高于河流水位，地下水补给河水。区内渠道密集，为方便灌溉修建了众多灌溉渠。这些灌溉渠不仅为农业生产提供了水源，也在一定程度上改变了区域的水文循环。灌溉水的入渗增加了地下水的补给量，影响着地下水的水位和水质。此外，在丘陵地带分布着一些水库，水库的存在调节了地表径流，在丰水期蓄水，枯水期放水，对区域水资源的时空分布起到了调节作用，也会通过库塘入渗补给地下水，影响地下水的动态变化。气象水文条件共同作用，影响着水源区水资源的形成、分布和循环，是研究水文地球化学特征和氨氮污染的重要背景因素。2.2地质概况2.2.1地层岩性研究区地层岩性主要由第四系全新统（Q4）和上更新统（Q3）地层组成。第四系全新统（Q4）广泛分布于研究区，主要为河流冲洪积物，岩性以砂砾卵石层为主，夹有粉质粘土、粉细砂等薄层。这些砂砾卵石层颗粒较粗，分选性和磨圆度较好，孔隙较大，透水性强，是地下水的良好储存和运移介质，赋存着丰富的第四系松散岩类孔隙水。粉质粘土和粉细砂层的存在，在一定程度上影响着地下水的垂向渗透和水力联系，它们相对隔水，可形成局部的弱透水层，对地下水的径流和排泄产生一定的阻滞作用。上更新统（Q3）地层主要为冰水-流水堆积物，岩性为含泥砂砾卵石层、粉质粘土、粘土等互层。含泥砂砾卵石层具有一定的透水性，是地下水储存和运移的次要空间；粉质粘土和粘土层透水性差，为相对隔水层，对地下水的分布和运动起到隔水边界和层间阻隔的作用。在含水层结构上，不同地层岩性的组合形成了复杂的含水层系统。上部的全新统砂砾卵石层为强透水层，接受降水入渗、地表水补给等，地下水径流速度相对较快；下部的上更新统含泥砂砾卵石层与粉质粘土、粘土层互层，导致含水层的渗透性在垂向上存在差异，形成了多层结构的含水层，使得地下水在不同含水层之间的水力联系和运动规律变得复杂。这种地层岩性的分布特征对地下水的赋存和运移产生了重要影响。从赋存方面来看，砂砾卵石层的高孔隙度和强透水性，为地下水提供了充足的储存空间，使得研究区具有丰富的地下水资源。粉质粘土和粘土层的相对隔水作用，有助于地下水的保存，防止地下水过度渗漏和流失。在运移方面，强透水层为地下水的快速径流提供了通道，而弱透水层则改变了地下水的径流方向和速度，使得地下水在含水层中呈现出复杂的三维流动状态。地层岩性的差异还影响着地下水与外界的物质交换过程，不同岩性与地下水之间的化学反应不同，进而影响着地下水的化学成分和水文地球化学特征。2.2.2地质构造研究区位于扬子准地台四川台坳川西台陷的北部边缘，区域构造相对稳定，主要受龙门山断裂带和龙泉山断裂带的影响，但距离断裂带较远，地震活动相对较弱。区内主要发育有一些规模较小的褶皱和断裂构造，这些构造对地下水的流动和水文地球化学特征具有重要的控制作用。褶皱构造使得地层发生弯曲变形，改变了地层的原始产状。在背斜构造部位，地层向上拱起，含水层的厚度变薄，地下水的储存空间相对减小；同时，背斜顶部岩石受张力作用，裂隙发育，有利于地下水的径流和排泄，使得地下水水位相对较低，水力坡度较大，径流速度较快。向斜构造部位则相反，地层向下凹陷，含水层厚度相对较大，地下水储存空间丰富；向斜底部岩石受挤压作用，相对致密，裂隙不发育，地下水径流相对缓慢，水位较高，水力坡度较小。褶皱构造还可能导致不同岩性地层的接触关系发生变化，形成水力联系复杂的含水层系统，影响地下水的流动路径和补给、排泄条件。断裂构造是地下水流动的重要通道和边界。区内的断裂构造多为正断层和逆断层，断层破碎带岩石破碎，透水性强，成为地下水的优势运移通道。地下水可沿断层破碎带快速流动，使得断层附近的地下水水位、水质和水化学特征与周围地区存在明显差异。一些断裂构造还可能切断含水层，形成隔水边界，阻碍地下水的水平径流，导致地下水在断层一侧聚集，水位升高，另一侧则水位降低，影响地下水的分布格局。断裂构造还可能使深部地下水与浅部地下水发生水力联系，深部地下水的化学成分和温度等特征可能通过断裂带影响浅部地下水，改变浅部地下水的水文地球化学特征。地质构造对地下水流动和水文地球化学特征的控制作用，使得研究区地下水的分布和演化过程变得复杂。在进行水文地球化学研究和氨氮污染模拟预测时，必须充分考虑地质构造因素，准确把握地下水的流动规律和物质交换过程，才能更准确地揭示水文地球化学特征的形成机制，预测氨氮污染的迁移趋势。三、水文地质条件分析3.1地下水类型与含水岩组3.1.1第四系松散岩类孔隙水第四系松散岩类孔隙水在旌阳区天元镇水源区分布广泛，主要赋存于第四系全新统（Q4）和上更新统（Q3）地层的松散堆积物孔隙中。全新统地层主要为河流冲洪积物，岩性以砂砾卵石层为主，夹有粉质粘土、粉细砂等薄层。这些砂砾卵石层颗粒较粗，分选性和磨圆度较好，孔隙较大，透水性强，是第四系松散岩类孔隙水的主要储存和运移空间。粉质粘土和粉细砂层相对隔水，在一定程度上影响着孔隙水的垂向渗透和水力联系，可形成局部的弱透水层。上更新统地层为冰水-流水堆积物，岩性为含泥砂砾卵石层、粉质粘土、粘土等互层。含泥砂砾卵石层具有一定的透水性，也是孔隙水储存和运移的空间之一，但透水性相对全新统砂砾卵石层较弱。粉质粘土和粘土层透水性差，为相对隔水层，对孔隙水的分布和运动起到隔水边界和层间阻隔的作用。第四系松散岩类孔隙水主要接受大气降水入渗补给，降水通过地表土壤孔隙下渗，补充到孔隙含水层中。在河流附近，河水与孔隙水存在密切的水力联系，河水在洪水期水位较高时，会补给孔隙水；在枯水期，孔隙水则可能补给河水。此外，灌溉水的入渗也是孔隙水的重要补给来源之一，区内密集的灌溉渠使得大量灌溉水渗入地下，增加了孔隙水的补给量。孔隙水的排泄方式主要有蒸发排泄和侧向径流排泄。在地势低洼且潜水埋藏较浅的地区，孔隙水通过蒸发作用排泄到大气中，导致盐分在地表积累，可能引发土壤盐渍化等问题。侧向径流排泄则是孔隙水在水力坡度的作用下，向地势较低的方向流动，最终排泄到河流、湖泊等地表水体中，或者补给其他含水层。第四系松散岩类孔隙水与其他含水层之间存在一定的水力联系。在一些区域，孔隙水可通过弱透水层的越流补给深部的基岩裂隙水；在基岩出露区，孔隙水与基岩裂隙水可能通过断裂、裂隙等通道直接连通，发生水力交换。这种水力联系使得不同含水层之间的水量和水质相互影响，对整个地下水系统的动态变化和水文地球化学特征产生重要作用。3.1.2第四系松散岩类孔隙水富水性分布第四系松散岩类孔隙水的富水性在空间上存在明显的分布差异，主要受含水层岩性、厚度、地形地貌以及补给条件等多种因素的综合影响。在含水层岩性方面，以全新统砂砾卵石层为主的区域，由于其颗粒粗、孔隙大、透水性强，能够储存和传输大量的地下水，富水性较好。在河流主流线附近及漫滩I级阶地，砂砾卵石层厚度较大且连续性好，单井涌水量一般较大，可达[X]m³/d以上，属于富水区域。而在一些夹有较多粉质粘土、粉细砂层的地段，含水层的透水性变差，孔隙水的储存和运移受到阻碍，富水性相对较弱，单井涌水量可能在[X]m³/d以下。含水层厚度对富水性的影响也较为显著。在沉积作用较强、地层堆积较厚的区域，如河流的中下游冲积平原，第四系松散堆积物厚度较大，相应的孔隙含水层厚度也大，能够容纳更多的地下水，富水性好。而在地形较高、沉积作用较弱的地区，含水层厚度较薄，富水性较差。在研究区的边缘地带，由于地势相对较高，第四系地层厚度较薄，孔隙水的富水性明显低于中心的平原区域。地形地貌对孔隙水富水性分布起着重要的控制作用。平原地区地势平坦，有利于大气降水的入渗和地表水体的补给，且水力坡度较小，地下水径流缓慢，有利于地下水的储存，使得平原区的孔隙水富水性普遍较好。在河流阶地分布区，不同阶地的富水性存在差异。漫滩I级阶地靠近河流，接受河水补给充足，含水层岩性好，富水性强；上更新统冰水-流水堆积的II级阶地，岩性相对较细，透水性略差，富水性次之。在丘陵地区，地形起伏较大，降水易形成地表径流快速流失，入渗补给量少，且含水层厚度变化较大，富水性相对较差。补给条件是影响孔隙水富水性的关键因素之一。降水充沛、地表水丰富且补给条件良好的区域，孔隙水能够得到充足的水源补充，富水性较好。研究区属四川盆地亚热带湿润季风气候区，降水充沛，为孔隙水提供了丰富的补给来源。河流、灌溉渠等地表水体与孔隙水之间的水力联系密切，也增加了孔隙水的补给量。在一些远离河流和灌溉水源的地区，补给量相对较少，富水性较弱。3.1.3基岩裂隙水基岩裂隙水主要赋存于研究区基岩的裂隙之中，其赋存条件与基岩的岩性、地质构造以及风化作用密切相关。区内基岩主要为各类沉积岩和变质岩，这些岩石在长期的地质作用过程中，受到构造应力、风化作用等影响，产生了各种裂隙，为基岩裂隙水的赋存提供了空间。构造裂隙是基岩裂隙水赋存的重要场所。区内的褶皱和断裂构造使得基岩产生了大量的构造裂隙。在褶皱的轴部和翼部，由于岩石受力变形，裂隙较为发育，有利于基岩裂隙水的储存和运移。背斜顶部岩石受张力作用，裂隙张开度较大，连通性较好，是基岩裂隙水的良好储存空间；向斜部位虽然岩石相对致密，但在翼部也存在一定数量的裂隙，能够赋存一定量的地下水。断裂构造，尤其是导水断裂，是基岩裂隙水的优势运移通道。断层破碎带岩石破碎，裂隙密集，透水性强，地下水可沿断层快速流动，使得断层附近的基岩裂隙水水位、水质和水化学特征与周围地区存在明显差异。风化裂隙也是基岩裂隙水的重要赋存空间。地表基岩在风化作用下，形成了风化裂隙层。风化裂隙在地表呈壳状分布，厚度一般为[X]m，这些裂隙密集均匀，相互连通，构成了风化裂隙含水系统。风化裂隙水一般为潜水，受大气降水补给，其动态变化与降水密切相关。在风化作用强烈的地区，风化裂隙发育良好，基岩裂隙水相对丰富；而在风化作用较弱的地区，风化裂隙较少，基岩裂隙水含量较低。基岩裂隙水与第四系松散岩类孔隙水之间存在一定的水力联系。在基岩出露区或基岩与第四系地层接触部位，孔隙水可通过裂隙、断裂等通道补给基岩裂隙水；反之，在一定条件下，基岩裂隙水也可向上补给孔隙水。这种水力联系使得两个含水层之间的水量和水质发生交换，影响着整个地下水系统的动态平衡和水文地球化学特征。基岩裂隙水的动态变化相对较为复杂，除了受降水补给影响外，还受到地质构造、含水层渗透性以及与其他含水层水力联系等多种因素的制约。在导水断裂附近，基岩裂隙水的水位变化可能较为迅速，且受区域地下水流动系统的影响较大；而在风化裂隙发育区，基岩裂隙水的水位变化则主要受降水的季节性变化影响。三、水文地质条件分析3.2地下水补给径流排泄特征3.2.1补给方式大气降水补给：大气降水是旌阳区天元镇水源区地下水的重要补给来源之一。研究区属四川盆地亚热带湿润季风气候区，多年平均降水量较为丰富。降水通过地表土壤孔隙、岩石裂隙等途径渗入地下，补充地下水。降水对地下水的补给量受多种因素影响，降水强度是重要因素之一，当降水强度较小时，雨水有足够的时间下渗，能够有效地补给地下水；而当降水强度过大时，地表容易形成径流，大部分雨水会随径流流失，入渗补给地下水的量相对减少。地形地貌也对降水补给有显著影响，在平原地区，地势平坦，地表径流流速较慢，有利于降水的入渗，使得大气降水能够充分补给地下水；而在丘陵或山区，地形起伏较大，降水易形成地表径流快速流走，入渗补给量相对较少。包气带岩性同样影响着降水补给，包气带岩性为砂砾石等透水性强的物质时，降水能够迅速下渗，补给地下水；若包气带为粉质粘土、粘土等透水性差的物质，降水下渗会受到阻碍，补给量会相应减少。降水对地下水的补给具有明显的季节性变化，在雨季，降水量大，地下水补给量增加，水位上升；旱季时，降水量少，地下水补给量减少，水位下降。通过对研究区降水与地下水位动态变化的监测分析，发现地下水位的变化与降水过程具有一定的相关性，在降水集中的时段后，地下水位往往会出现明显的上升趋势。地表水补给：地表水与地下水之间存在密切的水力联系，地表水是地下水的重要补给源之一。研究区内河流主要有绵远河、石亭江等，这些河流流量较大，且在洪水期水位较高，河水能够通过河床渗漏、侧向渗透等方式补给地下水。河流对地下水的补给强度受多种因素制约，河床岩性是关键因素之一，当河床岩性为砂砾石等透水性强的物质时，河水渗漏量大，对地下水的补给作用显著；若河床为粉质粘土、粘土等相对隔水层，河水渗漏量小，补给作用较弱。河流水位与地下水位的高差也影响着补给强度，高差越大，水力梯度越大，河水补给地下水的速度越快，补给量也越大。在洪水期，河流水位大幅上升，高于地下水位，河水快速补给地下水，使得地下水水位迅速抬升；枯水期时，河流水位下降，可能低于地下水位，此时地下水反而会补给河水。区内渠道密集，灌溉水的入渗也是地下水的重要补给方式。在农业灌溉期间，大量的灌溉水通过渠道渗漏和田间入渗进入地下，增加了地下水的补给量。对灌溉前后地下水位的监测表明，灌溉后地下水位明显上升，说明灌溉水对地下水的补给作用明显。此外，一些水库、池塘等地表水体也会在一定条件下补给地下水，调节地下水的水量平衡。侧向径流补给：在研究区的边界地带或不同含水层之间，地下水存在侧向径流补给的情况。当相邻区域的地下水位高于研究区时，地下水会在水力梯度的作用下，通过含水层的孔隙、裂隙等通道，以侧向径流的方式流入研究区，补给研究区的地下水。在研究区与周边山区的交界处，山区的地下水水位相对较高，会向地势较低的天元镇水源区侧向径流补给。这种侧向径流补给受含水层的渗透性和水力梯度的控制，含水层渗透性越好，水力梯度越大，侧向径流补给量就越大。在一些断层或裂隙发育的区域，地下水的侧向径流补给可能更为明显，因为这些区域的岩石破碎，透水性强，为地下水的流动提供了良好的通道。侧向径流补给对研究区地下水的水量和水质都有一定的影响，它不仅增加了研究区地下水的补给量，还可能带来不同化学组成的地下水，改变研究区地下水的水质。3.2.2径流与天然排泄方式径流方向与速度：地下水的径流方向主要受地形地貌和地质构造的控制。在旌阳区天元镇水源区，总体上地下水由地势较高的区域向地势较低的区域径流。在平原地区，由于地势相对平坦，地下水水力坡度较小，径流方向较为平缓，大致与地形坡降方向一致。在河流阶地分布区，地下水的径流方向可能会受到阶地地形和含水层结构的影响，在漫滩I级阶地，地下水可能会沿着阶地的倾斜方向向河流或更低地势处径流；上更新统II级阶地，由于岩性和结构的差异，地下水的径流方向可能会发生一定的变化。在基岩出露区，地质构造对地下水径流方向起着重要的控制作用。断层、裂隙等构造会改变地下水的流动路径，使地下水沿着构造通道流动。在一些背斜构造部位，地下水可能会向背斜顶部汇聚，然后通过断层或裂隙向其他区域径流；向斜构造部位则相反，地下水会在向斜底部汇聚，并向周边区域径流。地下水的径流速度受多种因素影响，含水层的渗透系数是关键因素之一。第四系松散岩类孔隙含水层中，以砂砾卵石层为主的区域，渗透系数大，地下水径流速度相对较快；而在夹有较多粉质粘土、粉细砂层的地段，渗透系数小，径流速度较慢。基岩裂隙含水层中，裂隙发育程度好、连通性强的区域，地下水径流速度较快；裂隙不发育或被充填的区域，径流速度较慢。水力坡度也对径流速度有重要影响，水力坡度越大，地下水的驱动力越大，径流速度越快。在靠近河流或地形变化较大的区域，水力坡度相对较大，地下水径流速度可能会加快。天然排泄方式与途径：地下水的天然排泄方式主要有蒸发排泄和侧向径流排泄。在地势低洼且潜水埋藏较浅的地区，地下水通过蒸发作用排泄到大气中。蒸发排泄受气候条件和潜水埋深的影响，在炎热干旱的季节，蒸发量大，地下水的蒸发排泄量也会相应增加；潜水埋深越浅，越容易发生蒸发排泄。在研究区的一些农田、湿地等区域，由于潜水埋藏较浅，在夏季高温时段，地下水的蒸发排泄较为明显，可能导致土壤盐分积累，影响土壤质量和农作物生长。侧向径流排泄是地下水在水力坡度的作用下，向地势较低的方向流动，最终排泄到河流、湖泊等地表水体中，或者补给其他含水层。在研究区的边界地带，地下水会向周边区域侧向径流排泄，与相邻区域的地下水进行水量和水质的交换。在河流附近，地下水会侧向径流补给河流，维持河流的水量平衡。此外，在一些岩溶发育地区，地下水还可能通过岩溶管道排泄到溶洞、暗河等岩溶系统中。在研究区若存在岩溶地质条件，岩溶排泄也是地下水的一种重要天然排泄途径，岩溶管道的发育和连通性对地下水的排泄速度和排泄量有重要影响。3.3地下水动态特征地下水动态是指地下水的水位、水量、水温、水化学成分在各种因素综合影响下，随时间作有规律的变化。研究区地下水动态特征对于深入了解地下水的形成、循环规律以及水资源的合理开发利用具有重要意义。3.3.1水位动态变化通过对研究区多个地下水监测点的长期水位监测数据进行分析，发现地下水水位呈现出明显的季节性变化特征。在雨季（5月-9月），由于大气降水充沛，地下水补给量增加，水位迅速上升。根据监测数据，雨季期间地下水位平均上升幅度可达[X]m，其中部分靠近河流或降水入渗条件较好的区域，水位上升幅度甚至超过[X]m。在旱季（10月-次年4月），降水减少，地下水补给量相应减少，同时蒸发排泄和侧向径流排泄持续进行，导致地下水位逐渐下降。旱季期间地下水位平均下降幅度约为[X]m。这种季节性的水位变化与当地的气候条件密切相关，大气降水的季节性分布是控制地下水位动态变化的主要因素之一。除了季节性变化外，地下水水位还存在多年变化趋势。从较长时间序列的监测数据来看，过去[X]年间，研究区地下水位整体呈现出波动下降的趋势。这可能是由于多种因素共同作用的结果，随着区域经济的发展，工农业用水量不断增加，地下水开采量逐渐增大，导致地下水的排泄量超过补给量，从而引起地下水位下降。研究区周边的城市化进程加快，地表硬化面积增加，降水入渗量减少，也对地下水补给产生了不利影响，进一步加剧了地下水位的下降趋势。不同含水层的水位动态变化也存在一定差异。第四系松散岩类孔隙水水位变化相对较为敏感，对降水和地表水补给的响应迅速，水位波动幅度较大。在降水集中的时段，孔隙水水位能够在短时间内快速上升；而在旱季，水位下降也较为明显。基岩裂隙水水位变化相对较为平缓，其水位动态除了受降水影响外，还受到地质构造、裂隙发育程度以及与其他含水层水力联系等多种因素的制约。在一些构造复杂、裂隙连通性较差的区域，基岩裂隙水水位变化相对较小，且与降水的相关性不如孔隙水明显。3.3.2水量动态变化地下水水量动态变化主要受补给和排泄因素的影响。在补给方面，大气降水和地表水的补给量直接影响着地下水的水量。如前所述，雨季降水补给增加，使得地下水水量相应增加；旱季补给减少，水量随之减少。地表水与地下水的水力联系也对水量动态产生重要影响，在河流丰水期，河水补给地下水，增加了地下水的水量；枯水期则相反，地下水补给河水，导致地下水水量减少。在排泄方面，除了自然排泄方式（如蒸发排泄和侧向径流排泄）外，人类活动对地下水排泄的影响日益显著。随着研究区工农业的发展，地下水开采量不断增加，成为地下水排泄的重要途径之一。大量抽取地下水用于灌溉、工业生产和居民生活，导致地下水水量减少，尤其是在农业灌溉季节，地下水开采量大幅增加，对地下水水量动态产生了明显的影响。不合理的地下水开采还可能引发一系列环境问题，如地面沉降、地裂缝等，进一步影响地下水的赋存和水量动态。通过对不同年份地下水开采量和水位变化的对比分析发现，地下水开采量与地下水位之间存在明显的负相关关系。当开采量增加时，地下水位下降明显；减少开采量后，地下水位有一定程度的回升。在某一时期，由于农业灌溉需求增加，地下水开采量大幅上升，导致研究区部分区域地下水位在短时间内急剧下降；而在采取节水措施，减少地下水开采量后，地下水位逐渐趋于稳定，并在一定程度上有所回升。这表明合理控制地下水开采量对于维持地下水水量动态平衡和保护地下水资源具有重要意义。四、水文地球化学特征及成因分析4.1数据采集与分析方法为深入研究旌阳区天元镇水源区的水文地球化学特征，本研究开展了系统的数据采集工作。在水源区内，综合考虑地形地貌、地层岩性、地下水类型以及污染源分布等因素，科学合理地布置了[X]个地下水采样点和[X]个地表水采样点。其中，地下水采样点涵盖了第四系松散岩类孔隙水和基岩裂隙水分布区域，确保能够全面反映不同类型地下水的水文地球化学特征；地表水采样点则主要分布在绵远河、石亭江等主要河流以及区内的灌溉渠和水库等地表水体中。水样采集严格按照相关标准和规范进行操作，以保证样品的代表性和准确性。对于地下水采样，使用专业的地下水采样设备，通过钻孔或井管采集水样，确保采集到的水样能够真实反映含水层的水质情况。在采样前，先对采样设备进行清洗和消毒，避免设备本身对水样造成污染。地表水采样时，根据水体的深度和流速，选择合适的采样方法和工具，如使用聚乙烯塑料桶采集表层水样，对于较深水体则采用直立式采水器采集不同深度的水样。在样品保存与运输方面，采取了严格的措施。采集后的水样立即加入相应的保存剂，以防止水样中的化学成分发生变化。对于测定金属离子的水样，加入硝酸使其pH值小于2，以保持金属离子的溶解状态；对于测定氨氮、硝酸盐等的水样，加入硫酸调节pH值至2左右，抑制微生物的生长和化学反应的进行。水样保存后，尽快将其放入低温冷藏箱中，在4℃左右的温度下运输至实验室进行分析测试，确保水样在运输过程中的稳定性。实验室检测分析运用了先进的仪器和方法，以获取准确的水样化学组成信息。采用离子色谱仪测定水样中的主要阴离子（如HCO₃⁻、SO₄²⁻、Cl⁻等）和阳离子（如Ca²⁺、Mg²⁺、Na⁺、K⁺等）浓度，离子色谱仪具有分析速度快、灵敏度高、分离效果好等优点，能够准确测定水样中各种离子的含量。利用电感耦合等离子体质谱仪（ICP-MS）测定微量元素（如Fe、Mn、Cu、Zn等）含量，ICP-MS能够检测出极低浓度的微量元素，具有极高的灵敏度和准确性。采用同位素质谱仪分析水样中的氢氧同位素（如δD、δ¹⁸O）、碳同位素（如δ¹³C）等，通过同位素分析可以揭示地下水的补给来源、径流路径以及与其他水体的相互关系。在数据处理与分析过程中，运用多种数理统计方法对检测数据进行深入分析。首先进行描述性统计分析，计算各项化学指标的平均值、中位数、最大值、最小值、标准差等统计参数，以了解数据的集中趋势、离散程度和分布特征。通过描述性统计，可以初步掌握水源区地下水和地表水的化学组成概况，发现数据中的异常值和变化趋势。开展相关性分析，研究不同化学指标之间的相关性，判断它们之间是否存在内在联系以及这种联系的紧密程度。通过相关性分析，可以揭示水文地球化学过程中的相互作用机制，确定影响水文地球化学特征的主要因素。在研究中发现，地下水中的Ca²⁺和HCO₃⁻含量存在显著的正相关关系，这表明它们可能受到相同的水文地球化学过程控制，如碳酸盐岩的溶解作用。运用主成分分析（PCA）和因子分析（FA）等多元统计方法，将多个相关变量转化为少数几个综合因子，提取数据中的主要信息，简化数据结构，更清晰地识别影响水文地球化学特征的主要因素。主成分分析通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量，即主成分，这些主成分能够反映原始数据的主要特征。因子分析则是从多个变量中提取出公共因子，每个公共因子代表了原始变量之间的一种潜在关系，通过对公共因子的分析可以深入了解数据的内在结构和影响因素。四、水文地球化学特征及成因分析4.2水文地球化学特征分析4.2.1描述性统计分析对采集的水样中主要离子浓度进行描述性统计分析，结果如表4-1所示。地下水中Ca²⁺浓度范围为[Ca²⁺最小值]-[Ca²⁺最大值]mg/L，平均值为[Ca²⁺平均值]mg/L，表明地下水中Ca²⁺含量存在一定的变化范围，且平均值反映了其总体水平。Mg²⁺浓度范围为[Mg²⁺最小值]-[Mg²⁺最大值]mg/L，平均值为[Mg²⁺平均值]mg/L，说明Mg²⁺浓度也有一定波动。Na⁺浓度范围较广，为[Na⁺最小值]-[Na⁺最大值]mg/L，平均值为[Na⁺平均值]mg/L，反映出地下水中Na⁺的来源和影响因素较为复杂。K⁺浓度相对较低，范围是[K⁺最小值]-[K⁺最大值]mg/L，平均值为[K⁺平均值]mg/L。在阴离子方面，HCO₃⁻浓度范围为[HCO₃⁻最小值]-[HCO₃⁻最大值]mg/L，平均值为[HCO₃⁻平均值]mg/L，是地下水中的主要阴离子之一。SO₄²⁻浓度范围为[SO₄²⁻最小值]-[SO₄²⁻最大值]mg/L，平均值为[SO₄²⁻平均值]mg/L，其含量变化可能与岩石矿物的溶解、人类活动等因素有关。Cl⁻浓度范围为[Cl⁻最小值]-[Cl⁻最大值]mg/L，平均值为[Cl⁻平均值]mg/L，Cl⁻浓度的变化可能受到海水入侵、工业污染、生活污水排放等因素的影响。地表水的离子浓度也呈现出一定的变化特征。Ca²⁺、Mg²⁺、Na⁺、K⁺等阳离子以及HCO₃⁻、SO₄²⁻、Cl⁻等阴离子的浓度范围和平均值与地下水存在一定差异。地表水的Ca²⁺浓度平均值相对地下水略低，这可能是由于地表水与大气、土壤等的接触更为频繁，其化学成分受到更多外界因素的影响。地表水的SO₄²⁻浓度平均值有时会高于地下水，可能是因为地表受到更多工业排放、农业面源污染等影响，导致SO₄²⁻含量增加。通过对主要离子浓度的描述性统计分析，可以初步了解水源区地下水和地表水的化学组成特征和数据分布规律。浓度范围反映了离子含量的变化幅度，标准差则衡量了数据的离散程度，较大的标准差表明离子浓度在不同采样点之间的差异较大，可能受到不同的地质条件、补给来源、人类活动等因素的影响。平均值则代表了离子浓度的总体水平，为后续分析提供了基础数据。对数据的分布形态进行分析，发现部分离子浓度可能呈现正态分布或对数正态分布，这有助于进一步选择合适的统计方法和模型来深入研究水文地球化学特征。[此处插入表4-1，表中应包含地下水和地表水的主要离子（Ca²⁺、Mg²⁺、Na⁺、K⁺、HCO₃⁻、SO₄²⁻、Cl⁻等）浓度的最小值、最大值、平均值、中位数、标准差等统计参数]4.2.2地下水化学类型运用Piper三线图来确定地下水化学类型，Piper三线图由两个等边三角形和一个菱形组成，能够直观地展示地下水中主要离子的相对含量和水化学类型特征。通过计算水样中Ca²⁺、Mg²⁺、Na⁺、K⁺、HCO₃⁻、SO₄²⁻、Cl⁻等主要离子的毫克当量百分数，将其投点到Piper三线图中，结果如图4-1所示。从图中可以看出，研究区地下水化学类型主要为Ca-HCO₃型和Ca・Mg-HCO₃型。在Ca-HCO₃型区域，Ca²⁺在阳离子中占主导地位，HCO₃⁻在阴离子中占主导地位，这表明该类型地下水主要受碳酸盐岩的溶解作用影响。在研究区，地层中广泛分布的碳酸盐岩（如石灰岩、白云岩等）在地下水的溶蚀作用下，释放出Ca²⁺和HCO₃⁻，使得地下水中这两种离子含量较高，从而形成Ca-HCO₃型水。Ca・Mg-HCO₃型水则是Ca²⁺和Mg²⁺在阳离子中均占有一定比例，HCO₃⁻仍为主要阴离子。这种类型的形成除了碳酸盐岩的溶解作用外，还可能与白云岩等含镁碳酸盐岩的溶解有关，使得地下水中Mg²⁺含量增加。部分采样点的水化学类型呈现出过渡特征，可能是由于不同来源的地下水混合，或者受到多种水文地球化学过程的共同影响。在靠近河流的采样点，可能受到河水补给的影响，导致水化学类型发生一定变化；在人类活动密集区域，工业废水、生活污水的排放可能改变地下水的化学成分，进而影响水化学类型。[此处插入Piper三线图4-1，图中应清晰标注各采样点的位置和对应的水化学类型，三角形和菱形区域应标注明确，坐标轴刻度清晰可读]4.2.3主要因子空间分布特征利用克里金插值法绘制主要离子浓度的空间分布图，以分析其空间变化规律。图4-2为Ca²⁺浓度空间分布图，从图中可以看出，Ca²⁺浓度在研究区呈现出一定的空间差异。在北部和西部部分区域，Ca²⁺浓度相对较高，可能是因为这些区域地层中碳酸盐岩含量较高，且地下水径流相对缓慢，有利于碳酸盐岩的溶解，使得Ca²⁺不断富集。在东部和南部的一些区域，Ca²⁺浓度相对较低，可能与这些区域的地层岩性、地下水补给来源以及人类活动影响有关。东部靠近河流，河水的稀释作用可能导致地下水中Ca²⁺浓度降低；南部可能受到人类活动排放的影响，引入了其他化学成分，改变了地下水中离子的相对含量。图4-3为SO₄²⁻浓度空间分布图，SO₄²⁻浓度在空间上也有明显的变化。在中部和南部的部分区域，SO₄²⁻浓度较高，这可能与工业活动排放的含硫废水、农业面源污染中含硫肥料的使用以及地层中含硫矿物的溶解有关。在北部和西部的一些区域，SO₄²⁻浓度相对较低，可能是因为这些区域受工业和农业污染影响较小，且地层中含硫矿物含量较少。[此处插入Ca²⁺浓度空间分布图4-2和SO₄²⁻浓度空间分布图4-3，图中应标注清楚地理方位、采样点位置以及浓度等值线，等值线间隔合理，能够清晰展示离子浓度的空间变化趋势]通过对主要离子浓度空间分布特征的分析，可以了解不同离子在研究区内的富集和扩散规律，以及其与地质条件、人类活动等因素的关系。这对于深入理解水文地球化学过程、识别潜在的污染源以及评估地下水质量的空间变化具有重要意义。空间分布特征还可以为进一步的数值模拟研究提供基础数据，帮助准确刻画地下水中化学成分的迁移转化过程。4.3成因分析4.3.1控制因素分析地质构造对水文地球化学特征起着重要的控制作用。研究区位于扬子准地台四川台坳川西台陷的北部边缘，区内发育的褶皱和断裂构造影响着地下水的流动路径和水-岩相互作用过程。在褶皱构造部位，背斜顶部岩石受张力作用，裂隙发育，地下水径流速度较快，水-岩相互作用时间相对较短，导致地下水中的化学成分相对简单。向斜部位岩石受挤压作用，相对致密，地下水径流缓慢，水-岩相互作用时间长，使得地下水中的溶解物质增多，化学成分更为复杂。断裂构造作为地下水的优势运移通道，会使不同含水层之间的水发生混合，改变地下水的化学组成。在一些断层附近，地下水中的离子浓度和水化学类型可能会发生明显变化，这是由于断层沟通了不同的含水层，使得具有不同化学组成的地下水相互混合。岩石类型是控制水文地球化学特征的关键因素之一。研究区主要地层岩性为第四系全新统和上更新统的松散堆积物以及基岩，不同岩石类型与地下水之间的化学反应不同，从而影响地下水的化学成分。第四系松散堆积物中的砂砾卵石层，透水性强，与地下水的接触面积大，有利于地下水与岩石颗粒表面的物质发生交换。地层中的碳酸盐岩（如石灰岩、白云岩等），在地下水的溶蚀作用下，会释放出Ca²⁺、Mg²⁺、HCO₃⁻等离子，使得地下水中这些离子含量升高，形成以Ca-HCO₃型或Ca・Mg-HCO₃型为主的水化学类型。基岩中的矿物成分也会对地下水化学成分产生影响，如基岩中的长石、云母等矿物在风化作用下，会释放出K⁺、Na⁺等阳离子，增加地下水中这些离子的含量。人类活动对水文地球化学特征的影响日益显著。随着研究区经济的发展，工业废水、生活污水的排放以及农业面源污染等，都对地下水和地表水的化学成分产生了重要影响。工业活动中，一些企业排放的含重金属、酸碱废水等，会导致地下水中重金属离子（如Cu²⁺、Pb²⁺、Zn²⁺等）浓度升高，pH值发生变化。生活污水中含有大量的有机物、氮、磷等污染物，这些污染物进入水体后，会消耗水中的溶解氧，促进微生物的生长繁殖，导致水体富营养化，同时也会改变地下水中的离子组成。农业面源污染方面，大量使用的化肥、农药等，通过地表径流和淋溶作用进入水体，使得地下水中的NO₃⁻、SO₄²⁻、Cl⁻等阴离子浓度增加，以及农药残留对地下水水质产生潜在危害。在一些农业灌溉区，长期不合理的灌溉方式导致地下水位上升，盐分在地表积累，引发土壤盐渍化，进而影响地下水的化学组成。4.3.2相关性分析对地下水中主要离子进行相关性分析，结果如表4-2所示。从表中可以看出，Ca²⁺与HCO₃⁻之间存在显著的正相关关系，相关系数达到[具体数值]，这表明它们可能主要来源于碳酸盐岩的溶解作用。在地下水的溶蚀过程中，碳酸盐岩（如CaCO₃）与水中的CO₂和H₂O发生反应，生成Ca²⁺和HCO₃⁻，使得这两种离子在地下水中的含量密切相关。Na⁺与Cl⁻之间也呈现出一定的正相关关系，相关系数为[具体数值]，说明它们可能有共同的来源或受到相似的水文地球化学过程影响。Na⁺和Cl⁻可能来源于岩盐的溶解，或者是受到海水入侵、工业废水排放等因素的影响，使得它们在地下水中的含量同步变化。Mg²⁺与Ca²⁺之间存在一定的相关性，相关系数为[具体数值]，这可能是因为它们在岩石矿物中的存在形式和溶解机制有相似之处，都与碳酸盐岩等矿物的溶解有关。而SO₄²⁻与其他主要离子之间的相关性相对较弱，这表明SO₄²⁺的来源可能较为复杂，除了受到岩石矿物溶解的影响外，还可能受到工业活动、农业面源污染等多种因素的影响。在工业生产中，一些含硫矿物的开采和加工会产生大量的含SO₄²⁻废水，这些废水的排放会增加地下水中SO₄²⁻的含量；农业生产中，含硫肥料的使用也会导致SO₄²⁻进入水体。[此处插入表4-2，表中应包含主要离子（Ca²⁺、Mg²⁺、Na⁺、K⁺、HCO₃⁻、SO₄²⁻、Cl⁻等）之间的相关系数，以及对应的显著性水平]通过相关性分析，可以初步揭示地下水中离子之间的内在联系，为进一步探讨水文地球化学过程和主要离子来源提供依据。但相关性分析只能反映变量之间的线性关系，对于一些复杂的非线性关系，还需要结合其他分析方法进行深入研究。4.3.3主要离子来源分析运用多种方法确定主要离子来源，进一步阐述其形成机制。通过矿物平衡计算法，对研究区地层中的主要矿物进行分析，计算矿物溶解后可能释放到地下水中的离子含量。结果表明，碳酸盐岩（如方解石CaCO₃、白云石CaMg(CO₃)₂等）的溶解是地下水中Ca²⁺、Mg²⁺、HCO₃⁻的主要来源。以方解石溶解为例，其化学反应方程式为：CaCO₃+CO₂+H₂O⇌Ca²⁺+2HCO₃⁻，在研究区的地下水环境中，由于CO₂的存在，促进了方解石的溶解，使得地下水中Ca²⁺和HCO₃⁻含量升高。岩石的风化作用也是主要离子的重要来源。研究区的岩石在长期的风化作用下，矿物逐渐分解，释放出各种离子。长石类矿物（如钾长石KAlSi₃O₈、钠长石NaAlSi₃O₈等）在风化过程中，会发生水解反应，释放出K⁺、Na⁺等阳离子。以钾长石水解为例：2KAlSi₃O₈+2H₂O+2CO₂→2K⁺+2HCO₃⁻+H₂Al₂Si₂O₈+4SiO₂，从而增加了地下水中K⁺的含量。人类活动对主要离子来源的影响也不容忽视。工业废水排放中含有大量的重金属离子和酸根离子，如电镀行业排放的废水中含有Cu²⁺、Zn²⁺等重金属离子，化工行业排放的废水中含有SO₄²⁻、Cl⁻等酸根离子，这些废水未经有效处理直接排放，会导致地下水中相应离子含量升高。生活污水排放中含有大量的有机物、氮、磷等污染物，其中的氮元素在微生物的作用下，会转化为氨氮、硝酸盐等形式，增加地下水中NO₃⁻等离子的含量。农业面源污染方面，化肥的大量使用使得地下水中的NO₃⁻、SO₄²⁻、Cl⁻等阴离子含量显著增加。在一些蔬菜种植区，由于过量施用氮肥，导致地下水中NO₃⁻浓度严重超标，对地下水水质造成了严重威胁。综合以上分析，研究区地下水中主要离子的来源是多种因素共同作用的结果，地质因素（如岩石类型、风化作用）是主要离子的自然来源基础，而人类活动则在一定程度上改变了离子的组成和含量，对水文地球化学特征产生了重要影响。在进行水资源保护和管理时，需要充分考虑这些因素，采取有效的措施减少人类活动对地下水水质的负面影响。4.4盐度危害分析水源区盐度对水资源利用和生态环境存在潜在危害，需进行深入分析评估。过高的盐度会影响水源区地下水和地表水的水质，使其不适宜作为生活饮用水和工农业用水。当盐度超过一定标准时，用于生活饮用可能会对人体健康产生不良影响，如导致心血管疾病、肾脏负担加重等问题。在农业灌溉方面，高盐度的水会使土壤盐分积累，造成土壤盐渍化，影响土壤的物理化学性质，降低土壤肥力，抑制农作物生长，导致农作物减产甚至绝收。在工业生产中，高盐度的水可能会对生产设备造成腐蚀，降低设备使用寿命，增加生产成本，影响产品质量。对生态环境而言，盐度变化会对水源区的水生生态系统产生显著影响。过高的盐度会改变水体的渗透压，影响水生生物的生理功能和生存环境，导致水生生物种类和数量减少，破坏水生生态系统的平衡。在一些盐度较高的水体中，可能只有耐盐性较强的生物能够生存，而许多不耐盐的生物则会逐渐消失，使得生物多样性降低。盐度的变化还可能引发水体富营养化问题，进一步恶化生态环境。通过对水源区水样的盐度分析，结合相关标准和实际情况，评估盐度危害程度。研究区部分区域的地下水和地表水盐度已超出适宜范围，存在一定的盐度危害风险。在靠近工业排放源和农业面源污染严重的区域，盐度相对较高，危害程度较大；而在远离污染源、受人类活动影响较小的区域，盐度相对较低，危害程度较小。为降低盐度危害，应采取相应的应对建议。加强对工业废水和生活污水的处理，确保达标排放，减少盐分等污染物的排放。在农业生产中，推广合理的灌溉技术和施肥方式，减少化肥的使用量，避免因不合理灌溉和施肥导致盐分在土壤和水体中积累。建立完善的水源区水质监测体系，实时监测盐度等水质指标的变化，及时发现盐度异常情况，并采取有效的治理措施。还可以通过工程措施，如建设排盐渠道、改良土壤等，降低土壤和水体中的盐分含量，减轻盐度危害。4.5污染源分布通过详细的实地调查与资料收集，识别出旌阳区天元镇水源区存在多种潜在污染源，这些污染源的分布具有一定特征，对水源区的水质产生了不同程度的影响。工业污染源方面，区内分布着一些小型制造业企业，如机械加工、建材生产等工厂。这些企业主要集中在镇工业园区以及部分交通便利的区域，靠近主要交通干线，便于原材料的运输和产品的输出。机械加工企业在生产过程中会产生含重金属（如铜、锌、铅等）的废水，这些重金属废水若未经有效处理直接排放，可能会通过地表径流或渗漏的方式进入地表水和地下水系统，导致水源区水中重金属含量超标，对人体健康和生态环境造成危害。建材生产企业排放的废水中可能含有大量的悬浮物、碱性物质以及氟化物等，会改变水体的酸碱度和化学组成，影响水生生物的生存环境，还可能与地下水中的化学成分发生反应，进一步影响水文地球化学特征。生活污染源主要包括居民生活污水和垃圾填埋场。随着天元镇人口的增长，居民生活污水排放量不断增加。生活污水主要通过污水管网收集，但部分老旧小区和农村地区的污水管网不完善，存在污水直接排放的情况。这些生活污水中含有大量的有机物、氮、磷等污染物，如蛋白质、碳水化合物、尿素等，会导致水体富营养化，引发藻类等浮游生物大量繁殖，消耗水中的溶解氧，使水体发黑发臭，影响水源区的水质和景观。垃圾填埋场位于镇边缘的低洼地带，填埋的生活垃圾在微生物的分解作用下会产生大量的渗滤液。渗滤液中含有高浓度的氨氮、重金属、有机物等污染物，若垃圾填埋场的防渗措施不到位，渗滤液可能会渗漏到地下水中，对地下水水质造成严重污染。在垃圾填埋场周边的地下水监测点，已检测出氨氮和重金属等污染物的含量明显高于其他区域，表明垃圾填埋场对周边地下水水质产生了较大影响。农业面源污染在水源区也较为普遍。天元镇是农业生产区，大量农田分布在水源区范围内。农业生产中广泛使用的化肥和农药是主要的污染源之一。化肥的主要成分包括氮、磷、钾等，过量使用化肥会导致土壤中养分过剩，通过地表径流和淋溶作用，这些养分进入地表水和地下水，增加水中的氮、磷含量，引发水体富营养化。农药的使用虽然有助于防治病虫害，但也会残留在土壤和水体中，对水生生物和人体健康产生潜在危害。在一些农田附近的沟渠和河流中，已检测出农药残留，如有机磷农药、拟除虫菊酯类农药等，说明农业面源污染对水源区水体的影响不容忽视。畜禽养殖也是农业面源污染的重要来源。养殖场产生的畜禽粪便和污水中含有大量的有机物、氨氮、病原体等污染物。部分养殖场的环保设施不完善，畜禽粪便和污水未经处理直接排放到周边环境中，会对地表水和地下水造成污染。在一些畜禽养殖场附近的水体中，氨氮含量严重超标，且存在大量的大肠杆菌等病原体，对水源区的水质安全构成了威胁。五、氨氮污染现状分析5.1氨氮污染监测数据收集与整理为全面了解旌阳区天元镇水源区氨氮污染现状，在水源区内合理设置了多个监测点。考虑到不同区域的水文地质条件、污染源分布以及用水功能等因素，共设置了[X]个地下水监测点和[X]个地表水监测点。地下水监测点涵盖了第四系松散岩类孔隙水和基岩裂隙水分布区域，包括在河流附近、工业集中区周边、农业灌溉区以及居民生活区内等不同位置，以确保能够监测到不同来源和不同类型的地下水氨氮污染情况。地表水监测点主要分布在绵远河、石亭江等主要河流，以及区内的灌溉渠和水库等地表水体中，在河流的上中下游、支流与干流交汇处、灌溉渠取水口和出水口等关键位置设置监测点，以全面掌握地表水中氨氮的分布和变化情况。监测频率为每月一次，在每月固定的时间进行采样监测，以保证数据的可比性和连续性。在特殊时期，如暴雨后、农业灌溉高峰期、工业生产旺季等，会增加监测次数，以便及时了解氨氮污染在特殊情况下的变化趋势。在一次暴雨过后，立即对各监测点进行了加密监测，发现部分地表水监测点的氨氮浓度出现了明显的上升，这表明暴雨可能会将地表的污染物冲刷进入水体，导致氨氮污染加剧。对收集到的氨氮监测数据进行系统整理。首先，对数据进行初步审核，检查数据的完整性和准确性，剔除异常值和错误数据。在审核过程中，发现个别监测点的氨氮浓度数据明显偏离其他监测点同期数据，经核实是由于采样过程中的操作失误导致数据错误，因此将该数据剔除。然后，将整理后的数据按照监测点、监测时间等信息进行分类存储，建立数据库，以便后续分析使用。运用专业的数据处理软件，对数据进行统计分析，计算氨氮浓度的平均值、中位数、最大值、最小值、标准差等统计参数，绘制氨氮浓度随时间和空间变化的图表，直观展示氨氮污染的变化趋势和分布特征。通过数据整理和分析，发现水源区部分区域的氨氮浓度呈现出逐渐上升的趋势，尤其是在工业集中区和农业灌溉区附近，氨氮浓度相对较高，这为进一步深入研究氨氮污染的原因和制定防治措施提供了重要依据。5.2氨氮污染现状评估通过对监测数据的深入分析，绘制氨氮浓度空间分布图（如图5-1所示），清晰展示氨氮浓度在水源区内的空间分布特征。从图中可以看出，氨氮浓度在空间上呈现出明显的不均匀分布。在工业集中区和农业灌溉区附近，氨氮浓度相对较高。在某工业集中区周边的地下水监测点，氨氮浓度高达[X]mg/L，远高于其他区域。这是由于工业生产过程中排放的含氨氮废水，以及农业生产中过量使用的氮肥，通过地表径流和渗漏等方式进入地下水，导致这些区域氨氮污染较为严重。而在远离污染源的山区和河流上游部分区域，氨氮浓度相对较低，一般在[X]mg/L以下，这些区域受人类活动影响较小，生态环境相对较好，氨氮污染程度较轻。[此处插入氨氮浓度空间分布图5-1，图中应标注清楚地理方位、监测点位置以及氨氮浓度等值线，等值线间隔合理，能够清晰展示氨氮浓度的空间变化趋势]依据国家相关水质标准，对氨氮浓度的超标情况进行统计分析。在地下水监测点中，有[X]%的监测点氨氮浓度超过了《地下水质量标准》（GB/T14848-2017）中III类水标准（氨氮浓度限值为0.5mg/L）。在某村庄附近的地下水监测点，氨氮浓度达到0.8mg/L，超过标准限值，这可能是由于村庄生活污水排放未经有效处理，渗入地下导致地下水氨氮超标。在地表水监测点中，也有部分监测点氨氮浓度超过了《地表水环境质量标准》（GB3838-2002）中相应功能区的标准限值。在一些灌溉渠的出水口，氨氮浓度较高，这是因为灌溉过程中，农田中的含氮污染物随灌溉水流入渠道，使得渠道水中氨氮含量升高。氨氮超标会对水体生态系统和人类健康产生严重危害，过高的氨氮含量会消耗水中的溶解氧，导致水体富营养化，引发藻类等浮游生物大量繁殖，破坏水生生态平衡；人类长期饮用氨氮超标的水，可能会引发一系列健康问题。对氨氮浓度的时间变化趋势进行分析，绘制氨氮浓度随时间变化的折线图（如图5-2所示）。从图中可以看出，氨氮浓度在不同时间段呈现出一定的波动变化。在农业灌溉高峰期（一般为春季和夏季），由于大量使用氮肥，地表径流增加，氨氮浓度会出现明显的上升。在2023年夏季，某农业灌溉区附近的地表水监测点氨氮浓度从0.3mg/L上升至0.6mg/L。而在冬季，农业活动减少，氨氮浓度相对较低。工业生产的季节性变化也会对氨氮浓度产生影响，一些工业企业在生产旺季排放的废水较多，导致周边水体氨氮浓度升高。从长期趋势来看，部分区域的氨氮浓度呈现出逐渐上升的趋势，这表明氨氮污染问题有逐渐加重的态势，需要引起高度重视。[此处插入氨氮浓度随时间变化的折线图5-2，图中横坐标为时间，纵坐标为氨氮浓度，应清晰标注不同监测点或区域的折线，以及时间刻度和氨氮浓度刻度]5.3氨氮污染来源解析运用多元统计分析、同位素分析以及污染源调查等多种方法，对氨氮污染来源进行深入