Next Article in Journal
Bioaugmentation of Aerobic Granular Sludge with Dye-Decolorizing Yeast for Textile Industrial Wastewater
Next Article in Special Issue
Study on Instability Characteristics of the Directional Borehole on the Coal-Seam Roof: A Case Study of the Tingnan Coal Mine
Previous Article in Journal
Influence of Sodium New Houttuyfonate as a New EGFR-TK Inhibitor on the Apoptosis and Autophagy of MCF-7 Cells and Its Toxicity to Caenorhabditis elegans
Previous Article in Special Issue
A Case Study on Strong Strata Behaviors Mechanism of Mining Reserved Roadway and Its Prevention Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Geophysical Prediction of Goaf and Water Accumulation in Pingshuo Dong Open-Cut Mine with Ultrashallow and High Drops

1
Department of Earth Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
2
Shanxi Institute of Geological Survey, Taiyuan 030006, China
3
No.2 Oil Production Plant, Xinjiang Oilfield Company, PetroChina, Karamay 834008, China
4
Coal Geological Geophysical Exploration Surveying & Mapping Institute of Shanxi Province, Jinzhong 030600, China
*
Author to whom correspondence should be addressed.
Processes 2023, 11(6), 1653; https://doi.org/10.3390/pr11061653
Submission received: 16 March 2023 / Revised: 15 May 2023 / Accepted: 19 May 2023 / Published: 29 May 2023
(This article belongs to the Special Issue Process Safety in Coal Mining)

Abstract

:
To address safety problems caused by goaf and water accumulation in open-pit mines, the shallow three-dimensional seismic method and transient electromagnetic method were integrated and applied to detect the mining goaf distribution scope and the water accumulation conditions. In view of the special topographic conditions of an ultrashallow layer and high drop in an open-pit mine, we proposed utilizing bin homogenization and multidomain joint denoising to improve the reliability of seismic data. By using seismic-sensitive attributes to predict the goaf, the transient electromagnetic method was employed to further predict the water accumulation in the goaf. The results show that the shallow seismic method could clearly reflect the reflected-wave variation features of the goaf. The features of a junction between a normal seam and goaf vary obviously, and the prediction effect of the goaf boundary with high resolution is in line with the actual situation. Furthermore, taking the goaf scope prediction with the shallow three-dimensional seismic method as a base, targeted transient electromagnetic detection was deployed, with a detailed analysis of the survey-line repeated areas of the shallow seismic and transient electromagnetic methods. Making full use of the advantages of the shallow seismic and transient electromagnetic methods, we propose a reasonable data interpretation method in combination with the special topographic conditions of open-pit mines, which greatly improves the prediction effect of goaf and water accumulation conditions.

1. Introduction

Since the founding of New China in 1949, nationwide construction has spurred demand for coal and the exploitation of coal resources. Currently, coal production in China accounts for almost half of the world’s total production [1]. Since 2005, to improve the efficient utilization and safe mining of coal resources, China has carried out the effective integration of coal resources [2]. However, most coal mines with integrated resources exhibit problems, such as a serious lack of data on the old goaf, unclear distribution of the boundary area, unclear water area and quantity in the goaf, and difficulty carrying out geophysical prospecting work underground in small closed-pit coal mines, which have exposed serious hidden dangers to coal mine safety. The Pingshuo Dong open-cast coal mine is a typical example. Due to long-term disorderly mining, no-filling mining, rich mining, and poor abandonment of surrounding small coal mines, many underground goafs exist in the Pingshuo Dong open-cut coal mine, and most of these goafs have no recorded specific information, such as location and size. If the roof thickness of the goaf is less than the minimum safe thickness, a roof collapse can occur, which can cause surface subsidence or cracking, seriously threatening the safety of personnel and equipment working above the roof [3,4].
Geophysical prospecting methods commonly used in goaf and water accumulation include the shallow seismic method [5,6], transient electromagnetic method [7,8,9], high-density electrical method [10,11], magnetotelluric method [12,13], and ground-penetrating radar [14,15]. Various geophysical methods have achieved certain effects in detecting goaf and water accumulation, but each geophysical method has its own applicable conditions and advantages. The shallow seismic exploration method allows continuous measurement, and the range of goaf detected is more accurate and reliable, but it can hardly reflect the accumulated water of the goaf area [16]. The transient electromagnetic method is especially sensitive to low resistance and has a good effect in detecting goafs with accumulated water, and the detection depth is relatively large, but it is greatly affected by electromagnetic interference. The high-density resistivity method has a better detection effect on goaf water accumulation, but it is affected by topography and grounding resistance [17]. The magnetotelluric method is effective at detecting water accumulation in goafs, and the detection depth is large, but it can be affected by static effects and near-field effects. Because a single geophysical method is limited, the geophysical interpretation often exhibits multiple solutions and low reliability [18,19]. In recent years, many studies and applications have been carried out for integrated geophysical methods [20,21,22,23]. Integrated geophysical prospecting can obtain more physical information, verify and supplement other methods, reduce the multiple solutions of a single geophysical prospecting method, and improve the resolution of exploration and the reliability of interpretation results. However, for open-pit mining areas, the terrain conditions are more complicated, the coal seams are very shallow, and there are high-drop steps (transport lines). For goaf distribution and water accumulation prediction, the high-density electric method cannot meet the condition of the line layout, the ground-penetrating radar detection depth is insufficient, and the shallow seismic acquisition coverage is uneven. How to achieve higher resolution and better exploration results has not been systematically discussed in depth. This paper takes the prediction of the goaf distribution and water accumulation of the Pingshuo Dong open-cast coal mine as an example. Aiming at this type of problem, this paper proposes combining small-array high-density acquisition, multidomain joint denoising, and bin homogenization, weakening the influence of high-drop steps on shallow seismic exploration, and then using the multiple-attribute method to predict the goaf. Based on the electromagnetic environment analysis and topographical features of the study area, the transient electromagnetic method was used to identify goafs. The practical application shows that the high-density 3D seismic and transient electromagnetic comprehensive interpretation can effectively identify the goaf and water accumulation characteristics of an ultrashallow high-altitude open-pit mine. This proves the feasibility and application of the proposed method and can provide a meaningful reference for research in similar areas.

2. Geological Background and Engineering Geological Conditions

The study area is located in the northeastern margin of the Ningwu coalfield in Shanxi, China. Most of the surface is covered by Cenozoic strata, which are typical loess hilly landforms. There are sporadic strata bedrock outcrops near the coal-seam outcrops in the northeastern part of this area and at the bottoms of major gullies in this area. The study area is located on the southwest side of the eastern open-cast coal mine and its surrounding area, and the exploration area is approximately 0.35 km2, the scope is delineated by a red circled. as shown in Figure 1.
The Pingshuo Dong open-cast coal mine is located in the low hills of the Pingshuo Platform in the Shanxi Loess Plateau. The whole area is mostly covered with loess. The loess plateau in this area has been strongly eroded and cut. The valley is relatively developed in a “V” shape with a cutting depth of 40–70 m. The highest point in the area is located in the northern part of Shifeng, with an altitude of 1537.99 m. The lowest point is located at the Maguan River in western Xuegaodeng village, with an altitude of 1213.8 m, and the maximum relative height difference is 324.19 m. The ground surface in the study area has been stripped, showing a step-like shape, with an elevation of 1190–1380 m, a minimum of 1190 m in the northeast corner, a height of 1380 m in the south, and a relative height difference of 190 m. Some areas are the main roads in the southern part of the mine, and the traffic is more intensive, as shown in Figure 2.

3. Methods

3.1. Attribute-Based Predicted Multiparameter Analysis for Goaf Prediction

The geophysical characteristics of the goaf and its upper strata have changed compared with those before mining, including differences in density, velocity and other elastic parameters, which provide the physical basis for detecting goafs by the seismic method [24]. A large amount of data proves that the existence of goafs changes the characteristics of coal-seam reflection waves, which is manifested in the distortion, dislocation, and even disappearance of coal-seam reflection-wave groups, the reduction in seismic amplitude, frequency, phase reversal, and waveform changes. These characteristics are the main signs for identifying goafs. Currently, seismic attributes are among the most effective ways to identify goafs.
The extraction of different seismic attributes is essentially the selection of different attribute parameters. Different geological characteristics correspond to different physical parameters. Therefore, selecting reasonable seismic attributes in a wide range of seismic attributes depends on the characteristics of geological bodies or geological phenomena. To improve the accuracy of identifying goafs, we analyzed the geological characteristics of goaf development and its seismic response and selected several layer attributes that are extremely sensitive to goafs and small faults from more than 20 attributes, including average amplitude energy, instantaneous frequency, dominant frequency, instantaneous phase, and waveform-clustering attributes.
Vakman compared various methods for extracting instantaneous parameters of signals and concluded that the analytical signal method is the most reasonable [25]. The analytical signal method uses the Hilbert transform to obtain the corresponding analytical signal to form a complex signal and then extracts the instantaneous parameters of the signal.
Assuming that the real signal is x ( t ) , the corresponding imaginary part can be obtained by the Hilbert transform [26]:
y ( t ) = H T [ x ( t ) ] = 1 π x ( τ ) t τ d τ
H T [ ] represents the Hilbert transform. Then, we can construct an analytic signal s ( t ) by x ( t ) and y ( t ) :
s ( t ) = x ( t ) + i y ( t ) = A ( t ) e i θ ( t )
According to the analytic signal, the instantaneous parameters of the original signal x ( t ) can be extracted, including
A ( t ) = x 2 ( t ) + y 2 ( t )
ϕ ( t ) = t a n 1 y ( t ) / x ( t )
ω ( t ) = 1 2 π d ϕ ( t ) d t = 1 2 π x ( t ) y ( t ) y ( t ) x ( t ) x 2 ( t ) + y 2 ( t )
where A(t), ϕ ( t ) , and ω ( t ) represent the instantaneous amplitude, instantaneous phase, and instantaneous frequency, respectively.
The instantaneous amplitude is a measure of the reflective intensity, which is the square root of the total energy of the real and imaginary parts of the analytical signal. It mainly reflects the change in energy and can highlight the change in special rock strata. The instantaneous phase is a measure of the continuity of seismic events in seismic profiles, and its phase can be displayed regardless of the intensity of energy, even if the weak-amplitude effective wave can be well displayed on the instantaneous phase diagram. When the wave propagates in an anisotropic homogeneous medium, its phase is continuous when the wave is in an abnormal medium. Its phase will change significantly in the abnormal position, and it is obviously discontinuous in the profile. The instantaneous frequency is the rate of time variation of the phase, which can reflect the lithological changes in the stratum and help to identify the stratum.
For the same detection object, the obvious change in three instantaneous pieces of information in the same location may reflect a change in the physical properties of the detection object. Among the three parameters, the instantaneous phase spectrum has the highest resolution. The change in the instantaneous frequency and amplitude spectrum is also more intuitive and can be used to determine the general location of underground anomalies and the instantaneous phase to determine the stratified boundary.
The C3 algorithm is a coherence algorithm based on feature structure analysis [27]. The calculation of coherence values in the algorithm can use multichannel seismic records. When calculating the coherence value at the sampling point (x, y, t), first define a data body containing the point (controlled by a two-dimensional spatial window and a one-dimensional time window), assuming that the time window contains J channels. For the seismic record, the number of sampling points in the time window is 2 N + 1; let d n T ( x , y , t ) = ( d n 1 , d n 2 , , d n J ) , n = N , N + 1 , , N 1 , N . The seismic trace in the data body can be represented as matrix.
D ( x , y , t ) = d N , d N + 1 , , d N T
The correlation matrix is calculated as follows:
C ( x , y , t ) = D T D = n = N n = N d n d n T = n = N n = N d n 1 2 n = N n = N d n 1 d n 2 n = N n = N d n 1 d n J n = N n = N d n 2 d n 1 n = N n = N d n 2 2 n = N n = N d n 2 d n J n = N n = N d n J d n 1 n = N n = N d n J d n 2 n = N n = N d n J 2
The coherence value Ec (x, y, t) at point (x, y, t) is defined as
E c ( x , y , t ) = λ max T r ( C ) = λ max j = 1 J C j j = λ max j = 1 J λ j
where λmax is the maximum value in λj.
It is easy to know from Equation (7) that the correlation matrix is positive definite or semipositive definite, so the value of Ec falls in the range of 1/J~J. When the waveforms of different traces are the same, Ec is equal to 1; when the waveforms are different, Ec is less than 1, and the greater the difference between waveforms, the smaller the value of Ec. The difference between waveforms can reflect the changes in geological bodies, so Ec can depict the discontinuities of geological bodies.

3.2. TEM-Based Apparent Resistivity Prediction for Water Accumulation Identification

The transient electromagnetic method (TEM) is a geophysical method that uses the principle of electromagnetic induction [28,29]. It emits pulsed electromagnetic waves as an excitation field source (called ’primary field’) through an ungrounded return line or ground-line-source electric dipole. It also observes and studies the spatial distribution and time characteristics of the electromagnetic field (called ’secondary field’) of the stratum or geological target body (underground goaf) generated by the excitation field (called ’primary field’), which can be used to infer and explain the spatial structure and physical characteristics of the stratum or goaf water accumulation.
Transient electromagnetic measurements are carried out by using the secondary eddy-current field excited by a step-pulse current in a large fixed source loop [30,31]. A rectangular emitter loop is placed on the surface of a homogeneous isotropic ground with magnetic conductivity and conductivity. The area of the loop is S, and the inner loop provides a step-pulse current. Before disconnecting the current, the emission current produces a stable magnetic field in the space around the loop and in the ground, as shown in Figure 3a.
When the current is turned off, the magnetic field generated by the emission current also disappears immediately. This change is quickly transmitted to the ground around the return line in the primary field with the aid of air and underground conductive medium and induces a current in the earth to maintain the magnetic field before the emission current is disconnected, so that the magnetic field of the space does not immediately disappear. However, because of the ohmic loss of the medium, the magnetic field generated by it will decay rapidly with the decay of the induced current. The rapid decay of the magnetic field will induce a new weaker eddy current in the surrounding medium. This process continues until the earth’s ohmic losses have exhausted the magnetic field energy (see Figure 3b), forming a geodetic electromagnetic field as the target of detection.

4. Application and Discussion

The high-density three-dimensional seismic and transient electromagnetic combined detection method is employed to identify the goafs in the detection area. Three-dimensional seismic design has 5843 survey points in the whole area, including 2399 shot points and 3444 detection points. The exploration area is 0.25 km2, and the coverage area of 24 times is 0.2 km2. The normal distance of the transient electromagnetic method is 26~35 m. The engineering survey line is arranged along the design roadway and its surroundings. The point distance of key areas above the roadway is 10 m, some areas are densified to 5 m, and the point distance of other areas is 20 m. The survey lines are arranged along the east–west axis, and a total of six survey lines are arranged. Two survey lines are arranged above the No. 4 coal roadway, numbered N0 and N35 from south to north. Four survey lines are arranged along the No. 9 coal roadway, numbered l0, L33, l59, and L85 from south to north, and the survey points are numbered 0, 20, and 40 from west to east according to the point distance, as shown in Figure 4.
The buried depth of Coal Seam 4 is 0–150 m, with an average of 70 m. The 9th coal seam has a buried depth of 40–210 m, with an average of about 110 m. The extremely shallow burial depth of the target layer and the small survey area led to the limited length of survey lines. Therefore, seismic exploration adopts the method of small array and high-density acquisition, with high shot-point density and high coverage times. After the field collection test, the shot interval has changed from 24 m to 6 m, and the coverage times have changed from 24 to 64. Meanwhile, due to the restriction of the field acquisition and construction conditions (vehicle walking platform), the distribution of shot and detection points is extremely uneven (Figure 5a), resulting in extremely uneven coverage of the study area (Figure 5b).
Due to the limitation of terrain conditions, such as large drop and steep slopes between the platforms, the original coverage times are extremely uneven. Therefore, the surface homogenization method is used to eliminate the adverse effects of uneven coverage, including static bin homogenization and dynamic bin homogenization. Static bin homogenization is used to make only one track in each offset group in a bin. If there are N tracks in an offset group in this bin, the N tracks in the offset are weighted and superimposed to represent the offset group; if there are N tracks in an offset group in the bin, then the N tracks are weighted and superimposed, and the superimposed track represents the offset group. Dynamic surface equalization is mainly used to find the offset group in the adjacent element. Define a large rectangular bin centered on the current bin. Any missing offset set in the current bin searches for the corresponding offset set in the adjacent bin within a circle, where the radius of the circle is a linear function of the offset set and its range is limited to the large bin.
During the processing, the seismic trace is edited on the single-shot record to eliminate the strong interference channel and the invalid redundant channel. On this basis, dynamic panel homogenization is used to supplement the missing offset set, and static panel homogenization is used to make the coverage times and energy more balanced (Figure 5c).
Figure 6 shows the CDP gathers before and after bin homogenization, where Figure 6a,b show the processing effects of bin-homogenizing the number of insufficient superpositions traced to the average number of superpositions. Figure 6c, d show the processing effects of bin-homogenizing the number of excessive superpositions traced to the average number of superpositions. From the processing effects, we find that the basic shape of the seismic profile has not changed. The overall energy of the seismic section is more balanced.
Figure 7 shows the comparison of stacked gathers before and after bin homogenization. We can determine that the energy of the missing trace is recovered (at the red elliptical mark), and the amplitude of the homogenized bin can reflect the real situation of the underground media.
Due to the influence of engineering construction, noise has a great influence on shallow seismic exploration. Suppressing noise is very important for the accurate identification of goafs. Data analysis shows that the main noise sources include linear interference, surface waves, and abnormally strong amplitudes. In this study, a multidomain joint denoising method is employed through experiments to classify and eliminate different types of noise. This follows the law of progressive, step-by-step denoising and gradually improves the signal-to-noise ratio. For linear interference, the apparent speed and frequency range of linear interference is analyzed, and the noise attenuation method with less waveform modification is adopted to eliminate the noise. For the outliers and abnormal amplitudes, the regional anomalous amplitude attenuation method and the frequency-division abnormal amplitude attenuation method are used for suppression. The strategy of multiple-trace statistics, single-trace denoising, and frequency-division suppression is adopted to suppress random abnormal amplitude noise and surface waves. Strong energy interference in seismic records is identified automatically in different frequency bands, and the spatial position of noise is determined. According to the defined threshold value and attenuation coefficient, time-varying and space-varying methods are adopted to suppress random abnormal amplitude noise. Figure 8 shows the comparison of single-shot records before and after suppressing noise. Within the scope of the red rectangular box, we can determine that the surface waves and other noises are effectively suppressed, and the energy of the primary reflection wave is obviously enhanced. Figure 9 is the comparison of stacked gathers before and after suppressing noise. Within the scope of the red elliptical frame, it can be seen that the weak waveform after noise suppression is displayed, which can more effectively reflect the subtle changes in underground media.
This study identifies goaf and water accumulation characteristics by three steps: seismic forward modeling assistant interpretation, goaf seismic-wave-anomaly characteristics interpretation, and seismic- and electrical-method comprehensive-anomaly interpretation. First, a goaf geological model suitable for the exploration area is established, and wave forward modeling is carried out to study the seismic response characteristics of various goafs (Figure 10) and assist in anomaly interpretation.
Due to the air or water content in the goaf, the absorption coefficient of the goaf is significantly higher than that of the normal coal seam, resulting in no reflection or weak reflection at the bottom interface of the goaf and substratum [32]. There is a fault in the wave group; the reflected wave is characterized by energy enhancement, frequency reduction, and an obvious concave phase (Figure 11a), or the energy becomes weaker and the continuity becomes worse (Figure 11b).
Comparing the reflected waves, seismic attributes, and drilling verification results of coal seams in the study area, a total of four coal-seam goafs are explained. The goaf labeled ① is located at the northern boundary of the exploration area, and the goaf anomaly is irregular. Combined with the interpretation results of the northern transient electromagnetic method, the goaf is produced by mining from north to south in the Yuling coal mine, and the abnormal area in the seismic work area is approximately 1359 m2.
The goaf marked with ② is located outside the southeast exploration boundary, in front of the Xiaoxiyao mine roadway, and the boundary of the abnormal area is irregular. The goaf is caused by mining from east to west in the Xiaoxiyao mine, and the abnormal area in the area is approximately 7352 m2.
The goaf marked with ③ and ④ is located southeast of the exploration boundary, and the abnormal boundary is irregular. The goaf is formed by the westward mining of the Xiaoxiyao mine. However, the northern part of the abnormal area is greatly disturbed by the raw coal transportation channel, and the southern part may be affected by the deterioration of data quality due to its close proximity to the boundary of the detection area, so the reliability of the results needs to be verified by more drilling results.
There is a seismic anomaly zone of coal-seam structure change in coal seam 9, which is located in the west of the exploration area. The plane shape is an irregular polygon, and the east side is connected with the seismic anomaly of the collapse column. The scope is delineated by a dotted line in Figure 12. The drilling data confirmed that the reason was that coal seams 9–1 and 9–2 merged into one layer.
A collapse column was identified via the drilling and log data from coal seam 4 (Figure 13). The lithology below K2 sandstone, as the marker bed of the No. 4 coal roof, is breccia, accompanied by a small amount of coal chips and blocks, with poor cementation. The shape of the collapse column on the plane is an approximate ellipse, especially in the instantaneous frequency attribute, and the characteristics of the collapse column are the most obvious. The collapsed column is a good water channel and water-accumulating area, so it is very important to identify the development of the coal-seam collapse column.
Based on the seismic interpretation of the development characteristics of special geological bodies such as goafs, faults, and collapse columns, the transient electromagnetic method is employed to predict the water accumulation in the key areas of the study area. The transient electromagnetic survey lines covering this area are N0 and N35. Due to terrain limitations, we adopted small coils and multiple circles of central loop to detect water accumulation in the goaf using the sensitivity of electrical resistivity to water.
The length of survey line N35 is 910 m, and it is located on the 1320 platform. There are two small-scale low-resistance anomalies with relatively small ranges and amplitudes, which are located near point 310 and point 570 (at the red circle of the Figure 14). Among them, the low-resistivity anomaly of measuring point 570 is also shown near the corresponding position of line N0, and the anomaly is caused by the enhanced water yield of the formation near the fault. The abnormal position of point 310 is consistent with the position of the collapse column interpreted by the seismic attribute, which is speculated to be caused by the increased water abundances of the collapse column.
The length of survey line N0 is 790 m, which is located on the 1350 platform. The apparent resistivity profile east of survey point 530 (in the direction of the large point) shows a large low-resistivity anomaly (at the red circle of Figure 15). Combined with the comprehensive interpretation of adjacent survey line N35 and three-dimensional seismic data, it is speculated that the anomaly here is caused by the fragmentation of the formation near the fault and the enhancement in water abundances.
Integrating the results of the seismic interpretation of the goaf and the transient electromagnetic interpretation of the low-resistance abnormal area of each survey line, the scope of the goaf and water accumulation, collapse column, and fault water-rich anomaly displayed in the detection area is delineated on the plane. The subsequent drilling and testing results are basically consistent with the prediction results; in particular, the accurate prediction of the scope of the goaf plays a vital role in safe mining and personnel safety. Figure 16 shows the verification hole DLT-18-6 passing through goaf No. 4. The verification hole DLT-18-6 in the goaf with number 4, compared with the information of the adjacent well, revealed that the section where no coal core was found happened to be the depth section of coal 4 and 9. Due to the presence of goaf, the diameter value of the coal-seam section shows significant fluctuations, and the gamma–gamma value is significantly lower than that of the No.3 coal seam. Combined with the prediction of water accumulation, relevant coal-mining measures can be effectively formulated to avoid production risks (Figure 17).
The local part of the study area is particularly affected by the on-site construction conditions, such as the mining party’s heavy vehicles, gas stations, drilling areas, and blasting areas, which are close to the geophysical exploration work area or even overlap and have a large interference with the seismic signal. In addition, the local coverage times are low, and the signal-to-noise ratio of the data is poor, which directly leads to the poor reliability of the interpretation results. For example, the southwest corner of the study area is a typical disturbed area, with a large area of abnormal areas. However, we predicted the possible results in advance, and the results were confirmed by a very small amount of drilling. Therefore, the judgment of the reliability of data is a very important link.
To ensure safe and efficient mining, it is recommended to deploy ground drilling or advance drilling and geophysical exploration in the abnormal area and poorly controlled area in advance during the mining process to further determine the nature and location of the abnormality and ensure safety [33]. Combined with relevant geological data such as tunneling and drilling in different stages, comprehensive geological research should be strengthened to further improve the reliability and accuracy of interpretation results.

5. Conclusions

For the shallow-buried, deep, and large-gap open-pit mine goaf and water accumulation prediction problem, we adopt an integrated geophysical prospecting strategy. To solve the problems of uneven coverage and noise caused by ultrashallow layers and high drops, we propose a direct bin homogenization and multidomain joint denoising method. The results show that the shallow three-dimensional seismic method can clearly reflect the reflected-wave variation features of the seam. The features of a junction between a normal seam and goaf vary obviously, the effect of the set goaf boundary and scope is good, and the resolution is high. With a detailed analysis of the survey-line repeated area of the shallow three-dimensional seismic and transient electromagnetic methods, along with the apparent resistivity threshold scope of the transient electromagnetic method in the normal seam, goaf and goaf water-accumulation area are determined and can be applied to guide the explanation work. The transient electromagnetic method could especially reflect sensitivity to the low resistance, judge the water accumulation in the goaf, and effectively avoid the limit of the shallow three-dimensional seismic method to judge the water accumulation in the goaf. These results provided useful references for predicting goaf and water accumulation of open-pit mines with ultrashallow and high-drop conditions.

Author Contributions

Conceptualization, S.Z.; methodology, S.Z. and W.G.; validation, S.C. and Q.C.; investigation, Q.C. and Q.M.; writing—original draft preparation, S.Z. and Q.M.; writing—review and editing, S.Z. and S.C.; visualization, W.G. and Y.D.; supervision, S.C.; project administration, S.C. and Q.M.; funding acquisition, S.Z. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (42102212), the unveiling projects of the Department of Science and Technology of Shanxi Province, grant number (20191101018); and the Research and Development Project of Yangquan Coal Industry (Group) Co., Ltd. (GY18027), the National Natural Science Foundation of China (U21A20107).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gui, H.; Qiu, H.; Qiu, W.; Tong, S.; Zhang, H. Overview of goaf water hazards control in China coalmines. Arab. J. Geosci. 2018, 11, 49. [Google Scholar] [CrossRef]
  2. Zhang, J.; Fu, M.; Geng, Y.; Tao, J. Energy saving and emission reduction: A project of coal-resource integration in Shanxi Province, China. Energy Policy 2011, 39, 3029–3032. [Google Scholar] [CrossRef]
  3. He, H.; Yu, W.; Wang, J.; Wen, K. Application research of high density resistivity method to karst geological exploration in mining area. In Proceedings of the 3rd International Conference on Environmental & Engineering Geophysics, Wuhan, China, 15–18 June 2008; pp. 385–388. [Google Scholar]
  4. Hu, G.; Su, J.; Zhao, N. Application of three dimensional laser scanning technology in Pingshuo open-pit mine goaf detection. Opencast Min. Technol. 2015, 7, 61–63, (In Chinese with English abstract). [Google Scholar]
  5. Yang, D.; Guo, W.; Tan, Y. Application of magnetotelluric method to the detection of overburden failure height in shallow seam mining. Arab. J. Geosci. 2018, 11, 350. [Google Scholar] [CrossRef]
  6. Lu, T.; Liu, S.D.; Wang, B.; Wu, R.; Hu, X. A review of geophysical exploration technology for mine water disaster in China: Applications and trends. Mine Water Environ. 2017, 36, 331–340. [Google Scholar] [CrossRef]
  7. Chang, J.; Yu, J.; Liu, Z. Three-dimensional numerical modeling of full-space transient electromagnetic responses of water in goaf. Appl. Geophys. 2016, 13, 539–552. [Google Scholar] [CrossRef]
  8. Su, B.; Yu, J.; Sheng, C. On the borehole electromagnetic method for exploration of coal mining goaf. Elektron. Elektrotech. 2016, 22, 37–40. [Google Scholar] [CrossRef]
  9. Yan, S.; Xue, G.; Qiu, W.; Li, H.; Zhong, H. Feasibility of central loop TEM method for prospecting multilayer water-filled goaf. Appl. Geophys. 2016, 13, 587–597. [Google Scholar] [CrossRef]
  10. Dam, J.; Meulenkamp, J. Some results of the geo-electrical resistivity method in ground water investigations in The Netherlands. Geophys. Prospect. 2010, 15, 92–115. [Google Scholar] [CrossRef]
  11. Liang, K.; Zhu, Y.; Shen, L. Application of high density resistivity method in groundwater exploration in Weifang. Ground Water 2013, 35, 112–113. [Google Scholar]
  12. Ueda, T.; Mitsuhata, Y.; Uchida, T.; Marui, A.; Ohsawa, K. A new marine magnetotelluric measurement system in a shallow-water environment for hydrogeological study. J. Appl. Geophys. 2014, 100, 23–31. [Google Scholar] [CrossRef]
  13. Wang, N.; Wang, Z.; Sun, Q.; Hui, J. Coal Mine Goaf Interpretation: Survey, Passive Electromagnetic Methods and Case Study. Minerals 2023, 13, 422. [Google Scholar] [CrossRef]
  14. Francke, J.; Utsi, V.; Singh, K.K.K. The design of an intrinsically safe ground penetrating radar for detecting abandoned workings in underground coal mines. In Proceedings of the IEEE 14th International Conference on Ground Penetrating Radar (GPR), Shanghai, China, 4–8 June 2012; pp. 125–130. [Google Scholar]
  15. Kotyrba, A.; Stańczyk, K. Sensing underground coal gasification by ground penetrating radar. Acta Geophys. 2017, 65, 1185–1196. [Google Scholar] [CrossRef]
  16. Brocher, T.; Ewing, J. A comparison of high-resolution seismic methods for determining seabed velocities in shallow water. J. Acoust. Soc. Am. 1998, 79, 286–298. [Google Scholar] [CrossRef]
  17. Wilkinson, P.; Chambers, J.; Meldrum, P.; Ogilvy, R.; Mellor, C.; Caunt, S. A comparison of self-potential tomography with electrical resistivity tomography for the detection of abandoned mineshafts. J. Environ. Eng. Geophys. 2005, 10, 381–389. [Google Scholar] [CrossRef]
  18. Li, S.; Zhou, Z.; Ye, Z.; Li, L.; Zhang, Q.; Xu, Z. Comprehensive geophysical prediction and treatment measures of karst caves in deep buried tunnel. J. Appl. Geophys. 2015, 116, 247–257. [Google Scholar] [CrossRef]
  19. Zhang, S.; Huang, H.; Dong, Y.; Yang, X.; Wang, C.; Luo, Y. Direct estimation of the fluid properties and brittleness via elastic impedance inversion for predicting sweet spots and the fracturing area in the unconventional reservoir. J. Nat. Gas Sci. Eng. 2017, 45, 415–427. [Google Scholar] [CrossRef]
  20. Lelièvre, P.; Oldenburg, D. A comprehensive study of including structural orientation information in geophysical inversions. Geophys. J. Int. 2010, 178, 623–637. [Google Scholar] [CrossRef]
  21. Liu, B.; Liu, Z.; Li, S.; Nie, L.; Su, M.; Sun, H. Comprehensive surface geophysical investigation of karst caves ahead of the tunnel face: A case study in the xiaoheyan section of the water supply project from songhua river, Jilin, China. J. Appl. Geophys. 2017, 144, 37–49. [Google Scholar]
  22. Chang, S.; Qu, Z.; Zhang, S.; Fu, Y.; Chen, Q. Suppression of multiple reflection–refraction waves in the shallow low-velocity zone via vertical source combination. Acta Geophys. 2021, 69, 1799–1811. [Google Scholar] [CrossRef]
  23. Zhang, S.; Huang, H.; Zhu, B.; Li, H.; Zhang, L. Seismic facies-controlled prestack simultaneous inversion of elastic and petrophysical parameters for favourable reservoir prediction. Explor. Geophys. 2018, 49, 655–668. [Google Scholar] [CrossRef]
  24. Zhang, X.; Ma, J.; Li, L. Monitoring of coal-mine goaf based on 4D seismic technology. Appl. Geophys. 2020, 17, 54–66. [Google Scholar] [CrossRef]
  25. Vakman, D. On the analytic signal, the teager-kaiser energy algorithm, and other methods for defining amplitude and frequency. IEEE Trans. Signal Process. 2002, 44, 791–797. [Google Scholar] [CrossRef]
  26. Taner, M.; Koehler, F.; Sheriff, R. Complex seismic trace analysis. Geophysics 1979, 44, 1041–1063. [Google Scholar] [CrossRef]
  27. Gersztenkorn, A.; Marfurt, K.J. Eigenstructure-based coherence computations as an aid to 3-d structural and stratigraphic mapping. Geophysics 2012, 64, 1468–1479. [Google Scholar] [CrossRef]
  28. Fitterman, D. Transient electromagnetic sounding for groundwater. Geophysics 1986, 51, 995–1005. [Google Scholar] [CrossRef]
  29. Shi, R.; Darcherif, A.; Sabonnadiere, J. Computation of transient electromagnetic fields radiated by a transmission line: An exact model. IEEE Trans. Magn. 1995, 31, 2423–2431. [Google Scholar] [CrossRef]
  30. Chen, J.; Pi, S.; Zhang, Y.; Lin, T. Weak coupling technology with noncoplanar bucking coil in a small-loop transient electromagnetic system. IEEE Trans. Ind. Electron. 2021, 69, 3151–3160. [Google Scholar] [CrossRef]
  31. Qi, T.; Zhang, F.; Pei, X.; Feng, G.; Wei, H. Simulation research and application on response characteristics of detecting waterflled goaf by transient electromagnetic method. Int. J. Coal Sci. Technol. 2022, 9, 17–26. [Google Scholar] [CrossRef]
  32. Xue, G.; Pan, D.; Yu, J. Review the applications of geophysical methods for mapping coal-mine voids. Prog. Geophys. 2018, 33, 2187–2192. (In Chinese) [Google Scholar]
  33. Wang, B.; Li, W.; Ye, Q.; Ma, K. A nonparametric spectrum estimation method for dispersion and attenuation analysis of borehole acoustic measurements. Pet. Sci. 2022, 20, 241–248. [Google Scholar] [CrossRef]
Figure 1. Topographic map of the Pingshuo Dong open-cast coal mine.
Figure 1. Topographic map of the Pingshuo Dong open-cast coal mine.
Processes 11 01653 g001
Figure 2. Topographical photos of the study area, (a) high drop, (b) narrow platform, (c) high voltage line, (d) complex noise environment.
Figure 2. Topographical photos of the study area, (a) high drop, (b) narrow platform, (c) high voltage line, (d) complex noise environment.
Processes 11 01653 g002
Figure 3. Rectangular loop lines and equivalent eddy-current ring, (a) establishment of stable magnetic field; (b) attenuation of magnetic field (Tx: transmission line frame. t: decay time).
Figure 3. Rectangular loop lines and equivalent eddy-current ring, (a) establishment of stable magnetic field; (b) attenuation of magnetic field (Tx: transmission line frame. t: decay time).
Processes 11 01653 g003
Figure 4. Layout of 3D seismic and transient electrical detection.
Figure 4. Layout of 3D seismic and transient electrical detection.
Processes 11 01653 g004
Figure 5. The distribution of shot and receiver points (a), the actual coverage times (b), and coverage times after bin homogenization (c) in the study area.
Figure 5. The distribution of shot and receiver points (a), the actual coverage times (b), and coverage times after bin homogenization (c) in the study area.
Processes 11 01653 g005
Figure 6. Comparison of stacked gathers before and after bin homogenization, (a) before bin homogenization, (b) after bin homogenization, (c) before bin homogenization, (d) after bin homogenization.
Figure 6. Comparison of stacked gathers before and after bin homogenization, (a) before bin homogenization, (b) after bin homogenization, (c) before bin homogenization, (d) after bin homogenization.
Processes 11 01653 g006
Figure 7. Comparison of stacked gathers before and after bin homogenization.
Figure 7. Comparison of stacked gathers before and after bin homogenization.
Processes 11 01653 g007
Figure 8. Comparison of single-shot records before (left) and after (right) noise suppression.
Figure 8. Comparison of single-shot records before (left) and after (right) noise suppression.
Processes 11 01653 g008
Figure 9. Comparison of stacked gathers before (upper) and after (lower) suppressing noise.
Figure 9. Comparison of stacked gathers before (upper) and after (lower) suppressing noise.
Processes 11 01653 g009
Figure 10. Forward modeling and seismic profile of goaf of different types and in different positions.
Figure 10. Forward modeling and seismic profile of goaf of different types and in different positions.
Processes 11 01653 g010
Figure 11. Seismic response of goaf in actual seismic profiles. (a) Increased energy and low frequency, and (b) weak energy and worse continuity.
Figure 11. Seismic response of goaf in actual seismic profiles. (a) Increased energy and low frequency, and (b) weak energy and worse continuity.
Processes 11 01653 g011
Figure 12. Distribution plan of seismic-interpretation abnormal areas of coal seams 4 (left) and 9 (right), ①–④ are the number of the goafs.
Figure 12. Distribution plan of seismic-interpretation abnormal areas of coal seams 4 (left) and 9 (right), ①–④ are the number of the goafs.
Processes 11 01653 g012
Figure 13. Comparison of logging and drilling properties of collapsed columns.
Figure 13. Comparison of logging and drilling properties of collapsed columns.
Processes 11 01653 g013
Figure 14. One-dimensional inversion profile of the apparent resistivity of the N35 line.
Figure 14. One-dimensional inversion profile of the apparent resistivity of the N35 line.
Processes 11 01653 g014
Figure 15. One-dimensional inversion profile of the apparent resistivity of the N35 line.
Figure 15. One-dimensional inversion profile of the apparent resistivity of the N35 line.
Processes 11 01653 g015
Figure 16. Borehole Histogram of the DLT-18-6 (The red line is GGFR, the blue line is CAL, and the green line is CAL).
Figure 16. Borehole Histogram of the DLT-18-6 (The red line is GGFR, the blue line is CAL, and the green line is CAL).
Processes 11 01653 g016
Figure 17. Plan distribution map for predicting risk areas of goaf and water accumulation.
Figure 17. Plan distribution map for predicting risk areas of goaf and water accumulation.
Processes 11 01653 g017
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, S.; Guan, W.; Chang, S.; Meng, Q.; Dong, Y.; Chen, Q. Integrated Geophysical Prediction of Goaf and Water Accumulation in Pingshuo Dong Open-Cut Mine with Ultrashallow and High Drops. Processes 2023, 11, 1653. https://doi.org/10.3390/pr11061653

AMA Style

Zhang S, Guan W, Chang S, Meng Q, Dong Y, Chen Q. Integrated Geophysical Prediction of Goaf and Water Accumulation in Pingshuo Dong Open-Cut Mine with Ultrashallow and High Drops. Processes. 2023; 11(6):1653. https://doi.org/10.3390/pr11061653

Chicago/Turabian Style

Zhang, Sheng, Wenyuan Guan, Suoliang Chang, Qinggang Meng, Yinping Dong, and Qiang Chen. 2023. "Integrated Geophysical Prediction of Goaf and Water Accumulation in Pingshuo Dong Open-Cut Mine with Ultrashallow and High Drops" Processes 11, no. 6: 1653. https://doi.org/10.3390/pr11061653

APA Style

Zhang, S., Guan, W., Chang, S., Meng, Q., Dong, Y., & Chen, Q. (2023). Integrated Geophysical Prediction of Goaf and Water Accumulation in Pingshuo Dong Open-Cut Mine with Ultrashallow and High Drops. Processes, 11(6), 1653. https://doi.org/10.3390/pr11061653

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop