三维裂隙网络非稳定渗流数值分析

研究裂隙岩体中渗透水流随时间、边界条件的变化规律，从而客观反映渗透水流情况，具有重要的意义。本文针对分布特别稀疏的三维主干岩体裂隙网络，以渗流区内的水量平衡原理为依据建立基本方程，推导三维岩体裂隙网络非稳定渗流的基本方程；并将三维裂隙网络渗流理论及有限元数值方法运用于主干岩体裂隙网络渗流分析中。基于裂隙岩体网络非稳定渗流数学模型，编制相应的程序。对于工程实例，根据实测情况利用Monte-Carlo模拟出岩体中二级裂隙网络分布情况进行分析，然后进行了三维非稳定渗流分析。计算结果充分体现了岩体中裂隙网络的储水性，当边界条件发生变化时，裂隙岩体中水头分布变化存在一定的滞后性；并表明主干裂隙的强导水性，使得裂隙岩体中的主干裂隙的渗流分析成为渗流分析中的重点。

Numerical analysis of 3-D unsteady seepage through fracture network in rock mass

It is of great importance to study the fluid flow in the fracture network varying with time and the boundary condition, and to simulate the distribution of seepage in the rock mass. For the three-dimensional main fracture network that distributing sparsely in rock mass, the basic equation of three-dimensional unsteady seepage is deduced based on the principle of water balance. The corresponding finite element method is also proposed and applied to the analysis of three-dimensional unsteady seepage through main fracture network in rock mass. The codes are compiled following the mathematic model of the three-dimensional unsteady seepage in the fracture network. An engineering example is analyzed by using the codes. The distribution of the secondary fracture network in the rock mass is simulated by the Monte-Carlo method. The calculated results show that the main fracture networks dominate water flow within rock mass, and the hysteresis of the seepage head distribution can be seen in the example when the boundary condition is varying.