An assessment method and system for the roof stability of concealed voids in open-pit coal mines

Abstract

The present invention provides a method and a system for evaluating the stability of the roof of hidden voids in an open-pit coal mine. The method includes identifying the hidden voids in the open-pit coal mine and obtaining the void shape information of the hidden voids; dividing the roof area of the hidden voids by the grid method to obtain the stress information and displacement information of the roof rock strata of each grid; obtaining the roof evaluation coefficient of each grid based on the void shape information, stress information and displacement information; and obtaining the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the void roof. Using the method of the present invention can improve the accuracy of evaluating the stability of the roof of hidden voids.

Description

Technical Field

[0001] The present invention relates to the technical field of stability assessment, and particularly to a method and a system for assessing the roof stability of hidden goafs in open-pit coal mines. Background Art

[0002] Hidden goafs in open-pit coal mines refer to a large number of underground mined-out areas formed in the early stage when room-and-pillar mining methods and other open-stoping mining methods are used in underground coal mines. Due to incomplete or lost design data, the positions and boundaries of many of these mined-out areas cannot be accurately determined in the end, leaving a large number of unknown mined-out areas inside and outside the mining boundary of open-pit coal mines. The existence of hidden goafs poses a serious threat to the safe production of open-pit coal mines, and may lead to a series of disastrous events such as roof collapse, ground settlement, and mining equipment falling into pits, endangering the safety of personnel and equipment.

[0003] Although the existing technologies can determine the positions, shapes of goafs, stress changes at the roofs of goafs, etc., due to the large area of the roof and the irregularity of the roof caused by mining and stripping operations, the data obtained by a single monitoring device are limited and cannot comprehensively represent the stability information of the overall roof, resulting in poor accuracy of stability assessment. Summary of the Invention

[0004] The present invention aims to solve at least one of the technical problems in the related technologies to some extent.

[0005] To this end, a first object of the present invention is to propose a method for assessing the roof stability of hidden goafs in open-pit coal mines to improve the accuracy of assessing the roof stability of hidden goafs.

[0006] A second object of the present invention is to propose a system for assessing the roof stability of hidden goafs in open-pit coal mines.

[0007] A third object of the present invention is to propose an electronic device.

[0008] A fourth object of the present invention is to propose a computer-readable storage medium.

[0009] To achieve the above object, a first aspect of the present invention proposes a method for assessing the roof stability of hidden goafs in open-pit coal mines, including:

[0010] Identifying hidden goafs in an open-pit coal mine and obtaining the goaf shape information of the hidden goafs;

[0011] Dividing the roof area of the hidden goafs by a grid method to obtain the stress information and displacement information of the roof rock strata of each grid;

[0012] Obtaining a roof evaluation coefficient for each grid based on the goaf shape information, the stress information, and the displacement information;

[0013] The roof warning coefficient is obtained based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof.

[0014] In the method of the first aspect of the present invention, obtaining the roof evaluation coefficients of each grid based on the goaf shape information, the stress information, and the displacement information includes: calculating the stress dynamic risk coefficient of each grid based on the stress information; calculating the displacement fluctuation coefficient of each grid based on the displacement information; calculating the spatial characteristic coefficient of the roof area based on the goaf shape information; and obtaining the roof evaluation coefficient of each grid based on the stress dynamic risk coefficient, the displacement fluctuation coefficient, and the spatial characteristic coefficient.

[0015] In the method of the first aspect of the present invention, the roof evaluation coefficients of each grid satisfy:

[0016]

[0017] wherein, PG represents the roof evaluation coefficient of any grid, and α1, α2, and α3 are the proportionality coefficients of the stress dynamic risk coefficient, the displacement fluctuation coefficient, and the spatial characteristic coefficient of this grid respectively, YL dt is the stress dynamic risk coefficient of this grid, WY bd is the displacement fluctuation coefficient of this grid, and KJ tz is the spatial characteristic coefficient of the roof area.

[0018] In the method of the first aspect of the present invention, the stress information includes the horizontal stress and the overlying stress per unit time, and calculating the stress dynamic risk coefficient of each grid based on the stress information includes: calculating the horizontal stress deviation coefficient and the overlying stress deviation coefficient respectively based on the horizontal stress and the overlying stress per unit time of each grid, and then calculating the stress dynamic risk coefficient of the corresponding grid.

[0019] In the method of the first aspect of the present invention, the displacement information includes the initial position of the displacement sensor and the distance from the current position of the displacement sensor to the initial position per unit time. Calculating the displacement fluctuation coefficient of each grid based on the displacement information includes: obtaining the average value and the standard deviation of the corresponding distance based on the initial position of the displacement sensor and the distance of each grid, and then calculating the displacement fluctuation coefficient of the corresponding grid.

[0020] In the method of the first aspect of the present invention, the goaf shape information includes the span, height, and inclination angle of the concealed goaf. Calculating the spatial characteristic coefficient of the roof area based on the goaf shape information includes: obtaining the deviations corresponding to the span, height, and inclination angle based on the span, height, and inclination angle, and then calculating the spatial characteristic coefficient of the roof area.

[0021] In the method of the first aspect of the present invention, obtaining the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof includes: calculating to obtain the roof warning coefficient based on the roof evaluation coefficients of all grids and the roof evaluation coefficient threshold, and determining the stability of the goaf roof based on the roof warning coefficient and the roof warning coefficient threshold.

[0022] To achieve the above object, a second aspect of the present invention proposes an evaluation system for the stability of the concealed goaf roof in an open-pit coal mine, including:

[0023] A first acquisition module, configured to identify the concealed goaf in the open-pit coal mine and acquire the goaf shape information of the concealed goaf;

[0024] A second acquisition module, configured to divide the roof area of the concealed goaf by the grid method to obtain the stress information and displacement information of the roof rock stratum of each grid;

[0025] A calculation module, configured to obtain the roof evaluation coefficient of each grid based on the goaf shape information, the stress information, and the displacement information;

[0026] An evaluation module, configured to obtain the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof.

[0027] To achieve the above object, a third aspect of the present invention proposes an electronic device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer execution instructions; the processor executes the computer execution instructions stored in the memory to implement the method proposed in the first aspect of the present invention.

[0028] To achieve the above object, a fourth aspect of the present invention proposes a computer-readable storage medium, in which computer execution instructions are stored, and when the computer execution instructions are executed by a processor, they are used to implement the method proposed in the first aspect of the present invention.

[0029] The method, system, electronic device and storage medium for evaluating the roof stability of hidden goafs in open-pit coal mines provided by the present invention identify the hidden goafs in the open-pit coal mines and obtain the goaf shape information of the hidden goafs; divide the roof area of the hidden goafs by the grid method to obtain the stress information and displacement information of the roof rock strata of each grid; obtain the roof evaluation coefficient of each grid based on the goaf shape information, stress information and displacement information; and obtain the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the roof stability of the goaf. In this case, comprehensively dividing the roof area into grids and comprehensively analyzing the goaf shape information, stress information and displacement information to determine the roof evaluation coefficient of the grid, and then obtaining the roof warning coefficient to evaluate the roof stability of the goaf helps to improve the roof stability, thereby improving the accuracy of the evaluation of the roof stability of the hidden goaf.

[0030] Additional aspects and advantages of the present invention will be given in part in the following description, become apparent in part from the following description, or be learned through the practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0032] Figure 1 is a schematic flow chart of a method for evaluating the roof stability of hidden goafs in an open-pit coal mine provided by an embodiment of the present invention;

[0033] Figure 2 is a block diagram of a system for evaluating the roof stability of hidden goafs in an open-pit coal mine provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0034] Embodiments of the present invention will be described in detail below. Examples of the embodiments are shown in the accompanying drawings, in which the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to explain the present invention and should not be construed as limiting the present invention.

[0035] The method and system for evaluating the roof stability of hidden goafs in an open-pit coal mine according to embodiments of the present invention will be described below with reference to the accompanying drawings.

[0036] Embodiments of the present invention provide a method for evaluating the roof stability of hidden goafs in an open-pit coal mine to improve the accuracy of the evaluation of the roof stability of hidden goafs.

[0037] Figure 1 is a schematic flow chart of a method for evaluating the roof stability of hidden goafs in an open-pit coal mine provided by an embodiment of the present invention.

[0038] As Figure 1As shown, the method for evaluating the stability of the roof of the hidden goaf in the open-pit coal mine includes the following steps:

[0039] Step S101: Identify the hidden goaf in the open-pit coal mine and obtain the goaf shape information of the hidden goaf.

[0040] In step S101, technologies such as geophysical prospecting, three-dimensional laser scanning, and sonar detection are used to identify the hidden goaf in the open-pit coal mine.

[0041] In step S101, the obtained goaf shape information includes the span, height, and inclination angle of the hidden goaf.

[0042] It is easy to understand that the span refers to the maximum width of the hidden goaf (subsequently abbreviated as the goaf) in the horizontal direction. The inclination angle refers to the inclination degree of the goaf relative to the horizontal plane, usually expressed in degrees. The height refers to the distance from the bottom plate of the goaf to the roof of the goaf.

[0043] Step S102: Divide the roof area of the hidden goaf by the grid method to obtain the stress information and displacement information of the roof rock strata of each grid.

[0044] In step S102, the steps of dividing the roof area of the hidden goaf by the grid method specifically include:

[0045] 1) First, perform regional division. Divide the roof area of the hidden goaf into multiple small areas according to a certain scale to form a grid-like division structure; the division scale of the small areas can be determined according to specific circumstances, usually depending on geological conditions, rock stratum characteristics, and the required analysis accuracy.

[0046] 2) Then, establish grids. Establish grids in each small area of the roof area of the hidden goaf, usually using a regular rectangular or orthogonal grid structure. The grids in each small area can be numbered through a coordinate system or other identifiers for subsequent analysis.

[0047] In step S102, after dividing the roof area of the hidden goaf by the grid method, monitoring devices are reasonably distributed into each grid according to the divided grids. The monitoring devices include stress sensors and displacement sensors.

[0048] In step S102, the obtained stress information of the roof rock strata includes the horizontal stress and overlying stress per unit time at the roof rock strata. The stress information of the roof rock strata is collected by stress sensors.

[0049] In step S102, the obtained displacement information of the roof rock strata includes the initial position of the displacement sensor and the distance from the current position of the displacement sensor to the initial position per unit time. The distance from the current position to the initial position is collected by the displacement sensor.

[0050] Step S103: Obtain the roof evaluation coefficients of each grid based on the goaf shape information, stress information, and displacement information.

[0051] In step S103, since the goaf shape information, stress information, and displacement information obtained in steps S101 and S102 are all the main factors affecting the roof stability, in this step, the roof evaluation coefficients of each grid are obtained by comprehensively considering the goaf shape information, stress information, and displacement information to evaluate the stability changes of different grids of the roof, and further determine the potential risks of different grids of the roof.

[0052] The obtaining of the roof evaluation coefficients of each grid based on the goaf shape information, stress information, and displacement information in step S103 includes: calculating the stress dynamic risk coefficient of each grid based on the stress information; calculating the displacement fluctuation coefficient of each grid based on the displacement information; calculating the spatial characteristic coefficient of the roof area based on the goaf shape information; obtaining the roof evaluation coefficients of each grid based on the stress dynamic risk coefficient, displacement fluctuation coefficient, and spatial characteristic coefficient.

[0053] Among them, calculating the stress dynamic risk coefficient of each grid based on the stress information includes: calculating the horizontal stress deviation coefficient and the overlying stress deviation coefficient respectively based on the horizontal stress and the overlying stress per unit time of each grid, and then calculating the stress dynamic risk coefficient of the corresponding grid. Calculating the displacement fluctuation coefficient of each grid based on the displacement information includes: obtaining the average value and standard deviation of the corresponding distance based on the initial position and distance of the displacement sensor of each grid, and then calculating the displacement fluctuation coefficient of the corresponding grid. Calculating the spatial characteristic coefficient of the roof area based on the goaf shape information includes: obtaining the deviations corresponding to the span, height, and inclination angle based on the span, height, and inclination angle, and then calculating the spatial characteristic coefficient of the roof area.

[0054] Specifically, for the stress dynamic risk coefficient, the obtaining process of the stress dynamic risk coefficient of each grid includes:

[0055] Select any grid, determine the horizontal stress and the overlying stress per unit time according to the stress sensor in this grid, and represent the horizontal stress and the overlying stress per unit time by SP n and SF n respectively, where n = 1, 2, 3,..., N, N represents the number of horizontal stress and overlying stress collected per unit time, N is a positive integer, and n is the number of horizontal stress and overlying stress collected per unit time;

[0056] Obtain the maximum horizontal stress threshold and the maximum overlying stress threshold, and represent the maximum horizontal stress threshold and the maximum overlying stress threshold by SP yz and SF yzIt is shown that the horizontal stress and overburden stress greater than the maximum horizontal stress threshold and the maximum overburden stress threshold respectively within a unit time are obtained, and the horizontal stress greater than the maximum horizontal stress threshold and the overburden stress greater than the maximum overburden stress threshold are represented by SP m and SF i It is shown that, where m = 1, 2, 3, ……, M, M represents the number of horizontal stresses greater than the maximum horizontal stress threshold, M is a positive integer, m is the number of the horizontal stress greater than the maximum horizontal stress threshold, i = 1, 2, 3, ……, I, I represents the number of overburden stresses greater than the maximum overburden stress threshold, I is a positive integer, and i is the number of the overburden stress greater than the maximum overburden stress threshold;

[0057] Calculate the horizontal stress deviation coefficient, and the calculation formula is: where PC sp is the horizontal stress deviation coefficient of this grid;

[0058] Calculate the overburden stress deviation coefficient, and the calculation formula is: where PC sf is the overburden stress deviation coefficient of this grid;

[0059] Calculate the stress dynamic risk coefficient of this grid, and the calculation formula is: where YL dt is the stress dynamic risk coefficient of this grid.

[0060] For the stress dynamic risk coefficients of the remaining grids, refer to the above calculation process in turn, so as to obtain the stress dynamic risk coefficients of each grid.

[0061] Regarding the displacement fluctuation coefficient, the acquisition process of the displacement fluctuation coefficient of each grid includes:

[0062] Arbitrarily select a grid, set the initial positions of the corresponding displacement sensors according to the positions where the displacement sensors are installed in this grid, obtain the distance from the current position of the displacement sensor to the initial position of the displacement sensor within a unit time, and represent the distance from the current position of the displacement sensor to the initial position of the displacement sensor within a unit time by JL n ;

[0063] Obtain the average value and standard deviation of the distance from the current position of the displacement sensor to the initial position of the displacement sensor within a unit time, and represent the average value and standard deviation of the distance from the current position of the displacement sensor to the initial position of the displacement sensor within a unit time by JL avg and JL bzc respectively, where N1 represents the number of displacement sensors, N1 is a positive integer, and n1 represents the number of the displacement sensor;

[0064] Calculate the displacement fluctuation coefficient, and the calculation formula is as follows: where, WY bd is the displacement fluctuation coefficient of this grid.

[0065] For the displacement fluctuation coefficients of the remaining grids, refer to the above calculation process in sequence, so as to obtain the displacement fluctuation coefficients of each grid.

[0066] Regarding the spatial characteristic coefficient, the acquisition process of the spatial characteristic coefficient in the roof area includes:

[0067] Collect the spatial characteristics of the hidden goaf in the open-pit coal mine. The spatial characteristics of the hidden goaf in the open-pit coal mine include the span, height and inclination angle of the goaf. Denote the span of the goaf as KD, the height of the goaf as GD, and the inclination angle of the goaf as JD;

[0068] Obtain the standard thresholds of the span, height and inclination angle of the goaf. By comparing the span, height and inclination angle of the hidden goaf in the open-pit coal mine with the standard thresholds of the span, height and inclination angle of the goaf, obtain the deviations of the span, height and inclination angle of the goaf, and denote the span deviation of the goaf as KD pc denote, the height deviation of the goaf as GD pc denote, and the inclination angle deviation of the goaf as JD pc denote;

[0069] Calculate the spatial characteristic coefficient, and the calculation formula is as follows: where, KJ tz is the spatial characteristic coefficient of the roof area.

[0070] After calculating the stress dynamic risk coefficient, displacement fluctuation coefficient and spatial characteristic coefficient of each grid, establish a data analysis model for the stress dynamic risk coefficient, displacement fluctuation coefficient and spatial characteristic coefficient to generate a roof evaluation coefficient. The roof evaluation coefficients of each grid satisfy:

[0071]

[0072] In the formula, PG represents the roof evaluation coefficient of any grid, α1, α2, α3 are the proportionality coefficients of the stress dynamic risk coefficient, displacement fluctuation coefficient and spatial characteristic coefficient of this grid respectively, α1, α2, α3 are greater than 0, YL dt is the stress dynamic risk coefficient of this grid, WY bd is the displacement fluctuation coefficient of this grid, KJ tz is the spatial characteristic coefficient of the roof area.

[0073] In step S103, since each small area in the roof area of the concealed goaf includes multiple grids, and each grid has a corresponding roof evaluation coefficient. For the convenience of subsequent description, the roof evaluation coefficient of the q-th grid in the u-th small area is represented by PG u, where q = 1, 2, 3, ……, UQ, UQ is the number of grids in the u-th small area, UQ is a positive integer, and q is the grid number in the small area; u = 1, 2, 3, ……, U, U is the number of small areas, U is a positive integer, and u is the small area number. Based on the roof evaluation coefficients of each grid, the stability changes of different grids can be evaluated.

[0074] In step S103, obtain the roof evaluation coefficient threshold, and represent the roof evaluation coefficient threshold by PG yz Compare the roof evaluation coefficients of each grid with the roof evaluation coefficient threshold. If the roof evaluation coefficient is greater than the roof evaluation coefficient threshold, mark the grid as a risk area. If the roof evaluation coefficient is less than the roof evaluation coefficient threshold, no marking will be performed. Thus, the potential risks of different grids on the roof are evaluated.

[0075] Step S104, obtain the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof.

[0076] In step S104, obtaining the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof includes: calculating the roof warning coefficient based on the roof evaluation coefficients of all grids and the roof evaluation coefficient threshold, and determining the stability of the goaf roof based on the roof warning coefficient and the roof warning coefficient threshold. Various thresholds involved in the embodiments of the present invention are usually obtained through methods such as empirical values, numerical simulations, theoretical calculations, and on-site data calibration, or set by staff in the professional field to ensure the accuracy and reliability of the evaluation results

[0077] Specifically, the roof warning coefficient is used to quantify the risk level, and the roof warning coefficients of different small areas satisfy:

[0078]

[0079] In the formula, YJ db, is the roof warning coefficient of the u-th small area, PG u, is the roof evaluation coefficient of the q-th grid in the u-th small area, and UQ is the number of grids in the u-th small area.

[0080] After calculating the roof warning coefficients of each area, obtain the threshold of the roof warning coefficient, and compare each roof warning coefficient with the threshold of the roof warning coefficient. If there is a roof warning coefficient greater than the threshold of the roof warning coefficient, generate a warning signal, indicating that the stability of the roof of the concealed goaf in the open-pit coal mine is poor and risks may occur, and it is necessary to remind the staff to check in time. If all roof warning coefficients are less than the threshold of the roof warning coefficient, no warning signal is generated, indicating that the stability of the roof area is good and there are no risks temporarily.

[0081] To implement the above embodiments, the present invention also proposes a system for evaluating the stability of the roof of a concealed goaf in an open-pit coal mine.

[0082] Figure 2 It is a block diagram of a system for evaluating the stability of the roof of a concealed goaf in an open-pit coal mine provided by an embodiment of the present invention.

[0083] As Figure 2 shown, the system for evaluating the stability of the roof of the concealed goaf in the open-pit coal mine includes a first acquisition module 11, a second acquisition module 12, a calculation module 13, and an evaluation module 14, where:

[0084] The first acquisition module 11 is used to identify the concealed goaf in the open-pit coal mine and obtain the goaf shape information of the concealed goaf.

[0085] The second acquisition module 12 is used to divide the roof area of the concealed goaf by the grid method to obtain the stress information and displacement information of the roof rock stratum of each grid.

[0086] The calculation module 13 is used to obtain the roof evaluation coefficient of each grid based on the goaf shape information, stress information, and displacement information.

[0087] The evaluation module 14 is used to obtain the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof.

[0088] Further, in a possible implementation manner of the embodiment of the present invention, the calculation module 13 is specifically used to: calculate the stress dynamic risk coefficient of each grid based on the stress information; calculate the displacement fluctuation coefficient of each grid based on the displacement information; calculate the spatial characteristic coefficient of the roof area based on the goaf shape information; and obtain the roof evaluation coefficient of each grid based on the stress dynamic risk coefficient, displacement fluctuation coefficient, and spatial characteristic coefficient.

[0089] Further, in a possible implementation manner of the embodiment of the present invention, in the calculation module 13, the roof evaluation coefficient of each grid satisfies:

[0090]

[0091] In the formula, PG represents the roof evaluation coefficient of any grid, and α1, α2, and α3 are the proportionality coefficients of the stress dynamic risk coefficient, displacement fluctuation coefficient, and spatial characteristic coefficient of the grid respectively. YL dt is the stress dynamic risk coefficient of the grid, and WY bd is the displacement fluctuation coefficient of the grid, and KJ tz is the spatial characteristic coefficient of the roof area.

[0092] Furthermore, in a possible implementation manner of the embodiment of the present invention, the stress information includes the horizontal stress and overlying stress per unit time. In the calculation module 13, calculating the stress dynamic risk coefficient of each grid based on the stress information includes: calculating the horizontal stress deviation coefficient and the overlying stress deviation coefficient respectively based on the horizontal stress and overlying stress per unit time of each grid, and then calculating the stress dynamic risk coefficient of the corresponding grid.

[0093] Furthermore, in a possible implementation manner of the embodiment of the present invention, the displacement information includes the initial position of the displacement sensor and the distance from the current position of the displacement sensor to the initial position per unit time. In the calculation module 13, calculating the displacement fluctuation coefficient of each grid based on the displacement information includes: obtaining the average value and standard deviation of the corresponding distance based on the initial position and distance of the displacement sensor of each grid, and then calculating the displacement fluctuation coefficient of the corresponding grid.

[0094] Furthermore, in a possible implementation manner of the embodiment of the present invention, the goaf shape information includes the span, height, and inclination angle of the hidden goaf. In the calculation module 13, calculating the spatial characteristic coefficient of the roof area based on the goaf shape information includes: obtaining the deviations corresponding to the span, height, and inclination angle based on the span, height, and inclination angle, and then calculating the spatial characteristic coefficient of the roof area.

[0095] Furthermore, in a possible implementation manner of the embodiment of the present invention, the evaluation module 14 is specifically configured to: calculate the roof warning coefficient based on the roof evaluation coefficient of all grids and the roof evaluation coefficient threshold, and determine the stability of the goaf roof based on the roof warning coefficient and the roof warning coefficient threshold.

[0096] It should be noted that the foregoing explanation of the embodiment of the method for evaluating the stability of the roof of the hidden goaf in an open-pit coal mine also applies to the system for evaluating the stability of the roof of the hidden goaf in an open-pit coal mine in this embodiment, and will not be elaborated here.

[0097] In an embodiment of the present invention, by identifying the concealed goaf in an open-pit coal mine, the goaf shape information of the concealed goaf is obtained; by using the grid method to divide the roof area of the concealed goaf, the stress information and displacement information of the roof rock strata of each grid are obtained; based on the goaf shape information, stress information and displacement information, the roof evaluation coefficient of each grid is obtained; based on the roof evaluation coefficients of all grids, the roof warning coefficient is obtained to evaluate the stability of the goaf roof. In this case, the roof area is comprehensively divided into grids, and the goaf shape information, stress information and displacement information are comprehensively analyzed to determine the roof evaluation coefficient of the grid, and then the roof warning coefficient is obtained to evaluate the stability of the goaf roof, which is helpful for the stability of the goaf roof, thereby improving the accuracy of the evaluation of the stability of the concealed goaf roof.

[0098] The method and system of the present invention divide the roof area into grids, analyze the information within the grids separately, and comprehensively analyze according to the stress information at the roof rock strata within the grids, the displacement information at the roof rock strata, and the goaf shape information to determine the roof evaluation coefficient of the grid area, and based on the roof evaluation coefficients of all the grids of the roof, measure the degree of danger of the roof stability, which is helpful for analyzing the specific conditions of different areas of the roof in detail and can accurately and real-time evaluate the stability of the roof. It solves the limitation problem of the existing evaluation of the stability of the concealed goaf roof in open-pit coal mines.

[0099] To implement the above embodiment, the present invention also proposes an electronic device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method provided in the foregoing embodiment.

[0100] To implement the above embodiment, the present invention also proposes a computer-readable storage medium, in which computer-executable instructions are stored, and when the computer-executable instructions are executed by a processor, they are used to implement the method provided in the foregoing embodiment.

[0101] To implement the above embodiment, the present invention also proposes a computer program product, including a computer program, and when the computer program is executed by a processor, it implements the method provided in the foregoing embodiment.

[0102] In the descriptions of the foregoing embodiments, the descriptions referring to terms such as "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials, or characteristics described in connection 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 may be combined in any one or more embodiments or examples in a suitable manner. In addition, without contradiction, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0103] In addition, the terms "first" and "second" are used for descriptive purposes only and cannot be construed as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise specifically defined.

[0104] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment, or portion of code including one or more executable instructions for implementing a customized logic function or process, and the scope of the preferred embodiments of the present invention includes additional implementations, where the functions may be executed in a substantially simultaneous manner or in a reverse order according to the functions involved, rather than in the order shown or discussed, which should be understood by those skilled in the art to which the embodiments of the present invention pertain.

[0105] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a definite sequence list of executable instructions for implementing logical functions, and can be specifically implemented in any computer-readable medium for use by an instruction execution system, apparatus, or device (such as a computer-based system, a system including a processor, or other systems that can fetch and execute instructions from the instruction execution system, apparatus, or device), or in conjunction with these instruction execution systems, apparatus, or devices. For the 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 (non-exhaustive list) of computer-readable media include the following: electrical connection parts with one or more wirings (electronic devices), portable computer disk cartridges (magnetic devices), random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber devices, and portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium can even be paper or other suitable media on which the program can be printed, because the program can be obtained electronically, for example, by optically scanning the paper or other media, followed by editing, interpretation, or other suitable processing as necessary, and then stored in a computer memory.

[0106] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any one or a combination of the following techniques well known in the art can be used: discrete logic circuits having logic gate circuits for implementing logical functions on data signals, application-specific integrated circuits having suitable combinational logic gate circuits, programmable gate arrays (PGAs), field programmable gate arrays (FPGAs), etc.

[0107] Those of ordinary skill in the art of this technology can understand that all or part of the steps carried by the method of the above embodiments can be completed by instructing relevant 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 embodiments.

[0108] In addition, each functional unit in various embodiments of the present invention may be integrated into one processing module, may exist physically alone for each unit, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. When the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

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

Claims

1. A method for evaluating the stability of the roof of hidden void areas in open-pit coal mines, characterized in that, Including: Identifying the hidden void area of the open-pit coal mine and obtaining the void shape information of the hidden void area; Dividing the roof area of the hidden void area by the grid method to obtain the stress information and displacement information of the roof rock strata of each grid. Among them, the roof area of the hidden void area is divided into multiple small areas according to a certain scale to form a grid-like division structure; the scale of the division of small areas is determined according to geological conditions, rock strata characteristics and the accuracy to be analyzed. After dividing the roof area of the hidden void area by the grid method, monitoring equipment is reasonably distributed into each grid according to the divided grids. The monitoring equipment includes a stress sensor and a displacement sensor. Among them, the stress information of the roof rock strata includes the horizontal stress and overlying stress per unit time at the roof rock strata; Obtaining the roof evaluation coefficient of each grid based on the void shape information, the stress information and the displacement information; Obtaining a roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the void roof; The obtaining the roof evaluation coefficient of each grid based on the void shape information, the stress information and the displacement information includes: Calculating the stress dynamic risk coefficient of each grid based on the stress information; Calculating the displacement fluctuation coefficient of each grid based on the displacement information; Calculating the spatial characteristic coefficient of the roof area based on the void shape information; Obtaining the roof evaluation coefficient of each grid based on the stress dynamic risk coefficient, the displacement fluctuation coefficient and the spatial characteristic coefficient; The roof evaluation coefficients of each grid satisfy: Wherein, PG represents the roof evaluation coefficient of any grid, and α1, α2, and α3 are the proportionality coefficients of the stress dynamic risk coefficient, displacement fluctuation coefficient, and spatial characteristic coefficient of the grid, respectively, YL dt is the stress dynamic risk coefficient of the grid, WY bd is the displacement fluctuation coefficient of the grid, KJ tz is the spatial characteristic coefficient of the roof area; The obtaining the roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the void roof includes: Calculating to obtain the roof warning coefficient based on the roof evaluation coefficients of all grids and the roof evaluation coefficient threshold; Determining the stability of the void roof based on the roof warning coefficient and the roof warning coefficient threshold.

2. The method for evaluating the roof stability of concealed voids in open-pit coal mines according to claim 1, wherein The stress information includes the horizontal stress and overlying stress per unit time. The calculating the stress dynamic risk coefficient of each grid based on the stress information includes: Based on the horizontal stress and overlying stress per unit time of each grid, calculating the horizontal stress deviation coefficient and the overlying stress deviation coefficient respectively, and then calculating the stress dynamic risk coefficient of the corresponding grid.

3. The method for evaluating the roof stability of concealed voids in open-pit coal mines according to claim 1, wherein The displacement information includes the initial position of the displacement sensor and the distance from the current position of the displacement sensor to the initial position per unit time. The calculating the displacement fluctuation coefficient of each grid based on the displacement information includes: Based on the initial position of the displacement sensor and the distance of each grid, obtaining the average value and standard deviation of the corresponding distance, and then calculating the displacement fluctuation coefficient of the corresponding grid.

4. The method for evaluating the stability of the roof of concealed voids in surface coal mines according to claim 1, wherein The void shape information includes the span, height and inclination angle of the hidden void area. The calculating the spatial characteristic coefficient of the roof area based on the void shape information includes: Obtaining the deviation corresponding to the span, height and inclination angle based on the span, height and inclination angle, and then calculating the spatial characteristic coefficient of the roof area.

5. An evaluation system for the stability of the roof of hidden void areas in open-pit coal mines, characterized in that, Including: A first acquisition module for identifying the hidden void area of the open-pit coal mine and obtaining the void shape information of the hidden void area; A second acquisition module, configured to divide the roof area of the concealed goaf by the grid method to obtain the stress information and displacement information of the roof rock strata of each grid. The roof area of the concealed goaf is divided into multiple small areas according to a certain scale to form a grid-like division structure. The scale of the division of the small areas is determined according to geological conditions, rock stratum characteristics, and the accuracy required for analysis. After dividing the roof area of the concealed goaf by the grid method, monitoring devices are reasonably distributed into each grid according to the divided grids. The monitoring devices include stress sensors and displacement sensors. The stress information of the roof rock strata includes the horizontal stress and overlying stress per unit time at the roof rock strata. A calculation module, configured to obtain the roof evaluation coefficient of each grid based on the goaf shape information, the stress information, and the displacement information. An evaluation module, configured to obtain a roof warning coefficient based on the roof evaluation coefficients of all grids to evaluate the stability of the goaf roof. The calculation module is configured to calculate the stress dynamic risk coefficient of each grid based on the stress information. Calculate the displacement fluctuation coefficient of each grid based on the displacement information. Calculate the spatial characteristic coefficient of the roof area based on the goaf shape information. Obtain the roof evaluation coefficient of each grid based on the stress dynamic risk coefficient, the displacement fluctuation coefficient, and the spatial characteristic coefficient. The roof evaluation coefficients of each grid satisfy: Wherein, PG represents the roof evaluation coefficient of any grid, and α1, α2, and α3 are the proportionality coefficients of the stress dynamic risk coefficient, displacement fluctuation coefficient, and spatial characteristic coefficient of the grid, respectively, YL dt is the stress dynamic risk coefficient of the grid, WY bd is the displacement fluctuation coefficient of the grid, KJ tz is the spatial characteristic coefficient of the roof area; The evaluation module is further configured to calculate and obtain the roof warning coefficient based on the roof evaluation coefficients of all grids and a roof evaluation coefficient threshold. Determine the stability of the goaf roof based on the roof warning coefficient and the roof warning coefficient threshold.

6. An electronic device, characterized in that, It includes: A processor and a memory communicatively connected to the processor; The memory stores computer execution instructions; The processor executes the computer execution instructions stored in the memory to implement the method according to any one of claims 1-4.

7. A computer-readable storage medium, characterized in that, Computer execution instructions are stored in the computer-readable storage medium, and when the computer execution instructions are executed by the processor, they are used to implement the method according to any one of claims 1-4.

Images