分形修正立方定律在采动岩体渗流研究中应用.docxVIP

分形修正立方定律在采动岩体渗流研究中应用 摘要：本文利用岩石破断面粗糙度具有分形性质 的特性，通过对立方定律中粗糙度系数的分形修正，得到更 加精确的渗透系数的求解公式。结合离散元数值模拟对石嘴 山煤矿工作面附近采动岩体的最大渗流速度变化进行分析， 验证了公式的准确性，得到其渗流规律。此方法对研究煤矿 采空区渗流规律具有重要意义，为预防采动岩体裂隙导通而 产生突水事故，确保煤矿安全生产提供理论依据。 关键词：分形修正立方定律粗糙度 0引言 构造裂隙多为陡倾角裂隙，本文将其简化为与层面近似 正交的裂隙，近似为等距分布。裂隙间距与岩层的厚度之间 存在近似线性关系且因地而异。孙广忠对采动岩体等级序的 节理裂隙作了大量统计研究，指出采动岩体的节理裂隙等间 距性的广泛存在并与其所切穿的岩层厚度有关，节理间距大 体为切割岩体厚度的1/2,多数为l/lo吴占寿也对节理裂 隙的等间距分布作了研究，节理裂隙的等间距确实存在，并 指出影响节理裂隙间距大小的因素很多，对于具体地方所得 到的定量关系也不尽一致。综上所述，釆动岩体的渗流可看 作是层面裂隙和构造裂隙所组合的裂隙网络系统为渗透介 质的渗流问题。本文运用离散裂隙网络的方法对釆动岩体的 渗透特性进行了初步探讨。 1分形描述的节理粗糙度公式 裂隙的几何形态和空间展布特征造成了岩体渗透性具 有分布的高度非线性和明显的各向异性。裂隙的粗糙度和裂 隙水力开度同样对渗透性具有重要影响。 首先，假定采动应力场变化只引起裂隙开度发生变化， 而岩石变形导致的裂隙粗糙度变化可以忽略不计。 Louis假设水流服从Darcy定律，根据单相，无紊乱， 粘性不可压缩介质的Navier-Stokes方程，得到了立方定律 的修正式。 根据Louis的立方定律修正式，得到裂隙面相对粗糙度 修正系数C： C=l+8.(1) D■二■ bm——有效隙宽 根据杜时贵等推导的结构面粗糙度系数Barton直边法 简明公式 JRC=400H (2) 把公式(2)带入公式(1)： C=l+8.8HHHH (3) 周创兵对拉西瓦水电站坝址区花岗岩裂隙进行了分形 测量，得到分形维数与粗糙度系数之间的统计关系为： JRC=B (4) 并对各种经验关系式做了进一步的回归分析，得到 JRC=479. 369D-1B (5) 结构面粗糙度两个方向上分形维度Hurst维数可以用下 面的方法描述： H2D ( 0 ) = (alcos 0 +blsin 0 ) + (a2cos2 0 +b2sin2 0 ) (6) 其中，al=B, a2=B, bl=B, b2=B ■ 根据由分形理论可知，一个分形体在E维欧式空间的维 数一般可以表示为： D二E+H (7) 式中，E表示欧式维数。那么对于粗糙表面，E=2,则上 式应写为： D=2+H (8) 把公式(7)和公式(6)带入公式(5)得到： JRC二479. 369[1+ (alcos 0 +blsin 0 ) + (a2cos2 0 + b2sin2 0 ) ] 1.0566 (9) 把(9)带入(3)可以得到最终的结构面粗糙度分形修 正系数C*,把C*带回Louis立方定律公式得到一组平行裂 隙分形描述的立方定律公式： Kf=B (10) B——裂隙内连通面积与总面积之比（连通系数）； S——裂隙平均间距。 岩石的损伤断裂破坏同其内部微裂隙结构特征有密切 的关系，因此结构面粗糙度是岩石的固有属性，只依据岩性 不同而改变，不会因为形状、大小而发生变化。岩石断面具 有分形特征，结构面粗糙度分形修正系数很好地反映了岩石 结构面这一属性。 另外，岩石裂隙具有各向异性，而二维分形Hurst维数 定量的描述了沿两个方向都变化的整个裂隙面粗糙性，因 此，公式（10）对渗透性的计算更加精确。 2采动岩体渗流-应力耦合研究 对于给定的工程岩体，影响岩体渗流的主要是岩体裂隙 的隙宽。当外部应力环境发生变化时，将导致岩体裂隙发生 变化，从而对岩体的渗透特性产生较大的影响。 由于不同年代地质构造，岩体被大量结构面切割，这些 裂隙虽然杂乱无章，但却有一定规律，往往由几组平行的裂 隙所切割。所以可以将岩体裂隙假设为平行、等间距、等隙 宽的几组裂隙进行理论研究。由于岩块的渗透性远远小于裂 隙的渗透性，研究中常常忽略岩块的渗流，认为水仅仅在裂 隙中流动。采用这些假设后，根据公式（10）得到一组平行 裂隙的渗透系数： K0=H （11） 当垂直裂隙方向的应力变化为△ 0时，其产生的位移 TOC \o 1-5 \h \z Au=A e ■ + ■■ (12) 单裂隙水力渗透系数为： k=kHl+A e ■ + ■■■ (13) 由于裂隙宽度b与裂隙间距s相比无穷小，可以忽略， 上式变为［13］： k=kHl+A e ■■■ (14)