Application of Comprehensive Geophysical Prospecting Method in Water Accumulation Exploration of Multilayer Goaf in Integrated Mine

A coal mine in Datong is an integrated mine. At present, there is goaf in the upper and lower part of the mining coal seam. There is a lot of ponding in the goaf, which has great potential safety hazards for production. In order to find out the scope and location of ponding in goaf, the comprehensive geophysical exploration method combining transient electromagnetic method and high-density resistivity method is used to carry out the research. Firstly, the time-base, turn-off time, receiving delay, current, superposition times, and other parameters of the instrument are tested on the surface of known goaf to obtain the best instrument parameters, and the parameters are used to verify the feasibility of the research scheme; then, the transient electromagnetic method is used for large-area exploration on the surface of the mine, the suspected goaf ponding area is found through comprehensive analysis, and the high-density resistivity exploration is arranged in the suspected goaf ponding area. According to the obtained results, the scope and location of the goaf ponding area are accurately located through comprehensive analysis. The results show that there are two goaf ponding areas in the exploration area, which are located above the 8# coal seam currently mined; the range and location of goaf ponding area can be accurately obtained by using the comprehensive geophysical method of high-density electrical method and transient electromagnetic method. This method can provide reference for mine water prevention and control in Datong area and has great practical significance to ensure coal mine safety production.


Introduction
In the production of integrated mine, the lack of geological data will cause great hidden trouble to the production safety [1]. At the same time, when the mines are some small coal mines merged and reorganized, most of the mining technology of these small coal mines is relatively backward, and the preserved geological data have low reference value [2,3]. erefore, the urgent task of coal mine production and construction is to solve the water damage problem in the goaf of old coal mines [4][5][6]. On November 29, 2020, a water inrush accident occurred in Yuanjiangshan Coal Mine in Hunan Province, China, due to the ultradeep and cross-boundary mining of the steep coal seam, under the combined action of rock pressure and water pressure in the upper gob area. Similar incidents show that when hydrogeological conditions are not accurately grasped, blind mining of coal mines has great hidden dangers [7,8]. Zuoyun Caoduogou Coal Industry Co., Ltd., of Shanxi Coal Import and Export Group is an integrated mine with less historical geological data, so the investigation of mine flooding is the primary task. At present, there are many methods to choose for the specific problems of mine water damage. Different experts and scholars have conducted in-depth research from the perspective of different exploration technologies and methods. Zhang et al. [9,10] have detected suspicious and key areas through geophysical exploration methods to more accurately obtain the goaf and hydrogeological conditions. Wen et al. [11] have adopted high-density resistivity and transient electromagnetic method and other geophysical methods have carried out theoretical research and engineering practice on the problem of water damage in goaf and found out a set of effective means for the exploration of shallow hydrogeological conditions; Yi [12,13] conducted transient electromagnetic interference test, optimization test, and electrical response characteristic analysis for shallow buried deep goaf, which improved the accuracy of shallow hydrogeological transient electromagnetic method exploration. From the above research results, many experts and scholars have conducted a lot of research around geophysical exploration, but the research results focus on solving a specific problem. In this paper, by using the comprehensive geophysical exploration method, the transient electromagnetic method, and high-density resistivity method, we can find out the situation of mine goaf ponding, so as to make a reference for coal mining planning and safety production.

Transient Electromagnetic Method.
e principle of the transient electromagnetic method is as follows: supply power to an earth-free coil, this coil can produce a pulsed magnetic field, the medium in the detecting direction after excited will produce eddy current [14], when the pulse intervals, the eddy current will not disappear immediately, and as a result of the existence of eddy current, a new magnetic field will be generated in the medium, and the signal of the new magnetic field can be received by another earth-free coils or grounded electrodes. [15][16][17]. Due to the similarities and differences in the structure, moisture content, scale, resistivity, and other factors of the underground medium, the attenuation characteristics of the secondary induced electromagnetic field are different, so the properties and spatial distribution patterns of the underground medium can be judged [18,19]. e technical principles are shown in Figure 1.
Transient electromagnetic method is the most sensitive method in improving the detection depth and finding lowresistance geological bodies in high-resistance areas. It has the advantages of automatic elimination of main noise sources, no topographic influence, combined observation at the same point, optimal coupling with the detection target, strong anomaly response, simple shape, and strong resolution. It is often used in the investigation of water bearing geology, such as karst caves and channels, coal mine goaf, and deep irregular water bearing structure [20,21].

High-Density Resistivity Method.
e working principle of electrical prospecting is that the resistivity of each stratum is different due to the different lithology; when the current is transmitted to the underground, the change of the underground geological condition is determined by observing the change of the resistivity received by the measuring electrode [22,23]. By understanding the characteristics reflected by the electrical measurement curve, the structural state, electrical property, and scale of the geological structure can be analyzed [24]. High-density resistivity method is a combination of electrical profile method and electrical sounding method; the geological interpretation map of the measured geoelectric section can be given quickly and accurately through the processing of various parameters obtained by corresponding programs and automatic inversion imaging. It is often used in many engineering exploration fields, such as engineering geological survey of major sites, site selection of dam foundation and pier, gob, and ground fissure detection [25]. e advantage of high-density resistivity method is that by arranging a single arrangement the position of the measuring point can be changed by instrument control, and then by increasing the distance of the power supply electrode and the measuring electrode on the same measuring point, the detection distance can be deepened. erefore, the change law of apparent resistivity from near to far in the detection direction of the same measuring point can be mastered [26]. e technical principles are shown in Figure 2 e exploration area is a long strip area of 380 m × 1040 m, with an area of 0.3952 km 2 . As shown in Figure 3, the exploration area belongs to the gentle hilly region of the Loess plateau, and the terrain in the area is generally high in the west and low in the east.  Advances in Civil Engineering

Geophysical Characteristics.
e working conditions of the ground in the exploration area are complex. Although the relative height difference is only 68.2 m, there is a loess gully and six high voltage wires with irregular distribution with the measuring line. e geophysical conditions of the shallow surface are complex. e electrical parameters of each coal seam and surrounding rock are shown in Table 1, it can be seen that the apparent resistivity of rock formations with the same lithology is positively correlated with the rock particle size, the apparent resistivity of rock formations with different lithology is positively correlated with their density, the electrical properties of coal seam and rock stratum are significantly different, the apparent resistivity difference is as high as 2-3 times, and the coal seam is the highest.

Field Test and Engineering Application
e instruments used for detection are IGGETEM-30B transient electromagnetic instrument and WDJD-2 highdensity resistivity measurement system.

Instrument Parameter Test.
According to the surface geological conditions of the exploration area, the single point experiment is selected at the 40th measuring point of the third line in the exploration area (it is located above the known goaf of 8113 working face of the coal mine). e control variable method is used to compare the attenuation curves under different parameters and then determine the appropriate instrument parameters.
(1) Time-Base Test. Under the premise that the fixed transmitting wire frame was 15 m × 30 m, turn-off time was 75 µs, current was 7.5 A, superposition times was 512, and receiving delay was 250 µs and Log14, the time-base tests of 20 ms and 40 ms are carried out, respectively, and the results are shown in Figure 4. It can be seen that the V (t)/I attenuation curves of 20 ms and 40 ms are basically the same, but because there is some interference in the exploration area, and the detection of coal layers is deeper, the time-base selection is 40 ms.
(2) Turn-Off Time Test. Under the premise that fixed transmitting wire frame was 15 m × 30 m, time-base was 40 ms, current was 7.5 A, superposition times was 512 times, and receiving delay was 250 µs and Log14, three tests with turn-off time of 50 µs, 75 µs, and 90 µs are carried out, and the results are shown in Figure 5. According to the analysis in Figure 5(b), the V (t)/I attenuation curves of 75 µs and 90 µs in the late period were well responded, while the V (t)/I attenuation curves of 50us in the early period were significantly higher. e V (t)/I attenuation curves of 75 µs and 90 µs were basically the same; therefore, the turn-off time was determined to be 75 µs.
(3) Receiving Delay Test. Under the premise that fixed transmitting wire frame was 15 m × 30 m, turn-off time was 75 µs, time-base was 40 ms, current was 7.5 A, and Advances in Civil Engineering superposition times was 512 times and Log14, three tests with receiving delay time of 200 µs, 250 µs, and 300 µs are carried out, and the results are shown in Figure 6. It can be seen that the overall linearity of V (t)/I attenuation curve is good in the late period, and the V (t)/I attenuation curve of 200us in the early period appears as a straight section, mainly due to the composition of the primary field. e V (t)/I attenuation curves of 250 µs and 300 µs are basically the same, and the V (t)/I attenuation curve of 250 µs is better.
(4) Current Test. Under the premise that fixed transmitting wire frame was 15 m × 30 m, turn-off time was 75us, time-base was 40 ms, receiving delay was 250 µs, and superposition times was 512 times and Log14, two tests of transmitting current 3.7 A and 7.5 A are carried out, and the results are shown in Figure 7. It can be seen that because of the existence of interference 3.7 A current attenuation curve has obvious inflection point, and the V (t)/I attenuation curve of 7.5 A is better than 3.7 A. e test results are shown in Figure 8. Generally speaking, if the attenuation curve is to be smoother, it requires a high enough number of superpositions. Considering the working time, the superposition times were determined to be 512 times.
It can be found that the value of the collected voltage is closely related to the turn-off time and delay time but has little relationship with the power supply current. Finally, the instrument parameter coil adopts a 15 × 30 m '∞' shape loop, the emission current is 7.6 A, the time-base is 40 ms, the superposition number is 512, the turn-off time is 75 µs, the take-off time is Log14, and the receiving delay time is 250 µs.

Method Effectiveness Test.
Select the known goaf water area, set the test line as 300 m, instrument parameters using the parameters selected by the test, the wire frame using 15 × 30 m '∞' type overlapping loop, and the test results are shown in Figures 9 and 10. It can be seen that there is an obvious high value anomaly at 23-37 points (red frame), and the anomaly amplitude is large, which is consistent with the known location of goaf water accumulation area. In Figure 10, the contour lines at 22-36 points (red frame) show obvious low-resistance anomalies, and the water accumulation area is obviously reflected, indicating that the effect of using the test device and instrument parameters in this exploration area is good, and the water accumulation in goaf is obvious. At the same time, it can be seen that the inversion depth of the transient electromagnetic method combined with known data can reach more than 200 m, and the goaf layer of the coal seam can be further determined from the depth.

Engineering Applications.
A total of 20 measuring lines are designed in the exploration area; each line is 1040 m long and 20 m apart. 53 measuring points are designed on each measuring line, and the spacing is 20 m. e south side of the exploration area is line 1, arranged in turn north. e west side of each line is No. 1 measuring point, which is arranged in turn east. e layout of the line is shown in Figure 11.
Firstly, the transient electromagnetic method is used to carry out full coverage exploration in the exploration area to determine the suspected water accumulation area, and then the high-density resistivity method is carried out at the position of the suspected water accumulation area. e highdensity resistivity method uses α arrangement as shown in Figure 12. e distance between the observation points is 10 m, the number of measurement channels is 60, and the distance between the channels is 10 m. e cross section of the detection results obtained by this method is a flat 'U' type. e selection of time parameters is as follows: power supply time is 2000 ms, power outage delay is 200 ms, power supply voltage is about 500 V, and power supply current is about 400 mA.

Analysis of Typical Measuring
Lines. Several measuring lines with obvious electrical changes, such as 5, 6, 9, 13, and 20, were selected to analyze their multichannel voltage profile and apparent resistivity profile (see Figures 13 and 14). It can be seen that there are high-voltage and low-resistance anomalies (red ellipses) at 1-5 points of 5 measuring lines, 28-31 points of 6 measuring lines, 23-30 points of 9 measuring lines, 25-28 points of 13 measuring lines, and 19 and 47 points of 20 measuring lines. e low-resistance anomaly area is basically about 1170-1180 m, which is speculated to be the reflection of goaf water in 8 # coal seam. Among them, 15-25 points of 13 measuring lines (blue frame) are the industrial site of air pit, and there are lost data; therefore, the high-voltage and low-resistance anomaly area are no longer analyzed.

Analysis of Typical Plane Graph.
According to the elevation and thickness of the coal seam in the exploration area, the plan of equal voltage value of the 24th and 28th measuring channels (see Figure 15) is selected to analyze the voltage change. Combined with the analysis of the drawings of the typical measuring lines, it can be seen that there are mainly two abnormal areas of high voltage (more

Advances in Civil Engineering
By drawing the plan of apparent resistivity of 8 # coal seam (see Figure 16), it can be concluded that the low-resistance region is basically corresponding to the high voltage region of the plan of equal voltage value of 24th and 28th measuring channels. e difference is that the low-resistance region in the northwest corner is increased, and the analysis is mainly caused by the increase of the terrain in the northwest corner of the exploration area.

Data Analysis of High-Density Resistivity Method.
According to the suspected goaf water accumulation areas, the high-density resistivity method is arranged at 16-45 points of 6 measuring lines and 25 points of 5-25 measuring lines, and the inversion simulation apparent resistivity section with the topography of each measuring line is plotted (see Figure 17). It can be seen that there are large low resistivity anomalies at the contour line of    Combined with geological data and surface conditions, it is comprehensively inferred that these two areas are the goaf water area of 8 # coal seam, and the location of the water area is shown in Figure 18. At the same time, there is an industrial site of air pit in the exploration area, which cannot be measured; there are many interference factors such as ferromagnetic materials and high voltage lines around it, which may bring    10 Advances in Civil Engineering multiple solutions to the data processing and have a great influence on the accuracy of the results.

Conclusions
(1) Firstly, the transient electromagnetic method is used to carry out full coverage exploration in the exploration area, and then the high-density resistivity method is carried out for the suspected goaf water area, which can ensure the reliability of geophysical prospecting results. (2) By using the comprehensive geophysical prospecting method, two goaf water accumulation areas with large scope were measured in the mine with multilayer goaf. According to their positions and elevations, it was deduced that the two goaf water accumulation areas were both located in the 8 # coal seam, which was located in the upper part of the coal seam under mining, providing a reliable reference for coal mine safety construction. (3) When there are buildings such as the industrial site of air pit or high voltage lines in the exploration area, there is a certain range of electromagnetic field around them, which has a great influence on the above two geophysical methods and can be combined with other geophysical methods such as seismic exploration to make the exploration results more accurate.

Data Availability
e data used to support the findings of this research are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest.