粗糙节理网络模拟及裂隙岩体水力耦合特性研究
【摘要】:
众所周知,节理是控制岩体水力学行为的重要因素,研究节理的变形行为及渗透行为是对裂隙岩体工程安全评估的重要组成部分。随着地下开采、隧道、边坡及核废料地下处置库等岩体工程的兴建,对裂隙岩体水力耦合研究的需求越来越迫切。然而目前有关节理变形行为及渗透行为的描述和数值模拟,均局限于光滑、平行板模型,节理产状参数相互独立且节理处于常法向应力边界等假设;在使用裂隙岩体水力学模型时,常常忽略了其前提条件。针对目前研究的不足,结合国际合作项目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等人的研究结果对比得,本文的计算结果偏小,这是因为本文考虑了节理的粗糙性及非贯通节理的弱化作用的影响。
|
|
|
|
1 |
朱珍德,孙钧;裂隙岩体非稳态渗流场与损伤场耦合分析模型[J];水文地质工程地质;1999年02期 |
2 |
于青春,陈德基,薛果夫,大西有三;裂隙岩体一般块体理论初步[J];水文地质工程地质;2005年06期 |
3 |
范波;罗平平;;钻孔压水试验理论研究现状及展望[J];煤炭工程;2010年01期 |
4 |
朱珍德,孙钧;裂隙岩体的渗流场与损伤场耦合分析模型及其工程应用[J];长江科学院院报;1999年05期 |
5 |
杨立中,黄涛;初论环境地质中裂隙岩体渗流-应力-温度耦合作用研究[J];水文地质工程地质;2000年02期 |
6 |
张电吉,汤平,白世伟;节理裂隙岩体渗流与应力耦合分析[J];武汉化工学院学报;2004年04期 |
7 |
陈兴周,李建林,柴军瑞,王如宾;考虑应力状态的裂隙岩体渗透系数确定方法简述[J];地下水;2005年05期 |
8 |
陶振宇;唐方福;;岩体水力模型的比较研究[J];工程勘察;1987年05期 |
9 |
仵彦卿;;利用多孔压水试验资料计算裂隙岩体的渗透系数张量[J];地质灾害与环境保护;1993年01期 |
10 |
贺少辉,廖国华,李中林;裂隙岩体初始裂隙网络渗流模型研究[J];南方冶金学院学报;1995年02期 |
11 |
袁绍国,杨万根,刘占魁;渗透系数的各向异性对裂隙水渗流潜水面的影响[J];水文地质工程地质;1995年03期 |
12 |
陈志坚,孙英学;裂隙岩体力学参数的弱化处理[J];江苏地质;2000年01期 |
13 |
李新平;米健;张成良;王艳丽;;三维应力作用下岩体单个裂隙的渗流特性分析[J];岩土力学;2006年S1期 |
14 |
刘先珊;周创兵;;裂隙岩体非饱和水力应力耦合的不连续介质模型研究[J];岩石力学与工程学报;2007年07期 |
15 |
方涛;柴军瑞;徐文彬;;裂隙岩体渗流水-岩耦合数值分析[J];矿业快报;2007年08期 |
16 |
高俊丽;;降雨入渗对裂隙岩质边坡稳定性研究[J];西部探矿工程;2008年09期 |
17 |
李守巨;龙子祺;刘迎曦;冯文文;;基于分形网络模型的裂隙岩体渗流特性分析(英文)[J];科学技术与工程;2008年21期 |
18 |
张强林;王媛;;裂隙岩体等效热传导系数探讨[J];西安石油大学学报(自然科学版);2009年04期 |
19 |
李睿冉;张国玉;;裂隙岩体离散裂隙网络介质模型研究[J];山东农业大学学报(自然科学版);2010年02期 |
20 |
贺少辉,廖国华,李中林;裂隙岩体裂隙网络渗流场──应力场耦合的刚体元分析[J];南方冶金学院学报;1995年03期 |
|