Method and device for determining water inflow of mine and water inflow detection system

Abstract

The application provides a method and device for determining the water inflow of a mine and a water inflow detection system. The method comprises: obtaining a plurality of contour maps; obtaining the development heights of all water conducting fracture zones at a plurality of time points, wherein the time points and the development heights correspond to each other; and determining the water inflow of a to-be-mined part of a current working face according to the plurality of contour maps and the plurality of development heights. The scheme solves the problem that the fine prediction of the water inflow of a mine cannot be achieved in the prior art.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of mine water inflow prediction, in particular to a mine water inflow determination method and device, a computer program product and a water inflow detection system. BACKGROUND

[0002] Mine water disaster is one of the five disasters (gas, coal dust, water disaster, fire and roof accident) that endanger the safety production of coal mines. After coal mining, the surrounding rock of the working face is disturbed and damaged, and the underground water flows into the underground through the surrounding rock damage zone, forming the mine water inflow. Accurate prediction of mine water inflow is crucial to mine safety production.

[0003] Current research on mine water inflow prediction methods often focuses on: (1) comprehensively considering the index factors affecting mine water inflow, using existing or mature mathematical theory methods to establish mathematical models for water inflow prediction; (2) generalizing the working face as a well or an overall range area, and using the "big well method" or "water collection corridor method" to predict the water inflow. The current research has achieved rich research results and provides technical support for mine water disaster prevention to some extent. However, the current scheme cannot achieve fine prediction of mine water inflow. SUMMARY

[0004] The main purpose of the present application is to provide a mine water inflow determination method, device, computer program product and water inflow detection system to at least solve the problem that the fine prediction of mine water inflow cannot be achieved in the prior art.

[0005] In order to achieve the above-mentioned purpose, according to one aspect of the present application, a mine water inflow determination method is provided, comprising: obtaining a plurality of contour maps, wherein the contour map is used to represent the relationship between the development height of the water flowing fractured zone of the coal mine working face and the permeability coefficient of the affected aquifer and the thickness of the aquifer, the water flowing fractured zone is the fractured zone that connects the pore water or underground water in the rock layer and makes the water flow through the fractured zone to the mine, the development height of the water flowing fractured zone is the maximum vertical height reached by the fractured zone formed by the disturbance of mining, a plurality of contour maps include one or more of a first contour map, a second contour map and a third contour map, the first contour map is the contour map of the mined part of the current working face, the second contour map is the contour map of the adjacent mined working face of the current working face, and the third contour map is the contour map of the to-be-mined part of the current working face; obtaining the development height of all the water flowing fractured zones at a plurality of time points, wherein the time points and the development heights are one-to-one corresponding; determining the water inflow of the to-be-mined part of the current working face according to a plurality of contour maps and a plurality of development heights.

[0006] Optionally, the obtaining the plurality of contour maps comprises: obtaining pumping test information of the mined part of the current working face, wherein the pumping test information comprises one or more of a layer position of a first aquifer of borehole pumping of the pumping test, a first permeability coefficient of the first aquifer when each borehole pumps, and a first thickness of the first aquifer; obtaining hardness information of overburden of the mined part of the current working face; calculating a first development height of the water flowing fractured zone of the mined part of the current working face according to the hardness information; and generating the first contour map by using ArcGIS software according to the pumping test information and the first development height of the water flowing fractured zone.

[0007] Optionally, the obtaining the plurality of contour maps comprises: obtaining water discharge information of the adjacent mined working face of the current working face, wherein the water discharge information comprises one or more of a water outlet position, a final hole position, a water quantity, a pressure, a water discharge time, and water inrush quantity change information of a water drainage hole; calculating a second permeability coefficient corresponding to the water discharge information according to Darcy's law; determining a second thickness of a second aquifer according to the water discharge information and the second permeability coefficient; obtaining a second development height, wherein the second development height is a development height corresponding to hardness information of overburden of the adjacent mined working face of the current working face; and generating the second contour map by using ArcGIS software according to the second permeability coefficient, the second thickness, and the second development height.

[0008] Optionally, the determining the water inrush quantity of the to-be-mined part of the current working face according to the plurality of contour maps and the plurality of development heights comprises: obtaining a thematic map, wherein the thematic map is determined by screening according to the plurality of contour maps; determining all of the aquifers affected by the to-be-mined part of the current working face according to the plurality of development heights to obtain summary information; and determining the water inrush quantity according to at least the summary information, wherein the permeability coefficients and the thicknesses of different aquifers in the thematic map are different.

[0009] Optionally, the acquiring the thematic map comprises: comparing the first contour map and the third contour map to obtain a first similarity; in a case that the first similarity is greater than or equal to a first similarity threshold, fusing the first contour map and the third contour map to obtain the thematic map, wherein the fusion manner comprises at least one or more of weighted average, Bayesian, Kalman filtering; in a case that the first similarity is less than the first similarity threshold, comparing the second contour map and the third contour map to obtain a second similarity; in a case that the second similarity is greater than or equal to a second similarity threshold, determining the third contour map as the thematic map, wherein the first similarity threshold is greater than the second similarity threshold; in a case that the second similarity is less than the second similarity threshold, fusing the second contour map and the third contour map to obtain the thematic map.

[0010] Optionally, the determining the total water-bearing layer inundated by the to-be-mined part of the current working face according to the plurality of development heights comprises: acquiring a stratigraphic relationship corresponding to the development height of each time point, wherein the stratigraphic relationship is a relationship of the water-bearing layer inundated by the water-fractured zone when the water-fractured zone reaches the development height; determining all the water-bearing layers inundated according to all the corresponding stratigraphic relationships of all the time points to obtain the summary information.

[0011] Optionally, the determining the water inflow amount according to at least the summary information comprises: calculating a plurality of sub-water inflow amounts according to the permeability coefficients of all the water-bearing layers inundated and the thicknesses of all the water-bearing layers corresponding to the summary information by using a large well method, wherein the sub-water inflow amount and the water-bearing layer inundated are one-to-one corresponding; and summing up all the sub-water inflow amounts according to the summary information to obtain the water inflow amount.

[0012] According to another aspect of the present application, a device for determining a water inflow of a mine is provided, comprising: a first obtaining unit configured to obtain a plurality of contour maps, wherein the contour maps are used to represent a relationship between a development height of a water flowing fractured zone of a coal mining face, a permeability coefficient of a water-bearing stratum affected by the coal mining face, and a thickness of the water-bearing stratum, the water flowing fractured zone is a fractured zone through which pore water or underground water in a stratum is connected to enable water to flow to the mine through the fractured zone, the development height of the water flowing fractured zone is a maximum vertical height reached by the fractured zone formed due to mining disturbance, the plurality of contour maps comprises one or more of a first contour map, a second contour map, and a third contour map, the first contour map is the contour map of a mined-out part of a current mining face, the second contour map is the contour map of a neighboring mined-out mining face of the current mining face, and the third contour map is the contour map of a to-be-mined part of the current mining face; a second obtaining unit configured to obtain the development height of the water flowing fractured zone at all time points, wherein the time points and the development heights are in one-to-one correspondence; and a determining unit configured to determine a water inflow of the to-be-mined part of the current mining face according to the plurality of contour maps and the development heights.

[0013] According to still another aspect of the present application, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the steps of any of the methods for determining a water inflow of a mine.

[0014] According to yet another aspect of the present application, a water inflow detection system is provided, comprising one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs comprise instructions for performing any of the methods for determining a water inflow of a mine.

[0015] By applying the technical solution of the present application, the spatial characteristics of a coal mining face can be integrated to determine the changes in the permeability coefficient of a water-bearing stratum and the thickness of the water-bearing stratum in different spaces, and the temporal characteristics of the coal mining face can be integrated to determine the development height of a water flowing fractured zone at different time points, and then the water inflow of a to-be-mined part of a current mining face can be finely predicted according to the temporal characteristics and the spatial characteristics, which is more fine than the way of directly dividing into one whole to predict the water inflow in the prior art, and the spatial and temporal relationship between the development height of the water flowing fractured zone and the water-bearing stratum is used to finely predict the water inflow. BRIEF DESCRIPTION OF DRAWINGS

[0016] The accompanying drawings, which form a part of the present application, are intended to provide further understanding of the present application and are incorporated herein for a purpose of explanations and are not intended as improper limitations to the present application. In the drawings:

[0017] Figure 1 A hardware structure block diagram of a mobile terminal for performing a method for determining a water inflow of a mine according to an embodiment of the present application is shown;

[0018] Figure 2 A flowchart of a method for determining a water inflow of a mine according to an embodiment of the present application is shown;

[0019] Figure 3 A structure block diagram of a device for determining a water inflow of a mine according to an embodiment of the present application is shown.

[0020] Among the above drawings, the following reference signs are included:

[0021] 102, processor; 104, memory; 106, transmission device; 108, input / output device. DETAILED DESCRIPTION

[0022] It should be noted that the embodiments and features of the embodiments in the present application can be combined with each other without conflict. The technical solutions in the embodiments of the present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments. Obviously, the described embodiments are only a part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor should be within the scope of protection of the present application.

[0023] In order to enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor should be within the scope of protection of the present application.

[0024] It should be noted that the terms "first", "second" and the like in the specification and claims of the present application and the above drawings are used to distinguish similar objects, and do not necessarily describe a specific order or sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application described herein can be implemented. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device including a series of steps or units does not necessarily limit to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to these processes, methods, products or devices.

[0025] Existing studies on high-capacity mines present objective challenges such as extremely long and wide working faces, where the heterogeneity of the aquifer in the roof is particularly prominent. Generalizing the working face to a single well or an entire region is inappropriate. Furthermore, the height of the water-conducting fracture zone in the roof exhibits spatiotemporal characteristics during coal seam mining. As the working face advances, the height of the water-conducting fracture zone will extend through different aquifers. Existing studies that only consider a single aquifer are therefore inappropriate.

[0026] As described in the background section, the prior art cannot achieve precise prediction of mine water inflow. To solve the above problems, embodiments of this application provide a method, apparatus, computer program product, and water inflow detection system for determining mine water inflow.

[0027] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.

[0028] The methods and embodiments provided in this application can be executed on a mobile terminal, computer terminal, or similar computing device. Taking running on a mobile terminal as an example, Figure 1 This is a hardware structure block diagram of a mobile terminal for a method of determining the water inflow in a mine, according to an embodiment of the present invention. Figure 1 As shown, a mobile terminal may include one or more ( Figure 1 Only one is shown in the diagram. A processor 102 (which may include, but is not limited to, a microprocessor MCU or a programmable logic device FPGA, etc.) and a memory 104 for storing data are also shown. The mobile terminal may further include a transmission device 106 for communication functions and an input / output device 108. Those skilled in the art will understand that... Figure 1 The structure shown is for illustrative purposes only and does not limit the structure of the mobile terminal described above. For example, the mobile terminal may also include components that are more... Figure 1 The more or fewer components shown, or having the same Figure 1 The different configurations shown.

[0029] The memory 104 can be used to store computer programs, such as software programs of application software and modules, such as a computer program corresponding to the device information display method of the embodiments of the present application. The processor 102 can execute various functional applications and data processing, i.e., implement the above method, by running the computer program stored in the memory 104. The memory 104 can include a high-speed random access memory, and can further include a nonvolatile memory, such as one or more magnetic storage devices, flash memories, or other nonvolatile solid-state memories. In some examples, the memory 104 can further include a memory remotely disposed relative to the processor 102, which can be connected to the mobile terminal through a network. Examples of the network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof. The transmission device 106 is configured to receive or send data through a network. The specific examples of the network can include a wireless network provided by a communication provider of the mobile terminal. In one example, the transmission device 106 includes a network adapter (NIC), which can be connected to other network devices through a base station to communicate with the Internet. In one example, the transmission device 106 can be a radio frequency (RF) module, which is configured to communicate with the Internet in a wireless manner.

[0030] In the embodiments, a method for determining a water inflow of a mine is provided, which is run on a mobile terminal, a computer terminal or a similar computing device. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown herein.

[0031] Figure 2 FIG. 1 is a flowchart of a method for determining a water inflow of a mine according to an embodiment of the present application. As shown in FIG. 1, the method includes the following steps: Figure 2

[0032] ​Step S201, obtaining a plurality of contour maps, wherein the contour maps are used to represent the relationship between the development height of the water flowing fractured zone of the coal mine working face and the permeability coefficient of the affected aquifer and the thickness of the aquifer, the water flowing fractured zone is a fractured zone that connects the pore water or groundwater in the rock formation and allows water to flow into the mine through the fractures, the development height of the water flowing fractured zone is the maximum vertical height reached by the fractured zone formed by the disturbance of mining, the plurality of contour maps includes one or more of a first contour map, a second contour map and a third contour map, the first contour map is the contour map of the mined part of the current working face, the second contour map is the contour map of the adjacent mined working face of the current working face, and the third contour map is the contour map of the to-be-mined part of the current working face.

[0033] Specifically, obtaining a plurality of contour maps means collecting and drawing a series of graphs reflecting the relationship between the development height of the water flowing fractured zone of the mine working face, the permeability coefficient of the aquifer and the thickness of the aquifer. Contour map is a graphical representation method of visualizing spatially continuous parameters (such as permeability coefficient, aquifer thickness) as a series of contour lines, which is used to intuitively show the trend of these parameters in space.

[0034] The water flowing fractured zone refers to the fractured zone that can connect the pore water or groundwater in the rock formation and make water flow into the mine through the fractures during the process of coal mining. The development height of the water flowing fractured zone is a parameter describing the maximum extension range of this fractured zone in the vertical direction, reflecting the connectivity of the fractured zone formed by mining disturbance to the groundwater.

[0035] The first contour map refers to the contour map of the mined part of the current working face. This contour map covers the development height of the water flowing fractured zone, the permeability coefficient and the thickness of the aquifer in the mined area, and is drawn based on the actual mining conditions and the water drainage characteristic index data, which can reflect the hydrogeological characteristics of the mined area.

[0036] The second contour map refers to the contour map of the adjacent mined working face of the current working face. This map is used to analyze the development height of the water flowing fractured zone and the parameters of the aquifer in the adjacent area of the current working face. By comparing the characteristics of the current working face and the adjacent working face, the water hazard risk that may occur in the area to be mined can be inferred, providing a reference for prediction.

[0037] The third contour map refers to the contour map of the to-be-mined part of the current working face. Based on the data analysis and model prediction of the mined area, this map is used to predict the development height of the water flowing fractured zone, the permeability coefficient and the thickness of the aquifer in the area that has not been mined, providing a basis for water inflow prediction and risk assessment of the to-be-mined working face.

[0038] Step S202, obtaining all the development heights of the above-mentioned water flowing fractured zone at multiple time points, wherein the time points and the development heights correspond to each other;

[0039] Specifically, the "multiple time points" refer to multiple key time nodes selected in different stages of coal mining, i.e., from the beginning of mining in the working face to the end of mining. These time points can include the initial mining stage, the initial caving step, the first periodic step, the time point at which the maximum water flowing fractured zone development height is reached, and the duration after mining, etc. Each time point corresponds to a different stage and degree of development of the water flowing fractured zone.

[0040] Specifically, the development height of the water flowing fractured zone refers to the extension range of the fractured zone in the vertical direction during the coal mining process, which is formed due to the collapse and cracking of the roof rock stratum, can connect the pore water or groundwater in the rock stratum, and make these water bodies flow through the fractures into the mine. The size of this height parameter directly affects the size of the mine water inflow and is a key indicator for evaluating the risk of mine water disasters.

[0041] Specifically, the development height of the water flowing fractured zone changes with the advancement of the mining time, and this change process reflects the dynamic development of the fractured zone formed due to mining disturbance, including the formation, expansion, reaching of the maximum height, and possible gradual closure process of the fractured zone. By monitoring and recording the development height of the water flowing fractured zone at each key time point, the relationship between time and the development height of the fractured zone can be established, providing basic data for predicting the mine water inflow. Obtaining the development height data of the water flowing fractured zone at multiple time points is crucial for understanding the spatiotemporal evolution law of the fractured zone, evaluating the change trend of the mine water inflow, and formulating reasonable water disaster prevention measures.

[0042] Step S203, determining the water inflow of the to-be-mined part of the current working face according to the multiple above-mentioned contour maps and the multiple above-mentioned development heights.

[0043] Specifically, the contour maps contain the distribution of key hydrogeological parameters such as the permeability coefficient and the thickness of the aquifer affected by the maximum development height of the water flowing fractured zone. By analyzing these contour maps, the spatial variation law of the aquifer parameters and the correlation between these parameters and the development height of the water flowing fractured zone can be understood. The development height data record the vertical extension range of the water flowing fractured zone at different mining stages, and these data correspond to each time point, representing the dynamic changes of the fractured zone with the mining process. When predicting the water inflow of the to-be-mined part, the development height of the water flowing fractured zone recorded at multiple key time points needs to be combined to understand the extent of the fractured zone reaching the aquifer, and then to evaluate the connectivity of the fractured zone and the aquifer and the water flow path.

[0044] Specifically, the water inflow prediction is not only a static analysis of the aquifer parameters, but also needs to consider the spatial and temporal characteristics of the fracture zone development. This means that the prediction model needs to consider the expansion law of the fracture zone over time and the impact of this expansion on the distribution of aquifer parameters. By combining the contour map with the development height data, it can be more accurately determined which aquifer will be affected by the fracture zone at different mining stages, and how the aquifer parameters affect the water inflow.

[0045] Using the above analysis results, combined with specific hydraulic models (such as the "large well method") and geological conditions, the water inflow of the working face to be mined can be accurately predicted at different time points and different aquifer regions. This prediction result can provide real-time and dynamic guidance for mine water disaster prevention, help mine managers plan water drainage measures in advance, reduce the risk of water disasters during mining, and ensure the safety of mine production.

[0046] Through the embodiment, the spatial characteristics of the coal mine working face can be integrated to determine the changes of the permeability coefficient of the aquifer and the thickness of the aquifer in different spaces, and the temporal characteristics of the coal mine working face can be integrated to determine the development height of the water-conducting fracture zone at different time points. Furthermore, according to the time characteristics and the space characteristics, the water inflow of the to-be-mined part of the current working face can be accurately predicted. Compared with the existing technology of directly dividing into a whole for water inflow prediction, the present embodiment is more accurate, and realizes the accurate prediction of the water inflow based on the spatial and temporal relationship between the development height of the water-conducting fracture zone and the aquifer.

[0047] Specifically, the scheme of the present application effectively realizes the fine prediction of the roof water inflow of the mine from multiple angles such as the scale of the mining area and working face, the temporal and spatial relationship between the overburden damage and the aquifer horizon.

[0048] In the specific implementation process, the plurality of contour maps can be obtained by the following steps: obtaining pumping test information of a mined part of the current working face, wherein the pumping test information includes one or more of the horizon of a first aquifer of a borehole pumping of the pumping test, a first permeability coefficient of the first aquifer during each borehole pumping, and a first thickness of the first aquifer; obtaining hardness information of an overburden of the mined part of the current working face; calculating a first development height of the water-conducting fracture zone of the mined part of the current working face according to the hardness information; and generating the first contour map by using ArcGIS software according to the pumping test information and the first development height of the water-conducting fracture zone.

[0049] In the scheme, the first contour map is generated based on the actual data of the mined working face and GIS technology, which can provide the direct relationship between the development height of the water flowing fractured zone and the hydrogeological parameters of the aquifer, so that the water inflow can be predicted according to the accurate first contour map, which helps to improve the efficiency of mine water disaster prevention.

[0050] Specifically, pumping test data and other data of the mining area range where the working face to be mined is located can be obtained. The pumping test data of the mining area range where the working face to be mined is located specifically refers to the pumping aquifer horizon of the pumping test borehole in the mining area range where the working face to be mined is located, the permeability coefficient of each pumping aquifer of the borehole, the thickness of the aquifer and other data.

[0051] Specifically, the permeability coefficient and thickness contour map of one or more aquifers affected by the maximum development height of the water flowing fractured zone in the mining area range where the working face to be mined is located, i.e., the first contour map, can be drawn.

[0052] Specifically, the steps of drawing the permeability coefficient and thickness contour map of one or more aquifers affected by the maximum development height of the water flowing fractured zone in the mining area range where the working face to be mined is located are as follows:

[0053] (1) According to the "Code for Coal Pillar Setting and Coal Mining Under Buildings, Water Bodies, Railways and Main Roadways" and the like, the maximum development height of the water flowing fractured zone in coal mining is predicted (all taking +, taking the positive value in calculation), and specifically, the maximum development height of the water flowing fractured zone is the maximum distance of the roof rock fissure expanding upward after coal mining, which determines whether the water in the roof aquifer will flow into the mine. When predicting this height, the following formula is usually used for calculation:

[0054] When the overburden strata type of the coal roof is hard:

[0055]

[0056] When the overburden strata type of the coal roof is medium hard:

[0057]

[0058] When the overburden strata type of the coal roof is soft:

[0059]

[0060] In the formula, H Li — height of water flowing fractured zone (m); ∑M — cumulative mining thickness (m).

[0061] (2) According to the predicted maximum development height of the water flowing fractured zone, combined with the stratum structure characteristics and the characteristics of the aquiclude, one or more aquifers within the maximum development height range of the water flowing fractured zone in the mining area are determined; that is, once the maximum development height of the water flowing fractured zone is predicted, combined with the stratum structure characteristics and the characteristics of the aquiclude, it can be determined which aquifers within this height range may be affected. The stratum structure characteristics include the composition, sequence and thickness of the stratum, while the characteristics of the aquiclude mainly focus on the position, thickness and water-resisting performance of the aquiclude. Through this analysis, it can be accurately determined which aquifer water may enter the mine due to mining activities.

[0062] Specifically, the stratum structure characteristics include the depth of the coal seam, the properties of the overburden layer (such as lithology, thickness, sequence, etc.). First, according to the mining depth of the coal seam, the vertical range of the development of the water flowing fractured zone can be determined. Then, the lithology of the overburden layer is analyzed, such as sandstone, limestone, etc., and the fracture development and permeability of these rock layers are key factors affecting the development of the water flowing fractured zone.

[0063] Specifically, the aquiclude refers to the rock layer between aquifers or between an aquifer and a non-aquifer that can prevent or slow down water flow. Identifying the position and thickness of the aquiclude is crucial for determining whether the water flowing fractured zone can penetrate to the next aquifer and predicting the water flow path. The thicker the aquiclude, the better its water-resisting performance, and the less likely it is for the water flowing fractured zone to penetrate it.

[0064] Specifically, the calculated maximum development height of the water flowing fractured zone is compared with the stratum structure diagram to determine all possible aquifers within this height range. If the height of the water flowing fractured zone exceeds the top surface of a certain aquiclude, then this aquifer and the aquifers above it may be affected. For the determined aquifers, further collect hydrogeological parameters related to them, such as permeability coefficient, porosity, aquifer thickness, etc., to establish a model of the aquifer. These parameters can be obtained through historical drilling data, pumping tests, geophysical exploration, etc.

[0065] (3) Combine ArcGIS to draw the permeability coefficient and aquifer thickness contour maps of one or more aquifers affected by the maximum development height of the water flowing fractured zone within the mining area of the working face to be mined. That is, after determining the aquifers that may be affected by the water flowing fractured zone, the next step is to draw the permeability coefficient and aquifer thickness contour maps of these aquifers. The permeability coefficient is a parameter that measures the water permeability of the aquifer, and the aquifer thickness directly affects the water storage capacity of the aquifer. Using ArcGIS software, based on existing pumping test data and drilling water discharge characteristic indicators, etc., the contour maps of these parameters can be drawn through spatial interpolation, etc. The contour map can directly show the spatial distribution of the parameters, which helps to identify the high permeability and thick layer areas of the aquifer, which are potential high-risk areas of water inrush.

[0066] Specifically, hardness refers to the ability of rock to resist external forces such as mining, including compressive strength, shear strength, etc. The hardness information of overburden rock has a direct impact on the formation and expansion of the water-conducting fracture zone. Hard rock may inhibit the formation of cracks, while weak rock is prone to cracking. Based on the hardness information of the overburden rock, the maximum vertical extension range (first development height) of the water-conducting fracture zone in mining can be predicted.

[0067] Specifically, ArcGIS is a powerful geographic information system software that can process and analyze spatial data, generate various types of contour maps, distribution maps, etc. In this process, the pumping test information (including the layer position of the aquifer, the permeability coefficient and the thickness) is combined with the first development height of the water-conducting fracture zone, and the interpolation and spatial analysis functions of ArcGIS are used to generate the first contour map reflecting the relationship between the development height of the water-conducting fracture zone and the hydrogeological parameters of the aquifer. This map can intuitively show the distribution of the permeability coefficient and the thickness of the aquifer within the development height range of the water-conducting fracture zone in the mined-out area, providing key visual information for subsequent water inflow prediction.

[0068] In the specific implementation process, obtaining multiple contour maps can be achieved by the following steps: obtaining the drainage information of the adjacent mined-out area of the current working face, wherein the drainage information includes one or more of the water outlet point position, the final hole position, the water quantity, the pressure, the drainage time and the water inflow change information of the drainage hole; calculating the second permeability coefficient corresponding to the drainage information according to Darcy's law; determining the second thickness of the second aquifer according to the drainage information and the second permeability coefficient; obtaining the second development height, wherein the second development height is the development height corresponding to the hardness information of the overburden rock of the adjacent mined-out area of the current working face; generating the second contour map using ArcGIS software according to the second permeability coefficient, the second thickness and the second development height.

[0069] In this scheme, based on the drainage information of the adjacent mined-out area, the second permeability coefficient and the second thickness of the second aquifer are calculated, and combined with the second development height of the water-conducting fracture zone, the second contour map is generated. Using ArcGIS for data visualization and spatial analysis can make the complex changes of hydrogeological parameters intuitively presented, so that the water inflow can be predicted according to the accurate second contour map, which helps to improve the efficiency of mine water disaster prevention.

[0070] Specifically, the drilling dewatering feature indicators and water inflow of the adjacent mined-out working face of the working face to be mined can be acquired. The acquisition of the drilling dewatering feature indicators and water inflow of the adjacent mined-out working face of the working face to be mined specifically refers to the elevation of the water outflow point of the dewatering hole, the position of the final hole, the real-time monitored water inflow and pressure data, and the corresponding dewatering time, and the water inflow variation law during the advancing of the adjacent mined-out working face of the working face to be mined.

[0071] Specifically, the permeability coefficient and aquifer thickness contour map of one or more aquifers affected by the maximum development height of the water-conducting fracture zone of the adjacent mined-out working face of the working face to be mined, i.e., the second contour map, can be drawn.

[0072] Specifically, the drawing of the permeability coefficient and aquifer thickness contour map of one or more aquifers affected by the maximum development height of the water-conducting fracture zone of the adjacent mined-out working face of the working face to be mined specifically includes the following steps:

[0073] 1. The maximum development height of the water-conducting fracture zone of the mined-out working face is predicted according to the above formula one, formula two and formula three.

[0074] 2. According to the drilling dewatering feature indicators of the adjacent mined-out working face of the working face to be mined, the permeability coefficient of the dewatering hole position is inversely deduced according to Darcy's law. After calculating the permeability coefficients of multiple dewatering holes, the permeability coefficient distribution map of the entire region can be obtained by using the Arcgis interpolation method. That is, according to the dewatering test data of the mined-out working face, the permeability coefficient of the dewatering hole position can be inversely deduced using Darcy's law. Darcy's law describes the linear relationship between the flow velocity of water in porous media and the hydraulic gradient under laminar flow conditions, where the permeability coefficient is a key parameter reflecting the permeability of the aquifer. By analyzing the real-time monitoring data of multiple dewatering holes, the permeability coefficient of each hole can be calculated, and then using the interpolation function of ArcGIS software, these point data can be converted into a continuous permeability coefficient distribution map, which intuitively displays the permeability coefficient distribution of the entire mined-out working face region.

[0075] Specifically, according to Darcy's law, the groundwater flow velocity is proportional to the water head gradient and inversely proportional to the permeability of the underground medium. The mathematical expression of Darcy's law is: Q = K * A * (h1 - h2) / L, where Q represents the amount of groundwater flowing out per unit time, K represents the permeability coefficient, A represents the cross-sectional area, h1 and h2 represent the groundwater level height, and L represents the flow distance.

[0076] Specifically, after the position of the dewatering hole, the real-time monitored water inflow and pressure data, and the corresponding dewatering time are known, the change of the groundwater level height can be calculated from the monitoring data, and thus the groundwater flow velocity can be obtained. Combined with the groundwater flow distance, the permeability coefficient K can be inversely deduced.

[0077] Specifically, by calculating the permeability coefficients of multiple drainage holes, data on permeability coefficients at different locations can be obtained. Using the Arcgis interpolation method, these data can be interpolated to generate a permeability coefficient distribution map for the entire area.

[0078] Specifically, interpolation is a method of estimating unknown point values using known point data. Common interpolation methods include Kriging interpolation and inverse distance weighting interpolation. In Arcgis software, appropriate interpolation methods can be selected, and known permeability coefficient data points can be input into the software for interpolation to generate a permeability coefficient distribution map.

[0079] Specifically, the final permeability coefficient distribution map can help engineers better understand the permeability distribution of underground media and provide important reference for mining work.

[0080] 3. Analyze the water quantity changes during the drainage process by combining real-time monitoring of flow and water pressure data. Based on these changes, further calibrate the thickness of the aquifer. Finally, use the groundwater flow model to fit and analyze the drainage data, combined with the known permeability coefficient, to calculate the thickness information of the aquifer. That is, by analyzing the real-time monitoring of flow and water pressure data, the change rule of water quantity with time during the drainage process can be observed, which helps to calibrate the thickness of the aquifer. The thickness of the aquifer is closely related to the water inflow, and a thicker aquifer can store and release more water. Using the groundwater flow model (such as numerical model or analytical model) to fit and analyze the drainage data, combined with the known permeability coefficient, the thickness information of the aquifer can be calculated. This process relies on the principles of hydraulics and optimization of model parameters, aiming to improve the accuracy of the aquifer model.

[0081] Specifically, through real-time monitoring of flow and water pressure data, the change of water quantity during the drainage process can be obtained. At the beginning of drainage, the flow will gradually increase, and the water level will also rise; when the drainage ends, the flow will gradually decrease, and the water level will gradually decrease. By analyzing the trend of flow and water pressure data, the change of water quantity in the aquifer can be inferred.

[0082] Specifically, calibrating the thickness of the aquifer requires combining real-time monitoring of water level data and the groundwater flow model. According to the change of water level data, it can be inferred whether the thickness of the aquifer is consistent with the expectation. If the actual water level change is inconsistent with the model predicted water level change, the thickness parameter of the aquifer may need to be adjusted to make the model more consistent with the actual situation.

[0083] Specifically, by fitting and analyzing the groundwater flow model to the water release data, and combining this with the known permeability coefficient, the thickness information of the aquifer can be calculated. By comparing the actual observed flow rate, water pressure, and water level data with the model's predicted data, the optimal aquifer thickness parameters can be determined, thus obtaining the aquifer thickness information that best reflects the actual situation.

[0084] 4. Based on the water inflow patterns during the advancement of the mined working face and the contour maps of aquifer permeability and thickness within the mining area, verify the rationality of the distribution of permeability and aquifer thickness within the mined working face. This involves comparing the calculated permeability and aquifer thickness distributions with the actual water inflow patterns monitored during the advancement of the mined working face, and simultaneously comparing them with the contour maps of permeability and aquifer thickness within the mining area. The purpose of this comparative analysis is to verify the rationality of the model's predictions and whether they accurately reflect the hydrogeological characteristics of the aquifer within the mined working face. If the predictions match the actual monitoring data and parameter distribution within the mining area, the model is considered reasonable; otherwise, the model parameters need to be adjusted until the predictions match the actual data.

[0085] Specifically, plot the curve of water inflow at the mined working face as the mining progress progresses, observe the increasing and decreasing trends of water inflow at different mining stages, and whether there are obvious inflection points or peaks. Analyze the contour maps of permeability coefficient and aquifer thickness within the mining area to determine the parameter distribution characteristics of the area where the mined working face is located, such as high-permeability zones and thick aquifer zones.

[0086] Specifically, based on the water inflow data of the mined working face, as well as the permeability coefficient and thickness data of the aquifer, a hydraulic model (such as the "large well method") is used to simulate the changes in water inflow during the mining process. The model should consider the development of water-conducting fracture zones caused by mining. The water inflow trend simulated by the model is compared with the actual monitored water inflow curve to evaluate the degree of agreement between the model prediction and the actual water inflow. If there is a significant difference between the model prediction and the actual data, it may be necessary to adjust the parameters used in the model, such as the permeability coefficient and aquifer thickness, to better reflect the actual geological conditions. This usually involves sensitivity analysis of the model parameters to determine which parameters have the greatest impact on the results. After repeated adjustments and verification, if the model can accurately predict the trend of water inflow changes, then the distribution of permeability coefficient and aquifer thickness within the mining area can be considered reasonable. In addition, key factors affecting water inflow changes can be further analyzed, such as the height of the water-conducting fracture zone and the heterogeneity of the aquifer.

[0087] Specifically, during coal mining, to control the water inflow after working face mining, water drainage measures are usually implemented in the mined working face, including borehole water drainage. The water drainage information includes the water outflow point position, the final hole position, the real-time monitored water inflow, pressure, water drainage time, and the variation law of water inflow. The water outflow point position and the final hole position reflect the arrangement of the water drainage borehole and whether it effectively penetrates the aquifer; the water inflow and pressure data reveal the water storage capacity and water conductivity of the aquifer; the water drainage time and water inflow variation information reflect the dynamic interaction between the water drainage borehole and the aquifer, as well as the response characteristics of the aquifer.

[0088] Specifically, Darcy's law is a fundamental principle describing fluid flow in porous media, which states that the water flow rate is proportional to the water head gradient and the permeability coefficient of the medium under laminar flow conditions. By analyzing the water inflow, pressure, and water drainage time data of the water drainage borehole, the principle of Darcy's law can be applied to inversely calculate the permeability coefficient of the corresponding aquifer, i.e., the second permeability coefficient. This calculation process requires combining the measured water inflow, water pressure, and time data with the geometric parameters of the aquifer (such as porosity, effective thickness, etc.), and using numerical methods or theoretical models to solve.

[0089] Specifically, the thickness of the aquifer is another key parameter affecting water inflow. Based on the analysis of water drainage information, combined with the second permeability coefficient, the thickness of the second aquifer (i.e., the second thickness) can be further determined or calibrated. Specifically, by analyzing the water inflow changes of the water drainage borehole at different time points, the water storage capacity and water flow path of the aquifer can be inferred, and the thickness of the aquifer can be estimated. This process may require the use of groundwater flow models, combined with measured flow and water pressure data, for fitting analysis to obtain the most reasonable thickness estimate.

[0090] Specifically, the second development height refers to the maximum development height of the water-conducting fracture zone in the mined working face adjacent to the current working face. This is usually calculated based on the hardness information of the overlying rock in this area through Formula 1, Formula 2, and Formula 3. The second development height reflects the degree of fracture formation in the rock layer due to mining disturbance, and is an important indicator for evaluating mine water inflow.

[0091] Specifically, the second permeability coefficient, the second thickness, and the second development height and other data are imported into the ArcGIS software, and the spatial analysis and interpolation functions thereof are used to generate a second contour map reflecting the change trend of the hydrogeological parameters of the second aquifer. The specific steps include: data import: importing the second permeability coefficient, the second thickness, and the second development height and other data in the form of spatial distribution into ArcGIS. Spatial interpolation: using the interpolation tools of ArcGIS, such as Kriging interpolation or inverse distance weighted interpolation (IDW), to perform spatial interpolation on the second permeability coefficient and the second thickness data to estimate the parameter distribution of the entire region. Contour map generation: based on the interpolation results, draw the contour map of the second permeability coefficient and the second thickness, which intuitively shows the parameter distribution characteristics of the second aquifer in the adjacent mined-out working face.

[0092] In addition, obtaining a third contour map in the plurality of contour maps can obtain the mining condition-geological condition-hydrogeological condition-drilling water drainage feature index and other data of the working face to be mined. The above-mentioned obtaining the mining condition-geological condition-hydrogeological condition-drilling water drainage feature index and other data of the working face to be mined specifically refers to: the mining condition of the working face to be mined: mining method, coal seam mining thickness, working face strike length, working face inclination length; geological condition: stratum structure characteristics, overburden type; hydrogeological condition: aquiclude characteristics, roof aquifer hydrogeological parameters; drilling water drainage feature index: the position of the water outlet point of the water drainage hole with an upward angle, the position of the final hole, real-time monitoring of water quantity and pressure data, and corresponding water drainage time, etc.

[0093] Specifically, the permeability coefficient and aquifer thickness contour map of one or more aquifers affected by the maximum development height of the water flowing fractured zone of the working face to be mined, i.e., the third contour map, can be drawn.

[0094] Specifically, the above-mentioned drawing the permeability coefficient and aquifer thickness contour map of one or more aquifers affected by the maximum development height of the water flowing fractured zone of the working face to be mined has the same steps as drawing the second contour map, which will not be repeated here.

[0095] Specifically, all relevant data such as mining conditions, geological conditions, hydrogeological conditions, and drilling water drainage characteristic indicators are ensured to be collected. These data include the mining method of the working face, the coal seam mining thickness, the working face size, the stratum structure, the overburden type, the characteristics of the aquifer (permeability coefficient, thickness, etc.), and the measured data of the drilling water drainage (water quantity, water pressure, drainage time, etc.). Based on the real-time monitoring data of the water drainage holes of the mined working face, the permeability coefficient of the water drainage hole position is inversely deduced by applying Darcy's law. This step needs to put the water quantity, water pressure change, and hole depth information of the drilling into the formula of Darcy's law, and the permeability coefficient is obtained by inverse calculation. The data of the permeability coefficient of each water drainage hole position calculated by Darcy's law is imported into GIS software such as ArcGIS. In the software, a suitable interpolation method (such as inverse distance weighted interpolation IDW, Kriging interpolation, etc.) is used to generate a permeability coefficient distribution map of the entire mined working face area. The selection of interpolation parameters depends on the distribution of data points and the requirement for prediction accuracy. The thickness information of the aquifer is inversely deduced by analyzing the water quantity and water pressure change during the drilling water drainage process. These thickness data are also imported into GIS, and an interpolation method is used to draw a distribution map of the aquifer thickness. The maximum development height of the water flowing fractured zone is determined by using the aforementioned prediction formula (based on the overburden type), and then the distribution of the aquifer within this height range is analyzed. This needs to combine the drawn permeability coefficient and aquifer thickness distribution maps with the development height information of the water flowing fractured zone, analyze the change law of the aquifer parameters within this height range, and verify the accuracy of the prediction model. The change law of the water inflow of the mined working face is compared and analyzed with the distribution maps of the permeability coefficient and the aquifer thickness, and whether the model prediction result is consistent with the actual water inflow data is checked. Once the model prediction is consistent with the actual water inflow, these verified parameter distribution maps can be used to predict the water inflow of the to-be-mined working face. In the prediction process, the mining conditions, geological characteristics, and possible development of the water flowing fractured zone of the to-be-mined working face need to be considered, and the drawn parameter distribution maps are combined to make a refined prediction of the water inflow.

[0096] In the implementation process, the water inflow of the to-be-mined part of the current working face is determined according to a plurality of the above contour maps and a plurality of the above development heights, which can be achieved by the following steps: obtaining a thematic map, wherein the thematic map is determined by screening according to a plurality of the above contour maps; determining all the above aquifers affected by the to-be-mined part of the current working face according to a plurality of the above development heights to obtain summary information; and determining the above water inflow according to at least the above summary information, wherein the above permeability coefficient and the above thickness of different aquifers in the above thematic map are different.

[0097] In the scheme, the aquifers possibly affected by the to-be-mined part of the current working face are determined by using a plurality of contour maps and the development height of the water flowing fractured zone, and the water inflow is finally calculated, which can further improve the accuracy of the water inflow prediction.

[0098] Specifically, the working face to be mined can be zoned into one or more thematic maps of aquifer permeability coefficient and aquifer thickness; the development law of the height of the water-conducting fracture zone during the advancement of the working face to be mined can be simulated; and the partitioned and layered fine prediction of the roof water inflow during the advancement of the working face can be made.

[0099] Specifically, the thematic maps are determined by analyzing and screening a plurality of contour maps (such as the first contour map and the second contour map), which reflect the distribution of key hydrogeological parameters such as the permeability coefficient and thickness of the aquifer in different regions and different mining stages. The generation of the thematic maps aims to establish a more comprehensive and detailed aquifer parameter distribution model, which integrates the aquifer information in the mining area, the mined working face, and the working face to be mined, and provides basic data for subsequent water inflow calculation. This process may involve spatial superposition analysis of multiple contour maps to identify patterns and trends of parameter changes.

[0100] Specifically, based on the development height data of multiple key mining stages (such as initial caving, periodic caving, and reaching the maximum development height of the water-conducting fracture zone) of the current working face, the entire aquifer that can be affected by the part to be mined can be determined. This analysis takes into account the dynamic changes of the water-conducting fracture zone with mining time, ensuring that the prediction model can cover all aquifers that may affect water inflow. Since the development height of the water-conducting fracture zone varies with mining depth and mining method, it is necessary to combine the specific working face mining conditions and mining progress to update this information in real time.

[0101] Specifically, the summary information refers to the information set of organizing and summarizing the parameters (permeability coefficient and thickness) of the entire aquifer affected by the part to be mined of the current working face. This information set includes detailed parameters of each aquifer and their specific location and distribution in the thematic map, providing an important basis for the subsequent water inflow prediction. The acquisition of summary information may involve statistical analysis of the parameters of each aquifer, identification of the regularity of parameter changes, and correlation with the development height of the water-conducting fracture zone.

[0102] Specifically, the calculation of water inflow is based on the permeability coefficient and thickness data of different aquifers in the thematic map. Since the parameters of different aquifers may differ greatly, the uniqueness of each aquifer, i.e., the stronger the permeability and the greater the thickness, the greater the contribution to water inflow, needs to be considered when predicting water inflow. This step may use the "big well method", "catchment corridor method" or other hydrogeological models to calculate the water inflow of each aquifer in layers and zones, combined with the mining conditions of the current working face and the development height of the water-conducting fracture zone. Finally, by integrating the water inflow calculation results of all aquifers, the predicted value of the overall water inflow of the working face to be mined can be obtained.

[0103] For example, the water inflow of a certain working face (working face B) in a large coal mine is being predicted. First, the contour maps of the mining area where the working face is located, the adjacent mined working face (working face A), and the mined part of working face B are obtained, which reflect the distribution of the permeability coefficient and thickness of the aquifer in different areas. Then, the development height of the water flowing fractured zone at key stages such as initial caving and first periodic caving of working face B is analyzed, and the aquifers that may be affected at these stages are determined. Next, the data of the parameters of these aquifers, including permeability coefficient and thickness, and their positions in the thematic map, are sorted out to form the summary information.

[0104] Based on the summary information, the "large well method" model is applied to make a refined calculation of the water inflow of each aquifer. In this process, the differences between different aquifer parameters and their interaction with the development height of the water flowing fractured zone are fully considered to ensure the accuracy and reliability of the prediction results. Through model calculation, the predicted values of the water inflow that working face B may encounter at key stages such as initial caving and first periodic caving, as well as the total water inflow prediction during the entire mining period, are obtained.

[0105] In some embodiments, the thematic map can be obtained by comparing the first and third contour maps to obtain a first similarity, and fusing the first and third contour maps to obtain the thematic map when the first similarity is greater than or equal to a first similarity threshold, wherein the fusion method includes at least one or more of weighted average, Bayesian, and Kalman filtering; when the first similarity is less than the first similarity threshold, comparing the second and third contour maps to obtain a second similarity, and determining the third contour map as the thematic map when the second similarity is greater than or equal to a second similarity threshold, wherein the first similarity threshold is greater than the second similarity threshold; and fusing the second and third contour maps to obtain the thematic map when the second similarity is less than the second similarity threshold.

[0106] In this scheme, a more accurate thematic map can be obtained through multi-contour map fusion and screening, which can improve the accuracy of subsequent water inflow prediction.

[0107] Specifically, the division of the one or more aquifer permeability coefficient and aquifer thickness thematic maps of the roof of the working face to be mined is as follows:

[0108] ①, Combined with the aquifer permeability coefficient and thickness thematic map of the mining area, adjacent working face, and the working face to be mined, the variation law of the aquifer permeability coefficient and thickness contour of the working face to be mined and the aquifer permeability coefficient and thickness contour of the area where the working face to be mined is located within the mining area range is compared and analyzed.

[0109] ②, If the aquifer permeability coefficient and thickness contour of the working face to be mined is basically consistent with the aquifer permeability coefficient and thickness contour of the area where the working face to be mined is located within the mining area range, the average value of the two is taken as the aquifer permeability coefficient and thickness contour thematic map of the roof of the working face to be mined.

[0110] ③, If the aquifer permeability coefficient and thickness contour of the working face to be mined is quite different or inconsistent with the aquifer permeability coefficient and thickness contour of the area where the working face to be mined is located within the mining area range; further compare the variation law of the aquifer permeability coefficient and thickness contour of the working face to be mined and the aquifer permeability coefficient and thickness contour of the adjacent mined working face, and analyze the extension trend of the aquifer permeability coefficient and thickness contour of the adjacent mined working face to the working face to be mined, which is consistent with the trend of the aquifer permeability coefficient and thickness contour of the working face to be mined, then the aquifer permeability coefficient and thickness contour thematic map of the roof of the working face to be mined is selected and used.

[0111] ④, Combined with the natural fracture level zoning in ArcGIS spatial analysis, the aquifer permeability coefficient and thickness contour thematic map of the roof of the working face to be mined is divided into four level areas: I, II, III and IV. Among them, I level represents large aquifer permeability coefficient / thickness.

[0112] Specifically, the development law of the simulated water flowing fractured zone height during the advancing process of the working face to be mined is as follows:

[0113] I. According to the stratum structure characteristics of the working face to be mined, combined with a three-dimensional numerical simulation software, a numerical model of the working face scale is established;

[0114] II. Combined with the actual mining situation of the working face to be mined, the development height of the water flowing fractured zone at key nodes such as the initial caving step distance, the first cycle step distance, and the initial advancing distance to reach the maximum water flowing fractured zone height is obtained by simulating the step-by-step excavation of the coal seam;

[0115] III. Summarize the evolution law of the water flowing fractured zone development height during the mining process of the working face to be mined.

[0116] Specifically, the first similarity is quantifying the similarity degree of the two in the distribution of aquifer parameters by comparing the first and third contour maps. This can be achieved through various statistical or image processing methods, such as calculating the correlation coefficient, structural similarity (SSIM) index between the two maps, or using a machine learning model for feature matching.

[0117] Specifically, if the first similarity reaches or exceeds a preset first similarity threshold, it means that the hydrogeological parameters of the mined-out face can represent or predict the conditions of the to-be-mined face to some extent. At this time, the weighted average method can be used to fuse the parameter data of the first and third contour maps to generate a more accurate thematic map.

[0118] The Bayesian method is a probability-based fusion technique that can dynamically update the prediction of aquifer parameter distribution based on prior knowledge and observed data. In water inflow prediction, the Bayesian formula can be used to combine the data of the first and third contour maps to evaluate the most likely distribution of aquifer parameters.

[0119] Kalman filtering is an effective state estimation method for processing noisy measurement data. In this method, the first contour map can be regarded as the "observation value" and the third contour map as the "prediction value". Through the Kalman filtering algorithm, the information of the two is fused to obtain a more accurate estimation of the aquifer parameter distribution, thereby generating a thematic map.

[0120] Specifically, if the first similarity is lower than the set threshold, it indicates that the hydrogeological conditions of the mined-out face may differ significantly from the to-be-mined face, and data fusion cannot be directly performed. At this time, the similarity between the second contour map (aquifer parameter distribution map of the adjacent mined-out face) and the third contour map (second similarity) will be compared.

[0121] The second similarity threshold is lower than the first similarity threshold, meaning that even if the similarity between the second and third contour maps is low, it may still contain useful information for generating a thematic map. This is because the mining conditions and geological structure of adjacent faces may be closer, so their aquifer parameter distribution may be more valuable for predicting the water inflow of the to-be-mined face.

[0122] Specifically, if the second similarity reaches or exceeds the second similarity threshold, the third contour map can be directly used as a thematic map without additional fusion steps. This indicates that the aquifer parameter distribution of the to-be-mined face is consistent with the adjacent mined-out face, which can be directly used as the basis for water inflow prediction.

[0123] Specifically, if the second similarity is still lower than the second similarity threshold, it indicates that the hydrogeological conditions of the to-be-mined face differ significantly from the mined-out face, and at this time, the second and third contour maps need to be fused to generate a thematic map that integrates the information of both, for further water inflow prediction.

[0124] For example, in practical applications, assume that the predicted is the water inflow of working face C, and working faces A and B are its mined and adjacent mined working faces. First, the contour maps of aquifer parameter distribution of working faces A, B, and C are obtained respectively. It is found that the first contour map of working face A has a large difference in permeability coefficient and thickness distribution from the third contour map of working face C, and the first similarity is low, failing to reach the first similarity threshold. However, the second contour map of working face B shows a high similarity to the third contour map of working face C in some parameters, exceeding the lower second similarity threshold. Therefore, it is decided to use the weighted average method to fuse the second contour map of working face B and the third contour map of working face C to generate the final thematic map for water inflow prediction.

[0125] In order to perform weighted averaging, different weights are assigned to the parameters of working face B and working face C to reflect their relative importance for water inflow prediction. For example, if working faces B and C are very similar in geological structure and mining conditions, the parameters of working face B can be given a higher weight. The fused thematic map integrates the aquifer parameter information of the two working faces, providing a more comprehensive and accurate basis for water inflow prediction of working face C. Through this similarity analysis and data fusion strategy, mine managers can better understand the hydrogeological differences between different working faces and adjust water drainage measures accordingly to ensure safe production and effective management of water resources in the mine.

[0126] In practical applications, the first similarity threshold can be 90%, and the second similarity threshold can be 85%. Of course, the first similarity threshold can be 70%, and the second similarity threshold can be 60%. The specific values are not limited, and those skilled in the art can make adaptive extensions according to actual conditions.

[0127] In some embodiments, according to a plurality of development heights, all aquifers of a to-be-mined part of the current working face are determined to obtain summary information, which can be achieved by the following steps: obtaining a stratigraphic relationship corresponding to the development height at each time point, wherein the stratigraphic relationship is a relationship of the aquifer reached by the water-conducting fractured zone when the water-conducting fractured zone reaches the development height; and determining all aquifers reached according to all corresponding stratigraphic relationships of all time points to obtain the summary information.

[0128] In this scheme, by comprehensively analyzing the laws of aquifer parameters and time changes and considering the influence of the development of the water-conducting fractured zone on different aquifers during mining, the stratigraphic relationship corresponding to the development height at each time point can be obtained, and all aquifers reached can be determined according to these information to generate summary information, which can determine all aquifers reached during mining, and then subsequent fine prediction of water inflow can be performed.

[0129] Specifically, during the coal mining process, as the working face advances, the overlying rock strata of the coal seam will be subjected to different mining disturbances, leading to changes in the development height of the water-conducting fracture zone over time. These changes will directly affect whether the water-conducting fracture zone will spread and to which aquifers. Therefore, for each key time point during the mining process, we need to obtain the development height of the water-conducting fracture zone and determine the aquifers reached by the fracture zone at this height. The stratigraphic relationship refers to the spatial relationship between the development height of the water-conducting fracture zone and the mine roof aquifers at a specific time point, including whether the fracture zone penetrates the aquifers, reaches which aquifers, and the relative position with the aquifers.

[0130] Specifically, after analyzing the stratigraphic relationships at multiple key time points during the mining process, we need to organize and summarize this information to determine all the aquifers reached by the water-conducting fracture zone during the entire mining period. This includes the aquifers that may be encountered during the initial mining, initial caving, periodic caving, and post-mining stages. Determine the reached aquifers: By comparing the development heights at different time points with the stratigraphic positions of the aquifers, identify the specific aquifers reached by the water-conducting fracture zone at each mining stage, as well as the parameter changes of these aquifers, such as permeability coefficient and thickness.

[0131] Specifically, summarizing information refers to organizing and synthesizing a collection of hydrogeological parameters (such as permeability coefficient, thickness, etc.) of all the reached aquifers and their spatiotemporal change information during the mining process. This information collection is crucial for predicting water inflow, as it directly reflects the potential water inflow risks after the working face is mined. Generate summary information: Convert the stratigraphic relationships at each time point into data of specific aquifer parameters, and generate a summary information list containing all the reached aquifer parameters and their change trends based on the development law of the water-conducting fracture zone.

[0132] In some embodiments, the water inflow is determined based on the above-mentioned summary information, which can be achieved by the following steps: according to the permeability coefficients of all the aquifers reached by the summary information and the thicknesses of all the aquifers, a plurality of sub-water inflows are calculated using the large well method, wherein the sub-water inflows correspond one-to-one to the reached aquifers; according to the summary information, the sum of all the sub-water inflows is calculated to obtain the water inflow.

[0133] In this scheme, the sum of all the sub-water inflows can be calculated, which not only allows the prediction of sub-water inflows at different mining stages of the working face, but also allows the prediction of the overall water inflow. Since the spatiotemporal characteristics during the mining process are considered, the accuracy of the prediction results is further improved.

[0134] Specifically, the above-mentioned roof water inrush amount zoning and layering fine prediction during working face advancement has the following specific steps:

[0135] (i) Analyze the relationship between the height of the water-conducting fracture zone during the advancement of the working face to be mined and the layer position of the roof aquifer, focusing on analyzing the water-conducting fracture zone height corresponding to the initial caving step distance, the first periodic step distance, the maximum height of the water-conducting fracture zone reached initially, and the maximum height of the water-conducting fracture zone sustained (4 periods), and the relationship between the roof aquifer layer position and one or more aquifers within the water-conducting fracture zone height range during the 4 periods;

[0136] (ii) Using the "large well method", calculate the roof water inrush amount for each of the 4 periods, different roof aquifer layer positions, and different roof aquifer regions during the advancement of the working face, based on the zoning of the roof aquifer into four regions;

[0137] (iii) Statistically analyze the evolution law of the roof water inrush amount during the advancement of the working face, and after integrating the calculation results of the water inrush amount for each period, the total water inrush amount for the entire working face can be obtained.

[0138] Specifically, the large well method is a commonly used hydrogeological model for estimating the inflow of groundwater in a certain area. In the prediction of roof water inrush in a mine, the water-conducting fracture zone formed after the working face is mined can be considered as a large "well". By calculating the sub-inrush amount of each aquifer affected by the fracture zone, the water inrush risk of the entire working face can be evaluated. Sub-inrush amount: refers to the water inrush amount corresponding to each affected aquifer, which is directly related to the permeability coefficient K and thickness T of the aquifer.

[0139] Specifically, for each aquifer, based on its permeability coefficient K, thickness T, and the spatiotemporal characteristics of the water-conducting fracture zone affecting the layer, the large well method model is applied to calculate the sub-inrush amount. This usually involves using the groundwater flow equation in combination with the specific parameters of the aquifer to estimate the amount of water transported from the aquifer to the working face at different mining stages.

[0140] Specifically, after calculating the sub-inrush amount of all affected aquifers, the total inrush amount of the entire working face at a specific mining stage can be obtained by summing these sub-inrush amounts. This process ensures that the prediction result takes into account all factors that may affect the inrush amount, improving the comprehensiveness and accuracy of the prediction. Water inrush amount: is the total amount of water transported from the roof aquifer to the underground after the working face is mined, which is the sum of the sub-inrush amounts of all affected aquifers.

[0141] The following is a specific implementation case to better understand this process:

[0142] Case: Roof water inrush prediction for working face E in mine E.

[0143] Step A: Determine key time points and height of water-conducting fracture zone development;

[0144] Initial collapse stage: In the early stage of mining, as the coal seam is mined, the water-conducting fracture zone begins to develop, with a development height recorded as H1. During this stage, the water-conducting fracture zone affects the uppermost aquifer L1.

[0145] First cycle collapse stage: As the working face advances further, the development height of the water-conducting fracture zone increases to H2, affecting aquifers L1 and L2.

[0146] Stable mining stage: At a certain stable stage during mining, the development height of the water-conducting fracture zone stabilizes at H3, affecting aquifers L2, L3, and L4.

[0147] Final collapse stage: In the last stage before the completion of mining, the development height of the water-conducting fracture zone reaches a maximum of H4, affecting aquifers L4 and L5.

[0148] Step B: Analyze the stratigraphic relationship at each time point;

[0149] At the initial collapse stage, the development height of the water-conducting fracture zone H1 affects aquifer L1, and the parameters of L1 (such as permeability coefficient K1 and thickness T1) are recorded.

[0150] At the first cycle collapse stage, the development height of the water-conducting fracture zone H2 affects aquifers L1 and L2, and the parameters of L1 are updated (such as permeability coefficient K1 and thickness T1 may change due to the expansion of the fracture zone), and the parameters of L2 (such as permeability coefficient K2 and thickness T2) are recorded.

[0151] At the stable mining stage, the development height of the water-conducting fracture zone H3 affects aquifers L2, L3, and L4, and the parameters of L2 and L3 (such as permeability coefficient K2, thickness T2, permeability coefficient K3, and thickness T3) are recorded, as well as the parameters of L4 (such as permeability coefficient K4 and thickness T4).

[0152] At the final collapse stage, the development height of the water-conducting fracture zone H4 affects aquifers L4 and L5, and the parameters of L4 are updated, as well as the parameters of L5 (such as permeability coefficient K5 and thickness T5).

[0153] Step C: Integrate the stratigraphic relationship at all time points to form a summary information;

[0154] Integrate the stratigraphic relationship information of the initial collapse, first cycle collapse, stable mining stage, and final collapse stage to form a comprehensive summary information list. This list not only includes all the aquifers (L1, L2, L3, L4, L5) affected by the water-conducting fracture zone during mining, but also details the parameter changes of each aquifer at different mining stages, such as permeability coefficient and thickness.

[0155] The summary information list can be as follows:

[0156] Aquifer L1: Impacted in the initial caving stage, with permeability coefficient K1, thickness T1.

[0157] Aquifer L2: Impacted in the initial caving and first periodic caving stages, with permeability coefficient K2 (which can vary at different stages), thickness T2.

[0158] Aquifer L3: Impacted in the stable mining stage, with permeability coefficient K3, thickness T3.

[0159] Aquifer L4: Impacted in the first periodic caving, stable mining stage, and final caving stage, with permeability coefficient K4 (which can vary at different stages), thickness T4.

[0160] Aquifer L5: Impacted in the final caving stage, with permeability coefficient K5, thickness T5.

[0161] Step D: Water inflow prediction based on the summary information;

[0162] Using the above summary information, combined with the "large well method" or other hydrogeological models, for each impacted aquifer (L1, L2, L3, L4, L5), at different mining stages (initial caving, first periodic caving, stable mining stage, final caving stage), the water inflow is calculated respectively. Finally, by integrating the water inflow prediction results of all aquifers at each mining stage, the roof water inflow prediction of the entire working face E can be obtained, including the water inflow change trend at the early, middle, and late stages of mining.

[0163] Suppose the roof water inflow of working face E is predicted. In the initial caving stage, the water-conducting fracture zone impacts aquifer L1, with a permeability coefficient K1 = 1 x 10^-5 m / s and a thickness T1 = 20 m; in the first periodic caving stage, it impacts aquifers L1 and L2, with L2 having a permeability coefficient K2 = 2 x 10^-5 m / s and a thickness T2 = 25 m; in the stable mining stage, it impacts aquifers L2, L3, and L4, with L3 and L4 having a permeability coefficient K3 = 1.5 x 10^-5 m / s, a thickness T3 = 30 m, and a permeability coefficient K4 = 1.8 x 10^-5 m / s, a thickness T4 = 35 m, respectively; in the final caving stage, it impacts aquifers L4 and L5, with L5 having a permeability coefficient K5 = 2.2 x 10^-5 m / s and a thickness T5 = 40 m.

[0164] According to the above parameters, the sub-inrush water quantity of the initial caving stage L1, the sub-inrush water quantity of the first periodic caving stage L1 and L2, the sub-inrush water quantity of the stable mining stage L2, L3 and L4, and the sub-inrush water quantity of the final caving stage L4 and L5 can be calculated respectively. Assuming that in the initial caving stage, the sub-inrush water quantity of L1 is calculated as Q1=100m 3 / day; in the first periodic caving stage, the sub-inrush water quantities of L1 and L2 are calculated as Q1'=120m 3 / day, Q2=150m 3 / day respectively; in the stable mining stage, the sub-inrush water quantities of L2, L3 and L4 are calculated as Q2'=180m 3 / day, Q3=200m 3 / day, Q4=220m 3 / day respectively; in the final caving stage, the sub-inrush water quantities of L4 and L5 are calculated as Q4'=240m 3 / day, Q5=260m 3 / day respectively.

[0165] Finally, all the sub-inrush water quantities are summed up to obtain the roof inrush water quantity of the entire working face E. Taking the initial caving stage as an example, since only L1 is affected at this time, the total inrush water quantity is Q1=100m 3 / day; in the first periodic caving stage, the total inrush water quantity is Q1'+Q2=120+150=270m 3 / day; in the stable mining stage, the total inrush water quantity is Q2'+Q3+Q4=180+200+220=600m 3 / day; in the final caving stage, the total inrush water quantity is Q4'+Q5=240+260=500m 3 / day.

[0166] In this case, by carefully analyzing the development height of the water-conducting fracture zone and the wave-through relationship with the aquifer, and integrating the aquifer parameter information at the corresponding time points, a detailed, phased and layered parameter basis can be provided for the roof inrush water quantity prediction of the working face E of the mine. This prediction method not only considers the dynamic changes of the water-conducting fracture zone during mining, but also fully evaluates the contribution of different aquifers to the inrush water quantity, thereby improving the accuracy of the prediction and providing a scientific basis for mine water disaster prevention and inrush water quantity management.

[0167] In summary, the scheme is aimed at large-capacity mines, and the objective situation of super-long and super-wide working face, and the heterogeneity of the roof aquifer is particularly prominent. The application analyzes the zoning of the hydrogeological parameters such as the permeability coefficient and the thickness of one or more aquifers in the roof from the perspective of the working face angle of the working face to be mined; the height of the roof water-conducting fracture zone has the characteristics of space and time during the coal mining process, and the water-conducting fracture law of four time periods, including the initial caving step distance, the water-conducting fracture zone height corresponding to the first cycle step distance, the maximum height of the water-conducting fracture zone reached initially, and the maximum height of the water-conducting fracture zone, is comprehensively considered, and the space-time relationship between the water-conducting fracture law and the roof aquifer is analyzed, so that the fine prediction of the roof water inflow of the working face to be mined in the time dimension, the sectional dimension and the layered dimension is realized.

[0168] The embodiment of the present application also provides a device for determining the water inflow of a mine. It should be noted that the device for determining the water inflow of a mine in the embodiment of the present application can be used to execute the method for determining the water inflow of a mine provided by the embodiment of the present application. The device is used to realize the above-mentioned embodiment and preferred embodiment, and will not be described here. As used below, the term "module" can be a combination of software and / or hardware that realizes a predetermined function. Although the device described in the following embodiment is preferably realized in software, the realization of hardware or a combination of software and hardware is also possible and is conceived.

[0169] The device for determining the water inflow of a mine provided by the embodiment of the present application is introduced below.

[0170] Figure 3 is a structural block diagram of a device for determining the water inflow of a mine according to the embodiment of the present application. As shown in Figure 3 , the device comprises:

[0171] A first acquisition unit 10 is configured to acquire a plurality of contour maps, wherein the contour maps are used to represent the relationship between the development height of a water-conducting fracture zone of a coal mine working face and the permeability coefficient of the affected aquifer and the thickness of the aquifer, the water-conducting fracture zone is a fracture zone that connects the pore water or underground water in the rock stratum to make the water flow through the fracture zone to the mine, the development height of the water-conducting fracture zone is the maximum vertical height reached by the fracture zone formed due to the disturbance of mining, the plurality of contour maps comprises one or more of a first contour map, a second contour map and a third contour map, the first contour map is the contour map of the mined part of the current working face, the second contour map is the contour map of the adjacent mined working face of the current working face, and the third contour map is the contour map of the to-be-mined part of the current working face;

[0172] The second acquisition unit 20 is configured to acquire all of the development heights of the water flowing fractured zones at a plurality of time points, wherein the time points and the development heights are in one-to-one correspondence.

[0173] The determination unit 30 is configured to determine the water inflow of the to-be-mined part of the current working face according to the plurality of contour maps and the plurality of development heights.

[0174] In the implementation process, the first acquisition unit includes a first acquisition module, a second acquisition module, a first calculation module and a first generation module. The first acquisition module is configured to acquire pumping test information of the mined part of the current working face, wherein the pumping test information includes one or more of the following: the layer position of the first aquifer of borehole pumping in the pumping test, the first permeability coefficient of the first aquifer during each borehole pumping, and the first thickness of the first aquifer. The second acquisition module is configured to acquire hardness information of the overburden strata of the mined part of the current working face. The first calculation module is configured to calculate the first development height of the water flowing fractured zone of the mined part of the current working face according to the hardness information. The first generation module is configured to generate the first contour map by using ArcGIS software according to the pumping test information and the first development height of the water flowing fractured zone.

[0175] In this scheme, the first contour map is generated based on the actual data of the mined working face and GIS technology, which can provide the direct relationship between the development height of the water flowing fractured zone and the hydrogeological parameters of the aquifer, so that the water inflow can be predicted according to the accurate first contour map, which helps to improve the efficiency of mine water disaster prevention.

[0176] In the implementation process, the first acquisition unit includes a third acquisition module, a second calculation module, a first determination module, a fourth acquisition module and a second generation module. The third acquisition module is configured to acquire drainage information of the adjacent mined working face of the current working face, wherein the drainage information includes one or more of the following: the water outlet position of the water drainage hole, the final hole position, the water quantity, the pressure, the drainage time and the water inflow change information. The second calculation module is configured to calculate a second permeability coefficient corresponding to the drainage information according to Darcy's law. The first determination module is configured to determine a second thickness of the second aquifer according to the drainage information and the second permeability coefficient. The fourth acquisition module is configured to acquire a second development height, wherein the second development height is the development height corresponding to the hardness information of the overburden strata of the adjacent mined working face of the current working face. The second generation module is configured to generate the second contour map by using ArcGIS software according to the second permeability coefficient, the second thickness and the second development height.

[0177] In the scheme, the second permeability coefficient and the second thickness of the second aquifer are calculated based on the water discharge information of the adjacent mined working face, and the second development height of the water flowing fractured zone is combined to generate a second contour map. Data visualization and spatial analysis are performed using ArcGIS, which can intuitively present the changes of complex hydrogeological parameters, so that the water inflow can be predicted according to the accurate second contour map, and the efficiency of mine water disaster prevention and control can be improved.

[0178] In the implementation process, the determining unit includes a fifth acquisition module, a second determining module and a third determining module. The fifth acquisition module is configured to acquire a thematic map, wherein the thematic map is determined by screening a plurality of contour maps. The second determining module is configured to determine all aquifers affected by the to-be-mined part of the current working face according to a plurality of development heights, and obtain summary information. The third determining module is configured to determine the water inflow according to the summary information, wherein the permeability coefficients and the thicknesses of different aquifers in the thematic map are different.

[0179] In the scheme, the aquifers possibly affected by the to-be-mined part of the current working face are determined by using a plurality of contour maps and the development height of the water flowing fractured zone, and the water inflow is finally calculated, so that the accuracy of water inflow prediction can be further improved.

[0180] In some embodiments, the fifth acquisition module includes a first comparison submodule, a first fusion submodule, a second comparison submodule, a first determination submodule and a second fusion submodule. The first comparison submodule compares the first contour map and the third contour map to obtain a first similarity. The first fusion submodule is configured to fuse the first contour map and the third contour map to obtain the thematic map when the first similarity is greater than or equal to a first similarity threshold. The fusion method includes at least one or more of weighted average, Bayesian and Kalman filtering. The second comparison submodule is configured to compare the second contour map and the third contour map to obtain a second similarity when the first similarity is less than the first similarity threshold. The first determination submodule is configured to determine the third contour map as the thematic map when the second similarity is greater than or equal to a second similarity threshold. The first similarity threshold is greater than the second similarity threshold. The second fusion submodule is configured to fuse the second contour map and the third contour map to obtain the thematic map when the second similarity is less than the second similarity threshold.

[0181] In the scheme, the thematic map can be obtained by fusing and screening a plurality of contour maps, and the accuracy of subsequent water inflow prediction can be improved.

[0182] In some embodiments, the second determining module comprises an acquisition submodule and a second determining submodule. The acquisition submodule is configured to acquire a layer relationship corresponding to the development height at each of the time points, wherein the layer relationship is a relationship of the aquifer affected by the water-conducting fractured zone when the water-conducting fractured zone reaches the development height; and the second determining submodule is configured to determine all the aquifers affected according to all the corresponding layer relationships at all the time points to obtain the summary information.

[0183] In this scheme, by comprehensively analyzing the parameters of the aquifer and the time-varying law and considering the influence of the development of the water-conducting fractured zone on different aquifers during the mining process, the layer relationship corresponding to the development height at each time point can be acquired, and all the aquifers affected can be determined according to these information to generate the summary information. In this way, all the aquifers affected during the mining process can be determined, and the subsequent fine prediction of the water inflow can be performed.

[0184] In some embodiments, the third determining module comprises a calculation submodule and a statistical submodule. The calculation submodule is configured to calculate a plurality of sub-water inflows by using the large well method according to the permeability coefficients of all the aquifers affected corresponding to the summary information and the thicknesses of all the aquifers; and the statistical submodule is configured to obtain the water inflow by summing up all the sub-water inflows according to the summary information.

[0185] In this scheme, the sum of all the sub-water inflows can be calculated, so that the sub-water inflow prediction of the working face at different mining stages can be obtained, and the overall water inflow can also be obtained. Since the time and space characteristics during the mining process are considered, the accuracy of the prediction result is further improved.

[0186] The device for determining the water inflow of the mine comprises a processor and a memory. The first acquisition unit, the second acquisition unit, and the determining unit are all stored in the memory as program units, and the corresponding functions are realized by the processor executing the program units stored in the memory. The modules are located in the same processor, or the modules are located in different processors in any combination.

[0187] The processor comprises a core, and the core retrieves the corresponding program units from the memory. The core can be set to one or more, and the problem that the fine prediction of the water inflow of the mine cannot be realized in the prior art can be solved by adjusting the core parameters.

[0188] The memory can include a non-persistent memory in a computer readable medium, a random access memory (RAM), and / or a non-volatile memory such as a read-only memory (ROM) or a flash memory (flash RAM), and the memory comprises at least one memory chip.

[0189] The embodiment of the present application provides a computer readable storage medium, the computer readable storage medium comprises a stored program, wherein the computer readable storage medium controls a device where the computer readable storage medium is located to execute the mine water inflow determination method when the program runs.

[0190] The embodiment of the present application provides a processor, the processor is used for running a program, wherein the processor executes the mine water inflow determination method when the program runs.

[0191] The embodiment of the present application provides a device, the device comprises a processor, a memory and a program stored in the memory and capable of running on the processor, and the processor executes the program to realize at least the mine water inflow determination method steps. The device herein can be a server, a PC, a PAD, a mobile phone and the like.

[0192] A computer program product comprises a non-volatile computer readable storage medium, the non-volatile computer readable storage medium stores a computer program, and the computer program is executed by a processor to realize the steps of the mine water inflow determination method in various embodiments of the present application.

[0193] The present application also provides a water inflow detection system, comprising one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs comprise a program for executing any of the mine water inflow determination methods.

[0194] Obviously, those skilled in the art should understand that the modules or steps of the present application described above can be realized by general computing devices, which can be concentrated on a single computing device or distributed on a network composed of multiple computing devices, and they can be realized by program codes executable by the computing devices, so that they can be stored in storage devices and executed by the computing devices, and in some cases, the steps shown or described can be executed in different order, or they can be manufactured into individual integrated circuit modules, or multiple modules or steps can be manufactured into a single integrated circuit module. Therefore, the present application is not limited to any specific combination of hardware and software.

[0195] Those skilled in the art will appreciate that embodiments of the application can be readily used as software, hardware, or a combination of software and hardware. In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0196] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.

[0197] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.

[0198] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.

[0199] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0200] The memory can include non-persistent memory, random access memory (RAM), and / or non-volatile memory, etc. in the form of a computer-readable medium, such as read only memory (ROM) or flash memory. The memory is an example of computer-readable media.

[0201] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible to a computing device. According to the definition herein, computer-readable media does not include transitory media such as modulated data signals and carriers.

[0202] It should also be noted that the terms "comprising", "including" or any other variant thereof are intended to cover non-exclusive inclusion, so that processes, methods, articles or devices including a series of elements not only include those elements, but also include other elements not explicitly listed or inherent to such processes, methods, articles or devices. Without more limitations, the element defined by the statement "comprising a" does not exclude the presence of additional identical elements in the process, method, article or device including the element.

[0203] From the above description, it can be seen that the above-mentioned embodiments of the present application achieve the following technical effects:

[0204] 1) The mine water inflow determination method of the present application can integrate the spatial characteristics of the coal mine working face to determine the permeability coefficient of the aquifer in different spaces and the thickness variation of the aquifer, and can also integrate the time characteristics of the coal mine working face to determine the development height of the water flowing fractured zone at different time points, and then according to the time characteristics and the spatial characteristics, the water inflow of the to-be-mined part of the current working face is finely predicted. Compared with the way of directly dividing into a whole to predict the water inflow in the prior art, the present application is more fine, and realizes the fine prediction of the water inflow according to the space-time relationship between the development height of the water flowing fractured zone and the aquifer.

[0205] 2) The mine water inflow determination device of this application can integrate the spatial characteristics of the coal mine working face to determine the changes in the permeability coefficient and thickness of the aquifer in different spaces. It can also integrate the temporal characteristics of the coal mine working face to determine the development height of the water-conducting fracture zone at different time points. Then, based on the temporal and spatial characteristics, it can make a refined prediction of the water inflow of the part to be mined in the current working face. Compared with the existing technology, which directly divides the water inflow into a whole for prediction, it is more refined and realizes the precise prediction of water inflow based on the development height of the water-conducting fracture zone and the spatiotemporal relationship of the aquifer.

[0206] The above description is merely a preferred embodiment of this application and is not intended to limit this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method of determining the inflow of water in a mine, characterised by, include: Multiple contour maps are obtained, wherein the contour maps are used to represent the relationship between the development height of the water-conducting fracture zone of the coal mine working face and the permeability coefficient and thickness of the affected aquifer. The water-conducting fracture zone is a fracture zone in which fractures connect pore water or groundwater in the rock strata, allowing water to flow into the mine through the fractures. The development height of the water-conducting fracture zone is the maximum vertical height reached by the fracture zone formed by the rock strata due to mining disturbance. The multiple contour maps include one or more of a first contour map, a second contour map, and a third contour map. The first contour map is the contour map of the mined portion of the current working face. The second contour map is the contour map of the adjacent mined working face of the current working face. The third contour map is the contour map of the unmined portion of the current working face. The development height of all the water-conducting fracture zones at multiple time points is obtained, wherein the time points and the development heights correspond one-to-one; Based on the multiple contour maps and multiple development heights, the water inflow of the section to be mined in the current working face is determined.

2. The method of claim 1, wherein, Obtain multiple contour maps, including: Obtain pumping test information of the mined portion of the current working face, wherein the pumping test information includes one or more of the following: the first aquifer layer position of the borehole pumping test, the first permeability coefficient of the first aquifer layer when pumping water from each borehole; and the first thickness of the first aquifer layer. Obtain the hardness information of the overlying rock strata of the mined portion of the current working face; Based on the hardness information, calculate the first development height of the water-conducting fracture zone in the mined portion of the current working face; Based on the pumping test information and the first development height of the water-conducting fracture zone, the first contour map is generated using ArcGIS software.

3. The method of claim 2, wherein, Obtain multiple contour maps, including: Obtain water discharge information of adjacent mined working faces of the current working face, wherein the water discharge information includes one or more of the following: the location of the outlet point of the drainage hole, the location of the final hole, the water volume, the pressure, the water discharge time, and the change in the inflow rate. Calculate the second permeability coefficient corresponding to the water discharge information according to Darcy's law; Based on the water release information and the second permeability coefficient, the second thickness of the second aquifer is determined; Obtain a second development height, wherein the second development height is the development height corresponding to the hardness information of the overlying strata of the adjacent mined working face of the current working face; Based on the second permeability coefficient, the second thickness, and the second development height, the second contour map is generated using ArcGIS software.

4. The method of claim 1, wherein, Based on multiple contour maps and multiple development heights, the water inflow of the section to be mined in the current working face is determined, including: Obtain a thematic map, wherein the thematic map is determined by filtering multiple contour maps; Based on the multiple development heights, determine all the aquifers affected by the portion to be mined in the current working face, and obtain summary information; Determine the water inflow rate according to the summary information, wherein the different aquifers have different permeability coefficients and thicknesses in the thematic map.

5. The method of claim 4, wherein, Obtain a thematic map, comprising: Compare the first contour map and the third contour map to obtain a first similarity; In a case where the first similarity is greater than or equal to a first similarity threshold, fuse the first contour map and the third contour map to obtain the thematic map, wherein the fusion method at least includes one or more of weighted average, Bayesian, Kalman filtering; In a case where the first similarity is less than the first similarity threshold, compare the second contour map and the third contour map to obtain a second similarity; In a case where the second similarity is greater than or equal to a second similarity threshold, determine the third contour map as the thematic map, wherein the first similarity threshold is greater than the second similarity threshold; In a case where the second similarity is less than the second similarity threshold, fuse the second contour map and the third contour map to obtain the thematic map.

6. The method of claim 4, wherein, According to a plurality of the development heights, determine all the aquifers that are to be mined by the current working face, to obtain summary information, comprising: Obtain a horizon relationship corresponding to the development height of each time point, wherein the horizon relationship is a relationship of the aquifers that are to be mined by the water flowing fractured zone when the water flowing fractured zone reaches the development height; According to all the horizon relationships corresponding to all the time points, determine all the aquifers that are to be mined, to obtain the summary information.

7. The method of claim 4, wherein, Determine the water inflow rate according to the summary information, comprising: According to the permeability coefficients of all the aquifers that are to be mined and the thicknesses of all the aquifers corresponding to the summary information, calculate a plurality of sub-water inflow rates by using a large well method, wherein the sub-water inflow rate and the aquifer that is to be mined are one-to-one corresponding; According to the summary information, calculate a sum of all the sub-water inflow rates to obtain the water inflow rate.

8. An apparatus for determining the inflow of a mine, characterized in that Comprise: The first acquisition unit is used for acquiring a plurality of contour maps, wherein the contour map is used for representing the relationship between the development height of the water flowing fractured zone of the coal mine working face and the permeability coefficient of the aquifer that is to be mined and the thickness of the aquifer, the water flowing fractured zone is the fractured zone that the pore water or the underground water in the fractured connected rock layer flows to the mine through the fracture, the development height of the water flowing fractured zone is the maximum vertical height of the fractured zone that reaches the rock layer due to the disturbance of mining, a plurality of the contour maps include one or more of a first contour map, a second contour map and a third contour map, the first contour map is the contour map of the mined part of the current working face, the second contour map is the contour map of the adjacent mined working face of the current working face, and the third contour map is the contour map of the to-be-mined part of the current working face; The second acquisition unit is used for acquiring the development height of all the water flowing fractured zones at a plurality of time points, wherein the time point and the development height are one-to-one corresponding; A determining unit is configured to determine the water inflow of the to-be-mined part of the current working face according to the multiple contour maps and the multiple development heights.

9. A computer program product comprising a computer program, characterized in that, The computer program, when executed by a processor, implements the steps of the method for determining the water inflow of the mine according to any one of claims 1 to 7.

10. A water inflow detection system characterized by comprising: Comprise: One or more processors, memories, and one or more programs, wherein the one or more programs are stored in the memories and configured to be executed by the one or more processors, and the one or more programs comprise programs for executing the method for determining the water inflow of the mine according to any one of claims 1 to 7.

Images