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

岩体力学性质受结构面控制，表现出强烈的各向异性。同时工程地质师多年的实践认为岩石结构面是具有统计性的，因此岩体结构面网络模拟成为研究岩石力性质的新方向。其中裂隙网络渗流模拟，特别是三维裂隙网络的渗流计算是结构面网络模拟研究工作的重点和难点。又由于对于具有随机性裂隙网络渗流的计算，需要多次形成裂隙网络样本，并对每个样本进行多次计算，因此对渗流计算比确定性渗流提出更高的要求。本文基于将单个裂隙看作为一个不可压缩的二维各向同性多孔介质的假定，运用加权剩余法，推导了单裂隙渗流的边界单元法公式。通过将裂隙之间的交线看作为内部边界，推导了多裂隙相交时的边界单元法公式。利用FORTRAN 90语言编制了相应的计算程序BIEMNF。为验证程序的正确性，利用计算实例进行了考核。

Boundary element method for 3-D fracture network seepage flow and its programing

Rock hydraulics is a very important aspect of the rock mechanics. For a long time, the rock engineering problem has been solved with multi-porous media model. However, the behavior of rock shows strong anisotropy which is mainly caused by fractures in it. Geo-engineers have made it clear that the distribution of rock fractures is stochastic but of statistical feature. It is possible to study the problem of rock mechanics by means of simulation with fracture network method. 3-D fracture network seepage flow calculation is the most important and difficult part in fracture network flow simulation. When determining the rock mass conduction tensor, the stochastic samples are needed to form many times, and each sample should be calculated several times. So based on the postulation that single fracture is an incompressible plane of isotropy multi-porous media, the Boundary Element Method of single fracture is derived by using the Weighted Residual Method. Supposing that the intersection of multi-group fractures is treated as the internal boundary, a Boundary Element Method of fracture network seepage flow is obtained. To perform the procedure, the authors have developed a corresponding program named as BIEMNF in FORTRAN 90, and then check it with some examples.