Water-resisting layer reconstruction range determination method and device, electronic equipment and storage medium

Abstract

The invention provides a water-resisting layer reconstruction range determination method and device, electronic equipment and a storage medium, and the method comprises the steps: obtaining geological parameters of a target mine, determining a water filling source based on the geological parameters, and constructing a spatial position partition map corresponding to a soil layer weak region range in combination with a triaxial seepage algorithm; according to the method, a water-rich partition map constructed based on a transient electromagnetic algorithm and a goaf ponding partition map are combined, a grouting treatment area selection map is constructed through continuous superposition, and a proper water-resisting layer reconstruction range can be rapidly and accurately selected based on an area selection index algorithm for grouting and water plugging so as to realize water-resisting layer reconstruction, so that treatment is carried out while mining is carried out in real time, and underground safety production is guaranteed.

Description

Technical Field

[0001] The present application relates to the technical field of water-preserving coal mining, and in particular to a method, device, electronic device and storage medium for determining a reconstruction range of an aquiclude. Background Art

[0002] The high-intensity mining of coal resources has caused the water-conducting fracture zone to develop to the top of the bedrock aquifer. In areas where the key impermeable layer of the overlying soil layer is missing (weak), the shallow Quaternary aquifer with strong water richness directly (overflow) recharges the bedrock aquifer through the water-conducting fracture zone, thus becoming a stable source of mine water replenishment, which can easily cause roof water disasters and endanger underground production safety.

[0003] Based on the above situation, grouting can be performed in the vertical water passage of the Quaternary aquifer to construct a horizontal water-blocking curtain to reconstruct the aquifer, cut off the supply of Quaternary groundwater to the mine water, and achieve water reduction underground and water conservation on the shallow surface. Therefore, how to quickly and accurately select the reconstruction range of the aquifer during grouting and water blocking has become a technical problem that needs to be solved urgently. Summary of the invention

[0004] In view of this, the purpose of this application is to propose a method for determining the reconstruction range of an impermeable layer to solve the above-mentioned technical problems.

[0005] Based on the above purpose, the first aspect of the present application provides a method for determining a reconstruction range of an aquiclude, comprising:

[0006] Obtaining geological parameters of the target mine, and determining a water source for water filling based on the geological parameters;

[0007] Determine the scope of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and construct a spatial position zoning map corresponding to the scope of the weak area of ​​the soil layer;

[0008] Determine the distribution range of the horizontal water-rich area based on the water source by transient electromagnetic algorithm, and construct a horizontal water-rich zoning map corresponding to the distribution range of the horizontal water-rich area; determine the distribution range of the vertical water-rich area based on the water source by transient electromagnetic algorithm, and construct a vertical water-rich zoning map corresponding to the distribution range of the vertical water-rich area;

[0009] Determine the water inflow from the underground drainage point according to the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow from the underground drainage point, and construct a goaf water accumulation zoning map corresponding to the distribution range of water accumulation in the goaf;

[0010] The spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map are superimposed to construct a grouting treatment area selection map;

[0011] The grouting treatment selection map is processed by a selection index algorithm, and the reconstruction range of the aquiclude is determined from the grouting treatment selection map.

[0012] Optionally, the determining the scope of the weak area of ​​the soil layer by a triaxial seepage algorithm based on the geological parameters includes:

[0013] Obtaining a stratigraphic sample of the target mine;

[0014] Recording permeability parameters of the formation sample under a preset pressure, wherein the permeability parameters include a permeability coefficient;

[0015] The area with a permeability greater than a preset permeability threshold in the stratum area corresponding to the stratum sample is taken as the weak area range of the soil layer.

[0016] Optionally, the determining the distribution range of the horizontal water-rich area based on the water source by using a transient electromagnetic algorithm includes:

[0017] Determine the target formation area corresponding to the water source;

[0018] Determine the induced eddy currents generated in the target formation area under the horizontal pulsed magnetic field, and determine the resistivity of the target formation area based on the induced eddy currents;

[0019] The area in the target formation area where the resistivity is less than a preset resistivity threshold is used as the distribution range of the horizontal water-rich area.

[0020] Optionally, the determining the distribution range of the vertical water-rich area based on the water source by using a transient electromagnetic algorithm includes:

[0021] Determine the target formation area corresponding to the water source;

[0022] determining an induced eddy current generated in the target formation area under a vertical pulsed magnetic field, and determining the resistivity of the target formation area based on the induced eddy current;

[0023] The area in the target formation area where the resistivity is less than a preset resistivity threshold is taken as the distribution range of the vertical water-rich area.

[0024] Optionally, determining the water inflow at the underground drainage point according to the geological parameters includes:

[0025] determining a hydraulic gradient based on the geological parameters;

[0026] Determining the permeability coefficient and water flow cross-sectional area of ​​the target mine;

[0027] The hydraulic gradient, the permeability coefficient and the water flow cross-sectional area are multiplied to obtain the water inflow at the underground drainage point.

[0028] Optionally, the determining the distribution range of water accumulation in the goaf area based on the water inflow from the underground drainage point includes:

[0029] Collecting apparent resistivity values ​​at different depths underground of the target mine by high-density electrical algorithms;

[0030] Performing imaging processing based on the apparent resistivity value to obtain a resistivity distribution map;

[0031] The area in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is regarded as a low-resistance area;

[0032] The water gushing out from the underground drainage point is used to determine the distribution range of water accumulation in the goaf from the low-resistance area.

[0033] Optionally, the step of processing the grouting treatment selection map by using a selection index algorithm and determining the reconstruction range of the aquiclude from the grouting treatment selection map comprises:

[0034] Any partition map of the grouting treatment area selection map is used as a target partition map, each target sub-region map includes multiple spatial positions, and any spatial position in the target partition map is used as a target spatial position;

[0035] Based on the target spatial position and the preset weight corresponding to the target spatial position, the governance constituency index is determined by the following formula:

[0036]

[0037] Where EI represents the governance constituency index, W i (x, y) represents the preset weight corresponding to the target spatial position x, y in the i-th target partition map, f i (x, y) represents the value of the target spatial position x, y in the i-th target partition map, n represents the total number of target partition maps, and i represents the order of the target partition maps;

[0038] The target partition map where the governance area index is greater than or equal to the preset index threshold is used as the reconstruction range of the impermeable layer.

[0039] Based on the same inventive concept, the second aspect of the present application provides a device for determining a reconstruction range of an aquiclude, comprising:

[0040] A water filling source determination module is configured to obtain geological parameters of the target mine and determine the water filling source based on the geological parameters;

[0041] A first partition map construction module is configured to determine the range of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and to construct a spatial position partition map corresponding to the range of the weak area of ​​the soil layer;

[0042] The second zoning diagram construction module is configured to determine the distribution range of the horizontal water-rich area based on the water source by transient electromagnetic algorithm, and construct a horizontal water-rich zoning diagram corresponding to the distribution range of the horizontal water-rich area, and determine the distribution range of the vertical water-rich area based on the water source by transient electromagnetic algorithm, and construct a vertical water-rich zoning diagram corresponding to the distribution range of the vertical water-rich area;

[0043] A third partition map construction module is configured to determine the water inflow of the underground drainage point according to the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow of the underground drainage point, and construct a goaf water accumulation partition map corresponding to the distribution range of water accumulation in the goaf;

[0044] A treatment area selection map construction module is configured to perform superposition processing on the spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map to construct a grouting treatment area selection map;

[0045] The reconstruction range determination module is configured to process the grouting treatment selection map through a selection index algorithm, and determine the reconstruction range of the impermeable layer from the grouting treatment selection map.

[0046] Based on the same inventive concept, the third aspect of the present application provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executed by the processor, wherein the processor implements the method described in the first aspect above when executing the computer program.

[0047] Based on the same inventive concept, the fourth aspect of the present application provides a non-transitory computer-readable storage medium, which stores computer instructions, and the computer instructions are used to enable a computer to execute the method described in the first aspect above.

[0048] From the above, it can be seen that the method, device, electronic device and storage medium for determining the range of reconstruction of the aquiclude provided in the present application, by acquiring the geological parameters of the target mine and determining the water source based on the geological parameters, and combining the three-axis seepage algorithm, construct a spatial position zoning map corresponding to the range of the weak area of ​​the soil layer, and combine it with the water-rich zoning map constructed based on the transient electromagnetic algorithm, and the water accumulation zoning map of the goaf, and continuously superimpose and construct a grouting treatment selection area map. Based on the selection area index algorithm, the appropriate aquiclude reconstruction range can be selected quickly and accurately for grouting and water blocking to achieve aquiclude reconstruction, so as to ensure safe production underground in real time during mining. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions in the present application or related technologies, the drawings required for use in the embodiments or related technical descriptions are briefly introduced below. Obviously, the drawings described below are only embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0050] Figure 1 A flowchart of a method for determining a reconstruction range of an aquiclude according to an embodiment of the present application;

[0051] Figure 2A A schematic diagram of a missing (weak) zoning map of the coal seam roof soil layer according to an embodiment of the present application;

[0052] Figure 2B A schematic diagram of a water-rich zoning diagram of a weathered bedrock aquifer horizontally disposed on a coal seam roof according to an embodiment of the present application;

[0053] Figure 2C A schematic diagram of a vertical groundwater richness zoning diagram of a coal seam roof according to an embodiment of the present application;

[0054] Figure 2D A schematic diagram of a water accumulation zoning diagram of an underground goaf area according to an embodiment of the present application;

[0055] Figure 2E A schematic diagram of a water-retaining coal mining roof grouting treatment area selection diagram according to an embodiment of the present application;

[0056] Figure 3 This is a structural block diagram of a device for determining a reconstruction range of an impermeable layer according to an embodiment of the present application;

[0057] Figure 4 A schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0058] In order to make the objectives, technical solutions and advantages of the present application more clearly understood, the present application is further described in detail below in combination with specific embodiments and with reference to the accompanying drawings.

[0059] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present application should be the usual meanings understood by people with ordinary skills in the field to which the present application belongs. The "first", "second" and similar words used in the embodiments of the present application do not represent any order, quantity or importance, but are only used to distinguish different components. "Including" or "comprising" and similar words mean that the elements or objects appearing in front of the word cover the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0060] It is understandable that before using the technical solutions of each embodiment of the present application, the type, scope of use, usage scenarios, etc. of the personal information involved will be informed to the user in an appropriate manner, and the user's authorization will be obtained.

[0061] For example, in response to receiving an active request from a user, a prompt message is sent to the user to clearly remind the user that the operation requested to be performed will require obtaining and using the user's personal information. Thus, the user can independently choose whether to provide personal information to the electronic device, application, server, storage medium or other software or hardware that performs the operation of the technical solution of the present application according to the prompt message.

[0062] As an optional but non-limiting implementation, in response to receiving the user's active request, the prompt information may be sent to the user in the form of a pop-up window, in which the prompt information may be presented in text form. In addition, the pop-up window may also carry a selection control for the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0063] It is understandable that the above notification and the process of obtaining user authorization are merely illustrative and do not constitute a limitation on the implementation method of this application. Other methods that meet relevant laws and regulations may also be applied to the implementation method of this application.

[0064] The high-intensity mining of coal resources has caused the water-conducting fracture zone to develop to the top of the bedrock aquifer. In areas where the key impermeable layer of the overlying soil layer is missing (weak), the shallow Quaternary aquifer with strong water richness directly (overflow) recharges the bedrock aquifer through the water-conducting fracture zone, thus becoming a stable source of mine water replenishment, which can easily cause roof water disasters and endanger underground production safety.

[0065] Based on the above situation, the aquiclude can be recreated by grouting in the vertical water channel of the Quaternary aquifer to construct a horizontal waterproof curtain, cut off the supply of Quaternary groundwater to the mine water, and achieve water reduction underground and water conservation on the shallow surface.

[0066] The present application found that if the reconstruction of the aquiclude is achieved through horizontal curtain grouting in the entire area where the key aquiclude is missing (weak) in the soil layer, it will undoubtedly incur huge economic costs. On the other hand, since the formation of the stratum is affected by a variety of complex internal and external dynamic geological effects such as sedimentation, structure, and weathering, its water storage space, water richness, etc. have significant variation characteristics in spatial distribution. Therefore, how to use this characteristic of the stratum to quickly and accurately select the key range for aquiclude reconstruction is a problem that needs to be solved urgently.

[0067] An embodiment of the present application provides a method for determining the reconstruction range of an aquiclude, by acquiring the geological parameters of a target mine, determining the water source based on the geological parameters, and combining a three-axis seepage algorithm to construct a spatial position zoning map corresponding to the range of weak areas of the soil layer, and combining the water-rich zoning map constructed based on a transient electromagnetic algorithm, as well as the water accumulation zoning map of the goaf, to continuously superimpose and construct a grouting treatment selection area map, and based on the selection area index algorithm, a suitable aquiclude reconstruction range can be quickly and accurately selected for grouting and water blocking to achieve aquiclude reconstruction, so as to ensure real-time mining and treatment and ensure safe production underground.

[0068] like Figure 1 As shown, the method of this embodiment includes:

[0069] Step 101, obtaining geological parameters of a target mine, and determining a water source for water filling based on the geological parameters.

[0070] In this step, geological parameters are various data and information used to describe geological structures and geological and hydrological characteristics, such as the type, thickness, distribution of rock formations, geological structures (such as faults and folds), depth of groundwater levels, permeability of rocks, and other geological factors that may affect mine stability.

[0071] Obtaining geological parameters usually requires geological exploration, which includes geological surveys, geophysical exploration (such as seismic exploration, electrical exploration), hydrogeological exploration, etc. These exploration activities can provide detailed information about the geological structure of the mining area.

[0072] For example, conduct stratigraphic surveys and production status surveys on mines, collect mine geological and hydrogeological data, determine the boundaries of the study area, stratigraphic structure, aquifer type, hydrogeological characteristics, and hydraulic connections between aquifers (i.e. geological parameters).

[0073] The source of water is the source of water that may enter the mine. These sources may include groundwater, surface water (such as rainwater, rivers, lakes, etc.) or water from other mines.

[0074] By analyzing geological parameters, it is possible to determine which rock layers or geological structures are potential water-filling channels. For example, highly permeable rock layers or faults may be the main path for groundwater flow. In addition, the height of the groundwater table and the direction of groundwater flow are also important factors in determining the source of water filling.

[0075] For example, the coal-bearing strata of this coal mine in the western region are the Yan'an Formation of the Middle Jurassic System. The overall characteristics of the occurrence of coal seams and aquifers are the symbiosis of coal and water, with the aquifer on top and the coal seam below.

[0076] The aquifers in the well field are divided into Quaternary phreatic aquifer, Jurassic Anding Formation weathered bedrock aquifer, and Jurassic bedrock aquifer from top to bottom.

[0077] The aquiclude is mainly a Neogene laterite aquiclude (there are missing areas);

[0078] At present, the main mining area is Jurassic Yan'an Formation 2-2. The water-conducting fracture zone generated by the coal seam has developed to the weathered bedrock aquifer, resulting in the weathered bedrock groundwater on the coal seam roof and Quaternary groundwater being the main sources of water in weak and missing soil areas.

[0079] Step 102: determine the scope of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and construct a spatial position zoning map corresponding to the scope of the weak area of ​​the soil layer.

[0080] In this step, the triaxial seepage algorithm takes into account the permeability characteristics of the soil layer in three-dimensional space and the flow law of water in the soil layer. By inputting geological parameters, the triaxial seepage algorithm can calculate key indicators such as the permeability coefficient and seepage velocity of the soil layer under different conditions, and then evaluate the stability and permeability of the soil layer.

[0081] The triaxial seepage algorithm can be used to simulate the stress-strain state and water flow path of the soil layer under different seepage conditions. By analyzing the simulation results, areas with high permeability, fast water flow velocity and poor stress state in the soil layer can be identified. These areas are often regarded as the weak area of ​​the soil layer (i.e. the distribution location of the missing (weak) area of ​​the key water-proof layer). The existence of weak areas may cause the soil layer to be destroyed under external loads or water flow, thereby causing geological disasters or engineering problems.

[0082] After determining the scope of the weak area of ​​the soil layer, the geographic information system (GIS) and computer-aided design (CAD) can be used to convert it into an intuitive spatial location zoning map. By marking these weak areas on the map and giving them different colors or symbols to distinguish their severity, a clear and intuitive spatial location zoning map can be formed, which helps to quickly understand the stability of the soil layer.

[0083] For example, firstly, the distribution of key aquicludes in the soil layer is explored through drilling, geophysical exploration and other methods; based on the soil samples obtained by drilling, indoor triaxial seepage experiments and physical and mechanical property tests are carried out to delineate the weak area of ​​the soil layer and draw a missing (weak) zoning map of the coal seam roof soil layer (i.e., a spatial location zoning map), such as Figure 2A As shown in the figure, it is found that there are a lot of missing red soil distribution (weak areas) in the exploration range of the coal mine. The missing (weak) zoning map of the coal seam roof soil layer includes missing (weak) soil layer areas and complete soil layer areas, among which: Figure 2A The horizontal axis represents the advancement distance in the strike direction of the working face, and the vertical axis represents the length in the dip direction of the working face; 5m represents the 5m contour line of soil thickness, and 10m represents the 10m contour line of soil thickness. In this example, the 10m contour line is used as the threshold of the weak soil area.

[0084] Step 103: Based on the water source, the distribution range of the horizontal water-rich area is determined by the transient electromagnetic algorithm, and a horizontal water-rich zoning map corresponding to the distribution range of the horizontal water-rich area is constructed. Based on the water source, the distribution range of the vertical water-rich area is determined by the transient electromagnetic algorithm, and a vertical water-rich zoning map corresponding to the distribution range of the vertical water-rich area is constructed.

[0085] In this step, the transient electromagnetic method (TEM) is suitable for detecting underground low-resistance bodies, such as aquifers, water-rich areas, etc. This method uses the principle of electromagnetic induction to infer the electrical parameters of the underground medium by emitting a horizontal transient electromagnetic field and observing its propagation and attenuation characteristics in the underground medium, and then analyzes the distribution and water richness of groundwater.

[0086] Based on the water source, the transient electromagnetic algorithm can accurately identify the distribution range of underground aquifers, the relative strength of water content, and the spatial morphology of abnormal areas. This is because the water-rich area has a significant difference in electrical properties relative to the surrounding strata, that is, the resistivity is lower. Through the measurement of transient electromagnetic method, the apparent resistivity data of the underground medium can be obtained, and then the apparent resistivity profile or slice diagram can be drawn, so as to intuitively show the distribution characteristics of the water-rich area.

[0087] Then, by analyzing the changing characteristics of apparent resistivity, the low-resistance abnormal area, that is, the water-rich area, is identified.

[0088] Combined with geological and hydrogeological data, a comprehensive analysis was conducted on the distribution range of water-rich areas, the relative strength of water content, and the spatial morphology of abnormal areas. Then, a horizontal water-rich zoning map was drawn based on the distribution characteristics of the water-rich areas.

[0089] In the zoning map, different colors or symbols can be used to represent areas with different water-rich levels, so as to intuitively show the distribution range of water-rich areas.

[0090] At the same time, key geological and hydrogeological information, such as faults, folds, aquifers, etc., can be marked on the zoning map to provide a scientific basis for the subsequent development and management of groundwater resources.

[0091] For example, the coal mine is mainly threatened by the fourth groundwater and weathered bedrock groundwater. Since the water-richness of the weathered bedrock aquifer and the underlying bedrock aquifer is obviously heterogeneous in spatial distribution, the reconstruction of the coal seam aquifer must fully consider the water-richness of the weathered bedrock aquifer and the main water-rich areas. According to the geological and hydrogeological conditions of the coal mine, the transient electromagnetic exploration method is used to explore the water-richness of the weathered bedrock aquifer and the underlying bedrock aquifer. Based on the low-resistance anomaly area of ​​the weathered bedrock, the water-richness zoning map of the weathered bedrock aquifer (i.e., the horizontal water-rich zoning map) is identified, such as Figure 2B As shown in the figure, it was found that there were three low-resistance abnormal areas (water-rich areas) in the weathered bedrock aquifer in the coal mine exploration range. The water-rich zoning diagram of the weathered bedrock aquifer horizontally in the coal seam roof includes strong water-rich areas and weak water-rich areas, among which: Figure 2B The horizontal axis represents the advancement distance of the working face in the strike direction, the vertical axis represents the length of the working face in the dip direction, and the text in the gray area represents the location of the highly water-rich area of ​​the horizontal weathered bedrock aquifer interpreted by geophysical exploration.

[0092] The vertical electromagnetic response signal of the underground medium is obtained, and the collected data is filtered, time-depth converted, and other processing is performed to eliminate noise and extract useful information. Then the processed data is imaged using an inversion algorithm to obtain a resistivity distribution map of the underground medium. After that, low-resistance abnormal areas are identified based on the resistivity distribution map. These areas usually correspond to water-rich areas. Therefore, by analyzing the resistivity data at different depths, the distribution range of the vertical water-rich area can be determined.

[0093] The identified water-rich areas are sorted according to their vertical distribution range to obtain parameters such as the depth and thickness of each water-rich area.

[0094] According to the parameters of the water-rich area, it is divided into different zones, usually according to factors such as the strength of water richness and distribution range.

[0095] Using geographic information system (GIS) or related mapping software, each water-rich zone is plotted on the map to form a vertical water-rich zone map.

[0096] For example, transient electromagnetic exploration is used to investigate the water-rich distribution of each aquifer on the vertical profile of the coal seam roof of the coal mine, identify the distribution and changes of groundwater drop funnels, draw a vertical groundwater-rich zoning map of the coal seam roof (i.e., vertical water-rich zoning map), and find the water passage for Quaternary groundwater to flow into the mine, such as Figure 2C As shown in the figure, it was found that there were two water channels within the exploration range of the coal mine. The vertical groundwater richness zoning map of the coal seam roof includes vertical water-conducting areas and vertical non-water-conducting areas, among which: Figure 2C The horizontal axis represents the advancement distance in the direction of the working face, and the vertical axis represents the elevation of the top plate of the working face at different vertical positions.

[0097] Step 104, determining the water inflow of underground drainage points according to the geological parameters, determining the distribution range of water accumulation in the goaf based on the water inflow of the underground drainage points, and constructing a goaf water accumulation zoning map corresponding to the distribution range of water accumulation in the goaf.

[0098] In this step, based on these geological parameters and using the hydrogeological model, the water inflow from each drainage point (usually a pre-set facility for extracting water from the mine) can be estimated. The water inflow refers to the amount of water flowing in or out through the drainage point within a certain period of time.

[0099] Goaf refers to areas in a mine that have been mined and abandoned. These areas are prone to water accumulation because they have lost their original rock support.

[0100] By analyzing the water inflow data of underground drainage points, combined with the mining history, geological structure and groundwater flow patterns of the mine, the possible distribution range of water accumulation in the goaf can be inferred. This usually involves predicting the water flow path, water accumulation depth and water accumulation volume.

[0101] After determining the distribution range of water accumulation in the goaf, the next step is to make a goaf water accumulation zoning map. This map should clearly show the water accumulation in different areas of the goaf, including the depth, range and possible flow direction of the water accumulation.

[0102] Among them, zoning maps usually use different colors or symbols to represent areas with different degrees of water accumulation, so as to intuitively understand the water accumulation conditions in the goaf and formulate effective drainage and safety management measures.

[0103] For example, through the analysis of water inflow at underground drainage points, hydrogeological and geophysical exploration, etc., the location, range and scale of water accumulation in the underground goaf of the coal mine can be determined, and a water accumulation zoning map of the underground goaf can be drawn, such as Figure 2D As shown, it was found that there were three waterlogged areas in the coal mine. The underground goaf waterlogging map includes waterlogged areas and non-waterlogged areas, among which: Figure 2DThe horizontal axis represents the advancement distance of the working face in the strike direction, the vertical axis represents the length of the working face in the dip direction, and the gray shaded part represents the scope of the water accumulation area in the goaf.

[0104] Step 105, superimposing the spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map to construct a grouting treatment area selection map.

[0105] In this step, based on the results of the overlay process, a grouting treatment area selection map is constructed to determine which areas are most suitable for grouting treatment. Grouting is used to reinforce soil, fill cavities, reduce groundwater infiltration, etc.

[0106] For example, based on the composite functions of GIS, the spatial position zoning map, the vertical water-rich zoning map described in the horizontal water-rich zoning map and the goaf water accumulation zoning map are superimposed to generate a water-preserving coal mining roof grouting treatment selection map (i.e., a grouting treatment selection map).

[0107] In the process of coal mining, the water accumulation in the mining area and the mining working face can be monitored in real time by continuously repeating steps 102, 103, 104 and 105, and the soil layer missing (weak) area, weathered bedrock water-rich area, groundwater funnel area, and underground goaf water distribution area can be obtained. Figure 2A , Figure 2B , Figure 2C and Figure 2D The water-retaining coal mining roof grouting treatment area selection map is generated by superposition, such as Figure 2E As shown, Figure 2E The horizontal axis represents the advancement distance of the working face in the strike direction, the vertical axis represents the length of the working face in the dip direction (Figures 1, 2, and 4 from top to bottom) or the elevation of the working face roof at different vertical positions (Figure 3 from top to bottom), and the circled area represents the distribution position of the strong water-rich area of ​​the vertical groundwater in the coal seam roof, that is, the distribution position of the vertical abnormal funnel (water channel).

[0108] Step 106, processing the grouting treatment selection map by using a selection index algorithm, and determining the reconstruction range of the aquiclude from the grouting treatment selection map.

[0109] In this step, the selection index algorithm is used to evaluate, select and determine the method for a specific area (such as a grouting treatment area). The algorithm is based on a variety of geological, engineering and environmental factors, and through calculation and analysis, it derives one or more indices or indicators used to guide decision-making. These indices or indicators can reflect information on the feasibility, effectiveness and cost of grouting treatment.

[0110] When processing the grouting treatment selection map, the selection index algorithm is used to evaluate and compare the various areas in the map. Based on the selection index algorithm, the grouting treatment index or index of each area is calculated.

[0111] Analyze and compare the calculation results to determine which areas are suitable for grouting treatment, as well as the priority and plan of grouting treatment.

[0112] Based on the geological conditions and grouting treatment requirements, determine the specific scope of the aquiclude reconstruction, including length, width and depth.

[0113] According to the scope of reconstruction of the aquiclude and the geological conditions, a detailed reconstruction plan is formulated, including the selection of grouting materials, determination of grouting pressure, construction of grouting holes, etc.

[0114] Carry out grouting construction according to the reconstruction plan to ensure the effective reconstruction of the waterproof layer and the effect of grouting treatment.

[0115] For example, the water-conducting channel in the weathered bedrock aquifer is selected as the grouting plugging area, and the specific plan for grouting treatment is determined, suitable grouting materials are selected and injected into the selected grouting plugging area, so as to rebuild the coal seam roof aquifer and realize treatment as mining progresses.

[0116] Specifically, in the soil layer missing (weak) zoning map of coal seam roof, the soil layer missing (weak) area and the soil layer intact area are assigned values ​​of 1 and 0, and the weights are assigned values ​​of 1 and 0 respectively; in the horizontal weathered bedrock aquifer water-rich zoning map of coal seam roof, the strong water-rich area and the weak water-rich area are assigned values ​​of 1 and 0, and the weights are assigned values ​​of 1 and 0 respectively; in the vertical groundwater richness zoning map of coal seam roof, the vertical water-conducting area and the vertical non-water-conducting area are assigned values ​​of 1 and 0, and the weights are assigned values ​​of 1 and 0 respectively; in the water accumulation zoning map of underground goaf, the water-accumulated area and the non-water-accumulated area are assigned values ​​of 1 and 0, and the weights are assigned values ​​of 1 and 0 respectively.

[0117] The composite overlay raster calculation function of GIS is used to process the generated water-preserving coal mining roof grouting treatment area map based on the selection index model to determine the reconstruction range of the impermeable layer.

[0118] Through the above scheme, by obtaining the geological parameters of the target mine and determining the water source based on the geological parameters, combined with the three-axis seepage algorithm, a spatial position zoning map corresponding to the weak area of ​​the soil layer is constructed, and combined with the water-rich zoning map constructed based on the transient electromagnetic algorithm and the water accumulation zoning map of the goaf, a grouting treatment selection map is continuously superimposed and constructed. Based on the selection index algorithm, the appropriate aquiclude reconstruction range can be selected quickly and accurately for grouting and water blocking to achieve aquiclude reconstruction, so as to ensure real-time mining and treatment to ensure safe production underground.

[0119] In some embodiments, in step 102, determining the scope of the weak area of ​​the soil layer by a triaxial seepage algorithm based on the geological parameters includes:

[0120] Step A1, obtaining a stratum sample of the target mine.

[0121] Step A2, recording the permeability parameters of the formation sample under a preset pressure, wherein the permeability parameters include a permeability coefficient.

[0122] Step A3: taking the area in the stratum area corresponding to the stratum sample where the permeability is greater than a preset permeability threshold as the weak area range of the soil layer.

[0123] In the above scenario, stratigraphic samples are collected from the target mine. These samples represent the physical and chemical characteristics of the strata within the mine and are the basis for subsequent analysis.

[0124] These formation samples are then subjected to laboratory testing under preset pressure conditions to record and measure their permeability parameters.

[0125] Permeability parameters mainly include the key indicator permeability coefficient. The permeability coefficient is a physical quantity that measures the difficulty of fluid passing through porous media (such as soil or rock). The larger the permeability coefficient, the stronger the permeability of the porous medium (such as soil or rock).

[0126] Finally, according to the permeability coefficient, the area in the formation area corresponding to the formation sample where the permeability coefficient is greater than the preset permeability coefficient threshold is identified as the weak area of ​​the soil layer.

[0127] The preset permeability coefficient threshold is a standard value set in advance based on factors such as mine safety, geological conditions, and mining needs. It is used to distinguish which areas of soil are more susceptible to penetration and erosion by fluids (such as water and gas), which may become potential safety hazards or mining problems.

[0128] By testing and analyzing the permeability parameters of the formation samples and combining them with the preset permeability coefficient threshold, weak areas of the soil in the target mine can be quickly identified and defined.

[0129] In some embodiments, in step 103, determining the distribution range of the horizontal water-rich area based on the water source by using a transient electromagnetic algorithm includes:

[0130] Step B1, determining the target formation area corresponding to the water source.

[0131] Step B2, determining the induced eddy currents generated in the target formation area under the horizontal pulsed magnetic field, and determining the resistivity of the target formation area based on the induced eddy currents.

[0132] Step B3: taking the area in the target formation area where the resistivity is less than a preset resistivity threshold as the distribution range of the horizontal water-rich area.

[0133] In the above scheme, geological layers or areas of formations where groundwater may be present are identified. This is usually based on geological surveys, geological maps, previous hydrogeological data, or preliminary results of geophysical surveys. The goal is to narrow the search area to focus on those formations that are most likely to contain groundwater.

[0134] By applying a horizontal pulsed magnetic field at or near the surface, the magnetic field will penetrate the ground and induce electric currents (eddy currents) in the conductive strata. The size and distribution of these eddy currents depend on the conductivity of the strata (i.e. the inverse of the resistivity).

[0135] The generation of eddy currents is a direct result of Faraday's law of electromagnetic induction. When the magnetic field changes over time, an electromotive force is generated in the conductor, which in turn drives the current to flow.

[0136] By measuring and analyzing the eddy currents induced by the pulsed magnetic field, the resistivity of the formation can be inferred. Resistivity is a measure of the degree to which a material hinders the flow of electric current. For groundwater-saturated formations, the resistivity is usually low due to the conductivity of water. Forward modeling and inversion analysis of electromagnetic fields are used to extract resistivity information from the observed data.

[0137] Finally, based on the resistivity measurements, areas with resistivity below a preset threshold are identified as areas of horizontal water abundance. This threshold is set based on knowledge of the local geological and hydrogeological conditions and analysis of resistivity data from similar areas.

[0138] By using the principle of electromagnetic induction and the difference in conductivity of the strata to indirectly detect the distribution of groundwater, non-invasive and efficient groundwater exploration can be achieved.

[0139] In some embodiments, in step 103, determining the distribution range of the vertical water-rich area based on the water source by using a transient electromagnetic algorithm includes:

[0140] Step C1, determining the target formation area corresponding to the water source.

[0141] Step C2, determining the induced eddy currents generated in the target formation area under the vertical pulse magnetic field, and determining the resistivity of the target formation area based on the induced eddy currents.

[0142] Step C3: taking the area of ​​the target formation area whose resistivity is less than a preset resistivity threshold as the distribution range of the vertical water-rich area.

[0143] In the above scheme, geological layers or areas of formations where groundwater may be present are identified. This is usually based on geological surveys, geological maps, previous hydrogeological data, or preliminary results of geophysical surveys. The goal is to narrow the search area to focus on those formations that are most likely to contain groundwater.

[0144] By applying a horizontal pulsed magnetic field at or near the surface, the magnetic field will penetrate the ground and induce electric currents (eddy currents) in the conductive strata. The size and distribution of these eddy currents depend on the conductivity of the strata (i.e. the inverse of the resistivity).

[0145] The generation of eddy currents is a direct result of Faraday's law of electromagnetic induction. When the magnetic field changes over time, an electromotive force is generated in the conductor, which in turn drives the current to flow.

[0146] By measuring and analyzing the eddy currents induced by the pulsed magnetic field, the resistivity of the formation can be inferred. Resistivity is a measure of the degree to which a material hinders the flow of electric current. For groundwater-saturated formations, the resistivity is usually low due to the conductivity of water. Forward modeling and inversion analysis of electromagnetic fields are used to extract resistivity information from the observed data.

[0147] Finally, based on the resistivity measurements, the vertical water-rich areas are identified as those with resistivity below a preset threshold, which is set based on knowledge of the local geological and hydrogeological conditions and analysis of resistivity data from similar areas.

[0148] By using the principle of electromagnetic induction and the difference in conductivity of the strata to indirectly detect the distribution of groundwater, non-invasive and efficient groundwater exploration can be achieved.

[0149] In some embodiments, in step 104, determining the water inflow of the underground drainage point according to the geological parameters includes:

[0150] Step D1, determining the hydraulic gradient according to the geological parameters.

[0151] Step D2, determining the permeability coefficient and water flow cross-sectional area of ​​the target mine.

[0152] Step D3, multiplying the hydraulic gradient, the permeability coefficient and the water flow cross-sectional area to obtain the water inflow at the underground drainage point.

[0153] In the above scheme, the hydraulic gradient describes the degree to which water flows in a certain direction under the action of gravity, usually expressed as the water level difference per unit length (such as how many centimeters or millimeters drop per meter).

[0154] Geological parameters may include ground tilt, groundwater level changes, soil type, etc., which will affect the calculation of hydraulic gradients.

[0155] By analyzing these geological parameters, the specific direction and speed of water flow in the mine can be determined, and the hydraulic gradient can be calculated.

[0156] Permeability (also known as hydraulic conductivity) is a physical quantity that describes the ability of soil or rock to allow water to pass through. It indicates the amount of water that passes through a unit area per unit time under a unit hydraulic gradient.

[0157] The cross-sectional area of ​​water flow refers to the cross-sectional area of ​​the specific channel (such as tunnels, cracks, pipes, etc.) through which water flows.

[0158] These two parameters are crucial in assessing the flow capacity of water in a mine.

[0159] The hydraulic gradient, the permeability coefficient and the water flow cross-sectional area are multiplied to obtain the water inflow at the underground drainage point:

[0160] This procedure is based on Darcy's Law, which describes the relationship between the rate of fluid flow in porous media and the hydraulic gradient, permeability and flow channel area.

[0161] Specifically, the water inflow (Q) can be calculated by the following formula: Q = K*I*A, where K is the permeability coefficient, I is the hydraulic gradient, and A is the cross-sectional area of ​​the water flow.

[0162] Through this formula, geological and hydrological parameters can be converted into specific water yield values, and the water yield of underground drainage points can be determined quickly and accurately.

[0163] In some embodiments, in step 104, determining the distribution range of water accumulation in the goaf based on the water inflow of the underground drainage point includes:

[0164] Step E1, collecting apparent resistivity values ​​at different underground depths of the target mine using a high-density electrical algorithm.

[0165] Step E2: performing imaging processing based on the apparent resistivity value to obtain a resistivity distribution map.

[0166] Step E3: The area in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is regarded as a low-resistance area.

[0167] Step E4, using the water inflow from the underground drainage point to determine the distribution range of water accumulation in the goaf from the low-resistance area.

[0168] In the above scheme, high-density electrical method is a geophysical exploration method that arranges dense electrode arrays on the surface or underground and measures the resistivity changes between different electrode combinations to infer the electrical structure of the underground medium.

[0169] In this process, data will be collected at different depths underground in the target mine to obtain the apparent resistivity value at each depth. Apparent resistivity is a physical quantity that reflects the electrical conductivity of underground media, and its value is affected by many factors such as media composition, water content, and temperature.

[0170] Based on the acquired apparent resistivity values, they are processed through a specific imaging algorithm (such as an inversion algorithm) to generate a resistivity distribution map.

[0171] The resistivity distribution map can intuitively display the resistivity distribution of different areas underground and is the basis for subsequent analysis.

[0172] In the resistivity distribution map, a preset apparent resistivity threshold is set. This threshold is usually determined based on experience or previous geological surveys to distinguish areas with different resistivity characteristics.

[0173] The area with apparent resistivity value less than this threshold value in the figure is regarded as low-resistance area. In mine detection, low-resistance area is often related to geological phenomena such as aquifer and water accumulation in goaf.

[0174] After the low-resistance areas are determined, these low-resistance areas are further analyzed in combination with the known water inflow information of the underground drainage points.

[0175] By analyzing the spatial relationship between water inflow and low-resistance areas, the specific distribution range of water accumulation in the goaf can be inferred. Areas with large water inflow often indicate areas where water accumulation is more concentrated or flows more actively.

[0176] In summary, this process uses high-density electrical algorithms to collect data, perform imaging processing, set thresholds to determine low-resistance areas, and finally combines water inflow information to quickly and accurately determine the distribution range of water accumulation in the goaf, providing important technical support for mine safety production and water accumulation control.

[0177] In some embodiments, step 106 includes:

[0178] Step F1, taking any partition map of the grouting treatment area selection map as a target partition map, each target sub-region map includes multiple spatial positions, and taking any spatial position in the target partition map as a target spatial position.

[0179] Step F2, based on the target spatial position and the preset weight corresponding to the target spatial position, determine the governance constituency index by the following formula:

[0180]

[0181] Where EI represents the governance constituency index, W i (x, y) represents the preset weight corresponding to the target spatial position x, y in the i-th target partition map, fi (x, y) represents the value of the target spatial position x, y in the i-th target partition map, n represents the total number of target partition maps, and i represents the order of the target partition maps.

[0182] Step F3, taking the target partition map where the governance area index is greater than or equal to the preset index threshold as the aquiclude reconstruction range.

[0183] In the above scheme, the grouting treatment area map is divided into multiple zone maps, each zone map represents a potential treatment or grouting area.

[0184] These partition maps are further subdivided into multiple spatial locations, each of which represents a specific point or area within the partition.

[0185] In this process, any partition map is selected as the target partition map for analysis. In the target partition map, any spatial position is selected as the target spatial position for further calculation.

[0186] Based on the target spatial location and the preset weight corresponding to the location, the governance constituency index is calculated using the following formula:

[0187]

[0188] Where EI represents the governance constituency index, W i (x, y) represents the preset weight corresponding to the target spatial position x, y in the i-th target partition map, f i (x, y) represents the value of the target spatial position x, y in the i-th target partition map of the evaluation unit, n represents the total number of target partition maps, and i represents the order of the target partition maps.

[0189] The preset weights corresponding to the target spatial positions may be based on various factors such as geological conditions, grouting requirements, and cost-effectiveness.

[0190] The governance constituency index is calculated comprehensively by considering the weight and location information of the target spatial location in the target zoning map.

[0191] The calculated governance constituency index is compared with the preset index threshold.

[0192] If the governance constituency index of a target zoning map is greater than or equal to the preset index threshold, the zoning map is considered to be suitable for aquiclude reconstruction.

[0193] For example, the area with a governance selection index of 4 is the water-retaining coal mining roof grouting governance area (i.e., the scope of reconstruction of the impermeable layer), and the area with a governance selection index less than 4 does not need to be governed. In addition, the scope of the governance area can be appropriately expanded based on the actual project to compensate for the calculation errors caused by the exploration accuracy of various exploration methods.

[0194] By comprehensively considering the weights and location information of the spatial locations in different zoning maps, it is possible to quickly and accurately determine which areas are most in need of grouting treatment to achieve aquiclude reconstruction.

[0195] It should be noted that the method of the embodiment of the present application can be performed by a single device, such as a computer or server. The method of this embodiment can also be applied to a distributed scenario and completed by multiple devices cooperating with each other. In the case of such a distributed scenario, one of the multiple devices can only perform one or more steps in the method of the embodiment of the present application, and the multiple devices will interact with each other to complete the described method.

[0196] It should be noted that the above describes some embodiments of the present application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the above embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0197] Based on the same inventive concept, corresponding to any of the above-mentioned embodiment methods, the present application also provides a device for determining the reconstruction range of an impermeable layer.

[0198] refer to Figure 3 , the device for determining the reconstruction range of the water-blocking layer comprises:

[0199] The water filling source determination module 301 is configured to obtain geological parameters of the target mine and determine the water filling source based on the geological parameters;

[0200] A first partition map construction module 302 is configured to determine the range of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and to construct a spatial position partition map corresponding to the range of the weak area of ​​the soil layer;

[0201] The second zoning diagram construction module 303 is configured to determine the distribution range of the horizontal water-rich area based on the water source by using the transient electromagnetic algorithm, and construct a horizontal water-rich zoning diagram corresponding to the distribution range of the horizontal water-rich area, and determine the distribution range of the vertical water-rich area based on the water source by using the transient electromagnetic algorithm, and construct a vertical water-rich zoning diagram corresponding to the distribution range of the vertical water-rich area;

[0202] The third partition map construction module 304 is configured to determine the water inflow of the underground drainage point according to the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow of the underground drainage point, and construct a goaf water accumulation partition map corresponding to the distribution range of water accumulation in the goaf;

[0203] The treatment selection area map construction module 305 is configured to perform superposition processing on the spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map to construct a grouting treatment selection area map;

[0204] The reconstruction range determination module 306 is configured to process the grouting treatment selection map through a selection index algorithm, and determine the reconstruction range of the impermeable layer from the grouting treatment selection map.

[0205] In some embodiments, the first partition map construction module 302 is specifically configured to:

[0206] Obtaining a stratigraphic sample of the target mine;

[0207] Recording permeability parameters of the formation sample under a preset pressure, wherein the permeability parameters include a permeability coefficient;

[0208] The area with a permeability greater than a preset permeability threshold in the stratum area corresponding to the stratum sample is taken as the weak area range of the soil layer.

[0209] In some embodiments, the second partition map construction module 303 is specifically configured to:

[0210] Determine the target formation area corresponding to the water source;

[0211] Determine the induced eddy currents generated in the target formation area under the horizontal pulsed magnetic field, and determine the resistivity of the target formation area based on the induced eddy currents;

[0212] The area in the target formation area where the resistivity is less than a preset resistivity threshold is used as the distribution range of the horizontal water-rich area.

[0213] In some embodiments, the second partition map construction module 303 is specifically configured to:

[0214] Determine the target formation area corresponding to the water source;

[0215] determining an induced eddy current generated in the target formation area under a vertical pulsed magnetic field, and determining the resistivity of the target formation area based on the induced eddy current;

[0216] The area in the target formation area where the resistivity is less than a preset resistivity threshold is taken as the distribution range of the vertical water-rich area.

[0217] In some embodiments, the third partition map construction module 304 is specifically configured to:

[0218] determining a hydraulic gradient based on the geological parameters;

[0219] Determining the permeability coefficient and water flow cross-sectional area of ​​the target mine;

[0220] The hydraulic gradient, the permeability coefficient and the water flow cross-sectional area are multiplied to obtain the water inflow at the underground drainage point.

[0221] In some embodiments, the third partition map construction module 304 is specifically configured to:

[0222] Collecting apparent resistivity values ​​at different depths underground of the target mine by high-density electrical algorithms;

[0223] Performing imaging processing based on the apparent resistivity value to obtain a resistivity distribution map;

[0224] The area in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is regarded as a low-resistance area;

[0225] The water gushing out from the underground drainage point is used to determine the distribution range of water accumulation in the goaf from the low-resistance area.

[0226] In some embodiments, the reconstruction range determination module 306 is specifically configured to:

[0227] Any partition map of the grouting treatment area selection map is used as a target partition map, each target sub-region map includes multiple spatial positions, and any spatial position in the target partition map is used as a target spatial position;

[0228] Based on the target spatial position and the preset weight corresponding to the target spatial position, the governance constituency index is determined by the following formula:

[0229]

[0230] Where EI represents the governance constituency index, W i (x, y) represents the preset weight corresponding to the target spatial position x, y in the i-th target partition map, f i (x, y) represents the value of the target spatial position x, y in the i-th target partition map, n represents the total number of target partition maps, and i represents the order of the target partition maps;

[0231] The target partition map where the governance area index is greater than or equal to the preset index threshold is used as the reconstruction range of the impermeable layer.

[0232] For the convenience of description, the above device is described in terms of functions divided into various modules. Of course, when implementing the present application, the functions of each module can be implemented in the same or multiple software and / or hardware.

[0233] The device of the above embodiment is used to implement the corresponding method for determining the reconstruction range of the impermeable layer in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be repeated here.

[0234] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments, the present application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for determining the reconstruction range of the impermeable layer described in any of the above embodiments is implemented.

[0235] Figure 4 A more specific schematic diagram of the hardware structure of an electronic device provided in this embodiment is shown, and the device may include: a processor 401, a memory 402, an input / output interface 403, a communication interface 404, and a bus 405. The processor 401, the memory 402, the input / output interface 403, and the communication interface 404 are connected to each other in communication within the device through the bus 405.

[0236] The processor 401 can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0237] The memory 402 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 402 can store an operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program codes are stored in the memory 402 and called and executed by the processor 401.

[0238] The input / output interface 403 is used to connect the input / output module to realize information input and output. The input / output module can be configured in the device as a component (not shown in the figure), or it can be externally connected to the device to provide corresponding functions. The input device may include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device may include a display, a speaker, a vibrator, an indicator light, etc.

[0239] The communication interface 404 is used to connect a communication module (not shown) to realize communication interaction between the device and other devices. The communication module can realize communication through a wired mode (such as USB, network cable, etc.) or a wireless mode (such as mobile network, WIFI, Bluetooth, etc.).

[0240] The bus 405 comprises a pathway for transmitting information between the various components of the device (eg, the processor 401 , the memory 402 , the input / output interface 403 , and the communication interface 404 ).

[0241] It should be noted that, although the above device only shows the processor 401, the memory 402, the input / output interface 403, the communication interface 404 and the bus 405, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, it can be understood by those skilled in the art that the above device may also only include the components necessary for implementing the embodiments of the present specification, and does not necessarily include all the components shown in the figure.

[0242] The electronic device of the above embodiment is used to implement the corresponding method for determining the reconstruction range of the water-proof layer in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be described in detail here.

[0243] Based on the same inventive concept, corresponding to any of the above-mentioned embodiment methods, the present application also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method for determining the reconstruction range of the impermeable layer as described in any of the above embodiments.

[0244] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. Information can be computer-readable instructions, data structures, modules of programs, 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, read-only compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape 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.

[0245] The computer instructions stored in the storage medium of the above embodiment are used to enable the computer to execute the method for determining the reconstruction range of the impermeable layer as described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0246] A person skilled in the art should understand that the discussion of any of the above embodiments is merely illustrative and is not intended to imply that the scope of the present application is limited to these examples. In line with the concept of the present application, the technical features in the above embodiments or different embodiments may be combined, the steps may be implemented in any order, and there are many other variations of the different aspects of the embodiments of the present application as described above, which are not provided in detail for the sake of simplicity.

[0247] In addition, to simplify the description and discussion, and in order not to make the embodiments of the present application difficult to understand, the known power supply / ground connection with the integrated circuit (IC) chip and other components may or may not be shown in the provided drawings. In addition, the device can be shown in the form of a block diagram to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that the details of the implementation of these block diagram devices are highly dependent on the platform to be implemented in the embodiments of the present application (that is, these details should be fully within the scope of understanding of those skilled in the art). In the case of elaborating specific details (e.g., circuits) to describe exemplary embodiments of the present application, it is obvious to those skilled in the art that the embodiments of the present application can be implemented without these specific details or when these specific details are changed. Therefore, these descriptions should be considered to be illustrative rather than restrictive.

[0248] Although the present application has been described in conjunction with specific embodiments of the present application, many replacements, modifications and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may use the embodiments discussed.

[0249] The embodiments of the present application are intended to cover all such substitutions, modifications and variations that fall within the broad scope of the present application. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present application should be included in the protection scope of the present application.

Claims

1. A method for determining a reconstruction range of an aquiclude, characterized in that: include: Obtaining geological parameters of the target mine, and determining a water source for water filling based on the geological parameters; Determine the scope of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and construct a spatial position zoning map corresponding to the scope of the weak area of ​​the soil layer; Determine the distribution range of the horizontal water-rich area based on the water source by transient electromagnetic algorithm, and construct a horizontal water-rich zoning map corresponding to the distribution range of the horizontal water-rich area; determine the distribution range of the vertical water-rich area based on the water source by transient electromagnetic algorithm, and construct a vertical water-rich zoning map corresponding to the distribution range of the vertical water-rich area; Determine the water inflow from the underground drainage point according to the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow from the underground drainage point, and construct a goaf water accumulation zoning map corresponding to the distribution range of water accumulation in the goaf; The spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map are superimposed to construct a grouting treatment area selection map; The grouting treatment selection map is processed by a selection index algorithm, and the reconstruction range of the aquiclude is determined from the grouting treatment selection map.

2. The method according to claim 1, characterized in that The determining of the weak area of ​​the soil layer by a triaxial seepage algorithm based on the geological parameters includes: Obtaining a stratigraphic sample of the target mine; Recording permeability parameters of the formation sample under a preset pressure, wherein the permeability parameters include a permeability coefficient; The area with a permeability greater than a preset permeability threshold in the stratum area corresponding to the stratum sample is taken as the weak area range of the soil layer.

3. The method according to claim 1, characterized in that The method of determining the distribution range of the water-rich area horizontally based on the water source by using a transient electromagnetic algorithm includes: Determine the target formation area corresponding to the water source; Determine the induced eddy currents generated in the target formation area under the horizontal pulsed magnetic field, and determine the resistivity of the target formation area based on the induced eddy currents; The area in the target formation area where the resistivity is less than a preset resistivity threshold is used as the distribution range of the horizontal water-rich area.

4. The method according to claim 1, characterized in that: The method of determining the distribution range of the vertical water-rich area based on the water source by using a transient electromagnetic algorithm includes: Determine the target formation area corresponding to the water source; determining an induced eddy current generated in the target formation area under a vertical pulsed magnetic field, and determining the resistivity of the target formation area based on the induced eddy current; The area in the target formation area where the resistivity is less than a preset resistivity threshold is taken as the distribution range of the vertical water-rich area.

5. The method according to claim 1, characterized in that Determining the water inflow at the underground drainage point according to the geological parameters includes: determining a hydraulic gradient based on the geological parameters; Determining the permeability coefficient and water flow cross-sectional area of ​​the target mine; The hydraulic gradient, the permeability coefficient and the water flow cross-sectional area are multiplied to obtain the water inflow at the underground drainage point.

6. The method according to claim 1, characterized in that The determining of the distribution range of water accumulation in the goaf area based on the water inflow from the underground drainage point includes: Collecting apparent resistivity values ​​at different depths underground of the target mine by high-density electrical method; Performing imaging processing based on the apparent resistivity value to obtain a resistivity distribution map; The area in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is regarded as a low-resistance area; The water gushing out from the underground drainage point is used to determine the distribution range of water accumulation in the goaf from the low-resistance area.

7. The method according to claim 1, characterized in that The step of processing the grouting treatment selection map by using a selection index algorithm and determining the reconstruction range of the aquiclude from the grouting treatment selection map includes: Any partition map of the grouting treatment area selection map is used as a target partition map, each target sub-region map includes multiple spatial positions, and any spatial position in the target partition map is used as a target spatial position; Based on the target spatial position and the preset weight corresponding to the target spatial position, the governance constituency index is determined by the following formula: Where EI represents the governance constituency index, W i (x, y) represents the preset weight corresponding to the target spatial position x, y in the i-th target partition map, f i (x, y) represents the value of the target spatial position x, y in the i-th target partition map, n represents the total number of target partition maps, and i represents the order of the target partition maps; The target partition map where the governance area index is greater than or equal to the preset index threshold is used as the reconstruction range of the impermeable layer.

8. A device for determining a water-blocking layer reconstruction range, characterized in that: include: A water filling source determination module is configured to obtain geological parameters of the target mine and determine the water filling source based on the geological parameters; A first partition map construction module is configured to determine the range of the weak area of ​​the soil layer through a triaxial seepage algorithm based on the geological parameters, and to construct a spatial position partition map corresponding to the range of the weak area of ​​the soil layer; The second zoning diagram construction module is configured to determine the distribution range of the horizontal water-rich area based on the water source by transient electromagnetic algorithm, and construct a horizontal water-rich zoning diagram corresponding to the distribution range of the horizontal water-rich area, and determine the distribution range of the vertical water-rich area based on the water source by transient electromagnetic algorithm, and construct a vertical water-rich zoning diagram corresponding to the distribution range of the vertical water-rich area; A third partition map construction module is configured to determine the water inflow of the underground drainage point according to the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow of the underground drainage point, and construct a goaf water accumulation partition map corresponding to the distribution range of water accumulation in the goaf; A treatment area selection map construction module is configured to perform superposition processing on the spatial position zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map and the goaf water accumulation zoning map to construct a grouting treatment area selection map; The reconstruction range determination module is configured to process the grouting treatment selection map through a selection index algorithm, and determine the reconstruction range of the impermeable layer from the grouting treatment selection map.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the program, the method according to any one of claims 1 to 7 is implemented.

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

Images