CN105631102B - 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 - Google Patents
一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 Download PDFInfo
- Publication number
- CN105631102B CN105631102B CN201510982092.2A CN201510982092A CN105631102B CN 105631102 B CN105631102 B CN 105631102B CN 201510982092 A CN201510982092 A CN 201510982092A CN 105631102 B CN105631102 B CN 105631102B
- Authority
- CN
- China
- Prior art keywords
- drilling
- stress
- rock
- tunnel
- release
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明公开了一种深部高应力巷道钻孔卸压参数的数值模拟确定方法,模拟前先取试验巷道围岩岩样在实验室中进行加卸载试验,拟合得到岩样强度参数与损伤变量的衰减关系,并将其嵌入FLAC3D内置应变软化本构模型中,以试验巷道实测矿压显现规律作为已知特征值,反演岩体的数值计算模型参数。建立深部巷道钻孔卸压数值计算模型,模拟分析卸压钻孔方位、卸压时机及钻孔参数(直径、长度、间排距)对巷道围岩稳定性的影响规律,提出上述影响卸压效果的各因素的确定方法,完善钻孔卸压技术体系。本发明在确定卸压钻孔参数时考虑了卸压程度与围岩稳定性的相互作用关系,设计的卸压钻孔可在保证围岩稳定的前提下发挥最大卸压效果,方法简单,具有很强的应用价值。
Description
技术领域
本发明涉及深部巷道围岩应力转移及围岩变形控制数值模拟的分析领域,具体为尤其是一种考虑围岩峰后强度参数衰减影响的深部巷道应力转移及围岩变形控制的数值模拟确定钻孔卸压参数的分析方法。
背景技术
近年来,我国煤炭资源需求量及开采强度的增加,导致浅部煤炭资源日益枯竭,许多矿区已相继转入深部开采阶段。进入深部开采后,煤岩体的组织结构、基本行为特征和工程响应均发生了根本性变化,巷道掘出后围岩易表现出软岩力学特征,碎胀、扩容等大变形现象严重,巷道维护效果极差,严重制约了深部煤炭资源的安全高效生产。
钻孔卸压技术通过人为钻孔的方法在巷道浅部围岩内部形成一个弱化区,将其周边高应力转移至深部稳定围岩中,同时为围岩膨胀变形提供补偿空间,减小巷道变形。作为巷内应力转移技术的一种,钻孔卸压技术具有工艺简单、施工速度快、工程量小等优点,但制约该技术推广应用的主要原因是钻孔卸压参数的确定,目前对于钻孔卸压参数的设计多是参考一些经验,忽略了钻孔对开挖卸荷后煤岩体的再次扰动及弱化效应,导致巷道冒顶、片帮事故时常发生。总结钻孔卸压参数设计主要存在的问题有:
1)采用数值模拟确定钻孔卸压参数时,对深部巷道围岩强度衰减规律认识不足,很少考虑岩石峰后强度参数的衰减过程,或将峰后衰减过程做简单的线性简化,二者均与实际情况偏差较大,缺少室内试验与现场工程岩体参数的反演与校验过程,导致岩体参数的选取具有一定的盲目性,影响了运算结果的输出;
2)设计钻孔卸压参数时,仅考虑了直径、长度等常规参数,忽略了钻孔直径与间排距的相互作用关系,且对于其它一些影响卸压效果的钻孔参数(如卸压时机、方位等)并未全局考虑,导致钻孔卸压参数的确定无法形成完整的技术体系,制约了钻孔卸压技术的推广应用。
因此,提供一种能够能够为克服上述现有技术中的不足而通过数据模拟确定钻孔卸压参数的涉及提供正确而准确的分析方法,从而为钻孔卸压参数的设定提供精准保障,已经是一个值得研究的问题。
发明内容
为了解决上述现有技术的不足,本发明提供了一种深部高应力巷道钻孔卸压参数确定的数值模拟方法。
本发明技术解决方案:
深部高应力巷道钻孔卸压参数的数值模拟确定方法,包括以下步骤:
步骤1)在工程现场调研阶段:对试验巷道进行①围岩变形破坏特征实时监测和②完整岩样取芯两项工作;工作①分别监测巷道不同深度处围岩位移变化规律、应力演化规律、支护体受力变化规律及围岩节理裂隙、破裂区的扩展演化规律等内容,工作②采用取芯钻对试验巷道完整岩样进行取芯;
步骤2)对步骤1)所取岩样进行实验室加卸载试验:将岩样加载至至峰后应变软化段设计目标点后卸载,获取不同损伤程度的初始损伤岩样,测定初始损伤岩样的力学参数,拟合得到初始损伤岩样力学参数与损伤变量的函数关系,建立深部巷道围岩强度衰减模型;
步骤3)采用FISH编程语言对FLAC3D软件进行二次开发:将步骤2)获得的围岩强度衰减模型嵌入软件内置应变软化本构关系中;建立深部巷道开挖数值计算模型,以步骤1)现场实测获得的真实矿压显现规律作为已知特征值,反演得到岩体的数值计算模型参数;
步骤4)建立深部巷道钻孔卸压数值计算模型:采用步骤3)获得的岩体参数对模型赋参,通过改变卸压钻孔开挖方位、开挖时间、开挖尺寸及开挖密度等工程施工参数,分析各工况下深部巷道围岩应力演化规律及变形破坏特征,选择合理的指标评价巷道的卸压程度,提出巷道钻孔卸压参数的确定方法;
步骤5)采用步骤4)确定的钻孔参数:对试验巷道进行工程应用,验证钻孔卸压参数确定方法的合理性。
所述的步骤1)主要监测设备及监测内容如下:
a)采取布置巷道表面位移测站和多基点位移计的方法,监测试验巷道不同深度处围岩位移变化规律;
b)采用安设钻孔测力计的方法,监测试验巷道开挖后围岩应力演化规律;
c)采用安设压力测试装置的方法,监测试验巷道支护体受力变化规律;
d)采用电磁辐射、地质雷达、钻孔窥视仪等现场测试仪器监测顶板和两帮节理、裂隙、破裂区的扩展演化规律。
所述的步骤2)通过实验室加卸载试验获取损伤变量与损伤岩样力学参数函数关系的方法,具体如下:
步骤a)在工程现场调研阶段,对试验巷道完整岩样进行取芯,在实验室内先进行单轴、三轴压缩试验,测定岩样的峰值强度σ1、轴向应变ε1及环向应变ε3,拟合得到σ1与围压σ3的函数关系,根据Mohr-Coulomb强度准则:σ1=Aσ3+B,式中A和B均为强度准则参数,计算式分别为:和计算得到岩样的内聚力c和内摩擦角
步骤b)以步骤a)获得的相应围压下(巷道真实应力环境)岩样峰值强度σ1的70%-80%作为卸围压点,将完整岩样加载至该点时卸除围压σ3,继续加载σ1,设计9个目标卸载点,分别由σ1的90%逐级递减至10%,加载至目标卸载点后进行卸载,测定卸载点岩样的轴向应变ε1和环向应变ε3,获取不同损伤程度的初始损伤岩样;
步骤c)采用单块法测定初始损伤岩样的力学参数,围压σ3等级为5MPa,逐级递增至巷道实际受力状态,测定初始损伤岩样不同σ3等级下的峰值强度σ1,拟合得到σ1与σ3的函数关系,以步骤a)中Mohr-Coulomb强度准则表达式,计算得到初始损伤岩样的内聚力c和内摩擦角
步骤d)选取塑性剪切应变γp作为损伤变量,由下式γp=|ε1-ε3|计算得到初始损伤岩样的塑性剪切应变,结合步骤a)和步骤c)试验结果,拟合得到γp和c、间的函数关系,建立深部巷道围岩强度衰减模型。
所述的步骤3)反演岩体数值计算模型参数的方法,具体如下:
e)步骤d)将围岩强度衰减模型中γp与FLAC3D应变软化本构模型的参数εps按式进行替换推导,式中ψ为岩样剪胀角,得到εps、c和表示的深部巷道围岩强度衰减模型;
步骤f)采用FLAC3D内置FISH编程语言将步骤e)替换推导得到的围岩强度衰减模型编写成运算代码,嵌入软件内置应变软化模型中,数值模拟过程直接调用二次开发后的应变软化模型;
步骤g)建立单位尺寸(1m×1m×1m)的FLAC3D运算模型,按模型按步骤a)获得的岩样参数赋参,模型下边界设置位移约束,水平方向施加载荷模拟围压,围压等级与室内试验时相同,输出模型的全应力-应变曲线,与室内试验结果对比,校验围步骤d)岩强度衰减模型的合理性,校验无误后继续步骤h:;
步骤h)建立深部巷道开挖数值计算模型,模型侧边界及底边界位移固定,按试验巷道真实应力环境施加应力边界条件,按照按照步骤a)获得的岩样参数进行初始赋参,开挖相同断面尺寸的巷道,运算平衡后输出巷道围岩塑性区及位移数据,与现场实测数据进行比较,若不吻合,采用迭代法逐步修正初始赋参,直至二者基本吻合为止,此时输出的围岩力学参数即为校验获取的岩体数值计算模型参数。
所述的步骤4)提出深部巷道钻孔卸压参数的确定方法,具体如下:
采用分类法建立2种不同的深部巷道钻孔卸压数值计算模型:I相同直径不同方位和II相同方位不同直径,模型I用以确定卸压钻孔方位,模型II用以确定开挖时间、开挖尺寸及开挖密度,模型侧边界及底边界位移固定,巷道开挖尺寸与现场一致,采用反演得到岩体参数对模型赋参;
i)对模型I分别施加静水应力场、垂直应力场和水平应力场3种应力边界条件,并对不同应力环境下巷道开挖不同方位的卸压钻孔,即巷帮水平钻孔、肩角倾斜钻孔、顶板垂直钻孔及不同方位钻孔交互开挖,以无钻孔时作为对比模型,每种应力场8种方案,共24种方案,分析不同应力环境下开挖不同方位钻孔对巷道卸压效果的影响规律,判定不同方案下巷道卸压程度,提出卸压钻孔方位的确定方法;
j)步骤i)完成后,按照试验巷道真实应力环境建立模型II,固定卸压钻孔长度、直径和间排距等参数,模拟卸压钻孔开挖时间对巷道卸压效果的影响规律,提出卸压时机的确定方法;
k)步骤j)完成后,采用模型II,固定卸压时机、钻孔直径和间排距等参数,模拟卸压钻孔长度对巷道卸压效果的影响规律,提出卸压钻孔长度的确定方法;
l)步骤k)完成后,采用模型II,固定卸压时机、钻孔长度等参数,模拟卸压钻孔密度(直径、间排距)对巷道卸压效果的影响规律,提出卸压钻孔直径、间排距的确定方法。
作为一种评价方法,步骤4)所述巷道卸压程度的评价具体如下:以开挖卸压钻孔后巷道周边高应力的转移效果和围岩变形控制效果作为评价卸压效果的直接指标,并依据不同卸压钻孔参数对应的卸压程度不同,提出卸压程度3种分类标准——①非充分卸压、②充分卸压和③过度卸压。
作为一种评价方法,巷道周边应力转移效果的评价具体如下:选取围岩应力增高区内的一个重要特征值——应力峰值,表征巷道周边高应力,以巷道无卸压钻孔时围岩应力峰值(原应力峰值)σpeak及峰值位置L(σpeak)作为已知特征值,通过对不同卸压钻孔参数下巷道围岩产生的新应力峰值σ′peak和位置L(σ′ peak)、以及原应力峰值位置L(σpeak)的应力变化规律综合对比分析,综合评判巷道周边高应力的转移效果。
作为一种评价方法,巷道围岩变形控制效果的评价具体如下:以巷道无卸压钻孔时围岩表面位移量作为已知特征值,通过对不同卸压钻孔参数下巷道围岩表面位移数值模拟结果的对比分析,评判巷道围岩变形控制效果;
所述的巷道卸压程度的分类具体如下:
①非充分卸压:卸压钻孔参数不足以转移巷道周边高应力,或在应力转移效果不明显的情况下,同时增加了巷道围岩变形量,称为非充分卸压;
②充分卸压:卸压钻孔参数在有效转移巷道周边高应力的同时,又能对巷道围岩变形起到一定的控制效果,称为充分卸压;
③过度卸压:卸压钻孔参数的改变对于巷道周边应力的转移作用较小,且由于卸压程度过高,导致巷道围岩变形量的急剧增加,围岩已不能保持自稳状态,称为过度卸压。
所述的卸压钻孔参数的确定具体如下:
①卸压钻孔方位:以巷道所处应力环境作为判定依据,合理的卸压钻孔方位应垂直于巷道最大主应力方向布置;巷道处于垂直应力场时,应以开挖巷帮水平钻孔为主,巷道处于水平应力场时,应以开挖顶板垂直钻孔为主,巷道处于静水应力场时,同时开挖巷帮水平钻孔和顶板垂直钻孔。
②卸压时机:卸压钻孔滞后巷道开挖时间越短卸压效果越好,滞后时间越长,钻孔对近乎趋于稳定的浅部围岩结构产生二次扰动,不利于巷道维护,现场工程应用中,卸压钻孔应尽量紧跟巷道迎头施工,尽可能实现巷道掘进与卸压钻孔开挖的平行作业;
③卸压钻孔长度:采用卸压钻孔长度L与无钻孔时围岩应力峰值距巷道表面距离L(σpeak)的比值L/L(σpeak)作为卸压钻孔长度的确定指标,L/L(σpeak)<1时,卸压钻孔无法有效转移巷道周边高应力,且破坏了浅部围岩结构的完整性,不利于巷道维护;L/L(σpeak)≥2后,虽能一定程度转移围岩内部的应力峰值,但对原应力峰值区L(σpeak)的应力作用并不明显,且延长了巷道不稳定时间,增加了钻孔工程量,同样不利于巷道的维护,因此,合理的卸压钻孔长度应在1≤L/L(σpeak)<2范围内;
④卸压钻孔密度(直径、间排距):卸压钻孔直径与间排距对卸压效果的影响具有一定的相关性,钻孔直径是确定间排距的基础,从减少钻孔工程量角度,应尽可能增大钻孔直径,进而增加钻孔间排距;采用巷道原应力峰值位置L(σpeak)相邻两孔间的应力分布状态作为判定指标,当相邻两孔间的应力峰值相互叠加呈“单峰”曲线状态分布,且峰值强度不低于巷道围岩的原岩应力时,巷道处于充分卸压状态。
积极有益效果:本发明方法采用实验室加卸载试验的方法,拟合得到了初始损伤岩样力学参数与损伤变量的函数关系,建立深部巷道围岩强度衰减模型,采用FISH编程语言将获得的围岩强度衰减模型嵌入FLAC3D软件内置应变软化本构关系中,反演得到岩体的数值计算模型参数,为数值计算提供了可靠的试验数据,使输出结果更为准确;同时,在确定卸压钻孔参数时考虑了卸压程度与围岩稳定性的相互作用关系,设计的卸压钻孔可在保证围岩稳定的前提下发挥最大卸压效果,且方法简单,具有很强的创新意义,拥有广泛的工程应用前景。
附图说明
图1为本发明实施例1的钻孔卸压参数数值模拟确定方法的流程图;
图2为岩样单轴、三轴压缩试验拟合得到σ1与围压σ3的函数关系图;
图3a为J1-J3岩样峰值强度与围压的拟合曲线;
图3b为J4-J6岩样峰值强度与围压的拟合曲线;
图3c为J7-J9岩样峰值强度与围压的拟合曲线;
图4a为初始损伤岩样内聚力与塑性剪切应变之间的拟合关系;
图4b为初始损伤岩样内摩擦角与塑性剪切应变之间的拟合关系;
图5为单轴、三轴压缩试验数值模型图;
图6为实验室试验与数值模拟输出全应力-应变曲线对比图;
图7为深部巷道开挖数值计算模型图;
图8a为数值模拟输出巷道围岩塑性区分布图;
图8b为现场实测巷道围岩松动圈范围图;
图9为巷道不同深度处围岩变形的数值模拟与现场实测对比图;
图10a为相同钻孔直径不同钻孔方位的数值计算模型图;
图10b为相同钻孔方位不同钻孔直径数值计算模型图;
图11a为不同卸压钻孔方位下巷帮垂直应力分布曲线图;
图11b为不同卸压钻孔方位下顶板水平应力分布曲线图;
图12a为静水应力场下不同卸压钻孔方位巷道围岩变形量图;
图12b为垂直应力场下不同卸压钻孔方位巷道围岩变形量图;
图12c为水平应力场下不同卸压钻孔方位巷道围岩变形量图;
图13a为不同钻孔卸压时机下巷帮垂直应力分布曲线图;
图13b为不同钻孔卸压时机下巷道围岩变形曲线图;
图14a为不同钻孔长度下巷帮垂直应力分布曲线图;
图14b为不同钻孔长度下巷道围岩变形曲线图;
图15a为不同钻孔排距下巷帮垂直应力分布曲线图;
图15b为不同卸压钻孔排距下巷道围岩变形曲线图;图15c为不同卸压钻孔排距下邻近孔间垂直应力分布曲线图;
图16a为不同钻孔间距下巷帮垂直应力分布曲线图;
图16b为不同卸压钻孔间距下巷道围岩变形曲线图
图16c为不同卸压钻孔间距下邻近孔间垂直应力分布曲线图;
图17a为不同卸压钻孔直径下巷帮垂直应力分布曲线图;
图17b为不同卸压钻孔直径下巷道围岩变形曲线图;
图17c为不同卸压钻孔间距下邻近孔间垂直应力分布曲线图;
图18为采用设计卸压钻孔参数施工后巷道围岩变形监测曲线图。
具体实施方式
下面结合附图及实施例,对本发明做进一步的说明:
实施例1
实施例以徐矿集团张双楼煤矿-1000m西大巷的生产地质条件及巷道支护工程实践为例,该工程巷道平均埋深1030m,巷道掘进断面为直墙半圆拱形,断面尺寸(宽×高)为4.8m×4.4m,现采用本发明所述的技术对巷道钻孔卸压参数进行设计。本实例的钻孔卸压参数的数值模拟确定方法,包括以下几个步骤:
1)根据本实施例的工程方案,对试验巷道完整围岩进行钻孔取样,并在实验室中进行加卸载试验:
分别对岩样进行单轴和三轴压缩试验,测定岩样的峰值强度σ1、轴向应变ε1及环向应变ε3,拟合得到σ1与围压σ3的函数关系,如图2所示。根据Mohr-Coulomb强度准则,计算得到完整岩样的内聚力c和内摩擦角如表1所示。
表1不同阶段岩样力学参数
阶段 | 内聚力/MPa | 内摩擦角/(°) |
峰值阶段 | 13.2597 | 37.324 |
残余阶段 | 2.9669 | 37.652 |
巷道真实应力环境约为25MPa,确定相应围压岩样峰值强度σ1的80%(102.8MPa)作为卸围压点,将完整岩样加载至该点时卸除围压σ3,继续加载σ1至峰后卸载点,卸载点共9个,分别由σ1的90%逐级递减至10%,获取初始损伤岩样,分别以J1-J9编号表示,测定卸载点岩样的轴向应变ε1和环向应变ε3。
采用单块法测定初始损伤岩样的力学参数,围压σ3等级为5MPa,逐级递增至25MPa,测定初始损伤岩样不同σ3等级下的峰值强度σ1,拟合得到σ1与σ3的函数关系,如图3所示,图3a为J1-J3岩样峰值强度与围压的拟合曲线图,图3b为J4-J6岩样峰值强度与围压的拟合曲线图,图3b为J7-J9岩样峰值强度与围压的拟合曲线图。
依据Mohr-Coulomb度准则表达式,计算得到初始损伤岩样的内聚力c和内摩擦角如表2所示。
表2初始损伤岩样力学参数计算结果
选取岩样的塑性剪切应变γp作为损伤变量,由下式γp=|ε1-ε3|计算得到初始损伤岩样的塑性剪切应变,计算结果如表3所示。
表3初始损伤岩样塑性剪切应变计算结果
依据试验结果,拟合得到γp和c、间的函数关系,建立深部巷道围岩强度衰减模型,如图4所示。图4a为初始损伤岩样内聚力与塑性剪切应变之间的拟合关系图,图4b为初始损伤岩样内摩擦角与塑性剪切应变之间的拟合关系图。
2)围岩强度衰减模型嵌入FLAC3D软件,反演岩体数值计算模型参数:
将强度参数衰减关系式中γp与FLAC3D应变软化本构模型的参数εps按进行替换推导,式中ψ为岩样剪胀角,取12°,得到以εps、c和表示的围岩强度参数衰减函数关系,如下式所示:
(R2=0.8295)
(R2=0.8002)
采用FLAC3D内置FISH编程语言将上式编写成运算代码,嵌入软件内置应变软化模型中,数值模拟过程直接调用二次开发后FISH语言。
建立单位尺寸(1m×1m×1m)的FLAC3D运算模型,图5为单轴、三轴压缩试验数值模型图,模型按表1赋参,模型下边界设置位移约束,水平方向施加载荷模拟围压,围压等级取0MPa、5MPa、15MPa和25MPa,输出模型的全应力-应变曲线,与室内试验结果对比,校验围岩强度衰减模型的合理性,图6为实验室试验与数值模拟输出全应力-应变曲线对比图。
建立深部巷道开挖数值计算模型,如图7所示,模型侧边界及底边界位移固定,上表面施加25MPa载荷,按表1对模型进行初始赋参,开挖宽×高为4.8m×4.4m的巷道断面,输出巷道围岩塑性区与位移分布云图,与现场实测数据进行比较,若不吻合,采用迭代法逐步修正初始赋参,直至二者基本吻合为止,此时输出的围岩力学参数即为校验获取的岩体数值计算模型参数,图8a为数值模拟输出巷道围岩塑性区分布图,图8b为现场实测巷道围岩松动圈范围,图9为巷道不同深度处围岩变形的数值模拟与现场实测对比图,表4为数值模拟反演岩体基本力学参数。
表4数值模拟反演岩体的基本力学参数
3)采用分类法建立2种深部巷道钻孔卸压数值计算模型,如图10所示。图10a为相同钻孔直径不同钻孔方位的数值计算模型图(模型I),图10b为相同钻孔方位不同钻孔直径数值计算模型图(模型II),模型I用以确定卸压钻孔方位,模型II用以确定钻孔卸压时机、开挖尺寸及开挖密度,模型侧边界及底边界位移固定,上表面施加25MPa载荷,按表4对模型进行赋参,开挖宽×高为4.8m×4.4m的巷道断面。
3.1)卸压钻孔方位的确定方法
对模型I施加静水应力场,即取侧压系数λ=1,分别开挖不同方位的卸压钻孔,即巷帮水平钻孔、肩角倾斜钻孔、顶板垂直钻孔及不同方位钻孔交互开挖,以无钻孔时作为对比模型,共8种方案。图11a为不同卸压钻孔方位下巷帮垂直应力分布曲线图,图11b为不同卸压钻孔方位下顶板水平应力分布曲线图。
图8a显示,卸压钻孔方位对巷帮垂直应力转移效果影响的主次顺序为:巷帮水平钻孔>肩角倾斜钻孔>顶板垂直钻孔;当开挖水平钻孔时,对转移巷帮垂直应力效果最好,在水平钻孔开挖的基础上,同时开挖肩角和顶板钻孔对于垂直应力转移效果的改善并不明显。
图8b显示,卸压钻孔方位对顶板水平应力转移效果影响的主次顺序为:顶板垂直钻孔>肩角倾斜钻孔>巷帮水平钻孔;开挖顶板钻孔对转移顶板水平应力效果最好,在此基础上,同时施工肩角钻孔和巷帮钻孔对顶板水平应力峰值及位置的影响并不明显。
分别对模型I施加静水应力场、垂直应力场和水平应力场3种应力边界条件,并开挖不同方位的卸压钻孔,每种应力条件下计算8组模型。图12a为静水应力场下不同卸压钻孔方位巷道围岩变形量,图12b为垂直应力场下不同卸压钻孔方位巷道围岩变形量,图12c为水平应力场下不同卸压钻孔方位巷道围岩变形量。
图12显示,巷道处于垂直应力场时,巷帮变形量大于顶底板,对巷帮开挖水平卸压钻孔,即可有效转移巷帮围岩周边的高应力,控制围岩变形,且钻孔工程量最小;与之相反,巷道处于水平应力场时,合理的卸压钻孔方位应以布置顶板垂直钻孔为主;而巷道处于静水应力环境时,围岩的稳定性受垂直应力和水平应力共同作用,同时开挖巷帮及顶板钻孔可有效转移围岩内部高应力,控制围岩变形。因此,确定卸压钻孔方位应以巷道所处应力环境作为判定依据,合理的卸压钻孔方位应垂直于巷道最大主应力方向布置;巷道处于垂直应力场时,应以开挖巷帮水平钻孔为主,巷道处于水平应力场时,应以开挖顶板垂直钻孔为主,巷道处于静水应力场时,同时开挖巷帮水平钻孔和顶板垂直钻孔。试验巷道所处应力环境为垂直应力场,因此,应以布置巷帮水平钻孔。
3.2)卸压时机的确定方法
采用模型II,固定卸压钻孔长度、直径和间排距等参数,模拟卸压钻孔开挖时间对巷道卸压效果的影响规律。图13a为不同钻孔卸压时机下巷帮垂直应力分布曲线,图13b为不同钻孔卸压时机下巷道围岩变形曲线。
图13显示,卸压钻孔滞后巷道开挖时间越短,越能较早参与围岩应力调整过程,巷道控制效果越好;当在巷道围岩应力调整末期开挖钻孔,虽可一定程度上起到转移围岩高应力的作用,但钻孔的开挖同时破坏了近乎趋于稳定的浅部围岩结构,巷道将继续产生变形,反而对其维护不利。合理的钻孔卸压时机应在巷道开挖后至围岩应力调整趋于稳定之前,即在巷道应力调整剧烈期内开挖卸压钻孔较为合理,现场工程应用中,卸压钻孔应尽量紧跟巷道迎头施工,尽可能实现巷道掘进与卸压钻孔开挖的平行作业,在保证巷道一次支护完成的前提下,越早开挖卸压钻孔,对于控制巷道围岩变形越有利。
3.3)卸压钻孔长度的确定方法
采用模型②,固定卸压时机、钻孔直径和排距等参数,模拟卸压钻孔长度对巷道卸压效果的影响规律。图14a为不同钻孔长度下巷帮垂直应力分布曲线,图14b为不同钻孔长度下巷道围岩变形曲线。
图14显示,当L/L(σpeak)<1时,巷道围岩变形量不减反增,主要由于钻孔长度不能有效转移巷道围岩周边的高应力,且钻孔开挖破坏了浅部围岩结构的完整性,导致巷道浅部围岩承载能力的降低,巷道变形量逐渐增加;当1≤L/L(σpeak)<2时,钻孔可有效转移巷道周边高应力,同时可有效控制巷道变形,卸压效果最好;当L/L(σpeak)≥2后,增加钻孔长度仅体现在改变巷道应力峰值位置方面,对于巷道原应力峰值处L(σpeak)的应力并无明显改善,钻孔长度的增加无疑延长了巷道围岩应力的调整时间,反而不利于控制围岩变形,且增加了巷道钻孔工程量。因此,巷道合理的卸压钻孔长度的取值范围应为1≤L/L(σpeak)<2,对于试验巷道,取L/L(σpeak)=1.5,即钻孔长度为10m。
3.4)卸压钻孔直径与间排距的确定方法
采用模型②,固定卸压时机、钻孔长度等参数,模拟卸压钻孔密度(直径、间排距)对巷道卸压效果的影响规律。图15a为不同钻孔排距下巷帮垂直应力分布曲线,图15b为不同卸压钻孔排距下巷道围岩变形曲线,图15c为不同卸压钻孔排距下邻近孔间垂直应力分布曲线,图16a为不同钻孔间距下巷帮垂直应力分布曲线,图16b为不同卸压钻孔间距下巷道围岩变形曲线,图16c为不同卸压钻孔间距下邻近孔间垂直应力分布曲线,图17a为不同卸压钻孔直径下巷帮垂直应力分布曲线,图17b为不同卸压钻孔直径下巷道围岩变形曲线,图17c为不同卸压钻孔间距下邻近孔间垂直应力分布曲线。
图15显示,卸压钻孔排距>1.8m时,巷道处于非充分卸压状态;0.6m≤钻孔排距≤1.8m时,巷道处于充分卸压状态;钻孔排距<0.6m时,巷道处于过度卸压状态。当巷道处于充分卸压状态时,原应力峰值位置L(σpeak)相邻两孔间的应力相互叠加,形成单峰值曲线,且峰值量值大于巷道围岩的原岩应力。
图16显示,卸压钻孔排距>1.8m时,巷道处于非充分卸压状态;1.2m≤钻孔排距≤1.8m时,巷道处于充分卸压状态;钻孔排距<1.2m时,巷道处于过度卸压状态。当巷道处于充分卸压状态时,原应力峰值位置L(σpeak)相邻两孔间的应力同样相互叠加,形成单峰值曲线,且峰值量值大于巷道围岩的原岩应力。
图17显示,巷道处于充分卸压状态时,减小钻孔直径可使卸压效果向非充分卸压转变,同样,钻孔直径的增大,将导致卸压状态向过度卸压转变,钻孔直径是确定卸压钻孔间排距的基础,卸压钻孔直径越大,其单孔作用半径将显著增加,此时,保证巷道处于充分卸压状态所需的间排距可适当加大,反之则减小。因此,确定卸压钻孔间排距之前,应首先确定卸压钻孔的直径。
卸压钻孔的直径与间排距是影响卸压效果的三个关键因素,三者对卸压效果的影响具有一定的相关性,卸压钻孔直径的改变影响着钻孔间排距的确定,钻孔直径是确定间排距的基础,从减少钻孔工程量角度,应尽可能增大钻孔直径,进而增加钻孔间排距。确定卸压钻孔间排距时,采用巷道原应力峰值位置L(σpeak)相邻两孔间的应力分布状态作为判定指标,当相邻两孔间的应力峰值相互叠加呈“单峰”曲线状态分布,且峰值强度不低于巷道围岩的原岩应力时,巷道处于充分卸压状态。
受张双楼煤矿钻孔施工机具功率限制,卸压钻孔最大钻进直径为115mm,为使巷道处于充分卸压状态,应适当减小钻孔的间排距,最终确定卸压钻孔间排距分别为1000mm、800mm。
4)卸压效果现场工程验证
图18为采用设计卸压钻孔参数施工后巷道围岩变形监测曲线,图中监测曲线显示,卸压钻孔施工100d后,巷道两帮及顶底板移近量分别为79mm和111mm,并趋于稳定,变形量仅为未施工卸压钻孔同期变形的30%左右,围岩维护情况较好,表明设计卸压参数合理。
本发明方法采用实验室加卸载试验的方法,拟合得到了初始损伤岩样力学参数与损伤变量的函数关系,建立深部巷道围岩强度衰减模型,采用FISH编程语言将获得的围岩强度衰减模型嵌入FLAC3D软件内置应变软化本构关系中,反演得到岩体的数值计算模型参数,为数值计算提供了可靠的试验数据,使输出结果更为准确;同时,在确定卸压钻孔参数时考虑了卸压程度与围岩稳定性的相互作用关系,设计的卸压钻孔可在保证围岩稳定的前提下发挥最大卸压效果,且方法简单,具有很强的创新意义,拥有广泛的工程应用前景。
以上实施案例仅用于说明本发明的优选实施方式,但本发明并不限于上述实施方式,在所述领域普通技术人员所具备的知识范围内,本发明的精神和原则之内所作的任何修改、等同替代及改进等,均应视为本申请的保护范围。
Claims (1)
1.一种深部高应力巷道钻孔卸压参数的数值模拟确定方法,其特征在于包括以下步骤:
1)在工程现场调研阶段,对试验巷道进行①围岩矿压数据实时监测和②岩样取芯两项工作;工作①分别采取布置巷道表面位移测站,安设多基点位移计、钻孔测力计、压力测试装置矿压设备,监测巷道不同深度处围岩位移变化规律、应力演化规律及支护体受力变化规律内容,采用电磁辐射、地质雷达、钻孔窥视仪现场测试仪器监测顶板和两帮节理、裂隙、破裂区的扩展演化规律,为数值模拟岩体参数的校验提供现场数据;工作②采用取芯钻对巷道完整岩样进行取芯,为实验室加卸载试验提供测试岩样;
2)对步骤1)获得的取芯岩样进行实验室加卸载试验,包括①将岩样加载至峰后应变软化段各目标点后卸载,获取不同损伤程度的初始损伤岩样,并测定初始损伤岩样的力学参数;②选取损伤变量评价岩样的损伤程度,拟合得到初始损伤岩样力学参数与损伤变量的函数关系,建立深部巷道围岩强度衰减模型,具体如下:
a)对试验巷道完整岩样进行取芯,对完整岩样进行单轴压缩试验和常规三轴压缩试验,测定完整岩样的峰值强度σ1、轴向应变ε1及环向应变ε3,拟合得到σ1与围压σ3的函数关系,根据Mohr-Coulomb强度准则:σ1=Aσ3+B,式中A和B均为强度准则参数,计算式分别为:和计算得到岩样的内聚力c和内摩擦角
b)根据试验巷道应力环境计算围岩受力状态,确定室内试验初始应力加载环境,以步骤a)获得的相应围压下岩样峰值强度σ1的70%-80%作为卸围压点,将岩样加载至该点时卸除围压σ3,模拟巷道开挖对围岩应力环境的扰动,继续加载σ1,设计9个目标卸载点,分别由σ1的90%逐级递减至10%,岩样加载至目标卸载点后卸载,测定卸载点岩样的轴向应变ε1和环向应变ε3,获取不同损伤程度的初始损伤岩样;
c)采用单块法测定初始损伤岩样的力学参数,围压σ3等级为5MPa,逐级递增至巷道实际受力状态,测定损伤岩样不同σ3等级下的峰值强度σ1,拟合得到σ1与σ3的函数关系,以步骤a)中Mohr-Coulomb强度准则表达式,计算得到初始损伤岩样的内聚力c和内摩擦角
d)选取塑性剪切应变γp作为损伤变量,由下式γp=|ε1-ε3|计算得到初始损伤岩样的塑性剪切应变,结合步骤a)和步骤c)试验结果,拟合得到γp和c、间的函数关系,建立深部巷道围岩强度衰减模型;
3)采用FISH编程语言对FLAC3D软件进行二次开发,将步骤2)获得的围岩强度衰减模型嵌入FLAC3D软件内置应变软化本构关系中;建立深部巷道开挖数值计算模型,以步骤1)现场实测获得的真实矿压显现规律作为已知特征值,反演得到岩体的数值计算模型参数,所采用的具体方法是:
e)将步骤d)获得的围岩强度衰减模型中γp与FLAC3D应变软化本构模型的参数εps按式进行替换推导,式中ψ为岩样剪胀角,得到以εps、c和表示的深部巷道围岩强度衰减模型;
f)采用FLAC3D内置FISH编程语言将步骤e)替换推导得到的围岩强度衰减模型编写成运算代码,嵌入软件内置应变软化模型中,数值模拟过程直接调用二次开发后的应变软化模型;
g)建立1m×1m×1m的单位尺寸FLAC3D运算模型,模型按步骤a)获得的岩样参数赋参,模型下边界设置位移约束,水平方向施加载荷模拟围压,围压等级与步骤a)相同,输出模型的全应力-应变曲线,与室内试验结果对比,校验步骤d)得到的围岩强度衰减模型的合理性,校验无误后继续步骤h);
h)建立深部巷道开挖数值计算模型,模型侧边界及底边界位移固定,按试验巷道真实应力环境施加应力边界条件,按照步骤a)获得的岩样参数进行初始赋参,开挖相同断面尺寸的巷道,运算平衡后输出巷道围岩塑性区及位移数据,与步骤1)现场实测数据进行比较,若不吻合,采用迭代法逐步修正初始赋参,直至二者吻合为止,此时输出的围岩力学参数即为校验获取的岩体数值计算模型参数;
4)建立深部巷道钻孔卸压数值计算模型,采用步骤3)获得的岩体参数对模型进行赋参,通过改变钻孔卸压开挖方位、开挖时间、开挖尺寸及开挖密度工程因素的施工参数,分析各工况下深部巷道围岩应力演化规律及变形破坏特征,选择指标评价巷道的卸压程度,得到卸压参数对深部巷道围岩稳定性的影响规律,从而提出巷道钻孔卸压参数的确定方法,所述的选择指标评价巷道卸压程度的具体方法是,以开挖钻孔卸压后巷道周边高应力的转移效果和围岩变形控制效果作为评价卸压效果的直接指标,并依据不同钻孔卸压参数对应的卸压程度不同,提出卸压程度3种分类标准——①非充分卸压、②充分卸压和③过度卸压;评价巷道周边高应力转移效果的具体方法是,选取围岩应力增高区内的一个重要特征值——应力峰值,表征巷道周边高应力,以巷道无钻孔时围岩应力峰值,即原应力峰值σpeak及峰值位置L(σpeak)作为已知特征值,通过对不同钻孔卸压参数下巷道围岩产生的新应力峰值σ′peak和位置L(σ′peak)、以及原应力峰值位置σpeak的应力变化规律综合对比分析,综合评判巷道周边高应力的转移效果;评价巷道围岩变形控制效果的具体方法是,以巷道无钻孔时围岩表面位移量作为已知特征值,通过对不同钻孔卸压参数下巷道围岩表面位移的数值模拟对比分析,评判巷道围岩变形控制效果;评价巷道卸压程度的3种分类标准定义方法为:
①非充分卸压:钻孔卸压参数不足以转移巷道周边高应力,或在应力转移效果不明显的情况下,同时增加了巷道围岩变形量,称为非充分卸压;
②充分卸压:钻孔卸压参数在有效转移巷道周边高应力的同时,又能对巷道围岩变形起到控制效果,称为充分卸压;
③过度卸压:钻孔卸压参数的改变对于巷道周边应力的转移作用小,且由于卸压程度高,导致巷道围岩变形量的增加,围岩已不能保持自稳状态,称为过度卸压;建立深部巷道钻孔卸压数值计算模型,模拟不同钻孔方位、开挖时间、开挖尺寸及开挖密度因素对深部巷道围岩应力演化及变形规律影响时,采用分类法建立2种不同的深部巷道钻孔卸压数值计算模型:I相同直径不同方位和II相同方位不同直径,模型I用以确定钻孔卸压方位,模型II用以确定开挖时间、开挖尺寸及开挖密度,模型侧边界及底边界位移固定,巷道开挖尺寸与现场一致,采用步骤h)反演得到岩体参数对模型赋参,模拟具体方法是:
i)对模型I分别施加静水应力场、垂直应力场和水平应力场3种应力边界条件,并对不同应力环境下巷道开挖不同方位的钻孔卸压,即巷帮水平钻孔、肩角倾斜钻孔、顶板垂直钻孔及不同方位钻孔交互开挖,以无钻孔时作为对比模型,每种应力场8种方案,共24种方案,分析不同应力环境下开挖不同方位钻孔对巷道卸压效果的影响规律,提出钻孔卸压方位的确定方法;
j)步骤i)完成后,按照试验巷道真实应力环境建立模型II,固定钻孔卸压长度、直径和间排距参数,模拟钻孔卸压开挖时间对巷道卸压效果的影响规律,提出卸压时机的确定方法;
k)步骤j)完成后,采用模型II,固定卸压时机、钻孔直径和间排距参数,模拟钻孔卸压长度对巷道卸压效果的影响规律,提出钻孔卸压长度的确定方法;
l)步骤k)完成后,采用模型II,固定卸压时机、钻孔长度参数,模拟钻孔卸压密度对巷道卸压效果的影响规律,提出钻孔卸压直径、间排距的确定方法,其中钻孔密度指泄压钻孔直径和间排距;钻孔卸压参数的确定方法具体如下:
①钻孔卸压方位:以巷道所处应力环境作为判定依据,合理的钻孔卸压方位应垂直于巷道最大主应力方向布置;巷道处于垂直应力场时,应以开挖巷帮水平钻孔为主,巷道处于水平应力场时,应以开挖顶板垂直钻孔为主,巷道处于静水应力场时,同时开挖巷帮水平钻孔和顶板垂直钻孔;
②卸压时机:钻孔卸压滞后巷道开挖时间越短卸压效果越好,滞后时间越长,钻孔卸压对近乎趋于稳定的浅部围岩结构产生二次扰动,不利于巷道维护,现场工程应用中,钻孔卸压应紧跟巷道迎头施工,实现巷道掘进与钻孔卸压开挖的平行作业;
③钻孔卸压长度:采用钻孔卸压长度L与无钻孔时围岩应力峰值距巷道表面距离L(σpeak)的比值L/L(σpeak)作为钻孔卸压长度的确定指标,L/L(σpeak)<1时,钻孔卸压无法有效转移巷道周边高应力,且破坏了浅部围岩结构的完整性,不利于巷道维护;L/L(σpeak)≥2后,虽能转移围岩内部的应力峰值,但对原应力峰值区L(σpeak)的应力作用并不明显,且延长了巷道不稳定时间,增加了钻孔工程量,同样不利于巷道的维护,因此,合理的钻孔卸压长度应在1≤L/L(σpeak)≥2范围内;
④钻孔卸压密度,即钻孔卸压直径和间排距:钻孔卸压直径与间排距对卸压效果的影响具有相关性,钻孔直径是确定间排距的基础,从减少钻孔工程量角度,应增大钻孔直径,进而增加钻孔间排距;采用巷道原应力峰值位置L(σpeak)相邻两孔间的应力分布状态作为确定钻孔间排距的判定指标,当相邻孔间应力峰值相互叠加呈“单峰”曲线分布,且峰值强度不低于巷道围岩的原岩应力时,巷道处于充分卸压状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510982092.2A CN105631102B (zh) | 2015-12-24 | 2015-12-24 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510982092.2A CN105631102B (zh) | 2015-12-24 | 2015-12-24 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631102A CN105631102A (zh) | 2016-06-01 |
CN105631102B true CN105631102B (zh) | 2018-12-04 |
Family
ID=56046031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510982092.2A Active CN105631102B (zh) | 2015-12-24 | 2015-12-24 | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631102B (zh) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106547942B (zh) * | 2016-09-26 | 2019-05-14 | 昆明理工大学 | 一种外倾顺层节理岩体被动侧向岩石压力的计算方法 |
CN106528963B (zh) * | 2016-10-21 | 2019-06-21 | 河南理工大学 | 一种卸压钻孔间排距的设计方法 |
CN107491578B (zh) * | 2017-06-27 | 2020-07-24 | 四川建筑职业技术学院 | 软岩大变形失稳条件下松动区及塑性区快速识别的方法 |
CN108060942B (zh) * | 2017-12-11 | 2019-05-24 | 中国矿业大学(北京) | 一种巷道围岩破坏监测的锚杆索支护方法 |
CN108593875B (zh) * | 2018-03-12 | 2021-04-13 | 徐一蜜 | 煤壁片帮实时预测装置及预测方法 |
CN108562895B (zh) * | 2018-03-12 | 2021-03-02 | 周胡琴 | 煤壁片帮实时预测方法 |
CN108843331B (zh) * | 2018-07-06 | 2019-10-11 | 山东科技大学 | 一种等效钻孔面积下小直径钻孔联合卸压布置方法 |
CN109271738B (zh) * | 2018-10-16 | 2023-05-12 | 山东科技大学 | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 |
CN111324939A (zh) * | 2018-12-14 | 2020-06-23 | 江西省煤田地质局二二四地质队 | 一种基于flac3d数值模拟矿产开采的地表下沉量预算方法 |
CN109655400B (zh) * | 2018-12-26 | 2021-02-26 | 天地科技股份有限公司 | 巷道塑性区分析方法 |
CN109812276B (zh) * | 2019-01-22 | 2020-06-02 | 河南理工大学 | 一种基于动载应力场确定邻空巷道合理位置的方法 |
CN110984968A (zh) * | 2019-12-16 | 2020-04-10 | 山东科技大学 | 一种随钻卸压监测方法 |
CN111175121B (zh) * | 2020-01-21 | 2022-04-29 | 山东科技大学 | 巷道围岩钻孔卸压相似模拟试验系统及使用方法 |
CN111365077B (zh) * | 2020-03-19 | 2021-02-12 | 中国矿业大学(北京) | 一种用于深部巷道煤帮持续大变形的卸压与监控方法 |
CN111400926B (zh) * | 2020-03-27 | 2023-10-10 | 安徽理工大学 | 一种基于flac3d数值软件的巷道围岩热物理参数反演方法 |
CN111767593A (zh) * | 2020-05-21 | 2020-10-13 | 中国矿业大学 | 一种基于地应力场类型确定井下巷道合理布置方位的方法 |
CN112611805B (zh) * | 2020-12-11 | 2024-04-05 | 武汉理工大学 | 一种基于衰减系数的评价围岩松动圈范围的方法 |
CN112554953A (zh) * | 2020-12-18 | 2021-03-26 | 山东科技大学 | 一种窗口介入法判断向斜轴部采场冲击危险性的方法 |
CN112483178A (zh) * | 2020-12-18 | 2021-03-12 | 山东科技大学 | 大窗口破坏法判断深部向斜轴部采场冲击危险性的方法 |
CN113158561B (zh) * | 2021-04-14 | 2023-02-03 | 山东大学 | 适用于多种岩体条件的tbm操作参数优化方法及系统 |
CN114396252B (zh) * | 2021-12-31 | 2023-09-22 | 中国矿业大学 | 一种基于钻孔应力监测的大直径钻孔合理间距确定方法 |
US11492891B1 (en) | 2022-02-16 | 2022-11-08 | Liaoning University | Method and apparatus of controlling drilling for rock burst prevention in coal mine roadway |
CN115508548B (zh) * | 2022-11-22 | 2023-04-07 | 山东科技大学 | 应力环境和围岩结构相结合的支护结构的确定方法及系统 |
CN116595809B (zh) * | 2023-07-18 | 2023-10-03 | 中国矿业大学(北京) | 地下工程围岩钻进卸压-探测评价方法 |
CN116933588B (zh) * | 2023-07-18 | 2024-03-29 | 安徽建筑大学 | 不同岩性深部巷道钻孔卸压参数分析方法 |
CN117074171B (zh) * | 2023-07-21 | 2024-03-08 | 中国矿业大学(北京) | 一种深部煤岩应力场耦合诱发煤岩失稳预警方法 |
CN117027928B (zh) * | 2023-09-03 | 2024-03-08 | 山东科技大学 | 基于冲击地压巷道分类的卸压钻孔回填方法 |
CN117172063A (zh) * | 2023-09-06 | 2023-12-05 | 山东能源集团有限公司 | 一种深孔爆破改性邻空巷道悬臂结构效果的评价方法 |
CN117290928B (zh) * | 2023-09-25 | 2024-06-14 | 西南交通大学 | 基于随钻参数的隧道围岩力学参数的反演方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102220866A (zh) * | 2011-04-17 | 2011-10-19 | 山东科技大学 | 深部煤巷冲击地压卸压与加固协同防治方法 |
CN102425432A (zh) * | 2011-12-15 | 2012-04-25 | 湖南科技大学 | 深部高应力碎胀围岩的适时让压支护系统 |
CN104331533A (zh) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | 一种封闭不良钻孔采动响应数值模拟方法 |
-
2015
- 2015-12-24 CN CN201510982092.2A patent/CN105631102B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102220866A (zh) * | 2011-04-17 | 2011-10-19 | 山东科技大学 | 深部煤巷冲击地压卸压与加固协同防治方法 |
CN102425432A (zh) * | 2011-12-15 | 2012-04-25 | 湖南科技大学 | 深部高应力碎胀围岩的适时让压支护系统 |
CN104331533A (zh) * | 2014-09-18 | 2015-02-04 | 山东科技大学 | 一种封闭不良钻孔采动响应数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
深部厚顶煤巷道围岩稳定性相似模型试验研究;肖同强 等;《煤炭学报》;20140703;第39卷(第6期);第1016-1022页 * |
深部巷道围岩钻孔卸压与围岩控制技术研究;郑贺 等;《矿业安全与环保》;20140922;第41卷(第5期);第51-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105631102A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631102B (zh) | 一种深部高应力巷道钻孔卸压参数的数值模拟确定方法 | |
CN110390152B (zh) | 一种模拟巷道围岩裂隙演化的离散元方法 | |
Ma et al. | Experimental Study on Vibration Reduction Technology of Hole‐by‐Hole Presplitting Blasting | |
Yang et al. | Study on presplitting blasting the roof strata of adjacent roadway to control roadway deformation | |
CN114251103B (zh) | 一种定向切缝压裂顶板大巷防冲护巷方法及安全采矿方法 | |
Hu et al. | Study on failure depth of coal seam floor in deep mining | |
Yang et al. | An analysis method for evaluating the safety of pressure water conveyance tunnel in argillaceous sandstone under water-weakening conditions | |
Li et al. | The breaking span of thick and hard roof based on the thick plate theory and strain energy distribution characteristics of coal seam and its application | |
Tao et al. | Analysis of the Critical Safety Thickness for Pretreatment of Mined‐Out Areas Underlying the Final Slopes of Open‐Pit Mines and the Effects of Treatment | |
He et al. | Overburden Damage Degree‐Based Optimization of High‐Intensity Mining Parameters and Engineering Practices in China’s Western Mining Area | |
Sakhno et al. | Numerical studies of floor heave control in deep mining roadways with soft rocks by the rock bolts reinforcement technology | |
Liu et al. | Stability analysis and monitoring method for the key block structure of the basic roof of noncoal pillar mining with automatically formed gob‐side entry | |
Zou et al. | Numerical investigation on influential factors for quality of smooth blasting in rock tunnels | |
Wu et al. | Research on Plastic Zone Evolution Law of Surrounding Rock of Gob‐Side Entry Retaining under Typical Roof Conditions in Deep Mine | |
Guo et al. | A three-dimensional supporting technology, optimization and inspiration from a deep coal mine in China | |
Sakhno et al. | Numerical modeling of controlling a floor heave of coal mine roadways with a method of reinforcing in wet soft rock | |
Li et al. | Study on the influence of in situ stress distribution on the stability of roadway surrounding rock | |
Lou et al. | A control approach of the roof in no‐pillar roadway formed by roof cutting and pressure releasing | |
Cao et al. | Analysis of cumulative damage for shared rock in a neighborhood tunnel under cyclic blasting loading using the ultrasonic test | |
Qi et al. | Stability Control Mechanism of High‐Stress Roadway Surrounding Rock by Roof Fracturing and Rock Mass Filling | |
CN114526073B (zh) | 一种两侧充分采动采区上山防冲煤柱设计方法及采矿方法 | |
Sun | Research on the deformation mechanism of mining roadway stratiform surrounding rock with nonuniform stress field | |
Yang et al. | Information‐Based Construction of High‐Speed Railway Tunnel | |
Wang et al. | An improved numerical simulation approach for the failure of rock bolts subjected to tensile load in deep roadway | |
Yang et al. | Orthogonal analysis and numerical simulation of rock mechanics parameters in stress field of shaft heading face |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |