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

Abstract

The invention provides a mine water inflow determination method and device and a water inflow detection system. The method comprises the following steps: acquiring a plurality of contour maps; the development heights of all the water-conducting fractured zones at multiple time points are obtained, and the time points and the development heights are in one-to-one correspondence; and determining 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. According to the scheme, the problem that in the prior art, fine prediction of the mine water inflow cannot be achieved is solved.

Description

Technical Field

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

[0002] Mine water damage is one of the five major disasters that endanger coal mine safety production (gas, coal dust, flood, fire and roof accident). After coal mining, the surrounding rock of the working face is disturbed and damaged, and groundwater flows into the mine through the surrounding rock damage zone, forming mine water inflow. Accurately predicting the mine water inflow is crucial for mine safety production.

[0003] Current research on methods for predicting water inflow in mines often focuses on: (1) comprehensively considering the index factors that affect the water inflow in mines, using existing or mature mathematical theoretical methods, establishing mathematical models, and predicting water inflow; (2) generalizing the working face into a well or an overall range area, and using the "big well method" or "water collection corridor method" to predict water inflow. Although the above current research has achieved relatively rich research results and provided technical support for mine water hazard prevention and control to a certain extent, the current scheme cannot achieve precise prediction of mine water inflow. Summary of the invention

[0004] The main purpose of the present application is to provide a method, device, computer program product and water inflow detection system for determining the water inflow of a mine, so as to at least solve the problem that the prior art cannot achieve accurate prediction of the water inflow of a mine.

[0005] In order to achieve the above-mentioned purpose, according to one aspect of the present application, a method for determining the water inflow of a mine is provided, comprising: obtaining a plurality of contour maps, 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 of the affected aquifer and the thickness of the aquifer, the water-conducting fracture zone is a fracture zone in which the fractures connect the pore water or groundwater in the rock formation so that water flows 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 formation due to mining disturbance, and the plurality of contour maps include a first contour map; One or more of a first contour map, a second contour map and a third contour map, wherein 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 heights of all the water-conducting fracture zones at multiple time points, wherein the time points and the development heights correspond one to one; and determining 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.

[0006] Optionally, multiple contour maps are obtained, including: obtaining pumping test information of the mined portion of the current working face, wherein the pumping test information includes one or more of the layer position of the first aquifer pumped from the borehole in the pumping test, the first permeability coefficient of the first aquifer when each borehole is pumped, and the first thickness of the first aquifer; obtaining hardness information of the overlying rock strata of the mined portion of the current working face; calculating the first development height of the water-conducting fracture zone of the mined portion of the current working face based on the hardness information; and generating the first contour map using ArcGIS software based on the pumping test information and the first development height of the water-conducting fracture zone.

[0007] Optionally, multiple contour maps are obtained, including: obtaining water discharge information of the mined working face adjacent to the current working face, wherein the water discharge information includes one or more of the water outlet position, end hole position, water volume, pressure, water discharge time and water inflow change information of the drainage hole; calculating the second permeability coefficient corresponding to the water discharge information according to Darcy's law; determining the second thickness of the second aquifer according to the water discharge 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 overlying rock stratum of the mined working face adjacent to 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.

[0008] Optionally, the water yield of the to-be-mined portion of the current working face is determined based on a plurality of the contour maps and a plurality of the development heights, including: obtaining a thematic map, wherein the thematic map is determined by screening based on a plurality of the contour maps; determining all the aquifers affected by the to-be-mined portion of the current working face based on the plurality of the development heights, and obtaining summary information; determining the water yield at least based on the summary information, wherein different aquifers have different permeability coefficients and thicknesses in the thematic map.

[0009] Optionally, obtaining the thematic map includes: comparing the first contour map and the third contour map to obtain a first similarity; when 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 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 contour map and the third contour map to obtain a second similarity; when the second similarity is greater than or equal to a second similarity threshold, determining that the third contour map is the thematic map, wherein the first similarity threshold is greater than the second similarity threshold; when 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, based on multiple development heights, all the aquifers affected by the to-be-mined part of the current working face are determined to obtain summary information, including: obtaining the stratigraphic relationship corresponding to the development height at each time point, wherein the stratigraphic relationship is the relationship between the aquifers affected when the water-conducting fracture zone reaches the development height; based on all the stratigraphic relationships corresponding to all the time points, all the aquifers affected are determined to obtain the summary information.

[0011] Optionally, the water yield is determined at least based on the summary information, including: according to the permeability coefficients of all the aquifers affected and the thicknesses of all the aquifers corresponding to the summary information, a large well method is used to calculate multiple sub-water yields, wherein the sub-water yields correspond one-to-one to the aquifers affected; according to the summary information, the sum of all the sub-water yields is calculated to obtain the water yield.

[0012] According to another aspect of the present application, a device for determining the water inflow of a mine is provided, comprising: a first acquisition unit, for acquiring a plurality of contour maps, 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 of the affected aquifer and the thickness of the aquifer, the water-conducting fracture zone is a fracture zone in which the fractures connect the pore water or groundwater in the rock formation so that water flows 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 formation due to mining disturbance, and the plurality of contour maps include a first contour map, a second contour map One or more of a figure 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; a second acquisition unit is used to acquire the development heights of all the water-conducting fracture zones at multiple time points, wherein the time points and the development heights correspond one to one; a determination unit is used 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.

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

[0014] According to 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 are configured to be executed by the one or more processors, and the one or more programs include a method for determining the water inflow of any one of the mines.

[0015] By applying the technical solution of the present application, the spatial characteristics of the coal mine working face can be integrated to determine the changes in the permeability coefficient and thickness of the aquifer in different spaces. The temporal characteristics of the coal mine working face can also be integrated to determine the development height of the water-conducting fracture zone at different time points, and then the water inflow of the to-be-mined part of the current working face can be finely predicted based on the temporal and spatial characteristics. This is more refined than the method of directly dividing the work face into a whole to predict the water inflow in the prior art, and realizes the fine prediction of the water inflow based on the spatial-temporal relationship between the development height of the water-conducting fracture zone and the aquifer. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The drawings constituting part of the present application are used to provide a further understanding of the present application. The exemplary embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

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

[0018] Figure 2 A schematic flow chart of a method for determining water inflow in a mine provided in an embodiment of the present application is shown;

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

[0020] The above drawings include the following reference numerals:

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

[0022] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0023] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0024] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0025] In existing studies, for large-capacity mines, there are objective conditions of ultra-long and ultra-wide working faces, and the heterogeneity of the aquifer in the roof of the working face is particularly prominent. It is inappropriate to generalize the working face into a well or an overall range area. During coal mining, the height of the water-conducting fracture zone in the roof has temporal and spatial characteristics. As the working face advances, the height of the water-conducting fracture zone will penetrate different aquifers. It is inappropriate to only consider a single aquifer in existing studies.

[0026] As introduced in the background technology, it is impossible to achieve precise prediction of mine water inflow in the prior art. To solve the above problem, the embodiments of the present application provide a method, device, computer program product and water inflow detection system for determining the water inflow of a mine.

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

[0028] The method embodiments provided in the embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking running on a mobile terminal as an example, Figure 1 FIG. 1 is a hardware structure block diagram of a mobile terminal for determining a method for determining water inflow in a mine according to an embodiment of the present invention. Figure 1 As shown, the mobile terminal may include one or more ( Figure 1 Only one is shown in the figure) a processor 102 (the processor 102 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA) and a memory 104 for storing data, wherein the mobile terminal may also include a transmission device 106 and an input / output device 108 for communication functions. It can be understood by those skilled in the art that Figure 1 The structure shown is for illustration only and does not limit the structure of the mobile terminal. Figure 1 More or fewer components as shown, or with Figure 1 Different configurations are shown.

[0029] The memory 104 can be used to store computer programs, for example, software programs and modules of application software, such as the computer program corresponding to the display method of device information in the embodiment of the present invention. The processor 102 executes various functional applications and data processing by running the computer program stored in the memory 104, that is, the above method is implemented. The memory 104 may include a high-speed random access memory, and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include a memory remotely arranged relative to the processor 102, and these remote memories can be connected to the mobile terminal via a network. Examples of the above-mentioned 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 used to receive or send data via a network. The above-mentioned specific examples of the network may include a wireless network provided by a communication provider of the mobile terminal. In one example, the transmission device 106 includes a network adapter (Network Interface Controller, referred to as NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission device 106 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0030] In this embodiment, a method for determining the water inflow of a mine running on a mobile terminal, a computer terminal or a similar computing device is provided. 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 may be executed in an order different from that shown here.

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

[0032] Step S201, obtaining a plurality of contour maps, wherein the contour maps are used to indicate the relationship between the development height of the water-conducting fracture zone of the 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 in which the fractures connect the pore water or groundwater in the rock formation so that water flows to 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 formation due to mining disturbance, the plurality of contour maps include one or more of the first contour map, the second contour map and the 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 mined working face adjacent to 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 multiple contour maps refers to collecting and drawing a series of maps reflecting the relationship between the height of the water-conducting fracture zone in the working face of the mine, the permeability coefficient of the aquifer, and the thickness of the aquifer. A contour map is a graphical representation method that visualizes parameters that continuously change in space (such as permeability coefficient, aquifer thickness) as a series of contour lines, which is used to intuitively display the changing trends of these parameters in space.

[0034] The water-conducting fracture zone refers to the fracture zone formed during coal mining due to the collapse and cracking of the roof rock layer, which can connect the pore water or groundwater in the rock layer and allow water to flow into the mine through these fractures. The development height of the water-conducting fracture zone is a parameter that describes the maximum extension range of this fracture zone in the vertical direction, reflecting the degree of connectivity of the fracture zone formed by the rock layer due to 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 height of the water-conducting fracture zone, the permeability coefficient and thickness of the aquifer in the mined area. It is drawn based on the actual mining situation and drainage characteristic index data, and 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 height of the water-conducting fracture zone and the aquifer parameters in the area adjacent to the current working face. By comparing the characteristics of the current working face with 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 unmined part of the current working face. Based on data analysis and model prediction of the mined area, this map is used to predict the height of the water-conducting fracture zone, the permeability coefficient and thickness of the aquifer in the unmined area, providing a basis for the prediction of water inflow and risk assessment of the working face to be mined.

[0038] Step S202, obtaining the development heights of all the water-conducting fracture zones at multiple time points, wherein the time points correspond to the development heights one by one;

[0039] Specifically, "multiple time points" refer to multiple key time nodes selected in different stages of coal mining, that is, in the entire process from the beginning of mining to the end of mining. These time points may include the initial mining period, the first collapse step, the first cycle step, the time point when the maximum height of the water-conducting fracture zone is reached, and the duration after mining. Each time point corresponds to a different stage and degree of development of the water-conducting fracture zone.

[0040] Specifically, the height of the water-conducting fracture zone refers to the vertical extension of the fracture zone that can connect the pore water or groundwater in the rock layer and allow these water bodies to flow into the mine through the fractures due to the collapse and cracking of the roof rock layer during coal mining. The size of this height parameter directly affects the size of the mine water inflow and is a key indicator for assessing the risk of mine water hazards.

[0041] Specifically, the height of the water-conducting fracture zone changes with the progress of mining time. This change process reflects the dynamic development of the fracture zone formed by the mining disturbance, including the formation, expansion, reaching the maximum height and possible gradual closure of the fracture zone. By monitoring and recording the height of the water-conducting fracture zone at each key time point, the relationship between time and the height of the fracture zone can be established, providing basic data for predicting the amount of water inflow in the mine. Obtaining the height data of the water-conducting fracture zone at multiple time points is crucial to understanding the spatiotemporal evolution of the fracture zone, evaluating the changing trend of the amount of water inflow in the mine, and formulating reasonable water hazard prevention and control measures.

[0042] Step S203, determining the water inflow of the to-be-mined portion of the current working face based on the plurality of contour maps and the plurality of development heights.

[0043] Specifically, the contour map contains the distribution of key hydrogeological parameters such as the permeability coefficient and aquifer thickness of the aquifer affected by the maximum development height of the water-conducting fracture zone. By analyzing these contour maps, we can understand the spatial variation of aquifer parameters and the correlation between these parameters and the development height of the water-conducting fracture zone. The development height data records the vertical extension range of the water-conducting fracture zone at different mining stages. These data correspond to time points one by one and represent the dynamic changes of the fracture zone with the mining process. When predicting the water inflow of the part to be mined, it is necessary to combine the development height of the water-conducting fracture zone recorded at multiple key time points to understand the extent to which the fracture zone reaches the aquifer, and then evaluate the connectivity between the fracture zone and the aquifer and the water flow path.

[0044] Specifically, water yield prediction is not just a static analysis of aquifer parameters, but requires a comprehensive consideration of the spatiotemporal characteristics of fracture zone development. This means that the prediction model needs to consider the expansion of fracture zones over time and the impact of this expansion on the distribution of aquifer parameters. By combining contour maps with development height data, it is possible to more accurately determine which aquifers will be affected by fracture zones at different stages of exploitation, and how aquifer parameters affect water yield.

[0045] By using the above analysis results, combined with specific hydraulic models (such as the "big well method") and geological conditions, the water inflow of the working face to be mined can be accurately predicted at different time points and in different aquifer areas. This prediction result can provide real-time and dynamic guidance for mine water hazard prevention and control, help mine managers plan water drainage measures in advance, reduce water hazard risks during mining, and ensure safe production in mines.

[0046] Through this embodiment, the spatial characteristics of the coal mine working face can be integrated to determine the changes in the permeability coefficient and thickness of the aquifer in different spaces. The temporal characteristics of the coal mine working face can also be integrated to determine the development height of the water-conducting fracture zone at different time points, and then the water inflow of the to-be-mined part of the current working face can be finely predicted based on the temporal and spatial characteristics. This is more refined than the method of directly dividing the working face into a whole to predict the water inflow in the prior art, and realizes the fine prediction of the water inflow based on the temporal and spatial relationship between the development height of the water-conducting fracture zone and the aquifer.

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

[0048] In the specific implementation process, obtaining multiple contour maps can be achieved through the following steps: obtaining pumping test information of the mined part of the above-mentioned current working face, wherein the above-mentioned pumping test information includes one or more of the layer position of the first aquifer pumped by the borehole in the pumping test, the first permeability coefficient of the above-mentioned first aquifer when each borehole is pumped, and the first thickness of the above-mentioned first aquifer; obtaining the hardness information of the overlying rock strata of the mined part of the above-mentioned current working face; calculating the first development height of the above-mentioned water-conducting fracture zone of the mined part of the above-mentioned current working face based on the above-mentioned hardness information; and generating the above-mentioned first contour map using ArcGIS software based on the above-mentioned pumping test information and the above-mentioned first development height of the above-mentioned water-conducting fracture zone.

[0049] In this solution, the first contour map is generated based on the actual data of the mined working face and GIS technology, which can provide a direct relationship between the development height of the water-conducting fracture zone and the hydrogeological parameters of the aquifer, so that the water inflow can be predicted based on the accurate first contour map in the future, which helps to improve the efficiency of mine water hazard prevention and control.

[0050] Specifically, data such as pumping test data of the mining area where the working face to be mined is located can be obtained. The above-mentioned acquisition of pumping test data of the mining area where the working face to be mined is located specifically refers to: the position of the aquifer layer pumped by the pumping test borehole within the mining area where the working face to be mined is located, the permeability coefficient of each borehole pumping aquifer, the thickness of the aquifer and other data.

[0051] Specifically, a contour map of permeability coefficient and aquifer thickness of one or more aquifers affected by the maximum development height of the water-conducting fracture zone within the mining area where the working face to be mined is located, that is, a first contour map, can be drawn.

[0052] Specifically, the above-mentioned 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 within the mining area where the working face to be mined is located, the specific steps are as follows:

[0053] (1) According to the "Specifications for the Establishment of Coal Pillars and Compressed Coal Mining in Buildings, Water Bodies, Railways and Main Wells and Lanes", the maximum development height of the water-conducting fracture zone in coal seam mining is estimated (the + sign is used, and the positive value is used in the calculation). Specifically, the maximum development height of the water-conducting fracture zone is the maximum distance that the roof rock cracks extend upward after coal seam mining. This height determines whether the water in the roof aquifer will flow into the mine. When estimating this height, the following formula is usually used for calculation:

[0054] When the overlying rock type of the coal seam roof is hard:

[0055]

[0056] When the overlying rock type of the coal seam roof is medium hard:

[0057]

[0058] When the overlying rock type of the coal seam roof is weak:

[0059]

[0060] Where: H Li ——Height of water-conducting fracture zone (m); ∑M——Cumulative mining thickness (m).

[0061] (2) Based on the predicted maximum height of the water-conducting fracture zone, combined with the stratigraphic structural characteristics and the characteristics of the aquiclude, determine one or more aquifers within the maximum height of the water-conducting fracture zone within the mining area; that is, once the maximum height of the water-conducting fracture zone is predicted, combined with the stratigraphic structural characteristics and the characteristics of the aquiclude, it can be determined which aquifers within this height range may be affected. The stratigraphic structural characteristics include the composition, sequence and thickness of the strata, while the characteristics of the aquiclude mainly focus on the location, thickness and water-proof performance of the aquiclude. Through this analysis, it is possible to accurately locate which aquifers may have water entering the mine due to mining activities.

[0062] Specifically, the stratigraphic structural characteristics include the depth of the coal seam and the properties of the overburden (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-conducting fracture zone can be determined. Then, the lithology of the overburden, such as sandstone and limestone, is analyzed. The fracture development and permeability of these rock layers are key factors affecting the development of the water-conducting fracture zone.

[0063] Specifically, an aquiclude refers to a rock layer between aquifers or between an aquifer and a non-aquifer that can block or slow the flow of water. Identifying the location and thickness of an aquiclude is crucial for determining whether a water-conducting fracture zone can penetrate to the next aquifer and predicting the water flow path. The thicker the aquiclude, the better the water-blocking performance, and the less likely it is that a water-conducting fracture zone will penetrate it.

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

[0065] (3) Use ArcGIS to draw contour maps of the permeability coefficient and thickness of one or more aquifers affected by the maximum development height of the water-conducting fracture zone within the mining area where the working face to be mined is located. That is, after determining the aquifers that may be affected by the water-conducting fracture zone, the next step is to draw contour maps of the permeability coefficient and thickness of these aquifers. The permeability coefficient is a parameter that measures the water permeability of the aquifer, and the thickness of the aquifer directly affects the water storage capacity of the aquifer. Using ArcGIS software, contour maps of these parameters can be drawn through spatial interpolation and other methods based on existing pumping test data and borehole drainage characteristic indicators. The contour map can intuitively show the distribution law of the parameters in space, which helps to identify the high permeability and thick layer areas of the aquifer, which are potential high-risk areas for 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 the overburden has a direct impact on the formation and expansion of the water-conducting fracture zone. Hard rock formations may inhibit the formation of cracks, while weak rock formations are prone to cracks. Based on the hardness information of the overburden, the maximum vertical extension range (first development height) of the water-conducting fracture zone during mining can be predicted.

[0067] Specifically, ArcGIS is a powerful geographic information system software that can process and analyze spatial data and generate various types of contour maps, distribution maps, etc. In this process, the pumping test information (including the layer position, permeability and thickness of the aquifer) 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 permeability and thickness of the aquifer within the development height range of the water-conducting fracture zone in the mined working face area, providing key visualization information for subsequent water inflow prediction.

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

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

[0070] Specifically, data such as the characteristic indicators of water drainage and drainage of the drilled holes and the water inflow of the mined working face adjacent to the working face to be mined can be obtained. The above-mentioned acquisition of the characteristic indicators of water drainage and drainage and the water inflow of the drilled holes adjacent to the mined working face to be mined specifically refers to: the elevation angle of the water outlet point of the drainage hole, the terminal hole position, the real-time monitored water volume and pressure data and the corresponding water discharge time, and the law of change of water inflow during the advancement of the mined working face adjacent to the working face to be mined.

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

[0072] Specifically, the above-mentioned 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 mined working face adjacent to the working face to be mined is carried out in the following specific steps:

[0073] 1. According to the above formulas 1, 2 and 3, the maximum development height of the water-conducting fracture zone in the mined working face is estimated.

[0074] 2. According to the characteristic indicators of the drilling drainage of the adjacent mined working face to be mined, the permeability coefficient of the drainage hole position can be inferred according to Darcy's law. The permeability coefficient distribution map of the entire area can be obtained by calculating the permeability coefficients of multiple drainage holes and drawing using ArcGIS interpolation method. That is, using Darcy's law, the permeability coefficient of the drainage hole position can be inferred based on the drilling drainage test data of the mined working face. Darcy's law describes that under laminar flow conditions, the flow velocity of water in porous media is linearly related to the hydraulic gradient, among which the permeability coefficient is the key parameter, reflecting the permeability of the aquifer. By analyzing the real-time monitoring data of multiple drainage holes, the permeability coefficient of each hole can be calculated, and then the interpolation function of ArcGIS software can be used to convert these point data into a continuous permeability coefficient distribution map, which intuitively displays the permeability coefficient distribution of the entire mined working face area.

[0075] Specifically, according to Darcy's law, the groundwater flow rate is proportional to the 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 outflow per unit time, K represents the permeability coefficient, A represents the cross-sectional area, h1 and h2 represent the groundwater level, and L represents the flow distance.

[0076] Specifically, after knowing the location of the water discharge hole, the real-time monitored water volume and pressure data, and the corresponding water discharge time, the change in groundwater level can be calculated through the monitoring data, thereby obtaining the groundwater flow velocity. Combined with the groundwater flow distance, the permeability coefficient K can be inferred.

[0077] Specifically, by calculating the permeability coefficients of multiple drainage holes, the permeability coefficient data at different locations can be obtained. Using 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 using known point data to infer unknown point values. Common interpolation methods include Kriging interpolation, inverse distance weighted interpolation, etc. In ArcGIS software, you can select a suitable interpolation method, input the known permeability coefficient data points into the software, perform interpolation processing, and generate a permeability coefficient distribution map.

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

[0080] 3. Combine the real-time monitored flow and water pressure data to analyze the changes in water volume during the water release process. Based on these changes, further calibrate the thickness of the aquifer. Finally, use the groundwater flow model to fit and analyze the water release data, and combine it with the known permeability coefficient to obtain the thickness information of the aquifer through inversion calculation. That is, by analyzing the real-time monitored flow and water pressure data, it is possible to observe the change pattern of water volume over time during the water release process, which helps to calibrate the thickness of the aquifer. The thickness of the aquifer is closely related to the amount of water inflow, and thicker aquifers can store and release more water. Use a groundwater flow model (such as a numerical model or an analytical model) to fit and analyze the water release data, and combine it with the known permeability coefficient to invert and calculate the thickness information of the aquifer. This process relies on the optimization of hydraulic principles and model parameters, aiming to improve the accuracy of the aquifer model.

[0081] Specifically, the real-time monitoring of flow and water pressure data can be used to determine the changes in water volume during the release process. When the release begins, the flow rate will gradually increase and the water level will rise accordingly; when the release ends, the flow rate will gradually decrease and the water level will gradually drop. By analyzing the changing trends of flow and water pressure data, the changes in water volume in the aquifer can be inferred.

[0082] Specifically, calibrating the thickness of the aquifer requires combining real-time monitored water level data with a groundwater flow model. Based on the changes in the water level data, it can be inferred whether the thickness of the aquifer is consistent with expectations. If the actual water level changes are inconsistent with the water level changes predicted by the model, the aquifer thickness parameters may need to be adjusted to make the model more consistent with the actual situation.

[0083] Specifically, the water release data is fitted and analyzed by the groundwater flow model, and combined with the known permeability coefficient, the thickness information of the aquifer can be inverted and calculated. By comparing the actual observed flow, water pressure and water level data with the data predicted by the model, the optimal aquifer thickness parameters can be determined, thereby obtaining the aquifer thickness information that best meets the actual situation.

[0084] 4. Combined with the water inflow law during the advancement of the mined working face and the contour lines of the permeability coefficient and aquifer thickness within the mining area, verify and analyze the rationality of the distribution of the permeability coefficient and aquifer thickness within the mined working face. That is, the permeability coefficient and aquifer thickness distribution calculated above are compared with the water inflow law actually monitored during the advancement of the mined working face, and compared with the contour lines of the permeability coefficient and aquifer thickness within the mining area. The purpose of this comparative analysis is to verify whether the prediction results of the model are reasonable and whether they can accurately reflect the hydrogeological characteristics of the aquifer of the mined working face. If the prediction results are consistent with the actual monitoring data and the parameter distribution within the mining area, it means that the establishment of the model is reasonable; if not, the model parameters need to be adjusted until the prediction results are consistent with the actual data.

[0085] Specifically, a curve is drawn showing the change of water inflow from the mined working face with the mining progress, observing the increase and decrease trend of water inflow at different mining stages, and whether there is an obvious turning point or peak. The contour map of permeability coefficient and aquifer thickness within the mining area is analyzed to determine the parameter distribution characteristics of the area where the mined working face is located, such as high permeability area, thick aquifer area, etc.

[0086] Specifically, based on the water inflow data of the mined working face, as well as the permeability and thickness data of the aquifer, a hydraulic model (such as the "big well method") is used to simulate the change in water inflow during the mining process of the working face. The model should take into account the development of the water-conducting fracture zone 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 consistency between the model prediction results 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 such as the permeability and aquifer thickness used in the model to make it more consistent with the actual geological conditions. This usually involves a sensitivity analysis of the model parameters to determine which parameters have the greatest impact on the results. After repeated adjustments and verifications, if the model can accurately predict the trend of water inflow, then it can be considered that the distribution of permeability and aquifer thickness within the mining area is reasonable. In addition, the key factors of water inflow changes, such as the height of the water-conducting fracture zone and the heterogeneity of the aquifer, can be further analyzed.

[0087] Specifically, in the process of coal mining, in order to control the amount of water inflow after the working face is mined, drainage measures are usually implemented on the mined working face, including drilling drainage. The drainage information includes the water outlet position, the end hole position, the real-time monitored water volume, pressure, drainage time, and the change law of water inflow. The water outlet position and the end hole position reflect the layout of the drainage holes and whether they effectively penetrate the aquifer; the water volume and pressure data reveal the water storage capacity and water conductivity of the aquifer; the drainage time and water inflow change information reflect the dynamic interaction between the drainage holes and the aquifer, as well as the response characteristics of the aquifer.

[0088] Specifically, Darcy's law is the basic principle for describing fluid flow in porous media. It states that the flow rate of water under laminar flow is proportional to the head gradient and the permeability coefficient of the medium. By analyzing the water volume, pressure and water discharge time data of the drainage hole, the principle of Darcy's law can be applied to reversely calculate the permeability coefficient of the corresponding aquifer, that is, the second permeability coefficient. This calculation process requires combining the measured water discharge volume, water pressure and time data with the geometric parameters of the aquifer (such as porosity, effective thickness, etc.) and solving them using numerical methods or theoretical models.

[0089] Specifically, the thickness of the aquifer is another key parameter that affects the amount of water inflow. Based on the analysis of the water release 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 changes in water volume at the water release hole at different time points, the water storage capacity and water flow path of the aquifer can be inferred, and then the thickness of the aquifer can be estimated. This process may require the use of a groundwater flow model, 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 formation in the area through formulas 1, 2, and 3. The second development height reflects the degree of rock formation crack formation caused by mining disturbance and is an important indicator for evaluating the amount of water inflow in a mine.

[0091] Specifically, the data of the second permeability coefficient, the second thickness and the second development height are imported into ArcGIS software, and its spatial analysis and interpolation functions are used to generate a second contour map reflecting the changing trend of the hydrogeological parameters of the second aquifer. The specific steps include: Data import: Import the data of the second permeability coefficient, the second thickness and the second development height into ArcGIS in the form of spatial distribution. Spatial interpolation: Use ArcGIS's interpolation tools, 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 area. Contour map generation: Based on the interpolation results, draw the contour map of the second permeability coefficient and the second thickness. This figure intuitively shows the parameter distribution characteristics of the second aquifer in the adjacent mined working face.

[0092] In addition, by obtaining the third contour map among multiple contour maps, data such as mining conditions, geological conditions, hydrogeological conditions, and borehole drainage characteristics can be obtained. The above-mentioned data such as mining conditions, geological conditions, hydrogeological conditions, and borehole drainage characteristics of the working face to be mined specifically refer to: mining conditions of the working face to be mined: mining method, coal seam mining thickness, working face strike length, and working face dip length; geological conditions: stratigraphic structural characteristics, overburden type; hydrogeological conditions: aquiclude characteristics, roof aquifer hydrogeological parameters; borehole drainage characteristics: drainage hole elevation water outlet position, terminal hole position, real-time monitored water volume and pressure data, and corresponding drainage time, etc.

[0093] Specifically, a contour map of permeability coefficient and aquifer thickness of one or more aquifers affected by the maximum development height of the water-conducting fracture zone of the working face to be mined, ie, a third contour map, can be drawn.

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

[0095] Specifically, ensure that all relevant data such as mining conditions, geological conditions, hydrogeological conditions, and borehole drainage characteristics are collected. These data include the mining method of the working face, the thickness of the coal seam, the size of the working face, the stratigraphic structure, the overburden type, the characteristics of the aquifer (permeability, thickness, etc.), and the measured data of the borehole drainage (water volume, water pressure, drainage time, etc.). Based on the real-time monitoring data of the drainage holes of the mined working face, the permeability coefficient of the drainage hole position is reversed using Darcy's law. This step requires substituting the water volume, water pressure changes, and hole depth of the borehole into the formula of Darcy's law to invert and calculate the permeability coefficient. The permeability coefficient data of each drainage hole position calculated by Darcy's law is imported into GIS software, such as ArcGIS. In the software, an appropriate interpolation method (such as inverse distance weighted interpolation IDW, Kriging interpolation, etc.) is used to generate a distribution map of the permeability coefficient of the entire mined working face area. The selection of interpolation parameters depends on the distribution of data points and the requirements for prediction accuracy. By analyzing the changes in water volume and water pressure during the process of drilling and draining water, the thickness information of the aquifer can be inferred. These thickness data are also imported into GIS, and the distribution map of the aquifer thickness is drawn using the interpolation method. Using the prediction formula mentioned above (based on the overburden type), the maximum development height of the water-conducting fracture zone is determined, and then the distribution of the aquifer within this height range is analyzed. This requires combining the drawn permeability coefficient and aquifer thickness distribution map with the water-conducting fracture zone development height information, and analyzing the change law of the aquifer parameters within this height range to 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 map of the permeability coefficient and aquifer thickness to check whether the model prediction results are consistent with the actual water inflow data. 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 working face to be mined. In the prediction process, it is necessary to consider the mining conditions, geological characteristics and possible development of the water-conducting fracture zone of the working face to be mined, and combine the drawn parameter distribution map to make a refined prediction of the water inflow.

[0096] In the specific implementation process, the water yield of the to-be-mined part of the above-mentioned current working face is determined according to the multiple above-mentioned contour maps and the multiple above-mentioned development heights, which can be achieved through the following steps: obtaining a thematic map, wherein the above-mentioned thematic map is screened and determined based on the multiple above-mentioned contour maps; determining all the above-mentioned aquifers affected by the to-be-mined part of the above-mentioned current working face according to the multiple above-mentioned development heights, and obtaining summary information; determining the above-mentioned water yield at least based on the above-mentioned summary information, wherein different aquifers have different above-mentioned permeability coefficients and above-mentioned thicknesses in the above-mentioned thematic map.

[0097] In this scheme, multiple contour maps and the development height of water-conducting fracture zones are used to determine the aquifers that may be affected by the part to be mined at the current working face, and ultimately the water inflow is calculated, which can further improve the accuracy of water inflow prediction.

[0098] Specifically, the permeability coefficient and thickness of one or more aquifers on the roof of the working face to be mined can be mapped; the development law of the height of the water-conducting fracture zone as the working face to be mined advances can be simulated; and the water inflow from the roof can be precisely predicted by zoning and stratification during the advancement of the working face.

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

[0100] Specifically, based on the development height data of multiple key mining stages of the current working face (such as initial collapse, periodic collapse, reaching the maximum development height of the water-conducting fracture zone, etc.), all aquifers that may 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 the water inflow. Since the development height of the water-conducting fracture zone varies with mining depth and mining method, it is necessary to update this information in real time in combination with the specific mining conditions and mining progress of the working face.

[0101] Specifically, summary information refers to a collection of information that organizes and summarizes the parameters (permeability and thickness) of all aquifers affected by the current working face to be mined. This information set includes detailed parameters of each aquifer, as well as their specific location and distribution in the thematic map, which provides an important basis for the subsequent water yield prediction. The acquisition of summary information may involve statistical analysis of the parameters of each aquifer, identification of the law of parameter changes, and correlation with the height of the development of the water-conducting fracture zone.

[0102] Specifically, the calculation of water inflow is based on the permeability and thickness data of different aquifers in the thematic map. Since the parameters of different aquifers may vary greatly, the uniqueness of each aquifer needs to be taken into account when predicting water inflow, that is, the stronger the permeability and the thicker the aquifer, the greater the water inflow it contributes. This step may adopt the "big well method", "water collection gallery method" or other hydrogeological models, combined with the current working face mining conditions and the height of the water-conducting fracture zone, to calculate the water inflow of each aquifer in layers and zones. Finally, by integrating the calculation results of the water inflow 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, 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. These maps reflect the distribution of permeability and thickness of aquifers in different areas. Then, the height of the water-conducting fracture zone development at the key stages of the first collapse and the first cycle collapse of working face B is analyzed, and the aquifers that may be affected in these stages are determined. Then, the data of the parameters of these aquifers, including permeability and thickness, and their positions in the thematic map are sorted out to form summary information.

[0104] Based on the summary information, the "big well method" model was applied to perform a detailed calculation of the water inflow of each aquifer. In this process, the differences in parameters of different aquifers and their interaction with the height of the development of water-conducting fracture zones were fully considered to ensure the accuracy and reliability of the prediction results. Through model calculations, the predicted values ​​of water inflow that working face B may encounter in key stages such as the initial collapse stage and the first cycle collapse stage, as well as the total water inflow prediction during the entire mining cycle were obtained.

[0105] In some embodiments, obtaining the thematic map can be specifically achieved through the following steps: comparing the first contour map and the third contour map to obtain a first similarity; when 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 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 contour map and the third contour map to obtain a second similarity; when the second similarity is greater than or equal to a second similarity threshold, determining that the third contour map is the thematic map, wherein the first similarity threshold is greater than the second similarity threshold; when 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.

[0106] In this scheme, a more accurate thematic map can be obtained by merging and screening multiple contour maps, which can improve the accuracy of subsequent water inflow prediction.

[0107] Specifically, the above-mentioned zoning of the permeability coefficient and thickness of one or more aquifers on the top plate of the working face to be mined is carried out in the following specific steps:

[0108] ①. Combined with the thematic maps of permeability coefficient and thickness of one or more aquifers in the mining area, adjacent working faces, and the top plate of the working face to be mined, focus on comparing and analyzing the changing patterns of the permeability coefficient and thickness contour lines of the aquifer in the working face to be mined and the permeability coefficient and thickness contour lines of the aquifer in the area where the working face to be mined is located within the mining area.

[0109] ②. If the permeability coefficient and thickness contour lines of the aquifer at the working face to be mined are basically consistent with the permeability coefficient and thickness contour lines of the aquifer in the area where the working face to be mined is located within the mining area, take the average value of the two as the thematic map of the permeability coefficient and thickness contour lines of the top plate aquifer of the working face to be mined.

[0110] ③. If the permeability coefficient and thickness contour lines of the aquifer at the working face to be mined are significantly different from or inconsistent with the permeability coefficient and thickness contour lines of the aquifer in the area where the working face to be mined is located within the mining area; and then further compare the permeability coefficient and thickness contour lines of the aquifer at the working face to be mined with the permeability coefficient and thickness contour lines of the adjacent mined working face, and by analyzing the extension trend of the permeability coefficient and thickness contour lines of the adjacent mined working face toward the working face to be mined, which is consistent with the trend of the permeability coefficient and thickness contour lines of the aquifer at the working face to be mined, then the thematic map of the permeability coefficient and thickness contour lines of the top plate aquifer at the working face to be mined is selected and adopted.

[0111] ④. Combined with the natural fracture level zoning in ArcGIS spatial analysis, the thematic maps of the permeability coefficient and thickness contours of the aquifer on the top of the working face to be mined are divided into four level areas: I, II, III and IV. Among them, level I means that the permeability coefficient / thickness of the aquifer is large.

[0112] Specifically, the above-mentioned simulation of the development law of the height of the water-conducting fracture zone as the working face to be mined advances is carried out in the following specific steps:

[0113] I. According to the stratum structural characteristics of the working face to be mined, combined with 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, simulate the step-by-step excavation of the coal seam to obtain the height of the water-conducting fracture zone at key nodes such as the initial collapse step, the first cycle step, and the initial advancement distance to reach the maximum height of the water-conducting fracture zone;

[0115] III. Summarize the evolution law of the development height of water-conducting fracture zones during the mining process of the coal seam to be mined.

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

[0117] Specifically, if the first similarity reaches or exceeds the preset first similarity threshold, it means that the hydrogeological parameters of the mined working face can represent or predict the situation of the working face to be mined to a certain extent. At this time, the weighted average method can be used to merge the data of the parameters 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 the distribution of aquifer parameters based on prior knowledge and observational data. In the prediction of water inflow, 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 can be regarded as the "prediction value". Through the Kalman filtering algorithm, the information of the two can be integrated to obtain a more accurate estimation of the distribution of aquifer parameters, 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 working face may be quite different from those of the working face to be mined, and data fusion cannot be performed directly. At this time, the similarity (second similarity) between the second contour map (aquifer parameter distribution map of the adjacent mined working face) and the third contour map will be compared.

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

[0122] Specifically, if the second similarity reaches or exceeds the second similarity threshold, the third contour map can be used directly as a thematic map without the need for an additional fusion step. This indicates that the distribution of aquifer parameters of the working face to be mined is relatively consistent with that of the adjacent mined working face, and 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, this indicates that there are significant differences in the hydrogeological conditions of the working face to be mined and those of the mined working face. At this time, it is necessary to merge the second and third contour maps to generate a thematic map that integrates the information of the two for further water inflow prediction.

[0124] For example, in practical applications, it is assumed that the water inflow of working face C is predicted, and working faces A and B are its mined and adjacent mined working faces. First, the contour maps of the aquifer parameter distribution of working faces A, B, and C were obtained respectively. It was found that there were large differences in the permeability coefficient and thickness distribution between the first contour map of working face A and the third contour map of working face C, and the first similarity was low, failing to reach the first similarity threshold. However, the second contour map of working face B and the third contour map of working face C showed a high similarity in some parameters, exceeding the lower second similarity threshold. Therefore, it was 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 faces B and 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, then the parameters of working face B may 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 the water inflow prediction of working face C. Through this strategy based on similarity analysis and data fusion, mine managers can better understand the hydrogeological differences between different working faces, and adjust the drainage measures accordingly to ensure the safe production of the mine and the effective management of water resources.

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

[0127] In some embodiments, based on the multiple development heights, all the aquifers affected by the to-be-mined part of the current working face are determined to obtain summary information, which can be specifically achieved through the following steps: obtaining the stratigraphic relationship corresponding to the development height at each of the above-mentioned time points, wherein the stratigraphic relationship is the relationship between the aquifers affected when the water-conducting fracture zone reaches the above-mentioned development height; determining all the aquifers affected based on all the stratigraphic relationships corresponding to all the above-mentioned 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 impact of the development of water-conducting fracture zones on different aquifers during the exploitation process, the stratigraphic relationship corresponding to the development height at each time point can be obtained, and all the aquifers affected can be determined based on this information, and summary information can be generated. In this way, all the aquifers affected by the exploitation process can be determined, and then the water inflow can be finely predicted later.

[0129] Specifically, during coal mining, as the working face advances, the rock strata above the coal seam will be subject to different mining disturbances, causing the development height of the water-conducting fracture zone to change over time. These changes will directly affect whether the water-conducting fracture zone is affected and which aquifers are affected. Therefore, for each key time point in the mining process, we need to obtain the development height of the water-conducting fracture zone and determine the aquifer affected 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 aquifer at a specific time point, including whether the fracture zone penetrates the aquifer, which aquifers are affected, and its relative position to the aquifer.

[0130] Specifically, after analyzing the stratigraphic relationships at multiple key time points during the mining process, it is necessary to organize and summarize this information to determine all aquifers affected by the water-conducting fracture zone during the entire mining cycle. This includes aquifers that may be encountered at different stages such as the initial mining, initial collapse, periodic collapse, and the later stages of mining. Determine the affected aquifers: By comparing the development heights at different time points with the stratigraphic positions of the aquifers, identify the specific aquifers affected by the water-conducting fracture zone at each mining stage, as well as the changes in the parameters of these aquifers, such as permeability and thickness.

[0131] Specifically, summary information refers to a collection of information on the hydrogeological parameters of all affected aquifers (such as permeability, thickness, etc.) and their spatiotemporal changes during the mining process. This information set is crucial for water inflow prediction because it directly reflects the water inflow risk that may be encountered after the working face is mined. Generate summary information: Convert the stratigraphic relationship at each time point into specific aquifer parameter data, and combine it with the development law of the water-conducting fracture zone to generate a summary information list containing all affected aquifer parameters and their change trends.

[0132] In some embodiments, the water yield is determined at least based on the summary information, which can be specifically achieved through the following steps: according to the permeability coefficients of all the aquifers affected and the thicknesses of all the aquifers corresponding to the summary information, a large well method is used to calculate a plurality of sub-water yields, wherein the sub-water yields correspond one-to-one to the aquifers affected; according to the summary information, the sum of all the sub-water yields is calculated to obtain the water yield.

[0133] In this scheme, the sum of all sub-water inflows can be counted, which not only can obtain the sub-water inflow prediction of the working face at different mining stages, but also the overall water inflow. Since the temporal and spatial characteristics of the mining process are taken into consideration, the accuracy of the prediction results is further improved.

[0134] Specifically, the water inflow from the roof during the advancement of the above-mentioned working face is precisely estimated by zoning and layering. The specific steps are as follows:

[0135] (i) Analyze the relationship between the height of the water-conducting fracture zone and the stratum of the roof aquifer during the advancement of the working face to be mined, focusing on the relationship between the height of the water-conducting fracture zone corresponding to the initial collapse step and the first cycle step, the initial maximum height of the water-conducting fracture zone development and the maximum continuous height of the water-conducting fracture zone (4 periods) and the stratum of one or more aquifers in the roof, and determine the aquifers affected by the height range of the water-conducting fracture zone in the 4 periods;

[0136] 2. Combined with the four zones of the zoning level of the roof aquifer of the working face to be mined, the “big well method” is used to calculate the roof water inflow in four periods, different roof aquifer levels, and different roof aquifer areas during the advancement of the working face;

[0137] 3. Statistically analyze the evolution law of roof water inflow during the advancement of the working face; and after statistically integrating the calculation results of the water inflow of each section, the water inflow of the entire working face can be obtained.

[0138] Specifically, the big well method is a commonly used hydrogeological model to estimate the inflow of groundwater in a certain area. In the prediction of water inflow from the mine roof, the water-conducting fracture zone formed after the working face is mined can be regarded as a large "well". The water inflow risk of the entire working face can be evaluated by calculating the sub-water inflow of each aquifer affected by the fracture zone. Sub-water inflow: refers to the water inflow 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, the sub-water yield is calculated using the large well method model based on its permeability coefficient K, thickness T, and the spatiotemporal characteristics of the water-conducting fracture zone affecting the layer. This usually involves using the groundwater flow equation, combined with the specific parameters of the aquifer, to estimate the amount of water delivered by the aquifer to the working face at different stages of exploitation.

[0140] Specifically, after calculating the sub-inflows of all affected aquifers, these sub-inflows are summed up to obtain the total inflow of the entire working face at a specific mining stage. This process ensures that the prediction results take into account all factors that may affect the inflow, improving the comprehensiveness and accuracy of the prediction. Inflow: It is the total amount of water transported from the roof aquifer to the well after the working face is mined. It is the accumulation of the sub-inflows of all affected aquifers.

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

[0142] Case: Prediction of water inflow at working face E in the mine.

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

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

[0145] The first cycle collapse stage: With the further advancement of the working face, the height of the water-conducting fracture zone increases to H2, at which time it affects the aquifers L1 and L2.

[0146] Stable mining stage: At a certain stable stage in the mining process, the water-conducting fracture zone develops highly stably at H3, which then affects the aquifers L2, L3 and L4.

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

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

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

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

[0151] During the stable exploitation stage, the water-conducting fracture zone development height H3 affects the 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), as well as the parameters of L4 (such as permeability coefficient K4 and thickness T4) are recorded.

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

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

[0154] The information on the stratigraphic relationships of the initial collapse, the first cycle collapse, the stable mining phase, and the final collapse phase are integrated to form a comprehensive summary information list. This list not only includes all aquifers (L1, L2, L3, L4, L5) affected by the water-conducting fracture zone during the mining process, but also records in detail the changes in parameters of each aquifer at different mining stages, such as permeability and thickness.

[0155] A summary listing might look like this:

[0156] Aquifer L1: Affected in the initial collapse stage, permeability K1, thickness T1.

[0157] Aquifer L2: Affected by the initial collapse and the first cycle collapse, with a permeability coefficient of K2 (which may vary in different stages) and a thickness of T2.

[0158] Aquifer L3: Affected during the stable exploitation stage, with a permeability coefficient of K3 and a thickness of T3.

[0159] Aquifer L4: Affected during the first cycle collapse, stable mining stage and final collapse stage, with a permeability coefficient of K4 (which may vary in different stages) and a thickness of T4.

[0160] Aquifer L5: Affected in the final collapse stage, permeability K5, thickness T5.

[0161] Step D: forecasting water inflow based on the summary information;

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

[0163] Assume that the prediction is the water inflow from the roof of working face E. In the initial collapse stage, the water-conducting fracture zone affected the aquifer L1, whose parameters were permeability coefficient K1 = 1x10^-5m / s and thickness T1 = 20m; in the first cycle collapse stage, the aquifers L1 and L2 were affected, of which the parameters of L2 were permeability coefficient K2 = 2x10^-5m / s and thickness T2 = 25m; in the stable mining stage, the aquifers L2, L3 and L4 were affected, and the parameters of L3 and L4 were permeability coefficient K3 = 1.5x10^-5m / s, thickness T3 = 30m and permeability coefficient K4 = 1.8x10^-5m / s, thickness T4 = 35m; in the final collapse stage, the aquifers L4 and L5 were affected, and the parameters of L5 were permeability coefficient K5 = 2.2x10^-5m / s and thickness T5 = 40m.

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

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

[0166] In this case, by carefully analyzing the relationship between the height of the water-conducting fracture zone and the aquifer sweep, and integrating the aquifer parameter information at the corresponding time point, a detailed, phased and layered parameter basis can be provided for the prediction of the roof water inflow of the E working face of the mine. This prediction method not only takes into account the dynamic changes of the water-conducting fracture zone during the mining process, but also fully evaluates the contribution of different aquifers to the water inflow, thereby improving the accuracy of the prediction and providing a scientific basis for mine water hazard prevention and water inflow management.

[0167] In summary, this scheme is aimed at large-capacity mines with the objective situation of super-long and super-wide working faces, and the heterogeneity of the aquifer on the working face roof is particularly prominent. The present invention analyzes the hierarchical zoning of hydrogeological parameters such as permeability coefficient and thickness of one or more aquifers on the roof from the mining area where the working face to be mined is located to the working face; the height of the water-conducting fracture zone in the roof has temporal and spatial characteristics during coal seam mining, and comprehensively considers the water-conducting fracture zone laws in four time periods, including the initial collapse step, the height of the water-conducting fracture zone corresponding to the first cycle step, the initial maximum height of the water-conducting fracture zone development, and the maximum continuous height of the water-conducting fracture zone, and analyzes its temporal and spatial relationship with the roof aquifer, thereby realizing the refined prediction of the roof water inflow in the time dimension, segmentation dimension and layering dimension of the working face to be mined.

[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 in the embodiment of the present application. The device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware for a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived.

[0169] The following is an introduction to the device for determining the water inflow of a mine provided in an embodiment of the present application.

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

[0171] The first acquisition unit 10 is used to acquire a plurality of contour maps, wherein the contour maps are used to indicate the relationship between the development height of the water-conducting fracture zone of the 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 in which the fractures connect the pore water or groundwater in the rock formation so that water flows to 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 formation due to mining disturbance, the plurality of contour maps include one or more of the first contour map, the second contour map and the 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 mined working face adjacent to 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] A second acquisition unit 20 is used to acquire the development heights of all the water-conducting fracture zones at multiple time points, wherein the time points correspond to the development heights one by one;

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

[0174] In the specific 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 acquires the pumping test information of the mined part of the current working face, wherein the pumping test information includes one or more of the layer position of the first aquifer pumped by the borehole in the pumping test, the first permeability coefficient of the first aquifer when each borehole is pumped, and the first thickness of the first aquifer; the second acquisition module acquires the hardness information of the overlying rock formation of the mined part of the current working face; based on the hardness information, the first calculation module is used to calculate the first development height of the water-conducting fracture zone of the mined part of the current working face; the first generation module is used to generate the first contour map using ArcGIS software based on the pumping test information and the first development height of the water-conducting fracture zone.

[0175] In this solution, the first contour map is generated based on the actual data of the mined working face and GIS technology, which can provide a direct relationship between the development height of the water-conducting fracture zone and the hydrogeological parameters of the aquifer, so that the water inflow can be predicted based on the accurate first contour map in the future, which helps to improve the efficiency of mine water hazard prevention and control.

[0176] In the specific 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 used to obtain the drainage information of the above-mentioned mined working face adjacent to the above-mentioned current working face, wherein the above-mentioned drainage information includes one or more of the water outlet position, the end hole position, the water volume, the pressure, the drainage time and the water inflow change information of the drainage hole; according to Darcy's law, the second calculation module is used to calculate the second permeability coefficient corresponding to the above-mentioned drainage information; the first determination module is used to determine the second thickness of the second aquifer according to the above-mentioned drainage information and the above-mentioned second permeability coefficient; the fourth acquisition module is used to obtain the second development height, wherein the above-mentioned second development height is the above-mentioned development height corresponding to the hardness information of the overlying rock stratum of the above-mentioned mined working face adjacent to the above-mentioned current working face; the second generation module is used to generate the above-mentioned second contour map according to the above-mentioned second permeability coefficient, the above-mentioned second thickness and the above-mentioned second development height using ArcGIS software.

[0177] In this solution, based on the water release information of the adjacent mined working face, the second permeability coefficient and second thickness of the second aquifer are calculated, and combined with the second development height of the water-conducting fracture zone, a second contour map is generated. Using ArcGIS for data visualization and spatial analysis can make the complex changes in hydrogeological parameters intuitively presented, so that the subsequent water inflow can be predicted based on the accurate second contour map, which helps to improve the efficiency of mine water hazard prevention and control.

[0178] During the specific implementation process, the determination unit includes a fifth acquisition module, a second determination module and a third determination module. The fifth acquisition module is used to obtain the thematic map, wherein the thematic map is screened and determined based on the multiple contour maps; the second determination module is used to determine all the aquifers affected by the part to be mined of the current working face based on the multiple development heights, and obtain summary information; the third determination module is used to determine the water yield based on at least the summary information, wherein different aquifers have different permeability coefficients and thicknesses in the thematic map.

[0179] In this scheme, multiple contour maps and the development height of water-conducting fracture zones are used to determine the aquifers that may be affected by the part to be mined at the current working face, and ultimately the water inflow is calculated, which can further improve the accuracy of water inflow prediction.

[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 voice compares the first contour map and the third contour map to obtain a first similarity; the first fusion submodule is used to fuse the first contour map and the third contour map when the first similarity is greater than or equal to a first similarity threshold to obtain the thematic map, wherein the fusion method includes at least one or more of weighted average, Bayesian, and Kalman filtering; the second comparison submodule is used to compare the second contour map and the third contour map when the first similarity is less than the first similarity threshold to obtain a second similarity; the first determination submodule is used to determine that the third contour map is 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; the second fusion submodule is used to fuse the second contour map and the third contour map when the second similarity is less than the second similarity threshold to obtain the thematic map.

[0181] In this scheme, a more accurate thematic map can be obtained by merging and screening multiple contour maps, which can improve the accuracy of subsequent water inflow prediction.

[0182] In some embodiments, the second determination module includes an acquisition submodule and a second determination submodule, the acquisition submodule is used to obtain the stratigraphic relationship corresponding to the above-mentioned development height at each of the above-mentioned time points, wherein the above-mentioned stratigraphic relationship is the relationship between the above-mentioned aquifers affected when the above-mentioned water-conducting fracture zone reaches the above-mentioned development height; the second determination submodule is used to determine all of the above-mentioned aquifers affected according to all of the above-mentioned stratigraphic relationships corresponding to all of the above-mentioned time points, and obtain the above-mentioned summary information.

[0183] In this scheme, by comprehensively analyzing the laws of aquifer parameters and time changes, and considering the impact of the development of water-conducting fracture zones on different aquifers during the exploitation process, the stratigraphic relationship corresponding to the development height at each time point can be obtained, and all the aquifers affected can be determined based on this information, and summary information can be generated. In this way, all the aquifers affected by the exploitation process can be determined, and then the water inflow can be finely predicted later.

[0184] In some embodiments, the third determination module includes a calculation submodule and a statistical submodule. The calculation submodule is used to calculate multiple sub-water yields using the large well method according to the above-mentioned permeability coefficients of all the above-mentioned aquifers and the above-mentioned thicknesses of all the above-mentioned aquifers corresponding to the above-mentioned summary information, wherein the above-mentioned sub-water yields correspond one-to-one to the above-mentioned aquifers; the statistical submodule is used to calculate the sum of all the above-mentioned sub-water yields according to the above-mentioned summary information to obtain the above-mentioned water yield.

[0185] In this scheme, the sum of all sub-water inflows can be counted, which not only can obtain the sub-water inflow prediction of the working face at different mining stages, but also the overall water inflow. Since the temporal and spatial characteristics of the mining process are taken into consideration, the accuracy of the prediction results is further improved.

[0186] The above-mentioned device for determining the water inflow of a mine includes a processor and a memory, and the above-mentioned first acquisition unit, the second acquisition unit and the determination unit are all stored in the memory as program units, and the processor executes the above-mentioned program units stored in the memory to realize corresponding functions. The above-mentioned modules are all located in the same processor; or, the above-mentioned modules are located in different processors in any combination.

[0187] The processor includes a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the problem that the fine prediction of the mine water inflow cannot be achieved in the prior art can be solved by adjusting the kernel parameters.

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

[0189] An embodiment of the present invention provides a computer-readable storage medium, which includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the method for determining the water inflow of the mine.

[0190] An embodiment of the present invention provides a processor, and the processor is used to run a program, wherein the method for determining the water inflow of the mine is executed when the program is run.

[0191] An embodiment of the present invention provides a device, the device includes a processor, a memory, and a program stored in the memory and executable on the processor, and when the processor executes the program, at least the steps of the method for determining the water inflow of a mine are implemented. The device in this article can be a server, a PC, a PAD, a mobile phone, etc.

[0192] A computer program product includes a non-volatile computer-readable storage medium, wherein the non-volatile computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for determining the water inflow of the above-mentioned mine in each embodiment of the present application are implemented.

[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 are configured to be executed by the one or more processors, and the one or more programs include a method for determining the water inflow of any of the above-mentioned mines.

[0194] Obviously, those skilled in the art should understand that the above modules or steps of the present invention can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order than here, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. Thus, the present invention is not limited to any specific combination of hardware and software.

[0195] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0196] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0197] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0198] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

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

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

[0201] Computer readable media include 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 disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0202] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

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

[0204] 1) The method for determining the water inflow of a mine in the present 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, and 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, and then make a fine prediction of the water inflow of the to-be-mined part of the current working face based on the temporal and spatial characteristics. This method is more refined than the method of directly dividing the working face into a whole to predict the water inflow in the prior art, and realizes the fine prediction of the water inflow based on the temporal and spatial relationship between the development height of the water-conducting fracture zone and the aquifer.

[0205] 2) The device for determining the water inflow of a mine in the present 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, and then make a fine prediction of the water inflow of the to-be-mined part of the current working face based on the temporal and spatial characteristics. This is more refined than the method of directly dividing the working face into a whole to predict the water inflow in the prior art, and realizes the fine prediction of the water inflow based on the temporal and spatial relationship between the development height of the water-conducting fracture zone and the aquifer.

[0206] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for determining the water inflow of a mine, characterized in that: include: Acquire multiple contour maps, 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 of the affected aquifer and the thickness of the aquifer, the water-conducting fracture zone is a fracture zone in which the fractures connect the pore water or groundwater in the rock formation so that water flows to 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 formation 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 part of the current working face, the second contour map is the contour map of the mined working face adjacent to 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 heights of all the water-conducting fracture zones at multiple time points, wherein the time points and the development heights correspond one to one; Based on the multiple contour maps and the multiple development heights, the water inflow of the to-be-mined portion of the current working face is determined.

2. The method according to claim 1, characterized in that Get multiple contour plots, including: Acquire pumping test information of the mined portion of the current working face, wherein the pumping test information includes one or more of the layer position of the first aquifer pumped from the borehole in the pumping test, the first permeability coefficient of the first aquifer when each borehole is pumped, and the first thickness of the first aquifer; Obtaining hardness information of the overlying rock formation of the mined portion of the current working face; Calculating the first development height of the water-conducting fracture zone of the mined portion of the current working face according to the hardness information; The first contour map is generated using ArcGIS software based on the pumping test information and the first development height of the water-conducting fracture zone.

3. The method according to claim 2, characterized in that Get multiple contour plots, including: Acquire the drainage information of the mined working face adjacent to the current working face, wherein the drainage information includes one or more of the water outlet position, end hole position, water volume, pressure, drainage time and water inflow change information of the drainage hole; According to Darcy's law, calculating a second permeability coefficient corresponding to the water release information; determining a second thickness of a second aquifer according to the water release information and the second permeability coefficient; Acquire a second development height, wherein the second development height is the development height corresponding to the hardness information of the overlying rock formation of the mined working face adjacent to the current working face; The second contour map is generated using ArcGIS software according to the second permeability coefficient, the second thickness and the second development height.

4. The method according to claim 1, characterized in that: Determining the water inflow of the to-be-mined portion of the current working face according to the plurality of contour maps and the plurality of development heights includes: Acquire a thematic map, wherein the thematic map is determined by screening based on a plurality of the contour maps; According to the plurality of development heights, all the aquifers affected by the to-be-mined portion of the current working face are determined to obtain summary information; The water yield is determined at least based on the summary information, wherein different aquifers have different permeability coefficients and thicknesses in the thematic map.

5. The method according to claim 4, characterized in that Get thematic maps, including: Comparing the first contour map and the third contour map to obtain a first similarity; When 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 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 contour map with the third contour map to obtain a second similarity; In the case where the second similarity is greater than or equal to a second similarity threshold, determining the third contour map to be the thematic map, wherein the first similarity threshold is greater than the second similarity threshold; When the second similarity is less than the second similarity threshold, the second contour map and the third contour map are fused to obtain the thematic map.

6. The method according to claim 4, characterized in that According to the plurality of development heights, all the aquifers affected by the part to be mined of the current working face are determined to obtain summary information, including: Obtaining the stratigraphic relationship corresponding to the development height at each time point, wherein the stratigraphic relationship is the relationship of the aquifer affected when the water-conducting fracture zone reaches the development height; According to all the corresponding stratigraphic relationships at all the time points, all the affected aquifers are determined to obtain the summary information.

7. The method according to claim 4, characterized in that Determining the water inflow at least based on the summary information includes: According to the permeability coefficients of all the aquifers affected and the thicknesses of all the aquifers corresponding to the summary information, a plurality of sub-flow rates are calculated using the large well method, wherein the sub-flow rates correspond to the affected aquifers one by one; According to the summary information, the sum of all the sub-water inflows is calculated to obtain the water inflow.

8. A device for determining the water inflow of a mine, characterized in that: include: A first acquisition unit is used to acquire a plurality of contour maps, wherein the contour map is used to indicate the relationship between the development height of the water-conducting fracture zone of the 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 in which the fractures connect the pore water or groundwater in the rock formation so that water flows to 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 formation due to mining disturbance, the 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 mined working face adjacent to the current working face, and the third contour map is the contour map of the to-be-mined part of the current working face; A second acquisition unit is used to acquire the development heights of all the water-conducting fracture zones at multiple time points, wherein the time points correspond to the development heights one by one; A determination unit is used to determine the water inflow of the to-be-mined portion of the current working face based on the multiple contour maps and the multiple development heights.

9. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method for determining the water inflow of a mine as described in any one of claims 1 to 7 are implemented.

10. A water inflow detection system, characterized in that: include: One or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include a method for determining the water inflow of a mine as described in any one of claims 1 to 7.

Images