FAC方法在地下水数值模拟中的应用.pdfVIP

水 利 学 报 年 月 %’()( *’+,-. 第 卷 第 期 !# $$ /0 $$ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） DD313/D !# $$1 $/4312 !# 方法在地下水数值模拟中的应用 $ $ ! / 江思珉 ，朱国荣 ，王浩然 ，孙斌堂 （ 南京大学 地球科学系，江苏 南京 ； 南京大学 软件学院，江苏 南京 ； $E !$3/ ! E !$3/ / E 中国冶金地质工程勘查总局山东局，山东 济南 !D$2 ） 摘要：为了解决有限单元法和有限差分法在模拟条带状水文地质体时可能产生海量数据和耗费大量计算机时的困 难，本文提出采用基于区域分解原理的快速自适应组合网格方法 法（ ），讨论了 F-G F7=H -I7JH8KL GA6JA=8HL MB8I F-G 法建立和求解地下水模型的基本原理，并将其应用于一个矿山边坡防渗墙体防渗性能的研究。算例采用三层网格 建立模型并求解，准确描述出防渗帷幕墙体内部地下水运动的规律，为即将开始的防渗墙设计提供了参考依据。 关键词：地下水模型；快速自适应组合网格方法；区域分解 中图分类号： 文献标识码： N#2 $ - $ 研究背景 有限单元法（ ）和有限差分法（ ）在地下水模型研究中是比较成熟的模拟技术，模拟一般性 F+O FPO ［ ］ 的地下水运动时，都可以比较准确地再现模拟对象的细节 $ Q / 。然而，当研究区内存在诸如断裂带、岩 墙、软弱夹层等线状展布的地质体（以下统称线形地质体）时，用 F+O 或 FPO 模拟其中的地下水运动是 比较困难的。计算中，对于可视为隔水边界的线形地质体，当其中的地下水运动能够忽略的情况下，可 ［，］ / 2 以采用结点双编号或内部隔水边界等处理方法来解决 ，否则就应对这些线形地质体精细离散，并为 ［］ D 达到精度要求提供必须的精细网格 。这样，为达到对线性地质体研究精度的要求，整个区域的离散工 作量将会大大增加，同时计算所需机时也将大幅度增长。 如果能够只对这些线形地质体进行必要的精细离散，同时对其他研究区域采用粗网格，则有可能避 免上述的问题。然而这种非拟一致网格（ ）的离散处理可能产生求解的困难。首先， RA;1CA;SAB68; MB8I 空间的非拟一致离散违背了各种地下水数值方法的离散宗旨，其次方程的求解方法可能因离散尺度相 对变化过大而