采矿数值模拟研究中采动岩体力学参数确定研究.doc

/week114 采矿数值模拟研究中采动岩体力学参数确定研究1 陈晓祥 1，谢文兵 2 ，魏文政 3，宋勤发 3 1 中国矿业大学建筑工程学院（221008） 2 中国矿业大学安全与工程学院（221008） 3 郑州矿业集团公司(450000) E-mail：chenxxiang@163.com 摘 要：在实验室和现场观测数据的基础上，结合前人的一些经验公式，建立了一种用以 确定数值模型中煤岩体物理力学参数的较为实用有效的方法，应用此方法可以使得数值模拟 参数合理化，使模拟结果接近真实的反映煤岩层变形和破坏行为. 关键词：采动岩体，数值模拟，力学参数，UDEC2D 1.引 言 众所周知，岩体是由岩块和结构面共同组成的有机复合体，由于岩体本身性质的差异和 受力条件的不同，岩体破坏可以分为脆性破坏和塑性破坏两种形式，但其强度特征严格受到 岩块和结构面强度、岩体结构等因素控制[1]。另外，由于岩体赋存环境的变异性，不能期望 得到岩体参数的精确值，只能通过实验室试验或通过对岩体宏观特性的概括和统计来预测或 估算岩体强度和变形的可能范围。经过试验对比，一般都认为诸如弹性模量、粘聚力和抗拉 强度等煤岩体力学性质的参数取值往往只有煤岩块相应参数值的 1/5~1/3，有的差别可能更 大，比值达到 1/20~1/10，而煤岩体的泊松比一般为煤岩块泊松比的 1.2~1.4 倍[2]。 岩石与岩体力学参数的关系，岩体力学参数与围压的关系，是岩石力学中尚未很好解决 的难题。如何正确选取围岩力学参数一直是数值模拟分析法中具有共性的问题。岩体参数选 取的正确与否，直接关系到最终计算结果的正确性 [3~5] 逼近法、反分析法等一些经验分析方法，这很难使实验室测定的岩石力学与现场的围岩力学 特性联系起来。因此迫切需要研究实验室与现场测定相结合的选取岩体力学参数的科学适用 的方法[6]。 2.岩体地质力学分级方法[7,8] 岩体地质力学分级方法是Bieniawski提出来的，该分级方法认为，影响岩体工程质量的 因素有①岩块的单轴抗压强度、②钻孔岩芯质量RQD、③节理间距、④节理走向及倾角、 ⑤节理面状况以及⑥地下水渗流条件等。由于该方法考虑的因素比较全面，因此，该方法提 出的力学参数估算公式在岩体工程实践中得到较广泛的应用。通过实验室和现场试验，就可 以获得煤岩体RMR岩体分类体系6个方面的相关评判数据，最后得出煤岩体的RMR。得出煤 岩体RMR后，利相关公式就可以求出煤岩体的弹性模量、粘聚力、抗拉强度等参数。岩石 的RMR值可以在表1查出来，也可以参考文献[9]中的详细的RMR实例表。 3.确定采动岩体力学参数[12~14] 3.1 岩体弹性模量的确定 岩体地质力学分级方法其建议的变形模量 Em 计算公式为： Em?? (2RMR?? 100), GPa （1）  1  本课题得到国家自然科学基金重大项目“深部岩体力学基础研究与应用”资助，高等学校博 士学科点专项科研基金资助（20030290017） /week114 当 RMR50 时，则无法计算出 Em，因此，Serafim 和 Pereira 按 RMR 系统提出了另一个 公式，扩大了式（1）的应用范围，使之可以应用与整个 RMR 值范围： RMR??10 E m?? 10 40 , GPa （2） Nicholson 与 Bieniawski 在实验室实验及 RMR 值的基础上又提出了一个由岩块的弹性模 量 Eint 计算岩体的弹性模量 Erm 的公式,两者的比值称为煤岩体弹性模量的减小因子： RF?? Erm Eint ? ? RMR?? )  （3） 而 Mitri 提出的公式（4）也给出了岩块的弹性模量和岩体的弹性模量之间的减小因子。 R F?? E rm E int ? ? RMR 100 ? ) ? ? （4）  表1  ，11]  3.2 确定岩体的体积模量 K 和剪切模量 G 弹性模量 E 和泊松比ν是表征材料力学属性的两个十分重要参数，但在一些情况下，E 和ν并不能十分有效地反映材料的力学行为，如变形等。因此在一些数值模拟系统中（如 FLAC 和 UDEC），一般采用体积模量 K 和剪切模量 G。K 和 G 均由弹性模量 E 和泊松比ν 转化而来，如下面两式所示： K?? G?? E 3(1?? 2? ) E 2(1???? ) (5) (6) 需要注意，当ν接近 0.5 时，K 可能会趋于无穷大，此时不能盲目地采用式（4.3）进行 计算，应根据力学实验或 P 波波速进行估算。  /week114 此外，岩体的剪切模量还可由下式获得： 1 G ? 1 Gr ? 1