Method and device for determining stoping speed of coal seam mining in impact danger

Abstract

The invention provides an impact danger coal seam mining speed determination method, which comprises the following steps of: constructing a three-dimensional geometric model reflecting a coal mine geological structure and a coal seam form according to geological data and engineering parameter data of a coal seam, and performing grid division to generate a discretization model for numerical calculation; the discretization numerical model sets boundary conditions and a loading mode of the discretization model according to set solving parameters in combination with mechanical properties of coal rock in geological data and a coal seam environment, and simulates mechanical response data of the coal rock; according to the mechanical response data, calculating the maximum elastic energy accumulated in the coal rock unit body when the coal seam has rock burst; according to the maximum elastic energy, the safe stoping speed of coal seam mining is calculated. Therefore, by preventing and controlling the occurrence of rock burst disasters and calculating the safe recovery speed, the elastic strain energy accumulated in the coal rock unit body can be reduced, and the occurrence probability of rock burst is reduced, so that the safe production of a coal mine is guaranteed.

Description

Technical Field

[0001] The present invention relates to the technical field of coal seam mining, and in particular to a method and device for determining the recovery speed of impact-hazardous coal seam mining. Background Art

[0002] Rock burst is a serious natural disaster encountered during coal mining. It typically occurs when a large amount of elastic energy accumulates within the coal and rock mass. When this elastic energy is suddenly released, it can cause severe damage to the working face. Mining rate is a key factor influencing rock burst. As mining rate increases, the peak maximum principal stress of the overburden increases exponentially, while the release of elastic energy in the mine fluctuates significantly. This increases the risk of rock burst and compromises safe coal mine production. Therefore, properly determining the mining rate is crucial for controlling rock burst. Summary of the Invention

[0003] The present invention aims to solve one of the technical problems in the related art at least to a certain extent.

[0004] To this end, the first purpose of the present invention is to propose a method for determining the mining speed of impact-hazardous coal seams. By preventing and controlling the occurrence of impact ground pressure disasters and calculating the safe mining speed, the elastic strain energy accumulated in the coal rock unit body can be reduced, and the probability of impact ground pressure occurrence can be reduced, thereby ensuring the safe production of coal mines.

[0005] The second purpose of the present invention is to provide a device for determining the recovery speed of impact-hazardous coal seam mining.

[0006] A third object of the present invention is to provide an electronic device.

[0007] A fourth object of the present invention is to provide a non-transitory computer-readable storage medium storing computer instructions.

[0008] To achieve the above-mentioned purpose, a first embodiment of the present invention provides a method for determining the recovery speed of coal seam mining with high impact risk, the method comprising:

[0009] Collect geological data and engineering parameter data of coal seams;

[0010] Based on the geological data and engineering parameter data, a three-dimensional geometric model reflecting the geological structure of the coal mine and the coal seam morphology is constructed by simulation calculation software, and the three-dimensional geometric model is meshed to generate a discretized model for numerical calculation;

[0011] Based on the mechanical properties of coal rock in the geological data of the coal seam, the mechanical parameters of the coal rock in the discretization model are defined, and the boundary conditions and loading methods of the discretization model are set according to the coal seam environment;

[0012] The discretized numerical model simulates the mechanical response data of the coal rock according to the set solution parameters and in combination with the mechanical parameters, boundary conditions and loading methods;

[0013] Calculating the maximum elastic energy accumulated in the coal rock unit body when the coal seam experiences rock burst based on the mechanical response data of the coal rock;

[0014] The safe recovery speed of coal seam mining is calculated based on the maximum elastic energy accumulated in the coal rock unit.

[0015] To achieve the above-mentioned purpose, a second embodiment of the present invention provides a device for determining the recovery speed of a coal seam with a high impact risk, the device comprising:

[0016] Acquisition module, used to collect geological data and engineering parameter data of coal seams;

[0017] A construction module is used to construct a three-dimensional geometric model reflecting the geological structure and coal seam morphology of the coal mine through simulation calculation software based on the geological data and engineering parameter data, and to mesh the three-dimensional geometric model to generate a discretized model for numerical calculation;

[0018] The definition module is used to define the mechanical parameters of the coal rock in the discretization model according to the mechanical properties of the coal rock in the geological data of the coal seam, and to set the boundary conditions and loading methods of the discretization model according to the coal seam environment;

[0019] A simulation module is used for simulating the mechanical response data of coal and rock according to the set solution parameters of the discretized numerical model and in combination with the mechanical parameters, boundary conditions and loading methods;

[0020] a calculation module for calculating the maximum elastic energy accumulated in the coal rock unit body when the coal seam experiences rock burst according to the mechanical response data of the coal rock;

[0021] The calculation module is used to calculate the safe recovery speed of coal seam mining based on the maximum elastic energy accumulated in the coal rock unit.

[0022] To achieve the above-mentioned purpose, the third aspect embodiment of the present invention proposes an electronic device, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method described in the first aspect.

[0023] In order to achieve the above-mentioned objectives, an embodiment of the fourth aspect of the present invention proposes a non-transitory computer-readable storage medium storing computer instructions, where the computer instructions are used to enable the computer to execute the method described in the first aspect.

[0024] The method, device, electronic device, and storage medium for determining the recovery rate of coal seam mining with rock burst hazards in the embodiments of the present invention construct a three-dimensional geometric model reflecting the geological structure and morphology of the coal mine based on the geological data and engineering parameter data of the coal seam, perform grid division, and generate a discretized model for numerical calculation. The discretized numerical model simulates the mechanical response data of the coal rock based on the set solution parameters and the mechanical properties of the coal rock in the geological data and the coal seam environment. Based on the mechanical response data, the maximum elastic energy accumulated in the coal rock unit body when rock burst occurs in the coal seam is calculated. Based on the maximum elastic energy, the safe recovery rate of coal seam mining is calculated. Therefore, by preventing and controlling the occurrence of rock burst disasters and calculating the safe recovery rate, the elastic strain energy accumulated in the coal rock unit body can be reduced, the probability of rock burst occurrence can be reduced, and the safe production of coal mines can be ensured.

[0025] Additional aspects and advantages of the present invention will be set forth in part in the description which follows and, in part, will be obvious from the description which follows, or may be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0027] Figure 1 A flow chart of a method for determining a recovery speed in impact-hazardous coal seam mining provided by an embodiment of the present invention;

[0028] Figure 2 A diagram showing the relationship between the working face advance tunnel and stress mapping provided by an embodiment of the present invention;

[0029] Figure 3 A schematic structural diagram of a device for determining the recovery speed of impact-hazardous coal seam mining provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0030] The following describes embodiments of the present invention in detail, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present invention, and are not to be construed as limiting the present invention.

[0031] It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.) and signals involved in the present invention are all authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the relevant laws, regulations and standards of relevant countries and regions.

[0032] The acquisition, transmission, storage, use, and processing of data in the technical solution of the present invention comply with the relevant provisions of relevant laws and regulations.

[0033] It should be noted that in the embodiments of the present disclosure, certain software, components, models, and other existing solutions in the industry may be mentioned. They should be regarded as exemplary and their purpose is only to illustrate the feasibility of implementing the technical solution of the present invention, but it does not mean that the inventor has or will necessarily use the solution.

[0034] The following describes a method and apparatus for determining the recovery speed of impact-hazardous coal seam mining according to an embodiment of the present invention with reference to the accompanying drawings.

[0035] Figure 1 A flow chart of a method for determining the recovery speed of impact-hazardous coal seam mining provided by an embodiment of the present invention.

[0036] like Figure 1 As shown, the method includes the following steps:

[0037] Step 101: Collect geological data and engineering parameter data of the coal seam.

[0038] In some possible implementations, geological data include mechanical properties of coal rock, geological structure, coal seam thickness and inclination; engineering parameter data include coal seam mining method, support method, and mining depth.

[0039] Step 102: Based on the geological data and engineering parameter data, a three-dimensional geometric model reflecting the geological structure of the coal mine and the coal seam morphology is constructed by simulation calculation software, and the three-dimensional geometric model is meshed to generate a discretized model for numerical calculation.

[0040] In some possible implementations, when the simulation calculation software is the three-dimensional numerical analysis software FLAC3D, a three-dimensional geometric model (numerical model) reflecting the geological structure of the coal mine and the morphology of the coal seam is constructed based on geological data and engineering parameter data using three methods: writing command streams, FLAC3D modeling tools, and external software.

[0041] Specifically, the model is constructed and grouped by writing command streams; it is constructed using the modeling tools provided by FLAC3D; the model is constructed using external software to obtain a three-dimensional geometric model that reflects the geological structure of the coal mine and the morphology of the coal seam, and meshing is performed, which is then imported into FLAC3D for grouping operations to generate a discretized model for numerical calculation.

[0042] Among them, common external (modeling) software includes but is not limited to Rhino (RHINO), multi-physics field simulation components (ANSYS), and finite element analysis software (ABAQUS).

[0043] Optionally, when meshing a three-dimensional geometric model, the computational accuracy and computational efficiency of the model should be balanced to achieve dual optimization of accuracy and efficiency.

[0044] Step 103 , based on the mechanical properties of the coal rock in the geological data of the coal seam, the mechanical parameters of the coal rock in the discretization model are defined, and the boundary conditions and loading mode of the discretization model are set according to the coal seam environment.

[0045] Optionally, the mechanical properties of the coal rock include the lithology, thickness, bulk density, unidirectional compressive strength, unidirectional tensile strength, elastic modulus, Poisson's ratio, density, cohesion, and internal friction angle of the coal rock.

[0046] In some possible implementations, setting the boundary conditions of the discretized model according to the coal seam environment includes: applying appropriate boundary conditions and loading methods to the discretized model according to the actual coal seam environment conditions at the coal mine site. For example, part of the boundary of the discretized model can be fixed to simulate constraint conditions, or a gravity field can be applied to simulate the influence of the earth's gravity.

[0047] Specifically, the boundary conditions of the discretization model are as follows: constraints along the X-axis are imposed on the boundaries at both ends of the X-axis of the discretization model, that is, the displacement of the boundary in the X direction is zero; constraints along the Y-axis are imposed on the boundaries at both ends of the Y-axis of the discretization model, that is, the displacement of the boundary in the Y direction is zero; the bottom boundary of the discretization model is fixed, that is, the displacement of the bottom boundary in the X, Y, and Z directions are all zero, and the top of the discretization model is a free boundary, where the X, Y, and Z axes represent the horizontal direction, vertical direction, and depth or telescopic direction of the discretization model, respectively.

[0048] In step 104 , the discretized numerical model simulates the mechanical response data of the coal rock according to the set solution parameters and in combination with the mechanical parameters, boundary conditions and loading methods.

[0049] In some possible implementations, the solution parameters include the solution time, the number of iteration steps, and the convergence criterion of the discretized model.

[0050] Optionally, the FLAC3D solver can be used to perform numerical calculations (solving the discretized numerical model) to obtain the mechanical response data (stress, strain, and displacement change data) of the discretized numerical model under given conditions (solving parameters, and combining mechanical parameters, boundary conditions, and loading methods).

[0051] Step 105: Calculate the maximum elastic energy accumulated in the coal rock unit body when rock burst occurs in the coal seam based on the mechanical response data of the coal rock.

[0052] In some possible implementations, the maximum elastic energy accumulated in a coal rock unit when a rock burst occurs in a coal seam is calculated based on the mechanical response data of the coal rock, including: the mechanical response data of the coal rock includes stress, strain, and displacement change data of the coal rock; when the mechanical response data of the coal rock is stress data, the maximum elastic energy U accumulated in a coal rock unit when a rock burst occurs in a coal seam is calculated based on the three-dimensional stress of the coal rock simulated by a discretized model within each grid cell in a three-dimensional geometric model. Among them, σ i is the three-dimensional stress in the X, Y, and Z directions, i = 1, 2, 3, E is the elastic modulus, and μ is the Poisson's ratio.

[0053] Step 106, calculating the safe recovery speed of coal seam mining based on the maximum elastic energy accumulated in the coal rock unit.

[0054] In some possible implementations, rock bursts require certain elastic energy conditions to occur. For example, rock bursts often occur within 50 m ahead of the return air tunnel, and single pillars (coal rock units) often break 0 to 5 m ahead of the working face. Statistical analysis shows that the unit elastic energy of microseismic events is often on the order of 10 5 to 10 7 J. Using the embedded programming language (fish) and visualization tool (TEC360 software) built into FLAC3D software, the three-dimensional stress of each grid cell in the three-dimensional geometric model was extracted and the maximum elastic energy U of the coal rock unit was calculated. U was then imported into the TEC360 software to draw a stress distribution cloud map 5 meters ahead of the working face advance tunnel. When the working face advances at a speed of 2.4 m / d, the maximum elastic energy accumulated in the coal rock unit is 9.10×10 4 J, compared with previous experience, the elastic energy conditions for rock burst are not met; when the working face advances at a speed of 3.6 m / d, the maximum elastic energy of the coal-rock unit is 1.07×10 5J, slightly exceeding the elastic energy requirement for rock burst, is considered a critical state. As the working face advances further, the elastic energy accumulated within the coal and rock units reaches the threshold for rock burst, and the rate of elastic energy increase increases with increasing working face speed. Based on the above stress and energy analysis, the safe advance speed of the advanced working face (the safe recovery speed for coal seam mining) should be less than 3.6 m / d.

[0055] The method for determining the recovery rate of coal seam mining with rock burst hazards in an embodiment of the present invention constructs a three-dimensional geometric model reflecting the geological structure and coal seam morphology of the coal mine based on the geological data and engineering parameter data of the coal seam, performs grid division, and generates a discretized model for numerical calculation. The discretized numerical model simulates the mechanical response data of the coal rock based on the set solution parameters and the mechanical properties of the coal rock in the geological data and the coal seam environment. Based on the mechanical response data, the maximum elastic energy accumulated in the coal rock unit body when rock burst occurs in the coal seam is calculated. Based on the maximum elastic energy, the safe recovery rate of coal seam mining is calculated. Therefore, by preventing and controlling the occurrence of rock burst disasters and calculating the safe recovery rate, the elastic strain energy accumulated in the coal rock unit body can be reduced, the probability of rock burst occurrence can be reduced, and the safe production of coal mines can be ensured.

[0056] In summary, a discretized model (numerical model) with dimensions of 600 m (length × width × height) × 496 m (height × height) × 110.6 m was established. The distribution and mechanical parameters of the coal seam's roof and floor strata are shown in Table 1. This discretized model uses the Moore-Coulomb model. To ensure simulation results are consistent with actual conditions, the mechanical parameters of each coal rock in the coal seam were determined based on laboratory measurements of coal samples collected from the field and previous engineering experience.

[0057] The boundary conditions for the numerical model were set as follows: constraints were applied along the X-axis at both ends of the X-axis, meaning that the boundary displacement in the X direction was zero; constraints were applied along the Y-axis at both ends of the Y-axis, meaning that the boundary displacement in the Y direction was zero; the bottom boundary was fixed, meaning that the displacements in the X, Y, and Z directions were all zero; and the top boundary was a free boundary. The initial stress conditions were: a gradient stress of 8.18 MPa to 6.35 MPa was applied in the X-axis; a gradient stress of 13.40 MPa to 10.40 MPa was applied in the Y-axis; and a gradient stress of 9.45 MPa to 7.33 MPa was applied in the Z-axis. The upper rock layer of the numerical model was simulated using an equivalent load, with a deadweight load of 6.23 MPa applied in the Z-axis. The coal mining drum cutting depth was 0.6m, and four schemes of daily working face advancement speed (coal seam mining recovery speed) were simulated at 2.4m / d, 3.6m / d, 4.8m / d, and 6.0m / d. The distribution of coal seam top and bottom strata (rock type, thickness m, bulk density (kN·m-3 ), uniaxial compressive strength MPa, uniaxial tensile strength MPa,) and mechanical parameters are shown in Table 1:

[0058] Table 1 Distribution and mechanical parameters of coal seam roof and floor strata

[0059]

[0060] Therefore, the stress data of coal and rock simulated by the numerical model are calculated under the working face advancement speed on each day. When the working face advancement speed is 60m, the stress data within 100m of the working face ahead roadway are plotted into a curve, as shown in Figure 2. Figure 2 When the working face advance speed is 2.4 m / d, the leading stress peak is 15.1 MPa, the peak point is about 10 m away from the coal wall, and the leading influence range (working face leading tunnel) is 40 m; when the working face advance speed is 3.2 m / d, the leading stress peak is 15.7 MPa, the peak point is about 10 m away from the coal wall, and the leading influence range is 50 m; when the working face advance speed is 4.8 m / d, the leading stress peak is 17.4 MPa, the peak point is about 8 m away from the coal wall, and the leading influence range is 68 m; when the working face advance speed is 6.0 m / d, the leading stress peak is 20.1 MPa, the peak point is about 6 m away from the coal wall, and the leading influence range is 84 m. When the working face advance speed exceeds 3.2 m / d, the leading stress peak and the leading stress influence range both increase significantly. At a working face advance speed of 6.0 m / d, the leading stress peak increases by 33% compared to 2.4 m / d, and the leading stress influence range expands by 2.1 times. This shows that there is a critical value for working face advance speed. When this critical value is exceeded, the stress state of the coal body changes significantly, and the leading stress peak and its influence range both increase significantly, which is consistent with the theoretical analysis results mentioned above.

[0061] In order to realize the above embodiment, the present invention also proposes a device for determining the recovery speed of impact-hazardous coal seam mining.

[0062] Figure 3 A schematic structural diagram of a device for determining the recovery speed of impact-hazardous coal seam mining provided by an embodiment of the present invention.

[0063] like Figure 3 As shown, the recovery speed determination device 30 for impact-hazardous coal seam mining includes: an acquisition module 31 , a construction module 32 , a definition module 33 , a simulation module 34 , an estimation module 35 , and a calculation module 36 .

[0064] The acquisition module 31 is used to acquire geological data and engineering parameter data of the coal seam;

[0065] A construction module 32 is used to construct a three-dimensional geometric model reflecting the geological structure and coal seam morphology of the coal mine through simulation calculation software based on the geological data and engineering parameter data, and to mesh the three-dimensional geometric model to generate a discretized model for numerical calculation;

[0066] A definition module 33 is used to define the mechanical parameters of the coal rock in the discretization model according to the mechanical properties of the coal rock in the geological data of the coal seam, and to set the boundary conditions and loading mode of the discretization model according to the coal seam environment;

[0067] The simulation module 34 is used for simulating the mechanical response data of the coal rock according to the set solution parameters of the discretized numerical model and in combination with the mechanical parameters, boundary conditions and loading mode;

[0068] The calculation module 35 is used to calculate the maximum elastic energy accumulated in the coal rock unit body when the coal seam produces rock burst according to the mechanical response data of the coal rock;

[0069] The calculation module 36 is used to calculate the safe recovery speed of coal seam mining based on the maximum elastic energy accumulated in the coal rock unit.

[0070] Furthermore, in a possible implementation of the embodiment of the present invention, wherein:

[0071] The geological data include mechanical properties of coal rock, geological structure, coal seam thickness and dip;

[0072] The mechanical properties of the coal rock include the lithology, thickness, bulk density, uniaxial compressive strength, uniaxial tensile strength, elastic modulus, Poisson's ratio, density, cohesion, and internal friction angle of the coal rock.

[0073] The solution parameters include the solution time, number of iteration steps, and convergence criterion of the discretized model.

[0074] Furthermore, in a possible implementation of an embodiment of the present invention, when the simulation calculation software is the three-dimensional numerical analysis software FLAC3D, a three-dimensional geometric model reflecting the geological structure of the coal mine and the morphology of the coal seam is constructed according to the geological data and engineering parameter data using three methods: writing command streams, FLAC3D modeling tools, and external software.

[0075] Furthermore, in a possible implementation of an embodiment of the present invention, the boundary conditions of the discretization model are as follows: constraints along the X-axis are imposed on the boundaries at both ends of the X-axis of the discretization model, that is, the displacement of the boundary in the X direction is zero; constraints along the Y-axis are imposed on the boundaries at both ends of the Y-axis of the discretization model, that is, the displacement of the boundary in the Y direction is zero; the bottom boundary of the discretization model is fixed, that is, the displacement of the bottom boundary in the X, Y, and Z directions are all zero, and the top of the discretization model is a free boundary, wherein the X, Y, and Z axes respectively represent the horizontal direction, vertical direction, and depth or telescopic direction of the discretization model.

[0076] Furthermore, in a possible implementation of the embodiment of the present invention, the calculation module 35 is specifically configured to:

[0077] The mechanical response data of coal rock include stress, strain and displacement change data of coal rock. When the mechanical response data of coal rock is stress data, the maximum elastic energy U accumulated in the coal rock unit body when the coal seam produces rock burst is calculated based on the three-dimensional stress of coal rock simulated by the discretized model in each grid unit in the three-dimensional geometric model.

[0078] Among them, σ i is the three-dimensional stress in the X, Y, and Z directions, i = 1, 2, 3, E is the elastic modulus, and μ is the Poisson's ratio.

[0079] The device for determining the recovery rate for mining a coal seam with rock burst danger according to an embodiment of the present invention constructs a three-dimensional geometric model reflecting the geological structure and morphology of the coal mine based on the geological data and engineering parameter data of the coal seam, performs grid division, and generates a discretized model for numerical calculation. The discretized numerical model simulates the mechanical response data of the coal rock based on the set solution parameters and the mechanical properties of the coal rock in the geological data and the coal seam environment. Based on the mechanical response data, the maximum elastic energy accumulated in the coal rock unit body when rock burst occurs in the coal seam is calculated. Based on the maximum elastic energy, the safe recovery rate for coal seam mining is calculated. Thus, by preventing and controlling the occurrence of rock burst disasters and calculating the safe recovery rate, the elastic strain energy accumulated in the coal rock unit body can be reduced, the probability of rock burst occurrence can be reduced, and the safe production of coal mines can be ensured.

[0080] In order to implement the above embodiment, the present invention further provides an electronic device, including:

[0081] at least one processor; and

[0082] a memory communicatively connected to the at least one processor; wherein,

[0083] The memory stores instructions that can be executed by the at least one processor. The instructions are executed by the at least one processor to enable the at least one processor to perform the aforementioned method.

[0084] In order to implement the above embodiment, the present invention further proposes a non-transitory computer-readable storage medium storing computer instructions, where the computer instructions are used to enable the computer to execute the above method.

[0085] In the description of this specification, the reference terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art can combine and combine different embodiments or examples described in this specification and features of different embodiments or examples without contradiction.

[0086] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of the technical features being referred to. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.

[0087] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may be performed out of the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention pertain.

[0088] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include the following: an electrical connection with one or more wires (electronic devices), a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and a portable compact disc read-only memory (CDROM). Furthermore, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing it in another suitable manner if necessary, and then storing it in a computer memory.

[0089] It should be understood that various parts of the present invention can be implemented using hardware, software, firmware, or a combination thereof. In the above-described embodiments, multiple steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0090] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0091] In addition, the functional units in the various embodiments of the present invention may be integrated into a single processing module, or each unit may exist physically separately, or two or more units may be integrated into a single module. The aforementioned integrated modules may be implemented in the form of hardware or in the form of software functional modules. If the integrated modules are implemented in the form of software functional modules and sold or used as independent products, they may also be stored in a computer-readable storage medium.

[0092] The storage medium mentioned above may be a read-only memory, a magnetic disk, or an optical disk, etc. Although the embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and are not to be construed as limiting the present invention. Persons skilled in the art may make changes, modifications, substitutions, and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for determining the recovery speed of impact-hazardous coal seam mining, characterized in that: The method comprises: Collect geological data and engineering parameter data of coal seams; Based on the geological data and engineering parameter data, a three-dimensional geometric model reflecting the geological structure of the coal mine and the coal seam morphology is constructed by simulation calculation software, and the three-dimensional geometric model is meshed to generate a discretized model for numerical calculation; Based on the mechanical properties of coal rock in the geological data of the coal seam, the mechanical parameters of the coal rock in the discretization model are defined, and the boundary conditions and loading methods of the discretization model are set according to the coal seam environment; The discretized numerical model simulates the mechanical response data of the coal rock according to the set solution parameters and in combination with the mechanical parameters, boundary conditions and loading methods; Calculating the maximum elastic energy accumulated in the coal rock unit body when the coal seam experiences rock burst based on the mechanical response data of the coal rock; The safe recovery speed of coal seam mining is calculated based on the maximum elastic energy accumulated in the coal rock unit.

2. The method for determining the recovery speed of impact-hazardous coal seam mining according to claim 1 is characterized in that: in: The geological data include mechanical properties of coal rock, geological structure, coal seam thickness and dip; The mechanical properties of the coal rock include lithology, thickness, bulk density, uniaxial compressive strength, uniaxial tensile strength, elastic modulus, Poisson's ratio, density, cohesion, and internal friction angle of the coal rock; The solution parameters include the solution time, number of iteration steps, and convergence criterion of the discretized model.

3. The method for determining the recovery speed of impact-hazardous coal seam mining according to claim 1 is characterized in that: When the simulation calculation software is the three-dimensional numerical analysis software FLAC3D, based on the geological data and engineering parameter data, three methods of writing command streams, FLAC3D modeling tools and external software are used to construct a three-dimensional geometric model reflecting the geological structure of the coal mine and the morphology of the coal seam.

4. The method for determining the recovery speed of impact-hazardous coal seam mining according to claim 1 is characterized in that: The boundary conditions of the discretized model are as follows: constraints along the X axis are imposed on the boundaries at both ends of the X axis of the discretized model, that is, the displacement of the boundary in the X direction is zero; constraints along the Y axis are imposed on the boundaries at both ends of the Y axis of the discretized model, that is, the displacement of the boundary in the Y direction is zero; the bottom boundary of the discretized model is fixed, that is, the displacement of the bottom boundary in the X, Y, and Z directions are all zero, and the top of the discretized model is a free boundary, wherein the X, Y, and Z axes represent the horizontal direction, vertical direction, and depth or telescopic direction of the discretized model, respectively.

5. The method for determining the recovery speed of impact-hazardous coal seam mining according to claim 4 is characterized in that: The calculation of the maximum elastic energy accumulated in the coal rock unit body when the coal seam rock burst occurs based on the mechanical response data of the coal rock includes: The mechanical response data of coal rock include stress, strain and displacement change data of coal rock. When the mechanical response data of coal rock is stress data, the maximum elastic energy U accumulated in the coal rock unit body when the coal seam produces rock burst is calculated based on the three-dimensional stress of coal rock simulated by the discretized model in each grid unit in the three-dimensional geometric model. Among them, σ i is the three-dimensional stress in the X, Y, and Z directions, i = 1, 2, 3, E is the elastic modulus, and μ is the Poisson's ratio.

6. A device for determining the recovery speed of impact-hazardous coal seam mining, characterized in that: The device comprises: Acquisition module, used to collect geological data and engineering parameter data of coal seams; A construction module is used to construct a three-dimensional geometric model reflecting the geological structure and coal seam morphology of the coal mine through simulation calculation software based on the geological data and engineering parameter data, and to mesh the three-dimensional geometric model to generate a discretized model for numerical calculation; The definition module is used to define the mechanical parameters of the coal rock in the discretization model according to the mechanical properties of the coal rock in the geological data of the coal seam, and to set the boundary conditions and loading methods of the discretization model according to the coal seam environment; A simulation module is used for simulating the mechanical response data of coal and rock according to the set solution parameters of the discretized numerical model and in combination with the mechanical parameters, boundary conditions and loading methods; a calculation module for calculating the maximum elastic energy accumulated in the coal rock unit body when the coal seam experiences rock burst according to the mechanical response data of the coal rock; The calculation module is used to calculate the safe recovery speed of coal seam mining based on the maximum elastic energy accumulated in the coal rock unit.

7. The device for determining the recovery speed of impact-hazardous coal seam mining according to claim 6, characterized in that: in: The geological data include mechanical properties of coal rock, geological structure, coal seam thickness and dip; The mechanical properties of the coal rock include lithology, thickness, bulk density, uniaxial compressive strength, uniaxial tensile strength, elastic modulus, Poisson's ratio, density, cohesion, and internal friction angle of the coal rock; The solution parameters include the solution time, number of iteration steps, and convergence criterion of the discretized model.

8. The device for determining the recovery speed of impact-hazardous coal seam mining according to claim 6, characterized in that: When the simulation calculation software is the three-dimensional numerical analysis software FLAC3D, based on the geological data and engineering parameter data, three methods of writing command streams, FLAC3D modeling tools and external software are used to construct a three-dimensional geometric model reflecting the geological structure of the coal mine and the morphology of the coal seam.

9. The device for determining the recovery speed of impact-hazardous coal seam mining according to claim 6, characterized in that: The boundary conditions of the discretized model are as follows: constraints along the X axis are imposed on the boundaries at both ends of the X axis of the discretized model, that is, the displacement of the boundary in the X direction is zero; constraints along the Y axis are imposed on the boundaries at both ends of the Y axis of the discretized model, that is, the displacement of the boundary in the Y direction is zero; the bottom boundary of the discretized model is fixed, that is, the displacement of the bottom boundary in the X, Y, and Z directions are all zero, and the top of the discretized model is a free boundary, wherein the X, Y, and Z axes represent the horizontal direction, vertical direction, and depth or telescopic direction of the discretized model, respectively.

10. The device for determining the recovery speed of impact-hazardous coal seam mining according to claim 9, characterized in that: The calculation module is specifically used to: The mechanical response data of coal rock include stress, strain and displacement change data of coal rock. When the mechanical response data of coal rock is stress data, the maximum elastic energy U accumulated in the coal rock unit body when the coal seam produces rock burst is calculated based on the three-dimensional stress of coal rock simulated by the discretized model in each grid unit in the three-dimensional geometric model. Among them, σ i is the three-dimensional stress in the X, Y, and Z directions, i = 1, 2, 3, E is the elastic modulus, and μ is the Poisson's ratio.

11. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to perform the method according to any one of claims 1 to 5.

12. A non-transitory computer-readable storage medium storing computer instructions, characterized in that: The computer instructions are used to cause the computer to execute the method according to any one of claims 1 to 5.

Images