三维裂隙网络渗流计算的边界元法和程序.pdfVIP

第 4 卷 　第 2 期 中国水利水电科学研究院学报 Vol4 　No 2 2006 年 6 月 Journal of China Institute of Water Resources and Hydropower Research June ,2006 ( ) 文章编号 2006 三维裂隙网络渗流计算的边界元法及程序 李新强 ,陈祖煜 ( 中国水利水电科学研究院 岩土工程研究所,北京　100044) 摘要 :岩体力学性质受结构面控制 ,表现出强烈的各向异性。同时工程地质师多年的实践认为岩石结构面是具有 统计性的 ,因此岩体结构面网络模拟成为研究岩石力性质的新方向。其中裂隙网络渗流模拟 ,特别是三维裂隙网 络的渗流计算是结构面网络模拟研究工作的重点和难点。又由于对于具有随机性裂隙网络渗流的计算 ,需要多次 形成裂隙网络样本 ,并对每个样本进行多次计算 ,因此对渗流计算比确定性渗流提出更高的要求。本文基于将单 个裂隙看作为一个不可压缩的二维各向同性多孔介质的假定 ,运用加权剩余法 ,推导了单裂隙渗流的边界单元法 公式。通过将裂隙之间的交线看作为内部边界 ,推导了多裂隙相交时的边界单元法公式。利用 FORTRAN 90 语言 编制了相应的计算程序 BIEMNF 。为验证程序的正确性 ,利用计算实例进行了考核。 关键词 :岩石水力学 ;渗流 ;边界元法 ;裂隙网络 中图分类号:P641135 文献标识码 :A 由于岩体的工程地质特性、以及所处的复杂多变的地质条件等, 国内外学者和地质师非常重视现场 试验测试岩体的渗透性。然而 ,各种试验方法不仅耗费很大 ,时间长 ,而且独自难以准确有效测定裂隙 岩体的渗透性。裂隙岩体由结构面和岩块组成 ,大量的试验和工程实践使人们认识到岩块渗透性相对 于结构面而言相当微小 ,岩体渗流决定于结构面状态。通过对结构面的观测和统计分析证实 ,岩体内绝 大多数结构面的分布不是杂乱无章的 ,而是服从一定的统计规律 , 因此可以利用概率统计的方法来处 理。使用岩体力学和工程地质学科内新兴的岩体结构面网络模拟技术 ,可以将这些结构面反映到岩体 空间中 ,生成结构面网络模型。这种模型使研究随机裂隙岩体的力学性质成为可能。 利用随机结构面网络模拟的手段 ,辅之以现场试验的方法 ,联合确定裂隙岩体的渗透性 ,可较为有 效解决裂隙岩体渗透性的定量确定问题。对于随机形成的裂隙网络 ,利用抽样分析的方法 ,对每一样本 的渗流计算是必须和必要的工作 , 同时是渗流网络模拟的关键环节。通过对现有包括有限元法、线素 法、图论法、有限差分法、离散元法 ,DDA 法和流形元法等多种数值方法在求解裂隙渗流的应用情况的 对比分析[1 ] ,可知各方法均能有效求解已知裂隙网络的渗流问题 ,不过在裂隙网络的处理工作量和难度 上存在差别 , 由于所作假定条件的合理程度不同,求解的结果也会有所差别。引入一种高效、快速和易 于对数据自动准备和处理的计算方法 ,是决定利用离散裂隙网络渗流评价裂隙岩体渗透性的实用性基 础工作 ,需要在各细节问题上作深入研究。 1 　裂隙渗流的数学模型 11 　裂隙网络渗流的基本假定 ( ) ( ) 1 不考虑垂直于裂隙壁面的流动 ;沿裂隙面渗流各向同性 ,并符合达西定律。 2 裂隙隙宽为随机 ( ) ( ) 值 ,且隙宽为等宽的。 3 忽略水的压缩性 ,暂不考虑裂隙变形对渗流的影响。 4 不同组结构面之间的 交叉处引起的局部水头损失忽略不计。 12 　渗流控制方程和定解条件 　与连续介质相似 ,对于一域为 Ω,边界为 Γ的标量场 h ,水流满足 Laplace 方程 收稿日期 ( ) 作者简介 :李新强 1972 - , 山东滕州人,博士 ,工程师 ,主要研究方向为