CN112115599A - 一种密集钻孔弱化顶板孔间距计算方法 - Google Patents
一种密集钻孔弱化顶板孔间距计算方法 Download PDFInfo
- Publication number
- CN112115599A CN112115599A CN202010945806.3A CN202010945806A CN112115599A CN 112115599 A CN112115599 A CN 112115599A CN 202010945806 A CN202010945806 A CN 202010945806A CN 112115599 A CN112115599 A CN 112115599A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- hole
- mass around
- stress
- mpa
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Earth Drilling (AREA)
Abstract
本发明公开了一种密集钻孔弱化顶板孔间距计算方法,首先根据钻孔周围岩体塑性区的边界推导得到单个钻孔周围塑性区发育半径的理论计算公式;在塑性区发育计算公式中考虑到实际工程对塑性区发育的影响,得到钻孔周围岩体原位岩石力学性能;同时对钻孔周围的围岩破碎程度进行分级以确定放大系数,进而计算得到密集钻孔塑性区发育半径;最后可以确定密集钻孔弱化顶板孔间距。合理的密集钻孔孔间距可以使得钻孔的卸压区相互叠加、互相连通,从而形成人工岩体结构弱化带。本发明在保证了弱化顶板效果的同时,可有效的控制密集钻孔的孔间距大小,极大的减轻了施工作业工作量,进一步提高了矿井安全系数及生产效率。
Description
技术领域
本发明属于矿山顶板岩层控制技术领域,尤其涉及一种密集钻孔弱化顶板孔间距计算方法。
背景技术
弱化顶板通过改变顶板岩体的物理力学性质,减小顶板悬露面积,防止或减弱大面积顶板来压,以达到工作面安全、高效生产,是矿山综采的重要手段。
目前弱化顶板的常用方法主要有聚能***和水力压裂。但聚能***难以保证***孔孔间切缝率和***裂缝方向,存在***可控性低、对巷道围岩扰动影响严重等问题,且不适用于高瓦斯矿井。水力压裂虽然对巷道围岩扰动较小,适用性强,但是裂缝角度、裂缝高度、裂缝贯穿率有限,现场操作复杂,学习成本较高。针对以上问题有学者提出密集钻孔弱化顶板切顶卸压技术能够有效弥补上述方法的部分缺陷,能够有效解决煤矿深部开采沿空留巷过程中小煤柱开采时冲击地压过大的问题;也能够针对采空区坚硬顶板无法及时垮落、顶板内弹性势能大量积聚的情况,及时切顶卸压。但目前该方法对密集钻孔施工时孔间距的具体参数并未有科学计算方法。钻孔间距作为该技术的关键参数,若选取的孔间距太大,则不能满足顶板弱化的要求,无法达到卸压的效果;若选取的孔间距太小,虽然可以达到顶板弱化卸压效果,但劳动成本会大大增加。
发明内容
发明目的:针对以上问题,本发明提出一种密集钻孔弱化顶板孔间距计算方法,可针对不同围岩条件进行钻孔施工间距计算,满足实际工程需要。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种密集钻孔弱化顶板孔距计算方法,包括如下步骤:
S1,根据钻孔围岩受力情况,设定整个钻孔径向的受力作用一致,参考图6,以钻孔中心为圆心,建立二维极坐标系;
S2,根据弹性力学,极坐标下钻孔周围岩体的应力分量的公式为:
式中,ρ为钻孔周围岩体任意一点的极坐标值;a为钻孔半径,mm;P为钻孔所在位置的水平地应力,MPa;σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;τθ为钻孔周围岩体任意一点处的环向剪应力,MPa;
S3,根据平衡条件和莫尔强度条件,基于步骤S2所述参数σρ、σθ和τρ,得到钻孔周围岩体中任意一点所在平面状态的极值应力公式为:
式中,σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;
进一步地,列出钻孔周围岩体在该平面上所受剪应力值τn和钻孔周围岩体在该平面上所受正应力σn的常用求解公式为:
式中,τn为钻孔周围岩体在该平面上所受剪应力值,MPa;σn为钻孔周围岩体在该平面上所受正应力,MPa;σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;
同时根据莫尔—库伦屈服条件,即岩体某个平面上的剪应力τn达到极限值时,岩体发生屈服;而该极限值与钻孔周围岩体损伤后的内聚力C、钻孔周围岩体的内摩擦角及钻孔周围岩体在该平面上所受正应力σn有关,公式为:
进一步地,公式(2)、公式(3)和公式(4)联立可得到钻孔周围岩体任意一点的屈服条件式:
式中,σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;C为钻孔周围岩体损伤后的内聚力,MPa;为钻孔周围岩体的内摩擦角,度。
S4,根据钻孔周围岩体的应力分布规律和弹性理论,钻孔周围岩体的应力状态满足屈服条件时得到钻孔塑性区的发育半径,即联立步骤S2所述公式(1)和步骤S3所述公式(5),得到钻孔周围岩体塑性区的发育半径的计算公式:
式中:R为钻孔周围岩体弹塑性区边界线上的任意一点的极坐标值;ρ为钻孔周围岩体任意一点的极坐标值;a为钻孔半径,mm;C为钻孔周围岩体损伤后的内聚力,MPa;为钻孔周围岩体的内摩擦角,度;P为钻孔所在位置的水平地应力,MPa;
S5,基于步骤S4所述钻孔塑性区的发育半径,得到密集钻孔弱化顶板孔间距的计算公式为:
L=2Ru (7)
式中,L为密集钻孔弱化顶板孔间距,mm;R为钻孔周围岩体弹塑性区边界线上的任意一点的极坐标值;u为钻孔交界半径放大系数,可根据地质条件的不同取1.2-1.6;
有益效果:与现有技术相比,本发明的技术方案具有以下有益的技术效果:
本发明提出了一种密集钻孔弱化顶板孔间距的计算方法,基于岩石的弹-脆-塑模型及摩尔库伦准则,计算单孔的塑性区发育情况,并综合考虑多重影响因素,计算得到合适的孔间距,使得密集钻孔岩体的塑性区能相互叠加、互相连通,形成人工岩体结构弱化带,科学有效的解决煤矿深部开采过程中,开采时冲击地压过大的问题;并能够用于放顶煤强制放顶等强矿压问题。
附图说明
图1是密集钻孔弱化顶板切顶卸压方法的钻孔布置平面图;
图2是图1中I-I向剖视图;
图3是钻孔成像仪示意图;
图4是原位岩石钻孔剪切仪示意图;
图5是水压致裂原位地应力测试***示意图;
图6是单个钻孔的围岩的弹塑性区发育图;
图7是密集钻孔弱化顶板构成人工弱化带示意图;
图中,1-上区段巷道;2-上区段工作面;3-采空区;4-下区段工作面;5-巷道切顶卸压密集钻孔;6-直接顶;7-强制放顶密集钻孔。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。其中,所描述的实施例是本发明一部分实施例,而不是全部的实施例。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。
本发明所述的一种密集钻孔弱化顶板孔间距计算方法是密集钻孔顶板弱化的重要参数计算方法。密集钻孔弱化顶板这一技术可应用于坚硬顶板切顶卸压及放顶煤强制放顶等,通过密集钻孔形成人工弱化带,在岩体自重或采动的不平衡高应力影响下,顶板能够沿此弱化带切落或顶煤开始垮落。下面以某矿的地质条件为例对密集钻孔的孔间距进行计算:
S1,回采工作面布置方式参考图1和图2,为了使回采工作面的直接顶6及时垮落,降低严空巷道的维护难度,选择在上区段工作面巷道沿煤柱一侧的顶板布置巷道切顶卸压密集钻孔5,密集钻孔半径为30mm,巷道切顶卸压密集钻孔5沿着巷道的延伸方向,垂直于巷道直接顶6;若要对工作面放顶煤强制放顶可在回采工作面靠近煤壁处布设强制放顶密集钻孔7。
S2,对相关顶板围岩在开采扰动影响下岩石力学参数进行准确测量,在超前开切眼5-10m内布设原位测试孔;参考图4,采用井下原位测试***实测获得巷道开挖后工作面回采影响下,顶板围岩钻孔周围岩体损伤后的内聚力C为1.1MPa、钻孔周围岩体的内摩擦角为30°;参考图5,采用水压致裂原位地应力测试***得到该矿顶板围岩的地应力为14MPa;参考图3,通过钻孔窥视仪对单一钻孔后一定时间内围岩裂隙发育程度进行分析,参考图7,考虑密集钻孔的群体效应,根据围岩裂隙发育情况及密集钻孔相互耦合作用选定放大系数u取1.5。
S3,参考图6,根据钻孔围岩受力情况,设定整个钻孔径向的受力作用一致,以钻孔中心为圆心,建立二维极坐标系;根据弹性力学得到极坐标下钻孔周围岩体中的应力分量公式为:
式中,ρ为钻孔周围岩体任意一点的极坐标值;a为钻孔半径,mm;P为钻孔所在位置的水平地应力,MPa;σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;τθ为钻孔周围岩体任意一点处的环向剪应力,MPa;
S4,根据平衡条件和莫尔强度条件,基于步骤S3所述参数σρ、σθ和τρ,得到钻孔周围岩体中任意一点所在平面状态的极值应力公式为:
式中,σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;
进一步地,列出钻孔周围岩体在该平面上所受剪应力值τn和钻孔周围岩体在该平面上所受正应力σn的常用求解公式为:
式中,τn为钻孔周围岩体在该平面上所受剪应力值,MPa;σn为钻孔周围岩体在该平面上所受正应力,MPa;σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;
进一步地,公式(2)、公式(3)和公式(4)联立可得到钻孔周围岩体任意一点的屈服条件式:
式中,σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;C钻孔周围岩体损伤后的内聚力,MPa;为钻孔周围岩体的内摩擦角,度。
S5,根据钻孔周围岩体的应力分布规律和弹性理论,钻孔周围岩体的应力状态满足屈服条件时得到钻孔塑性区的发育半径,即联立步骤S3所述公式(1)和步骤S4所述公式(5),得到钻孔周围岩体塑性区的发育半径的计算公式:
式中:R为钻孔周围岩体弹塑性区边界线上的任意一点的极坐标值;ρ为钻孔周围岩体任意一点的极坐标值;a为钻孔半径,mm;C钻孔周围岩体损伤后的内聚力,MPa;为钻孔周围岩体的内摩擦角,度;P为钻孔所在位置的水平地应力,MPa;
S6,基于步骤S5所述钻孔塑性区的发育半径,得到密集钻孔弱化顶板孔间距的计算公式为:
L=2Ru (7)
式中,L为密集钻孔弱化顶板孔间距,mm;R为钻孔周围岩体弹塑性区边界线上的任意一点的极坐标值;u为钻孔交界半径放大系数,本式u取值1.5;
基于步骤S5所述参数R和步骤S6所述参数u,计算得到密集钻孔弱化顶板孔间距L=119mm,由于工程施工中存在一定误差,考虑到容差,密集钻孔选择110-130mm为钻孔间距。
Claims (5)
1.一种密集钻孔弱化顶板孔间距计算方法,其特征在于,所述方法具体包括如下步骤:
S1:根据钻孔围岩受力情况,设定整个钻孔径向的受力作用一致,以钻孔中心为圆心,建立二维极坐标系;
S2:根据弹性力学,得到极坐标下钻孔周围岩体的应力分量的公式;
S3:根据平衡条件和莫尔强度条件得到钻孔周围岩体任意一点所在平面状态的极值应力公式,根据钻孔周围岩体在满足莫尔—库伦屈服条件时,结合钻孔周围岩体的原岩应力,计算得到钻孔周围岩体任意一点的屈服条件式;
S4:根据弹性理论,由满足屈服条件时的钻孔周围岩体的应力状态计算得到塑性区的发育半径的计算公式;
S5:基于步骤S4所述塑性区的发育半径,计算得到密集钻孔弱化顶板孔间距的数值。
3.根据权利要求2所述的一种密集钻孔弱化顶板孔间距计算方法,其特征在于:步骤S3的具体操作步骤如下:
S3.1:根据平衡条件和莫尔强度条件,得到钻孔周围岩体中任意一点所在平面状态的极值应力公式为:
式中,σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;σρ为钻孔周围岩体任意一点处的径向应力,MPa;σθ为钻孔周围岩体任意一点处的环向应力,MPa;τρ为钻孔周围岩体任意一点处的径向剪应力,MPa;
S3.2:根据莫尔—库伦屈服条件,即岩体在某个平面上的剪应力达到极限值时,岩体发生屈服;得到该极值关于钻孔周围岩体损伤后的内聚力C、钻孔周围岩体的内摩擦角及钻孔周围岩体在该平面上所受正应力σn的公式为:
S3.3:根据钻孔周围岩体任意一点处平面内的主应力计算得到岩体在该平面上所受剪应力值τn和钻孔周围岩体在该平面上所受正应力σn的公式为:
式中,τn为岩体在该平面上所受剪应力值,MPa;σn为钻孔周围岩体在该平面上所受正应力,MPa;σ1为钻孔周围岩体任意一点处平面内的最大主应力,MPa;σ2为钻孔周围岩体任意一点处平面内的最小主应力,MPa;
S3.4:基于步骤S3.1、S3.2和S3.3所述公式,联立得到钻孔周围岩体任意一点的屈服条件式为:
5.根据权利要求4所述的一种密集钻孔弱化顶板孔间距计算方法,其特征在于:步骤S5中所述密集钻孔弱化顶板孔间距的计算公式为:
L=2Ru
式中,L为密集钻孔弱化顶板孔间距,mm;R为钻孔周围岩体弹塑性区边界线上的任意一点的极坐标值;u为钻孔交界半径放大系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010945806.3A CN112115599B (zh) | 2020-09-10 | 2020-09-10 | 一种密集钻孔弱化顶板孔间距计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010945806.3A CN112115599B (zh) | 2020-09-10 | 2020-09-10 | 一种密集钻孔弱化顶板孔间距计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112115599A true CN112115599A (zh) | 2020-12-22 |
CN112115599B CN112115599B (zh) | 2022-09-30 |
Family
ID=73803140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010945806.3A Active CN112115599B (zh) | 2020-09-10 | 2020-09-10 | 一种密集钻孔弱化顶板孔间距计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112115599B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040040717A1 (en) * | 2002-07-08 | 2004-03-04 | Hill Gilman A. | Method for upward growth of a hydraulic fracture along a well bore sandpacked annulus |
CN204419145U (zh) * | 2014-12-29 | 2015-06-24 | 西安科技大学 | 煤层液态二氧化碳压裂装置 |
CN106528963A (zh) * | 2016-10-21 | 2017-03-22 | 河南理工大学 | 一种卸压钻孔间排距的设计方法 |
CN107067333A (zh) * | 2017-01-16 | 2017-08-18 | 长沙矿山研究院有限责任公司 | 一种高寒高海拔高陡边坡稳定性监控方法 |
CN108920851A (zh) * | 2018-07-10 | 2018-11-30 | 山东科技大学 | 一种基于目标支承压力的卸压钻孔间距确定方法 |
CN108951940A (zh) * | 2018-07-20 | 2018-12-07 | 中国地震局工程力学研究所 | 一种采用应力自释放连接装置的框架填充墙及作法 |
CN110219592A (zh) * | 2019-05-27 | 2019-09-10 | 中国矿业大学(北京) | 一种实现厚煤层充分卸压的立体网络化卸压钻孔布置方法 |
CN111520183A (zh) * | 2020-05-13 | 2020-08-11 | 安徽理工大学 | 厚层砂岩下煤层群开采***切顶卸压增透瓦斯治理方法 |
CN111577232A (zh) * | 2020-05-21 | 2020-08-25 | 重庆市能源投资集团科技有限责任公司 | 一种煤矿井下控制压裂安全保障方法 |
-
2020
- 2020-09-10 CN CN202010945806.3A patent/CN112115599B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040040717A1 (en) * | 2002-07-08 | 2004-03-04 | Hill Gilman A. | Method for upward growth of a hydraulic fracture along a well bore sandpacked annulus |
CN204419145U (zh) * | 2014-12-29 | 2015-06-24 | 西安科技大学 | 煤层液态二氧化碳压裂装置 |
CN106528963A (zh) * | 2016-10-21 | 2017-03-22 | 河南理工大学 | 一种卸压钻孔间排距的设计方法 |
CN107067333A (zh) * | 2017-01-16 | 2017-08-18 | 长沙矿山研究院有限责任公司 | 一种高寒高海拔高陡边坡稳定性监控方法 |
CN108920851A (zh) * | 2018-07-10 | 2018-11-30 | 山东科技大学 | 一种基于目标支承压力的卸压钻孔间距确定方法 |
CN108951940A (zh) * | 2018-07-20 | 2018-12-07 | 中国地震局工程力学研究所 | 一种采用应力自释放连接装置的框架填充墙及作法 |
CN110219592A (zh) * | 2019-05-27 | 2019-09-10 | 中国矿业大学(北京) | 一种实现厚煤层充分卸压的立体网络化卸压钻孔布置方法 |
CN111520183A (zh) * | 2020-05-13 | 2020-08-11 | 安徽理工大学 | 厚层砂岩下煤层群开采***切顶卸压增透瓦斯治理方法 |
CN111577232A (zh) * | 2020-05-21 | 2020-08-25 | 重庆市能源投资集团科技有限责任公司 | 一种煤矿井下控制压裂安全保障方法 |
Non-Patent Citations (3)
Title |
---|
JIANGWEI LIU 等: "Mechanisms of crack initiation and propagation in dense linear multihole directional hydraulic fracturing", 《SHOCK AND VIBRATION》 * |
刘天啸: "高应力巷道钻孔卸压机理及让压支护技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
杨敬轩 等: "围岩孔裂隙充水承压***过程分析", 《中国矿业大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112115599B (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111270987B (zh) | 一种煤矿井下远程区域精准防控冲击地压方法 | |
CN112780340B (zh) | 一种煤矿井上下区域超前防治冲击地压的方法 | |
CN111305876A (zh) | 深部巷道锚固-劈裂注浆-水力压裂卸压协同控制方法 | |
CN110685689B (zh) | 一种采动应力影响下大变形巷道卸压方法 | |
CN110966002B (zh) | 一种基于密集钻孔的切顶卸压方法 | |
CN111255463B (zh) | 坚硬顶板梳状长钻孔分段压裂巷道变形源头治理方法 | |
CN110792419B (zh) | 一种煤矿冲击地压井上下超前预控方法 | |
CN108361011B (zh) | 从井下巷道定向水力压裂控制采动影响范围的方法 | |
CN104913696A (zh) | 边坡开挖光面***法施工方法 | |
CN110206547B (zh) | 适应于深部地层及高应力区域竖井掘进与支护的方法 | |
CN114251103B (zh) | 一种定向切缝压裂顶板大巷防冲护巷方法及安全采矿方法 | |
CN110344831A (zh) | 切顶卸压无煤柱沿空自成巷留巷方法 | |
CN112096383B (zh) | 一种沿空巷道脉冲切顶卸压方法 | |
CN112160792A (zh) | 一种井下坚硬顶板的分段水力压裂工作方法 | |
CN111985101A (zh) | 一种深井冲击危险巷道卸支耦合防冲方法 | |
CN112983418A (zh) | 一种煤矿井下采煤工作面回撤通道水力压裂卸压的方法 | |
CN110454164B (zh) | 冲击矿压巷道缓冲吸能带水力预置方法 | |
CN114673497A (zh) | 一种基于坚硬顶板控向致裂的悬顶处理方法 | |
CN113216912A (zh) | 降低矿震等级用射孔枪及水压致裂降低矿震等级方法 | |
CN112115599B (zh) | 一种密集钻孔弱化顶板孔间距计算方法 | |
CN110043262A (zh) | 一种煤矿坚硬顶板水平井压裂裂缝井上下联合监测方法 | |
CN117167011A (zh) | 破碎围岩基于顶板长水平孔的串-并联全长卸压布孔方法 | |
CN115030719B (zh) | 水力压裂厚硬岩层与煤层卸压相结合的冲击矿压防治方法 | |
CN113847029B (zh) | 一种地下矿山高大点柱的围柱回采方法 | |
CN114526073A (zh) | 一种两侧充分采动采区上山防冲煤柱设计方法及采矿方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |