分形修正的立方定律在采动岩体渗流分析中的应用_赵兴强

1、3兰泽全，王宝德，马汉鹏，徐景德.浅谈如何提高煤矿安全培训1兰泽全，王宝德，马汉鹏，徐景德.浅谈如何提高煤矿安全培训的质量J.煤炭工程，2007(05).2陈莉.浅谈如何完善煤矿职工安全教育培训J.科教导刊(中旬参考文献：矿山天地分形修正的立方定律在采动岩体渗流分析中的应用赵兴强 （神 宁煤集 公司石沟 煤 分公司）摘要：本文利用岩石破断面粗糙度具有分形性质的特性，通过对立方定律中粗糙度系数的分形修正，得到更加精确的渗透系数的求解公式。结合离散元数值模拟对石嘴山煤矿工作面附近采动岩体的最大渗流速度变化进行分析，验证了公式的准确性，得到其渗流规律。此方法对研究煤矿采空区渗流规律具有重要意义，为预

2、防采动岩体裂隙导通而产生突水事故，确保煤矿安全生产提供理论依据。关键词：分形修正立方定律粗糙度0 引言定律的修正式。根据 Louis 的立方定律修正式，得到裂隙面相对粗糙度修正系数 C：C=1+8.8 !em 1.5 (1)DHDH = b2mbm有效隙宽构造裂隙多为陡倾角裂隙，本文将其简化为与层面近根据杜时贵等推导的结构面粗糙度系数 Barton 直边似正交的裂隙，近似为等距分布。裂隙间距与岩层的厚度法简明公式之间存在近似线性关系且因地而异 孙广忠对采动岩体等e。JRC=400m(2)级序的节理裂隙作了大量统计研究，指出采动岩体的节理L裂隙等间距性的广泛存在并与其所切穿的岩层厚度有关，把公式

3、(2)带入公式(1)：节理间距大体为切割岩体厚度的 1/2，多数为 1/1。吴占寿1.51.5C=1+8.8 ! 1 !JRCbL (3)也对节理裂隙的等间距分布作了研究，节理裂隙的等间距200m确实存在，并指出影响节理裂隙间距大小的因素很多，对周创兵对拉西瓦水电站坝址区花岗岩裂隙进行了分于具体地方所得到的定量关系也不尽一致。综上所述，采形测量，得到分形维数与粗糙度系数之间的统计关系为：动岩体的渗流可看作是层面裂隙和构造裂隙所组合的裂JRC=D- 0.97148(4)隙网络系统为渗透介质的渗流问题。本文运用离散裂隙网- 35.80710络的方法对采动岩体的渗透特性进行了初步探讨。并对各种经验关

4、系式做了进一步的回归分析，得到1 分形描述的节理粗糙度公式!1.0566(5)裂隙的几何形态和空间展布特征造成了岩体渗透性JRC=479.369 D- 1结构面粗糙度两个方向上分形维度 Hurst 维数可以具有分布的高度非线性和明显的各向异性。裂隙的粗糙度用下面的方法描述：和裂隙水力开度同样对渗透性具有重要影响。（）（22 ）(6)首先，假定采动应力场变化只引起裂隙开度发生变H2D()= a1cos+b1sin + a2cos +b2sin 其中，a1=H0 - H180，a2=H0 +H180，b1=H90 - H270，b2=化，而岩石变形导致的裂隙粗糙度变化可以忽略不计。222Louis

5、 假设水流服从 Darcy 定律，根据单相，无紊乱， H90 +H270粘性不可压缩介质的 NavierStokes 方程，得到了立方2（上接第 117 页）目标的实现有机结合、教师个人成就与集体荣誉有机结标激励、情感激励、物质激励等激励机制，激发教师的积极合，建成一支富有凝聚力、战斗力的教师团队。性和责任感。3.4 建立考核机制，培养创优意识 培训机构在师资3.7 加强培养培训，不断提升教师素质 大力开展教队伍建设中，建立起了教师业务技能水平和常规工作两大 师自身培训，加强对青年教师的培养，促进青年教师快速考核评估体系。对教师出勤情况班级信息员作专门的记录 成长。精心安排教师培训计划和培训内

6、容，注重培训实效。考核。在校内实行考核考评结果公示制，实行奖惩分离，奖 通过加强培训，努力提高教师的师德师风、德育育人能力、惩公开。新制度激发争先创优意识，通过考核评估机制的 理论讲授能力和实践指导能力。形成和实施，在教师队伍中形成了“班长抓班子、班子带队总之，师资队伍建设是一项复杂的系统工程。培训中伍、队伍促发展”、人人“创新争优”的良好局面。心既要建设一支学历高、能力强、结构合理的文化理论教3.5 创新机制，确立和落实科学人才观，推进人才机 师队伍，更要配备一支既有理论基础又有实践经验、掌握制的创新 科学定编，因需设岗，因事择人，合理调配教师， 高新技术、具备较高操作技能的专业课教师队伍，

7、只有这使专业课教师和文化课教师的比例尽快达到国家规定的 样，才能适应培养高技能人才的需要。要求和职教改革实际的需要。不合格师资限期通过进修、培训，达到合格师资水平。建立岗位责任制，确定教师的岗位工作目标、任务和职责；竞争上岗，建立全员聘任制。3.6 改革分配制度，建立激励机制 改革现行的工资刊)，2011(07).分配制度，建立有效的分配激励机制。注重教师的工作成绩和实际贡献，充分调动教师的工作积极性。同时建立目的质量J.煤炭工程，2007(05).118矿山天地图 1结构面粗糙度二维 Hurst 维数图根据由分形理论可知，一个分形体在 E 维欧式空间的维数一般可以表示为：D=E+H(7)式中

8、，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 立方定律公式得到一组平行裂隙分形描述的立方定律公式：3Kf=gb(10)*12SC裂隙内连通面积与总面积之比(连通系数)；S裂隙平均间距。岩石的损伤断裂破坏同其内部微裂隙结构特征有密切的关系，因此结构面粗糙度是岩石的固有属性，只依据岩性不同而改变，不会因为形状、大小而发生变化。岩

9、石断面具有分形特征，结构面粗糙度分形修正系数很好地反映了岩石结构面这一属性。另外，岩石裂隙具有各向异性，而二维分形 Hurst 维数定量的描述了沿两个方向都变化的整个裂隙面粗糙性，因此，公式(10)对渗透性的计算更加精确。2 采动岩体渗流应力耦合研究对于给定的工程岩体，影响岩体渗流的主要是岩体裂隙的隙宽。当外部应力环境发生变化时，将导致岩体裂隙发生变化，从而对岩体的渗透特性产生较大的影响。由于不同年代地质构造，岩体被大量结构面切割，这些裂隙虽然杂乱无章，但却有一定规律，往往由几组平行的裂隙所切割。所以可以将岩体裂隙假设为平行、等间距、等隙宽的几组裂隙进行理论研究。由于岩块的渗透性远远小于裂隙的

10、渗透性，研究中常常忽略岩块的渗流，认为水仅仅在裂隙中流动。采用这些假设后，根据公式(10)得到一组平行裂隙的渗透系数：3K0= gb(11)12sc*当垂直裂隙方向的应力变化为 时，其产生的位移为：- 1u=!KEn + 1s (12)单裂隙水力渗透系数为：3k=k0 #1+ !knEb0 + bs0 -1 $ (13)由于裂隙宽度 b 与裂隙间距 s 相比无穷小，可以忽略，上式变为13：3k=k0 #1+ !knEb0 -1 $ (14)经化简，上式变为：3k=k #1+ b $ (15)0b02.1 自重力引起的裂隙变形 自重应力引起的在裂隙尖端 x=c 处产生的张拉位移为：42 221-

11、 v姨a - c1 =(16)E自重应力引起的裂隙变形，考虑单一裂隙在采动作用下的位移特征。此问题可以应用 S.L.Crouch 提出的位移不连续法。节理裂隙两侧相对位移可视为不连续位移，裂隙位移可表示为：u1= n(17)Kn2.1.1 裂隙的隙宽变化 综合考虑前面讨论的两种情况，可得由开采应力引起的裂隙在垂直裂隙面方向的位移为：42 n221- v姨a - cb1 =1 +u1 =+(18)EKn2.1.2 裂隙渗透系数变化 将式(18)带回式(15)，可得采动应力引起的渗透系数变化值。342n22K =K1+1- v姨a - c+(19)EbbK 0!n当采动应力造成的裂隙尖端破坏程度较

12、小时，ac，上述公式可变为：K=K01+ n (20)bKn具体到采动岩体，原岩应力 n=h，假定裂隙面与水平面夹角为 时，则不同埋深的岩层渗透系数为：h0!322K =K1-h cos +sin (21)bKn由式(21)可以看出，随着埋深的增加，渗透系数逐渐减小，这与现实情况是相符的，由于岩体在地表的渗透系数很容易获得，所以通过上述公式可以求得地下任意深度岩层的渗透系数值。2.2 水压引起的裂隙变形 在采动过程中，饱和岩体裂隙中存在的水压力将会对裂隙面上的每一点形成一个均匀的压力荷载，在这种水压力作用下，裂隙将会有一种119矿山天地张开的趋势，同样利用 Paris 位移公式求解：4!2 2

13、21- v姨a - cw =(22)E可以看出，水压力越大，裂隙尖端的张开位移越大。在裂隙内部法向作用有水压 p 时，水压力引起的裂隙张开位移为：uw= p(23)Kn则综合考虑前面讨论的两种情况，得到由水压引起的裂隙张开位移为：模拟。根据数值模拟结果，采动渗流岩体基本可以分为四个区域：贮水区、渗流屏蔽区、定向渗流区和完全渗流区。渗流屏蔽区由渗流关键层组成，随开采渗流关键层的破断，渗流屏蔽区可能转化为定向渗流区。而完全渗流区后方随开采裂隙的逐渐愈合会转化为渗流屏蔽区。根据分形修正的立方定律结合离散元数值模拟，得到工作面附近采动岩体的最大渗流速度变化如图 2 所示。4p!2 p22bw =w +

14、uw =1- v姨a - c +(24)EKn将式(24)带回式(15)，可得采动应力引起的渗透系数变化值。34!2 p姨22K =K1+1- vp姨a - c+(25)EbbKw0姨n当采动应力造成的裂隙尖端破坏程度较小时，ac，上述公式可变为：3Kw =K0 姨1+ bKpn 姨 (26)根据式(26)可以看出，随着水压的增加，渗透系数增加的幅度很大，室内试验表明，随着水压的增加，裂隙岩样的渗透性明显地增加。2.3 裂隙岩体渗流特征与应力的关系 综合考虑开采应力和水压的联合作用，裂隙岩体的渗透系数为：!2!p+n姨2 2041- vp+姨K=K +K =K1+Eba- c +bK(27)w

15、姨n具体到采动岩体，在自重力影响下，不考虑裂隙端部破坏引起的裂隙位移，则有：h0姨3!22姨K =K1+p- h cos +sin (28)bKn3 算例以石嘴山煤矿为例，用分形修正的立方定律计算采动覆岩的渗透系数，分析采动岩体中的渗流行为。根据对石嘴山煤矿覆岩不同性质岩样的室内渗流试验结果，并应用分形修正的立方定律计算得到石嘴山煤矿采动覆岩的渗透系数见表 1。表 1石嘴山煤矿覆岩渗透系数计算参数岩性渗透系数 K0分形维数粗糙度系数分形渗透系数(10- 6m3/s)H2D()JRCK0(10- 6m3/s)细砂岩0.210.76580.207中砂岩0.210.73040.243粗砂岩0.220

16、.59020.251泥岩0.200.965160.152砂质泥岩0.210.865100.184裂隙初始开度 b=0.01mm，裂隙水密度取 w=1.0 103kg/m3，水的动力黏滞系数 =1.0103。根据石嘴山煤矿柱状图，计算不同深度岩层的渗透性，并用离散元数值图 2覆岩最大渗流速度图4 结论本文利用岩石破断面粗糙度具有分形性质的特性，通过对立方定律中粗糙度系数的分形修正，得到更加精确的渗透系数的求解公式。4.1 分形修正的渗透系数在岩石裂隙面的两个方向上描述其渗透性，相对于以往理论解的渗透系数更加精确。4.2 通过应力与渗透性的耦合，得到不同深度条件下的岩层渗透系数，结合离散元数值模拟，得到采动岩体的渗流规律：随着工作面远离切眼，开切眼附近围岩的裂隙导水能力逐渐减小，最终达到稳定，但渗透能力只是初始渗透能力的几分之一。也就是说，在煤层开采工程中，只要开切眼在初始没有突水，那么在以后的开采过程中，在切眼一侧发生顶板突水的可能性更小。随着工作面推进，工作面附近的采动覆岩的导水能力逐渐逐级递增，直到主关键层破断，发生顶板突水的可能性最大。工作面继续推进，工作面附近采动覆岩的导水系数又开始逐渐降低，最终达到稳定。当导水裂隙发育达