一种处理地下水流模拟中自由水面问题的算法研究.pdf

第=篇技术方法 一种处理地下水流模拟中自由水面问题的算法 林琳杨金忠 (武汉大学水资源与水电工程科学目家重点实验室武汉430072) 摘要本文采用t—y二维有限元分析和垂向有隈差分析相结合的方法建立地下水三维渗流筒化模型，提出有 效的数值算法解决自由水面、含水层琉干和重新湿润的问题。做抽水井抽未和回水算例．多方面验证了简化模型 及教值算法的准确性和稳定性。水量平衡计算增加了计算结果的可信性。 关键词地下水流摸型 自由水面问题网椿变动 1引言 二维和三维地下水流数学模型都需要解决潜水面变动问题，大多数计算方法是根据要解决的特定实际问 题提出不同的求解策略，如大坝渗流、区域地下水流模拟和非饱和土柱含水量变化”一J。目前，完全地下水 流模型在实际应用中仍有很多缺点和问题”，。 三维地下水流有限差分模拟软件VisualMODFLOW把计算单元按性质分为无渗流计算单元、变水头单元 和定水头单元三类，由用户事先指定”J。当一个计算单元中的水头低于该单元底面标高时，该单元处于非饱 和状态，MODFLOW将其定义为无效单元，不参与计算，这种处理方法不是很完善，因而一般建议尽量避免出 现干牯单元。Visual 指定单元的二次湿润参数，判断这个单元是否参与下一时刻水流计算，但这常产生不收敛或不稳定的解。 孙讷正”1曾提出一种地下水三维流动简单算法用于矿床涌水量的预测计算，用一个长方体把渗流区域 包括起来，剖分成网格，在迭代过程中采用抛弃单元的做法，追踪动边界的准确位置，渗流区域模拟计算的 节点和单元数目总在变化．这种方法理论上比较严谨，但在程序中处理起来很麻烦，水流方程求解复杂。 了移动网格技术，可以在复杂地层和非均质介质中计算区域位置，但仅允许网格上酃变动。FEFLOW【．1采用 graphicDam)方法处理潜水面引起的网格移动。 本文在z—y二维方向上采用有限单元法分析水量平衡，在：向上则采用有限差分法，建立地下水流三 维简化模型，并提出一种简单的处理方法解决地下水流区域随潜水面变动的问题。 2地下水三维渗流简化模型 地下水三维渗流方程为： i (-)¨J P-i等=蠹(丘蓑)一s2瓦I～瓦J” 式中：n为比储水系数；H为总水头；‘为时间；≈、耳为空间坐标z、hz；E为饱和水力传导度；s为源 汇项，垂向坐标=向下为正。 2，3，4)表示各分层底板及顶板垂向坐标。 程，为： p。e等=磊0、K，e等)+毒(墨口筹)+El。一v,I。，一s (：) 第十章数学模拟技术 地表毋 潜水面H 底板日 图l地下含水层渗流系统 一个垂向分层的兰维水均衡单元，如图2。丘^和M。分别表示水均衡单元的长度、宽度和垂向厚度；虬i-‘、 口：2、g：2、口：。、订。、％t-。分别表示水均衡单元左、右、前、后、上、下边界水流通量，质量守恒方程为： (3) 金有限元法进行空间离散，得到维数相同的各垂向分层有限元矩阵，矩阵形式如图3中工。，…．“，同时层 与层之间在垂直方向上对应的节点用有限差方法分析．物理意义上可解释为节点控制区域的水量均衡，即在 限元矩阵如图3，“O”代表层与层间联系，空白区域为零。简化模型较完全三维模型的优点是程序编制容 易、节省内存、对离散方程组存在有效的解法，便于处理潜水面的位置变化。 A…AO 惑， i Ll； 、 △…A O 7lq【五‘