CN-115952910-B - Roof pressure display prediction method and device based on microseismic characteristics

CN115952910BCN 115952910 BCN115952910 BCN 115952910BCN-115952910-B

Abstract

The invention provides a method and a device for predicting roof pressure display based on microseismic characteristics. The method comprises the steps of calculating a theoretical value of a running step distance according to a rock stratum structure value of a top plate, correcting the theoretical value of the running step distance to obtain a predicted value of the running step distance, obtaining a daily actual stoping progress parameter of a working face according to schedule data of a production site, calculating the ratio of the predicted value of the running step distance to the actual stoping progress parameter to obtain a predicted running display time attribute, and predicting the running display of the top plate and triggering an early warning device when the predicted running display time attribute is greater than or equal to the current running display time attribute. The invention can analyze and predict the top plate pressing display of the working surface, accurately judge the top plate pressing step distance and solve the problem of inaccurate prediction result in the prior art.

Inventors

ZHAO LISONG
ZHANG SHUJIAO
GUO XUETING
JIA JING
SHANG BIN
YANG BO
WANG PENG
ZHOU JINYAN
LI JIA
GAO GANG
LIU RUI

Assignees

河北煤炭科学研究院有限公司

Dates

Publication Date: 20260508
Application Date: 20230103

Claims (10)

1. A roof pressing display prediction method, comprising: Acquiring a rock stratum structure value of a top plate, and calculating a theoretical value of a step-pressing distance according to the rock stratum structure value; Correcting the theoretical value of the incoming pressure step distance to obtain a predicted value of the incoming pressure step distance; Acquiring schedule data of a production site, and acquiring actual extraction schedule parameters of a working face every day according to the schedule data; Calculating the ratio of the predicted value of the step distance to the actual stoping progress parameter to obtain a predicted step showing time attribute; and comparing the predicted coming pressure display time attribute with the current coming pressure display time attribute, and when the predicted coming pressure display time attribute is greater than or equal to the current coming pressure display time attribute, realizing the prediction of the top plate before coming pressure display and triggering the early warning device.
2. The roof pressure appearance prediction method of claim 1, wherein correcting the theoretical value of the step distance to obtain the predicted value of the step distance includes: and acquiring an actual measurement incoming pressure step average value of the field bracket, and correcting the incoming pressure step theoretical value according to the actual measurement incoming pressure step average value to obtain the incoming pressure step predicted value.
3. The roof depression prediction method as claimed in claim 2, wherein the calculation formula of the depression step pitch prediction value is: Wherein, the Predicted value of step distance for coming and pressing; A theoretical value for the step distance; The step average is measured for the stent provided.
4. A roof crush show prediction method as claimed in any one of claims 1 to 3 wherein the calculating a theoretical value of the crush step comprises: Obtaining a roof stratum structure value, and calculating a primary running step pitch theoretical value according to the roof stratum structure value and a preset first formula, wherein the roof stratum structure value comprises the width of a fixed cantilever beam, the tensile strength of rock and the load of the unit length of the cantilever beam; wherein the first formula comprises: Wherein, the The theoretical value of the initial step distance is the theoretical value of the initial step distance; Fixing the width of the cantilever beam for one end of the ith layer; Is the tensile strength of the rock; The load of the cantilever beam unit length is the self weight of the cantilever beam as the load; And carrying out coefficient correction on the initial step pitch theoretical value to obtain the step pitch theoretical value.
5. The roof crush display prediction method according to claim 4, wherein the performing coefficient correction on the initial crush step theoretical value to obtain the crush step theoretical value comprises: carrying out coefficient correction on the primary pressing step distance theoretical value by using a preset second formula to obtain a pressing step distance theoretical value; the second formula includes: Wherein, the The value range is that the theoretical value of the step distance is calculated, mu is the correction coefficient Between them.
6. The roof impression prediction method of claim 1, further comprising, prior to comparing the predicted impression time attribute with a current impression time attribute: acquiring an acquired microseismic characteristic data table to obtain a data table with horizon height attribute; acquiring a data table within a preset range as a first roof event group; classifying according to the spatial attribute of the first roof event group to obtain a second roof event group; and performing correlation attribute analysis on the second roof event group to obtain roof pressure event, wherein the microseismic characteristic data table has horizon height attribute and spatial attribute.
7. The roof impression prediction method of claim 6, further comprising, prior to said comparing said predicted impression time attribute with a current impression time attribute: Acquiring data of the first roof event group; Classifying according to the energy of the first roof event group to obtain a third roof event group with energy attribute; performing correlation attribute analysis on the roof pressing event and the third roof event group to obtain a first correlation time node group; and calculating the average time of the first correlation time node group to obtain the current incoming call presentation time attribute.
8. A roof pressing appearance prediction device, comprising: the first calculation module is used for obtaining the rock stratum structure value of the top plate and calculating a theoretical value of the step-pressing distance according to the rock stratum structure value; the data correction module is used for correcting the theoretical value of the incoming pressure step distance to obtain a predicted value of the incoming pressure step distance; the data acquisition module is used for acquiring schedule data of a production site and acquiring actual extraction progress parameters of a working face every day according to the schedule data; The second calculation module calculates the ratio of the predicted value of the step length to the actual stoping progress parameter to obtain a predicted step length showing time attribute; And the early warning module is used for comparing the predicted incoming call display time attribute with the current incoming call display time attribute, and realizing the prediction and triggering of the early warning device before the roof incoming call display when the predicted incoming call display time attribute is greater than or equal to the current incoming call display time attribute.
9. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the roof compression presentation prediction method as claimed in any one of the preceding claims 1 to 7 when the computer program is executed.
10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the roof pressure representation prediction method of any one of the preceding claims 1 to 7.

Description

Roof pressure display prediction method and device based on microseismic characteristics Technical Field The invention relates to the technical field of coal mine safety exploitation, in particular to a roof pressure display prediction method and device based on microseismic characteristics. Background The coal seam occurrence conditions in China are complex and various, and the hard roof working face has the characteristics of large roof thickness, high strength and strong integrity, so that the mine pressure display degree of the working face is stronger. The method for predicting the working face roof pressure display mainly comprises a manual observation method and a liquid support working resistance monitoring method, wherein the manual observation method belongs to a method based on experience judgment, the influence of subjective factors of people is large, the prediction result is not objective and accurate enough, the liquid support working resistance monitoring method belongs to a judgment method for analyzing based on generated data, and the method also has the defect of insufficient accuracy when being used as an early warning mode before roof pressure. Disclosure of Invention The embodiment of the invention provides a roof pressure display prediction method and device based on microseismic characteristics, which are used for solving the problem of insufficient accuracy of a prediction result of roof pressure prediction in the prior art. In a first aspect, an embodiment of the present invention provides a method for predicting a roof pressure display based on microseismic characteristics, including: Acquiring a rock stratum structure value of a top plate, and calculating a theoretical value of a step-pressing distance according to the rock stratum structure value; Correcting the theoretical value of the incoming pressure step distance to obtain a predicted value of the incoming pressure step distance; Acquiring schedule data of a production site, and acquiring actual extraction schedule parameters of a working face every day according to the schedule data; Calculating the ratio of the predicted value of the step distance to the actual stoping progress parameter to obtain a predicted step showing time attribute; and comparing the predicted coming pressure display time attribute with the current coming pressure display time attribute, and when the predicted coming pressure display time attribute is greater than or equal to the current coming pressure display time attribute, realizing the prediction of the top plate before coming pressure display and triggering the early warning device. In one possible implementation manner, correcting the theoretical value of the incoming pressure step distance to obtain the predicted value of the incoming pressure step distance comprises the steps of obtaining an actual measurement incoming pressure step distance average value of a field bracket, and correcting the theoretical value of the incoming pressure step distance according to the actual measurement incoming pressure step distance average value to obtain the predicted value of the incoming pressure step distance. In one possible implementation manner, the calculation formula of the step pitch prediction value is: Wherein L y is a predicted value of the step length of the step, L 0 is a theoretical value of the step length of the step, and K is an average value of the step length of the support after actual measurement. In one possible implementation, the calculating the theoretical step size value includes: Obtaining a roof stratum structure value, and calculating a primary running step pitch theoretical value according to the roof stratum structure value and a preset first formula, wherein the roof stratum structure value comprises the width of a fixed cantilever beam, the tensile strength of rock and the load of the unit length of the cantilever beam; wherein the first formula comprises: wherein L i is a theoretical value of the initial step distance, h i is the width of a fixed cantilever beam at one end of the ith layer; q i is the load of the cantilever beam unit length, and the dead weight of the cantilever beam is used as the load; And carrying out coefficient correction on the initial step pitch theoretical value to obtain the step pitch theoretical value. In one possible implementation manner, the performing coefficient correction on the initial step distance theoretical value to obtain the step distance theoretical value includes: And carrying out coefficient correction on the primary initial step distance theoretical value by using a preset second formula to obtain the initial step distance theoretical value, wherein the second formula comprises the following steps: L0＝*Li wherein L 0 is the theoretical value of the step pitch of the incoming pressure, mu is the correction coefficient, and the value range is that Between them. In one possible implementation manner, before comparing the predicte