›› 2014, Vol. 35 ›› Issue (8): 2394-2400.

• 数值分析 • 上一篇    下一篇

岩体裂隙网络等效渗透系数方向性的数值计算

刘日成1, 2,蒋宇静1,李 博1,王肖珊1,徐帮树2   

  1. 1.长崎大学 工学研究科,日本 长崎 8528521;2.山东大学 土建与水利学院 城市地下空间工程系,济南 250061
  • 收稿日期:2014-03-31 出版日期:2014-08-12 发布日期:2014-08-14
  • 通讯作者: 蒋宇静,男,1962年生,教授,博导,主要从事岩石力学方面的教学工作。E-mail:jiang@nagasaki-u.ac.jp E-mail:my1122002006@126.com
  • 作者简介:刘日成,男,1986年生,博士研究生,主要从事岩体裂隙网络渗流特性研究。
  • 基金资助:

    国家自然科学基金面上项目(No. 51279097);日本学术振兴会青年基金(B)(No. 24760690)。

Numerical calculation of directivity of equivalent permeability of fractured rock masses network

LIU Ri-cheng1, 2, JIANG Yu-jing1, LI Bo1, WANG Xiao-shan1, XU Bang-shu2   

  1. 1. Graduate School of Engineering, Nagasaki University, Nagasaki 8528521, Japan; 2. Department of Underground Space, School of Civil and Hydraulic Engineering, Shandong University, Jinan 250061, China
  • Received:2014-03-31 Online:2014-08-12 Published:2014-08-14

摘要: 采用渗流力学理论并结合Monte Carlo方法描述岩体裂隙的随机分布,研究渗流模型的尺寸效应并确定表征单元体积(REV),得到了3种开口度分布形式的等效渗透系数椭圆曲线,建立了等效渗透系数方向性的判别标准。离散裂隙网络(DFN)模型假定流体只在岩体裂隙内部流动,而不通过岩体本身渗流。基于二维离散元程序UDEC并进行二次开发,建立DFN模型,通过改变流体的流动方向,得到不同流动方向下岩体裂隙网络的等效渗透系数,并分析不同的开口度分布形式对岩体裂隙网络等效渗透系数方向性的影响。计算结果表明,表征单元体积存在的条件是等效渗透系数保持稳定且渗透椭圆比较光滑。等效渗透系数的方向性受开口度分布形式的影响很大:当开口度-长度关联分布时,等效渗透系数各向异性;当开口度对数正态分布时,等效渗透系数各向同性;当开口度恒定分布时,等效渗透系数的特性介于二者之间。变化系数(CV)是否大于5%是判定岩体裂隙网络渗透系数是否具有方向性的判别标准。

关键词: 渗透系数方向性, 开口度分布, 离散裂隙网络, 离散元方法, 裂隙岩体

Abstract: Directivity of permeability of fractured rock masses is investigated through numerical models with varying flow directions and 3 kinds of aperture distributions. Besides constant aperture distribution and distribution of apertures correlated with trace length of fractures, a log-normal distribution of apertures, which is verified through experiments and is closer to reality due to its allowance of aperture deviations to some extents, is further introduced into numerical models to evaluate their influences on the directivities of permeability. The basic assumptions are that fluid flow only take place in fractures, and that the intact rock is impermeable and linearly elastic. A large number of stochastic discrete fracture network (DFN) models of varying sizes and varying fracture properties are established to examine the existence of representative elementary volume (REV) and to determine the directivities of equivalent permeability, using a discrete element method (DEM), which provides intact rock masses as an assemblage of discrete blocks and the discontinuities/fractures as interfaces between blocks. These numerical simulation results show that REV exists if the equivalent permeability holds steady and the deviations of coefficient of variance (CV) are less than 10%. The distributions of fracture apertures can influence the directivities of equivalent permeability greatly. The directivity of equivalent permeability is remarkable when apertures are correlated with trace lengths, however, the equivalent permeability seems to be isotropic when apertures are log-normally distributed. The reasons may be that when fracture apertures are correlated with trace lengths, the apertures could vary significantly along different fracture sets with different trace lengths, causing differential permeability values in different directions. The inclined angle of maximum permeability and minimum permeability is about 90° when the curves of directional permeability are smooth. The ‘CV’ value is a crucial factor for identifying the directivities of permeability. When the ‘CV’ value is less than 5%, for example with log-normally distributed fracture apertures, the equivalent permeability in different directions is coincident with the mean permeability, which means the rock masses are hydraulically isotropic and can be treated as a continuum. When the ‘CV’ value varies from 5% to 10%, for example with a constant aperture distribution, the permeability trends to be directional for different rotation models. If the ‘CV’ value exceeds 10%, for example with length correlated apertures, the directivities of permeability become remarkable, resulting in heterogeneous characteristics of rock masses. The coefficient of variance (CV) with the value of varying from 5% to 10%, is the critical condition to identify the directivity of permeability of a fracture network.

中图分类号: 

  • TU 451
[1] 肖瑶, 邓华锋, 李建林, 支永艳, 熊雨. 长期浸泡作用下灌浆加固裂隙岩体劣化效应[J]. 岩土力学, 2019, 40(S1): 143-151.
[2] 支永艳, 邓华锋, 肖瑶, 段玲玲, 蔡佳, 李建林. 微生物灌浆加固裂隙岩体的渗流特性分析[J]. 岩土力学, 2019, 40(S1): 237-244.
[3] 李 崴, 王者超, 毕丽平, 刘 杰, . 辐射流条件下裂隙岩体渗透性表征单元体尺寸 与等效渗透系数[J]. 岩土力学, 2019, 40(2): 720-727.
[4] 肖思友, 苏立君, 姜元俊, 李丞, 刘振宇, . 坡度对碎屑流冲击立式拦挡墙力学特征的影响[J]. 岩土力学, 2019, 40(11): 4341-4351.
[5] 刘艳章, 郭赟林, 黄诗冰, 蔡原田, 李凯兵, 王刘宝, 李 伟, . 冻融作用下裂隙类砂岩断裂特征与强度损失研究[J]. 岩土力学, 2018, 39(S2): 62-71.
[6] 李东奇,李宗利,吕从聪. 考虑裂隙附加水压的岩体断裂强度分析[J]. , 2018, 39(9): 3174-3180.
[7] 杨 昊,张晋勋,单仁亮,武福美,郭志明,. 冻结饱水单裂隙岩体力学特性试验研究[J]. , 2018, 39(4): 1245-1255.
[8] 黄诗冰,刘泉声,程爱平,刘艳章, . 低温裂隙岩体水-热耦合模型研究及数值分析[J]. , 2018, 39(2): 735-744.
[9] 周 喻, 孙 铮, 王 莉, 王 宇, 丁银平, . 单侧限压缩下预制裂隙试样力学特性及板裂化 机制细观研究[J]. 岩土力学, 2018, 39(12): 4385-4394.
[10] 陈庆发,牛文静,郑文师,刘俊广,尹庭昌,范秋雁,. 裂隙岩体块体化程度评价方法的若干问题修正[J]. , 2018, 39(10): 3727-3734.
[11] 崔 臻,盛 谦, . 裂隙岩体等效力学参数结构效应的计算机模拟试验研究[J]. , 2018, 39(10): 3830-3840.
[12] 杨 超,黄 达,蔡 睿,黄润秋,. 张开穿透型单裂隙岩体三轴卸荷蠕变特性试验[J]. , 2018, 39(1): 53-62.
[13] 黄诗冰,刘泉声,程爱平,刘艳章,. 低温岩体裂隙冻胀力与冻胀扩展试验初探[J]. , 2018, 39(1): 78-84.
[14] 郑 卓,刘人太,李术才,王晓晨,张连震,王洪波. 基于Hoek-Brown准则的岩体注浆加固机制研究[J]. , 2017, 38(S1): 409-417.
[15] 符贵军,张思渊,张玉军. 一种双重孔隙-裂隙岩体流变模型及其在地下洞室有限元分析中的应用[J]. , 2017, 38(2): 601-609.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 孙树林,李 方,谌 军. 掺石灰黏土电阻率试验研究[J]. , 2010, 31(1): 51 -55 .
[2] 李英勇,张顶立,张宏博,宋修广. 边坡加固中预应力锚索失效机制与失效效应研究[J]. , 2010, 31(1): 144 -150 .
[3] 李 晶,缪林昌,钟建驰,冯兆祥. EPS颗粒混合轻质土反复荷载下变形和阻尼特性[J]. , 2010, 31(6): 1769 -1775 .
[4] 梁健伟,房营光,谷任国. 极细颗粒黏土渗流的微电场效应分析[J]. , 2010, 31(10): 3043 -3050 .
[5] 王丽艳,姜朋明,刘汉龙. 砂性地基中防波堤地震残余变形机制分析与液化度预测法[J]. , 2010, 31(11): 3556 -3562 .
[6] 李秀珍,王成华,邓宏艳. DDA法和Fisher判别法在潜在滑坡判识中的应用比较[J]. , 2011, 32(1): 186 -192 .
[7] 谷拴成,苏培莉,王建文,王宏科. 烧变岩体特性及其注浆扩散行为研究[J]. , 2009, 30(S2): 60 -63 .
[8] 吉武军. 黄土隧道工程问题调查分析[J]. , 2009, 30(S2): 387 -390 .
[9] 陈力华 ,林 志 ,李星平. 公路隧道中系统锚杆的功效研究[J]. , 2011, 32(6): 1843 -1848 .
[10] 郑 刚 张立明 刁 钰. 开挖条件下坑底工程桩工作性状及沉降计算分析[J]. , 2011, 32(10): 3089 -3096 .