# Application of Comprehensive Geophysical Prospecting Method in the Exploration of Coal Mined-Out Areas

Authors: Xiao Luo; Shuang Gong; Zhonggang Huo; Hongyan Li; Xinpin Ding
Journal: Advances in Civil Engineering
 (2019)

---

## Abstract

Due to the influence of mining technology, some coal mines will leave mined-out areas (or goafs) in the rock mass after mining. If the distribution law and influence scope of mined-out areas cannot be accurately explored and controlled, it will bring great potential safety hazards to subsequent coal mining and cause serious harm to adjacent houses, roads, bridges, railways, and so on. The unmanned aerial vehicle (UAV) aerial survey is used to accurately acquire the three-dimensional model of the slope in the key area of Yuanbaoshan open-pit coal mine, which provides basic information for the real-time visualization of the open-pit mine and the establishment of the engineering geological research profile. And the distribution of mined-out areas and water accumulation in the coal mine are investigated in detail by using the comprehensive geophysical prospecting method of the transient electromagnetic method (TEM) and controlled source audio-frequency magnetotellurics (CSAMT). Then, the borehole detection method is used to verify the goaf range obtained from comprehensive geophysical prospecting, and the results show good consistency. In addition, based on the results of borehole exploration and the comprehensive analysis of the borehole image, the instantaneous water level elevation in goaf and the development height of the water-conducting fracture zone and caving zone are concluded. Finally, the distribution of coal mined-out areas in the exploration area is obtained, and the results show that through the use of the transient electromagnetic method and controlled source audio-frequency magnetotellurics method, combined with drilling detection methods, the purpose of accurate exploration of mined-out areas in coal mines can be achieved.


---

## Body

## 1. Introduction

In recent years, the large-scale exploitation and utilization of coal resources in China has brought huge economic and social benefits. At the same time, the overexploitation of coal has formed a large number of mined-out areas (or goafs) underground. The coal mined-out area has led to large-scale surface subsidence. Ground buildings, transportation facilities, and living infrastructure and land have been seriously damaged [1–5]. It has brought great threat to people’s lives and production in the mining area and has caused serious impact and harm to China’s economic construction and social stability. Therefore, finding out the specific situation of the mined-out area and avoiding the occurrence of disastrous geological phenomena such as collapse, water inrush, and roof caving have become urgent problems to be solved. Furthermore, the open-pit mining method in coal mining, with its advantages of full utilization of resources, high production efficiency, low cost, and good working conditions, is taking up an increasing proportion in China’s coal mining [6, 7]. However, for a long time, many small coal mines in China have adopted room-and-pillar mining methods, leaving many goafs when they are closed. The size, location, and depth of these holes are uncertain, and with the progress of open-pit mining above the concealed goaf and the passage of time, the holes gradually collapse. It is difficult to grasp the actual void stratum. Because the roof collapse of underground goaf is concealed, there is almost no warning in advance. Once the collapse occurs, it may cause serious safety accidents. The operation of super-large equipment in the open-pit mine has great potential safety hazards in the upper part of goaf. Once the equipment falls into the goaf, it may cause casualties and equipment damage.Quantitative evaluation of goaf stability requires accurate geological data, and determining the exact location of goaf is the primary task of engineering geological investigation. The conventional survey method is to detect the core by arranging more boreholes [8, 9]. This method is essential to determine the location of the goaf. However, the implementation of drilling often costs a lot of time and money, especially for some projects under construction. Even if there is a guarantee of drilling funds, but the exploration time is limited, it may not be allowed to implement more boreholes or the exploration site also has limitations on drilling, which requires the study of rapid and accurate exploration technology of goafs or mined-out areas. Compared with the borehole detection method, the geophysical method is widely used in detection because of its high efficiency and nondestructive characteristics [10–14]. At present, there are many geophysical techniques. Each detection method has certain effect but also has limitations. For example, ground penetrating radar has strong resolution but limited depth [15, 16]. Seismic and Rayleigh wave-layering ability is strong, but the response to water filling is not obvious [17]. The direct current method is greatly affected by grounding conditions and topographic fluctuations [18, 19]. The transient electromagnetic method is less affected by topography and has better detection effect on water accumulation in goaf of the coal mine, but it has no obvious response to high-resistance anomaly [20–23]. It can be seen that the previous comprehensive geophysical methods are mostly applied to hydrogeological conditions, landslides, and groundwater exploration. However, there are few studies on the exploration of goaf and water-bearing areas in coal mines. And the feasibility and validity of the comprehensive geological exploration method applied to the underground goaf are rarely reported. According to the physical characteristics of different mined-out areas, the main technical way to solve these problems is to adopt the comprehensive geophysical prospecting method with necessary drilling verification. However, its actual exploration effect needs further research and verification.In the present study, the unmanned aerial vehicle (UAV) aerial survey is used to accurately acquire the three-dimensional model of the slope in the key area of Yuanbaoshan open-pit coal mine, which provides basic information for the real-time visualization of the open-pit mine and the establishment of the engineering geological research profile. And the distribution of mined-out areas and water accumulation in the coal mine are investigated in detail by using the comprehensive geophysical prospecting method of transient electromagnetic method (TEM) and controlled source audio-frequency magnetotellurics (CSAMT). Then, the borehole detection method is used to verify the goaf range obtained from comprehensive geophysical prospecting. In addition, based on the results of borehole exploration and the comprehensive analysis of the borehole image, the instantaneous water level elevation in goaf and the development height of water-conducting fracture zone and caving zone are concluded. Finally, the comprehensive results and the predicted goaf range are analyzed and discussed.


## 2. General Situation of Exploration Area

### 2.1. General Geological Situation

Yuanbaoshan open-pit coal mine is located 35 km northeast of Chifeng City, Inner Mongolia Autonomous Region, China, as shown in Figure1. The main coal seams in the mine are large in thickness, shallowly buried, and there are many production mines and small coal pits nearby. There are five closed small coal mines in the first mining area. Because all these small coal mines are roof-and-pillar mining, a large number of underground void areas or mined-out areas are left behind, which poses a serious threat to the safety of production on the spot. In recent two years, accidents of shovels and trucks falling into the goaf have occurred successively, as shown in Figures 2(a) and (2b). Therefore, in order to avoid accidents and realize safe mining, it is necessary to conduct the comprehensive exploration of the goafs (mined-out areas) distributed within the boundary of the open-pit coal mine. It also provides basic data and information for the stability evaluation and effective control of the proven goaf and water accumulation areas.Figure 1
Location of Yuanbaoshan open-pit coal mine.Figure 2
Safety accident scene in the open-pit mine due to underground goaf. (a) Due to the insufficient safety thickness of the roof of unknown goaf, heavy truck fell into the underground goaf 30 m deep. (b) During the mining process of the open-pit coal mine, the track on the right side of shovel fell into the collapsed underground goaf.
(a)
(b)According to the geological data collected in the past and the urgency of solving problems, the research area of this paper is divided into two parts. As shown in Figure3, the area of exploration area A is 0.26 km2, the area of exploration area B is 0.11 km2, and the total area is about 0.37 km2.Figure 3
Zoning map of the comprehensive geophysical prospecting study area.pi In recent years, unmanned aerial vehicle (UAV) surveying and mapping technology has gradually become the main tool of surveying and mapping due to its advantages of automation, intelligence, low cost, and high precision. The UAV aerial survey is used to accurately acquire the three-dimensional model of the slope in the key area of Yuanbaoshan open-pit coal mine, which provides basic information for the real-time visualization of the open-pit mine and the establishment of the engineering geological research profile. Some aerial survey results of UAV are shown in Figure4. The aerial survey results show that the surface elevation of the open-pit mining area is +472 m to +490 m. The topographic relief patterns of the two exploration areas are shown in Figure 5.Figure 4
Aerial survey achievements of unmanned aerial vehicles. (a) Top view of the working area of the bucket-wheel excavator. (b) Top view of the eastern area of the working slope. (c) Three-dimensional model of the surface local area in the middle part of the working slope. (d) Three-dimensional surface model of the eastern area of the working slope.
(a)
(b)
(c)
(d)Figure 5
Surface topographic map of exploration area. (a) Topography of exploration area A. (b) Topography of exploration area B.
(a)
(b)


### 2.2. Geophysical Characteristics

In order to avoid the multisolution of the single geophysical exploration method, improve the accuracy of exploration, and have a qualitative understanding of the geophysical characteristics of the exploration area, the resistivity logging method [24, 25] is used to study the resistivity distribution law of different depth strata in the exploration area. The column diagram of the resistivity logging method in exploration area B is shown in Figure 6. It can be seen that shallow strata over −70 m present high resistivity characteristics. The strata near −70 m to −180 m in the middle of the column diagram are characterized by low resistivity. About −180 m to −380 m is the reflection of coal-bearing strata, showing high resistivity. The deep strata below −380 m are interbedded with medium and fine sandstones, showing low resistivity as a whole. In summary, the strata in the exploration area show the electrical characteristics of “high resistivity to low resistivity to high resistivity to low resistivity” from shallow to deep vertically.Figure 6
Column diagram of borehole resistivity in exploration area B.


## 2.1. General Geological Situation

Yuanbaoshan open-pit coal mine is located 35 km northeast of Chifeng City, Inner Mongolia Autonomous Region, China, as shown in Figure1. The main coal seams in the mine are large in thickness, shallowly buried, and there are many production mines and small coal pits nearby. There are five closed small coal mines in the first mining area. Because all these small coal mines are roof-and-pillar mining, a large number of underground void areas or mined-out areas are left behind, which poses a serious threat to the safety of production on the spot. In recent two years, accidents of shovels and trucks falling into the goaf have occurred successively, as shown in Figures 2(a) and (2b). Therefore, in order to avoid accidents and realize safe mining, it is necessary to conduct the comprehensive exploration of the goafs (mined-out areas) distributed within the boundary of the open-pit coal mine. It also provides basic data and information for the stability evaluation and effective control of the proven goaf and water accumulation areas.Figure 1
Location of Yuanbaoshan open-pit coal mine.Figure 2
Safety accident scene in the open-pit mine due to underground goaf. (a) Due to the insufficient safety thickness of the roof of unknown goaf, heavy truck fell into the underground goaf 30 m deep. (b) During the mining process of the open-pit coal mine, the track on the right side of shovel fell into the collapsed underground goaf.
(a)
(b)According to the geological data collected in the past and the urgency of solving problems, the research area of this paper is divided into two parts. As shown in Figure3, the area of exploration area A is 0.26 km2, the area of exploration area B is 0.11 km2, and the total area is about 0.37 km2.Figure 3
Zoning map of the comprehensive geophysical prospecting study area.pi In recent years, unmanned aerial vehicle (UAV) surveying and mapping technology has gradually become the main tool of surveying and mapping due to its advantages of automation, intelligence, low cost, and high precision. The UAV aerial survey is used to accurately acquire the three-dimensional model of the slope in the key area of Yuanbaoshan open-pit coal mine, which provides basic information for the real-time visualization of the open-pit mine and the establishment of the engineering geological research profile. Some aerial survey results of UAV are shown in Figure4. The aerial survey results show that the surface elevation of the open-pit mining area is +472 m to +490 m. The topographic relief patterns of the two exploration areas are shown in Figure 5.Figure 4
Aerial survey achievements of unmanned aerial vehicles. (a) Top view of the working area of the bucket-wheel excavator. (b) Top view of the eastern area of the working slope. (c) Three-dimensional model of the surface local area in the middle part of the working slope. (d) Three-dimensional surface model of the eastern area of the working slope.
(a)
(b)
(c)
(d)Figure 5
Surface topographic map of exploration area. (a) Topography of exploration area A. (b) Topography of exploration area B.
(a)
(b)


## 2.2. Geophysical Characteristics

In order to avoid the multisolution of the single geophysical exploration method, improve the accuracy of exploration, and have a qualitative understanding of the geophysical characteristics of the exploration area, the resistivity logging method [24, 25] is used to study the resistivity distribution law of different depth strata in the exploration area. The column diagram of the resistivity logging method in exploration area B is shown in Figure 6. It can be seen that shallow strata over −70 m present high resistivity characteristics. The strata near −70 m to −180 m in the middle of the column diagram are characterized by low resistivity. About −180 m to −380 m is the reflection of coal-bearing strata, showing high resistivity. The deep strata below −380 m are interbedded with medium and fine sandstones, showing low resistivity as a whole. In summary, the strata in the exploration area show the electrical characteristics of “high resistivity to low resistivity to high resistivity to low resistivity” from shallow to deep vertically.Figure 6
Column diagram of borehole resistivity in exploration area B.


## 3. Comprehensive Geophysical Prospecting Method

### 3.1. Transient Electromagnetic Method

The transient electromagnetic method (TEM) belongs to the time-domain electromagnetic induction method. Before 1990s, it was mainly used to search for metal deposits, coalfields, petroleum, gas, and other fields. Its basic principle is to send a pulse magnetic field (or electric field) to the underground by using an ungrounded loop (or a couple source), that is, to send a current pulse square wave on the transmitting loop, and the square wave will dropinstantaneously. Generate a magnetic field that can propagate underground. Under the excitation of the primary magnetic field, the geological body will produce eddy current, whose size depends on the conductivity of the geological body. After primary magnetic field disappears, the eddy current can not disappear immediately. It will have an excessive process (attenuation process), which produces attenuated secondary magnetic field to propagate to the surface. The secondary magnetic field will be received by the receiving line or magnetic probe of the surface instrument. The variation of the secondary magnetic field will reflect the electrical distribution of underground geological bodies. Secondary induction electromotive force (SIEMF) is measured with different delay times in the receiver, so the attenuation characteristics of the secondary magnetic field with time are obtained. The pseudosection map of resistivity inversion can be obtained by calculating program, and then the geological condition of deep underground can be judged. The advantages of TEM are easy handling, strong adaptability to surface sites, and high efficiency [20].


### 3.2. Controlled Source Audio-Frequency Magnetotellurics

For a long time, people can use various exploration methods to solve the engineering geology and hydrogeology problems concerned about the shallow surface geology of the engineering area, but there has been a lack of effective methods to understand the deep geological conditions. In geophysical exploration, electrical prospecting has been paid attention to because of its variability of field sources and effectiveness of methods. Controlled source audio-frequency magnetotellurics (CSAMT) is a new electrical method developed in the late 1980s [26]. Based on electromagnetic wave propagation theory and Maxwell equations, Cagniard’s apparent resistivity [27, 28] equation is derived as follows:(1)ρaC=1μωExHy2,where E and H are the orthogonal electric and magnetic fields, respectively. ω and μ denote angular frequency and magnetic permeability. If the electric field is measured in millivolt/km and magnetic field is measured in nT (nanotesla), the above equation becomes(2)ρaC=15fExHy2.According to the skin effect theory of the electromagnetic wave, the formula of skin depth can be derived as follows:(3)δ=2ωμσ=503ρf,where ρ is the electrical resistivity and f is the frequency.It has been shown that when the surface resistivity is fixed, the propagation depth (or detection depth) of the electromagnetic wave is inversely proportional to the frequency. In high frequency, the detection depth is shallow; in low frequency, the detection depth is deep. Therefore, the detection depth can be changed by adjusting the transmission frequency to achieve the purpose of detection.The general principle of the CSAMT method is to use controllable artificial field source and measure the Cagniard’s apparent resistivity by the ratio of electric field to magnetic field. At the same time, based on the principle of skin depth of the electromagnetic wave, electrical sounding of different depths is carried out by changing frequency. Therefore, the CSAMT method has the following advantages in its work: it is suitable for various sites, high efficiency, and can receive seven electric fieldsE and one magnetic field H at the same time in a single launch; it has strong anti-interference ability; it has a large exploration depth range, up to 1,000 m; it can detect the resistivity discrepancies in the geological background, has small high-resistance shielding effect, and can penetrate high-resistance strata. It is precisely because of these advantages that it was rapidly popularized and applied in engineering survey all over the world in the late 1990s.V8 multifunctional electromagnetic prospecting system is used, and the main technical parameters of the instrument are as follows: input impedance >100 MΩ, compensation error <0.12‰, resolution: 1μV, power: 2800 W, maximum output current: 40 A, and maximum supply voltage: 140 V.


### 3.3. Data Processing Method

The data collected on site must be further processed before they can be transformed into the required geophysical and geological results. The main parameters of TEM field observation are sampling time, emission current, and normalized induction secondary field, etc. Data processing begins with preprocessing, that is, data cataloguing, quality inspection, rejection of unqualified data, data quality statistics and classification, and sorting out the sequence and format required for special software. Secondly, the data are filtered to remove the suppressed interference signals and restore the original characteristics of geological signals. Then, the parameters are converted; that is, the apparent resistivity and apparent depth are obtained by software conversion. On this basis, topographic correction and other processing are done according to survey, geology, and drilling data. Finally, various sectional diagrams are drawn for data analysis.There are large bucket excavator system and material conveying tape in exploration area A. Its working voltage is about 6000 v, which makes the data near the operation area subject to certain electromagnetic interference. In addition, other large-scale operation vehicles in the exploration area will also have a certain impact on data acquisition. The data affected by electromagnetic interference will have abrupt maximum or minimum values, and even the phenomenon of local data uplifts or decreases. Figure7 is the contrast chart of the single-point quadratic field attenuation curve in double logarithmic coordinates. The attenuation curve without electromagnetic interference tends to decay in turn, and the data points basically do not show sudden change. The first part of the attenuation curve of the quadratic field of the measured point under electromagnetic interference basically maintains the trend of attenuation in turn. However, after 5 ms, the data show obvious “jitter and jump” phenomenon and the original attenuation law is broken.Figure 7
Comparison of single-point attenuation curves of transient electromagnetic interference.Correction is needed for the disturbed data. Hundreds of sampling points were extracted, and the interfered data were processed by combining the method of limit difference, fitting, and multipoint smoothness. The comparison of resistivity sections before and after processing is shown in Figure8, in which Figure 8(a) is the apparent resistivity section of data subject to electromagnetic interference and Figure 8(b) is the apparent resistivity section after interference correction. After correction, the change rule of affected data is consistent with that of data not subject to electromagnetic interference.Figure 8
Comparison of pseudosection of resistivity (a) before and (b) after interference correction.
(a)
(b)The measured TEM data are the comprehensive reflection of the underground electrical layers. Vertically, the original curve or the original apparent resistivity section can only show the comprehensive electrical change trend of several sets of strata. To interpret the distribution of the real strata or the closest to the real strata and geological anomalies through the measured data, the data should be inverted to obtain the subsurface resistivity distribution.Figure9 shows the different data processing methods for transient electromagnetic exploration. From Section 2.2, it can be seen that the strata in the exploration area present the electrical law of “high-low-high-low” from shallow to deep. By comparing (a) and (b) in Figure 9, it can be concluded that the section formed by inversion results reflects the electrical property law of the formation from shallow to deep “low-high-low-high,” while the section of the late formula shows the change characteristic of “high-low” from shallow to deep, which is not consistent with the actual electrical property of the formation. In addition, although the late formula has the advantages of simplicity and rapidity, the algorithm is suitable for the late time. There is some distortion in the early calculation, which results in unreal high resistivity response in the shallow section, and the inversion calculation restores the actual electrical property of the formation. Therefore, the inversion method is chosen as the data processing method of TEM exploration.Figure 9
Different data processing methods for transient electromagnetic exploration. (a) Late formula for calculating the resistivity section. (b) Inversion method for calculating the resistivity section.
(a)
(b)The data processing flow of controlled source audio-frequency magnetotellurics (CSAMT) is shown in Figure10. Detailed data processing can be referred to [26–28].Figure 10
Data processing flow chart of the CSAMT method.


### 3.4. Results of Integrated Geophysical Prospecting

By using the integrated geophysical prospecting method (TEM and CSAMT), the mined-out area and water accumulation area of 6 coal seams in the exploration area were comprehensively explored. If the stratum is not affected by the goaf or the water-conducting structure, the resistivity of coal-bearing stratum changes orderly, and the isoline changes steadily in the resistivity section, showing a layered distribution. When there is a goaf water accumulation area or water-conducting structure, the resistivity of anomalous areas decreases, and the isoline distribution is distorted, deformed, or densely striped. Typical results for exploration areas A and B are shown in Figures11–15, respectively. The location of the selected exploration line is shown in Figure 3.Figure 11
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 5 of the TEM method. (b) Inversion of the resistivity profile from line 5 of the CSAMT method.
(a)
(b)Figure 12
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 9 of the TEM method. (b) Inversion of the resistivity profile from line 9 of the CSAMT method.
(a)
(b)Figure 13
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 13 of the TEM method. (b) Inversion of the resistivity profile from line 13 of the CSAMT method.
(a)
(b)Figure 14
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 31 of the TEM method. (b) Inversion of resistivity profile from line 31 of CSAMT method.
(a)
(b)Figure 15
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 33 of the TEM method. (b) Inversion of the resistivity profile from line 33 of the CSAMT method.
(a)
(b)Figure11 shows the resistivity sections of line 5 obtained by different methods. The survey line is located in the southern boundary of the exploration area. The length of the survey line is 1180 m, and the surface roughly flat. Coal seams gradually deepen from west to east. It can be seen that the vertical resistivity of strata in the range of exploration depth generally shows the trend of “high resistivity-low resistivity-high resistivity-low resistivity,” which coincides with the electrical characteristics of the actual stratum described in Section 2.2. The horizontal resistivity isoline of strata near the coal seam is obviously distorted between No. 7 and No. 33 survey point, showing obvious characteristics of low-resistance anomaly. In addition, there is a normal fault F4 near No. 27 survey point. The fault distance is 30 m, which is not suitable for mining. Therefore, the anomaly can be inferred as a fractured water-rich zone. The horizontal resistivity isoline of strata near the coal seam is obviously distorted between No. 39 and No. 43 survey point, also showing obvious characteristics of low-resistance anomaly. As shown in Figure 11(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method. It is also shown that the electrical characteristics and the resistivity anomalies delineated by both methods are in agreement.Figure12 shows the resistivity sections of line 9 obtained by different methods. The survey line is located in the middle of exploration area A. The length of the survey line is 840 m, and the topography is relatively flat. The coal seam is gradually deepening from west to east from No. 25 to No. 47, and the whole coal seam is slightly uplifted from No. 48 to No. 66. From Figure 12(a), it can be seen that the vertical resistivity of the stratum in the exploration depth generally shows a trend of “high resistivity to low resistivity to high resistivity to low resistivity,” which coincides with the electrical characteristics of the actual stratum described in Section 2.2. The cross-sectional map is bounded by No. 50 survey point, and the rise and fall of the target electrical layer are obvious. From the trend of the coal seam, the trend of the left isoline is basically consistent with the trend of the actual stratum. Laterally, between No. 25 and No. 43 survey point, the resistivity isoline of the seam near the No. 6 coal seam is obviously distorted, showing obvious characteristics of low resistivity anomaly (the magenta ellipse in the picture). It is inferred that the low-resistance anomaly may be related to water accumulation in goaf formed by crossborder mining of 6 coals in adjacent coal mines. The isoline on the right side of the section (i.e., between the 50–63 survey points) changes obviously, the vertical high resistivity layer rises obviously upward, and the lower part of the section shows obvious low-resistance anomaly. It is inferred that it is the electrical reflection of No. 6 coal mined out and accumulated water. The upper high-resistance anomaly is similar to the electrical reflection of the nonfilling cracks formed by the caving of overlying strata after coal mining. As shown in Figure 12(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method.Figure13 shows the resistivity sections of line 13 obtained by different methods. The survey line is located in the north of exploration area A. The length of the survey line is 720 m. The surface is nearly flat, and the undulation is small. From Figure 13(a), it can be seen that the vertical resistivity of the strata in the exploration depth generally shows a trend of “high resistivity to low resistivity to high resistivity to low resistivity,” which is consistent with the electrical characteristics of the actual strata mentioned above. The resistivity isoline of the stratum near the coal seam is obviously distorted and shows obvious low-resistance anomaly between No. 36 and No. 46 as well as between No. 54 and No. 65 survey point transversely. It is inferred that the low-resistance anomaly may be related to water accumulation in goaf formed by crossborder mining of 6 coals in adjacent coal mines. As shown in Figure 13(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method, and the anomalies resistivity variations obtained by both methods are in agreement.The analysis of comprehensive results of some survey lines in exploration area A is provided earlier, and integrated geophysical exploration techniques are also carried out in exploration area B. Typical results are shown in Figures14 and 15. The low resistivity anomalous area near the coal seam can be obtained by comparing the modelling results of the two methods. It can be seen that the electrical characteristics and abnormal reactions obtained by the two geophysical methods are basically consistent. The results of deduced water-filled goaf and water-rich fracture area in exploration areas A and B are shown in Figure 3.


### 3.5. Verification of Borehole Television Method

In order to further verify the position and shape of goaf obtained by the comprehensive geophysical prospecting method, the borehole television method [29] is used to observe the possible goaf in exploration areas A and B. FB580 four-picture borehole television is used for borehole television observation. The borehole television has a high-definition camera in front of the probe and three high-definition cameras in the side of the probe. It can record the situation in the hole in all directions under the condition of water. Borehole television can also observe the development height of the caving zone and fissure zone above goaf and the water level in the hole. FB580 borehole television and construction site are shown in Figure 16.Figure 16
(a) FB580 borehole television; (b) construction site.
(a)
(b)Figure17 shows the television observation results of boreholes in exploration area A. Figure 17(a) is a borehole television image with the probe lowered to a depth of 26.13 m. It can be seen that there is no water in the borehole. Figure 17(b) shows the image of the probe lowered to a depth of 49.83 m. It can be seen that there is water leaking in the hole. Figure 17(c) is the image of the probe lowered to a depth of 55.22 m. It can be seen that the water in the hole increases obviously. Figure 17(d) is the image of the probe lowered to a depth of 56.82 m. It can be seen that there is obvious water jet phenomenon in the hole and the water pressure is large. Figure 17(e) is the image of the probe lowered to a depth of 92.28 m. The water surface in the hole can be clearly seen. Figure 17(f) is the image of the probe lowered to a depth of 93.49 m. The reflection image of the probe light source on the water surface can be clearly seen. Therefore, the probe depth at this time can be determined as the water level of the goaf, and the water level elevation is +371.21 m. Figure 17(g) is the image of the probe lowered to a depth of 94.23 m. When the probe enters the water, the lateral camera of the probe starts to work, and the bedding surface of the rock mass in the hole wall can be seen. Figure 17(h) is the image of the probe lowered to a depth of 115.54 m. The bedding plane of rock mass in the hole wall can be clearly seen.Figure 17
Television observation results of boreholes in exploration area A.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure17(i) is an image of the probe lowered to a depth of 130.35 m, with the elevation of +334.35 m. It is evident that the strata separating fracture can be seen. From the conclusion of borehole television observation in the literature [29], it can been seen that microcracks had been progressively generated in the wall area. The majority of these cracks were not initiated from the preexisting mineral filled joints but appeared at locations where no joints had been present earlier. In addition, the site drilling does not return slurry at this elevation. It is inferred that the location is the upper bound of the water conducting fracture zone above the goaf, and the height of the water conducting fracture zone is about 44 m. Figure 17(j) is the image of the probe lowered to a depth of 142.33 m. There are irregular fractures in the rock strata of the borehole. Figure 17(k) is an image of the probe lowered to a depth of 156.16 m, with the elevation of +308.54 m. Local block fall occurs in the rock strata in boreholes. Therefore, it is inferred that the location is the upper limit of the caving zone above the goaf, and the height of the caving zone is about 18 m. Figure 17(l) is the image of the probe lowered to a depth of 162.84 m. The rock blocks piled up in disorder, gradually unable to see the surrounding situation in the borehole.According to the existing mining geological data, the floor elevation of the goaf is +284 m at the drilling position, and the maximum mining height of the coal seam is 6 m. Therefore, the coal seam roof elevation is about +290 m.Figure18 shows the television observation results of boreholes in exploration area B. Figure 18(a) is a borehole television image with the probe lowered to a depth of 51.96 m. It can be seen that there is no water in the borehole. Figure 18(b) shows the image of the probe lowered to a depth of 60.01 m. It can be seen that there is water leaking in the hole.Figure 18
Television observation results of boreholes in exploration area B.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure18(c) is the image of the probe lowered to a depth of 62.05 m. It can be seen that the water in the hole increases obviously. Figure 18(d) is the image of the probe lowered to a depth of 69.07 m. It can be seen that there is obvious water jet phenomenon in the hole and the water pressure is large. Figure 18(e) is the image of the probe lowered to a depth of 88.67 m. The water surface in the hole can be clearly seen. Figure 18(f) is the image of the probe lowered to a depth of 91.61 m. The reflection image of the probe light source on the water surface can be clearly seen. Therefore, the probe depth at this time can be determined as the water level of the goaf, and the water level elevation is +371.49 m. Figure 18(g) is the image of the probe lowered to a depth of 105.32 m. When the probe enters the water, the lateral camera of the probe starts to work, and the bedding surface of the rock mass in the hole wall can be seen. Figure 18(i) is an image of the probe lowered to a depth of 160.21 m, with the elevation of +302.89 m. It is evident that the strata separating fracture can be seen. The site drilling does not return slurry at this elevation. It is inferred that the location is the upper bound of the water-conducting fracture zone above the goaf, and the height of the water-conducting fracture zone is about 65 m. Figure 18(j) is the image of the probe lowered to a depth of 165.32 m. There are irregular vertical secondary fissures in the strata of boreholes. Figure 18(k) is an image of the probe lowered to a depth of 198.86 m, with the elevation of +264.24 m. Local block fall occurs in the rock strata in boreholes. Therefore, it is inferred that the location is the upper limit of the caving zone above the goaf, and the height of the caving zone is about 26 m. Figure 18(l) is the image of the probe lowered to a depth of 207.78 m. The rock blocks piled up in disorder, gradually unable to see the surrounding situation in the borehole.According to the existing mining geological data, the floor elevation of the goaf is +230 m at the drilling position, and the maximum mining height of the coal seam is 8 m. Therefore, the coal seam roof elevation is about +238 m.


## 3.1. Transient Electromagnetic Method

The transient electromagnetic method (TEM) belongs to the time-domain electromagnetic induction method. Before 1990s, it was mainly used to search for metal deposits, coalfields, petroleum, gas, and other fields. Its basic principle is to send a pulse magnetic field (or electric field) to the underground by using an ungrounded loop (or a couple source), that is, to send a current pulse square wave on the transmitting loop, and the square wave will dropinstantaneously. Generate a magnetic field that can propagate underground. Under the excitation of the primary magnetic field, the geological body will produce eddy current, whose size depends on the conductivity of the geological body. After primary magnetic field disappears, the eddy current can not disappear immediately. It will have an excessive process (attenuation process), which produces attenuated secondary magnetic field to propagate to the surface. The secondary magnetic field will be received by the receiving line or magnetic probe of the surface instrument. The variation of the secondary magnetic field will reflect the electrical distribution of underground geological bodies. Secondary induction electromotive force (SIEMF) is measured with different delay times in the receiver, so the attenuation characteristics of the secondary magnetic field with time are obtained. The pseudosection map of resistivity inversion can be obtained by calculating program, and then the geological condition of deep underground can be judged. The advantages of TEM are easy handling, strong adaptability to surface sites, and high efficiency [20].


## 3.2. Controlled Source Audio-Frequency Magnetotellurics

For a long time, people can use various exploration methods to solve the engineering geology and hydrogeology problems concerned about the shallow surface geology of the engineering area, but there has been a lack of effective methods to understand the deep geological conditions. In geophysical exploration, electrical prospecting has been paid attention to because of its variability of field sources and effectiveness of methods. Controlled source audio-frequency magnetotellurics (CSAMT) is a new electrical method developed in the late 1980s [26]. Based on electromagnetic wave propagation theory and Maxwell equations, Cagniard’s apparent resistivity [27, 28] equation is derived as follows:(1)ρaC=1μωExHy2,where E and H are the orthogonal electric and magnetic fields, respectively. ω and μ denote angular frequency and magnetic permeability. If the electric field is measured in millivolt/km and magnetic field is measured in nT (nanotesla), the above equation becomes(2)ρaC=15fExHy2.According to the skin effect theory of the electromagnetic wave, the formula of skin depth can be derived as follows:(3)δ=2ωμσ=503ρf,where ρ is the electrical resistivity and f is the frequency.It has been shown that when the surface resistivity is fixed, the propagation depth (or detection depth) of the electromagnetic wave is inversely proportional to the frequency. In high frequency, the detection depth is shallow; in low frequency, the detection depth is deep. Therefore, the detection depth can be changed by adjusting the transmission frequency to achieve the purpose of detection.The general principle of the CSAMT method is to use controllable artificial field source and measure the Cagniard’s apparent resistivity by the ratio of electric field to magnetic field. At the same time, based on the principle of skin depth of the electromagnetic wave, electrical sounding of different depths is carried out by changing frequency. Therefore, the CSAMT method has the following advantages in its work: it is suitable for various sites, high efficiency, and can receive seven electric fieldsE and one magnetic field H at the same time in a single launch; it has strong anti-interference ability; it has a large exploration depth range, up to 1,000 m; it can detect the resistivity discrepancies in the geological background, has small high-resistance shielding effect, and can penetrate high-resistance strata. It is precisely because of these advantages that it was rapidly popularized and applied in engineering survey all over the world in the late 1990s.V8 multifunctional electromagnetic prospecting system is used, and the main technical parameters of the instrument are as follows: input impedance >100 MΩ, compensation error <0.12‰, resolution: 1μV, power: 2800 W, maximum output current: 40 A, and maximum supply voltage: 140 V.


## 3.3. Data Processing Method

The data collected on site must be further processed before they can be transformed into the required geophysical and geological results. The main parameters of TEM field observation are sampling time, emission current, and normalized induction secondary field, etc. Data processing begins with preprocessing, that is, data cataloguing, quality inspection, rejection of unqualified data, data quality statistics and classification, and sorting out the sequence and format required for special software. Secondly, the data are filtered to remove the suppressed interference signals and restore the original characteristics of geological signals. Then, the parameters are converted; that is, the apparent resistivity and apparent depth are obtained by software conversion. On this basis, topographic correction and other processing are done according to survey, geology, and drilling data. Finally, various sectional diagrams are drawn for data analysis.There are large bucket excavator system and material conveying tape in exploration area A. Its working voltage is about 6000 v, which makes the data near the operation area subject to certain electromagnetic interference. In addition, other large-scale operation vehicles in the exploration area will also have a certain impact on data acquisition. The data affected by electromagnetic interference will have abrupt maximum or minimum values, and even the phenomenon of local data uplifts or decreases. Figure7 is the contrast chart of the single-point quadratic field attenuation curve in double logarithmic coordinates. The attenuation curve without electromagnetic interference tends to decay in turn, and the data points basically do not show sudden change. The first part of the attenuation curve of the quadratic field of the measured point under electromagnetic interference basically maintains the trend of attenuation in turn. However, after 5 ms, the data show obvious “jitter and jump” phenomenon and the original attenuation law is broken.Figure 7
Comparison of single-point attenuation curves of transient electromagnetic interference.Correction is needed for the disturbed data. Hundreds of sampling points were extracted, and the interfered data were processed by combining the method of limit difference, fitting, and multipoint smoothness. The comparison of resistivity sections before and after processing is shown in Figure8, in which Figure 8(a) is the apparent resistivity section of data subject to electromagnetic interference and Figure 8(b) is the apparent resistivity section after interference correction. After correction, the change rule of affected data is consistent with that of data not subject to electromagnetic interference.Figure 8
Comparison of pseudosection of resistivity (a) before and (b) after interference correction.
(a)
(b)The measured TEM data are the comprehensive reflection of the underground electrical layers. Vertically, the original curve or the original apparent resistivity section can only show the comprehensive electrical change trend of several sets of strata. To interpret the distribution of the real strata or the closest to the real strata and geological anomalies through the measured data, the data should be inverted to obtain the subsurface resistivity distribution.Figure9 shows the different data processing methods for transient electromagnetic exploration. From Section 2.2, it can be seen that the strata in the exploration area present the electrical law of “high-low-high-low” from shallow to deep. By comparing (a) and (b) in Figure 9, it can be concluded that the section formed by inversion results reflects the electrical property law of the formation from shallow to deep “low-high-low-high,” while the section of the late formula shows the change characteristic of “high-low” from shallow to deep, which is not consistent with the actual electrical property of the formation. In addition, although the late formula has the advantages of simplicity and rapidity, the algorithm is suitable for the late time. There is some distortion in the early calculation, which results in unreal high resistivity response in the shallow section, and the inversion calculation restores the actual electrical property of the formation. Therefore, the inversion method is chosen as the data processing method of TEM exploration.Figure 9
Different data processing methods for transient electromagnetic exploration. (a) Late formula for calculating the resistivity section. (b) Inversion method for calculating the resistivity section.
(a)
(b)The data processing flow of controlled source audio-frequency magnetotellurics (CSAMT) is shown in Figure10. Detailed data processing can be referred to [26–28].Figure 10
Data processing flow chart of the CSAMT method.


## 3.4. Results of Integrated Geophysical Prospecting

By using the integrated geophysical prospecting method (TEM and CSAMT), the mined-out area and water accumulation area of 6 coal seams in the exploration area were comprehensively explored. If the stratum is not affected by the goaf or the water-conducting structure, the resistivity of coal-bearing stratum changes orderly, and the isoline changes steadily in the resistivity section, showing a layered distribution. When there is a goaf water accumulation area or water-conducting structure, the resistivity of anomalous areas decreases, and the isoline distribution is distorted, deformed, or densely striped. Typical results for exploration areas A and B are shown in Figures11–15, respectively. The location of the selected exploration line is shown in Figure 3.Figure 11
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 5 of the TEM method. (b) Inversion of the resistivity profile from line 5 of the CSAMT method.
(a)
(b)Figure 12
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 9 of the TEM method. (b) Inversion of the resistivity profile from line 9 of the CSAMT method.
(a)
(b)Figure 13
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 13 of the TEM method. (b) Inversion of the resistivity profile from line 13 of the CSAMT method.
(a)
(b)Figure 14
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 31 of the TEM method. (b) Inversion of resistivity profile from line 31 of CSAMT method.
(a)
(b)Figure 15
Resistivity sections obtained by different methods. (a) Inversion of the resistivity profile from line 33 of the TEM method. (b) Inversion of the resistivity profile from line 33 of the CSAMT method.
(a)
(b)Figure11 shows the resistivity sections of line 5 obtained by different methods. The survey line is located in the southern boundary of the exploration area. The length of the survey line is 1180 m, and the surface roughly flat. Coal seams gradually deepen from west to east. It can be seen that the vertical resistivity of strata in the range of exploration depth generally shows the trend of “high resistivity-low resistivity-high resistivity-low resistivity,” which coincides with the electrical characteristics of the actual stratum described in Section 2.2. The horizontal resistivity isoline of strata near the coal seam is obviously distorted between No. 7 and No. 33 survey point, showing obvious characteristics of low-resistance anomaly. In addition, there is a normal fault F4 near No. 27 survey point. The fault distance is 30 m, which is not suitable for mining. Therefore, the anomaly can be inferred as a fractured water-rich zone. The horizontal resistivity isoline of strata near the coal seam is obviously distorted between No. 39 and No. 43 survey point, also showing obvious characteristics of low-resistance anomaly. As shown in Figure 11(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method. It is also shown that the electrical characteristics and the resistivity anomalies delineated by both methods are in agreement.Figure12 shows the resistivity sections of line 9 obtained by different methods. The survey line is located in the middle of exploration area A. The length of the survey line is 840 m, and the topography is relatively flat. The coal seam is gradually deepening from west to east from No. 25 to No. 47, and the whole coal seam is slightly uplifted from No. 48 to No. 66. From Figure 12(a), it can be seen that the vertical resistivity of the stratum in the exploration depth generally shows a trend of “high resistivity to low resistivity to high resistivity to low resistivity,” which coincides with the electrical characteristics of the actual stratum described in Section 2.2. The cross-sectional map is bounded by No. 50 survey point, and the rise and fall of the target electrical layer are obvious. From the trend of the coal seam, the trend of the left isoline is basically consistent with the trend of the actual stratum. Laterally, between No. 25 and No. 43 survey point, the resistivity isoline of the seam near the No. 6 coal seam is obviously distorted, showing obvious characteristics of low resistivity anomaly (the magenta ellipse in the picture). It is inferred that the low-resistance anomaly may be related to water accumulation in goaf formed by crossborder mining of 6 coals in adjacent coal mines. The isoline on the right side of the section (i.e., between the 50–63 survey points) changes obviously, the vertical high resistivity layer rises obviously upward, and the lower part of the section shows obvious low-resistance anomaly. It is inferred that it is the electrical reflection of No. 6 coal mined out and accumulated water. The upper high-resistance anomaly is similar to the electrical reflection of the nonfilling cracks formed by the caving of overlying strata after coal mining. As shown in Figure 12(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method.Figure13 shows the resistivity sections of line 13 obtained by different methods. The survey line is located in the north of exploration area A. The length of the survey line is 720 m. The surface is nearly flat, and the undulation is small. From Figure 13(a), it can be seen that the vertical resistivity of the strata in the exploration depth generally shows a trend of “high resistivity to low resistivity to high resistivity to low resistivity,” which is consistent with the electrical characteristics of the actual strata mentioned above. The resistivity isoline of the stratum near the coal seam is obviously distorted and shows obvious low-resistance anomaly between No. 36 and No. 46 as well as between No. 54 and No. 65 survey point transversely. It is inferred that the low-resistance anomaly may be related to water accumulation in goaf formed by crossborder mining of 6 coals in adjacent coal mines. As shown in Figure 13(b), the inversion resistivity profile obtained by the CSAMT method is basically consistent with that obtained by the TEM method, and the anomalies resistivity variations obtained by both methods are in agreement.The analysis of comprehensive results of some survey lines in exploration area A is provided earlier, and integrated geophysical exploration techniques are also carried out in exploration area B. Typical results are shown in Figures14 and 15. The low resistivity anomalous area near the coal seam can be obtained by comparing the modelling results of the two methods. It can be seen that the electrical characteristics and abnormal reactions obtained by the two geophysical methods are basically consistent. The results of deduced water-filled goaf and water-rich fracture area in exploration areas A and B are shown in Figure 3.


## 3.5. Verification of Borehole Television Method

In order to further verify the position and shape of goaf obtained by the comprehensive geophysical prospecting method, the borehole television method [29] is used to observe the possible goaf in exploration areas A and B. FB580 four-picture borehole television is used for borehole television observation. The borehole television has a high-definition camera in front of the probe and three high-definition cameras in the side of the probe. It can record the situation in the hole in all directions under the condition of water. Borehole television can also observe the development height of the caving zone and fissure zone above goaf and the water level in the hole. FB580 borehole television and construction site are shown in Figure 16.Figure 16
(a) FB580 borehole television; (b) construction site.
(a)
(b)Figure17 shows the television observation results of boreholes in exploration area A. Figure 17(a) is a borehole television image with the probe lowered to a depth of 26.13 m. It can be seen that there is no water in the borehole. Figure 17(b) shows the image of the probe lowered to a depth of 49.83 m. It can be seen that there is water leaking in the hole. Figure 17(c) is the image of the probe lowered to a depth of 55.22 m. It can be seen that the water in the hole increases obviously. Figure 17(d) is the image of the probe lowered to a depth of 56.82 m. It can be seen that there is obvious water jet phenomenon in the hole and the water pressure is large. Figure 17(e) is the image of the probe lowered to a depth of 92.28 m. The water surface in the hole can be clearly seen. Figure 17(f) is the image of the probe lowered to a depth of 93.49 m. The reflection image of the probe light source on the water surface can be clearly seen. Therefore, the probe depth at this time can be determined as the water level of the goaf, and the water level elevation is +371.21 m. Figure 17(g) is the image of the probe lowered to a depth of 94.23 m. When the probe enters the water, the lateral camera of the probe starts to work, and the bedding surface of the rock mass in the hole wall can be seen. Figure 17(h) is the image of the probe lowered to a depth of 115.54 m. The bedding plane of rock mass in the hole wall can be clearly seen.Figure 17
Television observation results of boreholes in exploration area A.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure17(i) is an image of the probe lowered to a depth of 130.35 m, with the elevation of +334.35 m. It is evident that the strata separating fracture can be seen. From the conclusion of borehole television observation in the literature [29], it can been seen that microcracks had been progressively generated in the wall area. The majority of these cracks were not initiated from the preexisting mineral filled joints but appeared at locations where no joints had been present earlier. In addition, the site drilling does not return slurry at this elevation. It is inferred that the location is the upper bound of the water conducting fracture zone above the goaf, and the height of the water conducting fracture zone is about 44 m. Figure 17(j) is the image of the probe lowered to a depth of 142.33 m. There are irregular fractures in the rock strata of the borehole. Figure 17(k) is an image of the probe lowered to a depth of 156.16 m, with the elevation of +308.54 m. Local block fall occurs in the rock strata in boreholes. Therefore, it is inferred that the location is the upper limit of the caving zone above the goaf, and the height of the caving zone is about 18 m. Figure 17(l) is the image of the probe lowered to a depth of 162.84 m. The rock blocks piled up in disorder, gradually unable to see the surrounding situation in the borehole.According to the existing mining geological data, the floor elevation of the goaf is +284 m at the drilling position, and the maximum mining height of the coal seam is 6 m. Therefore, the coal seam roof elevation is about +290 m.Figure18 shows the television observation results of boreholes in exploration area B. Figure 18(a) is a borehole television image with the probe lowered to a depth of 51.96 m. It can be seen that there is no water in the borehole. Figure 18(b) shows the image of the probe lowered to a depth of 60.01 m. It can be seen that there is water leaking in the hole.Figure 18
Television observation results of boreholes in exploration area B.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure18(c) is the image of the probe lowered to a depth of 62.05 m. It can be seen that the water in the hole increases obviously. Figure 18(d) is the image of the probe lowered to a depth of 69.07 m. It can be seen that there is obvious water jet phenomenon in the hole and the water pressure is large. Figure 18(e) is the image of the probe lowered to a depth of 88.67 m. The water surface in the hole can be clearly seen. Figure 18(f) is the image of the probe lowered to a depth of 91.61 m. The reflection image of the probe light source on the water surface can be clearly seen. Therefore, the probe depth at this time can be determined as the water level of the goaf, and the water level elevation is +371.49 m. Figure 18(g) is the image of the probe lowered to a depth of 105.32 m. When the probe enters the water, the lateral camera of the probe starts to work, and the bedding surface of the rock mass in the hole wall can be seen. Figure 18(i) is an image of the probe lowered to a depth of 160.21 m, with the elevation of +302.89 m. It is evident that the strata separating fracture can be seen. The site drilling does not return slurry at this elevation. It is inferred that the location is the upper bound of the water-conducting fracture zone above the goaf, and the height of the water-conducting fracture zone is about 65 m. Figure 18(j) is the image of the probe lowered to a depth of 165.32 m. There are irregular vertical secondary fissures in the strata of boreholes. Figure 18(k) is an image of the probe lowered to a depth of 198.86 m, with the elevation of +264.24 m. Local block fall occurs in the rock strata in boreholes. Therefore, it is inferred that the location is the upper limit of the caving zone above the goaf, and the height of the caving zone is about 26 m. Figure 18(l) is the image of the probe lowered to a depth of 207.78 m. The rock blocks piled up in disorder, gradually unable to see the surrounding situation in the borehole.According to the existing mining geological data, the floor elevation of the goaf is +230 m at the drilling position, and the maximum mining height of the coal seam is 8 m. Therefore, the coal seam roof elevation is about +238 m.


## 4. Results and Discussion

The distribution of mined-out areas and water accumulation in the coal mine are investigated in detail by using the comprehensive geophysical prospecting method of the transient electromagnetic method (TEM) and controlled source audio-frequency magnetotellurics (CSAMT). Finally, the distribution of coal mined-out areas in the exploration area is obtained, and the zoning map of comprehensive geophysical prospecting achievements in exploration areas A and B are shown in Figure3. Based on the results of borehole exploration and the comprehensive analysis of the borehole image, it is concluded that the instantaneous water level elevation in goaf of exploration area A is about +371 m, the development height of the water-conducting fracture zone above goaf is about 44 m, and the height of the caving zone is about 18 m. In addition, the instantaneous water level elevation in goaf of exploration area B is about +371 m, the development height of the water-conducting fracture zone above goaf is about 65 m, and the height of the caving zone is about 26 m. The results of various geophysical methods are mutually corroborated, and the results are in good agreement.In addition, combined with the conclusions of similar detection methods in the past [30–32], it is found that the resolution and sensitivity of TEM and CSAMT are both high, and they can meet the necessary technical parameters for exploration and detection of underground goaf. Its exploration effect depends on the device type, receiving frequency, stacking times, and other related technical parameters. The underground goaf does show low-resistance anomaly according to the length of development time of the underground cavity and the degree of water accumulation [30]. It shows remarkable characteristics on the apparent resistivity section. Under the guarantee of suitable construction parameters and combined with existing geological data, the development characteristics of underground voids can be more accurately detected, and the survey effect is obvious [30, 31]. It provides a good example for the application of the comprehensive geophysical prospecting method of TEM and CSAMT in the exploration of the underground cavity development area. And it has strong reference significance for the investigation and prevention of collapse geological environment disasters caused by such underground cavity.


## 5. Conclusions

The distribution of mined-out areas and water accumulation in the coal mine are investigated in detail by using the comprehensive geophysical prospecting method. Firstly, the unmanned aerial vehicle (UAV) aerial survey is used to accurately acquire the three-dimensional model of the slope in the key area of Yuanbaoshan open-pit coal mine, which provides basic information for the real-time visualization of the open-pit mine and the establishment of the engineering geological research profile. Secondly, the distribution of goaf in exploration areas A and B is detected by using the transient electromagnetic method (TEM) and controlled source audio-frequency magnetotellurics (CSAMT) combined geophysical methods. The results show that the electrical characteristics and low-resistance anomaly reflected by the two geophysical methods are in good agreement with each other. Finally, the borehole detection method is used to verify the goaf range obtained from comprehensive geophysical prospecting, and the results show good consistency. In addition, based on the results of borehole exploration and the comprehensive analysis of the borehole image, it is concluded that the instantaneous water level elevation in goaf of exploration area A is about +371 m, the development height of the water conducting fracture zone above goaf is about 44 m, and the height of the caving zone is about 18 m. In addition, the instantaneous water level elevation in goaf of exploration area B is about +371 m, the development height of the water-conducting fracture zone above goaf is about 65 m, and the height of the caving zone is about 26 m. Overall, the results of comprehensive exploration show that through the use of the TEM and CSAMT method, combined with drilling detection methods, the purpose of accurate exploration of mined-out areas in coal mines can be achieved. Therefore, this comprehensive geophysical exploration research method provides a feasible and effective technical path for open-pit mines facing similar problems. It can be used for reference and extended to other mines or underground engineering exploration research.



---

*Source: 2368402-2019-08-08.xml*