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

Abstract

The application provides a water-resisting layer reconstruction range determination method and device, electronic equipment and a storage medium. The geological parameters of a target mine are obtained, and a water filling source is determined based on the geological parameters. A spatial position zoning map corresponding to the weak area range of the soil layer is constructed by combining a triaxial seepage algorithm. A water-rich zoning map and a goaf water accumulation zoning map are constructed based on a transient electromagnetic algorithm, and a grouting treatment selection map is continuously constructed by superimposition. A suitable water-resisting layer reconstruction range can be quickly and accurately selected by using a selection index algorithm to realize grouting water plugging and water-resisting layer reconstruction, so as to realize real-time mining and treatment and ensure safe production in the mine.

Description

Technical Field

[0001] This application relates to the field of water-conserving coal mining technology, and in particular to a method, apparatus, electronic device and storage medium for determining the reconstruction range of a water-retaining layer. Background Technology

[0002] The high-intensity mining of coal resources has led to the development of water-conducting fracture zones up to the top of the bedrock aquifer. In areas where the key aquitard layer of the overlying soil is missing (or weak), the shallow Quaternary aquifer with strong water-bearing capacity directly (overflows) supplies the bedrock aquifer through the water-conducting fracture zone, thus becoming a stable source of mine water supply. This can easily cause roof water disasters and endanger underground safe production.

[0003] Based on the above, a horizontal aquitard can be constructed by grouting into the vertical water passage of the Quaternary aquifer to recreate the aquitard, cutting off the supply of mine water from the Quaternary groundwater and achieving water reduction underground and water retention at the surface. Therefore, how to quickly and accurately select the reconstruction range of the aquitard during grouting and water plugging has become an urgent technical problem to be solved. 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 the waterproof layer in order to solve the above-mentioned technical problems.

[0005] To achieve the above objectives, the first aspect of this application provides a method for determining the reconstruction range of an impermeable layer, comprising:

[0006] Obtain the geological parameters of the target mine, and determine the water source based on the geological parameters;

[0007] Based on the geological parameters, the range of the weak soil layer is determined by the triaxial seepage algorithm, and a spatial location zoning map corresponding to the range of the weak soil layer is constructed.

[0008] 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.

[0009] The water inflow at the underground drainage point is determined based on the geological parameters. The distribution range of water accumulation in the goaf is determined based on the water inflow at the underground drainage point. A goaf water accumulation zoning map corresponding to the distribution range of water accumulation in the goaf is constructed.

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

[0011] The grouting treatment selection map is processed using the selection area index algorithm to determine the reconstruction range of the waterproof layer from the grouting treatment selection map.

[0012] Optionally, determining the extent of the weak soil zone based on the geological parameters using a triaxial seepage algorithm includes:

[0013] Obtain formation samples from the target mine;

[0014] Record the permeability parameters of the formation sample under a preset pressure, the permeability parameters including the permeability coefficient;

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

[0016] Optionally, determining the distribution range of the horizontally rich water zone based on the water source using a transient electromagnetic algorithm includes:

[0017] Determine the target geological region corresponding to the water source;

[0018] The induced eddy current generated in the target stratum region under a horizontal pulsed magnetic field is determined, and the resistivity of the target stratum region is determined based on the induced eddy current;

[0019] The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the horizontally water-rich zone.

[0020] Optionally, determining the distribution range of the vertically rich water zone based on the water source using a transient electromagnetic algorithm includes:

[0021] Determine the target geological region corresponding to the water source;

[0022] The induced eddy current generated in the target stratum region under a vertical pulsed magnetic field is determined, and the resistivity of the target stratum region is determined based on the induced eddy current;

[0023] The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the vertical water-rich zone.

[0024] Optionally, determining the water inflow at the downhole drainage point based on the geological parameters includes:

[0025] The hydraulic gradient is determined based on the geological parameters.

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

[0027] The water inflow at the downhole drainage point is obtained by multiplying the hydraulic gradient, the permeability coefficient, and the cross-sectional area of ​​the water flow.

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

[0029] Apparent resistivity values ​​at different depths underground in the target mine were collected using a high-density electrical resistivity algorithm.

[0030] Imaging processing is performed based on the apparent resistivity value to obtain a resistivity distribution map;

[0031] The region in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is defined as a low-resistivity region.

[0032] The distribution range of water accumulation in the goaf is determined from the low-resistivity area by using the water inflow rate at the underground drainage point.

[0033] Optionally, the step of processing the grouting treatment selection map using the selection area index algorithm to determine the reconstruction range of the impermeable layer from the grouting treatment selection map includes:

[0034] Any partition map of the grouting treatment selection map is taken as the target partition map. Each target partition map includes multiple spatial locations. Any spatial location in the target partition map is taken as the target spatial location.

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

[0036]

[0037] Where EI represents the Governance Constituency Index, W i (x, y) represents the preset weights corresponding to the spatial locations x and y in the i-th target partition map, f i (x, y) represents the values ​​at target spatial positions x and 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 zoning map whose governance selection area index is greater than or equal to the preset index threshold is used as the reconstruction range of the waterproof layer.

[0039] Based on the same inventive concept, a second aspect of this application provides a device for determining the reconstruction range of a waterproof layer, comprising:

[0040] The water source determination module is configured to acquire the geological parameters of the target mine and determine the water source based on the geological parameters.

[0041] The first zoning map construction module is configured to determine the range of weak soil layers based on the geological parameters using a triaxial seepage algorithm, and to construct a spatial location zoning map corresponding to the range of weak soil layers.

[0042] The second partition map construction module is configured to determine the distribution range of the horizontal water-rich area based on the water-filling source using a transient electromagnetic algorithm, construct a horizontal water-rich partition map 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-filling source using a transient electromagnetic algorithm, construct a vertical water-rich partition map corresponding to the distribution range of the vertical water-rich area.

[0043] The third zoning map construction module is configured to determine the water inflow of the underground drainage point based on 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 zoning map corresponding to the distribution range of water accumulation in the goaf.

[0044] The treatment zone map construction module is configured to overlay the spatial location zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map, and the goaf water accumulation zoning map to construct the grouting treatment zone map;

[0045] The reconstruction range determination module is configured to process the grouting treatment selection map using a selection area index algorithm to determine the reconstruction range of the waterproof layer from the grouting treatment selection map.

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

[0047] Based on the same inventive concept, a fourth aspect of this application provides a non-transitory computer-readable storage medium that stores computer instructions for causing a computer to perform the method described in the first aspect above.

[0048] As can be seen from the above, the method, device, electronic equipment and storage medium for determining the reconstruction range of the aquitard provided in this application obtain the geological parameters of the target mine and determine the water source based on the geological parameters. Combined with the triaxial seepage algorithm, a spatial location zoning map corresponding to the range of the weak soil layer is constructed. Combined with the water-rich zoning map constructed based on the transient electromagnetic algorithm and the goaf water accumulation zoning map, a grouting treatment selection map is continuously superimposed. Based on the selection index algorithm, a suitable aquitard reconstruction range can be quickly and accurately selected for grouting and water plugging to achieve aquitard reconstruction, so as to treat the aquitard in real time and ensure safe production underground. Attached Figure Description

[0049] To more clearly illustrate the technical solutions in this application or related technologies, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the drawings described below are only embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0050] Figure 1 This is a flowchart of the method for determining the reconstruction range of the waterproof layer according to an embodiment of this application;

[0051] Figure 2A This is a schematic diagram of a zoning map showing the missing (weak) soil layers in the roof of a coal seam, according to an embodiment of this application.

[0052] Figure 2B This is a schematic diagram of the water-bearing zoning of the coal seam roof horizontally towards the weathered bedrock aquifer, according to an embodiment of this application.

[0053] Figure 2C This is a schematic diagram of the vertical groundwater zoning of the coal seam roof according to an embodiment of this application;

[0054] Figure 2D This is a schematic diagram of the water accumulation zone in the underground goaf according to an embodiment of this application;

[0055] Figure 2E This is a schematic diagram of the selected area for grouting treatment of the roof of water-retaining coal mining, according to an embodiment of this application.

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

[0057] Figure 4 This is a schematic diagram of an electronic device according to an embodiment of this application. Detailed Implementation

[0058] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with specific embodiments and the accompanying drawings.

[0059] It should be noted that, unless otherwise defined, the technical or scientific terms used in the embodiments of this application should have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms "first," "second," and similar terms used in the embodiments of this application do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" mean that the element or object preceding the word encompasses the elements or objects listed after the word and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" 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 understood that before using the technical solutions of the various embodiments in this application, users will be informed of the type, scope of use, and usage scenarios of the personal information involved in an appropriate manner, and user authorization will be obtained.

[0061] For example, upon receiving a user's active request, a prompt message is sent to the user to explicitly inform them that the requested operation will require the acquisition and use of the user's personal information. This allows the user to independently choose, based on the prompt message, whether to provide personal information to the software or hardware such as electronic devices, applications, servers, or storage media performing the operations described in this application.

[0062] As an optional but not limited implementation, in response to a user's active request, sending a prompt message to the user can be done via a pop-up window, where the prompt message can be presented in text format. Furthermore, the pop-up window can also include a selection control allowing the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0063] It is understood that the above notification and user authorization process is merely illustrative and does not limit the implementation of this application. Other methods that comply with relevant laws and regulations may also be applied to the implementation of this application.

[0064] The high-intensity mining of coal resources has led to the development of water-conducting fracture zones up to the top of the bedrock aquifer. In areas where the key aquitard layer of the overlying soil is missing (or weak), the shallow Quaternary aquifer with strong water-bearing capacity directly (overflows) supplies the bedrock aquifer through the water-conducting fracture zone, thus becoming a stable source of mine water supply. This can easily cause roof water disasters and endanger underground safe production.

[0065] Based on the above, a horizontal water-resistant curtain can be constructed by grouting in the vertical water passage of the Quaternary aquifer to recreate the water-resistant layer, cut off the supply of Quaternary groundwater to mine water, and achieve water reduction underground and water retention on the surface.

[0066] This application finds that if horizontal curtain grouting were used to reconstruct the aquitard in all areas of the soil layer where the critical aquitard is missing (weak), it would undoubtedly incur enormous economic costs. On the other hand, due to the complex internal and external dynamic geological processes such as sedimentation, tectonics, and weathering involved in the formation of strata, their water-bearing capacity and water-bearing properties exhibit significant spatial variability. Therefore, how to leverage this characteristic of the strata to quickly and accurately select key areas for aquitard reconstruction is a pressing problem that needs to be solved.

[0067] The embodiments of this application provide a method for determining the reconstruction range of an aquitard. By obtaining the geological parameters of the target mine and determining the water source based on the geological parameters, a spatial location zoning map corresponding to the range of weak soil layers is constructed by combining a triaxial seepage algorithm. Combined with a water-rich zoning map constructed based on a transient electromagnetic algorithm and a goaf water accumulation zoning map, a grouting treatment selection map is continuously superimposed to construct a grouting treatment selection map. Based on the selection index algorithm, a suitable aquitard reconstruction range can be quickly and accurately selected for grouting and water plugging to achieve aquitard reconstruction, so as to treat the aquitard in real time and ensure safe production underground.

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

[0069] Step 101: Obtain the geological parameters of the target mine and determine the water source based on the geological parameters.

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

[0071] Obtaining geological parameters typically requires geological exploration, which includes methods such as geological surveys, geophysical exploration (e.g., seismic exploration, electrical resistivity tomography), and hydrogeological exploration. These exploration activities can provide detailed information about the geological structure of the mining area.

[0072] For example, conducting geological exploration and production status surveys of mines, collecting geological and hydrogeological data of the mines, determining the boundary range of the study area, the stratigraphic structure, aquifer types, hydrogeological characteristics, and the hydraulic connections between aquifers (i.e., geological parameters).

[0073] Water sources refer to the sources of water that may enter a 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 identify which rock strata or geological structures are potential water-bearing channels. For example, highly permeable rock strata or faults may be the main pathways for groundwater flow. In addition, the height of the groundwater level and the direction of groundwater flow are also important factors in determining the source of water recharge.

[0075] For example, in the western region, the coal-bearing strata of this coal mine are the Middle Jurassic Yan'an Formation. The general characteristic of the occurrence of coal seams and aquifers is that coal and water coexist, with the aquifer above and the coal seam below.

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

[0077] The impermeable layer is mainly a Neogene laterite impermeable layer (with missing areas);

[0078] Currently, the main coal seam being mined is the Yan'an Formation 2-2 of the Jurassic system. The water-conducting fracture zone generated by the coal seam has developed to the weathered bedrock aquifer, resulting in groundwater from the weathered bedrock of the coal seam roof and Quaternary unconfined water being the main water sources in areas with weak or missing soil layers.

[0079] Step 102: Based on the geological parameters, determine the range of the weak soil layer using a triaxial seepage algorithm, and construct a spatial location zoning map corresponding to the range of the weak soil layer.

[0080] In this step, the triaxial seepage algorithm considers the permeability characteristics of the soil layer in three-dimensional space, as well as the flow pattern of water within 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, thereby assessing the stability and permeability of the soil layer.

[0081] Using triaxial seepage algorithms, the stress-strain state and water flow path of soil layers under different seepage conditions can be simulated. By analyzing the simulation results, regions in the soil layer with high permeability, high water flow velocity, and poor stress state can be identified. These regions are often considered as weak areas of the soil layer (i.e., the locations of areas where the critical impermeable layer is missing (weak)). The existence of weak areas may lead to soil failure under external loads or water flow, thereby triggering geological disasters or engineering problems.

[0082] After identifying the areas of weakness in the soil, Geographic Information System (GIS) and Computer-Aided Design (CAD) can be used to transform them into an intuitive spatial zoning map. By marking these weak areas on the map and assigning them different colors or symbols to distinguish their severity, a clear and intuitive spatial zoning map can be created, which helps to quickly understand the stability of the soil layer.

[0083] For example, firstly, the distribution of key water-impermeable layers in the soil is investigated using methods such as drilling and geophysical exploration; based on soil samples obtained from boreholes, indoor triaxial seepage experiments and physical and mechanical property tests are conducted to delineate the range of weak areas in the soil layer and draw a spatial location zoning map of the missing (weak) soil layers in the coal seam roof. Figure 2A As shown, numerous areas of missing (weak) laterite distribution were found within the exploration area of ​​the coal mine. The aforementioned zoning map of missing (weak) soil layers in the coal seam roof includes areas of missing (weak) soil layers and areas of intact soil layers. Figure 2A The horizontal axis represents the distance advanced in the direction of the working face, and the vertical axis represents the length in the direction of the working face dip. 5m represents the 5m soil layer thickness contour line, and 10m represents the 10m soil layer thickness contour line. In this example, the 10m contour line is used as the threshold for the weak soil layer.

[0084] Step 103: Based on the water source, determine the distribution range of the horizontal water-rich area using a transient electromagnetic algorithm, construct a horizontal water-rich zoning map corresponding to the distribution range of the horizontal water-rich area, and based on the water source, determine the distribution range of the vertical water-rich area using a transient electromagnetic algorithm, construct a vertical water-rich zoning map corresponding to the distribution range of the vertical water-rich area.

[0085] In this step, transient electromagnetic method (TEM) is suitable for detecting low-resistivity underground bodies, such as aquifers and water-rich areas. This method utilizes 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, thereby analyzing the distribution and water-richness of groundwater.

[0086] Based on a water-bearing source, transient electromagnetic algorithms can accurately identify key information such as the distribution range of underground aquifers, the relative strength of water content, and the spatial morphology of anomalous zones. This is because water-rich areas have a significant electrical difference compared to the surrounding strata, i.e., lower resistivity. Through transient electromagnetic measurements, the apparent resistivity data of the underground medium can be obtained, and then apparent resistivity profiles or slices can be drawn, thus visually displaying the distribution characteristics of water-rich areas.

[0087] Then, by analyzing the variation characteristics of apparent resistivity, the low-resistivity anomaly region, i.e. the water-rich region, was identified.

[0088] Based on 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 anomalous areas. Subsequently, a horizontal water-rich zoning map was drawn based on the distribution characteristics of the water-rich areas.

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

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

[0091] For example, this coal mine is primarily threatened by fourth-level groundwater and weathered bedrock groundwater. Due to the significant spatial heterogeneity in the water-bearing capacity of the weathered bedrock aquifer and the underlying bedrock aquifer, the reconstruction of the coal seam aquitard must fully consider the water-bearing capacity of the weathered bedrock aquifer and its main water-bearing areas. Based on the geological and hydrogeological conditions of the coal mine, transient electromagnetic exploration was used to investigate the water-bearing capacity of the weathered bedrock aquifer and the underlying bedrock aquifer. Based on the low-resistivity anomaly zone of the weathered bedrock, a water-bearing capacity zoning map of the weathered bedrock aquifer (i.e., a horizontal water-bearing zoning map) was identified, such as... Figure 2B As shown, three low-resistivity anomaly zones (water-rich areas) were found in the weathered bedrock aquifer within the exploration area of ​​the coal mine. The horizontal water-rich zoning map of the weathered bedrock aquifer in the coal seam roof includes strongly water-rich areas and weakly water-rich areas, among which... Figure 2B The horizontal axis represents the distance advanced in the direction of the working face, the vertical axis represents the length in the direction of the working face dip, and the gray text 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 subsurface medium is acquired, and the collected data is processed through filtering and time-depth conversion to eliminate noise and extract useful information. Then, an inversion algorithm is used to image the processed data to obtain a resistivity distribution map of the subsurface medium. Subsequently, based on the resistivity distribution map, low-resistivity anomaly zones are identified. These areas typically correspond to water-rich areas; therefore, by analyzing resistivity data at different depths, the distribution range of vertically water-rich areas can be determined.

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

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

[0095] Using a Geographic Information System (GIS) or related mapping software, each water-rich zone is plotted on a map to form a vertical water-rich zone map.

[0096] For example, transient electromagnetic exploration methods can be used to investigate the water-bearing distribution of each aquifer on the vertical profile of the coal seam roof, identify the distribution and changes of groundwater drawdown cones, draw a vertical groundwater water-bearing zoning map of the coal seam roof (i.e., a vertical water-bearing zoning map), and locate the water passages for Quaternary groundwater to enter the mine. Figure 2C As shown, two water-passing channels were found within the exploration area of ​​the coal mine. The vertical groundwater zoning map of the coal seam roof includes vertically water-conducting zones and vertically non-water-conducting zones, among which... Figure 2C The horizontal axis represents the distance advanced in the direction of the working face, and the vertical axis represents the elevation of the top plate of the working face at different positions.

[0097] Step 104: Determine the water inflow at the underground drainage point based on the geological parameters, determine the distribution range of water accumulation in the goaf based on the water inflow at the underground drainage point, and construct 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, a hydrogeological model can be used to estimate the water inflow at each drainage point underground (usually a pre-installed facility for extracting water from the mine). The water inflow refers to the amount of water flowing into or out of the drainage point within a given time period.

[0099] Goaf areas refer to areas in mines that have been completely mined out and abandoned. These areas are prone to water accumulation because they have lost their original rock support.

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

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

[0102] In particular, zoning maps typically use different colors or symbols to represent areas with different levels of water accumulation, so as to intuitively understand the water accumulation situation in the goaf and thus formulate effective drainage and safety management measures.

[0103] For example, by analyzing the water inflow at underground drainage points and conducting hydrogeological and geophysical exploration, the location, extent, and scale of water accumulation in the underground goaf of the coal mine can be determined, and a zoning map of the water accumulation in the underground goaf can be drawn. Figure 2D As shown, three waterlogged areas were found underground in the coal mine. The aforementioned underground goaf waterlogging map includes waterlogged areas and non-waterlogged areas, among which... Figure 2DThe horizontal axis represents the distance the working face advances in the direction of the strike, the vertical axis represents the length of the working face in the dip direction, and the gray shaded area represents the extent of the existing goaf and water accumulation area.

[0104] Step 105: Overlay the spatial location 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 zone selection map.

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

[0106] For example, based on GIS composite functions, by overlaying spatial location zoning maps, horizontal water-rich zoning maps, vertical water-rich zoning maps, and goaf water accumulation zoning maps, a water-conserving coal mining roof grouting treatment selection map (i.e. grouting treatment selection map) can be generated.

[0107] During coal mining, the water accumulation in the mining area and the mining face can be monitored in real time by repeatedly performing steps 102, 103, 104, and 105. This allows for the identification of areas with missing (weak) soil layers, water-rich areas in weathered bedrock, groundwater funnel areas, and water distribution areas in underground goafs. Utilizing... Figure 2A , Figure 2B , Figure 2C and Figure 2D The selected area map for grouting treatment of water-retaining coal mining roof is generated by overlaying, as shown in the figure. Figure 2E As shown, where, Figure 2E The horizontal axis represents the distance advanced in the direction of the working face, and 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 vertical elevation of different positions of the working face roof (Figure 3 from top to bottom). The circled areas represent the distribution of areas with strong water-rich groundwater in the vertical direction of the coal seam roof, that is, the distribution of vertical abnormal funnels (water passages).

[0108] Step 106: Process the grouting treatment selection map using the selection area index algorithm to determine the reconstruction range of the waterproof layer from the grouting treatment selection map.

[0109] In this step, the selection index algorithm is used to evaluate, select, and determine specific areas (such as grouting treatment areas). Based on various geological, engineering, and environmental factors, the algorithm calculates and analyzes to derive one or more indices or indicators to guide decision-making. These indices or indicators reflect information regarding the feasibility, effectiveness, and cost of grouting treatment.

[0110] When processing grouting treatment zoning maps, the zoning index algorithm is used to evaluate and compare various regions on the map. The grouting treatment index or indicator for each region is calculated based on the zoning index algorithm.

[0111] The calculation results are analyzed and compared to determine which areas are suitable for grouting treatment, as well as the priority and plan for grouting treatment.

[0112] Based on geological conditions and grouting treatment requirements, the specific scope of the aquitard reconstruction is determined, including its length, width, and depth.

[0113] Based on the scope and geological conditions of the aquitard reconstruction, a detailed reconstruction plan is formulated, including the selection of grouting materials, the determination of grouting pressure, and the construction of grouting holes.

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

[0115] For example, water-conducting channels in weathered bedrock aquifers are selected as grouting and sealing areas, and specific grouting treatment plans are determined. Suitable grouting materials are selected and injected into the selected grouting and sealing areas to rebuild the water-resistant layer of the coal seam roof, achieving treatment as mining progresses.

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

[0117] Using the composite overlay raster calculation function of GIS, and based on the selection area index model, the generated selection area map for grouting treatment of water-retaining coal mining roof is processed to determine the reconstruction range of the water-retaining layer.

[0118] The above scheme obtains the geological parameters of the target mine and determines the water source based on these parameters. Combining the triaxial seepage algorithm, a spatial location zoning map corresponding to the weak soil layer is constructed. This is combined with a water-rich zoning map constructed based on the transient electromagnetic algorithm and a goaf water accumulation zoning map. By continuously superimposing these maps, a grouting treatment selection map is constructed. Based on the selection index algorithm, a suitable water-resistant layer reconstruction range can be quickly and accurately selected for grouting and water plugging to achieve water-resistant layer reconstruction. This allows for real-time treatment as mining progresses, ensuring safe underground production.

[0119] In some embodiments, step 102, determining the extent of the weak soil zone based on the geological parameters using a triaxial seepage algorithm, includes:

[0120] Step A1: Obtain formation samples from the target mine.

[0121] Step A2: Record the permeability parameters of the formation sample under a preset pressure, including the permeability coefficient.

[0122] Step A3: The area in the stratigraphic region corresponding to the stratigraphic sample with a permeability coefficient greater than a preset permeability coefficient threshold is defined as the range of the weak soil layer.

[0123] In the above scheme, formation samples are collected from the target mine. These samples represent the physical and chemical properties of the formations within the mine and form the basis for subsequent analysis.

[0124] Under preset pressure conditions, these formation samples were subjected to laboratory tests to record and measure their permeability parameters.

[0125] The main permeability parameter includes the key indicator, the permeability coefficient. The permeability coefficient is a physical quantity that measures how easily a fluid passes through a porous medium (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, based on the permeability coefficient, areas in the formation corresponding to the formation sample with a permeability coefficient greater than a preset permeability coefficient threshold are identified as weak soil zones.

[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 permeation and erosion by fluids (such as water and air), which may become potential safety hazards or mining problems.

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

[0129] In some embodiments, step 103, determining the distribution range of the horizontally rich water zone based on the water source using a transient electromagnetic algorithm, includes:

[0130] Step B1: Determine the target geological region corresponding to the water source.

[0131] Step B2: Determine the induced eddy current generated in the target stratum region under the horizontal pulsed magnetic field, and determine the resistivity of the target stratum region based on the induced eddy current.

[0132] Step B3: The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the horizontally water-rich zone.

[0133] In the above approach, geological layers or stratigraphic regions that may contain groundwater are identified. This is typically based on preliminary results from geological surveys, geological maps, previous hydrogeological data, or geophysical exploration. The goal is to narrow the search area, focusing on those strata most likely to contain groundwater.

[0134] By applying a horizontally pulsed magnetic field at or near the Earth's surface, this field penetrates the ground and induces currents (eddy currents) in conductive strata. The size and distribution of these eddy currents depend on the conductivity of the strata (i.e., the reciprocal of its resistivity).

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

[0136] The resistivity of a formation can be inferred by measuring and analyzing the eddy currents induced by a pulsed magnetic field. Resistivity is a measure of how much a material impedes the flow of electric current. For formations saturated with groundwater, the resistivity is usually low due to the conductivity of water. Forward modeling and inverse analysis of the electromagnetic field are used to extract resistivity information from the observed data.

[0137] Finally, based on the resistivity measurement results, areas with resistivity below a preset threshold were identified as the distribution range of horizontally water-rich areas. This threshold was set based on an understanding of local geological and hydrogeological conditions and an analysis of resistivity data from similar areas.

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

[0139] In some embodiments, step 103, determining the distribution range of the vertically rich water zone based on the water source using a transient electromagnetic algorithm, includes:

[0140] Step C1: Determine the target geological region corresponding to the water source.

[0141] Step C2: Determine the induced eddy current generated in the target stratum region under the vertical pulsed magnetic field, and determine the resistivity of the target stratum region based on the induced eddy current.

[0142] Step C3: The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the vertical water-rich zone.

[0143] In the above approach, geological layers or stratigraphic regions that may contain groundwater are identified. This is typically based on preliminary results from geological surveys, geological maps, previous hydrogeological data, or geophysical exploration. The goal is to narrow the search area, focusing on those strata most likely to contain groundwater.

[0144] By applying a horizontally pulsed magnetic field at or near the Earth's surface, this field penetrates the ground and induces currents (eddy currents) in conductive strata. The size and distribution of these eddy currents depend on the conductivity of the strata (i.e., the reciprocal of its resistivity).

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

[0146] The resistivity of a formation can be inferred by measuring and analyzing the eddy currents induced by a pulsed magnetic field. Resistivity is a measure of how much a material impedes the flow of electric current. For formations saturated with groundwater, the resistivity is usually low due to the conductivity of water. Forward modeling and inverse analysis of the electromagnetic field are used to extract resistivity information from the observed data.

[0147] Finally, based on the resistivity measurement results, areas with resistivity below a preset threshold were identified as the distribution range of vertically water-rich zones. This threshold was set based on an understanding of local geological and hydrogeological conditions and an analysis of resistivity data from similar areas.

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

[0149] In some embodiments, step 104, determining the water inflow at the downhole drainage point based on the geological parameters, includes:

[0150] Step D1: Determine the hydraulic gradient based on the geological parameters.

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

[0152] Step D3: Multiply the hydraulic gradient, the permeability coefficient, and the cross-sectional area of ​​the water flow to obtain the water inflow at the downhole drainage point.

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

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

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

[0156] The permeability coefficient (also known as hydraulic conductivity) is a physical quantity that describes the ability of soil or rock to allow water to pass through. It represents 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 through which water flows (such as alleyways, cracks, pipes, etc.).

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

[0159] The product of the hydraulic gradient, the permeability coefficient, and the cross-sectional area of ​​the water flow is used to obtain the water inflow at the downhole drainage point:

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

[0161] Specifically, the flow rate (Q) can be calculated using 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 flow.

[0162] This formula can be used to convert geological and hydrological parameters into specific water inflow values, thereby enabling the rapid and accurate determination of the water inflow at the well drainage point.

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

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

[0165] Step E2: Based on the apparent resistivity value, perform imaging processing to obtain a resistivity distribution map.

[0166] Step E3: The region in the resistivity distribution map whose apparent resistivity value is less than the preset apparent resistivity threshold is designated as a low-resistivity region.

[0167] Step E4: Determine the distribution range of water accumulation in the goaf from the low-resistivity area using the water inflow from the downhole drainage point.

[0168] In the above scheme, high-density electrical resistivity tomography is a geophysical exploration method that uses dense arrays of electrodes arranged on or under the surface to measure the resistivity changes between different electrode combinations in order to infer the electrical structure of the subsurface medium.

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

[0170] Based on the acquired apparent resistivity values, a resistivity distribution map is generated through processing using a specific imaging algorithm (such as an inversion algorithm).

[0171] Resistivity distribution maps can intuitively show the resistivity distribution in different underground areas and are the basis for subsequent analysis.

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

[0173] The area in the diagram with an apparent resistivity value less than this threshold is considered a low-resistivity area. In mine exploration, low-resistivity areas are often associated with geological phenomena such as aquifers and water accumulation in goaf areas.

[0174] After identifying the low-resistivity regions, these regions were further analyzed using known information on water inflow from downhole drainage points.

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

[0176] In summary, this process, through high-density electrical algorithms to collect data, image processing, setting thresholds to determine low-resistivity areas, and finally combining water inflow information, can quickly and accurately determine the distribution range of water accumulation in goaf areas, providing important technical support for safe production and water accumulation control in mines.

[0177] In some embodiments, step 106 includes:

[0178] Step F1: Take any partition map of the grouting treatment selection map as the target partition map. Each target partition map includes multiple spatial locations. Take any spatial location in the target partition map as the target spatial location.

[0179] Step F2: Based on the target spatial location and the preset weight corresponding to the target spatial location, the governance selection index is determined using the following formula:

[0180]

[0181] Where EI represents the Governance Constituency Index, W i (x, y) represents the preset weights corresponding to the spatial locations x and y in the i-th target partition map, fi (x, y) represents the values ​​at target spatial positions x and 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: The target zoning map whose governance selection area index is greater than or equal to the preset index threshold is used as the reconstruction range of the waterproof layer.

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

[0184] These partition maps are further subdivided into multiple spatial locations, each representing a specific point or region within that partition.

[0185] In this process, any partition map is selected as the target partition map for analysis. Within the target partition map, any spatial location is then selected as the target spatial location for further calculations.

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

[0187]

[0188] Where EI represents the Governance Constituency Index, W i (x, y) represents the preset weights corresponding to the spatial locations x and y in the i-th target partition map, f i (x, y) represents the values ​​at target spatial positions x and 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 location may be based on a variety of factors such as geological conditions, grouting requirements, and cost-effectiveness.

[0190] The governance selection index is calculated by taking into account 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 selection index of a target zoning map is greater than or equal to a preset index threshold, then the zoning map is considered a suitable area for waterproofing layer reconstruction.

[0193] For example, areas with a treatment selection index of 4 are designated as the treatment zone for water-retaining coal mining roof grouting (i.e., the area for water-retaining layer reconstruction), while areas with a treatment selection index less than 4 do not require treatment. Furthermore, the treatment zone can be appropriately expanded based on the actual engineering situation to compensate for calculation errors caused by the accuracy of various exploration methods.

[0194] By comprehensively considering the weight and location information of spatial locations in different partition maps, it is possible to quickly and accurately determine which areas most need grouting treatment in order to rebuild the waterproof layer.

[0195] It should be noted that the method in this embodiment can be executed by a single device, such as a computer or server. The method can also be applied in a distributed scenario, where multiple devices cooperate to complete the task. In such a distributed scenario, one of these devices may execute only one or more steps of the method in this embodiment, and the multiple devices will interact with each other to complete the method described.

[0196] It should be noted that the above description describes some embodiments of this 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 a different order than that shown in the above embodiments and still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require a specific or sequential order to achieve the desired result. 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 embodiments, this application also provides a device for determining the reconstruction range of a waterproof layer.

[0198] refer to Figure 3 The device for determining the reconstruction range of the waterproof layer includes:

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

[0200] The first partition map construction module 302 is configured to determine the range of the weak soil layer based on the geological parameters using a triaxial seepage algorithm, and to construct a spatial location partition map corresponding to the range of the weak soil layer.

[0201] The second partition map construction module 303 is configured to determine the distribution range of the horizontal water-rich area based on the water-filling water source using a transient electromagnetic algorithm, construct a horizontal water-rich partition map 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-filling water source using a transient electromagnetic algorithm, construct a vertical water-rich partition map 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 based on 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 zone map construction module 305 is configured to overlay the spatial location zoning map, the horizontal water-rich zoning map, the vertical water-rich zoning map, and the goaf water accumulation zoning map to construct the grouting treatment zone map;

[0204] The reconstruction range determination module 306 is configured to process the grouting treatment selection map using a selection area index algorithm to determine the reconstruction range of the waterproof layer from the grouting treatment selection map.

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

[0206] Obtain formation samples from the target mine;

[0207] Record the permeability parameters of the formation sample under a preset pressure, the permeability parameters including the permeability coefficient;

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

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

[0210] Determine the target geological region corresponding to the water source;

[0211] The induced eddy current generated in the target stratum region under a horizontal pulsed magnetic field is determined, and the resistivity of the target stratum region is determined based on the induced eddy current;

[0212] The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the horizontally water-rich zone.

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

[0214] Determine the target geological region corresponding to the water source;

[0215] The induced eddy current generated in the target stratum region under a vertical pulsed magnetic field is determined, and the resistivity of the target stratum region is determined based on the induced eddy current;

[0216] The region with resistivity less than a preset resistivity threshold in the target stratum area is defined as the distribution range of the vertical water-rich zone.

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

[0218] The hydraulic gradient is determined based on the geological parameters.

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

[0220] The water inflow at the downhole drainage point is obtained by multiplying the hydraulic gradient, the permeability coefficient, and the cross-sectional area of ​​the water flow.

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

[0222] Apparent resistivity values ​​at different depths underground in the target mine were collected using a high-density electrical resistivity algorithm.

[0223] Imaging processing is performed based on the apparent resistivity value to obtain a resistivity distribution map;

[0224] The region in the resistivity distribution map where the apparent resistivity value is less than a preset apparent resistivity threshold is defined as a low-resistivity region.

[0225] The distribution range of water accumulation in the goaf is determined from the low-resistivity area by using the water inflow rate at the underground drainage point.

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

[0227] Any partition map of the grouting treatment selection map is taken as the target partition map. Each target partition map includes multiple spatial locations. Any spatial location in the target partition map is taken as the target spatial location.

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

[0229]

[0230] Where EI represents the Governance Constituency Index, W i (x, y) represents the preset weights corresponding to the spatial locations x and y in the i-th target partition map, f i (x, y) represents the values ​​at target spatial positions x and 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 zoning map whose governance selection area index is greater than or equal to the preset index threshold is used as the reconstruction range of the waterproof layer.

[0232] For ease of description, the above devices are described in terms of function, divided into various modules. Of course, in implementing this application, the functions of each module can be implemented in one or more software and / or hardware.

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

[0234] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this 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 the processor executes the program to implement the method for determining the reconstruction range of the waterproof layer as described in any of the above embodiments.

[0235] Figure 4 This illustration shows a more specific hardware structure diagram of an electronic device provided in this embodiment. 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, memory 402, input / output interface 403, and communication interface 404 are interconnected internally via the bus 405.

[0236] The processor 401 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, 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 the operating system and other applications. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program code is stored in the memory 402 and is called and executed by the processor 401.

[0238] Input / output interface 403 is used to connect input / output modules to realize information input and output. Input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Input devices may include keyboards, mice, touch screens, microphones, various sensors, etc., and output devices may include displays, speakers, vibrators, indicator lights, etc.

[0239] Communication interface 404 is used to connect a communication module (not shown in the figure) to enable communication between this device and other devices. The communication module can communicate via wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0240] Bus 405 includes a pathway for transmitting information between various components of the device (e.g., processor 401, memory 402, input / output interface 403, and communication interface 404).

[0241] It should be noted that although the above-described device only shows the processor 401, memory 402, input / output interface 403, communication interface 404, and bus 405, in specific implementations, the device may also include other components necessary for normal operation. Furthermore, those skilled in the art will understand that the above-described device may only include the components necessary for implementing the embodiments of this specification, and not necessarily all the components shown in the figures.

[0242] The electronic devices described above are used to implement the corresponding water-resistant layer reconstruction range determination method in any of the foregoing embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0243] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this application also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the method for determining the water-resistant layer reconstruction range 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, 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 technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by a computing device.

[0245] The computer instructions stored in the storage medium of the above embodiments are used to cause the computer to execute the method for determining the reconstruction range of the waterproof 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] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of this application is limited to these examples; under the concept of this application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of the embodiments of this application as described above, which are not provided in detail for the sake of brevity.

[0247] Additionally, to simplify the description and discussion, and to avoid obscuring the embodiments of this application, the well-known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. Furthermore, the apparatus may be shown in block diagram form to avoid obscuring the embodiments of this application, and this also takes into account the fact that the details of the implementation of these block diagram apparatuses are highly dependent on the platform on which the embodiments of this application will be implemented (i.e., these details should be fully understood by those skilled in the art). While specific details (e.g., circuits) have been set forth to describe exemplary embodiments of this application, it will be apparent to those skilled in the art that the embodiments of this application can be implemented without these specific details or with variations thereof. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0248] Although this application has been described in conjunction with specific embodiments thereof, many substitutions, 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 be used with the embodiments discussed.

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

Claims

1. A method for determining a reconstruction range of a water layer, characterized by, The method comprises: obtaining geological parameters of a target mine, and determining a roof water filling source based on the geological parameters, the roof water filling source comprising a quaternary aquifer recharging to a bedrock aquifer through a water conducting fractured zone; determining a roof soil layer weak zone range by a triaxial seepage algorithm based on the geological parameters, and constructing a spatial position zoning map corresponding to the roof soil layer weak zone range; determining a distribution range of a roof horizontal water enrichment zone by a transient electromagnetic algorithm based on the roof water filling source, constructing a horizontal water enrichment zoning map corresponding to the distribution range of the roof horizontal water enrichment zone, and determining a distribution range of a roof vertical water enrichment zone by the transient electromagnetic algorithm based on the roof water filling source, and constructing a vertical water enrichment zoning map corresponding to the distribution range of the roof vertical water enrichment zone; determining a water inflow of a downhole drainage point according to the geological parameters, determining a goaf water distribution range based on the water inflow of the downhole drainage point, and constructing a goaf water zoning map corresponding to the goaf water distribution range; superimposing the spatial position zoning map, the horizontal water enrichment zoning map, the vertical water enrichment zoning map and the goaf water zoning map to construct a roof grouting treatment selection zone map; processing the roof grouting treatment selection zone map by a selection index algorithm to determine a water-resisting layer reconstruction range from the roof grouting treatment selection zone map; the method for determining the goaf water distribution range based on the water inflow of the downhole drainage point comprises: collecting apparent resistivity values of different depths of the target mine by a high-density electrical method; performing imaging processing based on the apparent resistivity values to obtain a resistivity distribution map; regarding an area with an apparent resistivity value less than a preset apparent resistivity threshold value in the resistivity distribution map as a low-resistance area; determining the goaf water distribution range from the low-resistance area by using the water inflow of the downhole drainage point; the method for processing the roof grouting treatment selection zone map by the selection index algorithm to determine the water-resisting layer reconstruction range from the roof grouting treatment selection zone map comprises: regarding any zoning map of the roof grouting treatment selection zone map as a target zoning map, each target zoning map comprising a plurality of spatial positions, and regarding any spatial position in the target zoning map as a target spatial position; determining a treatment selection index based on the target spatial position and a preset weight corresponding to the target spatial position by the following formula: ； wherein, represents a governance constituency index, represents a preset weight corresponding to a target spatial position in the i-th target partition map, represents a preset weight corresponding to a target spatial position in the i-th target partition map, represents a value based on the spatial position partition map, the horizontal water-rich partition map, the vertical water-rich partition map, and the goaf water accumulation partition map at a target spatial position in the i-th target partition map. regarding a target zoning map with the treatment selection index greater than or equal to a preset index threshold value as the water-resisting layer reconstruction range.

2. The method of claim 1, wherein, the method for determining the roof soil layer weak zone range based on the geological parameters by the triaxial seepage algorithm comprises: obtaining a stratum sample of the target mine; recording a permeability parameter of the stratum sample under a preset pressure, the permeability parameter comprising a permeability coefficient; regarding an area with a permeability coefficient greater than a preset permeability coefficient threshold value in a stratum area corresponding to the stratum sample as the roof soil layer weak zone range.

3. The method of claim 1, wherein, the method for determining the distribution range of the roof horizontal water enrichment zone based on the roof water filling source by the transient electromagnetic algorithm comprises: determining a target stratum area corresponding to the roof water filling source; determining induced eddy current generated by the target stratum region under horizontal direction pulsed magnetic field, determining the resistivity of the target stratum region based on the induced eddy current; regarding the region of the resistivity of the target stratum region less than a preset resistivity threshold as the distribution range of the horizontal water enrichment zone of the roof.

4. The method of claim 1, wherein, the distribution range of the vertical water enrichment zone of the roof is determined by the transient electromagnetic algorithm based on the water source of the roof, determining the target stratum region corresponding to the water source of the roof; determining induced eddy current generated by the target stratum region under vertical direction pulsed magnetic field, determining the resistivity of the target stratum region based on the induced eddy current; regarding the region of the resistivity of the target stratum region less than a preset resistivity threshold as the distribution range of the vertical water enrichment zone of the roof.

5. The method of claim 1, wherein, the water inflow of the underground drainage point is determined according to the geological parameters, determining the hydraulic gradient according to the geological parameters; determining the permeability coefficient and the water flow section area of the target mine; the product of the hydraulic gradient, the permeability coefficient and the water flow section area is processed to obtain the water inflow of the underground drainage point.

6. A device for determining the reconstruction range of a waterproof layer, characterized in that, comprising: the water source determination module is configured to obtain the geological parameters of the target mine, and determine the water source of the roof based on the geological parameters, wherein the water source of the roof includes the fourth system aquifer supplying to the bedrock aquifer through the water conducting fractured zone; the first partition map construction module is configured to determine the range of the roof soil weak zone by the triaxial seepage algorithm based on the geological parameters, and construct the spatial position partition map corresponding to the range of the roof soil weak zone; the second partition map construction module is configured to determine the distribution range of the horizontal water enrichment zone of the roof by the transient electromagnetic algorithm based on the water source of the roof, construct the horizontal water enrichment partition map corresponding to the distribution range of the horizontal water enrichment zone of the roof, and determine the distribution range of the vertical water enrichment zone of the roof by the transient electromagnetic algorithm based on the water source of the roof, and construct the vertical water enrichment partition map corresponding to the distribution range of the vertical water enrichment zone of the roof; the 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 the goaf water based on the water inflow of the underground drainage point, and construct the goaf water partition map corresponding to the distribution range of the goaf water; the treatment selection area map construction module is configured to superimpose process the spatial position partition map, the horizontal water enrichment partition map, the vertical water enrichment partition map and the goaf water partition map, and construct the roof grouting treatment selection area map; the reconstruction range determination module is configured to process the roof grouting treatment selection area map by the selection area index algorithm, and determine the water-resisting layer reconstruction range from the roof grouting treatment selection area map; the third partition map construction module is specifically configured to: collect the apparent resistivity value of different depths of the target mine by the high density electrical algorithm; image processing is performed based on the apparent resistivity value to obtain the resistivity distribution map; regard the region of the apparent resistivity value less than a preset apparent resistivity threshold in the resistivity distribution map as a low resistance region; determine the distribution range of the goaf water from the low resistance region by using the water inflow of the underground drainage point; The reconstruction range determination module is specifically configured to: take any subarea map of the roof grouting treatment selection area map as a target subarea map, each target subarea map including a plurality of spatial positions, and take any spatial position in the target subarea map as a target spatial position; determine a treatment selection area index based on the target spatial position and a preset weight corresponding to the target spatial position according to the following formula: ； wherein, denotes a governance district index, denotes a preset weight corresponding to a target spatial position in the i-th target district map, denotes a preset weight corresponding to a target spatial position in the i-th target district map, denotes a value based on the spatial position district map, the horizontal water-rich district map, the vertical water-rich district map, and the goaf water accumulation district map at a target spatial position in the i-th target district map. take a target subarea map with a treatment selection area index greater than or equal to a preset index threshold as the water-resisting layer reconstruction range.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and running on the processor, characterized in that, The processor implements the method of any one of claims 1 to 5 when executing the program.

8. A non-transitory computer-readable storage medium storing computer instructions, wherein, The computer instructions are used to make the computer execute the method of any one of claims 1 to 5.

Images