分形修正立方定律在采动岩体渗流分析中应用.docVIP

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