A mine roadway exploration method, system, device and medium

Abstract

The application relates to a mine roadway exploration method, system, device and medium, the method comprising: using an equivalent counter-flux transient electromagnetic method to explore a to-be-explored area of a mine, and obtaining abnormal resistivity distribution; using a high-density resistivity method to determine profile resistivity distribution based on the abnormal resistivity distribution; determining a preliminary roadway position based on the abnormal resistivity distribution and the profile resistivity distribution; verifying the preliminary roadway position to obtain a target roadway position; and drilling and scanning the target roadway position to obtain target survey information. The application solves the problem that, in the prior art, a roadway position determined by directly surveying a surface environment of a roadway with complex physical characteristics is often relatively large in error, and source treatment is relatively difficult and low in efficiency.

Description

Technical Field

[0001] This application relates to the field of mining exploration technology, and in particular to a method, system, equipment and medium for mining roadway exploration. Background Technology

[0002] Abandoned mines, such as pyrite and coal mines, generate large amounts of acidic underground wastewater, which transforms into surface water at tunnel exits, adit exits, or shafts, discharging into rivers or drinking water sources. This water quality poses a significant threat to the surrounding environment and human health. Source control is an effective method for preventing and controlling the problem in such mines, and the key to source control lies in identifying the underground acidic water flow channels. Current technologies often involve directly surveying the surface environment of abandoned mines to determine the location of underground channels.

[0003] However, due to the long-term changes in the natural environment, there are more or less sinkholes on the surface of the roadways at the top of the roadways, which can cause rainwater to flow in directly, often forming old mine water. This makes the physical characteristics of the roadways more complex. The roadway locations determined by directly surveying the surface environment in the existing technology are often inaccurate, resulting in greater difficulty and lower efficiency in source control. Summary of the Invention

[0004] To overcome the problem that the location of a mine roadway is often determined by directly surveying the surface environment of a roadway with complex physical properties in the existing technology, which often results in large errors and makes source control difficult and inefficient, this application provides a mine roadway exploration method, system, equipment and medium.

[0005] Firstly, in order to solve the aforementioned technical problems, this application provides a method for mine roadway exploration, comprising: The equivalent reverse flux transient electromagnetic method was used to explore the area to be explored in the mine and obtain the distribution of abnormal resistivity. The high-density resistivity method is used to determine the profile resistivity distribution based on the abnormal resistivity distribution. Based on the abnormal resistivity distribution and the profile resistivity distribution, the initial location of the roadway was determined; The initial location of the tunnel was verified to obtain the location of the target tunnel; Drilling and scanning were performed on the target tunnel location to obtain target survey information.

[0006] Secondly, this application also provides a mine roadway exploration system, comprising: The first acquisition module is used to explore the area to be explored in the mine using the equivalent reverse flux transient electromagnetic method to obtain the abnormal resistivity distribution. The second acquisition module is used to determine the profile resistivity distribution based on the abnormal resistivity distribution using the high-density resistivity method. The determination module is used to determine the initial location of the roadway based on the abnormal resistivity distribution and the profile resistivity distribution. The verification module is used to verify the initial tunnel location and obtain the target tunnel location; The module is used to drill and scan the target tunnel location to obtain target survey information.

[0007] Thirdly, this application also provides a computing device, including a memory, a processor, and a program stored in the memory and running on the processor, wherein the processor executes the program to implement the steps of a mine roadway exploration method as described above.

[0008] Fourthly, this application also provides a computer-readable storage medium storing instructions that, when executed on a terminal device, cause the terminal device to perform the steps of a mine roadway exploration method.

[0009] The beneficial effects of this application are as follows: By using the equivalent reverse flux transient electromagnetic method and the high-density resistivity method, the abnormal resistivity distribution and profile resistivity distribution of the area to be explored in the mine can be obtained. This allows for understanding the resistivity characteristics of roadways with complex physical properties that may exist in the area to be explored, thus directly determining the possible initial roadway locations based on the two distributions. The initial roadway locations are then verified to obtain target roadway locations with satisfactory accuracy. Precise drilling and scanning are then performed on the target roadway locations to obtain target survey information near the roadway. In this way, by repeatedly determining and verifying the locations of roadways with complex physical properties in the area to be explored through multiple methods, precise roadway exploration is achieved. This improves the accuracy of the obtained target roadway locations, reduces the workload of drilling and scanning, and shortens the time required to obtain accurate target survey information for the roadways. This reduces the time required for roadway source management based on the target survey information, thereby reducing the difficulty of roadway source management and improving management efficiency. Attached Figure Description

[0010] Figure 1 This is a schematic flowchart illustrating a mine roadway exploration method as an exemplary embodiment of this application; Figure 2 This is a map showing the distribution of geological survey points obtained by implementing a mine tunnel exploration method in an exemplary embodiment of this application; Figure 3 This is a regional distribution map of another geological condition obtained by implementing a mine tunnel exploration method in an exemplary embodiment of this application; Figure 4 A comparison chart of test results distribution for different apparent resistivity testing devices; Figure 5This is a schematic flowchart illustrating the application of the provided mine roadway exploration method in an exemplary embodiment of this application; Figure 6 This is a schematic diagram of the structure of a mine roadway exploration system, which is an exemplary embodiment of this application. Detailed Implementation

[0011] Due to long-term natural environmental changes, the lithological and physical properties of the local strata in abandoned mines have changed, making the physical characteristics of the mine tunnels quite complex. Based on the strong electrical differences between different strata and the hydrogeological structure of the area, the lithological and physical properties of the local strata are analyzed. The metamorphic rock group is fragmented and incomplete due to the destruction by later granites and gabbro-diabases. The main lithology is graphite-bearing biotite plagioclase gneiss, with relatively high resistivity, reaching over 800 ohm-meters. Quaternary residual slope deposits are mainly distributed in valleys, consisting primarily of gravelly soil, with resistivity typically between 50 and 100 ohm-meters. The granites, mainly biotite granites and pegmatite dikes, exhibit strong high resistivity, reaching over 1000 ohm-meters. In tectonically fractured areas, the resistivity decreases; previous surveys have shown it to drop to around 200 ohm-meters. Meanwhile, the slag, due to the presence of hydrous metallic minerals or highly conductive water, typically has a resistivity less than 50 ohm-meters. Meanwhile, the type of backfill can influence the physical properties of roadways. Unfilled roadways, affected by high-resistivity, low-velocity air, exhibit strong high-resistivity, low-velocity anomalies, while water-filled and soil-filled roadways show low-resistivity, low-velocity anomalies. Furthermore, chemical testing of contaminated water reveals its strong conductivity, indicating that contaminated water backfilling significantly alters the electrical characteristics of roadways. Therefore, abandoned mines often exhibit significant physical property anomalies, and these anomalies provide valuable insights for geophysical exploration data interpretation.

[0012] Based on this, this application utilizes the equivalent reverse flux transient electromagnetic method, high-density resistivity method, micro-motion exploration method, charging method, and borehole ground-penetrating radar scanning to explore the tunnels of abandoned mines. This method can accurately locate abandoned tunnels (2m wide × 2m high) with a burial depth of less than 150m. It can also determine the tunnel survey information, including the development of fractures in the tunnel and its surrounding area, the flow path of underground acidic water, and the bottom boundary of hydrogeological conditions. Based on the tunnel locations and survey information obtained in this application, a hydrogeological conceptual model can be established, improving the feasibility of source control of underground acidic wastewater. This provides a referable technical methodology system for the accurate identification of hydrogeological conditions and ecological restoration of abandoned mines, reducing the difficulty and improving the efficiency of source control.

[0013] The following embodiments are further explanations and supplements to this application and do not constitute any limitation on this application.

[0014] The following describes a method, system, equipment, and medium for mine roadway exploration according to an embodiment of this application, with reference to the accompanying drawings.

[0015] This application discloses a method for exploring mine roadways. The method is applied to a terminal device. This application uses the terminal device as the executing entity to describe the solution. The terminal device is used to execute the steps of a method for exploring mine roadways.

[0016] Please see Figure 1 , Figure 1 A method for mine roadway exploration is illustrated in an exemplary embodiment of this application, such as... Figure 1 As shown, this application provides a method for mine roadway exploration, including: Step S11: The area to be explored in the mine is explored using the equivalent reverse flux transient electromagnetic method to obtain the distribution of abnormal resistivity. Step S12: Using the high-density resistivity method, determine the profile resistivity distribution based on the abnormal resistivity distribution. Step S13: Determine the initial location of the roadway based on the abnormal resistivity distribution and the profile resistivity distribution. Step S14: Verify the initial measured location of the tunnel to obtain the location of the target tunnel; Step S15: Drill and scan the target tunnel location to obtain target survey information.

[0017] This embodiment of a mine roadway exploration method uses the equivalent reverse flux transient electromagnetic method and the high-density resistivity method to obtain the anomalous resistivity distribution and profile resistivity distribution of the area to be explored in the mine. This allows for understanding the resistivity characteristics of roadways with complex physical properties that may exist in the area to be explored, thus directly determining the possible initial roadway locations based on the two distributions. The initial roadway locations are then verified to obtain target roadway locations with sufficient accuracy. Precise drilling and scanning are then performed on the target roadway locations to obtain target exploration information near the roadway. In this way, by repeatedly determining and verifying the locations of roadways with complex physical properties in the area to be explored through multiple methods, precise roadway exploration is achieved. This improves the accuracy of the obtained target roadway locations, reduces the workload of drilling and scanning, and shortens the time required to obtain accurate target exploration information. This reduces the time required for roadway source management based on target exploration information, thereby reducing the difficulty of roadway source management and improving management efficiency.

[0018] In this embodiment, the Time Domain Electromagnetic Method (TEM) is a time-domain electromagnetic exploration method that detects underground geological bodies by receiving changes in the secondary induced field. Currently, the Opposing Coils Transient Electromagnetic Method (OCTEM) has been developed for shallow exploration. It is a novel transient electromagnetic method for measuring the decay and diffusion of the transient electromagnetic field of equivalent reverse magnetic flux. It utilizes the spatiotemporal distribution of the electromagnetic field generated by the equivalent reverse magnetic flux when two identical sets of coils are energized with opposite currents. Two identical sets of coils, parallel and coaxial, are used as emission sources. The pure secondary field coupled to the center of the earth is measured on the zero-flux plane of the primary field synthesized by these two coil sources. A dual-coil system transmits transient pulse electromagnetic field signals on the ground. One set of coils is placed near the ground surface. At the instant the transient pulse is de-energized, the superimposed magnetic field near the ground surface is at its maximum. Therefore, under the same change time, the maximum surface of the induced eddy current is concentrated near the ground surface, and the magnetic field generated by the induced eddy current is the strongest. As the delay of the turn-off interval increases, the induced eddy current on the ground surface gradually decays and generates a new eddy current maximum surface, which gradually spreads towards the depth and edge away from the transmitting coil. This is the "smoke ring effect" of the transient electromagnetic method vividly described by MNNabighian.

[0019] The method described in this application is applicable to tunnel exploration in abandoned metal mines, non-metal mines, and abandoned mines with underground wastewater outflow. It offers high accuracy for tunnel exploration at depths exceeding 150m when tunnel data is lost or its location is unclear. The area to be explored in this application is delineated based on collected mining and geological data, hydrogeological surveys, and the determination of water outlet locations and flow rates. Simultaneously, based on the mining data, geological data, and hydrogeological survey results, geological information such as faults, veins, collected tunnels, water-bearing mine shafts, and slag heaps in the area to be explored can be obtained. Based on this geological information, a map is drawn as follows: Figure 2 As shown and as Figure 3 The map shows the geological distribution of the area to be explored. Figure 2 and Figure 3 The faults, veins, collection tunnels, water-bearing mine shafts, and slag heaps are displayed in different colors.

[0020] Optionally, the equivalent reverse flux transient electromagnetic method is used to explore the area to be explored in the mine to obtain the distribution of abnormal resistivity, including: Using the equivalent reverse flux transient electromagnetic method, the locations of multiple measuring points in the area to be detected in the mine are determined based on a preset distance requirement, and the apparent resistivity of each measuring point is obtained. Based on the transient relaxation inversion method, several abnormal resistivities among multiple apparent resistivities were determined; The distribution of multiple abnormal resistivity values ​​is obtained as the distribution of abnormal resistivity values ​​in the area to be detected.

[0021] In this embodiment, the equivalent reverse flux transient electromagnetic method is used to determine the locations of multiple measuring points in the target area of ​​the mine based on preset distance requirements, thereby achieving uniform surveying of the target area and obtaining the apparent resistivity at each measuring point. Then, based on the transient relaxation inversion method, anomalies between each apparent resistivity and its corresponding inversion result can be identified. Based on multiple anomalies, multiple abnormal resistivities can be determined. That is, when an anomaly exists, the corresponding apparent resistivity is the abnormal resistivity. The distribution of multiple abnormal resistivities is used as the abnormal resistivity distribution of the target area, facilitating the subsequent acquisition of target survey information based on the abnormal resistivity distribution. The area composed of multiple abnormal resistivities is a low-resistivity distribution area.

[0022] In this embodiment, when the electrical distribution of underground rocks is uneven (there are two or more types of rocks or ores with different conductivity) or the surface is uneven, the resistivity obtained by using the method and calculation formula for measuring uniform horizontal earth resistivity is called apparent resistivity, denoted by the symbol ρs, and the unit is the same as resistivity, which is Ω·m. Using the equivalent reverse flux transient electromagnetic method, the measuring points are arranged in a grid at preset distances on the surface of the area to be explored in the mine. This is mainly used to obtain the electrical structure of the area to be explored. At each measuring point, a resistivity detection device is used to perform ground scanning to obtain the apparent resistivity, which helps to clarify the hydrogeological characteristics. The designed grid spacing is 20m × 30m. Due to factors such as terrain and leveling, there are some minor adjustments to the actual measuring points, all within the design range, and the measuring points are uniformly distributed overall. Figure 2 As shown in the figure, the distribution of measuring points for the equivalent back flux transient electromagnetic method is illustrated by the distribution of the black square measuring points. The equivalent back flux transient electromagnetic method implemented in this application has high accuracy in identifying geological electrical conditions above 250m underground.

[0023] Based on the apparent resistivity at each measurement point, the transient relaxation inversion method is used to invert multiple apparent resistivities, identifying the anomalous and normal resistivity of the area to be detected. The distribution of linear anomalous bands and planar anomalous regions formed by multiple anomalous resistivities is then used as the distribution of anomalous resistivity in the area to be detected, forming an intuitive planar distribution map of anomalous points, such as... Figure 3 As shown. Figure 3 In the diagram, yellow and red planar or linear patches represent low resistivity areas, while blue represents high resistivity areas. Low resistivity areas are caused by residual metal veins or underground mining voids; that is, the distribution of low resistivity areas reflects the abnormal resistivity distribution of the area to be detected.

[0024] in addition, Figure 3 The distribution of measurement points for the equivalent reverse flux transient electromagnetic method is also shown. Furthermore, the distribution of measurement points for the high-density resistivity method, the distribution of the survey lines for the micro-motion exploration method, and the distribution of the survey lines for the charging method used subsequently in this application are also shown in [the document / document / etc.]. Figure 2 and Figure 3 It will be displayed in the middle. Figure 2 The paper also shows the distribution of anomalies obtained by subsequent micro-motion exploration and charging methods. Figure 2 and Figure 3 The various survey lines and points are displayed in different colors.

[0025] Optionally, based on the transient relaxation inversion method, several anomalous resistivities among multiple apparent resistivities are determined, including: Multiple apparent resistivitys are inverted using the transient relaxation inversion method to obtain multiple relative resistivitys; For each relative resistivity, obtain the difference between the relative resistivity and the corresponding apparent resistivity; When the absolute value of the difference is greater than the preset abnormal value, the corresponding apparent resistivity is taken as the abnormal resistivity.

[0026] In this embodiment, the transient relaxation inversion method is used to invert multiple apparent resistivitys, obtaining the relative resistivity after inversion for each apparent resistivity. For each relative resistivity, the difference between the relative resistivity and the corresponding apparent resistivity can indicate anomalies. That is, when the absolute value of the difference is greater than a preset anomaly value, it indicates that an anomaly exists between the two. Therefore, when the absolute value of the difference is greater than the preset anomaly value, the corresponding apparent resistivity is abnormal, and this apparent resistivity is regarded as an abnormal resistivity. This facilitates the subsequent determination of the abnormal resistivity distribution based on the abnormal resistivity, and the subsequent determination of target survey information based on the abnormal resistivity distribution.

[0027] In this embodiment, the transient relaxation inversion method is a formula derivation based on the equivalent back magnetic flux theory, an inversion method applicable to the equivalent back magnetic flux transient electromagnetic method. The transient relaxation inversion method is sensitive to karst and mining-induced defects, thus it can invert the "relative resistivity" in actual geoelectric conditions, that is, invert to obtain a "relative resistivity" that more closely matches the actual geoelectric conditions for each apparent resistivity. Since there is a certain difference between the relative resistivity and the corresponding apparent resistivity, the difference between the two can be used to determine whether the apparent resistivity is abnormal, thereby achieving anomaly detection and obtaining a more accurate distribution of tunnel locations (abnormal resistivity distribution) in the area to be detected.

[0028] Optionally, the profile resistivity distribution is determined using the high-density resistivity method based on the anomalous resistivity distribution, including: Based on the abnormal resistivity distribution, abnormal regions were identified. By intersecting the abnormal regions at large angles, high-density resistivity measurement lines can be obtained. The vertical profile of the high-density resistivity survey line was surveyed using the high-density resistivity method, and the resistivity of multiple profiles was obtained. The resistivity distribution of multiple profiles is obtained as the resistivity distribution of the area to be detected.

[0029] In this embodiment, based on the abnormal resistivity distribution, abnormal regions are delineated, and these regions are intersected at large angles to obtain high-density resistivity survey lines. These high-density resistivity survey lines are those with a high probability of containing underground tunnels. The high-density resistivity method is then used to survey the vertical profiles of these high-density resistivity survey lines, obtaining multiple profile resistivity values. The distribution of these multiple profile resistivity values ​​is then used as the profile resistivity distribution of the area to be explored, facilitating subsequent target survey information based on this profile resistivity distribution.

[0030] In an exemplary embodiment, the intersections between different high-density resistivity survey lines represent potential tunnel fractures, such as... Figure 2 As shown, the intersection of high-density resistivity survey line L2 and high-density resistivity survey line L6 may contain tunnel fissures.

[0031] In this embodiment, for the anomalous regions delineated by the anomalous resistivity distribution obtained through exploration using the equivalent reverse flux transient electromagnetic method, a high-density resistivity survey line (e.g., ...) is obtained by intersecting the anomalous regions at a large angle using the high-density resistivity method. Figure 2 As shown, a vertical profile (measured at a large angle to the long axis of the red patch) was surveyed, and a Winner and Schlumberger apparatus were used to jointly survey the measuring points of the vertical profile of the high-density resistivity measuring line. The joint survey involved a preliminary exploration followed by a detailed exploration to improve the accuracy of the joint survey results and further enhance the accuracy of the subsequently obtained target tunnel location. The measuring point spacing for the preliminary exploration was set at 10m, and the measuring point spacing for the detailed exploration was set at 5m. Implementing the high-density resistivity method of this application can determine the apparent resistivity characteristics of tunnels and fault profiles at depths shallower than 150m.

[0032] Conventional high-density resistivity methods rely on single-device field operations, resulting in limitations in acquiring subsurface electrical parameter signals. However, with advancements in big data fusion technology and data processing capabilities in recent years, multi-device joint surveying and inversion techniques have become a trend. Since different types of devices vary in terms of vertical and horizontal resolution, anti-interference capabilities, and timeliness, multi-device joint resistivity inversion technology, through data fusion and joint inversion calculations, can achieve refined detection of subsurface electrical structures, particularly suitable for detailed detection of complex geological structures. Therefore, this application's high-density resistivity method employs a one-time electrode placement, and its intelligent electrode-running working mode greatly improves the applicability of conventional DC resistivity methods. Furthermore, it utilizes joint surveying with various device types, including Winner, Schlumberger, and dipole-dipole, enabling the acquisition of apparent resistivity with different characteristics on the vertical profile of the high-density resistivity survey line. This provides a more comprehensive distribution of resistivity across multiple profiles, achieving refined detection of complex and anomalous geological conditions.

[0033] In an exemplary embodiment of this application, the obtained high-density resistivity survey line is longer than 600m. When using the Winner and Schlumberger devices for joint surveying, both devices are set up at each measuring point, and the apparent resistivity data of both devices are collected simultaneously. The apparent resistivity data collected by both devices are used as the profile resistivity.

[0034] Please see Figure 4 , Figure 4 A comparison chart of the distribution of test results from different apparent resistivity testing devices, such as... Figure 4 The diagram shows the resistivity distributions of multiple profiles measured using only the Wenner instrument, the Schlumberger instrument alone, and a combined survey using both instruments. Yellow and red patchy or linear areas represent low-resistivity regions, while blue areas represent high-resistivity regions. Based on these three distributions, it is clear that the combined survey using both instruments provides a more comprehensive and accurate representation of the resistivity distributions.

[0035] Optionally, based on the abnormal resistivity distribution and the profile resistivity distribution, the initial location of the roadway is determined, including: Based on the abnormal resistivity distribution, abnormal regions were identified. Based on the resistivity distribution of the profile, high-density regions are identified. Obtain the overlapping area between abnormal areas and high-density areas; Based on the overlapping areas, the initial location of the tunnel was obtained.

[0036] In this embodiment, the abnormal area corresponding to the abnormal resistivity distribution is the area where a roadway may exist, and the high-density area corresponding to the profile resistivity distribution is the area where a roadway is more likely to exist. Therefore, the possibility of a roadway existing in the overlapping area between the abnormal area and the high-density area is greater than that of the former two. Thus, the initial roadway location can be obtained based on this overlapping area, which can improve the accuracy of the initial roadway location and reduce the workload of subsequent drilling and scanning. This can shorten the time required to obtain accurate target survey information of the roadway, thereby shortening the time for roadway source management based on target survey information, reducing the difficulty of roadway source management, and improving management efficiency.

[0037] Optionally, the overlapping area includes multiple overlapping locations; the preliminary tunnel location is obtained based on the overlapping locations, including: Determine the sorting order of multiple overlapping positions according to the preset sorting requirements; By connecting multiple overlapping locations in sorted order, the initial location of the tunnel can be obtained.

[0038] In this embodiment, the preset sorting requirements include: selecting one of multiple overlapping positions as the starting position and marking the sequence number of the starting position as 1; finding the first and second closest overlapping positions that are closest to the starting position from the overlapping positions without marking the sequence number; marking the sequence number of the first overlapping position as the sequence number of the starting position minus 1, and marking the sequence number of the second overlapping position as the sequence number of the starting position plus 1; finding the third overlapping position that is closest to the first overlapping position from the overlapping positions without marking the sequence number, marking the sequence number of the third overlapping position as the sequence number of the first overlapping position minus 1, and finding the fourth overlapping position that is closest to the second overlapping position from the overlapping positions without marking the sequence number, marking the sequence number of the fourth overlapping position as the sequence number of the second overlapping position plus 1, until each overlapping position is marked with a sequence number; and sorting the multiple overlapping positions according to their sequence numbers from smallest to largest to obtain the corresponding sorting order. In this embodiment, by connecting multiple overlapping locations according to a sorting order determined based on preset sorting requirements, point-to-line conversion is achieved, which enables the understanding of the predicted direction of the tunnel. Based on this predicted direction, the initial tunnel location is obtained, which facilitates subsequent verification based on the initial tunnel location to obtain a more accurate target tunnel location.

[0039] Optionally, the initial measured location of the tunnel is verified to obtain the location of the target tunnel, including: The wave velocity of the initially measured roadway location was measured using the micro-motion exploration method to obtain the precise roadway location. The target roadway location is obtained by using the charging method to perform electrical measurements on the roadway location using rated power.

[0040] In this embodiment, the initial roadway location is verified using both micro-motion exploration and charging methods, further eliminating useless roadway locations. This ensures that the obtained target roadway location meets the accuracy requirements, facilitating subsequent precise drilling and scanning exploration of the roadway. This reduces the workload of drilling and scanning, thereby shortening the time required to obtain accurate target exploration information of the roadway. Consequently, the time required for roadway source management based on the target exploration information is reduced, thus lowering the difficulty of roadway source management and improving management efficiency. The rated power is greater than the preset power.

[0041] In an exemplary embodiment, when a mining tunnel is filled with water, a continuous water inflow will occur at a certain outcrop in the tunnel, and the water volume will be relatively stable. This phenomenon can be obtained through hydrogeological exploration, indicating that the tunnel is filled. Furthermore, when the tunnel has mined-out areas, is filled with water, or has cavities, its wave velocity will differ from the surrounding geological body. This is the basis for interpreting physical parameters using the micro-motion exploration method. Simultaneously, when the tunnel is filled with water, a large number of metal ions or salt-containing minerals will dissolve in the water, resulting in extremely high electrical conductivity and forming conductive channels along the tunnel and fissures. At this time, if an electrical measurement is performed at the tunnel outlet using a charging method with a power higher than the preset power, a significant electrical response will be observed on the surface potential anomaly detection within a certain distance range (1km to 3km) from the measurement point. Thus, based on the initial tunnel location and combined with the tunnel filling conditions obtained from the hydrogeological survey, this application shows that within a range of 1km to 3km from the tunnel outlet, both the equivalent reverse flux transient electromagnetic method and the high-density resistivity method indicate locations with low resistivity. Utilizing differences in electrical properties and wave velocity properties, the micro-motion exploration method and the charging method can be used for more precise verification exploration, obtaining the target tunnel location with the required accuracy. Specifically, the measurement points for the micro-motion exploration method and the charging method are the same as those used in the detailed exploration with the high-density resistivity method, i.e., the distance between the measurement points for both methods is 5m.

[0042] Please see Figure 5 , Figure 5 This is a schematic flowchart illustrating the application of the provided mine tunnel exploration method in an exemplary embodiment of this application. Figure 5As shown, for precise exploration of roadways in abandoned mines, the first step is to conduct on-site reconnaissance, data collection, and hydrogeological investigation to obtain geological information such as faults, veins, collected roadways, water-bearing tunnels, and slag heaps in the mining area (the area to be explored). Based on this geological information, the mine roadway exploration method of this application is applied, using the equivalent reverse flux transient electromagnetic method to explore the area to be explored, obtaining the abnormal resistivity distribution. Based on the abnormal resistivity distribution, the abnormal areas within the area to be explored can be delineated. Secondly, using the high-density resistivity method, combined with the abnormal resistivity distribution, the initial roadway location and the predicted location of possible roadway fractures are obtained. Then, based on the wave velocity and electrical differences in the roadway, the micro-motion exploration method and the charging method are used to verify the initial roadway location and the predicted location of roadway fractures, determining the target roadway location and the target fracture location related to the water outlet with the required accuracy. Finally, drilling verification was conducted based on the target tunnel and fracture locations to drill channels from the surface to these locations. Ground-penetrating radar was then used to scan these channels, obtaining target survey information within a 5m or 10m radius of the tunnel. This information pinpoints the precise location for subsequent tunnel curtain grouting and filling, allowing for the development of effective remediation measures in advance, thereby improving the efficiency and reducing the difficulty of source control. Anomalies included linear resistivity anomalies and planar anomalies; the initial tunnel and predicted fracture locations included depth information; and the target survey information included the distribution of surrounding tunnels and the development of tunnel fractures.

[0043] Please see Figure 6 , Figure 6 A mine roadway exploration system is illustrated as an exemplary embodiment of this application, such as Figure 6 As shown, this application provides a mine roadway exploration system, comprising: The first acquisition module is used to explore the area to be explored in the mine using the equivalent reverse flux transient electromagnetic method to obtain the abnormal resistivity distribution. The second acquisition module is used to determine the profile resistivity distribution based on the abnormal resistivity distribution using the high-density resistivity method. The determination module is used to determine the initial location of the roadway based on the abnormal resistivity distribution and the profile resistivity distribution. The verification module is used to verify the initial tunnel location and obtain the target tunnel location; The module is used to drill and scan the target tunnel location to obtain target survey information.

[0044] This embodiment of a mine roadway exploration system utilizes a first acquisition module and a second acquisition module to obtain the anomalous resistivity distribution and profile resistivity distribution of the area to be explored in the mine using the equivalent reverse flux transient electromagnetic method and the high-density resistivity method. This allows for understanding the resistivity characteristics of roadways with complex physical properties that may exist in the area to be explored. Based on these two distributions, a determination module can then preliminarily determine the possible initial roadway locations. A verification module then verifies the initial roadway locations to obtain target roadway locations with satisfactory accuracy. The acquisition module then performs precise drilling and scanning on the target roadway locations to obtain target exploration information near the roadway. By repeatedly determining and verifying the locations of roadways with complex physical properties in the area to be explored through multiple methods, precise roadway exploration is achieved. This improves the accuracy of the obtained target roadway locations, reduces the workload of drilling and scanning, and shortens the time required to obtain accurate target exploration information. This reduces the time required for roadway source management based on the target exploration information, thereby reducing the difficulty of roadway source management and improving management efficiency.

[0045] Optionally, the first obtaining module is specifically used for: Using the equivalent reverse flux transient electromagnetic method, the locations of multiple measuring points in the area to be detected in the mine are determined based on a preset distance requirement, and the apparent resistivity of each measuring point is obtained. Based on the transient relaxation inversion method, several abnormal resistivities among multiple apparent resistivities were determined; The distribution of multiple abnormal resistivity values ​​is obtained as the distribution of abnormal resistivity values ​​in the area to be detected.

[0046] Optionally, the first obtaining module is specifically used for: Multiple apparent resistivitys are inverted using the transient relaxation inversion method to obtain multiple relative resistivitys; For each relative resistivity, obtain the difference between the relative resistivity and the corresponding apparent resistivity; When the absolute value of the difference is greater than the preset abnormal value, the corresponding apparent resistivity is taken as the abnormal resistivity.

[0047] Optionally, the second obtaining module is specifically used for: Based on the abnormal resistivity distribution, abnormal regions were identified. By intersecting the abnormal regions at large angles, high-density resistivity measurement lines can be obtained. The vertical profile of the high-density resistivity survey line was surveyed using the high-density resistivity method, and the resistivity of multiple profiles was obtained. The resistivity distribution of multiple profiles is obtained as the resistivity distribution of the area to be detected.

[0048] Optionally, a module is defined, specifically for: Based on the abnormal resistivity distribution, abnormal regions were identified. Based on the resistivity distribution of the profile, high-density regions are identified. Obtain the overlapping area between abnormal areas and high-density areas; Based on the overlapping areas, the initial location of the tunnel was obtained.

[0049] Optionally, the overlapping area includes multiple overlapping locations; the determining module is specifically used for: Determine the sorting order of multiple overlapping positions according to the preset sorting requirements; By connecting multiple overlapping locations in sorted order, the initial location of the tunnel can be obtained.

[0050] Optionally, a module is obtained, specifically for: The wave velocity of the initially measured roadway location was measured using the micro-motion exploration method to obtain the precise roadway location. The target roadway location is obtained by using the charging method to perform electrical measurements on the roadway location using rated power.

[0051] A computing device according to an embodiment of this application includes a memory, a processor, and a program stored in the memory and running on the processor. When the processor executes the program, it implements some or all of the steps of the above-described mine roadway exploration method.

[0052] The computing device can be a computer, and the corresponding program is computer software. The parameters and steps of the computing device described above can be referred to the parameters and steps in the embodiment of the mine roadway exploration method above, and will not be repeated here.

[0053] This application embodiment provides a computer-readable storage medium storing instructions that, when executed, perform the steps of the aforementioned mine roadway exploration method.

[0054] The computer-readable storage medium may be a transient computer-readable storage medium or a non-transitory computer-readable storage medium.

[0055] The technical solutions of this disclosure can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes one or more instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of this disclosure. The aforementioned computer-readable storage medium can be a non-transitory computer-readable storage medium, including: USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks, and other media capable of storing program code; it can also be a transient computer-readable storage medium.

[0056] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. Each block in a flowchart or block diagram may represent a module, segment, or portion of code, which contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0057] Those skilled in the art will recognize that this application can be implemented as a system, method, or computer program product. Therefore, this disclosure can be implemented in the following forms: it can be entirely hardware, entirely software (including firmware, resident software, microcode, etc.), or a combination of hardware and software, generally referred to herein as a "module" or "system." Furthermore, in some embodiments, this application can also be implemented as a computer program product contained in one or more computer-readable media, which contains computer-readable program code. Computer-readable storage media can be, for example, but not limited to—electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatuses, or devices, or any combination thereof.

[0058] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0059] Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of this application.

Claims

1. A method for exploring mine roadways, characterized in that, include: The equivalent back flux transient electromagnetic method is used to explore the area to be explored in a mine and obtain the distribution of abnormal resistivity. This includes: using the equivalent back flux transient electromagnetic method, determining the locations of multiple measuring points in the area to be explored based on preset distance requirements. These measuring points are arranged in a grid pattern on the surface of the area to be explored, with a grid spacing of 20m × 30m, and obtaining the apparent resistivity at each measuring point; using a transient relaxation inversion method to invert the multiple apparent resistivityes to obtain multiple relative resistivityes; for each relative resistivity, obtaining the difference between the relative resistivity and the corresponding apparent resistivity; when the absolute value of the difference is greater than a preset abnormal value, the corresponding apparent resistivity is taken as the abnormal resistivity; and obtaining the distribution of the multiple abnormal resistivityes as the abnormal resistivity distribution of the area to be explored. The high-density resistivity method is used to determine the profile resistivity distribution based on the abnormal resistivity distribution, including: dividing the abnormal region based on the abnormal resistivity distribution; intersecting the abnormal region at a large angle to obtain a high-density resistivity survey line; using a Winner and Schlumberger device to jointly survey the vertical profile of the high-density resistivity survey line, first conducting an initial exploration with a measurement point spacing of 10m, and then conducting a detailed exploration with a measurement point spacing of 5m, obtaining multiple profile resistivity; and acquiring the distribution of the multiple profile resistivity as the profile resistivity distribution of the area to be explored. Determining the initial location of the roadway based on the abnormal resistivity distribution and the profile resistivity distribution includes: dividing the area into anomaly regions based on the abnormal resistivity distribution; dividing the area into high-density regions based on the profile resistivity distribution; obtaining the overlapping areas between the anomaly regions and the high-density regions, the overlapping areas including multiple overlapping positions; determining the sorting order of the multiple overlapping positions according to preset sorting requirements; and connecting the multiple overlapping positions according to the sorting order to obtain the initial location of the roadway. The initial tunnel location is verified to obtain the target tunnel location, including: using the micro-motion exploration method to measure the wave velocity of the initial tunnel location, identifying the tunnel location based on the difference between the wave velocity when the tunnel is filled with water or has a cavity and the surrounding geological body, thus obtaining the precise tunnel location; using the charging method, electrical measurements are performed at the tunnel outlet using a rated power higher than the preset power, and potential anomaly detection is performed on the ground surface within a range of 1km to 3km from the measurement point to obtain the target tunnel location; wherein, the measurement points of the micro-motion exploration method and the charging method are the same as the measurement points used in the detailed exploration using the high-density resistivity method, and the distance between the measurement points is 5m; Drilling and scanning are performed on the target tunnel location to obtain target survey information, which includes the development status of fractures in and around the tunnel, the flow path of underground acid water, and the bottom boundary of hydrogeological conditions. The preset sorting requirements include: selecting one of multiple overlapping positions as the starting position and marking the starting position with a number of 1; finding the first and second closest overlapping positions from the unmarked overlapping positions; marking the first overlapping position with a number equal to -1 of the starting position and the second overlapping position with a number equal to +1 of the starting position; finding the third overlapping position from the unmarked overlapping positions with a number equal to -1 of the first overlapping position and marking the third overlapping position with a number equal to -1 of the first overlapping position; finding the fourth overlapping position from the unmarked overlapping positions with a number equal to +1 of the second overlapping position and marking the fourth overlapping position with a number equal to +1 of the second overlapping position, until each overlapping position is marked with a number; and sorting the multiple overlapping positions according to their numbers from smallest to largest to obtain the corresponding sorting order.

2. A mine roadway exploration system, characterized in that, include: The first acquisition module is used to explore the area to be explored in the mine using the equivalent reverse flux transient electromagnetic method to obtain the abnormal resistivity distribution. The first acquisition module is specifically used for: using the equivalent reverse flux transient electromagnetic method, determining the locations of multiple measuring points in the area to be detected in the mine based on a preset distance requirement; the measuring points are arranged in a grid pattern on the surface of the area to be detected with a grid spacing of 20m×30m, and acquiring the apparent resistivity of each measuring point; using the transient relaxation inversion method to invert the multiple apparent resistivityes to obtain multiple relative resistivityes; and for each relative resistivity, acquiring the difference between the relative resistivity and the corresponding apparent resistivity. When the absolute value of the difference is greater than a preset abnormal value, the corresponding apparent resistivity is taken as the abnormal resistivity; the distribution of multiple abnormal resistivities is obtained as the abnormal resistivity distribution of the area to be detected. The second acquisition module is used to determine the profile resistivity distribution based on the abnormal resistivity distribution using the high-density resistivity method. Specifically, the second acquisition module is used to: delineate abnormal regions based on the abnormal resistivity distribution; intersect the abnormal regions at large angles to obtain high-density resistivity survey lines; conduct joint surveys of the vertical profiles of the high-density resistivity survey lines using a Winner and Schlumberger device, first performing an initial exploration with a measurement point spacing of 10m, then a detailed exploration with a measurement point spacing of 5m, obtaining multiple profile resistivity values; and acquire the distribution of the multiple profile resistivity values ​​as the profile resistivity distribution of the area to be detected. The determination module is used to determine the initial location of the roadway based on the abnormal resistivity distribution and the profile resistivity distribution. Specifically, the determination module is used to: divide the area into anomaly regions based on the abnormal resistivity distribution; divide the area into high-density regions based on the profile resistivity distribution; obtain overlapping areas between the anomaly regions and the high-density regions, the overlapping areas including multiple overlapping positions; determine the sorting order of the multiple overlapping positions according to preset sorting requirements; and connect the multiple overlapping positions according to the sorting order to obtain the initial location of the roadway. The verification module is used to verify the initial measured tunnel location to obtain the target tunnel location. Specifically, the verification module is used to: measure the wave velocity of the initial measured tunnel location using the micro-motion exploration method, identify the tunnel location based on the difference between the wave velocity and the surrounding geological body when the tunnel is filled with water or has a cavity, and obtain the precise measured tunnel location; and use the charging method to perform electrical measurements at the tunnel outlet using a rated power higher than the preset power, and detect potential anomalies on the ground surface within a range of 1km to 3km from the measurement point to obtain the target tunnel location. The measurement points for the micro-motion exploration method and the charging method are the same as the measurement points for detailed exploration using the high-density resistivity method, and the distance between the measurement points is 5m. The module is used to perform drilling and scanning on the target tunnel location to obtain target survey information; the target survey information includes the development status of fractures in and around the tunnel, the underground acid water flow path, and the bottom boundary of hydrogeological conditions; The preset sorting requirements include: selecting one of multiple overlapping positions as the starting position and marking the starting position with a number of 1; finding the first and second closest overlapping positions from the unmarked overlapping positions; marking the first overlapping position with a number equal to -1 of the starting position and the second overlapping position with a number equal to +1 of the starting position; finding the third overlapping position from the unmarked overlapping positions with a number equal to -1 of the first overlapping position and marking the third overlapping position with a number equal to -1 of the first overlapping position; finding the fourth overlapping position from the unmarked overlapping positions with a number equal to +1 of the second overlapping position and marking the fourth overlapping position with a number equal to +1 of the second overlapping position, until each overlapping position is marked with a number; and sorting the multiple overlapping positions according to their numbers from smallest to largest to obtain the corresponding sorting order.

3. A computing device, comprising a memory, a processor, and a program stored in the memory and running on the processor, characterized in that, When the processor executes the program, it implements the steps of the mine roadway exploration method as described in claim 1.

4. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores instructions that, when executed on a terminal device, cause the terminal device to perform the steps of the mine roadway exploration method as described in claim 1.

Images