水气二相渗流与双重介质变形流固耦合数学模型.pdfVIP

水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期 文章编号: 0559-9350(2005)04-0405-08 水气二相渗流与双重介质变形的流固耦合数学模型 1 2 1,3 4 刘晓丽，梁冰，王思敬 ，李宏艳 (1.清华大学 水利水电工程系，北京 100084；2.辽宁工程技术大学 力学与工程科学系，辽宁 阜新 123000；3.中国科学院 地质与地球物理研究所，北京 100029；4.北京师范大学 环境科学研究所，北京 100875) 摘要：为了真实反映地质岩体系统内孔隙、裂隙对发生于其中的物理过程的控制与影响，基于岩体渗流水力学和 多相渗流力学理论，将工程地质体简化为孔隙-裂隙双重介质，建立了水气二相渗流与双重介质变形的流固耦合数 学模型，采用Galerkin有限元方法推导了相应的有限元方程，编制出三维有限元计算程序。算例表明：双重介质 模型能够比较准确地描述地质体的结构；在进行注水井数值实验模拟时，考虑耦合与不考虑耦合时的孔隙压力的 差别随时间的推移会越来越大。 关键词：双重介质；水气二相流；流固耦合；有限元 中图分类号：TU441 文献标识码：A 近几年，环境岩体工程发展迅速。这是一门交叉学科，它覆盖了在大气圈、生命圈、水圈、岩石圈及 地质微生物圈等各种环境下土和岩石及其相互作用的问题，其中岩石圈与水圈相互作用问题的研究一直受 到国内外学者的关注[1～5] 。这方面的研究主要是应用岩石力学、渗流力学、环境流体力学的基本理论来解 决工程实际问题，它涉及到资源、环境、水利水电工程、地下工程、核废料储存、水库诱发地震、地球动 力学等诸多领域。马尔帕塞(Malpasset)薄拱坝、瓦依昂(Vajont)水库等诸多工程由于渗流作用而发生破 坏的事实，使人们认识到渗流分析对工程安全及经济合理性分析的重要性。地下流体和地壳构造运动的相 互联系、相互作用更是构成了地学研究的前缘和热点，这使得对包括多相流体(油、水、气)在内的地下流 固耦合问题的研究变得更为重要和迫切。 随着岩体工程规模越来越大、复杂程度越来越高，人们对岩石材料的认识也越来越深。岩体是富含一 系列软弱结构面和裂隙的天然地质体，成岩作用以及成岩环境的差异和随机性造成了岩体宏观结构和微观 结构上的非均质性、各向异性和不连续性，再加上不同地应力环境及其变化，形成了岩体中大量方向各异 的岩体结构面。在岩体中存在的结构面、节理、裂隙对岩体的强度和变形以及渗透性起着至关重要的作用。 岩体的破坏与失稳总是伴随着原有裂隙面的的演化、扩展和贯通而发生，因而，岩体强烈地受到内部节理 裂隙面性态及分布方式的影响、制约和控制。从工程力学参数上考虑，裂隙的存在与扩展使得岩体的变形 模量和强度参数降低、渗透系数增大、损伤加剧并呈强烈的各向异性[6]。 基于岩土体结构的复杂性，本文在渗流分析过程中，将孔隙和微裂隙视为孔隙结构，而将结构面、节 理、裂隙统一视为裂隙结构，建立了水、气二相渗流与孔隙-裂隙双重介质变形的流固耦合数学模型。 1 水、气二相渗流与双重介质变形的流固耦合模型 收稿日期：2004-06-29 基金项目：国家自然科学基金资助项目；国家“十五”攻关项目(2001BA803B0404) 作者简介：刘晓丽(1978-)，男，山东泰安人，博士生，主要从事地下流固耦合理论与应用及环境流体力学方面的研究。 E-mail:liuxl04@mails.tsinghua.edu.cn 1 水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期 岩土体中存在大量不连续的地质结构面，如节理、断层、裂隙、岩层层面、片理面等，由这些稀疏的 结构面形成的裂隙网络和