Roof-to-pressure appearance prediction method and device based on microseismic characteristics

Abstract

The application provides a roof weighting appearance prediction method and device based on microseismic characteristics. The method comprises the following steps: calculating a theoretical value of a weighting step distance according to a numerical value of a rock structure of a roof; correcting the theoretical value of the weighting step distance to obtain a predicted value of the weighting step distance; obtaining an actual mining progress parameter of a working face per day according to schedule data of a production site; calculating a ratio of the predicted value of the weighting step distance and the actual mining progress parameter to obtain a predicted weighting appearance time attribute; and predicting the roof weighting appearance and triggering a warning device when the predicted weighting appearance time attribute is greater than or equal to a current weighting appearance time attribute. The application can analyze and predict the roof weighting appearance of the working face, accurately determine the weighting step distance of the roof, and solve the problem of inaccurate prediction results in the prior art.

Description

Technical Field

[0001] This invention relates to the field of coal mine safety mining technology, and in particular to a method and device for predicting roof pressure manifestation based on microseismic characteristics. Background Technology

[0002] my country's coal seam occurrence conditions are complex and diverse. Working faces with hard roofs exhibit more pronounced signs of rock pressure due to their large roof thickness, high strength, and strong overall integrity. Currently, the main methods for predicting roof pressure indices are manual observation and liquid support resistance monitoring. Manual observation, being based on experience and subjective factors, is highly susceptible to influence, resulting in less objective and accurate predictions. Liquid support resistance monitoring, based on analysis of existing data, also suffers from insufficient accuracy as an early warning system for roof pressure. Summary of the Invention

[0003] This invention provides a method and apparatus for predicting roof inrush manifestation based on microseismic characteristics, in order to solve the problem of insufficient accuracy in the prediction results of roof inrush prediction in the prior art.

[0004] In a first aspect, embodiments of the present invention provide a method for predicting roof pressure manifestation based on microseismic characteristics, including:

[0005] Obtain the rock strata structure values ​​of the top plate, and calculate the theoretical value of the pressure step distance based on the rock strata structure values;

[0006] The theoretical value of the incoming pressure step distance is corrected to obtain the predicted value of the incoming pressure step distance.

[0007] Obtain the schedule data from the production site, and obtain the daily actual mining progress parameters of the working face based on the schedule data;

[0008] Calculate the ratio of the predicted pressure step distance to the actual mining progress parameter to obtain the predicted pressure manifestation time attribute;

[0009] By comparing the predicted pressure display time attribute with the current pressure display time attribute, when the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, the system can predict the pressure display on the top plate and trigger an early warning device.

[0010] In one possible implementation, correcting the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance includes: obtaining the average value of the measured pressure step distance of the support on site, correcting the theoretical value of the pressure step distance based on the average value of the measured pressure step distance, and obtaining the predicted value of the pressure step distance.

[0011] In one possible implementation, the formula for calculating the predicted pressure step distance is:

[0012]

[0013] Among them, L y L0 is the predicted step distance; L0 is the theoretical step distance; K is the average measured step distance of the provided support.

[0014] In one possible implementation, the calculation of the theoretical value of the pressure step distance includes:

[0015] Obtain the numerical values ​​of the roof rock strata structure, and calculate the theoretical value of the initial pressure step distance based on the numerical values ​​of the roof rock strata structure and the preset first formula. The numerical values ​​of the roof rock strata structure include: the width of the fixed cantilever beam, the tensile strength of the rock, and the load per unit length of the cantilever beam.

[0016] The first formula includes:

[0017]

[0018] Among them, L i The initial pressure step distance is the theoretical value; h i Define the width of the cantilever beam fixed at one end of the i-th layer; q represents the tensile strength of the rock. i The load per unit length of the cantilever beam is the load, with the self-weight of the cantilever beam as the load.

[0019] The initial theoretical value of the pressure step distance is corrected by a coefficient to obtain the theoretical value of the pressure step distance.

[0020] In one possible implementation, the step of correcting the initial theoretical value of the pressure step distance with a coefficient to obtain the theoretical value of the pressure step distance includes:

[0021] The initial theoretical value of the pressing step distance is corrected by a coefficient using a preset second formula to obtain the theoretical value of the pressing step distance; the second formula includes:

[0022] L0 = *L i

[0023] Where L0 is the theoretical value of the pressure step distance; μ is the correction coefficient, with a value ranging from... between.

[0024] In one possible implementation, before comparing the predicted time attribute of pressure manifestation with the current time attribute of pressure manifestation, the method further includes: acquiring a table of collected microseismic feature data to obtain a data table with layer height attributes;

[0025] Retrieve a data table within a preset range and use it as the first top plate event group;

[0026] The second roof event group is obtained by classifying the spatial attributes of the first roof event group.

[0027] Correlation attribute analysis was performed on the second roof event group to obtain the roof pressure event; wherein, the microseismic characteristic data table has layer height attribute and spatial attribute.

[0028] In one possible implementation, the method further includes, before comparing the current pressure display time attribute with the preset pressure display time attribute:

[0029] Obtain the data from the first top plate event group;

[0030] Based on the energy magnitude of the first roof event group, a third roof event group with energy attributes is obtained;

[0031] Correlation attribute analysis is performed on the top plate pressure event and the third top plate event group to obtain the first correlation time node group;

[0032] Calculate the average time of the first relevant time node group to obtain the current pressure display time attribute.

[0033] Secondly, embodiments of the present invention provide a top plate pressure manifestation prediction device based on microseismic characteristics, comprising:

[0034] The first calculation module is used to obtain the rock strata structure values ​​of the top plate and calculate the theoretical value of the pressure step distance based on the rock strata structure values.

[0035] The data correction module is used to correct the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance.

[0036] The data acquisition module is used to acquire the schedule data of the production site and obtain the actual daily mining progress parameters of the working face based on the schedule data.

[0037] The second calculation module is used to calculate the ratio of the predicted pressure step distance to the actual mining progress parameter, so as to obtain the predicted pressure manifestation time attribute.

[0038] The early warning module is used to compare the predicted pressure display time attribute with the current pressure display time attribute. When the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, it can predict the pressure display on the top plate and trigger the early warning device.

[0039] Thirdly, embodiments of the present invention provide a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the method as described in the first aspect or any possible implementation of the first aspect.

[0040] Fourthly, embodiments of the present invention provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method as described in the first aspect or any possible implementation thereof.

[0041] The advantages of this invention compared to existing technologies are as follows: In this embodiment, the theoretical value of the pressure step distance is calculated based on the obtained values ​​of the roof rock structure; the theoretical value of the pressure step distance is corrected to obtain the predicted value of the pressure step distance; the actual daily mining progress parameters of the working face are obtained based on the obtained production site schedule data; the ratio of the predicted pressure step distance value to the actual mining progress parameters is calculated to obtain the predicted pressure manifestation time attribute; when the predicted pressure manifestation time attribute is greater than or equal to the current pressure manifestation time attribute, the roof pressure manifestation is predicted and an early warning device is triggered. By obtaining the predicted pressure step distance value and comparing it with the daily actual mining progress parameters of the working face, the predicted pressure manifestation time attribute is obtained. When this predicted value is exceeded, the early warning module is triggered. This invention can analyze and predict the roof pressure display of the working face and accurately determine the roof pressure step distance, solving the problem of inaccurate prediction results in existing technologies. Attached Figure Description

[0042] To more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0043] Figure 1 This is a flowchart illustrating the implementation of the roof pressure manifestation prediction method based on microseismic characteristics provided in this embodiment of the invention.

[0044] Figure 2 This is a schematic diagram of the top plate pressure manifestation prediction device based on microseismic characteristics provided in an embodiment of the present invention;

[0045] Figure 3 This is a schematic diagram of a terminal device provided in an embodiment of the present invention. Detailed Implementation

[0046] In the following description, specific details such as particular system architectures and techniques are set forth for illustrative purposes and not for limitation, in order to provide a thorough understanding of the embodiments of the invention. However, those skilled in the art will understand that the invention can be implemented in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, apparatuses, circuits, and methods are omitted so as not to obscure the description of the invention with unnecessary detail.

[0047] To make the objectives, technical solutions, and advantages of the present invention clearer, specific embodiments will be described below in conjunction with the accompanying drawings.

[0048] Figure 1 This is a flowchart illustrating the implementation of the roof pressure manifestation prediction method based on microseismic characteristics provided in this embodiment of the invention. (Refer to...) Figure 1 The details are as follows:

[0049] In step 101, the rock strata structure values ​​of the top plate are obtained, and the theoretical value of the pressure step distance is calculated based on the rock strata structure values.

[0050] In one possible implementation, calculating the theoretical value of the pressure step distance includes:

[0051] Obtain the numerical values ​​of the roof rock strata structure. Based on the numerical values ​​of the roof rock strata structure and the preset first formula, calculate the theoretical value of the initial pressure step distance. The numerical values ​​of the roof rock strata structure include: the width of the fixed cantilever beam, the tensile strength of the rock, and the load per unit length of the cantilever beam.

[0052] The first formula includes:

[0053]

[0054] Among them, L o The initial pressure step distance is the theoretical value; h i Define the width of the cantilever beam fixed at one end of the i-th layer; q represents the tensile strength of the rock. i This is the load per unit length of the cantilever beam, with the self-weight of the cantilever beam as the load.

[0055] The initial theoretical value of the pressure step distance is corrected by a coefficient to obtain the theoretical value of the pressure step distance.

[0056] In one possible implementation, the theoretical value of the initial pressing step distance is corrected by a coefficient using a pre-defined second formula to obtain the theoretical value of the pressing step distance.

[0057] The second formula includes:

[0058] L0 = *L i

[0059] Where L0 is the theoretical value of the pressure step distance; μ is the correction coefficient, with a value ranging from... between.

[0060] In this embodiment of the invention, the structural values ​​of the top strata include cantilever beam length, width, rock tensile strength, unit load, and other strata structural parameters. The collected data is imported into a first formula to calculate the theoretical value of the initial pressure step distance. This theoretical value is the length of the fixed cantilever beam, calculated as a simply supported beam. If the maximum shear stress is used as the basis for rock fracture, the maximum shear force occurs at both ends of the beam. Therefore, when the maximum shear stress reaches the ultimate strength, the ultimate span is formed. Considering it is a simply supported beam, when the maximum tensile stress equals the ultimate tensile strength, the beam reaches the ultimate span, which is also the length limit of the cantilever beam. The initial pressure step distance of the top strata is comparable to the ultimate span of the initial fracture. Therefore, substituting the rock strata structural values ​​into the length limit of the cantilever beam yields the theoretical value of the initial pressure step distance of the top strata. Experiments show that the theoretical pressure step distance is usually the theoretical value of the initial pressure step distance. By multiplying the initial theoretical value of the indentation step by a factor, the theoretical value of the indentation step L0 can be obtained. The theoretical value of the indentation step is obtained by collecting data on the rock strata structure.

[0061] In step 102, the theoretical value of the pressure step distance is corrected to obtain the predicted value of the pressure step distance.

[0062] In one possible implementation, the average value of the measured pressure step distance on the field support is obtained, and the theoretical value of the pressure step distance is corrected based on the average value of the measured pressure step distance to obtain the predicted value of the pressure step distance.

[0063] In one possible implementation, the formula for calculating the predicted step distance is:

[0064]

[0065] Among them, L y L0 is the predicted step distance; L0 is the theoretical step distance; K is the average measured step distance of the provided support.

[0066] In this embodiment of the invention, the calculated theoretical value of the pressure step distance is corrected. The average measured pressure step distance using the support provided by the mine is used to correct the theoretical value. The ratio of the theoretical pressure step distance L0 to the average measured pressure step distance using the support provided by the mine is calculated as the corrected value of the theoretical pressure step distance. This correction, applied to the theoretical pressure step distance obtained from the rock structure data, yields the predicted pressure step distance, thus improving the accuracy of the predicted pressure step distance.

[0067] In step 103, the production schedule data is obtained, and the actual daily mining progress parameters of the working face are obtained based on the schedule data.

[0068] In this embodiment of the invention, the production schedule data is obtained by collecting roof data using a microseismic monitoring system, and the collected data is transmitted to a ground server. The daily actual mining progress parameter M of the working face is obtained from the collected production schedule data, and the unit of the actual mining progress parameter is meters per day.

[0069] In step 104, the ratio of the predicted pressure step distance to the actual mining progress parameter is calculated to obtain the predicted pressure manifestation time attribute.

[0070] In this embodiment, the predicted pressure step distance is calculated according to a formula, and the unit is meters. The ratio of the predicted pressure step distance to the actual mining progress parameter is in days, which yields the predicted pressure manifestation time attribute, thereby determining the positional relationship between the roof pressure manifestation and the mining progress parameter.

[0071] In step 105, the predicted pressure display time attribute and the current pressure display time attribute are compared. When the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, the system predicts the pressure display on the top plate and triggers the early warning device.

[0072] In this embodiment, the pressure-indicating time attribute is compared with the preset pressure-indicating time attribute. Considering that methods such as calculating the pressure-indicating step distance using a single working face support resistance have a certain lag, this invention analyzes and predicts the relationship between the pressure-indicating time attribute and the current time attribute to predict the pressure-indicating time of the roof slab. The current pressure-indicating time attribute is obtained through correlation attribute analysis based on data collected on-site. When the time attribute of the incoming pressure is displayed, the early warning module is triggered, which realizes prediction and early warning before the top plate pressure is displayed, and solves the current problem of analyzing and predicting the top plate pressure of the working face.

[0073] In summary, this embodiment of the invention calculates the theoretical value of the pressure step distance based on the obtained values ​​of the roof rock structure; corrects the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance; obtains the daily actual mining progress parameters of the working face based on the obtained production site schedule data; calculates the ratio of the predicted pressure step distance value to the actual mining progress parameters to obtain the predicted pressure manifestation time attribute; when the predicted pressure manifestation time attribute is greater than or equal to the current pressure manifestation time attribute, the roof pressure manifestation is predicted and an early warning device is triggered. By obtaining the predicted pressure step distance value and comparing it with the daily actual mining progress parameters of the working face, the predicted pressure manifestation time attribute is obtained. When compared with the current pressure manifestation time attribute, the early warning module is triggered when it exceeds the current value. This invention can analyze and predict the roof pressure display of the working face and accurately determine the roof pressure step distance, solving the problem of inaccurate prediction results in the prior art.

[0074] In one possible implementation, obtaining the event of the top plate pressing down on the layer includes:

[0075] The collected microseismic feature data table is obtained, resulting in a data table with layer height attributes;

[0076] Retrieve a data table within a preset range and use it as the first top plate event group;

[0077] The second roof event group is obtained by classifying the spatial attributes of the first roof event group.

[0078] Correlation attribute analysis was performed on the second roof event group to obtain the roof pressure event; among them, the microseismic characteristic data table has layer height attribute and spatial attribute.

[0079] In this embodiment of the invention, a microseismic monitoring system is used to collect data and transmit it to a ground server to read a microseismic characteristic data table containing layer height and spatial attributes. The data table with layer height attribute H is initially screened, and data with layer height attribute H ≥ 0 are integrated and designated as the first roof event group C0. The first roof event group C0 is then classified according to spatial attributes to obtain a second roof event group, which includes several events. Correlation attribute analysis is performed on the second roof event group within a specific event period to obtain the event group with strong inrush pressure, referred to as roof inrush layer event C1.

[0080] For example, in an embodiment of the present invention, the first top plate event group is classified according to spatial attributes to obtain two sets of layer event sets, which are labeled as group C2 and group C3 respectively; a correlation attribute analysis is performed on group C2 and group C3 with a period of 15 days to obtain events that are strongly correlated with group C0, which are labeled as group C1.

[0081] In one possible implementation, the process of comparing the current pressure-to-display time attribute with the preset pressure-to-display time attribute includes:

[0082] Retrieve data from the first top plate event group;

[0083] Based on the energy magnitude of the first top plate event group, a third top plate event group with energy attributes is obtained;

[0084] Correlation attribute analysis was performed on the roof pressure event and the third roof event group to obtain the first correlation time node group;

[0085] Calculate the average time of the first relevant time node group to obtain the current pressure display time attribute.

[0086] In this embodiment of the invention, microseismic data of the first roof event group C0 is acquired and initially screened. This data group has energy attributes, and within a specified interval, it is grouped according to energy levels into total energy, medium-low energy, and high energy, respectively labeled as group N0, group N1, and group N2. The third roof event group includes total energy N0, medium-low energy N1, and high energy N2. Correlation analysis is performed on the roof pressure event C3 with groups N0, N1, and N2, calculating the correlation weight attributes of the group data. Correlation attribute classification weight curves are analyzed, and multiple time points with strong correlations are identified to obtain the first correlation time node group, labeled as T1…Tn. Outliers are removed, and the average time T0 is obtained, yielding the pressure display time attribute.

[0087] It should be understood that the sequence number of each step in the above embodiments does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0088] The following are device embodiments of the present invention. For details not described in detail, please refer to the corresponding method embodiments described above.

[0089] Figure 2 A schematic diagram of the top plate pressure manifestation prediction device based on microseismic characteristics provided in an embodiment of the present invention is shown. For ease of explanation, only the parts related to the embodiment of the present invention are shown, and are described in detail below:

[0090] like Figure 2 As shown, the roof pressure manifestation prediction device 20 based on microseismic characteristics includes: a first calculation module 21, a data correction module 22, a data acquisition module 23, a second calculation module 24, and an early warning module 25.

[0091] The first calculation module 21 is used to obtain the rock strata structure values ​​of the top plate and calculate the theoretical value of the pressure step distance based on the rock strata structure values.

[0092] Data correction module 22 is used to correct the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance;

[0093] The data acquisition module 23 is used to acquire the schedule data of the production site and obtain the daily actual mining progress parameters of the working face based on the schedule data.

[0094] The second calculation module 24 is used to calculate the ratio of the predicted pressure step distance to the actual mining progress parameter, and obtain the predicted pressure manifestation time attribute.

[0095] The early warning module 25 is used to compare the predicted pressure display time attribute with the current pressure display time attribute. When the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, the module predicts the pressure display on the top plate and triggers the early warning device.

[0096] In one possible implementation, the data correction module 22 is used to correct the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance, including: obtaining the average value of the measured pressure step distance of the support on site, correcting the theoretical value of the pressure step distance based on the average value of the measured pressure step distance, and obtaining the predicted value of the pressure step distance.

[0097] In one possible implementation, the data correction module 22 is also used to calculate the predicted value of the pressure step distance:

[0098]

[0099] Among them, L y L0 is the predicted step distance; L0 is the theoretical step distance; K is the average measured step distance of the provided support.

[0100] In one possible implementation, the first calculation module 21, used to calculate the theoretical value of the pressure step distance, includes:

[0101] Obtain the numerical values ​​of the roof rock strata structure, and calculate the theoretical value of the initial pressure step distance based on the numerical values ​​of the roof rock strata structure and the preset first formula. The numerical values ​​of the roof rock strata structure include: the width of the fixed cantilever beam, the tensile strength of the rock, and the load per unit length of the cantilever beam.

[0102] The first formula includes:

[0103]

[0104] Among them, L o The initial pressure step distance is the theoretical value; h i Define the width of the cantilever beam fixed at one end of the i-th layer; q represents the tensile strength of the rock. i The load per unit length of the cantilever beam is the load, with the self-weight of the cantilever beam as the load.

[0105] The initial theoretical value of the pressure step distance is corrected by a coefficient to obtain the theoretical value of the pressure step distance.

[0106] In one possible implementation, the first calculation module 21 uses a preset second formula to correct the initial theoretical value of the pressing step distance by a coefficient, thereby obtaining the theoretical value of the pressing step distance.

[0107] The second formula includes:

[0108] L0 = *L i

[0109] Where L0 is the theoretical value of the pressure step distance; μ is the correction coefficient, with a value ranging from... between.

[0110] In one possible implementation, the second calculation module 24 is further configured to acquire the top plate pressing layer event, including:

[0111] The collected microseismic feature data table is obtained, resulting in a data table with layer height attributes;

[0112] Retrieve a data table within a preset range and use it as the first top plate event group;

[0113] The second roof event group is obtained by classifying the spatial attributes of the first roof event group.

[0114] Correlation attribute analysis was performed on the second roof event group to obtain the roof pressure event; among them, the microseismic characteristic data table has layer height attribute and spatial attribute.

[0115] In one possible implementation, the second calculation module 24 is further configured to acquire data of the first roof event group; classify the first roof event group according to its energy magnitude to obtain a third roof event group with energy attributes; perform correlation attribute analysis on the roof pressure event and the third roof event group to obtain a first correlation time node group; and calculate the average time of the first correlation time node group to obtain the current pressure manifestation time attribute.

[0116] Figure 3 This is a schematic diagram of a terminal device provided in an embodiment of the present invention. Figure 3 As shown, the terminal device 3 in this embodiment includes: a processor 30, a memory 31, and a computer program 32 stored in the memory 31 and executable on the processor 30. When the processor 30 executes the computer program 32, it implements the steps in the various image segmentation method embodiments described above, for example... Figure 1 Steps 101 to 104 are shown. Alternatively, when the processor 30 executes the computer program 32, it implements the functions of each module / unit in the above-described device embodiments, for example... Figure 2 The functions of modules 110 to 140 are shown.

[0117] The computer program 32 can be divided into one or more modules / units, which are stored in the memory 31 and executed by the processor 30 to complete the present invention. The one or more modules / units can be a series of computer program instruction segments capable of performing specific functions, which describe the execution process of the computer program 32 in the terminal device 3.

[0118] The terminal device 3 can be a desktop computer, laptop, handheld computer, or cloud server, etc. The terminal device may include, but is not limited to, a processor 30 and a memory 31. Those skilled in the art will understand that... Figure 3 This is merely an example of terminal device 3 and does not constitute a limitation on terminal device 3. It may include more or fewer components than shown, or combine certain components, or different components. For example, the terminal device may also include input / output devices, network access devices, buses, etc.

[0119] The processor 30 may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0120] The memory 31 can be an internal storage unit of the terminal device 3, such as a hard disk or memory of the terminal device 3. The memory 31 can also be an external storage device of the terminal device 3, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., equipped on the terminal device 3. Furthermore, the memory 31 can include both internal and external storage units of the terminal device 3. The memory 31 is used to store the computer program and other programs and data required by the terminal device. The memory 31 can also be used to temporarily store data that has been output or will be output.

[0121] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is used as an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. Furthermore, the specific names of the functional units and modules are only for easy differentiation and are not intended to limit the scope of protection of this application. The specific working process of the units and modules in the above system can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0122] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0123] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0124] In the embodiments provided by this invention, it should be understood that the disclosed apparatus / terminal devices and methods can be implemented in other ways. For example, the apparatus / terminal device embodiments described above are merely illustrative. For instance, the division of modules or units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the displayed or discussed mutual coupling or direct coupling or communication connection may be through some interfaces; the indirect coupling or communication connection of devices or units may be electrical, mechanical, or other forms.

[0125] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0126] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0127] If the integrated module / unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording media, USB flash drives, portable hard drives, magnetic disks, optical disks, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc. It should be noted that the content included in the computer-readable medium can be appropriately added or removed according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electrical carrier signals and telecommunication signals.

[0128] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A method for predicting the appearance of pressure from a top plate, characterized in that, include: Obtain the rock strata structure values ​​of the top plate, and calculate the theoretical value of the pressure step distance based on the rock strata structure values; The theoretical value of the incoming pressure step distance is corrected to obtain the predicted value of the incoming pressure step distance. Obtain the schedule data from the production site, and obtain the daily actual mining progress parameters of the working face based on the schedule data; Calculate the ratio of the predicted pressure step distance to the actual mining progress parameter to obtain the predicted pressure manifestation time attribute; By comparing the predicted pressure display time attribute with the current pressure display time attribute, when the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, the system can predict the pressure display on the top plate and trigger an early warning device.

2. The top plate pressure manifestation prediction method as described in claim 1, characterized in that, The theoretical value of the incoming pressure step distance is corrected to obtain the predicted value of the incoming pressure step distance, including: Obtain the average measured pressure step distance of the support on site, and correct the theoretical value of the pressure step distance based on the average measured pressure step distance to obtain the predicted value of the pressure step distance.

3. The method for predicting top plate pressure as described in claim 2, characterized in that, The formula for calculating the predicted pressure step distance is: in, This is the predicted value for the pressure step distance; This is the theoretical value of the pressure step distance; The average step distance was measured for the provided support bracket.

4. The method for predicting top plate pressure manifestation as described in any one of claims 1-3, characterized in that, The calculated theoretical value of the step distance includes: Obtain the numerical values ​​of the roof rock strata structure, and calculate the theoretical value of the initial pressure step distance based on the numerical values ​​of the roof rock strata structure and the preset first formula. The numerical values ​​of the roof rock strata structure include: the width of the fixed cantilever beam, the tensile strength of the rock, and the load per unit length of the cantilever beam. The first formula includes: in, This is the theoretical value for the initial pressure step distance; Define the width of the cantilever beam fixed at one end of the i-th layer; The tensile strength of the rock; The load per unit length of the cantilever beam is the load, with the self-weight of the cantilever beam as the load. The initial theoretical value of the pressure step distance is corrected by a coefficient to obtain the theoretical value of the pressure step distance.

5. The method for predicting top plate pressure as described in claim 4, characterized in that, The step of correcting the initial theoretical value of the pressure step distance to obtain the theoretical value of the pressure step distance includes: The theoretical value of the initial pressing step distance is corrected by a coefficient using a preset second formula to obtain the theoretical value of the pressing step distance. The second formula includes: in, The theoretical value of the pressure step distance is given; μ is a correction coefficient, with a value ranging from... between.

6. The method for predicting top plate pressure as described in claim 1, characterized in that, Before comparing the predicted pressure influx time attribute with the current pressure influx time attribute, the following is also included: The collected microseismic feature data table is obtained, resulting in a data table with layer height attributes; Retrieve a data table within a preset range and use it as the first top plate event group; The second roof event group is obtained by classifying the spatial attributes of the first roof event group. Correlation attribute analysis was performed on the second roof event group to obtain the roof pressure event; wherein, the microseismic characteristic data table has layer height attribute and spatial attribute.

7. The method for predicting top plate pressure as described in claim 6, characterized in that, The method further includes, prior to comparing the predicted pressure display time attribute with the current pressure display time attribute: Obtain the data from the first top plate event group; Based on the energy magnitude of the first roof event group, a third roof event group with energy attributes is obtained; Correlation attribute analysis was performed on the roof pressure event and the third roof event group to obtain the first correlation time node group; Calculate the average time of the first relevant time node group to obtain the current pressure display time attribute.

8. A top-plate pressure-based display prediction device, characterized in that, include: The first calculation module is used to obtain the rock strata structure values ​​of the top plate and calculate the theoretical value of the pressure step distance based on the rock strata structure values. The data correction module is used to correct the theoretical value of the pressure step distance to obtain the predicted value of the pressure step distance. The data acquisition module is used to acquire the schedule data of the production site and obtain the actual daily mining progress parameters of the working face based on the schedule data. The second calculation module calculates the ratio of the predicted pressure step distance to the actual mining progress parameter to obtain the predicted pressure manifestation time attribute. The early warning module is used to compare the predicted pressure display time attribute with the current pressure display time attribute. When the predicted pressure display time attribute is greater than or equal to the current pressure display time attribute, it can predict the pressure display on the top plate and trigger the early warning device.

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, When the processor executes the computer program, it implements the steps of the top plate pressure display prediction method as described in any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the top plate pressure prediction method as described in any one of claims 1 to 7.

Images