基于离散裂隙网络模型的节理岩体渗透张量及特性分析

Permeability tensor and seepage properties for jointed rock masses based on discrete fracture network model

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明:单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。

Joint systems in rock masses are geometrically complex; and investigation of the hydraulic properties is important for understanding the hydraulic behavior of jointed rock masses and stability analysis of rock engineering. A discrete fracture network (DFN) model was designed using C++language based on VC++6.0 platform to realize a fracture network for seepage analysis. The permeability coefficients, possibility of equivalent permeability tensors of the jointed rock masses were examined by simulating flow through DFN models of varying joint properties. A directionality coefficient was defined to quantitatively assess the principal flow direction in jointed rock masses. The results show that the permeability coefficients are anisotropic according to permeability tensor analysis for rock masses with only one set of joints. The degree of obliquities of joint plane affects obviously the principal direction of permeability. The permeability coefficients increase logarithmically with the increase of joints connectivity rates. In the scenarios where two sets of joints existed, the anisotropic properties vary with the joints geometrical distribution. The permeability anisotropy in a fracture network with two mutually perpendicular sets of joints still exists; and however, the directionality coefficients are examined smaller. For the dip angles whose probability distribution can be effectively by normal function, the permeability coefficients increase with increasing of the standard deviation. The principal direction of permeability runs similarly along the direction of angle bisector between two sets of joints.