CN111828095B - 一种顶板来压动态预测方法 - Google Patents
一种顶板来压动态预测方法 Download PDFInfo
- Publication number
- CN111828095B CN111828095B CN202010657512.0A CN202010657512A CN111828095B CN 111828095 B CN111828095 B CN 111828095B CN 202010657512 A CN202010657512 A CN 202010657512A CN 111828095 B CN111828095 B CN 111828095B
- Authority
- CN
- China
- Prior art keywords
- pressure
- cutting cycle
- mining face
- coal cutting
- fully mechanized
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000003245 coal Substances 0.000 claims abstract description 368
- 238000005065 mining Methods 0.000 claims abstract description 259
- 239000002243 precursor Substances 0.000 claims abstract description 81
- 238000005520 cutting process Methods 0.000 claims description 253
- 238000004364 calculation method Methods 0.000 claims description 73
- 238000012544 monitoring process Methods 0.000 claims description 21
- 238000012216 screening Methods 0.000 claims description 20
- 238000003825 pressing Methods 0.000 claims description 11
- 238000013461 design Methods 0.000 claims description 5
- 238000000547 structure data Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010219 correlation analysis Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0028—Force sensors associated with force applying means
- G01L5/0038—Force sensors associated with force applying means applying a pushing force
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Business, Economics & Management (AREA)
- Analytical Chemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geology (AREA)
- Chemical & Material Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种顶板来压动态预测方法,其特征在于:包括如下步骤:S1、采集综采工作面的地质基础信息;S2、记录顶板来压时的顶板来压信息;S3、搭建来压前兆指标库;S4、选取来压步距最大值Lmax,将综采工作面前方未开采的0~Lmax长度的区域设定为顶板来压初次预测区域;S5、提取顶板来压初次预测区域中综采工作面的基础信息;并从综采工作面后方的已开采区域中提取相近的基础信息,寻找与该相近基础信息相对应的5~8组顶板来压时的顶板来压信息,S6、计算平均值并作为首次预测结果输出;S7、计算综采工作面在回采过程中的来压前兆指标;S8、进行一致性检验;S9、寻找一致性程度最高的记录对下次来压信息进行动态预测更新。
Description
技术领域
本发明涉及矿山压力监测预测技术领域,具体涉及一种顶板来压动态预测方法。
背景技术
目前,煤矿井下综采工作面主要通过实时监测支架压力大小,进行顶板来压实时监测,但是,对综采工作面安全生产的时效性偏低,即综采工作面来压后才能监测到,不能做到事前预测,提前防备。因此,对综采工作面来压进行来压时间、来压地点、来压强度等信息的提前预测意义重大。但是,由于影响顶板来压的因素过多,且各因素间动态变化、相互影响,造成来压预测难度高、预测准确度低,因此,有必要研究一种准确率较高的顶板来压动态预测方法来解决上述问题。
发明内容
本发明目的是针对上述问题,提供一种有效预测顶板下次来压步距、顶板下次来压时间的顶板来压动态预测方法。
为了实现上述目的,本发明的技术方案是:
一种顶板来压动态预测方法,包括如下步骤:
S1、采集综采工作面的地质基础信息;同时根据采煤机的运行轨迹进行割煤循环的划分;所述地质基础信息包括几何尺寸数据、地质构造数据、顶板岩性数据、综采工作面左右两侧区域开采状态数据;综采工作面几何尺寸数据包括综采工作面倾向长度、埋深,地质构造数据包括断层位置、断层落差、断层类型,顶板岩性数据包括综采工作面各坐标点(以开切眼为x轴,上顺槽或下顺槽为y轴,搭建综采工作面平面坐标系,采用坐标确定综采工作面各点位置)的顶板厚度、顶板硬度,综采工作面左右两侧区域开采状态数据包括仅右翼采空、仅左翼采空、两翼均采空、两翼均未采空四种状态。
综采工作面倾向长度为综采工作面的长度,即上、下顺槽的平面中心距。综采工作面埋深即埋藏深度,为综采工作面对应的地面高程与综采工作面底板高程的差值。
S2、建立综采工作面顶板来压监测模型,实时监测综采工作面顶板是否来压并记录顶板来压时的综采工作面的地质基础信息及顶板来压信息;所述顶板来压信息包括来压步距、动载系数、来压部位;
S3、利用步骤S2建立的综采工作面顶板来压监测模型进行来压监测,当监测到顶板来压时,采集综采工作面的来压前兆指标,搭建综采工作面的来压前兆指标库;来压前兆指标库包括来压前兆指标及顶板来压信息;所述来压前兆指标为在顶板来压剩余0.8m、1.6m、2.4m……nm时,对应的支架增阻率ηi、顶板破碎率λi、活柱下缩率ωi、顶底板移近率εi、煤壁片帮率Φi、安全阀开启率ξi以及两次顶板来压之间的综采工作面推进速度v、初撑力利用率δ、额定阻力富余量κi;
S4、提取步骤S2中记录的顶板来压时的来压步距,从中选取出来压步距最大值Lmax,将综采工作面前方未开采的0~Lmax长度的区域设定为顶板来压初次预测区域;
S5、提取综采工作面顶板来压初次预测区域中的地质基础信息;并以预测区域地质基础信息为参照,从综采工作面后方的已开采的顶板来压区域中提取地质基础信息相近的来压数据,来压数据包含地质基础信息及顶板来压信息;
S6、基于步骤S4得到的初次预测区域,在历史综采工作面顶板来压监测模型中寻找与初次预测区域相近地质基础信息相对应的5~8组顶板来压信息,计算顶板来压信息各项数据的平均值,将各项顶板来压信息的平均值作为综采工作面顶板下一次来压的首次预测结果输出;基于当前前综采工作面的推进度,计算剩余来压距离。
S7、伴随综采工作面的不断推进,在产生的新的割煤循环后,在新产生割煤循环期间内,计算综采工作面在回采过程中的来压前兆指标;
S8、将新产生割煤循环期间内的来压前兆指标与步骤S3中采集的综采工作面的来压前兆指标进行一致性检验;
S9、通过一致性检验从综采工作面的来压前兆指标中匹配与当前综采工作面的割煤循环一致性程度最高的记录,以该记录对应来压剩余距离与顶板来压信息作为综采工作面的下次来压信息进行动态预测更新。
将各循环计算的各项指标同历史记录各项指标对比,求各项指标的均方差,最终利用各项指标均方差来判断同当前循环最相近的历史来压记录。并以历史该次记录的来压特征作为当前来压预测标准。
动态预测更新是在首次预测结果之后,综采工作面来压时将进行下一次来压的首次预测更新,之后随割煤循环的推进,将依据循环内指标进行动态更新。
进一步的,所述步骤S1中的综采工作面几何尺寸数据包括综采工作面倾向长度、埋深,地质构造数据包括断层位置、断层落差、断层类型,顶板岩性数据包括综采工作面各坐标点的顶板厚度、顶板硬度,综采工作面左右两侧区域开采状态数据包括仅右翼采空、仅左翼采空、两翼均采空、两翼均未采空四种状态。
进一步的,所述步骤S2中建立综采工作面顶板来压监测模型包括以下步骤:
S21、采集综采工作面的支架压力数据、采高数据,计算各割煤循环内各个支架的最大工作阻力Pmax、时间加权平均阻力Pt、活柱下缩率ω、动载系数m,并分别计算得到该数据的平均值Pmax-aver,Pt-aver,ω-aver,m-aver,将该数据的平均值与该数据的均方差τPmax,τPt,τω,τm进行比较,得到四个顶板来压判据标准;
所述最大工作阻力Pmax取各个支架在割煤循环期间内的最大压力;
所述时间加权平均阻力Pt的计算公式为:
Pt=(P1×T1+P2×T2+…+Pn×Tn)/(T1+T2+…+Tn);
其中,Pt为当前割煤循环内单个支架的时间加权平均阻力;Pn为当前割煤循环内单个支架的第n个压力值;Tn为当前割煤循环内单个支架第n个压力值的持续时间;
所述活柱下缩率ωi的计算公式为:
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1);
其中,ωi为当前割煤循环内单个支架活柱下缩率;hn为当前割煤循环内单个支架的第n个高度值;tn为当前割煤循环内单个支架的第n个高度值产生时间;
所述动载系数m的计算公式为:
m=Pt/Ptt;
其中,m为当前割煤循环内单个支架的动载系数;Pt为当前割煤循环内单个支架的时间加权平均阻力;Ptt为非来压期间内所有支架的时间加权平均阻力;
所述最大工作阻力的平均值Pmax-aver的计算公式为:
Pmax-aver=(Pmax1+Pmax2+…+Pmaxi+…+Pmaxn)/n;
其中,Pmax-aver为当前割煤循环内综采工作面所有支架最大工作阻力的平均值;Pmaxi为当前割煤循环内综采工作面第i个支架的最大工作阻力;n为当前割煤循环内综采工作面全部支架所有支架的最大工作阻力的个数;
所述时间加权平均阻力的平均值Pt-aver的计算公式为:
Pt-aver=(Pt1+Pt2+…+Pti+…+Ptn)/n;
其中,Pt-aver为当前割煤循环内综采工作面所有支架时间加权平均阻力的平均值;Pti为当前割煤循环内综采工作面第i个支架的时间加权平均阻力;n为当前割煤循环内综采工作面所有支架时间加权平均阻力的总数;
所述活柱下缩率的平均值ω-aver的计算公式为:
ω-aver=(ω1+ω2+…+ωi+…+ωn)/n;
其中,ω-aver为当前割煤循环综采工作面所有支架活柱下缩率的平均值;ωi为当前割煤循环第i个支架的活柱下缩率;n为当前割煤循环内活柱下缩率的个数;
所述动载系数的平均值m-aver的计算公式为:
m-aver=(m1+m2+…+mi+…+mn)/n;
其中,m-aver为当前割煤循环综采工作面所有支架动载系数的平均值;mi为当前割煤循环第i个支架的动载系数;n为当前割煤循环内动载系数的个数;
所述最大工作阻力指标均方差τPmax的计算公式为:
其中,τPmax为当前割煤循环综采工作面所有支架的最大工作阻力的均方差;J为综采工作面的支架总数;Pmax_i为当前割煤循环内第i个支架的最大工作阻力;Pmax-aver为当前割煤循环综采工作面所有支架的最大工作阻力的平均值;
所述时间加权平均阻力指标均方差τPt的计算公式为:
其中,τPt为当前割煤循环综采工作面所有支架的时间加权平均阻力指标均方差;J为综采工作面支架总数;Pt_i为当前割煤循环内第i个支架的时间加权平均阻力;Pt-aver为当前割煤循环综采工作面全部支架时间加权平均阻力的平均值;
所述活柱下缩率均方差τω的计算公式为:
其中,τω为当前割煤循环综采工作面所有支架的活柱下缩率均方差;N为当前割煤循环内活柱下缩率的个数;ωi为当前割煤循环内第i个支架的活柱下缩率;ω-aver为当前割煤循环内综采工作面全部支架活柱下缩率的平均值;
所述动载系数均方差τm的计算公式为:
其中,τm为当前割煤循环内综采工作面全部支架的动载系数均方差;N为当前割煤循环内的动载系数个数;mi为当前割煤循环内第i个支架的动载系数;m-aver为当前割煤循环内综采工作面所有支架动载系数的平均值;
判据标准1:最大工作阻力Pmax≥Pmax-aver+τPmax;
判据标准2:时间加权平均阻力Pt≥Pt-aver+τPt;
判据标准3:活柱下缩率ω≥ω-aver+τω;
判据标准4:动载系数m≥m-aver+τm;
S22、来压监测判断标准:伴随综采工作面不断推进,实时计算综采工作面在回采过程中各部位支架的最大工作阻力Pmax、时间加权阻力Pt、活柱下缩率ω、动载系数m;当各支架满足4个判据标准,则将满足标准的支架顶板,记为强烈来压;如果仅满足任意3个判据标准,则将满足标准的支架顶板,记为明显来压;如果其仅满足任意2个判据标准,则将满足标准的支架顶板,记为普通来压;如果其仅满足1个判据标准或1个判据标准都不满足,则将该区域的支架顶板,记为未来压。
进一步的,步骤S3中,所述支架增阻率ηi的计算公式为:
ηi=(P2-P1)/(t2-t1)+(P3-P2)/(t3-t2)+…+(Pn-Pn-1)/(tn-tn-1);
其中,ηi为当前割煤循环内单台支架增阻率;Pn为当前割煤循环内单台支架的第n个压力值;tn为当前割煤循环内单台支架第n个压力的产生时间;
所述顶板破碎率λi的计算公式为:
λi=d1/S1*100%;
其中,λi为当前割煤循环内综采工作面i号支架顶板破碎率;d1为当前割煤循环内端面无支护区冒顶的宽度,单位为m;S1为当前割煤循环内端面无支护区顶板的宽度,单位为m;
所述活柱下缩率ωi计算公式为:
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1);
其中,ωi为当前割煤循环内单台支架活柱下缩率;hn为当前割煤循环内单台支架的第n个高度值;tn为当前割煤循环内单台支架的第n个高度值的产生时间;
所述顶底板移近率εi的计算公式为:
εi=Lshift/(tq-tp);
其中,εi为当前割煤循环内顶底板移近率;Lshift为当前割煤循环内综采工作面顶底板移进变化量;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
所述煤壁片帮率Φi的计算公式为:
Φi=Jslide/J;
其中,Φi为当前割煤循环内煤壁片帮率;Jslide为当前割煤循环内煤壁片帮的支架个数;J为综采工作面的支架总数;
所述安全阀开启率ξi的计算公式为:
ξi=tsafe/(tq-tp);
其中,ξi为当前割煤循环内单台支架安全阀开启率;tsafe为当前割煤循环内单台支架压力处于支架额定安全阀开启压力之上的时长;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
所述两次顶板来压之间的综采工作面推进速度v的计算公式为:
v=ncycle×B/(tw-tm);
其中,v为两次顶板来压之间的综采工作面平均推进速度;ncycle为两次顶板来压之间的割煤循环数;B为采煤机的截割深度;tw为后一次顶板来压的开始时间;tm为前一次顶板来压的结束时间;
所述初撑力利用率δ的计算公式为:
δ=Pset/Prate;
其中,δ为当前割煤循环内单台支架的初撑力利用率;Pset为当前割煤循环单台支架的实际初撑力;Prate为支架的额定初撑力;
所述额定阻力富余量κi的计算公式为:
κi=Pt/P额工作阻力;
其中,κi为当前割煤循环内单台支架的额定阻力富余量;Pt为当前割煤循环单台支架的时间加权平均阻力;P额工作阻力为支架的额定工作阻力。
进一步的,所述步骤S3中煤壁片帮率Φi的计算公式为:
Φi=Jslide/J;
其中,Φi为煤壁片帮率;Jslide为当前割煤循环内煤壁片帮的支架个数;J为当前割煤循环内综采工作面的支架总个数。
进一步的,所述步骤S8中进行一致性检验,包括以下步骤:
S81、从步骤S3综采工作面的来压前兆指标库中,筛选各来压剩余距离下设定数量与当前综采工作面的地质基础信息最为相近的来压前兆指标。以当前综采工作面来压前兆指标为基础,设定各项指标的范围;从采集的来压前兆指标中,去除采集的来压前兆指标中不在设定范围的数据。获取筛选后,剩余各来压剩余距离下来压前兆指标数据。按照来压剩余距离分类,对相同来压剩余距离下的数据将其来压前兆指标求取平均值,生成各来压剩余距离下唯一对应的来压前兆指标;对相同来压剩余距离下的数据将顶板来压信息求取平均值,生成各来压剩余距离下唯一对应的顶板来压信息。
S82、获取通过步骤S7,实时计算当前割煤循环中综采工作面的各支架顶板的来压前兆指标;
S83、对来压前各设定距离下来压前兆指标的均方差进行计算,其中,设定距离为0.8m时支架增阻率的均方差的计算公式为:
其中,τ0.8mη为设定距离为0.8m时当前割煤循环中支架增阻率的均方差;ηi为当前割煤循环中综采工作面的第i个支架增阻率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中支架增阻率的平均值;N为当前割煤循环内综采工作面的支架增阻率总个数;
S84、依据步骤S83中的计算公式分别求解设定距离为0.8m时各项来压前兆指标的均方差,分别为分别为支架增阻率的均方差τ0.8mη、顶板破碎率的均方差τ0.8mλ、支架活柱下缩率的均方差τ0.8mω、顶底板移近率的均方差τ0.8mε、煤壁片帮率的均方差τ0.8mΦ、安全阀开启率的均方差τ0.8mξ、两次顶板来压之间的综采工作面推进速度τ0.8mv、初撑力利用率τ0.8mδ、额定阻力富余量τ0.8mκ,对来压剩余距离0.8m时各项来压前兆指标的均方差,求解各项来压前兆指标的均方差平均值τ0.8m;
支架顶板破碎率的均方差τ0.8mλ的计算公式为:
其中,τ0.8mλ为设定距离为0.8m时当前割煤循环中支架顶板破碎率的均方差;λi为当前割煤循环中综采工作面的第i个支架顶板破碎率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中支架顶板破碎率的平均值;Nλ为当前割煤循环内综采工作面的支架顶板破碎率总个数;
支架活柱下缩率的均方差τ0.8mω的计算公式为:
其中,τ0.8mω为设定距离为0.8m时当前割煤循环中支架活柱下缩率的均方差;ωi为当前割煤循环中综采工作面的第i个支架活柱下缩率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中支架活柱下缩率的平均值;Nω为当前割煤循环内综采工作面的支架活柱下缩率总个数;
顶底板移近率的均方差τ0.8mε的计算公式为:
其中,τ0.8mε为设定距离为0.8m时当前割煤循环中顶底板移近率的均方差;εi为当前割煤循环中综采工作面的第i个支架顶底板移近率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中支架顶底板移近率的平均值;Nε为当前割煤循环内综采工作面的顶底板移近率总个数;
煤壁片帮率的均方差τ0.8mΦ的计算公式为:
其中,τ0.8mΦ为设定距离为0.8m时当前割煤循环中煤壁片帮率的均方差;Φi为当前割煤循环中综采工作面的煤壁片帮率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中综采工作面煤壁片帮率的平均值;
安全阀开启率的均方差τ0.8mξ的计算公式为:
其中,τ0.8mξ为设定距离为0.8m时当前割煤循环中安全阀开启率的均方差;ξi为当前割煤循环中综采工作面的第i个支架安全阀开启率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中综采工作面安全阀开启率的平均值;Nξ为当前割煤循环内综采工作面的安全阀开启率总个数;
两次顶板来压之间的综采工作面推进速度的均方差τ0.8mv计算公式为:
其中,τ0.8mv为设定距离为0.8m时当前割煤循环中综采工作面推进速度的均方差;vi为当前割煤循环中综采工作面的综采工作面推进速度;为步骤S81确定来压前0.8m时,筛选基础前兆指标中综采工作面综采工作面推进速度的平均值;
初撑力利用率的均方差τ0.8mδ的计算公式为:
其中,τ0.8mδ为设定距离为0.8m时当前割煤循环中初撑力利用率的均方差;δi为当前割煤循环中综采工作面的第i个支架初撑力利用率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中综采工作面初撑力利用率的平均值;Nδ为当前割煤循环内综采工作面的初撑力利用率总个数;
额定阻力富余量的均方差τ0.8mκ的计算公式为:
其中,τ0.8mκ为设定距离为0.8m时当前割煤循环中额定阻力富余量的均方差;κi为当前割煤循环中综采工作面的第i个支架额定阻力富余量;为步骤S81确定来压前0.8m时,筛选基础前兆指标中综采工作面额定阻力富余量的平均值;Nκ为当前割煤循环内综采工作面的额定阻力富余量总个数;
通过0.8m时各项来压前兆指标的均方差求解均方差平均值τ0.8m,其公式如下:
τ0.8m=(τ0.8mη+τ0.8mλ+τ0.8mω+τ0.8mε+τ0.8mΦ+τ0.8mξ+τ0.8mv+τ0.8mδ+τ0.8mκ)/9;
其中,τ0.8m为来压剩余距离0.8m时各项来压前兆指标的均方差的平均值。
S85、重复步骤S83、S84,分别求解设定距离为0.8m、1.6m、2.4m…nm时各项来压前兆指标的均方差平均值τ0.8m、τ1.6m、τ2.4m…τnm;
S86、筛选τ0.8m、τ1.6m、τ2.4m…τnm中的最小值τxm,当τxm小于设定数值时,则认为xm为顶板下次来压的剩余来压距离;根据综采工作面割煤循环推进速度,计算顶板下次来压对应的日期;根据综采工作面累计推进距离,计算顶板下次来压对应的来压位置;根据顶板来压信息,预测顶板下次来压对应的强度。
与现有技术相比,本发明具有的优点和积极效果是:
本发明专利基于大数据统计与关联分析,对历史的顶板来压数据进行分类、各指标定量化处理、关联关系拟合等,建立顶板来压历史数据库,通过实时计算当前各指标特征值,并与顶板来压历史数据库进行对比,从而实时、动态为综采工作面来压提供来压预测结果。本发明有效提高了顶板下次来压信息预测准确性,给顶板来压预测工作作出了一定的贡献。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为煤机轨迹及其一次割煤循环划分图;
图2为来压监测、预警流程图;
图3为综采工作面起始布设示意图;
图4为综采工作面回采阶段示意图;
图5为综采工作面破碎状态示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
如图1至图5所示,本实施例中的顶板来压动态预测方法,包括以下步骤:
S1、采集综采工作面的地质基础信息:包括对综采工作面进行地质勘探,在综采工作面上、下顺槽内按照固定间距,综采工作面从开切眼至停产线等间隔设置勘测区域,各勘测区域的顶板岩性、厚度、硬度信息以及通过钻孔揭露的方式,统计综采工作面断层、陷落柱不稳定特殊区域位置信息。通过勘测区域顶板信息,依据岩层传递连续性原则,依据勘探出来的各点顶板信息,延展出综采工作面整体顶板状态;
综采工作面顶板岩性数据库
综采工作面地质构造数据库
S2、根据安装在支架上的红外传感器采集的煤机轨迹,实现割煤循环的划分;根据当前综采工作面割煤进刀工艺,选择设定区域支架号为综采工作面机头进刀区、选择设定区域支架号为综采工作面机尾进刀区。在综采工作面机头进刀区寻找煤机轨迹的峰值点,在综采工作面机尾进刀区寻找煤机轨迹的谷值点,将寻找的两种点定义为煤机轨迹的循环切割点,每相邻两个煤机轨迹循环切割点定义为一次割煤循环。
S3、采集综采工作面支架历史压力数据、采高数据,计算各割煤循环各支架最大工作阻力Pmax、时间加权阻力Pt、活柱下缩率ω、动载系数m,并计算四个指标的平均值Pmax-aver,Pt-aver,ω-aver,m-aver,与四个指标均方差τPmax,τPt,τω,τm,构成四个顶板来压判据:通过来压判据建立综采工作面顶板来压监测模型;
S31、割煤循环各支架最大工作阻力Pmax、时间加权阻力Pt、活柱下缩率ω、动载系数m,并计算四个指标的平均值Pmax-aver,Pt-aver,ω-aver,m-aver的确定。
确定割煤循环的时间期间,取综采工作面各个支架在期间内产生的最大压力Pmax,作为综采工作面各个支架的最大压力数据Pmax。
确定割煤循环内任意一台支架的时间加权平均阻力Pt,计算公式如下,
Pt=(P1×T1+P2×T2+…+Pn×Tn)/(T1+T2+…+Tn);
其中,Pt为当前割煤循环内单个支架的时间加权平均阻力;Pn为当前割煤循环内单个支架的第n个压力值;Tn为当前割煤循环内单个支架第n个压力值的持续时间;
所述活柱下缩率ωi的计算公式为:
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1);
其中,ωi为当前割煤循环内单个支架活柱下缩率;hn为当前割煤循环内单个支架的第n个高度值;tn为当前割煤循环内单个支架的第n个高度值产生时间;
确定割煤循环内任意一台支架的动载系数m,计算公式如下,
m=Pt/Ptt
其中,m为当前割煤循环内单个支架的动载系数;Pt为当前割煤循环内单个支架的时间加权平均阻力;Ptt为非来压期间内所有支架的时间加权平均阻力;
所述最大工作阻力的平均值Pmax-aver的计算公式为:
Pmax-aver=(Pmax1+Pmax2+…+Pmaxi+…+Pmaxn)/n;
其中,Pmax-aver为当前割煤循环内综采工作面所有支架最大工作阻力的平均值;Pmaxi为当前割煤循环内综采工作面第i个支架的最大工作阻力;n为当前割煤循环内综采工作面全部支架所有支架的最大工作阻力的个数;
所述时间加权平均阻力的平均值Pt-aver的计算公式为:
Pt-aver=(Pt1+Pt2+…+Pti+…+Ptn)/n;
其中,Pt-aver为当前割煤循环内综采工作面所有支架时间加权平均阻力的平均值;Pti为当前割煤循环内综采工作面第i个支架的时间加权平均阻力;n为当前割煤循环内综采工作面所有支架时间加权平均阻力的总数;
所述活柱下缩率的平均值ω-aver的计算公式为:
ω-aver=(ω1+ω2+…+ωi+…+ωn)/n;
其中,ω-aver为当前割煤循环综采工作面所有支架活柱下缩率的平均值;ωi为当前割煤循环第i个支架的活柱下缩率;n为当前割煤循环内活柱下缩率的个数;
所述动载系数的平均值m-aver的计算公式为:
m-aver=(m1+m2+…+mi+…+mn)/n;
其中,m-aver为当前割煤循环综采工作面所有支架动载系数的平均值;mi为当前割煤循环第i个支架的动载系数;n为当前割煤循环内动载系数的个数;
所述最大工作阻力指标均方差τPmax的计算公式为:
其中,τPmax为当前割煤循环综采工作面所有支架的最大工作阻力的均方差;J为综采工作面的支架总数;Pmax_i为当前割煤循环内第i个支架的最大工作阻力;Pmax-aver为当前割煤循环综采工作面所有支架的最大工作阻力的平均值;
所述时间加权平均阻力指标均方差τPt的计算公式为:
其中,τPt为当前割煤循环综采工作面所有支架的时间加权平均阻力指标均方差;J为综采工作面支架总数;Pt_i为当前割煤循环内第i个支架的时间加权平均阻力;Pt-aver为当前割煤循环综采工作面全部支架时间加权平均阻力的平均值;
所述活柱下缩率均方差τω的计算公式为:
其中,τω为当前割煤循环综采工作面所有支架的活柱下缩率均方差;N为当前割煤循环内活柱下缩率的个数;ωi为当前割煤循环内第i个支架的活柱下缩率;ω-aver为当前割煤循环内综采工作面全部支架活柱下缩率的平均值;
所述动载系数均方差τm的计算公式为:
其中,τm为当前割煤循环内综采工作面全部支架的动载系数均方差;N为当前割煤循环内的动载系数个数;mi为当前割煤循环内第i个支架的动载系数;m-aver为当前割煤循环内综采工作面所有支架动载系数的平均值;
S33、构成四个液压支架对应的顶板来压判据:
①判据1:最大工作阻力Pmax≥Pmax-aver+τPmax;
②判据2:时间加权平均阻力Pt≥Pt-aver+τPt;
③判据3:活柱下缩率ω≥ω-aver+τω;
④判据4:动载系数m≥m-aver+τm;
S34、来压监测判断标准:若各循环内支架满足①~④中4个判据标准,当各支架满足4个判据标准,则将满足标准的支架顶板,记为强烈来压;如果仅满足任意3个判据标准,则将满足标准的支架顶板,记为明显来压;如果其仅满足任意2个判据标准,则将满足标准的支架顶板,记为普通来压;如果其仅满足1个判据标准或1个判据标准都不满足,则将该区域的支架顶板,记为未来压。
S4、利用建立的来压监测模型,一旦监测到顶板来压,则采集并统计顶板来压前的历史数据,包括综采工作面在顶板来压前0.8m、1.6m、2.4m……nm来压支架对应的支架增阻率ηi、顶板破碎率λi、活柱下缩率ωi、顶底板移近率εi、煤壁片帮率Φi、安全阀开启率ξi,两次顶板来压之间的综采工作面推进速度v、初撑力利用率δ、额定阻力富余量κi,建立综采工作面来压前兆指标库;
来压前兆指标数据库
所述的支架增阻率ηi、顶板破碎率λi、活柱下缩率ωi、顶底板移近率εi、煤壁片帮率Φi、安全阀开启率ξi,两次顶板来压之间的综采工作面推进速度v、初撑力利用率δ、额定阻力富余量κi。
指标计算公式如下所示,
ηi=(P2-P1)/(t2-t1)+(P3-P2)/(t3-t2)+…+(Pn-Pn-1)/(tn-tn-1)
其中,ηi为当前割煤循环内单台支架增阻率;Pn为当前割煤循环内单台支架的第n个压力值;tn为当前割煤循环内单台支架第n个压力的产生时间;
支架的顶板破碎率λi计算手段为人工现场抽样或采用面向顶板空顶区域的激光三维扫描仪进行测量;
人工现场抽样:从1#支架支架开始,设定支架间隔选定支架号,现场实测选定支架中支架上部顶板破碎的信息,包括端面无支护区冒顶的宽度d1,m;端面无支护区顶板的宽度S1,m。
激光三维扫描仪测量:在支架前梁下方布设激光三维扫描仪,实时监测端面无支护区冒顶的宽度d1,m;端面无支护区顶板的宽度S1,m。
顶板破碎率λi计算方法为:
λi=d1/S1*100%
其中,λi为当前割煤循环内综采工作面顶板破碎率;d1为当前割煤循环内端面无支护区冒顶的宽度,单位为m;S1为当前割煤循环内端面无支护区顶板的宽度,单位为m;
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1)
其中,ωi为当前割煤循环内单台支架活柱下缩率;hn为当前割煤循环内单台支架的第n个高度值;tn为当前割煤循环内单台支架的第n个高度值的产生时间;
εi=Lshift/(tq-tp)
其中,εi为当前割煤循环内顶底板移近率;Lshift为当前割煤循环内综采工作面顶底板移进变化量;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
Φi=Jslide/J
其中,Φi为当前割煤循环内煤壁片帮率;Jslide为当前割煤循环内煤壁片帮的支架个数;J为综采工作面的支架总数;
ξi=tsafe/(tq-tp)
其中,ξi为当前割煤循环内安全阀开启率;tsafe为当前割煤循环内支架压力处于支架额定安全阀开启压力之上的时长;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
v=ncycle×B/(tw-tm);
其中,v为两次顶板来压之间的综采工作面平均推进速度;ncycle为两次顶板来压之间的割煤循环数;B为采煤机的截割深度;tw为后一次顶板来压的开始时间;tm为前一次顶板来压的结束时间;
所述初撑力利用率δ的计算公式为:
δ=Pset/Prate;
其中,δ为当前割煤循环内单台支架的初撑力利用率;Pset为当前割煤循环单台支架的实际初撑力;Prate为支架的额定初撑力;
κi=(Pt-P额工作阻力)/P额工作阻力;
其中,κi为当前割煤循环内单台支架的额定阻力富余量;Pt为当前割煤循环单台支架的时间加权平均阻力;P额工作阻力为支架的额定工作阻力。
S5、顶板来压初次预测区域制定;从建立的顶板来压信息数据库中,提取历史来压步距信息,从中选取最大值Lmax,将综采工作面前方未开采的0~Lmax区域作为顶板来压初次预测区域;
S6、以综采工作面基础信息数据库为基础,提取当前综采工作面顶板初次来压预测区域的基础信息,包括:顶板来压初次预测区域内的地质数据、顶板岩性数据、两侧区域开采状态数据;
综采工作面顶板来压的首次预测结果表
时间 | 预测来压步距 | 预测来压时间 | 预测动载系数 | 预测来压部位 |
T1 | L1 | Q1 | K1 | 距开切眼35m |
S7、基于初次来压预测区域内提取的相近的地质基础信息,包括顶板岩性数据、两侧区域开采状态的历史数据,在历史综采工作面顶板来压监测模型中寻找与初次预测区域相近地质基础信息相对应的5~8组顶板来压信息,计算顶板来压信息各项数据的平均值,将各项顶板来压信息的平均值作为综采工作面顶板下一次来压的首次预测结果输出;
S8、随综采工作面不断推进,在产生的新的割煤循环后,在新产生割煤循环期间内,计算来压前兆指标库内的各项前兆指标;
S9、将新产生割煤循环期间内的来压前兆指标与综采工作面来压前兆指标数据库进行数据一致性检验,主要包括如下步骤:
(1)从步骤S3综采工作面的来压前兆指标库中,筛选各来压剩余距离下设定数量与当前综采工作面的地质基础信息最为相近的来压前兆指标。以当前综采工作面来压前兆指标为基础,设定各项指标的范围;从采集的来压前兆指标中,去除采集的来压前兆指标中不在设定范围的数据。获取筛选后,剩余各来压剩余距离下来压前兆指标数据。按照来压剩余距离分类,对相同来压剩余距离下的数据将其来压前兆指标求取平均值,生成各来压剩余距离下唯一对应的来压前兆指标;对相同来压剩余距离下的数据将顶板来压信息求取平均值,生成各来压剩余距离下唯一对应的顶板来压信息。
(2)获取通过步骤S7实时计算当前割煤循环中综采工作面的各支架顶板的来压前兆指标;
(3)分别对来压前各设定距离下指标的均方差计算,以来压前0.8m的支架增阻率ηi计算为例,计算方式如下所示:
其中,τ0.8mη为来压剩余距离0.8m时,来压前兆指标中支架增阻率的均方差;ηi为当前割煤循环中综采工作面的第i个支架增阻率;为来压剩余距离为0.8m时的来压前兆指标中支架增阻率指标的平均值;N为当前割煤循环内综采工作面的支架增阻率总个数;
(4)依据上述计算公式分别求解来压剩余距离0.8m时,来压前兆指标中其余9项指标的均方差,分别为支架增阻率的均方差τ0.8mη、支架增阻率的均方差τ0.8mλ、支架活柱下缩率的均方差τ0.8mω、顶底板移近率的均方差τ0.8mε、煤壁片帮率的均方差τ0.8mΦ、安全阀开启率的均方差τ0.8mξ、两次顶板来压之间的综采工作面推进速度τ0.8mv、初撑力利用率τ0.8mδ、额定阻力富余量τ0.8mκ,对来压剩余距离0.8m时各项来压前兆指标的均方差,求解各项来压前兆指标的均方差平均值τ0.8m;
0.8m时各项来压前兆指标的均方差求解均方差平均值τ0.8m,公式如下:
τ0.8m=(τ0.8mη+τ0.8mλ+τ0.8mω+τ0.8mε+τ0.8mΦ+τ0.8mξ+τ0.8mv+τ0.8mδ+τ0.8mκ)/9
其中,τ0.8m为来压剩余距离0.8m时各项来压前兆指标的均方差的平均值。
(5)依据来压前0.8m均方差以及均方差均值τ0.8m的求解方式,分别求解来压前0.8m、1.6m、2.4m…nm的均方差均值τ0.8m、τ1.6m、τ2.4m…τnm。
(6)筛选τ0.8m、τ1.6m、τ2.4m…τnm中的最小值τxm,当τxm小于设定数值时,则认为xm为顶板下次来压的剩余来压距离;根据综采工作面割煤循环推进速度,计算顶板下次来压对应的日期;根据综采工作面累计推进距离,计算顶板下次来压对应的来压位置;根据顶板来压信息,预测顶板下次来压对应的强度。
S10、通过一致性检验,从综采工作面来压前兆指标数据库,匹配与当前综采工作面割煤循环一致性程度最高的记录,以该记录对应来压剩余距离与顶板来压信息作为综采工作面的下次来压信息进行动态预测更新。跟新方案如下:根据综采工作面割煤循环推进速度,计算顶板下次来压对应的日期;根据综采工作面累计推进距离,计算顶板下次来压对应的来压位置;根据顶板来压信息,预测顶板下次来压对应的强度。
本发明专利基于大数据统计与关联分析,对历史的顶板来压数据进行分类、各指标定量化处理、关联关系拟合等,建立历史综采工作面顶板来压监测模型,通过实时计算当前各指标特征值,并与历史综采工作面顶板来压监测模型进行对比,从而实时、动态进为综采工作面来压提供来压预测结果。本发明有效提高了顶板下次来压信息预测准确性,给顶板来压预测工作作出了一定的贡献。
Claims (6)
1.一种顶板来压动态预测方法,其特征在于:包括如下步骤:
S1、采集综采工作面的地质基础信息;同时根据采煤机的运行轨迹进行割煤循环的划分;所述地质基础信息包括几何尺寸数据、地质构造数据、顶板岩性数据、综采工作面左右两侧区域开采状态数据;
S2、建立综采工作面顶板来压监测模型,实时监测综采工作面顶板是否来压并记录顶板来压时的顶板来压信息;所述顶板来压信息包括来压步距、动载系数、来压部位;
S3、利用步骤S2建立的综采工作面顶板来压监测模型进行来压监测,当监测到顶板来压时,采集综采工作面的来压前兆指标及顶板来压信息,搭建综采工作面的来压前兆指标库;所述来压前兆指标为在顶板来压剩余0.8m、1.6m、2.4m……nm时,对应的支架增阻率ηi、顶板破碎率λi、活柱下缩率ωi、顶底板移近率εi、煤壁片帮率Φi、安全阀开启率ξi以及两次顶板来压之间的综采工作面推进速度v、初撑力利用率δ、额定阻力富余量κi;
S4、提取步骤S2中记录的顶板来压时的来压步距,从中选取出来压步距最大值Lmax,将综采工作面前方未开采的0~Lmax长度的区域设定为顶板来压初次预测区域;
S5、提取综采工作面顶板来压初次预测区域中的地质基础信息;并以该地质基础信息为参照,从综采工作面后方的已开采区域中提取与该地质基础信息相近的5~8组来压数据;
S6、基于步骤S5获取的得到的来压数据与当前综采工作面推进度,计算来压剩余距离以及顶板来压信息各项数据的平均值,将各项顶板来压信息的平均值作为综采工作面顶板下一次来压的首次预测结果输出;
S7、伴随综采工作面的不断推进,在产生的新的割煤循环后,在新产生割煤循环期间内,计算综采工作面在回采过程中的来压前兆指标;
S8、将新产生割煤循环期间内的来压前兆指标与步骤S3中采集的综采工作面的来压前兆指标进行一致性检验;
S9、通过一致性检验从综采工作面的来压前兆指标中匹配与当前综采工作面的割煤循环一致性程度最高的记录,以该记录对应来压剩余距离与顶板来压信息作为综采工作面的下次来压信息进行动态预测更新。
2.如权利要求1所述的一种顶板来压动态预测方法,其特征在于:所述步骤S1中的综采工作面几何尺寸数据包括综采工作面倾向长度、埋深,地质构造数据包括断层位置、断层落差、断层类型,顶板岩性数据包括综采工作面各坐标点的顶板厚度、顶板硬度,综采工作面左右两侧区域开采状态数据包括仅右翼采空、仅左翼采空、两翼均采空、两翼均未采空四种状态。
3.如权利要求2所述的一种顶板来压动态预测方法,其特征在于:所述步骤S2中建立综采工作面顶板来压监测模型包括以下步骤:
S21、采集综采工作面的支架压力数据、采高数据,计算各割煤循环内各个支架的最大工作阻力Pmax、时间加权平均阻力Pt、活柱下缩率ω、动载系数m,并分别计算得到该数据的平均值Pmax-aver,Pt-aver,ω-aver,m-aver,将该数据的平均值与该数据的均方差τPmax,τPt,τω,τm进行比较,得到四个顶板来压判据标准;
所述最大工作阻力Pmax取各个支架在割煤循环期间内的最大压力;
所述时间加权平均阻力Pt的计算公式为:
Pt=(P1×T1+P2×T2+…+Pn×Tn)/(T1+T2+…+Tn);
其中,Pt为当前割煤循环内单个支架的时间加权平均阻力;Pn为当前割煤循环内单个支架的第n个压力值;Tn为当前割煤循环内单个支架第n个压力值的持续时间;
所述活柱下缩率ωi的计算公式为:
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1);
其中,ωi为当前割煤循环内单个支架活柱下缩率;hn为当前割煤循环内单个支架的第n个高度值;tn为当前割煤循环内单个支架的第n个高度值产生时间;
所述动载系数m的计算公式为:
m=Pt/Ptt;
其中,m为当前割煤循环内单个支架的动载系数;Pt为当前割煤循环内单个支架的时间加权平均阻力;Ptt为非来压期间内所有支架的时间加权平均阻力;
所述最大工作阻力的平均值Pmax-aver的计算公式为:
Pmax-aver=(Pmax1+Pmax2+…+Pmaxi+…+Pmaxn)/n;
其中,Pmax-aver为当前割煤循环内综采工作面所有支架最大工作阻力的平均值;Pmaxi为当前割煤循环内综采工作面第i个支架的最大工作阻力;n为当前割煤循环内综采工作面所有支架最大工作阻力的个数;
所述时间加权平均阻力的平均值Pt-aver的计算公式为:
Pt-aver=(Pt1+Pt2+…+Pti+…+Ptn)/n;
其中,Pt-aver为当前割煤循环内综采工作面所有支架时间加权平均阻力的平均值;Pti为当前割煤循环内综采工作面第i个支架的时间加权平均阻力;n为当前割煤循环内综采工作面所有支架时间加权平均阻力的总数;
所述活柱下缩率的平均值ω-aver的计算公式为:
ω-aver=(ω1+ω2+…+ωi+…+ωn)/n;
其中,ω-aver为当前割煤循环综采工作面所有支架活柱下缩率的平均值;ωi为当前割煤循环第i个支架的活柱下缩率;n为当前割煤循环内活柱下缩率的个数;
所述动载系数的平均值m-aver的计算公式为:
m-aver=(m1+m2+…+mi+…+mn)/n;
其中,m-aver为当前割煤循环综采工作面所有支架动载系数的平均值;mi为当前割煤循环第i个支架的动载系数;n为当前割煤循环内动载系数的个数;
所述最大工作阻力指标均方差τPmax的计算公式为:
其中,τPmax为当前割煤循环综采工作面所有支架的最大工作阻力的均方差;J为综采工作面的支架总数;Pmax_i为当前割煤循环内第i个支架的最大工作阻力;Pmax-aver为当前割煤循环综采工作面所有支架的最大工作阻力的平均值;
所述时间加权平均阻力指标均方差τPt的计算公式为:
其中,τPt为当前割煤循环综采工作面所有支架的时间加权平均阻力指标均方差;J为综采工作面支架总数;Pt_i为当前割煤循环内第i个支架的时间加权平均阻力;Pt-aver为当前割煤循环综采工作面全部支架时间加权平均阻力的平均值;
所述活柱下缩率均方差τω的计算公式为:
其中,τω为当前割煤循环综采工作面所有支架的活柱下缩率均方差;N为当前割煤循环内活柱下缩率的个数;ωi为当前割煤循环内第i个支架的活柱下缩率;ω-aver为当前割煤循环内综采工作面全部支架活柱下缩率的平均值;
所述动载系数均方差τm的计算公式为:
其中,τm为当前割煤循环内综采工作面全部支架的动载系数均方差;N为当前割煤循环内的动载系数个数;mi为当前割煤循环内第i个支架的动载系数;m-aver为当前割煤循环内综采工作面所有支架动载系数的平均值;
判据标准1:最大工作阻力Pmax≥Pmax-aver+τPmax;
判据标准2:时间加权平均阻力Pt≥Pt-aver+τPt;
判据标准3:活柱下缩率ω≥ω-aver+τω;
判据标准4:动载系数m≥m-aver+τm;
S22、来压监测判断标准:伴随综采工作面不断推进,实时计算综采工作面在回采过程中各部位支架的最大工作阻力Pmax、时间加权阻力Pt、活柱下缩率ω、动载系数m;当各支架满足4个判据标准,则将满足标准的支架顶板,记为强烈来压;如果仅满足任意3个判据标准,则将满足标准的支架顶板,记为明显来压;如果其仅满足任意2个判据标准,则将满足标准的支架顶板,记为普通来压;如果其仅满足1个判据标准或1个判据标准都不满足,则将该区域的支架顶板,记为未来压。
4.如权利要求3所述的一种顶板来压动态预测方法,其特征在于:步骤S3中,所述支架增阻率ηi的计算公式为:
ηi=(P2-P1)/(t2-t1)+(P3-P2)/(t3-t2)+…+(Pn-Pn-1)/(tn-tn-1);
其中,ηi为当前割煤循环内单台支架增阻率;Pn为当前割煤循环内单台支架的第n个压力值;tn为当前割煤循环内单台支架第n个压力的产生时间;
所述顶板破碎率λi的计算公式为:
λi=d1/S1*100%;
其中,λi为当前割煤循环内综采工作面i号支架顶板破碎率;d1为当前割煤循环内端面无支护区冒顶的宽度,单位为m;S1为当前割煤循环内端面无支护区顶板的宽度,单位为m;
所述活柱下缩率ωi计算公式为:
ωi=(h2-h1)/(t2-t1)+(h3-h2)/(t3-t2)+…+(hn-hn-1)/(tn-tn-1);
其中,ωi为当前割煤循环内单台支架活柱下缩率;hn为当前割煤循环内单台支架的第n个高度值;tn为当前割煤循环内单台支架的第n个高度值的产生时间;
所述顶底板移近率εi的计算公式为:
εi=Lshift/(tq-tp);
其中,εi为当前割煤循环内顶底板移近率;Lshift为当前割煤循环内综采工作面顶底板移进变化量;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
所述煤壁片帮率Φi的计算公式为:
Φi=Jslide/J;
其中,Φi为当前割煤循环内煤壁片帮率;Jslide为当前割煤循环内煤壁片帮的支架个数;J为综采工作面的支架总数;
所述安全阀开启率ξi的计算公式为:
ξi=tsafe/(tq-tp);
其中,ξi为当前割煤循环内单台支架安全阀开启率;tsafe为当前割煤循环内单台支架压力处于支架额定安全阀开启压力之上的时长;tq为当前割煤循环的结束时间;tp为当前割煤循环的开始时间;
所述两次顶板来压之间的综采工作面推进速度v的计算公式为:
v=ncycle×B/(tw-tm);
其中,v为两次顶板来压之间的综采工作面平均推进速度;ncycle为两次顶板来压之间的割煤循环数;B为采煤机的截割深度;tw为后一次顶板来压的开始时间;tm为前一次顶板来压的结束时间;
所述初撑力利用率δ的计算公式为:
δ=Pset/Prate;
其中,δ为当前割煤循环内单台支架的初撑力利用率;Pset为当前割煤循环单台支架的实际初撑力;Prate为支架的额定初撑力;
所述额定阻力富余量κi的计算公式为:
κi=(Pt-P额工作阻力)/P额工作阻力;
其中,κi为当前割煤循环内单台支架的额定阻力富余量;Pt为当前割煤循环单台支架的时间加权平均阻力;P额工作阻力为支架的额定工作阻力。
5.如权利要求4所述的一种顶板来压动态预测方法,其特征在于:所述步骤S3中煤壁片帮率Φi的计算公式为:
Φi=Jslide/J;
其中,Φi为煤壁片帮率;Jslide为当前割煤循环内煤壁片帮的支架个数;J为当前割煤循环内综采工作面的支架总个数。
6.如权利要求5所述的一种顶板来压动态预测方法,其特征在于:所述步骤S8中进行一致性检验,包括以下步骤:
S81、从步骤S3综采工作面的来压前兆指标库中,筛选各来压剩余距离下设定数量与当前综采工作面的地质基础信息最为相近的来压前兆指标;对来压前兆指标进行筛选去除,对筛选后剩余的来压前各设定距离下的来压前兆指标求取平均值,将该平均值作为相近地质条件下来压前设定距离下的来压前兆指标;
S82、通过步骤S7实时计算当前割煤循环中综采工作面的各支架顶板的来压前兆指标;
S83、对来压前各设定距离下来压前兆指标的均方差进行计算,其中,设定距离为0.8m时支架增阻率的均方差的计算公式为:
其中,τ0.8mη为来压剩余距离0.8m时,来压前兆指标中支架增阻率的均方差;ηi为当前割煤循环中综采工作面的第i个支架增阻率;为步骤S81确定来压前0.8m时,筛选基础前兆指标中支架增阻率指标的平均值;N为当前割煤循环内综采工作面的支架增阻率总个数;
S84、依据步骤S83中的计算公式分别求解来压剩余距离0.8m时,来压前兆指标中其余9项指标的均方差,分别为支架增阻率的均方差τ0.8mη、支架增阻率的均方差τ0.8mλ、支架活柱下缩率的均方差τ0.8mω、顶底板移近率的均方差τ0.8mε、煤壁片帮率的均方差τ0.8mΦ、安全阀开启率的均方差τ0.8mξ、两次顶板来压之间的综采工作面推进速度τ0.8mv、初撑力利用率τ0.8mδ、额定阻力富余量τ0.8mκ,对来压剩余距离0.8m时各项来压前兆指标的均方差,求解各项来压前兆指标的均方差平均值τ0.8m;
0.8m时各项来压前兆指标的均方差求解均方差平均值τ0.8m,其计算公式为:
τ0.8m=(τ0.8mη+τ0.8mλ+τ0.8mω+τ0.8mε+τ0.8mΦ+τ0.8mξ+τ0.8mv+τ0.8mδ+τ0.8mκ)/9;
其中,τ0.8m为来压剩余距离0.8m时各项来压前兆指标的均方差的平均值;
S85、重复步骤S83、S84,分别求解设定距离为0.8m、1.6m、2.4m…nm时各项来压前兆指标的均方差平均值τ0.8m、τ1.6m、τ2.4m…τnm;
S86、筛选τ0.8m、τ1.6m、τ2.4m…τnm中的最小值τxm,当τxm小于设定数值时,则认为xm为顶板下次来压的剩余来压距离;根据综采工作面割煤循环推进速度,计算顶板下次来压对应的日期;根据综采工作面累计推进距离,计算顶板下次来压对应的来压位置;根据顶板来压信息,预测顶板下次来压对应的来压强度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010657512.0A CN111828095B (zh) | 2020-07-09 | 2020-07-09 | 一种顶板来压动态预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010657512.0A CN111828095B (zh) | 2020-07-09 | 2020-07-09 | 一种顶板来压动态预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111828095A CN111828095A (zh) | 2020-10-27 |
CN111828095B true CN111828095B (zh) | 2021-11-19 |
Family
ID=72901167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010657512.0A Active CN111828095B (zh) | 2020-07-09 | 2020-07-09 | 一种顶板来压动态预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111828095B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113323698B (zh) * | 2021-05-31 | 2023-11-03 | 国能神东煤炭集团有限责任公司 | 综采工作面推进距离计算方法、存储介质及电子设备 |
CN113482677B (zh) * | 2021-08-02 | 2023-12-22 | 中煤科工开采研究院有限公司 | 一种基于聚类算法的四柱式液压支架工况评估方法 |
CN115795773B (zh) * | 2022-01-20 | 2023-06-23 | 山东科技大学 | 一种煤体灾变时顶板弹性能贡献率影响因素的分析方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103208032A (zh) * | 2013-03-20 | 2013-07-17 | 辽宁工程技术大学 | 一种基于小波和混沌优化lssvm的周期来压预测方法 |
CN104794327A (zh) * | 2015-05-06 | 2015-07-22 | 西安科技大学 | 基于决策融合的多模型矿井顶板安全预警模型 |
CN105298542A (zh) * | 2015-11-18 | 2016-02-03 | 中国神华能源股份有限公司 | 一种用于监测综采工作面顶板的方法及系统 |
CN105756709A (zh) * | 2016-02-24 | 2016-07-13 | 中国矿业大学 | 一种工作面顶板来压及破断的监测方法 |
CN106919744A (zh) * | 2017-02-23 | 2017-07-04 | 陕西煤业化工技术研究院有限责任公司 | 浅埋近距离煤层群开采顶板初次来压支架支护阻力的确定方法 |
CN107559045A (zh) * | 2017-08-24 | 2018-01-09 | 中国矿业大学 | 煤矿井下工作面矿压预测预警系统及方法 |
CN109798149A (zh) * | 2019-01-28 | 2019-05-24 | 辽宁工程技术大学 | 一种坚硬顶板厚煤层采场来压强度分级预测方法 |
CN110889087A (zh) * | 2019-10-29 | 2020-03-17 | 青岛本末岩控技术有限公司 | 一种基于综采工作面的矿压大数据实时预测系统及方法 |
-
2020
- 2020-07-09 CN CN202010657512.0A patent/CN111828095B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103208032A (zh) * | 2013-03-20 | 2013-07-17 | 辽宁工程技术大学 | 一种基于小波和混沌优化lssvm的周期来压预测方法 |
CN104794327A (zh) * | 2015-05-06 | 2015-07-22 | 西安科技大学 | 基于决策融合的多模型矿井顶板安全预警模型 |
CN105298542A (zh) * | 2015-11-18 | 2016-02-03 | 中国神华能源股份有限公司 | 一种用于监测综采工作面顶板的方法及系统 |
CN105756709A (zh) * | 2016-02-24 | 2016-07-13 | 中国矿业大学 | 一种工作面顶板来压及破断的监测方法 |
CN106919744A (zh) * | 2017-02-23 | 2017-07-04 | 陕西煤业化工技术研究院有限责任公司 | 浅埋近距离煤层群开采顶板初次来压支架支护阻力的确定方法 |
CN107559045A (zh) * | 2017-08-24 | 2018-01-09 | 中国矿业大学 | 煤矿井下工作面矿压预测预警系统及方法 |
CN109798149A (zh) * | 2019-01-28 | 2019-05-24 | 辽宁工程技术大学 | 一种坚硬顶板厚煤层采场来压强度分级预测方法 |
CN110889087A (zh) * | 2019-10-29 | 2020-03-17 | 青岛本末岩控技术有限公司 | 一种基于综采工作面的矿压大数据实时预测系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111828095A (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111828095B (zh) | 一种顶板来压动态预测方法 | |
Wu et al. | Real-time rock mass condition prediction with TBM tunneling big data using a novel rock–machine mutual feedback perception method | |
AU2019312330B2 (en) | Method for predicting mine strata pressure behavior data of stoping tunnel | |
CN104533400B (zh) | 一种重构测井曲线的方法 | |
CN104747185A (zh) | 非均质油藏储层综合分类评价方法 | |
WO2020199495A1 (zh) | 岩层参数的确定方法及装置 | |
CN107194068A (zh) | 页岩气压裂过程井下异常工况实时预测预警方法及装置 | |
CN112614021B (zh) | 一种基于已建隧道信息智能识别的隧道围岩地质信息预测方法 | |
CN115422740A (zh) | 一种巨厚煤层分层综放开采导水裂隙带高度预测方法 | |
CN115293316A (zh) | 基于ssa-elman中深埋厚煤层导水裂隙带预测方法 | |
CN112085242B (zh) | 一种基于大数据深度学习方法的低阻层精准预测方法与装置 | |
CN114357750A (zh) | 一种采空区充水状态评估方法 | |
Qiu et al. | TBM tunnel surrounding rock classification method and real-time identification model based on tunneling performance | |
CN115438823A (zh) | 一种井壁失稳机制分析与预测方法及系统 | |
CN113323676A (zh) | 用主成分分析-长短记忆模型确定盾构机刀盘扭矩的方法 | |
CN109268003A (zh) | 一种地质工程一体化框架下的致密油增产数学建模方法 | |
CN111340275B (zh) | 基于随钻探测技术的隧道支护模式选型实时预测方法 | |
CN116432855A (zh) | 一种基于掘进数据的隧道塌方情况的预测方法 | |
CN117131971A (zh) | 基于Xgboost算法的储层优势渗流通道预测方法 | |
CN115964667A (zh) | 基于深度学习和重采样的河流-湖泊岩相测井识别方法 | |
CN112906143B (zh) | 一种考虑数据分布域适应的综采工作面矿压预测模型建立方法 | |
CN115046516A (zh) | 基于单滑面r型深孔测斜曲线的滑动面位置精准确定方法 | |
CN115049173A (zh) | 深度学习和Eaton法耦合驱动地层孔隙压力预测方法 | |
Blouin et al. | Automated facies prediction in drillholes using machine learning | |
CN115438090A (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 |