Quantitative Calculation of Aquifer Water Quantity Using TEM Data

the disadvantage is the high workload and low efficiency. Therefore, transient electromagnetic method (TEM) was proposed, and TEM data was processed to detect the water yield property of a targeted layer in the coal mine. Through a mine-oriented 3D Transient Electromagnetic Method observing system, the signal of induced voltage is obtainable. Transient Electromagnetic Method 3D data volume can be calculated through the calculation of all-time resistivity and time-depth conversion. After an appropriate apparent resistivity value is set, the spatial distribution range of an aquifer can be determined. Then, with the water-filling coefficient of the aquifer, its water quantity can be estimated. The water yield property detection results in the No.4 coal seam goaf of the No.80101 workface in Jude Mine of Shanxi, China, demonstrates that the apparent resistivity of this goaf is less than 3 Ω.m, and the projection area of low resistivity anomaly zone is 22,383 m2. By using the formula Q=KMS, we can estimate that the water volume is 33,574 m3. Three boreholes have been constructed for the next dredging and drainage project, which results in a total water yield of 33,089 m3. The error percentage of the predicted water quantity is less than 1.5%. It can thus be concluded that it is feasible to predict aquifer water content with TEM data. ABSTRACT

El agua al interior de las minas es una amenaza para la minería de carbón. La cantidad del agua en las minas se detecta al perforar en las vías del socavón, con la desventaja de la carga de trabajo y la baja eficacia que significan. Por esto, el presente trabajo propone el Método Electromagnético Transitorio (TEM, del inglés Transient Electromagnetic Method) y utiliza la información generada para detectar las propiedades de producción de agua en una capa específica de una mina de carbón. A través de un sistema de observación tridimensional TEM orientado a minas se puede obtener la señal del voltaje inducido. La capacidad de la información del sistema TEM tridimensional se puede calcular al obtener la resistividad y la conversión tiempo-profundidad. Tan pronto se establece el valor apropiado de la resistividad aparente se puede determinar el rango de distribución de un acuífero. Luego, con el coeficiente de saturación del acuífero se puede calcular la cantidad de agua. Los resultados de la detección de las propiedades de producción de agua en la mina abandonada número cuatro de frente de trabajo número 80101 de la mina Jude, en Shanxi, China, demostró que la resistividad aparente de esta excavación es de 3 Ω.m, y el área de proyección de la zona anómala de baja resistividad es 22.383 m 2 . Al usar la fórmula Q=KMS, se estimó que el volumen de agua es de 33,574 m 3 . Tres pozos se han construido para el próximo proyecto de dragado y drenaje, cuyos resultados de producción de agua es de 33,089 m 3 . El porcentaje de error de la predicción de la cantidad de agua es menor al 1.5 %. Se concluye, por ende, que es factible producir el contenido de agua de un acuífero con la información del TEM.

Introduction
China is the largest producer of coal in the world because coal plays an essential role in the economic development of China ( Figure 1). Coal mining accidents have become a frequent occurrence in recent years (Zhang et al., 2014;Wang et al., 2016). Although the death rate resulting from these accidents in China has declined, it is still 70 times higher than that in the United States and 17 times higher than that in South Africa. Against this backdrop, it is urgent to study how to find effective ways to detect water quantity in the aquifer and reduce coal mining accidents (Zhang et al., 2015). Among some geophysical methods currently available, TEM technology is one that has received much attention and has been widely used to reduce coal mining accidents. In a TEM application, an ungrounded loop-line first transmits downward a pulse-type primary field, which induces an eddy current. The spatial and temporal distribution of a secondary field caused by the eddy current can then be observed through the coil (Yu, 2007;Wang et al., 2011;Tanaka and Kunisada 2011). By measuring the variation rule of the secondary magnetic field over time during the timeoff period, the geoelectricity features at different depths can be obtained (Guillemoteau et al., 2011;Tuncer et al., 2014). From Figure 2 it can be seen that TEM receives induced signals during the time-off period, which will be attenuated over time.
By using the TEM method, we can solve some geological problems Hu et al., 2013). Due to its features such as light equipment, little lateral influence, and high resolution, TEM can be widely applied in water disaster prevention and control (Zhang et al., 2010;Mollidor et al., 2013;Xue et al., 2013;Tao et al., 2013). However, because of the focus on theory development and processing method restraints, past analyses have usually been qualitative, which cannot meet the standard of quantification, and, to certain extent, hinders the application of TEM on a wider basis (Xu et al., 2012;Wang et al., 2014). Hence, in this study, we will develop a theory and propose a method to calculate TEM responses. To do so, we will use equivalent substitution, interpret the diffusion law of TEM field from the perspective of physics, calculate physical parameters, and conduct 3D apparent resistivity imaging and quantitative analyses of detection results.

Figure 2
Schematic diagram of TEM survey.

Apparent Resistivity Calculation
The transmitting loop center induction electromotive force of horizontal round loop in uniformity whole space medium is: Where is the vacuum permeability, n is the turn number of transmitting loop, I is the respective emission current, S is the receiving loop similar area, 0 r is the transmitting loop radius, t is the observation time and is the true resistivity of the medium.
From Equation (3), we can obtain an equation as below: (4) If the loop forms a square about b meters on a side, and , then, with Equations (1) ~ (4), the all-time apparent resistivity can be calculated with the binary search algorithm. The calculation process is illustrated in Figure 3.
For any observation window t, u is set as [0,10], and the searching region of u can be seen as an absolute value, and among them, us is the initial value of the searching region, ue is final value, ue is a median value which is u value. u is calculated by a binary search algorithm. is actual measured magnetic value, and H(ti ) is:

Depth Conversion
The detection depth D of the TEM method is related to transmitting magnetic torque, all-time apparent resistivity, and the minimum resolutive voltage. The time-depth conversion equation for mine-oriented TEM method is: (10) (6) where is the conversion coefficient, which is about the heterogeneous conductor in the vicinity of the roadway and usually values between 0.6 and 1.8; is the factor of proportionality and , where is the corresponding resistivity value of moment i t . (11)

Interpolation Operation
Due to the unstable workload, the measured and computational data of experiment and field are always not accurate enough, so interpolation is introduced as a supplement.
( , , ) k k k K x y z is an arbitrary point in the three-dimensional space. There is given point which is the vicinity of K.The inverse distance weighting method can be used to attribute the value XK with the equation: Where i d is the distance between the interpolating point and the given point in the vicinity. YK and ZK can be calculated in the same way as XK .

3D Slice Extraction
3D visualization is a technique to reveal and describe spatial data, receive the subsurface geological structure and features, which provides information for accurate description of the 3D geological structure, and facilitate the exploration and development of coal mine.
According to the data received from different detection directions, it can be seen that 3D visualization can clearly display the spatial distribution of anomalies with different resistivity values that correspond to different threshold values, the scattered points, and the range of apparent resistivity. However, a cross-section map can be clearer in showing the anomaly areas at a particular depth on a profile. Cross-sections can be extracted from different depths of Plane XY, YZ, and ZX. Apparent resistivity data can be extracted from any depth or vertical profile.

Aquifer Water Quantity Estimation
Through the cross-section map, the aquifer water volume can be calculated with the following formula: (12) In this equation, Q is the aquifer water quantity, K is the water containing coefficient, M is the aquifer thickness, and S is the water-rich area in the aquifer.

Geological Setting
The No.80101 Workface of Jude Mine in Shanxi Province is under L1 limestone in the middle and lower part of Taiyuan formation, 60.6m beneath the 4# seam. The seam occurrence is stable with a thickness range of 1.8 to 3.2m, and a 2.8m thickness on average. Part of the area contains a 0.05 to 0.1m dirt band of mudstone or carbonolite (carbonaceous mudstone). At the bottom of the seam is a pyrite seam in forms of lamella and nodule. The lithological column is shown in Figure 4.

Field Detection
To detect the spatial distribution of the goaf and estimate the water quantity, MTEM was adopted with the YCS360 Electromagnetic System (as shown in Figure 5). The detection of the No.80101 workface starts at the intersection of the return airway and the concentrated rail roadway, then goes on from the return airway to the openoff cut, and ends at the intersection of intake airway and concentrated rail roadway. The total detection workload is 3470 m, including 1660 m along the return airway, 150 m in the openoff cut, and 1660 m along the intake airway. Monitoring points are set every 10 meters, and the detection is conducted in 6 directions, including 60°, 45°, and 30° to the external wall roof, standard to the ceiling, and 30° and 45° to the interior wall roof. The TEM detection direction diagram is shown      in Figure 6, and the arrangement diagram of MTEM is shown in Figure 7.

Data Interpretation and Result Analysis
The data is processed based on the standard flow chart ( Figure  8), and data interpretation is carried out from the perspective of geophysical characteristics. The features are shown in 3D images. Water quantity is higher in areas where the fractures are concentrated, indicated by higher electrical conductivity and lower resistivity (i.e., high potentiallow resistivity) anomalies. This characteristic identifies aquifer areas as those where roof resistivity is less than 3Ω.m. Figure 9 is a 3D spatial distribution graph of anomalies at different depths above the No.80101 workface. 40 meters above the roof is a water-rich stratum, and key preventions should be focused on Region YC1, YC7, and YC10. Anomaly area YC7 is located between 1085m and 1850 m along the return airway (between monitoring point KF13 to KF15), 0-120m along the inclination slope (120m within return airway and workface), and 40m to 60m above the roof (more details in Fig 9), with an anomaly-impacted area of 22,383 m2. The goaf water-filling coefficient is 0.25 to 0.35, and the detection result of the actual ratio is 0.3. It is indicated by the effect sketch that M is 5m, based on the equation Q=KMS. Thus we can deduce that Q=33,574m³. The verification of the calculated values based on field measurements is shown in Table 1. H19-2, H19-3 and H19-S1 boreholes have been constructed for the later next and drainage project, which results in a total water quantity of 33,089m3, and the error percentage of the predicted water volume is less than 1.5%.

Conclusions
TEM resistivity imaging boasts rapid imaging, high quantification level, and high resolution, especially when it comes to detecting low resistivity body with a high-resistivity overlayer. With the projection area S of low-resistivity anomaly zone, the water containing coefficient K and the aquifer thickness M, the aquifer water quantity can be estimated with an error percentage of the predicted water quantity below 2%. The TEM 3D presentation and analysis has offered a new method to determine the water amount in an aquifer.