CN104268364A - 一种确定煤壁滑移片帮深度及其滑移危险性的方法 - Google Patents

一种确定煤壁滑移片帮深度及其滑移危险性的方法 Download PDF

Info

Publication number
CN104268364A
CN104268364A CN201410354960.8A CN201410354960A CN104268364A CN 104268364 A CN104268364 A CN 104268364A CN 201410354960 A CN201410354960 A CN 201410354960A CN 104268364 A CN104268364 A CN 104268364A
Authority
CN
China
Prior art keywords
depth
coal
slippage
slip
rib
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
Application number
CN201410354960.8A
Other languages
English (en)
Other versions
CN104268364B (zh
Inventor
康天合
李晓坡
杨永康
李晨阳
李恒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN201410354960.8A priority Critical patent/CN104268364B/zh
Publication of CN104268364A publication Critical patent/CN104268364A/zh
Application granted granted Critical
Publication of CN104268364B publication Critical patent/CN104268364B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Solid Fuels And Fuel-Associated Substances (AREA)

Abstract

本发明公开了一种确定煤壁滑移片帮深度及其滑移危险性的方法,其所述方法是假定n个距离顶板高度yi(i=1,2,3,…,n)的圆心和n个滑移圆弧半径Rj(j=1,2,3,…,n),由yi和Rj定义n×n个滑移面Sij,并求其危险系数Kij;作n×n个滑移面Sij的滑移危险系数K曲面,用K=1平面与K曲面相交形成Kcr=1曲线;定义滑移面与顶板交点到煤壁距离为滑移面深度d,作n×n个滑移面Sij的滑移深度d曲面,将Kcr=1曲线在d曲面上投影形成临界滑移深度dcr曲线;在dcr曲线上找出最大值dcrmax(ycr,Rcr),由ycr和Rcr所确定的滑移面为最大深度临界滑移面,其dcrmax值为煤壁滑移片帮深度,在大采高工作面煤壁管理中将以此深度设计防治煤壁片帮的控制措施及其参数。

Description

一种确定煤壁滑移片帮深度及其滑移危险性的方法
技术领域
本发明属于煤炭开采回采工作面围岩控制技术领域,尤其涉及一种确定回采工作面煤壁滑移片帮深度及其滑移危险性的方法。
背景技术
厚及特厚煤层在我国煤炭资源中约占40%,目前大采高综采一次开采高度为3.5~7.0m,大采高综放一次开采厚度为5~20m,其中底层割煤高度3~5m。然而,实践证明大采高综采/放工作面由于煤壁高度的加大,矿压显现更加剧烈,特别是极易发生煤壁片帮。严重的煤壁片帮不仅影响工作面的正常生产,而且还威胁着矿工的生命安全。因此,控制煤壁片帮成为大采高综采/放开采的关键技术问题。
在煤壁片帮方面,现有的研究成果多集中在影响煤壁片帮的因素分析、片帮机理研究、片帮方式、片帮形式和类型以及防治煤壁片帮的技术措施上,尚无查到一种系统有效的预测煤壁片帮深度及其对滑移危险性评价的方法,这一现状直接限制了控制煤壁片帮方法及其参数的准确确定。文献《大采高综采工作面煤壁片帮机理及控制技术》(煤炭科学技术,2008,36(9):1-3,24)采用库仑-摩尔强度理论,得出的压剪式煤壁片帮最大深度的计算式为Δa=M tan(45°-φ/2)和滑落式煤壁片帮最大深度的计算式为Δb=M tan(90°-φ),可见该计算式中仅考虑了煤壁高度M和煤的内摩擦角φ,所考虑因素单一,并且其预测值与实际片帮深度相差较大。
发明内容
本发明的目的是针对现有技术存在的缺陷,提供一种能够准确确定煤壁滑移片帮深度及其滑移危险性的方法。所述方法的内容与步骤为:
1)作图给出煤层顶板CE和底板AF,连接AC形成高度为M的煤壁,在顶板CE之上距CE距离为yi(i=1,2,3,…,n)作与CE平行的直线Li,以A为圆心,以Rj(j=1,2,3,…,n)为半径作圆弧ωj,ωj与Li相交成n×n个交点Oij,再以Oij为圆心,Rj为半径作圆弧,所作圆弧与煤层底板交点均为A点,和煤层顶板相交成n×n个交点Bij,形成n×n个圆弧                                                 ij,定义圆弧 ij为假定的n×n个滑移面Sij,且定义点Bij到C点的距离CBij为滑移面Sij的滑移深度dij(n×n个);根据煤层采高和采场矿山压力显现情况,估计煤壁的最大深度滑移面的滑移深度dij max,根据煤壁片帮及煤壁加固的工程特点,取最大深度滑移面滑移深度dij max的预测精度为0.05~0.3m,则将dij划分为宽度bijm=0.05~0.3m的k个垂直煤条,将k个垂直煤条由外向里编号,计算或量出各煤条的高度hijm(m=1,2,…,k)和各煤条底部滑移面切线与顶底板平面的夹角 ijm(m=1,2,…,k)。
2)根据式(1)计算顶板作用在每个煤条上的垂直应力qijm
               (1) 
式(1)中,为顶板作用在煤条上的垂直应力,kPa;c为煤的内聚力,kPa;为煤的内摩擦角,°;p为支架护帮板的水平推力,在此可不予考虑,取p=0;xijm为煤条中心距煤壁边缘的距离;为广义Mises准则系数。
3)根据式(2)计算假定滑移面Sij的滑移危险系数Kij。如果Kij>1,所假定的滑移面是稳定的;Kij=1,所假定滑移面是临界滑移面;Kij<1,所假定的滑移面是不稳定的。
                        (2)
式(2)中,为假定滑移面Sij的滑移危险系数;为煤的内摩擦角,°;c为煤的内聚力,kPa;为顶板与煤条接触面上的法向作用力,kN,;bijm为煤条的宽度,m;为工作面沿推进方向的倾角,°,仰斜取正,俯斜取负;fijm为顶板对煤条的摩擦阻力,kN,;Wijm为煤条自重,kN,为煤的容重,kN/m3;xijm为圆心Oij到Fijm的力臂,m,
4)作由yi和Rj所定义的n×n个滑移面Sij(yi,Rj)的滑移危险系数K曲面,用K=1的平面与K曲面相交形成Kcr=1曲线;作n×n个滑移面Sij的滑移深度d曲面,将Kcr=1曲线在d曲面上投影形成临界滑移深度dcr曲线;在dcr曲线上找出最大值dcrmax(ycr,Rcr),由ycr和Rcr所确定的滑移面为最大深度临界滑移面,其dcrmax值为煤壁滑移片帮深度,在大采高工作面煤壁管理中将以此深度设计防治煤壁片帮的控制措施及其参数。
本发明提供的一种确定煤壁滑移片帮深度及其滑移危险性的方法,其特征在于提出一种考虑煤壁高度M、顶板压力F、煤的内聚力c和煤的内摩擦角等多因素的能够更加准确地预测煤壁滑移片帮深度及其滑移危险性的方法,与现有技术《大采高综采工作面煤壁片帮机理及控制技术》(煤炭科学技术,2008,36(9):1-3,24)相比,考虑因素较全面,预测结果更加准确。如:斜沟煤矿18107工作面采高6.2m,煤层的容重为14.4 kN/m3,残余粘聚力c为0.5MPa,内摩擦角为35°,采用现有技术预测该工作面压剪式和滑落式的片帮深度分别为3.23m和8.85m,采用本发明方法预测该工作面煤壁滑移片帮深度为2.2m,而现场实测结果表明,该工作面煤壁片帮深度最大为2.08m,因此,本发明对煤壁滑移片帮深度的预测更加接近实测结果,为控制煤壁片帮方法和参数的选取提供了一种简单、实用和预测结果更加准确的方法,具有重要的科学意义和工程实用价值。
附图说明
附图1是确定工作面煤壁和滑移面方法示意图。
附图2是确定煤壁滑移面垂直煤条划分及其参数方法示意图。
附图3是确定煤壁滑移面滑移危险系数和滑移深度的方法示意图。
附图4是确定最大深度临界滑移面的方法示意图。
附图5是确定某矿一个工作面煤壁-滑移面示意图。
附图6是确定某工作面煤壁滑移面垂直煤条划分示意图。
附图7是确定某工作面煤壁滑移面滑移危险系数和滑移面深度的示意图。
附图8是确定某工作面煤壁最大深度临界滑移面示意图。
具体实施方式
以下结合附图和实例对本发明作进一步的描述:
实施本发明所提供的一种确定煤壁滑移片帮深度及其滑移危险性的方法,其所述方法是按下列步骤进行的:
步骤一、如图5所示,CE为所采煤层顶板,AF为所采煤层底板,AC为工作面煤壁,高度6.2m,在顶板CE以上距CE距离为y3=4m处作平行于CE的直线L3,以A为圆心,半径R3=25m作圆弧ω3和直线L3相交于点O33,以O33为圆心,O33A=R3=25m为半径作圆弧和顶板CE交于点B33,圆弧 33为y3=4m和R3=25m所确定的滑移面S33,再取圆心距顶板高度y1=0m,y2=2m,y4=6m,y5=8m和滑移面半径R1=21m,R2=23m,R4=27m,R5=30m,由yi和Rj组合确定25个滑移面Sij(yi,Rj),见表1,25个滑移面Sij的深度dij(yi,Rj)见表2。根据该煤层采高和采场矿山压力显现情况,估计煤壁的最大深度滑移面的滑移深度dij max为2.5m,根据煤壁片帮及煤壁加固的工程特点,取最大深度滑移面滑移深度的预测精度为0.2m,则将dij划分为宽度bijm=0.2m的k个垂直煤条,将k个垂直煤条由外向里编号,计算或量出各煤条的高度hijm(m=1,2,…,k)和各煤条底部滑移面切线与顶底板平面的夹角 ijm(m=1,2,…,k)。例如图6为滑移面 33上滑移煤体AB33C的煤条划分示意图,共划分k=9个煤条,煤条宽度为0.2m,煤条几何参数见表3。
表1  yi和Rj确定的滑移面
表2  yi和Rj确定的滑移面深度dij/ m
表3  滑移面 33上煤条几何参数
煤条编号m 高度h33m/m 夹角 33m /(°)
1 5.90 67
2 5.37 68
3 4.86 69
4 4.31 71
5 3.72 72
6 3.08 74
7 2.36 75
8 1.54 77
9 0.57 79
步骤二、根据式(1)计算顶板作用在每个煤条上的垂直应力qijm。例如顶板作用在滑移面 33上9个煤条的垂直应力见表4(m=1,2,…,9)。
表4 顶板作用在滑移面 33上煤条的垂直应力
煤条编号m 1 2 3 4 5 6 7 8 9
垂直应力q33m /kPa 1900 2333 2837 3425 4110 4908 5839 6924 8188
步骤三、在工作面取完整性较好的煤样,测得其容重为14.4kN/m3、内聚力为500kPa和内摩擦角为35°;将顶板作用在滑移面 33上9个煤条的垂直应力、煤条的几何参数和物理力学参数代入式(2)中采用迭代法计算滑移面 33的滑移危险系数为0.972153;同理可计算出其余24个滑移面危险系数,整理由yi和Rj所定义的5×5个滑移面Sij的滑移危险系数Kij(yi,Rj)见表5。
表5  假定滑移面上的危险系数Kij
步骤四、如图7所示,作由yi和Rj所定义的5×5个滑移面Sij的滑移危险系数K曲面,用K=1的平面与K曲面相交形成Kcr=1曲线;作5×5个滑移面Sij(yi,Rj)的滑移深度d曲面,将Kcr=1曲线在d曲面上投影形成d曲线,在d曲线上找出最大值2.2m(5.71m,26.65m),由y=5.71m,R=26.65m所确定的滑移面为最大深度临界滑移面,见图8所示,其2.2m为煤壁滑移片帮深度,在大采高工作面煤壁管理中将以此深度设计防治煤壁片帮的控制措施及其参数。

Claims (1)

1.一种确定煤壁滑移片帮深度及其滑移危险性的方法,其特征在于提出一种考虑煤壁高度M、顶板压力F、煤的内聚力c和煤的内摩擦角                                                因素的能够更加准确地确定煤壁滑移片帮深度及其滑移危险性的方法,具体的技术方案为:
Ⅰ、作图给出煤层顶板CE和底板AF,连接AC形成高度为M的煤壁,在顶板CE之上距CE距离为yi(i=1,2,3,…,n)作与CE平行的直线Li,以A为圆心,以Rj(j=1,2,3,…,n)为半径作圆弧ωj,ωj与Li相交成n×n个交点Oij,再以Oij为圆心,Rj为半径作圆弧,所作圆弧与煤层底板交点均为A点,和煤层顶板相交成n×n个交点Bij,形成n×n个圆弧 ij,定义圆弧 ij为假定的n×n个滑移面Sij,且定义点Bij到C点的距离CBij为滑移面Sij的滑移深度dij(n×n个);根据煤层采高和采场矿山压力显现情况,估计煤壁的最大深度滑移面的滑移深度dij max,根据煤壁片帮及煤壁加固的工程特点,取最大深度滑移面滑移深度dij max的预测精度为0.05~0.3m,则将dij划分为宽度bijm=0.05~0.3m的k个垂直煤条,将k个垂直煤条由外向里编号,计算或量出各煤条的高度hijm(m=1,2,…,k)和各煤条底部滑移面切线与顶底板平面的夹角 ijm(m=1,2,…,k);
Ⅱ、根据式(1)计算顶板作用在每个煤条上的垂直应力qijm
           (1)
Ⅲ、根据式(2)计算假定滑移面Sij的滑移危险系数Kij
             (2)
Ⅳ、作由yi和Rj所定义的n×n个滑移面Sij(yi,Rj)的滑移危险系数K曲面,用K=1的平面与K曲面相交形成Kcr=1曲线;作n×n个滑移面Sij的滑移深度d曲面,将Kcr=1曲线在d曲面上投影形成临界滑移深度dcr曲线;在dcr曲线上找出最大值dcrmax(ycr,Rcr),由ycr和Rcr所确定的滑移面为最大深度临界滑移面,其dcrmax值为煤壁滑移片帮深度,在大采高工作面煤壁管理中将以此深度设计防治煤壁滑移片帮的控制措施及其参数。
CN201410354960.8A 2014-07-24 2014-07-24 一种确定煤壁滑移片帮深度及其滑移危险性的方法 Expired - Fee Related CN104268364B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410354960.8A CN104268364B (zh) 2014-07-24 2014-07-24 一种确定煤壁滑移片帮深度及其滑移危险性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410354960.8A CN104268364B (zh) 2014-07-24 2014-07-24 一种确定煤壁滑移片帮深度及其滑移危险性的方法

Publications (2)

Publication Number Publication Date
CN104268364A true CN104268364A (zh) 2015-01-07
CN104268364B CN104268364B (zh) 2018-05-08

Family

ID=52159885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410354960.8A Expired - Fee Related CN104268364B (zh) 2014-07-24 2014-07-24 一种确定煤壁滑移片帮深度及其滑移危险性的方法

Country Status (1)

Country Link
CN (1) CN104268364B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090125141A1 (en) * 2005-11-29 2009-05-14 Takahiro Noda Working apparatus and working method for circuit board
CN103225509A (zh) * 2013-04-09 2013-07-31 中国矿业大学 一种大采高综采采场顶板分类与支护阻力确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090125141A1 (en) * 2005-11-29 2009-05-14 Takahiro Noda Working apparatus and working method for circuit board
CN103225509A (zh) * 2013-04-09 2013-07-31 中国矿业大学 一种大采高综采采场顶板分类与支护阻力确定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHANG HUI 等: "Coal wall collapsing zone analysis and controlling technology", 《ELECTRONIC,COMMUNICATIONS AND CONTROL(ICECC),2011 INTERNATIONAL CONFERENCE ON》 *
刘少伟 等: "倾斜煤层回采巷道上帮煤体滑移危险分析与应用", 《中国矿业大学学报》 *
华心祝: "大采高综采工作面煤壁片帮机理及控制技术", 《煤炭科学技术》 *
徐兵: "大采高工作面煤壁片帮冒顶控制技术", 《辽宁工程技术大学学报(自然科学版)》 *

Also Published As

Publication number Publication date
CN104268364B (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
US20200232320A1 (en) Four-double working method for risk assessment and prediction of water inrush at roof aquifer of coal seam
CN107665285B (zh) 一种基于敏感因子确定遗留煤柱下煤巷合理位置分析方法
CN106801422B (zh) 一种露天矿边坡形状优化设计方法
CN103902780A (zh) 固体充填采煤地表变形预计方法
CN110489928B (zh) 浅埋煤层矿区导水裂隙带发育高度预计方法及系统
CN109255834A (zh) 矿床三维地质储量估算方法
CN112364484A (zh) 一种考虑三维空间效应的邻近既有地铁的基坑开挖影响计算方法
CN106638537A (zh) 一种有限土体土压力计算方法
CN108867666A (zh) 基于开挖变形的结构面控制边坡稳定性评价方法
CN101845815A (zh) 一种软土地层含承压水基坑突涌塑性破坏计算处理方法
CN104268364A (zh) 一种确定煤壁滑移片帮深度及其滑移危险性的方法
CN113221399B (zh) 一种铁路隧道装配式衬砌构件划分方案决策方法
CN107239673A (zh) 一种浅埋煤层群层间双关键层结构支架载荷的确定方法
CN106643610B (zh) 一种煤体塑性区宽度的计算方法
CN108416174B (zh) 大型露天矿山边坡设计安全系数取值方法
Salieiev et al. Development of a methodology for assessing the expediency of mine workings decommissioning based on the geomechanical factor
CN113094905B (zh) 一种适用于多中段连续空区支撑关键点的计算方法及系统
CN113772363B (zh) 一种刮板输送机位姿模型建立方法和系统
CN110096833B (zh) 一种适用于顺层偏压隧道的围岩荷载计算方法
CN109635508A (zh) 一种基于关键层结构的地表偏态下沉系数预判方法
CN110059365A (zh) 隧道与地下空间各种围岩条件下深浅埋临界深度计算方法
CN112943211B (zh) 一种应用于碳酸盐岩油藏的水平布井方法
CN116955957B (zh) 一种煤矿开采工作面顶板导水裂隙带发育高度的预测方法
Marek et al. The reliability analysis of sheet pile wall located in soil with random properties based on CPTu results
CN108875773B (zh) 大型露天矿山边坡稳定性分类方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180508