Roof decompression determination method, device and storage medium

Abstract

The present invention provides a method, device and storage medium for determining roof pressure relief, including determining a first effective microseismic event activity interval in a plane section of the coal seam roof and a second effective microseismic event activity interval in a vertical section of the coal seam roof within the three-dimensional space of the coal seam roof of an adjacent working face based on the energy and frequency of each microseismic event within the three-dimensional space of the coal seam roof of an adjacent working face; and determining the spatial position of the roof pressure relief of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face. The present invention accurately defines the effective microseismic event activity intervals in the plane direction and the vertical direction of the adjacent working face by analyzing the energy and frequency of each microseismic event within the three-dimensional space of the coal seam roof of an adjacent working face, and then determines the spatial position of the roof pressure relief of the current working face based on these effective microseismic event activity intervals of the adjacent working face, thereby improving the accuracy of determining the roof pressure relief position.

Description

Technical Field

[0001] The present invention relates to the technical field of coal mining, and in particular to a method, device and storage medium for determining roof pressure relief. Background Art

[0002] With the increasing intensity of coal mining, the depth of coal resources has expanded from shallow to deep, and the mining environment has become more complex and variable. Under the combined influence of high ground stress and strong mining stress, the overlying hard and thick roof easily causes significant mine pressure in the working face area.

[0003] To alleviate this problem, downhole short-hole fracturing, long horizontal hole fracturing, and surface fracturing are currently widely used to pre-break the coal seam roof. Therefore, precisely determining the spatial location of the fracturing is crucial, as it not only improves fracturing accuracy but also effectively reduces the workload of the fracturing project.

[0004] Although existing technologies can use geological drilling columns and lithology monitoring in different areas of the roof to determine the roof unloading layer, these methods still have limitations in accurately locating the specific spatial position of the roof unloading and cannot meet high-precision requirements. Summary of the Invention

[0005] The present invention provides a method, device and storage medium for determining roof pressure relief, which are used to solve the technical defect in the prior art that the specific spatial position of roof pressure relief cannot be accurately located.

[0006] The present invention provides a method for determining roof pressure relief, comprising the following steps.

[0007] Based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face are determined;

[0008] Based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face, the roof decompression space position of the current working face is determined.

[0009] According to a method for determining roof pressure relief provided by the present invention, the method determines a first effective microseismic event activity interval of a coal seam roof plane interval within the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event within the three-dimensional space of the coal seam roof adjacent to the working face, comprising:

[0010] Determine a first total energy and a first total frequency of microseismic events in each plane subinterval divided in a plane interval of a coal seam roof in a three-dimensional space of a coal seam roof adjacent to a working face;

[0011] Determining a plane subinterval energy threshold based on the first total energy of the microseismic events in each plane subinterval in the coal seam roof plane interval;

[0012] Determining a plane subinterval frequency threshold based on the first total frequency of microseismic events in each plane subinterval in the coal seam roof plane interval;

[0013] Determine a first target plane subinterval in which the first total energy is greater than the plane subinterval energy threshold and a second target plane subinterval in which the first total frequency is greater than the plane subinterval frequency threshold;

[0014] The first target plane sub-interval and the second target plane sub-interval are determined as first effective microseismic event activity intervals in the coal seam roof plane interval in the three-dimensional space of the coal seam roof adjacent to the working face.

[0015] According to a method for determining roof pressure relief provided by the present invention, the method further includes:

[0016] Based on the width of the adjacent working face transport lane, the adjacent working face width, the adjacent working face return air lane width and the adjacent working face return air lane side coal pillar width, the coal seam roof plane interval in the three-dimensional space of the adjacent working face coal seam roof is determined;

[0017] The coal seam roof plane interval is divided into a plurality of plane sub-intervals according to a first preset interval length.

[0018] According to a method for determining roof pressure relief provided by the present invention, the method determines a second effective microseismic event activity interval of a vertical interval of a coal seam roof in the three-dimensional space of a coal seam roof adjacent to a working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, including:

[0019] Determine the second total energy and the second total frequency of the microseismic events in each vertical sub-interval divided in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face;

[0020] determining a vertical sub-interval energy threshold based on the second total energy of the microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof;

[0021] determining a vertical sub-interval frequency threshold based on the second total frequency of microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof;

[0022] Determine a first target vertical subinterval in which the second total energy is greater than the vertical subinterval energy threshold and a second target vertical subinterval in which the second total frequency is greater than the vertical subinterval frequency threshold;

[0023] The first target vertical sub-interval and the second target vertical sub-interval are determined as second effective microseismic event activity intervals in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof of the adjacent working face.

[0024] According to a method for determining roof pressure relief provided by the present invention, the method further includes:

[0025] Based on the height of the coal seam roof fracture zone, the vertical interval of the coal seam roof in the three-dimensional space of the adjacent working face is determined;

[0026] The vertical interval of the coal seam roof is divided into a plurality of vertical sub-intervals according to a second preset interval length.

[0027] According to a method for determining roof decompression provided by the present invention, determining the spatial position of roof decompression of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face includes:

[0028] Determine a second interval in a plane interval of a coal seam roof within a three-dimensional space of a coal seam roof adjacent to a working face, excluding the first effective microseismic event activity interval;

[0029] Based on the second interval of the adjacent working surface, determining a roof pressure relief area in a plane direction of the working surface;

[0030] Based on the second effective microseismic event activity interval of the adjacent working face, the vertical roof decompression area of ​​the working face is determined.

[0031] The present invention also provides a device for determining roof pressure relief, comprising:

[0032] A first roof decompression determination module is configured to determine, based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof within the three-dimensional space of the coal seam roof adjacent to the working face;

[0033] The second roof pressure relief determination module is used to determine the roof pressure relief spatial position of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face.

[0034] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, any of the above-described methods for determining top plate decompression is implemented.

[0035] The method for determining roof pressure relief provided by the present invention includes determining a first effective microseismic event activity interval in a plane section of the coal seam roof and a second effective microseismic event activity interval in a vertical section of the coal seam roof within the three-dimensional space of the coal seam roof of the adjacent working face based on the energy and frequency of each microseismic event within the three-dimensional space of the coal seam roof of the adjacent working face; and determining the spatial position of the roof pressure relief of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face. The present invention accurately defines the effective microseismic event activity intervals in the plane direction and the vertical direction of the adjacent working face by analyzing the energy and frequency of each microseismic event within the three-dimensional space of the coal seam roof of the adjacent working face, and then determines the spatial position of the roof pressure relief of the current working face based on these effective microseismic event activity intervals of the adjacent working face, thereby improving the accuracy of determining the roof pressure relief position. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] In order to more clearly illustrate the technical solutions in the present invention or the prior art, a brief introduction is given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0037] Figure 1 It is a flow chart of the method for determining roof pressure relief provided by the present invention.

[0038] Figure 2 This is one of the scenario schematic diagrams of the method for determining roof pressure relief provided by the present invention.

[0039] Figure 3 This is the second scenario diagram of the method for determining roof pressure relief provided by the present invention.

[0040] Figure 4 This is the third scenario diagram of the method for determining roof pressure relief provided by the present invention.

[0041] Figure 5 It is a structural schematic diagram of a method and device for determining roof pressure relief provided by the present invention.

[0042] Figure 6 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0043] To make the objectives, technical solutions, and advantages of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0044] Figure 1 FIG. 1 is a flow chart of a method for determining roof pressure relief provided by the present invention, as shown in FIG. Figure 1 As shown, the method includes step 210 and step 220.

[0045] Step 210: Determine a first effective microseismic event activity interval in a plane section of the coal seam roof and a second effective microseismic event activity interval in a vertical section of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face;

[0046] By statistically analyzing the distribution of microseismic event energy and frequency within the three-dimensional space of the coal seam roof of adjacent working faces (mined sections), the first effective microseismic event activity interval of the adjacent working face can be determined. The first effective microseismic event activity interval of the adjacent working face is the interval in the coal seam roof plane where microseismic events are most concentrated and have higher energy.

[0047] Similarly, in the three-dimensional space of the coal seam roof of the adjacent working face, by statistically analyzing the distribution of the energy and frequency of microseismic events in the vertical interval of the coal seam roof of the adjacent working face, the second effective microseismic event activity interval of the adjacent working face can be determined.

[0048] Step 220: Determine the roof decompression space position of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face.

[0049] Finally, the roof decompression space position of the current working face (the section currently being mined) is determined by combining the first effective microseismic event activity interval in the plane direction of the adjacent working face and the second effective microseismic event activity interval in the vertical direction.

[0050] Specifically, determine a second interval in the plane interval of the coal seam roof within the three-dimensional space of the coal seam roof of the adjacent working face, excluding the first effective microseismic event activity interval; based on the second interval of the adjacent working face, determine the roof unloading area in the plane direction of the working face; based on the second effective microseismic event activity interval of the adjacent working face, determine the roof unloading area in the vertical direction of the working face.

[0051] It should be understood that microseismic activity in the area outside the first effective microseismic event activity interval in the planar direction of the coal seam roof is not significant. Therefore, in this embodiment, the planar roof pressure relief area of ​​the working face is determined based on the second interval of the coal seam roof planar interval of the adjacent working face, excluding the first effective microseismic event activity interval. In the vertical direction, in this embodiment, the vertical roof pressure relief area of ​​the working face is determined based on the second effective microseismic event activity interval of the adjacent working face.

[0052] Specifically, refer to Figure 3 As shown in the figure, the first effective microseismic event activity interval (i.e. Figure 3 After determining the blue dense interval range in the horizontal direction shown in the figure, determine the second interval in the plane interval of the coal seam roof of the adjacent working face, and then use this second interval as a reference to find the corresponding area in the plane interval of the coal seam roof of this working face and define it as the plane direction roof unloading area.

[0053] It should be understood that the coal seam roof plane intervals of this working face and the adjacent working face follow the same division principles. In addition, the horizontal position of the plane-direction roof unloading zone in this working face within the coal seam roof plane interval must be consistent with the horizontal position of the second interval in the adjacent working face within its coal seam roof plane interval, ensuring the relative consistency of the two in the horizontal dimension.

[0054] Similarly, the second effective microseismic event activity interval (such as Figure 4 After taking the second effective microseismic event activity interval as a benchmark, a corresponding interval is found and determined within the vertical interval of the coal seam roof of this working face, and it is marked as the vertical roof unloading area.

[0055] It should be noted that the vertical divisions of the coal seam roof vertical intervals are consistent, both in the current working face and in the adjacent working face. Furthermore, the vertical position of the roof pressure relief area determined in the current working face within its coal seam roof vertical interval matches the vertical position of the second effective microseismic event interval in the adjacent working face within the corresponding coal seam roof vertical interval, ensuring relative vertical consistency between the two.

[0056] The method for determining roof pressure relief provided by an embodiment of the present invention analyzes the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof of the adjacent working face, accurately defines the effective microseismic event activity intervals in the planar direction and vertical direction of the adjacent working face, and then determines the spatial position of the roof pressure relief of the current working face based on these effective microseismic event activity intervals of the adjacent working face, thereby improving the accuracy of determining the roof pressure relief position.

[0057] In some embodiments, determining a first effective microseismic event activity interval of a coal seam roof plane interval in the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face includes:

[0058] Determine a first total energy and a first total frequency of microseismic events in each plane subinterval divided in a plane interval of a coal seam roof in a three-dimensional space of a coal seam roof adjacent to a working face;

[0059] Determining a plane subinterval energy threshold based on the first total energy of the microseismic events in each plane subinterval in the coal seam roof plane interval;

[0060] Determining a plane subinterval frequency threshold based on the first total frequency of microseismic events in each plane subinterval in the coal seam roof plane interval;

[0061] Determine a first target plane subinterval in which the first total energy is greater than the plane subinterval energy threshold and a second target plane subinterval in which the first total frequency is greater than the plane subinterval frequency threshold;

[0062] The first target plane sub-interval and the second target plane sub-interval are determined as first effective microseismic event activity intervals in the coal seam roof plane interval in the three-dimensional space of the coal seam roof adjacent to the working face.

[0063] In the three-dimensional space of the coal seam roof adjacent to the working face, the plane interval of the coal seam roof adjacent to the working face is first divided into several plane subintervals. For each plane subinterval, the first total energy and first total frequency of all microseismic events occurring in the plane subinterval are calculated.

[0064] Afterwards, the plane sub-interval energy threshold of the coal seam roof plane interval is determined based on the first total energy of the microseismic events in each plane sub-interval of the coal seam roof plane interval; and the plane sub-interval frequency threshold of the coal seam roof plane interval is determined based on the first total frequency of the microseismic events in each plane sub-interval of the coal seam roof plane interval.

[0065] In one example, the plane subinterval energy threshold , where the average energy value of the plane subinterval is , the mean square error of the energy of the plane subinterval , is the first total energy of the microseismic event in the i-th plane subinterval.

[0066] In one example, the plane subinterval frequency threshold , where the plane subinterval frequency average , the mean square error of the energy of the plane subinterval , is the first total frequency of microseismic events in the i-th plane subinterval.

[0067] Compare the first total energy of the microseismic events within each plane subinterval with the plane subinterval energy threshold, identify plane subintervals whose first total energy is greater than the plane subinterval energy threshold, and mark these plane subintervals as first target plane subintervals. Similarly, compare the first total frequency of the microseismic events within each plane subinterval with the plane subinterval frequency threshold, identify plane subintervals whose first total frequency is greater than the plane subinterval frequency threshold, and mark these plane subintervals as second target plane subintervals.

[0068] Finally, the first target plane subinterval and the second target plane subinterval are merged to determine these subintervals as the first effective microseismic event activity interval in the coal seam roof plane interval in the three-dimensional space of the coal seam roof of the adjacent working face.

[0069] In some embodiments, it further includes:

[0070] Based on the width of the adjacent working face transport lane, the adjacent working face width, the adjacent working face return air lane width and the adjacent working face return air lane side coal pillar width, the coal seam roof plane interval in the three-dimensional space of the adjacent working face coal seam roof is determined;

[0071] The coal seam roof plane interval is divided into a plurality of plane sub-intervals according to a first preset interval length.

[0072] In this embodiment, before performing microseismic monitoring and analysis, the plane interval of the coal seam roof of the adjacent working face in the plane direction is first determined. Figure 2 As shown in the figure, the plane interval of the coal seam roof is determined based on the following parameters: the width of the coal pillar on the return air lane side of the adjacent working face L1 (the width of the coal pillar on the return air lane side of the adjacent working face), the width of the return air lane of the adjacent working face L2 (the width of the return air lane of the adjacent working face), the width of the adjacent working face L3 (the width of the adjacent working face) and the width of the transport lane of the adjacent working face L4 (the width of the transport lane of the adjacent working face).

[0073] Through these parameters, the plane interval L=L1+L2+L3+L4 of the coal seam roof of the adjacent working face can be determined. The coal seam roof plane interval L is then further divided into multiple smaller plane sub-intervals to facilitate a more detailed analysis of microseismic events. Specifically, it is further divided into multiple smaller plane sub-intervals according to a preset interval length (such as 10m). For example, the first plane sub-interval starts from the coal pillar on the return air lane side of the adjacent working face, extends horizontally for 10 meters in the direction of the transport lane of the adjacent working face, and then begins to divide the second plane sub-interval... and so on. Before dividing to the end of the transport lane of the adjacent working face, if the remaining plane width is less than 10 meters, then this area less than 10 meters will be discarded and will not be treated as an independent plane sub-interval.

[0074] The method for determining roof decompression provided by an embodiment of the present invention divides the plane interval of the coal seam roof in the three-dimensional space of the coal seam roof of the adjacent working face into multiple plane sub-intervals, and then accurately identifies the first effective microseismic event activity interval in the plane direction in the three-dimensional space of the coal seam roof of the adjacent working face based on the first total energy and the first total frequency of the microseismic events in each plane sub-interval.

[0075] In some embodiments, determining a second effective microseismic event activity interval of a vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face includes:

[0076] Determine the second total energy and the second total frequency of the microseismic events in each vertical sub-interval divided in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face;

[0077] determining a vertical sub-interval energy threshold based on the second total energy of the microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof;

[0078] determining a vertical sub-interval frequency threshold based on the second total frequency of microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof;

[0079] Determine a first target vertical subinterval in which the second total energy is greater than the vertical subinterval energy threshold and a second target vertical subinterval in which the second total frequency is greater than the vertical subinterval frequency threshold;

[0080] The first target vertical sub-interval and the second target vertical sub-interval are determined as second effective microseismic event activity intervals in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof of the adjacent working face.

[0081] In the three-dimensional space of the coal seam roof adjacent to the working face, the vertical interval of the coal seam roof is first divided into several vertical subintervals. For each vertical subinterval, the second total energy and second total frequency of all microseismic events occurring in the vertical subinterval are calculated.

[0082] Afterwards, the vertical sub-interval energy threshold is determined based on the second total energy of the microseismic events in each vertical sub-interval of the vertical interval of the coal seam roof; and the vertical sub-interval frequency threshold is determined based on the second total frequency of the microseismic events in each vertical sub-interval of the vertical interval of the coal seam roof.

[0083] In one example, the vertical subinterval energy threshold , where the average energy value of the vertical sub-interval is , vertical sub-interval energy mean square error , is the second total energy of the microseismic event in the i-th vertical subinterval.

[0084] In one example, the vertical subinterval frequency threshold , where the vertical sub-interval frequency average , vertical sub-interval energy mean square error , is the second total frequency of microseismic events in the i-th vertical subinterval.

[0085] Compare the first total energy of the microseismic events within each vertical subinterval with the vertical subinterval energy threshold, identify vertical subintervals whose second total energy is greater than the vertical subinterval energy threshold, and mark these vertical subintervals as first target vertical subintervals. Similarly, compare the second total frequency of the microseismic events within each vertical subinterval with the vertical subinterval frequency threshold, identify vertical subintervals whose second total frequency is greater than the vertical subinterval frequency threshold, and mark these vertical subintervals as second target vertical subintervals.

[0086] Finally, the first target vertical subinterval and the second target vertical subinterval are merged to determine these subintervals as the second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof of the adjacent working face.

[0087] In some embodiments, it further includes:

[0088] Based on the height of the coal seam roof fracture zone, determine the vertical interval of the coal seam roof in the three-dimensional space of the adjacent working face;

[0089] The vertical interval of the coal seam roof is divided into a plurality of vertical sub-intervals according to a second preset interval length.

[0090] Similarly, before conducting microseismic monitoring and analysis, the analysis interval of the coal seam roof in the vertical direction is determined first.

[0091] In this embodiment, the height of the coal seam roof fracture zone is first determined. This height serves as the basis for determining the vertical analysis interval for the adjacent working face. The coal seam roof fracture zone refers to the area in the roof rock strata where cracks and fractures occur due to mining activities. Therefore, the height of the coal seam roof fracture zone has a significant impact on roof stability. Specifically, the vertical interval determined typically begins at the coal seam roof and extends to the height of the fracture zone. This interval is the primary range for analyzing microseismic events in the vertical direction of the roof.

[0092] After determining the vertical interval of the coal seam roof adjacent to the working face, it is further divided into multiple smaller vertical subintervals to facilitate more detailed analysis of microseismic events. Specifically, it is further divided into multiple smaller vertical subintervals according to a preset interval length (such as 5 meters). For example, the first vertical subinterval starts from the coal seam roof and extends upward for 5 meters, then the second vertical subinterval is divided... and so on. Before dividing to the top of the fracture zone, if the remaining rock layer height is less than 5 meters, this area less than 5 meters will be discarded and not treated as a separate vertical subinterval.

[0093] The method for determining roof decompression provided in an embodiment of the present invention divides the vertical area within the three-dimensional space of the coal seam roof into multiple vertical sub-intervals, and then accurately identifies the second effective microseismic event activity interval in the vertical direction within the three-dimensional space of the coal seam roof based on the second total energy and second total frequency of the microseismic events in each vertical sub-interval.

[0094] The roof pressure relief determination method and device provided by the present invention are described below. The roof pressure relief determination method and device described below and the roof pressure relief determination method described above can be referenced to each other.

[0095] refer to Figure 5 In this embodiment, the top plate pressure relief determination device includes:

[0096] The first roof decompression determination module 510 is configured to determine, based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof within the three-dimensional space of the coal seam roof adjacent to the working face;

[0097] The second roof pressure relief determination module 520 is configured to determine the roof pressure relief spatial position of the current working face based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face.

[0098] The roof pressure relief determination device provided by an embodiment of the present invention analyzes the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof of the adjacent working face, accurately defines the effective microseismic event activity intervals in the planar direction and vertical direction of the adjacent working face, and then determines the spatial position of the roof pressure relief of the current working face based on these effective microseismic event activity intervals of the adjacent working face, thereby improving the accuracy of the roof pressure relief position determination.

[0099] Figure 6 An example of a physical structure diagram of an electronic device is shown below. Figure 6 As shown, the electronic device may include: a processor 610, a communications interface 620, a memory 630, and a communication bus 640, wherein the processor 610, the communications interface 620, and the memory 630 communicate with each other via the communication bus 640. The processor 610 may call logic instructions in the memory 630 to execute a roof depressurization determination method, which includes:

[0100] Based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face are determined;

[0101] Based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face, the roof decompression space position of the current working face is determined.

[0102] Furthermore, the logic instructions in the aforementioned memory 630 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0103] In another aspect, the present invention further provides a computer program product, comprising a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can perform the roof decompression determination method provided by the above methods, the method comprising:

[0104] Based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face are determined;

[0105] Based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face, the roof decompression space position of the current working face is determined.

[0106] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the method for determining roof depressurization provided by the above methods is implemented, the method comprising:

[0107] Based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face are determined;

[0108] Based on the first effective microseismic event activity interval and the second effective microseismic event activity interval of the adjacent working face, the roof decompression space position of the current working face is determined.

[0109] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of the present embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0110] Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented using software plus a necessary general-purpose hardware platform, or of course, hardware. Based on this understanding, the essence of the above technical solution, or the portion that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, or an optical disk, and includes a number of instructions for causing a computer device (such as a personal computer, server, or network device) to execute the methods described in each embodiment or certain portions of the embodiments.

[0111] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for determining roof pressure relief, characterized in that: include: Based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face are determined; Determine a second interval in a plane interval of a coal seam roof within a three-dimensional space of a coal seam roof adjacent to a working face, excluding the first effective microseismic event activity interval; Based on the second interval of the adjacent working surface, determining a roof pressure relief area in a plane direction of the working surface; Based on the second effective microseismic event activity interval of the adjacent working face, the vertical roof decompression area of ​​the working face is determined.

2. The method for determining roof pressure relief according to claim 1, wherein: The method of determining a first effective microseismic event activity interval of a coal seam roof plane interval in the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face comprises: Determine a first total energy and a first total frequency of microseismic events in each plane subinterval divided in a plane interval of a coal seam roof in a three-dimensional space of a coal seam roof adjacent to a working face; Determining a plane subinterval energy threshold based on the first total energy of the microseismic events in each plane subinterval in the coal seam roof plane interval; Determining a plane subinterval frequency threshold based on the first total frequency of microseismic events in each plane subinterval in the coal seam roof plane interval; Determine a first target plane subinterval in which the first total energy is greater than the plane subinterval energy threshold and a second target plane subinterval in which the first total frequency is greater than the plane subinterval frequency threshold; The first target plane sub-interval and the second target plane sub-interval are determined as first effective microseismic event activity intervals in the coal seam roof plane interval in the three-dimensional space of the coal seam roof adjacent to the working face.

3. The method for determining roof pressure relief according to claim 2, wherein: The method further comprises: Based on the width of the adjacent working face transport lane, the adjacent working face width, the adjacent working face return air lane width and the adjacent working face return air lane side coal pillar width, the coal seam roof plane interval in the three-dimensional space of the adjacent working face coal seam roof is determined; The coal seam roof plane interval is divided into a plurality of plane sub-intervals according to a first preset interval length.

4. The method for determining roof pressure relief according to claim 1, wherein: The method of determining a second effective microseismic event activity interval in a vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face comprises: Determine the second total energy and the second total frequency of the microseismic events in each vertical sub-interval divided in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof adjacent to the working face; determining a vertical sub-interval energy threshold based on the second total energy of the microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof; determining a vertical sub-interval frequency threshold based on the second total frequency of microseismic events in each vertical sub-interval in the vertical interval of the coal seam roof; Determine a first target vertical subinterval in which the second total energy is greater than the vertical subinterval energy threshold and a second target vertical subinterval in which the second total frequency is greater than the vertical subinterval frequency threshold; The first target vertical sub-interval and the second target vertical sub-interval are determined as second effective microseismic event activity intervals in the vertical interval of the coal seam roof in the three-dimensional space of the coal seam roof of the adjacent working face.

5. The method for determining roof pressure relief according to claim 4, wherein: The method further comprises: Based on the height of the coal seam roof fracture zone, the vertical interval of the coal seam roof in the three-dimensional space of the adjacent working face is determined; The vertical interval of the coal seam roof is divided into a plurality of vertical sub-intervals according to a second preset interval length.

6. A device for determining roof pressure relief, characterized in that: include: A first roof decompression determination module is configured to determine, based on the energy and frequency of each microseismic event in the three-dimensional space of the coal seam roof adjacent to the working face, a first effective microseismic event activity interval in the plane interval of the coal seam roof and a second effective microseismic event activity interval in the vertical interval of the coal seam roof within the three-dimensional space of the coal seam roof adjacent to the working face; The second roof decompression determination module is used to determine the second interval in the plane interval of the coal seam roof within the three-dimensional space of the coal seam roof of the adjacent working face, except the first effective microseismic event activity interval; based on the second interval of the adjacent working face, determine the plane direction roof decompression area of ​​the current working face; based on the second effective microseismic event activity interval of the adjacent working face, determine the vertical direction roof decompression area of ​​the current working face.

7. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for determining roof decompression according to any one of claims 1 to 5 is implemented.

Images