兴平市地下水动态特征解析与数值模拟探究

兴平市地下水动态特征解析与数值模拟探究一、绪论1.1研究背景与意义1.1.1研究背景兴平市地处关中平原中部，是连接西安与宝鸡的重要节点城市，地理位置优越，交通便利，经济发展迅速。然而，兴平市面临着严峻的水资源短缺问题，多年平均降水量为555.5mm，降水年内分配不均，年际变化较大，时空分布不均，且境内虽有渭河过境，但受地形和水质污染限制，渭河地表水难以利用。除农业灌溉部分引用宝鸡峡渠系水外，各部门用水主要依赖地下水。兴平市水资源总量为1.68亿m³，人均地下水资源占有量仅280m³，约为全国平均水平的13％，远低于国际公认的500m³的警戒线，水资源供需矛盾突出。随着兴平市经济社会的快速发展，尤其是四大工业园区建设的强力推进，大量规模工业企业落户，工业用水量增长迅速，加之人口不断增长，城市框架持续拉大，生活用水量也日益增加，对地下水的开采量不断攀升。长期以来，不合理的地下水开发利用导致了一系列问题，其中最突出的是“兴化漏斗区”的形成。由于兴平化肥厂等企业长期大量抽取地下水，使得该区域地下水位持续下降，形成了以采水点为中心，四周向中心呈梯度下降的漏斗形区域。据相关资料显示，兴平化肥厂水源地的水位下降幅度尤为明显，严重影响了当地的生态环境和水资源可持续利用。过度开采地下水还引发了地面沉降、建筑物开裂等地质灾害，对城市的基础设施和居民的生命财产安全构成了威胁。1.1.2研究意义本研究对兴平市地下水动态特征进行深入分析，并建立数值模拟模型，具有重要的现实意义。通过研究可以准确掌握兴平市地下水的动态变化规律，包括水位、水流流场和化学特征等，为合理开发利用地下水资源提供科学依据，实现水资源的可持续利用。对“兴化漏斗区”等问题进行研究，有助于制定针对性的防治措施，减缓地下水位下降速度，修复受损的生态环境，保护地质环境安全。水资源是经济社会发展的重要支撑，合理的地下水开发利用和保护能够保障兴平市工业、农业和生活用水需求，促进经济社会的可持续发展，为兴平市打造产业高地、推动经济腾飞提供有力的水资源保障。1.2研究内容与技术路线1.2.1研究内容本研究将对兴平市地下水动态特征进行全面分析。收集兴平市多年的降雨数据，分析降雨的年际变化、年内分配以及空间分布特征；研究过境河流渭河的径流特征，包括径流量的年际变化、年内分配等。通过对兴平市地下水监测井网的观测资料进行整理和分析，研究地下水位的年际变化和年内变化规律，探讨不同地貌单元地下水位的变化差异。利用数值模拟软件，绘制兴平市不同时期的地下水水流流场图，分析水流的运动方向、流速分布等特征，研究水流流场的变化规律及其影响因素。采集兴平市地下水水样，分析地下水的化学组成，包括主要离子成分、酸碱度、硬度等，研究地下水化学特征的空间分布规律及其与地质条件、人类活动的关系。本研究还将建立兴平市地下水水流数值模拟模型。根据兴平市的地形地貌、地质构造、水文地质条件等资料，建立地下水系统概念模型，明确模拟区域的范围、边界条件、含水层结构等。选用合适的数学模型来描述兴平市地下水的运动规律，确定模型的基本方程和定解条件。对兴平市地下水系统进行概化，包括含水层的结构特征概化、地下水流动特征概化以及边界条件的概化等，将复杂的实际问题简化为数学模型能够处理的形式。详细分析兴平市地下水的补给和排泄条件，确定补给量和排泄量的计算方法，为模型的建立提供准确的参数。利用数值模拟软件FEFLOW，根据前面建立的概念模型、数学模型和概化结果，建立兴平市地下水水流数值模拟模型，并编辑模型的各项参数，如渗透系数、给水度、储水系数等。通过将模型计算结果与实际观测数据进行对比，对模型进行识别和验证，调整模型参数，使模型能够准确地反映兴平市地下水的实际运动情况。运用建立好的数值模拟模型，对兴平市地下水的未来变化趋势进行预测，为水资源管理和保护提供决策依据。1.2.2技术路线本研究的技术路线如图1-1所示。首先，广泛收集兴平市的气象、水文、地质、地下水监测等相关资料，对资料进行整理和分析，了解研究区的自然环境和地下水开发利用现状。基于收集到的资料，对兴平市地下水的降雨径流特征、水位特征、水流流场特征和化学特征进行详细分析，为建立数值模拟模型提供基础数据和理论支持。根据地下水系统的特点和研究目的，建立地下水系统概念模型和数学模型，对研究区进行系统概化，确定模型的各项参数和边界条件。利用数值模拟软件FEFLOW建立兴平市地下水水流数值模拟模型，对模型进行识别和验证，确保模型的准确性和可靠性。最后，运用验证后的模型对兴平市地下水的动态变化进行模拟和预测，并根据模拟结果提出合理的水资源管理和保护建议。\begin{matrix}&\text{资料收集}&&\text{气象、水文、地质、地下水监测资料}&\\&\downarrow&&\downarrow&\\&\text{特征分析}&&\text{降雨径流、水位、水流流场、化学特征分析}&\\&\downarrow&&\downarrow&\\&\text{模型建立}&&\text{概念模型、数学模型、系统概化、补给排泄条件分析}&\\&\downarrow&&\downarrow&\\&\text{模型识别与验证}&&\text{对比模拟结果与实际观测数据，调整模型参数}&\\&\downarrow&&\downarrow&\\&\text{模型应用}&&\text{模拟和预测地下水动态变化，提出管理和保护建议}&\end{matrix}图1-1技术路线图二、地下水数值模拟研究现状及发展趋势2.1国外地下水数值模拟研究现状2.1.1地下水数值模拟方法国外在地下水数值模拟方法研究方面处于前沿地位，发展出了多种成熟且高效的模拟方法，以满足不同场景和研究目的的需求。水流模拟是地下水数值模拟的基础，主要用于模拟地下水的流向以及地下水水头与时间的关系。它基于达西定律，通过建立数学模型来描述地下水在多孔介质中的流动过程。例如，有限差分法将模拟区域离散为网格，把连续的偏微分方程转化为差分方程进行求解，能够较为准确地计算地下水在各网格节点的水头值和流量，从而分析水流的运动趋势，在平原地区含水层的水流模拟中应用广泛。有限元法则是将模拟区域划分为有限个单元，通过对单元进行插值和加权余量计算，得到近似解，其对复杂边界条件和不规则区域的适应性更强，在山区等地形复杂的地下水水流模拟中优势明显。地下水运移模拟主要关注地下水、热和溶质组分的运移速率。在实际应用中，要特别考虑“优先流”现象，即局部具有高连通性的渗透性区域，使得水、热、溶质组分在该处的运移速率快于周围地区。为了准确模拟这一现象，研究人员开发了多种模型，如基于随机理论的随机模型，考虑了介质的非均质性和不确定性对运移的影响，能够更真实地反映实际情况；多尺度模型则通过对不同尺度的介质特性进行分析和模拟，有效处理了不同尺度下的运移问题，在研究大区域地下水溶质运移时发挥了重要作用。反应模拟聚焦于模拟水中、气-水界面、水-岩界面所发生的物理、化学、生物反应。在研究地下水与土壤中重金属离子的相互作用时，通过建立化学平衡模型，结合离子交换、溶解沉淀等反应方程，能够预测地下水中重金属离子的浓度变化和迁移路径。生物反应模拟则考虑微生物的代谢活动对地下水化学组成的影响，在研究地下水生物修复过程中具有重要意义。反应运移模拟将地球化学模拟（包括动力学模拟）和溶质运移模拟（包括非饱和介质二维、三维流）有机结合，是当前地下水模拟的重要发展趋势。它能够模拟地下水运移过程中发生的溶解与沉淀、吸附与解吸、氧化与还原、配合、中和、生物降解等各种复杂反应。在模拟垃圾填埋场渗滤液对地下水的污染时，综合考虑渗滤液中各种污染物的运移以及与地下水和周围介质发生的化学反应，能够更准确地评估污染的范围和程度，为制定合理的防治措施提供依据。2.1.2地下水数值模拟软件国外开发了众多功能强大、应用广泛的地下水数值模拟软件，这些软件各具特色，适用于不同的研究和工程需求。FEFLOW（FiniteElementsubsurfaceFLOWsystem）是一款由德国WASY公司开发的地下水流与溶质运移模拟软件，在国际上享有很高的声誉。它采用有限元数值模拟方法，具备模拟多维的、非稳定流动和传输过程的能力，支持多种边界条件和复杂的地质模型。FEFLOW可以模拟饱和与非饱和流、稳定与非稳定流、溶质运移、热传递、密度驱动流等多种情况。在海水入侵模拟中，FEFLOW能够准确模拟海水与淡水的界面变化以及盐分在地下水中的扩散过程，为沿海地区水资源管理提供科学依据。其界面友好，操作相对简便，拥有丰富的前后处理功能，便于用户进行模型构建、参数设置和结果分析。Modflow（Modularthree-dimensionalfinite-differenceground-waterflowmodel）是美国地质调查局开发的基于有限差分法的孔隙介质中三维地下水流模拟程序，具有程序结构模块化、离散方法简单化和求解方法多样化等优点，是目前应用最为广泛的地下水流模拟程序之一。它能够模拟稳定和非稳定地下水流，适用于非均质含水层、简单的各项异性含水层，以及承压、潜水、承压-无压含水层的地下水流模拟。在研究大型平原地区的地下水资源时，Modflow通过对含水层结构、边界条件和源汇项的合理设置，能够准确预测地下水位的变化和水资源的分布情况。Modflow还不断进行更新和扩展，推出了多个版本，如MODFLOW88、92、6、2000和2005等，功能日益完善。MT3DMS是一款专门用于模拟地下水系统中对流、弥散和化学反应的三维溶质运移模型，常与Modflow一起使用。它能够考虑多种溶质的运移过程，以及溶质之间的化学反应，在研究地下水污染扩散方面应用广泛。在模拟化工园区地下水污染时，MT3DMS可以准确模拟污染物在地下水中的迁移路径和浓度变化，为污染治理提供有力支持。此外，还有一些其他的地下水数值模拟软件，如GMS（GroundwaterModelingSystem），它是一个综合性的地下水模拟软件平台，集成了多种地下水模型，如MODFLOW、FEMWATER、MT3DMS等，具有强大的前处理和后处理功能，可实现模型的自动化运行和可视化展示；SEAWAT是一款用于模拟变密度地下水流动和溶质运移的软件，特别适用于研究海水入侵、盐水-淡水界面移动等问题。这些软件在不同的领域和场景中发挥着重要作用，推动了地下水数值模拟技术的发展和应用。2.2FEFLOW软件2.2.1简介FEFLOW（FiniteElementsubsurfaceFLOWsystem）作为一款先进的地下水流与溶质运移模拟软件，具有卓越的功能和优势。它能够实现对地下水系统的全面模拟，涵盖饱和与非饱和流、稳定与非稳定流、溶质运移、热传递以及密度驱动流等多种复杂过程。FEFLOW采用有限元数值模拟方法，该方法在处理复杂地质条件和边界条件时具有独特的优势，能够将模拟区域划分为各种形状的有限单元，从而更精确地拟合实际的地质结构和水流状况。在功能方面，FEFLOW拥有丰富的模型库和参数设置选项，用户可以根据研究区域的具体情况灵活选择和调整模型参数，以确保模拟结果的准确性。软件内建了直观友好的图形用户界面，大大降低了用户的操作难度。在该界面下，用户能够方便地进行模型的构建，包括定义模拟区域的范围、边界条件和初始条件等；还可以轻松设置各种水文地质参数，如渗透系数、储水系数、给水度等；并且能够实时控制模拟的运行过程，随时查看模拟进度和中间结果。在模拟完成后，FEFLOW提供了强大的结果分析功能，帮助用户深入理解地下水系统的动态特性。FEFLOW在应用领域上十分广泛，涵盖了水文地质学、环境工程、矿业工程、地下水科学等多个学科领域。在水文地质学研究中，它可用于分析地下水的补给与排泄规律，研究地下水与地表水之间的相互作用关系；在环境工程领域，能够模拟污染物在地下水中的迁移转化过程，为地下水污染防治提供科学依据；在矿业工程中，可用于评估矿山开采对地下水水位和水质的影响，制定合理的水资源保护和利用方案；在地下水科学研究中，FEFLOW为深入探究地下水系统的演化机制和规律提供了有力工具。2.2.2应用领域FEFLOW在众多领域都有成功的应用案例，展现出其强大的模拟能力和实际应用价值。在水文地质领域，FEFLOW被广泛应用于研究地下水的补给与排泄过程。通过建立精确的水文地质模型，结合实际的气象、水文和地质数据，FEFLOW能够准确模拟不同补给源（如降雨入渗、河流渗漏等）对地下水水位和流量的影响，以及地下水的排泄方式（如蒸发蒸腾、向河流排泄等）。在某山区流域的水文地质研究中，利用FEFLOW模拟发现，该地区地下水的主要补给来源是降雨入渗，且在不同地形和地质条件下，补给量存在显著差异。通过对模拟结果的分析，为该地区水资源的合理开发和管理提供了重要依据。在水资源管理方面，FEFLOW可用于优化地下水的开采方案。通过模拟不同开采强度和开采布局下地下水水位的变化，预测可能出现的水资源问题，如地面沉降、海水入侵等，从而制定出科学合理的地下水开采计划。在沿海地区的某城市，为了避免过度开采地下水导致海水入侵，利用FEFLOW对不同开采方案进行模拟分析，最终确定了最优的开采方案，既满足了城市发展的用水需求，又保护了地下水资源的可持续性。在地下水污染防治领域，FEFLOW能够模拟污染物在地下水中的迁移过程及其时间空间分布模式。在某化工污染场地的治理中，通过FEFLOW模拟污染物的扩散路径和浓度变化，为制定污染修复方案提供了关键信息。根据模拟结果，确定了污染的范围和程度，从而有针对性地采取了地下水抽提、原位修复等措施，有效降低了污染物对地下水和周边环境的影响。在地热资源开发利用中，FEFLOW可用于模拟地下热水的流动和热量传递过程。通过建立热-流耦合模型，分析地热田的温度分布和热储层的特性，为地热资源的合理开发和高效利用提供技术支持。在某地热田的开发项目中，利用FEFLOW模拟不同开采方式下热储层的温度变化和压力分布，优化了地热井的布局和开采参数，提高了地热资源的开发效率。2.2.3系统输入特点FEFLOW对系统输入具有明确的要求和灵活的处理方式，以确保能够准确地反映实际的地下水系统。在数据格式方面，FEFLOW具有良好的兼容性，输入数据格式既可以是常见的ASCII码文件，方便用户进行数据的编辑和处理；也可以直接读取GIS地理信息系统文件，如ArcView的shp格式，充分利用GIS强大的空间分析和数据管理功能，实现数据与空间的统一管理，提高数据利用效率和模拟精度。这使得用户能够方便地整合各种来源的数据，包括地形数据、地质数据、水文数据等，为模型的建立提供全面准确的信息。对于数据精度，FEFLOW要求用户在输入数据时尽可能保证准确性和可靠性。在获取渗透系数等关键水文地质参数时，需要通过现场试验、实验室测试等多种方法进行精确测定。对于离散的空间抽样数据，FEFLOW提供了克里格（Kriging）、阿基玛（Akima）和距离反比加权（IDW）三种方法进行内插或外推，以提高数据的空间连续性和精度，从而更准确地反映地下介质的特性和水流状态。在边界条件和初始条件的定义上，FEFLOW给予用户充分的灵活性。用户可以视具体情况定义第一、第二和第三类边界，并且还可以对边界条件增加特定的限制条件，以满足复杂的实际情况。在模拟河流与地下水的相互作用时，可以根据河流的水位变化、流量大小以及河床的渗透特性等因素，准确设定河流边界条件及其限制条件，使得模拟结果更符合实际情况。FEFLOW还能方便地定义复合含水层中的分层开采井和混合开采井以及注水井，所有边界条件及其限制条件、渗透系数、补排量等参数既可设置为常数，也能定义为随时间变化的函数，以适应不同的水文地质条件和动态变化的地下水系统。2.2.4数值求解方法FEFLOW采用了先进的数值求解方法，以确保模拟结果的准确性和计算效率。其中，有限单元法是FEFLOW的核心求解方法之一。该方法的基本原理是将连续的求解区域离散为有限个单元，这些单元通过节点相互连接。对于每个单元，基于变分原理或加权余量法，将描述地下水流动或溶质运移的偏微分方程转化为一组代数方程。在处理二维地下水流动问题时，将模拟区域划分为三角形或四边形单元，通过对每个单元内的水头或溶质浓度进行插值近似，建立起单元的离散方程。然后，将所有单元的方程组装成一个整体的方程组，通过求解这个方程组得到整个模拟区域的解。有限单元法的优势在于能够灵活地处理复杂的几何形状和边界条件，对于不规则的模拟区域和具有复杂地质结构的情况具有很好的适应性。除了有限单元法，FEFLOW还采用了多种快速直接求解法，如PCG（预调共轭梯度法）、BICGSTAB（双共轭梯度稳定法）、CGS（共轭梯度平方算法）、GMRES（广义最小残差法）以及带预处理的再启动ORTHOMIN法等。这些求解方法针对不同类型的线性方程组具有各自的优势，能够根据具体问题选择最合适的求解算法，提高计算效率和收敛速度。在处理大规模的地下水模拟问题时，PCG方法能够有效地减少计算量和内存需求，快速收敛到准确的解；而GMRES方法在处理非对称矩阵时表现出色，能够稳定地求解复杂的地下水流动方程。FEFLOW还运用了灵活多变的up-wind技术，如流线up-wind和奇值捕捉法（Shockcapturing）。流线up-wind技术通过考虑水流的方向，对对流项进行更准确的离散处理，减少了数值弥散现象，提高了对溶质运移模拟的精度；奇值捕捉法则能够有效地捕捉到水流或溶质运移过程中的突变现象，如在模拟地下水与海水的界面移动时，能够准确地描述界面的变化情况，使得模拟结果更符合实际物理过程。2.2.5系统结果输出及显示FEFLOW提供了丰富多样的结果输出形式和直观的可视化方式，方便用户对模拟结果进行分析和理解。在结果输出方面，FEFLOW可以输出多种数据类型，包括地下水水位、流速、流量、溶质浓度、温度等随时间和空间的变化数据。这些数据以文本文件或二进制文件的形式保存，用户可以根据需要进行进一步的处理和分析。FEFLOW还支持输出模拟过程中的中间结果，如迭代过程中的残差变化、各时间步的计算结果等，为用户深入了解模拟过程和检查计算的稳定性提供了便利。在可视化方面，FEFLOW内建了强大的图形显示功能。用户可以通过图形界面直观地查看模拟结果，包括绘制二维和三维的等值线图、矢量图、剖面图等。通过等值线图，可以清晰地展示地下水水位、溶质浓度等变量在平面上的分布情况，直观地反映出高值区和低值区的位置和范围；矢量图则能够形象地表示地下水的流动方向和流速大小，帮助用户分析水流的运动趋势；剖面图可以展示不同深度处地下水的特性变化，对于研究含水层的结构和水流特征具有重要意义。FEFLOW还支持动画制作功能，用户可以将不同时间步的模拟结果制作成动画，动态地展示地下水系统的变化过程，更直观地观察地下水水位、水流流场和溶质运移等的动态变化。FEFLOW还允许用户将模拟结果导出到其他专业的绘图软件或数据分析软件中进行进一步的处理和展示，以满足不同用户的需求和个性化的分析要求。将模拟结果导出到ArcGIS中，利用其强大的地理信息分析和制图功能，结合地形、土地利用等其他地理信息数据，进行更深入的空间分析和可视化表达。2.3数值模拟在我国的发展与应用我国的地下水数值模拟研究起步相对较晚，但在过去几十年中取得了显著的进展。自20世纪70年代引入数值模拟技术以来，我国科研人员和工程技术人员积极开展相关研究和应用工作，不断推动数值模拟技术在我国的发展。早期，主要是对国外先进的数值模拟方法和软件进行学习和引进，结合国内的实际情况进行应用探索。随着我国科研实力的不断提升和计算机技术的飞速发展，国内在地下水数值模拟领域的研究逐渐深入，自主研发能力也不断增强。在不同地区，数值模拟技术得到了广泛的应用。在华北平原，由于长期面临着严重的水资源短缺和地下水超采问题，数值模拟技术被大量应用于地下水资源评价和管理。通过建立高精度的地下水数值模型，模拟不同开采方案下地下水位的变化趋势，评估水资源的可持续性，为制定合理的水资源管理政策提供了科学依据。在上海等沿海城市，为了应对地面沉降和海水入侵等问题，利用数值模拟技术研究地下水开采与地面沉降之间的关系，以及海水入侵的规律和影响因素，从而采取有效的防控措施，保护城市的地质环境和水资源安全。在不同领域，数值模拟也发挥着重要作用。在水利工程领域，数值模拟可用于水库、大坝等水利设施对地下水环境的影响评估。通过模拟水库蓄水后周边地下水位的变化，预测可能出现的浸没、渗漏等问题，为工程设计和运行管理提供参考。在矿山开采领域，数值模拟能够评估采矿活动对地下水的影响，包括矿井涌水、地下水污染等，为矿山的安全生产和环境保护提供技术支持。在农业灌溉领域，数值模拟可用于优化灌溉方案，通过模拟不同灌溉方式和灌溉量下土壤水分和地下水的动态变化，提高水资源利用效率，实现农业的可持续发展。近年来，我国在地下水数值模拟软件的研发方面也取得了一定的成果。一些科研机构和高校开发了具有自主知识产权的地下水数值模拟软件，这些软件在功能和性能上不断提升，逐渐在国内市场占据一席之地。同时，国内对数值模拟技术的研究也更加注重多学科交叉融合，将地下水数值模拟与地质、水文、环境等学科相结合，开展综合性的研究，为解决复杂的实际问题提供了更全面的方法和思路。2.4数值模拟软件的发展趋势及前景展望2.4.1数值模拟软件的发展趋势随着计算机技术、信息技术和人工智能技术的不断发展，数值模拟软件呈现出智能化、集成化和可视化的发展趋势。智能化是当前数值模拟软件发展的重要方向之一。未来的数值模拟软件将引入更多的人工智能元素，实现自动网格划分、智能参数优化和人工智能辅助决策等功能。自动网格划分功能能够根据模拟区域的几何形状和物理特性，自动生成高质量的网格，减少人工干预，提高建模效率和准确性。智能参数优化则利用机器学习算法，根据实际观测数据自动调整模型参数，使模拟结果更符合实际情况。人工智能辅助决策功能可以根据模拟结果和预设的目标函数，为用户提供决策建议，帮助用户制定最优三、研究区概况及地下水开发利用现状3.1自然环境3.1.1兴平位置兴平市地处陕西省关中平原中部，地理位置优越，是连接关中地区重要城市的关键节点。其地理坐标为东经108°17′49″至108°37′7″，北纬34°12′50″至34°26′53″。兴平市东连咸阳市秦都区，为其与咸阳的经济、文化交流提供了便利条件，促进了区域一体化发展；南临渭水，与周至县、鄠邑区隔河相望，独特的地理位置使得兴平市在水资源利用和区域协作方面具有重要意义；西邻武功，加强了区域间的产业合作与资源共享；北依莽山与礼泉、乾县接壤，复杂的地形地貌造就了多样的生态环境和丰富的自然资源。兴平市境内交通网络发达，陇海铁路、“312”国道和104省道中线、北线横贯东西，便捷的交通条件使其与周边地区的联系更加紧密，为经济发展和人员往来提供了有力支撑。兴平市还处于“关中—天水”经济区建设腹地以及“西咸阳经济圈”西大门的重要位置，作为一个以大中城市为依托的“卫星城”，兴平市能够充分利用周边大城市的资源和优势，实现自身的快速发展。东距西安航空港咸阳国际机场约20千米，便捷的航空交通为兴平市的对外交流和经济合作提供了广阔的空间，促进了区域间的人流、物流和信息流的快速流动。西经扶风县至宝鸡市162千米，使得兴平市在关中地区的交通网络中扮演着重要的角色，加强了与宝鸡及西部地区的经济联系。兴平市独特的地理位置使其在区域发展中具有重要的战略地位，为其经济、社会和生态的可持续发展奠定了坚实的基础。3.1.2气候兴平市属暖温带半湿润半干旱大陆性季风气候，这种气候类型使得兴平市四季分明，光照充足，雨热同季。春季，蒙古高压逐渐衰退，西风带移动性槽脊活动频繁，热带暖气团逐渐北进，大地回暖迅速，降水随之增多。然而，由于冷空气活动依旧频繁，天气变化无常，易出现寒潮、霜冻、大风等灾害性天气，且多浮尘，春旱现象时有发生。据统计，春季降水量占全年降水量的20%-25%，但由于蒸发量大，春旱对农业生产的影响较为严重，约有30%-40%的年份会出现不同程度的春旱。夏季，兴平市主要受副热带高压影响，气候炎热多雨，降水多以暴雨和雷阵雨形式出现，是一年中降水最多的季节，占全年降水量的40%-50%。降水集中且强度大，容易引发洪涝灾害。夏季也是农作物生长的关键时期，充沛的降水和充足的光照有利于农作物的生长，但高温天气也可能对部分农作物造成热害，影响作物的产量和品质。秋季，蒙古高压逐渐增强，副热带高压开始南撤，北方冷空气南下。受秦岭山系阻挡，锋区南下迟缓，常形成准静止锋，导致初秋出现连阴雨秋汛期，降水占全年降水量的20%-25%。这种天气会影响农作物的收获和晾晒，增加病虫害发生的风险，对农业生产造成一定的损失。若秋季极地高压迅速建立并南下，同时副热带高压仍维持在较高纬度，高低空高压发生重合现象，则会出现秋高气爽的天气，有利于农作物的后期生长和收获。冬季，兴平市主要受蒙古高压影响，气候寒冷干燥，雨雪稀少，常有冬旱发生，降水量仅占全年降水量的5%-10%。寒冷的天气对农作物的越冬和畜牧业的发展带来一定挑战，需要采取相应的防寒保暖措施，以保障农业生产和畜牧业的稳定。兴平市年平均气温11.6℃，极端最高气温达42.2℃（1966年6月15日），极端最低气温可达-18.6℃（1991年12月28日）。年降水量584.7毫米，降水年内分配不均，7-9月降水量最多，占全年降水量的47%以上，1、2及12月份降水量最少，占全年降水量的5%以下。无霜期平均218天左右，全年日照时数2065.2小时。气候条件对兴平市地下水的影响显著，降水是地下水的主要补给来源之一，降水的时空分布直接影响地下水的补给量和水位变化。夏季降水集中，入渗补给量大，地下水位会相应上升；而冬季降水稀少，蒸发量大，地下水位则会有所下降。气温的变化也会影响土壤水分的蒸发和植物的蒸腾作用，进而间接影响地下水的动态变化。3.1.3地形地貌兴平市地势呈现西北高、东南低的态势，地貌主要可划分为南部渭河阶地区和北部黄土台塬区两个类型区。南部渭河阶地区属于渭河平原的一部分，是渭河干流长期切割形成的冲积平原。其北起晁庄、界庄、前北铁、李家坡、板桥、北汤台、杏花村、庄头、马堪村、来祁寨、窦马村以南，西起武功县界，东到秦都区界，南止渭河，南北宽约18.5公里，东西长约32公里，面积达324.47平方公里，占全市总面积的63.94％。该区域发育完整，阶面平坦宽阔，河漫滩与一级、二级、三级阶地特征明显。河漫滩位于本市沿渭西桥、永流、建坊、龙兴、吴耳、郭公、高家、北张、马村、仪空、蔡东、蔡西、千佛、张耳、塔耳、何孔、马坊、侯村、湾里、南田村南部，发育不甚完整，南北宽0.2-2.5公里，东西长32公里，海拔约390米，是本市渭河阶地最低区。在区内北部吴耳村和村以南有宽约1.5公里，长6公里的漫滩阶地，高出漫滩0.5米，常被季节性洪水淹没，组成物质以全新统的中、细沙为主，沙质粘土和沙卵石次之，渗漏严重，肥性差，不耐旱涝，虽有部分耕种，但收成难以保障。一级阶地位于河漫滩地区以北，台面平坦，并与二级阶地以坎相接，高出河漫滩2-3米，台原宽窄不一，距地表2-3米以内均有含水层，上部为潜水，下部为承压水，地质层是第四系的上更新统，上部为沙质粘土，下部为沙及沙卵石层，厚约190米，土壤主要是斑斑黑油土，具有蓄水保肥、抗旱耐涝、耕种期长、适种性广等特性。二级阶地在一级阶地以北，台面平坦，在杏花村以西与黄土台原相接，台面宽9公里，在杏花村以东与三级阶地相连，台面宽8.5公里，高出一级阶地10-15米，地下水埋深5-10米，含水层厚约30米，地质层仍属第四系上更新统，上部为黄土质沙质粘土，下部为沙及沙石，沙质粘土、沙砾石覆盖层厚约156米，土壤主要是黑油土，土层深厚肥沃，适于农耕。三级阶地位于二级阶地以北，黄土台原区以南，在杏花村以东呈楔形分布，以西地貌界线不甚明确，宽约1公里，高出二级阶地约20-50米，地下水埋深10-15米，含水层厚10-15米，地质层是第四系的中更新统，系渭河沉积层，20米以上为黄土类，以下见沙、沙卵石层，土壤主要是红油土，适种小麦。北部黄土台塬区位于县境北部，属于风积平原，奠基于中更新世，形成于上更新世，属于地质第四系上更新统，海拔460-541.8米，平均比降16.63‰，在杏花村以东高出三级阶地20-41米。台原面广阔平坦，南缘坡度较陡，南北向冲沟多，沟壁陡直，有天生桥、再生谷等自然地质现象，断面成V形，大部分沟壑已整修为梯田。台原区面积182.96平方公里，占全市总面积的36.06％，是本市农业区划的粮食主要生产区。其基底构造与阶地区相反，属缓慢上升区，上部覆盖深厚的黄土层，可达数百米，钻孔167米深处见有薄的沙层，土壤主要是红油土，面积385612亩，占全市土地总面积的53.5％，适于小麦及多种作物生长，也是本市农业发展潜力较大地区。兴平市的地形地貌对地下水的储存和径流有着重要影响。渭河阶地区地势平坦，含水层分布广泛且厚度较大，有利于地下水的储存和侧向径流，地下水的补给主要来自降水入渗和渭河的渗漏补给，排泄方式主要为人工开采和向渭河排泄。黄土台塬区由于地势较高，地下水埋藏较深，含水层厚度相对较薄，降水入渗是其主要的补给方式，地下水径流缓慢，主要以蒸发和人工开采为排泄方式。不同地貌单元的地形地貌特征导致了地下水的赋存和运移条件的差异，进而影响了地下水的动态变化。3.1.4土壤植被兴平市土壤类型丰富多样，主要包括红油土、黑油土、斑斑黑油土等。红油土主要分布在北部黄土台塬区，其土体深厚，质地适中，通气透水性良好，保水保肥能力较强，富含钾、磷等矿物质养分，肥力较高，适合小麦、玉米、豆类等多种农作物生长。黑油土和斑斑黑油土主要分布在南部渭河阶地区，黑油土土层深厚肥沃，具有良好的团粒结构，保水保肥性能优越，是优质的农业土壤，适合种植蔬菜、水果等经济作物；斑斑黑油土则具有蓄水保肥、抗旱耐涝的特点，耕种期长，适种性广，能满足多种农作物的生长需求。兴平市植被覆盖情况良好，森林植被属暖温带落叶阔叶林带。在北部黄土台塬区，主要植被有刺槐、杨树、柏树等乔木，以及酸枣、荆条等灌木，这些植被在保持水土、涵养水源方面发挥着重要作用。南部渭河阶地区，除了农田植被外，还有一些人工栽培的果树，如苹果树、梨树、桃树等，以及河漫滩和湿地周边的芦苇、菖蒲等水生植物。农田植被主要以小麦、玉米、油菜等农作物为主，随着农业产业结构的调整，经济作物的种植面积逐渐扩大，如蔬菜、瓜果等。土壤和植被对地下水的补给和蒸发有着重要影响。土壤的质地和结构决定了降水的入渗能力，疏松多孔的土壤有利于降水快速入渗，增加地下水的补给量；而质地黏重的土壤则会阻碍降水入渗，减少地下水的补给。植被通过蒸腾作用影响地下水的蒸发，植被覆盖率高的地区，蒸腾作用强烈，会消耗大量的土壤水分，减少地下水的蒸发；同时，植被还能截留降水，减缓地表径流，增加降水入渗的时间和机会，从而促进地下水的补给。在北部黄土台塬区，植被的根系能够深入土壤，改善土壤结构，增加土壤孔隙度，有利于降水入渗，提高地下水的补给量。在南部渭河阶地区，农田植被的合理种植和管理能够减少土壤水分的蒸发，保持土壤湿度，间接影响地下水的动态变化。3.1.5河流水系兴平市境内河流主要为渭河，渭河沿本市境内南缘自西向东流，是兴平市最重要的过境河流。渭河发源于甘肃省渭源县鸟鼠山，流经陇西、漳县、武山、甘谷、天水、宝鸡、眉县、兴平、咸阳、西安、渭南、华县、华阴，至潼关港口注入黄河，全长818公里，流域面积134.934平方公里。在兴平市境内，渭河从桑镇的西桥村入境，流经永流、马家、建坊、吴耳、郭公、高家、北张、东马村南，从仪空、蔡西村开始向东北经蔡东、千佛、张耳、高王、塔耳、何孔、马坊、田南村南至段家村东出境，进入秦都区，市境内长30.5公里，流域面积473.20平方公里。渭河属于季节性多泥沙河流，年平均径流量55.67亿立方米，7-9月份径流量占全年的46%。根据咸阳水文站资料，多年平均流量176.0立方米/秒，最大流量7220立方米/秒（1954年8月18日），最小流量40立方米/秒（1941年7月5日）。年平均水位高程383.29米，最大高程385.79米（1954年8月18日），最小高程382.29米（1938年8月11日）。历史上大洪水年份有1898年的11500立方米/秒，1933年的5800立方米/秒，1954年的7220立方米/秒，基本上41年一遇。含沙量据咸阳水文站资料，年平均值26.8公斤/立方米，最大值596公斤/立方米（1960年8月7日），最小值0.045公斤/立方米（1960年元月13日）。年平均输沙量1.49亿吨（咸阳站），最大输沙量3.31亿吨（1937年），最小输沙量0.52亿吨（1942年）。渭河的河床属于淤积性沙质，河道纵比降1/550-1/2000，河槽宽500-2000米，河岸易被冲塌，河道变化大，30年南北摆动变幅最大4公里，以北移为主，有些河道已靠近一、二级阶地处。渭河冬季有结冰和流凌现象，冰层厚几厘米到十多厘米，冰凌块很少超过5平方米大，凌速为1.3米/秒。渭河与兴平市地下水存在密切的水力联系。在丰水期，渭河水位高于地下水位，渭河河水通过河床渗漏补给地下水，增加地下水的储量，抬高地下水位。在枯水期，地下水位高于渭河水位，地下水则向渭河排泄，维持渭河的基流。这种水力联系使得渭河成为兴平市地下水的重要补给源和排泄通道，对地下水的动态变化产生重要影响。渭河的水质状况也会影响地下水的质量，若渭河受到污染，其渗漏补给可能会导致地下水污染，威胁当地的水资源安全。3.1.6水文地质条件兴平市含水层结构较为复杂，主要包括第四系松散岩类孔隙含水层和新近系碎屑岩类裂隙孔隙含水层。第四系松散岩类孔隙含水层广泛分布于渭河阶地区和黄土台塬区，根据含水层的岩性、结构和水力特征，可进一步划分为潜水含水层和承压含水层。潜水含水层主要分布在渭河阶地区的河漫滩和一级阶地，岩性以中、细沙和沙质粘土为主，厚度一般为10-20米，富水性较好，单井出水量一般在1000-3000立方米/日。承压含水层主要分布在渭河阶地区的二级、三级阶地和黄土台塬区，岩性以沙及沙卵石层为主，厚度较大，一般为30-80米，富水性中等，单井出水量一般在500-1000立方米/日。新近系碎屑岩类裂隙孔隙含水层主要分布在北部黄土台塬区的深部，岩性为砂岩、泥岩互层，裂隙发育程度较低，富水性较差，单井出水量一般小于500立方米/日。兴平市地下水类型主要为孔隙水和裂隙孔隙水。孔隙水主要赋存于第四系松散岩类孔隙中，是兴平市地下水的主要类型，其分布广泛，水量丰富，是工农业生产和生活用水的主要水源。裂隙孔隙水主要赋存于新近系碎屑岩类裂隙孔隙中，由于其富水性较差，开采利用难度较大，目前开发利用程度较低。兴平市不同区域的富水性存在明显差异。渭河阶地区由于含水层厚度大、岩性颗粒粗，富水性较好，是兴平市地下水的主要开采区。黄土台塬区含水层厚度相对较薄，岩性颗粒较细，富水性中等，地下水开采主要集中在靠近渭河阶地的区域。北部山区由于含水层不发育，富水性差，地下水开采量较小。水文地质条件对兴平市地下水的开发利用具有重要影响，了解水文地质条件有助于合理规划地下水开采布局，优化水资源配置，实现地下水资源的可持续利用。3.2社会经济状况3.2.1区域及人口情况兴平市下辖5个街道、8个镇，市政府驻东城街道县门街东路27号。截至2022年，县境内户籍人口597048人，人口分布呈现出一定的特点。城区人口相对密集，东城街道和西城街道作为兴平市的政治、经济和文化中心，集中了大量的人口，人口密度较大，基础设施较为完善，商业、教育、医疗等资源丰富，吸引了众多居民居住和就业。而一些偏远的镇，如桑镇、汤坊镇等，人口相对较少，人口密度较低，这些地区以农业生产为主，经济发展相对滞后，就业机会有限，导致人口外流现象较为明显。兴平市人口的年龄结构也在发生变化，随着社会经济的发展和医疗卫生条件的改善，人口老龄化趋势四、地下水动态特征分析4.1降雨径流特征分析4.1.1降雨特征兴平市多年降雨数据的统计与分析对于了解其水资源状况和地下水补给机制具有重要意义。通过收集兴平市气象站多年来的降雨数据，运用统计分析方法，对降雨量、降雨频率和季节分布特征进行深入研究。兴平市多年平均降雨量为555.5mm，年际变化较大。在过去的几十年中，降雨量最大值出现在[具体年份1]，达到[X1]mm，而最小值出现在[具体年份2]，仅为[X2]mm，最大值与最小值之间的差值较为显著，反映了兴平市降雨量在年际间的不稳定。这种年际变化可能与全球气候变化、大气环流异常以及地形地貌等多种因素有关。全球气候变化导致极端天气事件增多，使得兴平市的降雨量在不同年份出现较大波动；大气环流异常会影响水汽的输送和降水的形成，从而导致年际降雨量的变化；兴平市的地形地貌也会对降雨产生影响，山区和平原地区的降雨分布可能存在差异。兴平市降雨频率也呈现出一定的特征。年降雨日数平均为[X3]天，其中小雨（日降雨量小于10mm）日数占比较大，约为[X4]%；中雨（日降雨量10-25mm）日数占比约为[X5]%；大雨（日降雨量25-50mm）及暴雨（日降雨量大于50mm）日数占比较小，分别约为[X6]%和[X7]%。小雨日数较多，虽然单次降雨量较小，但由于其出现频率高，对地下水的补给具有一定的累积效应。中雨和大雨在短时间内能够提供较大的降雨量，对地下水的补给作用更为明显。暴雨虽然日数较少，但因其降雨强度大，可能会引发洪涝灾害，同时也会在短时间内大量补给地下水，但这种补给方式可能会导致部分雨水来不及入渗而形成地表径流流失。从季节分布来看，兴平市降雨主要集中在夏季（6-8月），该时期降雨量约占全年的45%-50%。夏季受副热带高压影响，暖湿气流活跃，容易形成降雨天气。且夏季气温高，水汽蒸发量大，大气中的水汽含量充足，为降雨提供了丰富的物质基础。春季（3-5月）和秋季（9-11月）降雨量相对较少，分别约占全年的20%-25%。春季气温逐渐升高，大地回暖，但此时冷暖空气活动频繁，天气变化无常，降雨相对较少。秋季随着冷空气逐渐南下，副热带高压南撤，降雨也逐渐减少。冬季（12-2月）降雨量最少，仅占全年的5%-10%，冬季受蒙古高压控制，气候寒冷干燥，水汽来源少，不利于降雨的形成。降雨的季节分布对地下水的动态变化有着重要影响。夏季降雨集中，入渗补给量大，地下水位会相应上升；而冬季降雨稀少，蒸发量大，地下水位则会有所下降。4.1.2径流特征渭河作为兴平市最重要的河流，其径流变化趋势对地下水补给有着深远影响。通过收集渭河兴平段相关水文站的长期径流数据，运用数据分析方法，深入研究渭河径流的变化趋势及其对地下水补给的影响。渭河径流量年际变化较大，且总体呈下降趋势。以咸阳水文站为例，在过去的几十年中，渭河多年平均径流量为[X8]m³/s，但不同年份间径流量波动明显。如在[具体年份3]，径流量达到[X9]m³/s，而在[具体年份4]，径流量仅为[X10]m³/s。这种年际变化可能是由多种因素共同作用导致的。气候因素是影响渭河径流量的重要原因之一，全球气候变化导致降水分布不均，渭河上游地区降水量减少，使得渭河径流量相应减少。人类活动的影响也不容忽视，随着渭河上游地区经济的发展，用水量不断增加，包括农业灌溉、工业用水和生活用水等，导致渭河径流量减少。渭河上游地区的生态破坏，如植被破坏、水土流失等，也会影响河流的涵养水源能力，进一步减少渭河径流量。渭河径流的年内分配也不均匀，主要集中在汛期（7-10月）。汛期径流量约占全年径流量的60%-70%，这是因为汛期降雨集中，大量雨水汇入渭河，使得径流量大幅增加。非汛期径流量相对较小，约占全年径流量的30%-40%。这种年内分配特征对地下水补给产生了重要影响。在汛期，渭河水位升高，当渭河水位高于地下水位时，渭河河水通过河床渗漏补给地下水，增加地下水的储量，抬高地下水位。在非汛期，渭河水位下降，地下水则向渭河排泄，维持渭河的基流。这种补给和排泄关系使得渭河成为兴平市地下水的重要补给源和排泄通道，对地下水的动态变化产生重要影响。为了更准确地分析渭河径流对地下水补给的影响，还可以通过建立水文模型进行模拟分析。利用SWAT（SoilandWaterAssessmentTool）模型，结合兴平市的地形、土壤、植被等数据，模拟不同径流条件下渭河对地下水的补给量。模拟结果表明，在渭河径流量较大的年份，对地下水的补给量也相应增加；而在径流量较小的年份，补给量则减少。在丰水年，渭河对地下水的补给量可达[X11]m³，而在枯水年，补给量仅为[X12]m³。这进一步说明了渭河径流变化对地下水补给的重要影响。4.2地下水监测井网及观测资料4.2.1地下水监测井网兴平市地下水监测井网的合理布局和有效运行对于准确获取地下水动态信息至关重要。兴平市地下水监测井网由多个监测井组成，这些监测井分布在全市不同区域，包括渭河阶地区和黄土台塬区等不同地貌单元，以全面监测地下水的动态变化。监测井的密度根据不同区域的重要性和地质条件进行合理设置，在城区、工业集中区和农业灌溉区等地下水开采量大的区域，监测井密度相对较高，以更准确地监测地下水水位和水质的变化。在渭河阶地区，由于含水层较厚，地下水开采量大，监测井密度为每[X13]平方公里设置1口监测井；而在黄土台塬区，含水层相对较薄，地下水开采量较小，监测井密度为每[X14]平方公里设置1口监测井。监测井的监测项目包括水位、水温、水质等多个方面。水位监测是最基本的项目，通过定期测量监测井中的水位，了解地下水水位的变化情况。采用自动水位监测仪，能够实时记录水位数据，并通过无线传输技术将数据传输到监测中心，实现对水位的动态监测。水温监测可以反映地下水的热状况，对于研究地下水的运动和补给机制具有一定的参考价值。通过在监测井中安装温度传感器，定期测量地下水的温度。水质监测则包括对地下水的化学组成、酸碱度、硬度、溶解性固体等指标的检测，以评估地下水的质量状况和污染程度。采集地下水水样，送实验室进行化学分析，检测各项水质指标。为了确保监测井的正常运行和数据的准确性，还建立了完善的监测井维护和管理体系。定期对监测井进行检查和维护，包括清洗井管、更换损坏的设备等，保证监测井的畅通和设备的正常运行。对监测人员进行培训，提高其业务水平和操作技能，确保数据的准确采集和记录。4.2.2地下水监测资料兴平市多年来积累了丰富的地下水监测资料，这些资料是研究地下水动态特征的重要依据。监测资料的数据来源主要是兴平市地下水监测站的监测数据，该监测站负责对全市地下水监测井进行定期监测和数据采集。监测数据的时间跨度较长，涵盖了[起始年份]至[结束年份]的多年数据，具有较高的连续性和可靠性。为了保证监测资料的可靠性，在数据采集过程中，严格按照相关标准和规范进行操作。在水位测量时，使用经过校准的水位测量仪，确保测量精度在允许范围内。在水质采样时，遵循样品采集的规范要求，保证样品的代表性和完整性。对采集到的数据进行严格的审核和质量控制，检查数据的合理性和一致性，剔除异常数据。通过与历史数据进行对比，检查数据是否存在异常波动；对同一监测井不同时间的数据进行连续性检查，确保数据的可靠性。对于监测资料的处理，采用了科学的方法。将原始监测数据进行整理和分类，按照时间序列和监测井编号进行归档，方便后续查询和分析。运用数据分析软件，对监测数据进行统计分析，计算地下水水位的平均值、最大值、最小值、标准差等统计参数，了解地下水水位的总体变化特征。采用趋势分析方法，如线性回归分析、滑动平均法等，分析地下水水位的变化趋势，判断地下水位是上升、下降还是保持稳定。通过相关性分析，研究地下水水位与降雨、径流、开采量等因素之间的关系，揭示地下水动态变化的影响因素。4.3地下水位特征分析4.3.1年际水位特征兴平市多年地下水位变化趋势的分析对于了解地下水资源的动态变化和可持续利用具有重要意义。通过对兴平市地下水监测井多年水位数据的收集和整理，运用趋势分析方法，深入探讨地下水位的年际变化趋势及其影响因素。兴平市多年地下水位总体呈下降趋势。以[具体监测井1]为例，在[起始年份]至[结束年份]期间，地下水位从[起始水位1]m下降至[结束水位1]m，累计下降了[下降幅度1]m，平均每年下降[平均下降速率1]m。这种下降趋势在不同区域可能存在差异，在“兴化漏斗区”，由于长期大量抽取地下水，地下水位下降幅度更为明显。兴平化肥厂水源地附近的监测井，地下水位下降幅度达到[下降幅度2]m，平均每年下降[平均下降速率2]m，形成了明显的漏斗状区域。影响兴平市地下水位年际变化的因素主要包括降雨、开采量和渭河径流等。降雨是地下水的重要补给来源之一，当降雨量增加时，入渗补给量相应增加，地下水位会上升；反之，当降雨量减少时，地下水位会下降。根据相关性分析，兴平市地下水位与降雨量之间存在一定的正相关关系，相关系数为[相关系数1]。开采量是影响地下水位的关键因素，随着兴平市经济社会的发展，工业、农业和生活用水对地下水的开采量不断增加，导致地下水位持续下降。在工业集中区，由于大量工业企业抽取地下水用于生产，使得该区域地下水位下降明显。渭河径流对地下水位也有一定影响，在渭河丰水期，河水补给地下水，地下水位会上升；而在枯水期，地下水排泄到渭河，地下水位会下降。为了进一步分析地下水位下降的原因，还可以结合土地利用变化、水资源管理政策等因素进行综合研究。随着兴平市城市化进程的加快，土地利用类型发生了变化，大量耕地被转化为建设用地，导致地表植被覆盖减少，入渗能力下降，间接影响了地下水的补给。水资源管理政策的不完善，对地下水开采的监管不力，也加剧了地下水位的下降。4.3.2年内水位特征兴平市年内地下水位变化规律的研究有助于深入了解地下水的动态变化过程及其与人类活动和自然因素的关系。通过对地下水监测井年内水位数据的详细分析，结合降雨和用水情况，揭示年内水位变化的原因。兴平市年内地下水位呈现出明显的季节性变化规律。一般来说，地下水位在春季（3-5月）开始下降，这主要是因为春季气温回升，蒸发量增大，同时农业灌溉用水开始增加，对地下水的开采量增大，导致地下水位下降。在夏季（6-8月），随着降雨的增加，尤其是集中降雨的出现，入渗补给量增大，地下水位开始上升。夏季是兴平市的雨季，降雨量占全年的45%-50%，大量雨水入渗补给地下水，使得地下水位显著上升。秋季（9-11月），降雨逐渐减少，蒸发量仍然较大，农业灌溉用水也有所减少，但工业和生活用水相对稳定，地下水位开始缓慢下降。冬季（12-2月），降雨稀少，气温较低，蒸发量较小，地下水开采量也相对减少，地下水位相对稳定，但由于前期的开采和蒸发，地下水位仍处于较低水平。以[具体监测井2]为例，在春季3月份，地下水位为[水位1]m，到5月份下降至[水位2]m，下降了[下降幅度3]m；在夏季6月份，地下水位开始上升，到8月份上升至[水位3]m，上升了[上升幅度1]m；秋季9月份地下水位为[水位4]m，到11月份下降至[水位5]m，下降了[下降幅度4]m；冬季12月份地下水位为[水位6]m，到2月份基本保持稳定。通过对降雨和用水情况的进一步分析，可以更清晰地解释年内地下水位变化的原因。在春季，农业灌溉用水量占总用水量的[X15]%左右，大量抽取地下水用于灌溉，导致地下水位下降。而在夏季，虽然农业灌溉用水仍然较多，但降雨的入渗补给量大于开采量，使得地下水位上升。夏季降雨入渗补给量约为[补给量1]m³，而农业灌溉开采量为[开采量1]m³，补给量大于开采量[差值1]m³。秋季随着降雨减少和农业灌溉用水的减少，地下水位开始下降。冬季由于降雨稀少和开采量减少，地下水位相对稳定。4.4地下水水流流场特征绘制兴平市地下水水流流场图是研究地下水流动特征的重要手段。通过收集兴平市不同时期的地下水水位数据，利用专业的绘图软件，如Surfer、ArcGIS等，绘制出不同时期的地下水水流流场图。在绘制过程中，首先根据监测井的位置和水位数据，进行数据插值处理，生成连续的水位面。然后，根据水位面的变化，计算出地下水的水力坡度和流向，绘制出流线，从而得到地下水水流流场图。从流场图可以清晰地看出兴平市地下水的流动方向、速度和水力坡度等特征。兴平市地下水总体流动方向是从西北向东南流动，这与兴平市的地形地貌特征相符，西北地势较高，东南地势较低，地下水在重力作用下由高向低流动。在渭河阶地区，地下水流动速度相对较快，这是因为渭河阶地区含水层厚度较大，岩性颗粒较粗，透水性较好，有利于地下水的流动。根据计算，渭河阶地区地下水的平均流速约为[流速1]m/d。而在黄土台塬区，地下水流动速度相对较慢，黄土台塬区含水层厚度相对较薄，岩性颗粒较细，透水性较差，导致地下水流动速度较慢，平均流速约为[流速2]m/d。地下水的水力坡度也存在一定的空间分布特征。在“兴化漏斗区”，由于长期大量抽取地下水，形成了明显的水位降落漏斗，水力坡度较大。兴平化肥厂水源地附近，水力坡度达到[水力坡度1]，表明地下水在该区域的流动速度较快，向采水点汇聚。而在远离开采区的区域，水力坡度相对较小，地下水流动较为平缓。在黄土台塬区的一些偏远地区，水力坡度仅为[水力坡度2]。为了更准确地分析地下水水流流场特征，还可以结合数值模拟方法进行研究。利用FEFLOW软件建立兴平市地下水水流数值模拟模型，通过模拟不同工况下地下水的流动情况，进一步验证和深化对流场特征的认识。模拟结果表明，在增加开采量的情况下，“兴化漏斗区”的范围会进一步扩大，水力坡度会增大，地下水流动速度会加快；而在增加补给量的情况下，地下水位会上升，水力坡度会减小，地下水流动速度会减缓。4.5地下水化学分析对兴平市地下水化学成分的检测是评估地下水质量和研究地下水污染状况的重要基础。通过在兴平市不同区域采集地下水水样，送专业实验室进行化学分析，检测地下水的主要化学成分，包括阳离子（如Ca²⁺、Mg²⁺、Na⁺、K⁺等）、阴离子（如Cl⁻、SO₄²⁻、HCO₃⁻等）、酸碱度（pH值）、硬度等指标。检测结果表明，兴平市地下水化学组成存在一定的空间分布差异。在渭河阶地区，由于含水层主要为第四系松散岩类孔隙含水层，岩性以沙及沙卵石层为主，地下水与岩石之间的相互作用较强，导致地下水中Ca²⁺、Mg²⁺等阳离子含量相对较高，水的硬度较大。根据检测数据，渭河阶地区地下水的总硬度平均值为[硬度1]mg/L，属于硬水。而在黄土台塬区，含水层主要为黄土质沙质粘土，地下水中的阳离子含量相对较低，水的硬度相对较小，总硬度平均值为[硬度2]mg/L，属于中硬水。通过对检测结果的综合分析，对兴平市地下水水质进行评价。采用《地下水质量标准》（GB/T14848-2017）作为评价标准，对各项指标进行单项评价和综合评价。评价结果显示，兴平市大部分地区地下水水质较好，符合Ⅲ类标准，可作为生活饮用水水源。在一些工业集中区和城市周边，由于受到工业废水、生活污水和农业面源污染的影响，地下水水质存在一定程度的污染。部分区域地下水中的Cl⁻、SO₄²⁻含量超标，pH值偏低，表明可能受到了酸性废水的污染；一些区域地下水中的氨氮、硝酸盐氮含量超标，可能与农业面源污染有关，大量使用化肥和农药，导致氮素通过地表径流和入渗进入地下水。进一步分析地下水污染的来源和程度。工业污染是地下水污染的重要来源之一，兴平市的一些工业企业，如化工、电镀、造纸五、地下水水流数值模拟模型5.1数学模型的理论5.1.1建立模型注意的问题在建立兴平市地下水水流数值模拟模型时，需充分考虑多方面因素以确保模型的准确性和可靠性。数据准确性是建模的基础，其来源广泛，包括地质勘探数据、长期的地下水水位监测数据、抽水试验数据以及气象数据等。地质勘探数据能详细揭示研究区的地层结构、岩性分布等信息，为确定含水层和隔水层的位置、厚度及空间分布提供依据。地下水水位监测数据是分析地下水位动态变化的关键，其时间序列的完整性和准确性直接影响对地下水系统动态特征的把握。抽水试验数据则用于获取含水层的渗透系数、导水系数等重要水文地质参数，试验过程的规范性和数据测量的精度对参数的可靠性至关重要。气象数据，如降雨量、蒸发量等，对于分析地下水的补给和排泄条件不可或缺。在数据收集过程中，要严格把控数据质量，对数据进行筛选和验证。对于水位监测数据，需检查数据的连续性，剔除异常值，可通过与周边监测井数据对比、分析数据变化趋势等方式进行验证。对于抽水试验数据，要确保试验方法符合相关标准，对试验设备进行校准，以保证数据的准确性。边界条件的处理是建模的关键环节，直接影响模型的模拟精度。兴平市地下水系统的边界条件复杂多样，包括自然边界和人为边界。自然边界如渭河等河流边界，其水位和流量的变化对地下水系统有重要影响。在处理河流边界时，需要准确获取渭河的水位和流量数据，这些数据可通过水文站监测资料获取。根据渭河与地下水的水力联系，确定边界类型，若渭河水位高于地下水位，渭河河水补给地下水，可将其视为补给边界；反之，若地下水位高于渭河水位，地下水排泄到渭河，可视为排泄边界。人为边界如地下水开采边界，随着兴平市经济的发展，地下水开采量不断增加，开采井的分布和开采强度对地下水流场产生显著影响。需要详细调查开采井的位置、开采量和开采时间等信息，将开采量作为源汇项加入模型中，以准确模拟人为开采对地下水系统的影响。参数选取的合理性决定了模型能否准确反映地下水系统的实际情况。渗透系数、给水度、储水系数等是模型的关键参数，其取值需要综合多种方法确定。渗透系数可通过抽水试验、实验室测试以及经验公式估算等方法获取。在进行抽水试验时，要根据含水层的类型和特征选择合适的试验方法，如稳定流抽水试验或非稳定流抽水试验，并结合试验数据运用相应的计算公式计算渗透系数。实验室测试可对岩芯样品进行分析，获取岩石的渗透性能参数，但需注意样品的代表性。经验公式估算则是根据研究区的地质条件和已有的研究成果，选择合适的经验公式进行计算。给水度和储水系数的确定可参考相关的地质资料和前人的研究成果，同时结合研究区的实际情况进行调整。通过敏感性分析，评估不同参数对模拟结果的影响程度，对敏感性较高的参数进行更精确的确定，以提高模型的模拟精度。5.1.2地下水系统概念模型构建兴平市地下水系统概念模型是建立数值模拟模型的重要基础，它旨在对复杂的地下水系统进行合理简化和抽象，以便于后续的数学描述和数值模拟。兴平市地下水系统概念模型的构建基于对研究区地质、水文地质条件的深入了解。在确定模拟区域范围时，充分考虑兴平市的地形地貌、地质构造以及地下水的补径排条件。以兴平市的行政边界为基础，结合地形等高线和地质构造特征，将模拟区域范围确定为东至[具体边界位置1]，西至[具体边界位置2]，南至[具体边界位置3]，北至[具体边界位置4]。该范围涵盖了兴平市主要的含水层分布区域，包括渭河阶地区和黄土台塬区，能够全面反映兴平市地下水系统的特征。模拟区域的边界条件分为自然边界和人为边界。自然边界方面，渭河作为兴平市最重要的地表水体，是地下水系统的重要边界。渭河与地下水存在密切的水力联系，在丰水期，渭河水位高于地下水位，渭河河水通过河床渗漏补给地下水；在枯水期，地下水位高于渭河水位，地下水向渭河排泄。根据这一特点，将渭河边界概化为一类流量边界，通过监测渭河的水位和流量数据，确定边界上的补给或排泄量。研究区北部的黄土台塬区与周边地区的地下水存在水力联系，将其边界概化为二类流量边界，根据区域水文地质条件和地下水流向，确定边界上的流量值。人为边界方面，兴平市的地下水开采活动较为频繁，开采井分布广泛。将开采井作为源汇项处理，根据实际的开采量和开采时间，确定开采井在模型中的位置和开采强度。兴平市含水层结构较为复杂，主要包括第四系松散岩类孔隙含水层和新近系碎屑岩类裂隙孔隙含水层。第四系松散岩类孔隙含水层广泛分布于渭河阶地区和黄土台塬区，根据含水层的岩性、结构和水力特征，进一步将其划分为潜水含水层和承压含水层。潜水含水层主要分布在渭河阶地区的河漫滩和一级阶地，岩性以中、细沙和沙质粘土为主，厚度一般为10-20米，富水性较好，单井出水量一般在1000-3000立方米/日。承压含水层主要分布在渭河阶地区的二级、三级阶地和黄土台塬区，岩性以沙及沙卵石层为主，厚度较大，一般为30-80米，富水性中等，单井出水量一般在500-1000立方米/日。新近系碎屑岩类裂隙孔隙含水层主要分布在北部黄土台塬区的深部，岩性为砂岩、泥岩互层，裂隙发育程度较低，富水性较差，单井出水量一般小于500立方米/日。在概念模型中，对含水层进行分层处理，每层赋予相应的水文地质参数，如渗透系数、给水度、储水系数等，以准确反映含水层的特性。5.2地下水系统数学模型5.2.1潜水非稳定流连续性方程潜水非稳定流连续性方程是描述潜水含水层中地下水运动的重要方程，其推导基于质量守恒定律和达西定律。考虑一个在空间直角坐标系下的潜水含水层单元，该单元在x、y、z方向的长度分别为Δx、Δy、Δz，含水层的厚度为h，渗透系数在x、y、z方向分别为Kx、Ky、Kz，给水度为μ。根据达西定律，在x方向上，单位时间内通过垂直于x轴的单位面积的流量qx为：qx=-Kx\frac{\partialh}{\partialx}同理，在y方向和z方向上的流量qy和qz分别为：qy=-Ky\frac{\partialh}{\partialy}qz=-Kz\frac{\partialh}{\partialz}在Δt时间内，流入和流出该单元的水量差应等于单元内储存水量的变化。流入单元的水量为：Q_{in}=qx|_{x}\Deltay\Deltaz+qy|_{y}\Deltax\Deltaz+qz|_{z}\Deltax\Deltay流出单元的水量为：Q_{out}=qx|_{x+\Deltax}\Deltay\Deltaz+qy|_{y+\Deltay}\Deltax\Deltaz+qz|_{z+\Deltaz}\Deltax\Deltay单元内储存水量的变化为：\DeltaW=\mu\Deltax\Deltay\Deltah其中，\Deltah为在\Deltat时间内潜水位的变化。根据质量守恒定律，Q_{in}-Q_{out}=\DeltaW，将上述流量表达式代入可得：-\left(Kx\frac{\partialh}{\partialx}\right|_{x+\Deltax}-Kx\frac{\partialh}{\partialx}\left|_{x}\right)\Deltay\Deltaz-\left(Ky\frac{\partialh}{\partialy}\right|_{y+\Deltay}-Ky\frac{\partialh}{\partialy}\left|_{y}\right)\Deltax\Deltaz-\left(Kz\frac{\partialh}{\partialz}\right|_{z+\Deltaz}-Kz\frac{\partialh}{\partialz}\left|_{z}\right)\Deltax\Deltay=\mu\Deltax\Deltay\Deltah当\Deltax、\Deltay、\Deltaz、\Deltat趋于0时，对上式进行极限运算，利用偏导数的定义，得到潜水非稳定流连续性方程：\frac{\partial}{\partialx}\left(Kx\frac{\partialh}{\partialx}\right)+\frac{\partial}{\partialy}\left(Ky\frac{\partialh}{\partialy}\right)+\frac{\partial}{\partialz}\left(Kz\frac{\partialh}{\partialz}\right)+W=\mu\frac{\partialh}{\partialt}其中，W为单位体积含水层中源汇项的流量，如降雨入渗补给、蒸发排泄、人工开采等。该方程的物理意义为：方程左边第一项表示x方向上由于水力坡度引起的流量变化，第二项表示y方向上的流量变化，第三项表示z方向上的流量变化，这三项之和反映了地下水在三维空间中的流动情况；第四项W表示源汇项对流量的影响；方程右边表示单位时间内潜水含水层中储存水量的变化，与潜水位的变化率成正比。5.2.2潜水非稳定流的定解条件潜水非稳定流的定解条件包括初始条件和边界条件，它们是确定数学模型唯一解的必要条件。初始条件是指在初始时刻（t=0），潜水含水层中各点的水位分布情况，即：h(x,y,z,0)=h_0(x,y,z)其中，h_0(x,y,z)为已知的初始水位函数，可通过研究区在初始时刻的地下水位监测数据确定。在兴平市地下水数值模拟中，收集模拟起始时刻各监测井的水位数据，利用空间插值方法，如克里格插值法，得到整个模拟区域的初始水位分布，作为初始条件代入模型中。边界条件分为三类：一类边界条件（Dirichlet边界条件）是指在边界上给定水位值，即：h(x,y,z,t)|_{\Gamma_1}=h_1(x,y,z,t)其中，\Gamma_1为一类边界，h_1(x,y,z,t)为边界上已知的水位函数。在兴平市，若将渭河边界视为一类边界，当渭河水位已知时，可将渭河水位作为边界上的水位值。二类边界条件（Neumann边界条件）是指在边界上给定流量值，即：-K\frac{\partialh}{\partialn}|_{\Gamma_2}=q_2(x,y,z,t)其中，\Gamma_2为二类边界，K为边界处的渗透系数，\frac{\partialh}{\partialn}为水位在边界法向方向上的导数，q_2(x,y,z,t)为边界上已知的流量函数。在处理研究区北部黄土台塬区与周边地区的地下水水力联系边界时，若已知边界上的流量值，可将其作为二类边界条件。三类边界条件（Cauchy边界条件）是指在边界上给定水位和流量的线性组合，即：K\frac{\partialh}{\partialn}|_{\Gamma_3}+\alp