多尺度有限元法在非均质多孔介质地下水流问题的应用-港口、海岸及近海工程专业论文.docxVIP

万方数据 万方数据 独创声明 本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研 究成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他 人已经发表或撰写过的研究成果，也不包含其他教育机构的学位或证书使用过的 材料。与我一同工作的同志对本研究所做的任何贡献均己在论文中作了明确的说 明并表示谢意。 学位论文作者签名:在)说红签字日期:训 年 j 月 1 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定，.并同意以下 事项: 1、学校有权保留并向国家有关部门或机构送交论文的复印件和磁盘，允许论 文被查阅和借阅。 2 、学校可以将学位论文的全部或部分内容编入有关数据库进行检索，可以采 用影印、缩印或扫描等复制手段保存、汇编学位论文。同时授权清华大学中国 学术期刊(光盘版)电子杂志社用于出版和编入 CNKI ((中国知识资源总库))，授权 中国科学技术信息研究所将本学位论文收录到《中国学位论文全文数据库比(保 密的学位论文在解密后适用本授权书) 学位论文作者签名:离现红 导师签字:仰眼 签字日期 :}OI叶年 6月 1 日 签字日期:明仔年 j 月， I 日 多尺度有限元法在非均质多孔介质地下水流问题的应用 摘 要 利用数值模型对地下水流动和溶质运移问题进行模拟的方法以其有效性、灵 活性和相对廉价性逐渐成为地下水领域一种重要的研究方法，并受到越来越多的 重视和广泛的应用。然而，自然界中地下水含水系统由于地质条件及周围环境的 影响，其结构变得非常复杂，呈现出高度的非均质特性，而这种非均质特性通常 体现在渗透系数的多尺度波动上。在应用传统有限单元或有限差分等方法对这类 大区域多尺度问题进行数值模拟时，通常存在因占用内存过大或耗时过长而难以 求解的问题。因此，各国学者致力于寻求既可以减少计算资源，又能保证计算精 度的新数值方法来处理非均质含水层的多尺度问题。近年来提出的多尺度有限单 元法（MsFEM）就是其中十分有效的一种方法，它可以在较粗的网格单元取得问 题的合理解，克服了传统有限单元法在求解此类问题时因网格过小过密而导致的 自由度太大无法求解的难题。本文的研究也是以该多尺度有限单元法的基本原理 为基础的，主要内容如下： (1)本文首先系统的总结介绍了如今地下水数值模拟方面所使用的各类求解 方法以及相应的研究成果，通过分析发现虽然现在已经有许多的计算方法及商用 软件可以模拟地下水流问题，但是多尺度求解算法并没有形成系统完整的程序， 现有的软件也没有应用这种多尺度方法的，所以开发编写多尺度算法的程序代码 仍然是十分必要的； (2)其次，本文简要介绍了 OpenGeoSys 软件的一些相关内容，然后以 OGS 软 件为基础，用 C++语言编写了新的源程序，其中主要是设计了多尺度基函数的程序 实现过程，并将该多尺度有限单元法以新功能模块的形式引入到 OGS 源程序中； (3)最后，本文建立了各向异性非均质多孔介质中饱和二维稳定流的数值模型 以及各向同性非均质多孔介质中饱和稳定流的二维和三维数值模型，其中的水力 梯度是通过在边界上设置不同的水头值来实现的，然后用新编写的程序对模型进 行数值计算，并将计算结果与传统方法的计算结果进行了比较，结果表明新编写 的程序既能节省计算量又能保证计算精度。 本文通过对 OpenGeoSys 软件源程序进行拓展，构建了一种求解饱和非均质 多孔介质中地下水流问题的高效的计算程序，既节省计算量又具有较高的精度， 在多尺度非均质多孔介质地下水模拟中具有很大的应用潜力。开发利用多尺度有 限元方法求解地下水流数学模型的模拟软件系统，具有重要的理论意义和实用价 值。 关键词：多孔介质；数值模拟；多尺度有限元法；OpenGeoSys 软件 Application of Multiscale Finite Element Methods in Groundwater Flow of Heterogeneous Porous Media Abstract The numerical modles are becoming an important method for the simulation of the groundwater flow and solute transport problems ,and now it have received more and more attention and been used widely for its’ effectiveness,flexibility and cheapness.However the groundwater aquifer systems in the por