一种深部高应力巷道钻孔卸压参数的数值模拟确定方法
【技术领域】
[0001] 本发明设及深部巷道围岩应力转移及围岩变形控制数值模拟的分析领域，具体为 尤其是一种考虑围岩峰后强度参数衰减影响的深部巷道应力转移及围岩变形控制的数值 模拟确定钻孔卸压参数的分析方法。
【背景技术】
[0002] 近年来，我国煤炭资源需求量及开采强度的增加，导致浅部煤炭资源日益枯竭，许 多矿区已相继转入深部开采阶段。进入深部开采后，煤岩体的组织结构、基本行为特征和工 程响应均发生了根本性变化，巷道掘出后围岩易表现出软岩力学特征，碎胀、扩容等大变形 现象严重，巷道维护效果极差，严重制约了深部煤炭资源的安全高效生产。
[0003] 钻孔卸压技术通过人为钻孔的方法在巷道浅部围岩内部形成一个弱化区，将其周 边高应力转移至深部稳定围岩中，同时为围岩膨胀变形提供补偿空间，减小巷道变形。作为 巷内应力转移技术的一种，钻孔卸压技术具有工艺简单、施工速度快、工程量小等优点，但 制约该技术推广应用的主要原因是钻孔卸压参数的确定，目前对于钻孔卸压参数的设计多 是参考一些经验，忽略了钻孔对开挖卸荷后煤岩体的再次扰动及弱化效应，导致巷道冒顶、 片帮事故时常发生。总结钻孔卸压参数设计主要存在的问题有：
[0004] 1)采用数值模拟确定钻孔卸压参数时，对深部巷道围岩强度衰减规律认识不足， 很少考虑岩石峰后强度参数的衰减过程，或将峰后衰减过程做简单的线性简化，二者均与 实际情况偏差较大，缺少室内试验与现场工程岩体参数的反演与校验过程，导致岩体参数 的选取具有一定的盲目性，影响了运算结果的输出；
[0005] 2)设计钻孔卸压参数时，仅考虑了直径、长度等常规参数，忽略了钻孔直径与间排 距的相互作用关系，且对于其它一些影响卸压效果的钻孔参数(如卸压时机、方位等)并未 全局考虑，导致钻孔卸压参数的确定无法形成完整的技术体系，制约了钻孔卸压技术的推 广应用。
[0006] 因此，提供一种能够能够为克服上述现有技术中的不足而通过数据模拟确定钻孔 卸压参数的设及提供正确而准确的分析方法，从而为钻孔卸压参数的设定提供精准保障， 已经是一个值得研究的问题。
【发明内容】
[0007] 为了解决上述现有技术的不足，本发明提供了一种深部高应力巷道钻孔卸压参数 确定的数值模拟方法。
[000引本发明技术解决方案：
[0009] 深部高应力巷道钻孔卸压参数的数值模拟确定方法，包括W下步骤：
[0010] 步骤1)在工程现场调研阶段:对试验巷道进行①围岩变形破坏特征实时监测和② 完整岩样取忍两项工作;工作①分别监测巷道不同深度处围岩位移变化规律、应力演化规 律、支护体受力变化规律及围岩节理裂隙、破裂区的扩展演化规律等内容，工作②采用取忍 钻对试验巷道完整岩样进行取忍；
[0011]步骤2)对步骤1)所取岩样进行实验室加卸载试验:将岩样加载至至峰后应变软化 段设计目标点后卸载，获取不同损伤程度的初始损伤岩样，测定初始损伤岩样的力学参数， 拟合得到初始损伤岩样力学参数与损伤变量的函数关系，建立深部巷道围岩强度衰减模 型；
[0012]步骤3)采用FI細编程语言对FLAC3d软件进行二次开发:将步骤2)获得的围岩强度 衰减模型嵌入软件内置应变软化本构关系中；建立深部巷道开挖数值计算模型，W步骤1) 现场实测获得的真实矿压显现规律作为已知特征值，反演得到岩体的数值计算模型参数；
[0013] 步骤4)建立深部巷道钻孔卸压数值计算模型:采用步骤3)获得的岩体参数对模型 赋参，通过改变卸压钻孔开挖方位、开挖时间、开挖尺寸及开挖密度等工程施工参数，分析 各工况下深部巷道围岩应力演化规律及变形破坏特征，选择合理的指标评价巷道的卸压程 度，提出巷道钻孔卸压参数的确定方法；
[0014] 步骤5)采用步骤4)确定的钻孔参数:对试验巷道进行工程应用，验证钻孔卸压参 数确定方法的合理性。
[0015] 所述的步骤1)主要监测设备及监测内容如下：
[0016] a)采取布置巷道表面位移测站和多基点位移计的方法，监测试验巷道不同深度处 围岩位移变化规律；
[0017] b)采用安设钻孔测力计的方法，监测试验巷道开挖后围岩应力演化规律；
[0018] C)采用安设压力测试装置的方法，监测试验巷道支护体受力变化规律；
[0019] d)采用电磁福射、地质雷达、钻孔窥视仪等现场测试仪器监测顶板和两帮节理、裂 隙、破裂区的扩展演化规律。
[0020] 所述的步骤2)通过实验室加卸载试验获取损伤变量与损伤岩样力学参数函数关 系的方法，具体如下：
[0021] 步骤a)在工程现场调研阶段，对试验巷道完整岩样进行取忍，在实验室内先进行 单轴、S轴压缩试验，测定岩样的峰值强度01、轴向应变ει及环向应变63,拟合得到01与围压 化的函数关系，根据Mohr-Coulomb强度准则:〇ι = Α〇3+Β，式中A和Β均为强度准则参数，计算 式分别为
计算得到岩样的内聚力C和内摩擦角P ;
[0022] 步骤b)W步骤a)获得的相应围压下(巷道真实应力环境)岩样峰值强度σι的70%- 80 %作为卸围压点，将完整岩样加载至该点时卸除围压〇3,继续加载〇1，设计9个目标卸载 点，分别由σ 1的9 0 %逐级递减至10 %，加载至目标卸载点后进行卸载，现憶卸载点岩样的轴 向应变ε?和环向应变63,获取不同损伤程度的初始损伤岩样；
[0023] 步骤C)采用单块法测定初始损伤岩样的力学参数，围压03等级为5MPa，逐级递增 至巷道实际受力状态，测定初始损伤岩样不同03等级下的峰值强度01，拟合得到01与03的函 数关系，W步骤a)中Mohr-Coulomb强度准则表达式，计算得到初始损伤岩样的内聚力C和内 摩擦角φ ;
[0024] 步骤d)选取塑性剪切应变丫 Ρ作为损伤变量，由下式丫 ρ= I ε广ε31计算得到初始损 伤岩样的塑性剪切应变，结合步骤a)和步骤C)试验结果，拟合得到丫 Ρ和C、則'司的函数关 系，建立深部巷道围岩强度衰减模型。
[0025] 所述的步骤3)反演岩体数值计算模型参数的方法，具体如下：
[0026] e)步骤d)将围岩强度衰减模型中丫 P与化AC3d应变软化本构模型的参数Eps按式
，Φ为岩样剪胀角，得到ePs、c和Ψ 表示的深部巷道围岩强度衰减模型；
[0027] 步骤f)采用FLAC3d内置FI甜编程语言将步骤e)替换推导得到的围度衰减模型编 写成运算代码，嵌入软件内置应变软化模型中，数值模拟过程直接调用二次开发后的应变 软化模型；
[00%]步骤g)建立单位尺寸（ImXlmXlm)的FLAC3d运算模型，按模型按步骤a)获得的岩 样参数赋参，模型下边界设置位移约束，水平方向施加载荷模拟围压，围压等级与室内试验 时相同，输出模型的全应力-应变曲线，与室内试验结果对比，校验围步骤d)岩强度衰减模 型的合理性，校验无误后继续步骤h:;
[0029] 步骤h)建立深部巷道开挖数值计算模型，模型侧边界及底边界位移固定，按试验 巷道真实应力环境施加应力边界条件，按照按照步骤a)获得的岩样参数进行初始赋参，开 挖相同断面尺寸的巷道，运算平衡后输出巷道围岩塑性区及位移数据，与现场实测数据进 行比较，若不吻合，采用迭代法逐步修正初始赋参，直至二者基本吻合为止，此时输出的围 岩力学参数即为校验获取的岩体数值计算模型参数。
[0030] 所述的步骤4)提出深部巷道钻孔卸压参数的确定方法，具体如下：
[0031] 采用分类法建立巧巾不同的深部巷道钻孔卸压数值计算模型：1相同直径不同方位 和II相同方位不同直径，模型I用W确定卸压钻孔方位，模型II用W确定开挖时间、开挖尺 寸及开挖密度，模型侧边界及底边界位移固定，巷道开挖尺寸与现场一致，采用反演得到岩 体参数对模型赋参；
[0032] i)对模型I分别施加静水应力场、垂直应力场和水平应力场3种应力边界条件，并 对不同应力环境下巷道开挖不同方位的卸压钻孔，即巷