Sage Journals: Discover world-class research

Abstract

Once a coal seam has been weathered and oxidized, the quality of the coal is lost and the coal resources are wasted. At the same time, new risks are introduced to the safe operation of the coal mine. Therefore, it is critical to determine the extent of the weathered and oxidized zone in the coal seam. This paper depicts the electromagnetic characteristics of the weathered and oxidized zone of coal seams, briefly introduces the method of detection by transient electromagnetic method (TEM) and the data processing steps, and performs the numerical simulation for the simplified resistivity model of weathered and oxidized zone in a coal seam. It is observed that after the coal seam has been weathered and oxidized, the electromagnetic characteristics (secondary field potential and resistivity) will be dramatically affected, namely high secondary field potential and low resistivity. For experimental verification purposes, this study selects the test project of the Shixin Minefield in Shanxi, China as an example, and the accuracy for the detection of the weathered and oxidized zone of the coal seam using TEM is successfully verified.

Keywords

Transient electromagnetic method coal seams weathered and oxidized zone geophysical exploration electrical characteristics

Introduction

Due to crustal movement, coal seams buried deep in the crust are exposed to the surface or buried in the shallow crust A variety of intrinsic physical and chemical features of coal seams may be modified as a result of the interplay of temperature variations, gas, water, and biological activity, in turn producing a weathered zone, and this process is known as weathering. In the weathering process, in addition to the dissolution of water, there is also oxidation, thus forming an oxidized zone. A weathered zone and an oxidized zone are often collectively referred to as a weathered and oxidized zone. The weathered and oxidized zone of a coal seam generally develops in an outcrop or shallow buried area, and the depth of the weathered and oxidized zone differs due to the influence of various factors such as topography, rock layer properties, structure, and hydrogeological conditions (Cui et al., 2015; Yang, 2019; Zhan and Wang, 2016). The weathered and oxidized effect of coal seam greatly affects the effective mining of coal seams, thereby causing direct economic losses and waste of resources. It can even result in spontaneous combustion of coal seams, gas and coal dust explosions, smoke poisoning, or casualties (Wang and Liu, 2014).

At present, methods for detecting weathered and oxidized zones in coal seams include sonic logging, roadway high-density electrical approaches, roadway seismic investigation, and drilling. The time difference between the sound wave traveling from the sending probe through the formation to the receiving probe is recorded using sonic logging. The change in the time difference determines the degree of coal seam compactness, and then weathered and oxidized zone of the coal seam can be detected. This method is prone to being influenced by several factors: when the diameter of the well changes, the sonic time difference in the expansion section is reduced, thereby causing false anomalies in the time difference curve; when the layer thickness changes, then the acoustic velocity logging tool has poor resolution for thin formations which are smaller than the spacing; finally, if the stratum is loose or cracks are developed, and the stratum absorbs excessive sound wave energy, this causes the cycle jump to appear as an irregular jump on the sound velocity curve (You and Huang, 2015). The roadway high-density electrical technique measures the electrical differences between the weathered and oxidized zone and the coal seam. Strict building criteria must be adhered to when applying this approach in coal mine roadways (Huang, 2014; Zhang, 2008). Underground seismic exploration is generated and received in the coal-measure strata, absent from the influence of loose surface layer absorption or the shielding effect of the coal-overlying measure's high-speed layer. In addition, the complex working environment in the roadway, along with the irregularity of the distribution of weathered and oxidized zones, further increase the difficulty involved in this method. For example, when the thickness of the coal seam changes greatly and the roof and floor of the coal seam are undulating, then intermittent and short-axis reflection wave groups generally exist on the seismic profile. The continuity of the wave groups is often difficult to track, and stray reflections appear (Huang, 2016; Shan and Li, 2014; Zhan et al., 2015). Finally, as a direct detection method, drilling can effectively reveal the composition, depth, and degree of a weathered and oxidized zone, but it is costly to use, the network is sparse, and it is difficult to accurately locate the small area of weathered and oxidized zone. Therefore, this method is not suitable for large-scale implementation (Cheng and Shi, 2013).

The transient electromagnetic method (TEM) is an effective method that can be used to study the spatial or temporal distribution of electromagnetic fields based on the differences in electrical properties such as the conductivity of rocks or ore bodies in the Earth's crust (Cao et al., 2019; Chen, 2017; Fu et al., 2019; Guo et al., 2020; Lai et al., 2016; Li, 2019), so as to solve geological problems. TEM has the advantages of flexible field source, diverse methods, stability and high efficiency. It has been developed as an important method for detecting oil and gas (Anton et al., 2010; Cairns et al., 1996; Meyer et al., 2016), metal and non-metal minerals (Asten and Duncan, 2012; Cai et al., 2017; Chen et al., 2010; Gao et al., 2016; Qin et al., 2019; Xue et al., 2014; Xue et al., 2020; Xue and Yu, 2017; Yu et al., 2007; Zhi et al., 2016). TEM has been widely applied in the fields of geological structure research (Hisatoshi et al., 2011; Wu et al., 2020), and geothermal exploration (Cheng and Gao, 2011; Di et al., 2019; Tohru et al., 1998; Zhou et al., 2018). However, at present few researchers have attempted to apply TEM to the detection of coal seam weathered and oxidized zone.

In this paper, to verify the effectiveness of TEM in detecting the weathered and oxidized zone of a coal seam, we analyze the geophysical characteristics of weathered and oxidized zone of coal seam, propose implementation steps of data processing, and conduct numerical simulation for a simplified model, taking the project of the Shixin Coal Mine in Shanxi Province, China as an example, and satisfactory detection results are achieved.

Methodology

Geophysical characteristics of a coal seam weathered and oxidized zone

The resistivity of coal is quite high; however, as the degree of weathering and oxidation of coal increases, its resistivity and density exhibit a decreasing trend, while the natural gamma amplitude shows an increasing trend. In weathered coal, due to the increase in the proportion of shale content, the main characteristics are low resistivity, low density, and high natural gamma amplitude. Meanwhile, the main characteristics of oxidized coal show that the resistivity and density are slightly higher than those of weathered coal, which are between normal coal and weathered coal, and the natural gamma amplitude is also between that of normal and weathered coal (Sun and Hu, 2018).

In general, coal seams bear the characteristics of layered distribution, while their electrical characteristics are relatively uniform in the lateral direction, and gradually change according to a certain rule in the vertical direction. When the coal seam is dense and integrated, then the resistivity of the coal seam is relatively high (about 1800–3000 Ω·m, with an average of 2400 Ω·m). However, as the degree of weathering and oxidation of the coal seam increases, its resistivity also exhibits a downward trend. It can be inferred that the deeper the degrees of weathering and oxidation are, the lower the apparent resistivity will be. As an example, the resistivity of oxidized coal is about 1000–2000 Ω·m, with an average of 1500 Ω·m, which is roughly 40% lower than that of normal coal. In addition, the resistivity of weathered coal is about 0–1000 Ω·m, with an average of 500 Ω·m, which is about 80% lower than that of normal coal (Figure 1). This is the geophysical premise for coal seam weathered and oxidized zone detection using the electromagnetic method.

Figure 1.

Comparison of geophysical logging characteristics of (a) weathered coal, (b) oxidized coal and (c) normal coal (cited from Sun and Hu, 2018).

TEM methodology

Method selection basis

As mentioned above, the detection methods of the weathered and oxidized zone of coal seams mainly include sonic logging, roadway seismic exploration, roadway high-density electrical method, and drilling. These methods can indicate the characteristics and range of the coal seam weathered and oxidized zone from some aspects, but they all have certain limitations, and the detection accuracy requires improvement. In general, the resistivity of the weathered and oxidized zone is lower than that of the coal seam, and TEM bears a more obvious detection effect on the high-resistivity surrounding rock (since the coal seam is regarded as the high-resistivity surrounding rock), thus it can be used to detect the range of the coal seam weathered and oxidized zone.

Data processing

Based on the processed data, the slice map of the apparent resistivity for the target layer and the contour map of the secondary field potential of a certain channel can be drawn, and the data can be explained according to these two maps. The following method can be used for data processing (Bai et al., 2003; Chen and Xu, 2017; Di et al., 2016; Efthimios and Michael, 2000):

Transmitting the TEM data to the computer, and under the excitation of the step pulse current, the transient response of the ground half-space (vertical component of magnetic induction $B_{z} (t)$ and its rate of changing over time $\frac{d B_{Z} (t)}{d t}$ can be expressed as follows:

B_{Z} (t) = \frac{μ m}{4 π r^{3}} [(\frac{9}{2 Z^{2}} - 1) e r f (Z) - \frac{1}{\sqrt{π}} (\frac{9}{Z} + 4 Z) e^{- Z^{2}}]

(1)

\frac{d B_{Z} (t)}{d t} = - \frac{μ m}{2 π r^{5}} [- 9 e r f (z) + \frac{2 Z}{\sqrt{π}} (9 + 6 Z^{2} + 4 Z^{4}) e^{- Z^{2}}]

(2)

Z = \frac{r}{2} \sqrt{\frac{u}{ρ t}}

(3)

Where u is the permeability of a uniform conductive half space (approximately taken as 4π × 10⁻⁷ h/m); m is the emission magnetic moment m = NIS; N is the number of turns of the emission loop; I is the emission current; S is the area of the loop; r is the receiving-sending distance; ρ is the resistivity of the uniform earth medium; t is the transient delay; and

e r f (z)

is the error function, the expression of which is

e r f (Z) = \frac{2}{\sqrt{π}} \int_{0}^{Z} e^{- t^{2}} d t

The resistivity can be obtained from the transformation of Formula (3), as follows:

ρ = \frac{r^{2} μ}{4 t} \frac{1}{Z^{2}}

(4)

Next, the following definitions are made:

F_{B} (Z) = [(\frac{9}{2 Z^{2}} - 1) e r f (Z) - \frac{1}{\sqrt{π}} (\frac{9}{Z} + 4 Z) e^{- Z^{2}}]

(5)

F_{D} (Z) = \frac{1}{Z^{2}} [- 9 e r f (Z) + \frac{2 Z}{\sqrt{π}} (9 + 6 Z^{2} + 4 Z^{4}) e^{- Z^{2}}]

(6)

Then, we can obtain the following Formulas from Formulas (1) and (2):

B_{Z} (t) = \frac{μ m}{4 π r^{3}} F_{B} (Z)

(7)

\frac{d B_{Z} (t)}{d t} = - \frac{μ m}{8 π t r^{3}} F_{D} (Z)

(8)

F_{B} (Z)

is the kernel function of the vertical component of the magnetic induction intensity

B_{Z} (t)

; and

F_{D} (Z)

is the kernel function of the rate of changing over time

\frac{d B_{Z} (t)}{d t}

of the vertical component of the magnetic induction intensity

B_{Z} (t)

Using Formula (8), the secondary field potential at different time channel can be calculated. Then, by cutting the secondary field potential and the resistivity along the target coal seam into slices, the contour map of the secondary field potential and the slice map of the resistivity of the target coal seam can be obtained.

Numerical modeling

In order to verify the ability of TEM to identify weathered and oxidized zone of coal seam, numerical simulation analysis is carried out. The numerical simulation program is the open-source code SIMPEG written by the Geophysical Group of the University of British Columbia, and the forward method is the finite volume method, which is presented as follows.

Electromagnetic method based on the finite volume method

The Maxwell equations satisfied by the electromagnetic field in the time domain can be expressed as follows:

\nabla \times E_{s} + μ \frac{\partial H_{s}}{\partial t} = 0

(9)

\nabla \times H_{s} - σ E_{s} - (σ - σ_{b}) E_{b} - ε \frac{\partial E_{s}}{\partial t} = 0

(10)

Then, taking the curl on both sides of the Formula (9) at the same time, we can obtain:

\nabla \times \nabla \times E_{s} + μ \frac{\partial \nabla \times H_{s}}{\partial t} = 0

(11)

Substituting Formula (10) into Formula (11) and ignoring the displacement current, we can obtain:

\nabla \times \nabla \times E_{s} + μ σ \frac{\partial E_{s}}{\partial t} + μ (σ - σ_{b}) \frac{\partial E_{b}}{\partial t} = 0

(12)

Formula (12) is the double curl equation of the secondary electric field generated by the abnormal body. The finite volume method and backward Euler implicit scheme method can be used to discretize the equations, so as to achieve the numerical solution of the equation.

We can divide the calculation area into a series of non-repetitive control volumes in the finite volume method, so that there is a control volume around each calculation point, then integrate each control volume with the differential equation which is to be solved. In this paper, we use Yee's staggered grid (Yee, 1996) for spatial discretization. Taking the E y component of the electric field strength E at any point as an example, the control volume can be explained as follows (Figure 2):

Figure 2.

Control volume of E in the Yee grid.

Integrating each variable within its control volume, and then integrating the control volume V with the anomalous field diffusion equation in any calculation point, we can obtain:

\int \int \int_{V} \nabla \times \nabla \times E_{s} d V + \int \int \int_{V} μ σ \frac{\partial E_{s}}{\partial t} d V + \int \int \int_{V} μ (σ - σ_{b}) \frac{\partial E_{b}}{\partial t} d V = 0

(13)

Formula (13) is the electromagnetic integral equation of the time-domain finite volume method, which regards the secondary electric field generated by the three-dimensional geological body as being unknown. Expressing Formula (16) in discrete form, we can obtain:

G^{V} \frac{\partial E_{s}}{\partial t} + G^{\partial V} E_{s} = - G^{V r} \frac{\partial E_{b}}{\partial t}

(14)

Where

G^{V}

and

G^{\partial V}

are the discrete operator of the volume integral and the area integral of the secondary electric field containing conductivity parameters, respectively. In addition,

G^{V r}

is a discrete operator of the volume integral of the background electric field containing the conductivity kernel function.

The implicit Euler method can then be used, expressed as follows:

f^{'} (t) = \frac{f (t) - f (t - Δ t)}{Δ t}

(15)

According to Formula (15), the derivative term in Formula (14) can be discretized as:

\frac{\partial E_{s}}{\partial t} = \frac{1}{Δ t} (E_{s}^{n} - E_{s}^{n - 1})

(16)

Where n represents the current moment and n-1 represents the previous moment, respectively. After discretizing Formula (14), we can obtain the following:

{C E}_{s}^{n} + {D E}_{s}^{n - 1} = R (E_{b}^{n} - E_{b}^{n - 1})

(17)

Where C and D are the secondary electric field coefficient matrices at any point in space at different times, both of which are sparse symmetric matrices, and R is the sparse symmetric matrix composed of related parameters, such as background conductivity.

Numerical simulation result

Based on the above theoretical analysis, we design a simplified resistivity model and perform numerical simulation (Figure 3). The basic parameters of the TEM numerical calculation are as follows: the transmitting source is a single-turn square loop (200 × 200 m); the voltage responses at different moments are received in the center of the receiving loop; the transmitting waveform is a bipolar rectangular pulse; the time range is 0.48 to 28.7 ms; the number of time channels is 22; and the distance of the receiving point is 20 m.

Figure 3.

Simplified resistivity model of the geothermal system with the coal seam as the target body.

It can be seen that in the simulated resistivity profile, normal coal, weathered coal, and oxidized coal are all clearly reflected, and there are significant differences among them. Both the weathered and oxidized zones of the coal seam exhibit anomalous characteristics of low resistivity. At the distance of 100–200 m with the depth of about 130 m, the resistivity is as low as 40–80 (Ω·m), reflecting the existence of weathered coal; at the distance of 300–400 m, the resistivity is also lower than that of the normal coal with an range of about 80–120 (Ω·m), which is the manifestation of the coal seam being oxidized. Based on the numerical simulation results, the effect of TEM in detecting the weathered and oxidized zone of coal seam has been well proven theoretically.

Example analysis

Transient electromagnetic exploration was carried out at the Shixin Coal Mine, so as to assess the accuracy of the transient electromagnetic method in located the coal seam weathered and oxidized zone.

Geological setting

The Shixin Coal Mine is located in Gujiao City, Shanxi Province, China. The geographical coordinates are as follows: east longitude 112°15′02″–112°17′03″; north latitude 37°52′20″–37°53′40″ (Figure 4).

Figure 4.

Numerical simulation results of the simplified resistivity model.

Figure 5.

Location of the Shixin Coal Mine.

The strata in the coal mine, from old to new, are as follows: Middle Ordovician Fengfeng Formation (O₂f), Benxi Formation of Middle Carboniferous System (C₂b), Taiyuan Formation of Upper Carboniferous System (C₃t), Lower Permian Shanxi Formation (P₁s), Lower Permian Lower Shihezi Formation (P₁x), Upper Permian Upper Shihezi Formation (P₂s), and fourth series (Q).

Among them, the Lower Permian Shanxi Formation (P₁s) is the main coal-bearing formation. This group is continuously deposited on the underlying Taiyuan Formation. It is a set of continental coal-bearing strata. The thickness of the whole group is 43.86–64.74 m, with an average of 52.37 m. The lithology is mainly dark gray and gray-black mudstone, siltstone, and sandy mudstone interposed with gray-white sandstone, containing five layers of coal, of which the numbering from top to bottom is No. 02, 03, 1, 2, 4, among which the No. 2 coal seam is the main mineable coal seam. The No. 2 coal seam is located in the upper part of the Shanxi Formation, which is relatively stable, with a thickness of 0.75–3.00 m, averaging 2.08 m. It has a simple structure, and occasionally contains one layer of gangue. The roof of the No. 2 coal seam consists of fine-grained sandstone, partially sandy mudstone, and siltstone, while the bottom is gray-white, medium-grained quartz sandstone, and partially coarse-grained sandstone. The No. 2 coal seam is a stable mineable coal seam.

According to previous geological data, the weathered and oxidized zone the No. 2 coal seam of the Shixin Coal Mine is mainly located in the eastern and northern parts of the minefield. During the mining of the No. 2 coal seam, it was found that its roof was missing in some areas, and the roof fell during the coal cutting, making it difficult to provide support. The coal mine staff drilled the roof of the No. 2 coal seam at point A in transportation roadway 284 (Figure 6).

Figure 6.

Location of transportation roadway 284.

When drilling 2.8 m into the roof of the No. 2 coal seam, phenomena such as dripping water and softening of the rock occurred. It was thus determined that there was a weathered and oxidized zone on the roof of the No. 2 coal seam (Figure 7).

Figure 7.

Roof of transportation roadway 284.

However, the reliability of inferring the range of coal seam weathered and oxidized zone only by geological and exploration data is far from sufficient, and is of little significance for the safety production of coal mines. The scope of coal seam weathered and oxidized zone requires further exploration by geophysical methods. The exploration area is the entire minefield of the Shixin Coal Mine, with a total area of 6.3 km², and there are 61 survey lines laid out. The dimensions of the working network are 40 × 20 m, i.e. the line distance is 40 m and the point distance is 20 m (Figure 6).

Method validity test

In the data collection work, the GDP-32II multi-function workstation of ZONGE Company was used; the configuration of TEM was the center loop device, and data collecting area was within 1/9 of the area inside the loop (i.e. the side length of the data collecting area was 1/3 of the side length of the transmitting loop). The transmitting parameters were as follows: the transmitting loop size was 160 × 160 m, the transmitting current was 8 A, and the number of turns of the transmitting loop was 1. Receiving parameters: the acquisition probe matched with GDP-32II was used to collect the data, which had an equivalent acquisition area of 10,000 m², and the data sampling delay time was 220 μs.

The method validity test location is selected on line 137 (the green line in Figure 7). According to the underground mining data, the eastern part of this line is the junction of the coal seam weathered and oxidized zone and the normal coal seam. The team decided to test here to distinguish the electrical difference between the weathered and oxidized zone and the normal coal seam.

According to steps described in the previous part “Data Processing,” the fifth surveying channel is the target channel. Figure 8 shows the multi-channel potential profile and resistivity profile of the testing line. It can be seen that there is a high potential abnormality present in the 1360–1600 section in the multi-channel potential profile, and the potential value is greater than 10,060 μV/A (green grid area). This value is regarded as the abnormal potential threshold value, which is approximately 50% higher than the normal section potential value (about 7000 μV/A). In the resistivity profile, there are obvious low-resistivity characteristics (gray grid area) in “Introduction”, 360–1600, and the resistivity value is less than 55 Ω·m, which is roughly 45% smaller than the normal resistivity value (80 Ω·m), which is regarded as the abnormal threshold of resistivity. “Introduction”, 360–1600 is inferred as a coal seam weathered and oxidized zone, which is in good agreement with the results (blue grid grid) revealed by the underground mining data, thus indicating that the use of TEM to detect the scope of the coal seam weathered and oxidized zone is effective.

Figure 8.

Multi-channel potential profile (a) and resistivity profile of the testing line (b).

Application Effect

According to the previous steps, the contour map of the secondary field potential of the fifth surveying channel and the slice map of the resistivity of the No. 2 coal seam are drawn, as shown in Figure 9. It can be concluded that the weathered and oxidized zone of the coal seam has two forms of electrical characteristics, i.e. high-potential and low-resistivity.

Figure 9.

Comparison of (a) the secondary field potential contour map of the fifth surveying channel and (b) the resistivity slice map of the No. 2 coal seam.

Figure 9 shows the comparison between the contour map of the secondary field potential of the fifth surveying channel and the slice map of the resistivity of the No. 2 coal seam. The following observations can be made:

The secondary field potential in the exploration area is in the range of (9940–10,130) μV/A. According to the test results, the area where the secondary field potential is greater than the abnormal threshold of 10,600 μV/A is treated as a high-potential abnormal area. Two high-potential anomaly areas are delineated, which are located in the eastern and central parts of the exploration area, and the eastern anomaly area is relatively large (the gray range in Figure 9(a)).

In addition, the resistivity in the exploration area undergoes a large overall change, and the resistivity value is between 15–155 (Ω·m). According to the test results, the area where the resistivity value is less than the abnormal threshold of 50 Ω·m is classified as the low-resistivity abnormal area. A total of two low-resistivity anomaly areas are delineated, which are respectively located on the eastern and western sides of the exploration area, and the eastern anomaly area is relatively large (the gray range in Figure 9(b)).

According to the characteristics of secondary field potential and resistivity, combined with geological data, the eastern anomaly area is inferred as the coal seam weathered and oxidized zone.

Discussion

By the above analysis, the effectiveness of TEM to detect the distribution of weathered and oxidized zone of a coal seam has been well verified, and the detection mechanism can be discussed from the aspects of physics and chemistry (Gayle et al., 1984).

On the one hand, when a coal seam is weathered and oxidized, the physical processes include freezing and thawing, and wetting and drying. During freezing, fractures in the coal are extended and enlarged when the increased volume of ice forces unfrozen water farther into the coal. Physical breakdown from wetting and drying is caused by internal stresses that are set up from unequal volume changes that take place as water is absorbed and desorbed faster on the outside of a coal particle than that on the inside. Thus, wetting and drying create fractures that can be extended and enlarged during freezing and thawing; freezing and thawing force water into new areas where wetting and drying can take place, resulting in coal seams with higher void ratios and moisture content.

The primary chemical cause for the creation of coal seam weathered and oxidized zone, on the other hand, is that the roof of the coal measure strata and cracks around the structure emerge, which have strong water circulation conditions. Surface water transports oxygen and carbon dioxide into the earth. Under the leaching action of groundwater, the coal seams dissolve carbonate rocks and sulfate rocks (CaCO₃, MgCO₃, CaSO₄, MgSO₄, etc.), the water quality is mainly HCO₃-Ca·Mg or SO₄·HCO₃-Ca·Mg, which has stronger electrical conductivity compared to a normal coal seam (Zhao and Fu, 2008).

Based on these physical-chemical characteristics, it can be seen that after the coal seam is weathered and oxidized, the electrical characteristics will alter significantly, which is mainly reflected in the significant decrease in resistivity. According to these dramatically changes of electrical characteristics, the distribution range of the weathered and oxidized zone of a coal seam can be well detected by the transient electromagnetic method.

Conclusions

Through theoretical and engineering case analysis, some important conclusions are summarized from this study: The weathering and oxidation of coal seam will exert a great influence on coal quality and coal seam mining. Therefore, the identification of the weathered and oxidized zone in the coal seam is critical to the safe production of coal mines. The use of logging data to evaluate physical characteristics is critical. It shows that the resistivity value of the weathered and oxidized zone is at least 40% smaller than that of normal coal, which provides guidance for the creation of transient electromagnetic work, and can decrease work blindness. After the coal seam has been oxidized or weathered, it possesses the electrical properties of “high potential value (50% greater than the normal value) and low resistivity (45% smaller than the normal resistivity value.” The employment of a transient electromagnetic method in the Shixin Coal Mine to identify the coal seam weathered and oxidized zone range is possible and effective, thereby providing a more dependable geological basis for the safe production of coal mines.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This Research is sponsored by the National Natural Science Foundation of China (Nos. 41674075, 42062015, and 42174080), co-funded by the Natural Science Foundation of Guangxi Province (Nos. 2016GXNSFGA380004 and 2018AD19204), Research Start-up Foundation of Guilin University of Technology (No. RD2100002165), and China Postdoctoral Science Foundation Funded Project (No. 2021MD703820).

Data availability

The data which support the research are updated with the manuscript.

ORCID iDs

Yang Yang

Siqin Liu

References

Anton

Ronnie

David

, et al. (2010) Multi-transient electromagnetic repeatability experiment over the north sea harding field. Geophysical Prospecting 58: 1159–1176.

Asten

Duncan

(2012) Some consequences and advantages of using B-field sensors in time-domain EM measurement for mineral exploration and unexploded ordnance search. Geophysics 77: WB137–WB148.

Bai

Jian

, et al. (2003) Numerical calculation of all-time apparent resistivity for central loop transient electromagnetic method. Chinese Journal of Geophysics 46(5): 998–1010.

Cai

Xiong

(2017) Finite-element time-domain modelling of electromagnetic data in general dispersive medium using adaptive Padé series. Computers & Geosciences 109: 1945.

Cairns

Evans

Edwards

(1996) A time-domain electromagnetic survey of the TAG hydrothermal mound. Geophysical Research Letters 23(23): 3455–3458.

Cao

Kang

Wei

(2019) Exploration and practice of step-by-step inversion of high-density resistivity method: taking a geological section of Changde Ding Cheng development zone as an example. Progress in Geophysics (In Chinese) 34(6): 2398–2405.

Chen

(2017) The influence of topography on TEM data. Coal Geology & Exploration 45(6): 139–143.

Chen

(2017) Suitable algorithm of transient electromagnetic turn-off effects and whole-stage apparent resistivity. Coal Geology & Exploration 45(4): 131–134.

Chen

Liu

, et al. (2010) Integrated geophysical exploration for the Longtoushan Ag-Pb-Zn deposit in the southeast of the Da Xing’an Ling mountains, Inner Mongolia, northern China. Exploration Geophysics 41(4): 279–288.

10.

Cheng

Shi

(2013) Current status and development of coal geophysical technology in China. Progress in Geophysics 28(4): 2024–2032.

11.

Cheng

Gao

(2011) Application of tem in hydrological geology exploration. Coal Technology 30(7): 124–126.

12.

Cui

Zhu

, et al. (2015) Significance and research on depth of coalbed methane oxidation zone in Dacheng area of Hebei province. Coal Technology 34(9): 113–115.

13.

Lei

Wang

, et al. (2016) An integrated test of the multi-channel transient electromagnetic system. Chinese Journal of Geophysics 59(12): 4399–4407.

14.

Zhu

Xue

, et al. (2019) New development of the electromagnetic (EM) methods for deep exploration. Chinese Journal of Geophysics 62(6): 2128–2138.

15.

Efthimios

Michael

(2000) Fast imaging of TEM data based on S-inversion. Journal of Applied Geophysics 43: 15–32.

16.

Wang

, et al. (2019) Elimination of mutual inductance effect for small-loop transient electromagnetic devices. Geophysics 84(3): E143–E154.

17.

Gao

Han

Liu

(2016) Research on the forward model DC method of abnormal body in full space based on ANSYS. Progress in Geophysics 31(5): 2089–2094.

18.

Gayle

Donald

and Rimstidt (1984) Natural weathering of coal. Fuel 63(March): 292–296.

19.

Guo

Xue

Liu

, et al. (2020) Electromagnetic methods for mineral exploration in China: a review. Ore Geology Reviews 118: 103357.

20.

Hisatoshi

Toru

Akira

, et al. (2011) Further investigations of underground resistivity structures in coastal areas using grounded-source airborne electromagnetics. Earth, Planets and Space 63: e9–e12.

21.

Huang

(2014) Application of mine high-density resistivity method in detection of water-filled karst fissures. Mining Safety & Environmental Protection 41(5): 56–58 + 68.

22.

Huang

(2016) Advance detection of mine concealed water flowing structure by seismic prospecting method. Mining Safety & Environmental Protection 43(4): 79–81.

23.

Lai

Chen

Zhang

(2016) Application of shallow transient electromagnetic method in the detection of city road disease. Progress in Geophysics (In Chinese) 31(6): 2743–2746.

24.

(2019) Study on mine transient electromagnetic method inversion based on PSO-DLS combination algorithm. Coal Science and Technology (In Chinese) 47(9): 268–272.

25.

Meyer

Siemon

Becken

(2016) New airborne methods and procedures for the exploration of mineral resources: An overview of BGR activities. In: 86th Ann. International Mtg.Soc. Expi. Geophys, Expanded Abstracts.

26.

Qin

Wang

(2019) Fast resistivity imaging of transient electromagnetic using ANN. IEEE geoscience and Remote Sensing Letters 16(9): 1373–1377.

27.

Shan

(2014) Application of seismic multi-attribute analysis technology in small colliery gob detection. Coal Mining Technology 19(5): 23–25.

28.

Sun

(2018) Geological and geophysical features of weathered, oxidized coals. Coal Geology of China 30(s2): 70–73.

29.

Tohru

Yoshikazu

Ken’ichirou

, et al. (1998) Development of grounded electrical source airborne transient EM (GREATEM). Exploration Geophysics 29(1-2): 61–64.

30.

Wang

Liu

(2014) Research on the safety mining technology of the seam under wind oxidized zone in Xiqu Coal Mine. Shanxi Coking Coal Science & Technology 1: 8–13.

31.

Zhi

Deng

, et al. (2020) Exploration of deep geological structure of baiyun gold deposit in eastern Liaoning province with TEM. Earth Science 45(11): 4027–4037.

32.

Xue

Liu

Sakyi

, et al. (2014) Discovery of a hidden BIF deposit in Anhui province, China by integrated geological and geophysical investigations. Ore Geology Reviews 63: 470–477.

33.

Xue

(2017) New development of TEM research and application in coal mine exploration. Progress in Geophysics 32(1): 319–326.

34.

Xue

Zhang

Zhou

, et al. (2020) Development’s measurements of TEM sounding in China. Geological Journal 55(3): 1771–1780.

35.

Yang

(2019) Development mechanism of wind oxidation zone in large buried deep coal seam and its influence on mining. China Energy and Environmental Protection 41(8): 78–83.

36.

Yee

(1996) Numerical solution of initial boundary problem involving Maxwell’s equations in isotropic media. IEEE Transactions on Antennas and Propagation 14(3): 302–307.

37.

You

Huang

(2015) Acoustic velocity logging curve based compartmentalization for bedrock air oxidized zone - taking supplementary hydrogeological exploration first mining area in Jinjiaqu for example. Ningxia Engineering Technology 14(4): 289–292.

38.

Liu

Tang

, et al. (2007) Transient electromagnetic detecting technique for water hazard to the roof of fully M echanized sub-level caving face. Journal of China University of Mining and Technology 36(4): 542–546.

39.

Zhan

Wang

(2016) Comprehensive geophysical prospecting technology for wind oxidize belt of inner coal seam on driving face. Shanxi Coal 36(5): 53–56.

40.

Zhan

Wang

Liu

, et al. (2015) Application of high-resolution seismic wave reflection technology in coal mine oxidation zone detection. Coal Engineering 47(1): 53–56.

41.

Zhang

(2008) Prospecting coal face structure by mine high-density electrical method. Coal Technology 27(10): 112–114.

42.

Zhao

(2008) Neogene system hydro-geological features and mining research with retaining wind oxidation zone pillars. Coal Technology 27(11): 56–58.

43.

Zhi

Wang

(2016) Three-component interpretation technique of fixed source TEM and its experimental application in metallic ore district. Geophysical and Geochemical Exploration 40(4): 798–803.

44.

Zhou

Liu

(2018) Research on the 3D mimetic finite volume method for loop-source TEM response in biaxial anisotropic formation. Chinese Joumal of Geophysics 61(1): 368–378.

Exploring weathered and oxidized zones of coal seams using the transient electromagnetic method

Abstract

Keywords

Introduction

Methodology

Geophysical characteristics of a coal seam weathered and oxidized zone

TEM methodology

Method selection basis

Data processing

Numerical modeling

Electromagnetic method based on the finite volume method

Numerical simulation result

Example analysis

Geological setting

Method validity test

Application Effect

Discussion

Conclusions

Footnotes

Declaration of conflicting interests

Funding

Data availability

ORCID iDs

References