粗糙节理网络模拟及裂隙岩体水力耦合特性研究

【摘要】： 众所周知,节理是控制岩体水力学行为的重要因素,研究节理的变形行为及渗透行为是对裂隙岩体工程安全评估的重要组成部分。随着地下开采、隧道、边坡及核废料地下处置库等岩体工程的兴建,对裂隙岩体水力耦合研究的需求越来越迫切。然而目前有关节理变形行为及渗透行为的描述和数值模拟,均局限于光滑、平行板模型,节理产状参数相互独立且节理处于常法向应力边界等假设；在使用裂隙岩体水力学模型时,常常忽略了其前提条件。针对目前研究的不足,结合国际合作项目DECOVALEX,以核废料地下处置库安全评估为最终目的,在国内外研究的基础上,采用离散裂隙网络模型方法,通过离散元数值计算(DEM)分析了节理粗糙性、非贯通节理的弱化作用、产状参数相关性以及应力状态等多因素对裂隙岩体水力耦合特性的影响,为分析实际裂隙岩体的水力学特性提供研究方法,其具体内容如下： 1)提出了一种新的粗糙节理网络模拟方法—SAW法,将粗糙节理视为不断行走的SAW链,通过控制SAW链的一些基本参数(如键长、凸起高度等)来控制粗糙节理形态。本文详细讲述了SAW法的主要思想、实现程序、基本参数的取值以及其使用条件,最终给出了较为系统的粗糙节理网络模拟方法。且应用SAW,研究了不同取样长度下,描述节理粗糙性的各统计参数与JRC之间的相关关系。 2)考虑节理模型的广泛适用性以及其参数取值的方便性,引入适用于任何边界条件下的粗糙节理本构模型—SA模型,将其嵌入UDEC自定义节理本构模型,并通过数值单轴压缩试验和直剪实验验证了该自定义节理本构模型的正确性,为水、力学计算提供基础。 3)基于SAW,生成单条变隙宽、粗糙节理,对其进行数值渗透实验,分析了其水流规律,得出：节理等效水力隙宽与力学隙宽均值呈指数关系,与节理粗糙度系数JRC呈负指数关系。拟合三者的关系,将其与BB模型修正立方定律进行对比分析,得到了较好的相似性,从数值模拟的角度验证了BB模型的修正立方定律,并将修正的立方定律与SA节理模型联合使用,最终给出了粗糙节理水力学本构模型,为研究裂隙岩体水力学特性打下基础。 4)为考虑非贯通节理对裂隙岩体力学强度的弱化作用,引入一损伤变量D,将其定义为块体内所有非贯通节理的实际长度总和与延长至块体边界的长度总和之比。并给出了各块体内损伤变量D的实现程序,然后将含有非贯通节理的岩块弹性模量进行折减,分析非贯通节理对裂隙岩体水力学特性的影响规律,得出：非贯通节理的弱化作用对裂隙岩体力学REV的影响较小,对岩体等效弹性模量有较大的消减作用；其弱化作用对裂隙岩体水力耦合REV的影响较小,对裂隙岩体渗透特性有消弱作用；尤其对岩体渗透率的最大值影响较大,对渗透率的最小值和渗透主方向的影响较小。在实际应用中,应考虑非贯通节理弱化作用的影响。 5)研究不同分布规律的节理迹长与隙宽的相关性,重点讲述了迹长与隙宽均符合幂律分布时的相关性,在此基础上,对比分析了常隙宽情况下、变隙宽且不相关情况下以及变隙宽且相关情况下裂隙岩体的水力耦合特性,得出：考虑迹长与隙宽相关性对裂隙岩体的渗透REV、渗透率张量均有很大影响；只考虑隙宽的变化对裂隙岩体渗透REV的影响较小,但对其渗透率张量的影响较大。因此,在实际工程中,应根据地质调查结果对其分别进行研究。 6)基于上述研究的基础,结合国际合作项目DECOVALEXⅢBMT 2,计算了Sellafield处核废料处置库各个区域内裂隙岩体的等效渗透系数张量,并研究了不同地应力状态下,其渗透系数张量的变化规律,得出：A)区域1裂隙岩体的等效渗透系数张量主渗透系数K1和K2分别为3.89×10-6m/s和2.27×10-6m/s,其主渗透方向为-15.32°；B)区域2裂隙岩体的等效渗透系数张量主渗透系数K1和K2分别为5.76×10-6m/s和1.06×10-6m/s,其主渗透方向为-32.03°；C)断层带裂隙岩体的等效渗透系数张量主渗透系数K1和K2分别为5.75×10-6m/s和1.97×10-6m/s,其主渗透方向为-31.33°；D)区域1裂隙岩体等效渗透系数张量随深度的变化规律为：Kxx=8.209×10-6·d-0.536, Kxy=3.081×10-6·d-0905,Kyy=2.205×10-6·d-0.635；E)区域2裂隙岩体等效渗透系数张量随深度的变化规律为：Kxx=4.275×10-6·d-0.255,Kxy=5.299×10-6·d-0.325,Kyy=5.971×10-6·d-0.115；F)断层带裂隙岩体等效渗透系数张量随深度的变化规律为：Kxx=4.04×10-6·d-0.141,Kxy=3.048×10-6·d-0.246,Kyy=5.99×10-6·d-0.114；G)与Blum等人的研究结果对比得,本文的计算结果偏小,这是因为本文考虑了节理的粗糙性及非贯通节理的弱化作用的影响。