Next Article in Journal
Deciphering Double-Walled Corrugated Board Geometry Using Image Analysis and Genetic Algorithms
Previous Article in Journal
Simulation of Human Movement in Zero Gravity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 24
Issue 6
10.3390/s24061771
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Combined Method of Seismic Monitoring and Transient Electromagnetic Detection for the Evaluation of Hydraulic Fracturing Effect in Coal Burst Prevention

by
Jiang Bian
1,
Aixin Liu
1,2,
Shuo Yang
1,
Qiang Lu
1,2,
Bo Jia
1,
Fuhong Li
3,
Xingen Ma
4,
Siyuan Gong
1,* and
Wu Cai
2,*
1
School of Mines, China University of Mining and Technology, Xuzhou 221116, China
2
State Key Laboratory for Fine Exploration and Intelligent Development of Coal Resources, China University of Mining and Technology, Xuzhou 221116, China
3
Hetaoyu Coal Mine, Huaneng Qingyang Coal Power Co., Ltd., Qingyang 745300, China
4
Huaneng Coal Technology Research Co., Ltd., Beijing 100070, China
*
Authors to whom correspondence should be addressed.
Sensors 2024, 24(6), 1771; https://doi.org/10.3390/s24061771
Submission received: 20 February 2024 / Revised: 6 March 2024 / Accepted: 7 March 2024 / Published: 9 March 2024
(This article belongs to the Section Physical Sensors)
Browse Figures
Review Reports Versions Notes

Abstract

:
In order to mitigate the risk of roof-dominated coal burst in underground coal mining, horizontal long borehole staged hydraulic fracturing technology has been prevailingly employed to facilitate the weakening treatment of the hard roof in advance. Such weakening effect, however, can hardly be evaluated, which leads to a lack of a basis in which to design the schemes and parameters of hydraulic fracturing. In this study, a combined underground–ground integrated microseismic monitoring and transient electromagnetic detection method was utilized to carry out simultaneous evaluations of the seismic responses to each staged fracturing and the apparent resistivity changes before and after all finished fracturing. On this basis, the comparable and applicable fracturing effects on coal burst prevention were evaluated and validated by the distribution of microseismic events and their energy magnitude during the mining process. Results show that the observed mining-induced seismic events are consistent with the evaluation results obtained from the combined seismic-electromagnetic detection method. However, there is a limited reduction effect on resistivity near the fractured section that induces far-field seismic events. Mining-induced seismic events are concentrated primarily within specific areas, while microseismic events in the fractured area exhibit high frequency but low energy overall. This study validates the rationality of combined seismic-electromagnetic detection results and provides valuable insights for optimizing fracturing construction schemes as well as comprehensively evaluating outcomes associated with underground directional long borehole staged hydraulic fracturing.

1. Introduction

Hydraulic fracturing, especially within the horizontal long borehole staged fracturing technology (L-shape hydraulic fracturing), was initially employed in the field of oil and gas extraction. By injecting fluid into the borehole through high-pressure pumps, fractures are stimulated in the rock when the fluid pressure exceeds its strength. Such fractures serve as channels to improve the extraction efficiency of oil, gas, and geothermal resources [1]. Since the 1990s, hydraulic fracturing technology has been gradually applied in coalbed methane gas management and coal burst prevention in coal mines [2,3]. Its objectives include enhancing coal seam permeability [4,5], increasing the gas extraction rate [6], weakening hard roofs [7,8], reducing the initial pressure step distance [9,10], and minimizing coal burst tendency [11,12]. In recent years, the application of directional hydraulic fracturing in coal burst prevention has been demonstrated in a few coal mines. For instance, Su [13] took Hejiata Coal Mine in Shenfu Coalfield as the application object, and verified that after hydraulic fracturing pre-split roof, the roof can be stratified and caved successively, and the initial strength of the old roof can be weakened. Liu et al. [14] adopted the hydraulic fracturing construction technology of ultra-thick hard roof, judged the degree of crack development through drilling and peeking, and verified the pressure relief effect according to the peak working resistance and load increase coefficient of the support. Pang et al. [15] designed the fracturing scheme based on the key layer theory and Hubbert–Willis non-permeable formation rupture pressure calculation formula, and verified the good effect of directional hydraulic fracturing with long horizontal holes by using the periodic pressure step distance and advanced support pressure.
Currently, the monitoring techniques for hydraulic fracturing prevention and control technology primarily encompass borehole imaging, well temperature logging, transient electromagnetic detection [16], direct current detection, microseismic monitoring [17,18], and construction parameter analysis. Men et al. [19] adopted ARAMIS/E microseismic monitoring and ARES-5 ground sound monitoring to conduct the joint monitoring of hydraulic fracturing of ultra-long boreholes in the roof of ultra-thick coal seam. The results show that the fracture of the rock layer is mainly in a small area and has a high frequency during fracturing, which effectively reduces the energy accumulation of the overlying rock layer on the roof. Kim J et al. [20] used the crosswell electromagnetic (EM) geophysical methods to detect the fracture propagation induced by hydraulic fracturing; the resulting conductivity images provided a clearer understanding of the fractured reservoir, proving the EM can serve as an effective monitoring tool for monitoring the injected fluids during hydraulic fracturing. Kang et al. [21] used ground three-component microseismic monitoring to analyze the distribution of a fracture network, and used a mine pressure monitoring and warning platform, borehole stress gauge, underground microseismic monitoring, etc., to form a hydraulic fracturing analysis system for directional drilling in coal mines. Zhang et al. [22] analyzed the fracture area and fracturing influence range of coal and rock mass after fracturing by comparing the numerical simulation results with the transient electromagnetic detection. Peng et al. [23] used a vibration wave CT monitoring method and direct current method to reflect the extension degree of fractures around fracture holes and the connectivity and extension degree of primary fractures, and built a multi-variable and multi-scale detection method for the influence range of hydraulic fracturing.
In this study, the engineering background was provided by the fully mechanized top coal caving face longwall panel 1802 in Hetaoyu Coal Mine, where the application and analysis of combined seismic-electromagnetic monitoring technology for roof hydraulic fracturing effect evaluation were conducted. Based on the results obtained from simultaneous microseismic monitoring during fracturing and transient electromagnetic detection measuring apparent resistivity changes before and after fracturing, the effectiveness of the fracturing project was evaluated to establish a foundation for enhancing and optimizing the fracturing scheme. Finally, the effectiveness of this method is further verified by analyzing the pressure relief effect of each region through the microseismic events monitored during the mining process of the working face.

2. Project Overview

2.1. Geological Survey of the Mine

The Hetaoyu Mine is situated at Gansu Province, China, spanning approximately 14.728 km from east to west and 12.438 km from north to south, the mine covers an expansive area of about 183.189 km2. There are 3 coal seams (No. 2, No. 5, and No. 8) to be mined, in which the No. 8 coal seam was designed as the primary mining layer [24]. There are two fully mechanized caving faces in production at Hetaoyu mine: 2804 and 1802. The coal seam’s roof and floor rocks exhibit low solidity, weak rock strength, uneven local bump conditions, limited fissure development, expansion of both the roof and floor when encountering water, frequent floor occurrences, a high risk of roof collapse, as well as a significant presence of gas, making it a high-risk mine for explosion due to coal dust accumulation [25]. The three-dimensional geological survey of the Hetaoyu Mine is shown in Figure 1, within the location of drilling field 1 included.

2.2. Directional Long Borehole Staged Hydraulic Fracturing

In the context of fully mechanized mining of thick coal seams with hard roofs, the directional long borehole staged hydraulic fracturing is a technical measure for the regional advanced transformation of roof stress distribution.
The construction area of drilling field 1 of Hetaoyu Coal Mine studied in this paper is located in the initial mining area of the 1802 working face, and the mining of the working face passes through the square special structure area and the S5 anticline area, and the impact risk is high. In order to investigate the hydraulic fracturing relief effect of the key layer of overlying rock, three boreholes, No. 1, No. 2, and No. 3, were designed by directional drilling technology. According to the borehole column shape, the occurrence general situation of the strata above the coal seam was determined, and the fracturing plan of the key layer was designed, and the construction effect of the staged fracturing is shown in Figure 2. The hole opening position of drilling field 1 is about 530 m away from the open cut-off the working face 1802. All holes are drilled in single-hole mode. Table 1 describes the drilling parameters.
In the process of fracturing engineering implementation, the number of fracturing stages and the spacing of each fracturing stage should be adjusted according to the exposure of rock strata during drilling. The actual fracturing sequence of the first drilling field is borehole No. 3 (frac stages 1–9) → borehole No. 1 → borehole No. 2 → borehole No. 3 (frac stages 10–14). The hole depth of borehole No. 1 is about 531 m, and there are 12 stages of fracturing. The hole depth of borehole No. 2 is about 522 m and there are 14 stages of fracturing. Borehole No. 3 has a depth of approximately 573 m and a total of 14 stages of fracturing. Figure 3 shows the distribution of pressure fractures. Table 2 lists the length data of each fracturing stage.

3. Combined Seismic-Electromagnetic Detection Method

3.1. Underground–Ground Integrated Microseismic Monitoring System

The joint monitoring system utilized in this study primarily comprises the robust well vibration monitoring system and the subterranean SOS microseismic monitoring system. Building upon the pre-existing SOS microseismic monitoring system at Hetaoyu Coal Mine, this system meticulously designs and arranges a ground-based three-component strong seismograph to achieve the synchronized transmission of vibration data via the 4G network. Consequently, it establishes an all-encompassing three-dimensional envelope of the monitored area, enhancing the positioning accuracy of the microseismic monitoring system while reducing positioning errors in the z direction for individual underground pickups. This refined approach enables the more effective surveillance of vibration locations resulting from overlying rock fractures during hydraulic fracturing operations and facilitates the accurate assessment of fracture development levels [26,27]. The microseismic location diagram of the combined well and ground system is shown in Figure 4 below, and the microseismic parameters of underground and ground sensors are shown in Table 3.
For the purpose of achieving high accuracy in the source location, the selection of the longitudinal wave (P wave), which is more discernible, is typically favored for localization, as shown in Figure 5. Compared to other waves, the determination error of the first arrival time of the P wave is smaller and its positioning accuracy is higher. Our system uses the STA/LTA method to automatically identify the microseismic signal. Due to the low signal-to-noise ratio of the mining area, it is necessary to manually adjust the initial arrival time of the P wave, and the theoretical and actual calibration time of each sensor is marked. The deviation is controlled within ±20 ms. For signals with noise interference, the FFT algorithm within the low-pass filtering method is used for filtering and calibration. In practical applications, it is often assumed that the vibration wave propagates in a uniform and isotropic medium, making it very challenging to determine the propagation time using any propagation velocity. This assumption implies that the P wave maintains a constant velocity in all propagation directions. The formula below represents the minimum time taken for the vibration signal to travel from its source to the pickup location.
t i − t 0 = ( x 0 − x i ) 2 + ( y 0 − y i ) 2 + ( z 0 − z i ) 2 v x o , y 0 , z 0
In the formula, x 0 , y 0 , z 0 are the source coordinates; t 0 is the time of onset of seismic source;   x i , y i , z i are the coordinates of the seismic geophone;   t i is the arrival time of the P wave at the seismic geophone;   v x 0 , y 0 , z 0 is the velocity of P wave propagation in medium.
According to the fracturing area of drilling field 1 in the 1802 working face of Hetaoyu Coal Mine, considering the actual equipment installation and optimization of the microseismic monitoring effect during fracturing, a layout scheme comprising 13 sensors serving drilling field 1 was devised, as shown in Figure 6. Among them, a total of 7 SOS sensors were assembled downhole at the two roadways of 1802 working face; U19 is a large-range sensor, and its installation position is only 0.5 m away from the U18 sensor, which is used for the energy calculation of the large vibration. The rest are small-range sensors, which are used for the precise positioning of the microseismic source. And a total of 6 ET-GSY strong seismograph stations were assembled on the ground, which finally form an underground–ground envelope monitoring system for the fracturing area of drilling field 1.
By using the seismic wave passive CT inversion technology, the change in the impact dangerous area of drilling field 1 in 1802 working face before and after fracturing can be identified and analyzed. A microseismic monitoring system is utilized to calibrate the triggering time of downhole vibrations and determine the spatial distance between rupture positions and sensors. When vibration waves propagate through coal and rock masses, they induce changes in wave velocity due to stress effects [28].
The abnormal value of wave velocity can better reflect the stress distribution state of the detection area, thus providing a basis for the division of impact hazard area and the implementation of pressure relief prevention and control measures. The higher the abnormal value the wave velocity is, the stronger the stress concentration characteristics and the larger the concentration probability is, and the stronger the impact hazard is. The calculation method for determining the abnormal value of wave velocity is presented in Formula (2), while the correlation between the abnormal value of wave velocity and the probability of stress concentration is illustrated in Table 4.
A n = v p − v p a v p a
In the formula, v p is the P-wave velocity value at a point in the inversion region; v p a is the average of the model wave velocity.

3.2. Transient Electromagnetic Detection

The principle of mine transient electromagnetic detection is the principle of electromagnetic induction. A transmitting coil with a certain current is set up outside the detection area by using an ungrounded return line, and an induced current is generated in the conductive rock orebody around the detection area. The current establishes a stable magnetic field in the space around the transmitting coil. At any initial moment, when the original electromagnetic field is suddenly closed or changed, the underground free electrons cannot stop moving immediately because of inertia, so they will continue to produce secondary electromagnetic fields, which is the so-called “transient” process. By measuring the intensity and direction of the transient electromagnetic field, the position, shape, size, depth, and electrical parameters of the anomalous body can be obtained [29]. Figure 7 shows the detector host.
The primary objective of this study is to identify the development and expansion of hydrofractures subsequent to fracturing the roof of the working face. In rock masses characterized by good integrity and compactness, their resistivity tends to be relatively high. However, when fracturing fluid penetrates through the rock mass, leading to fracture generation and expansion, the resistivity in these fractured regions decreases due to fluid influx [30].
The design detection route is shown in Figure 8. Prior to and following fracturing in the drilling field 1, detection commenced from the 1802 cut hole, with a detection point every 10 m along the air return lane and transportation lane. To comprehensively monitor fissure development within the stratum, three detection angles were designated for each point of measurement, utilizing 2 m × 2 m electromagnetic coils to measure angles (15°, 45°, and 60°) relative to the coal wall of the working face. The detection schematic diagram is shown in Figure 9, and the field measurement is shown in Figure 10.

4. Application Results and Fracturing Effect Analysis

4.1. Analysis of Fracturing Construction Parameters

The accumulated fracturing fluid in borehole No. 1 is 458.9 m3, with the maximum construction pressure ranging from 22.4 MPa to 34.9 MPa and the maximum pressure drop ranging from 4.9 MPa to 12.8 MPa. The accumulated fracturing fluid in borehole No. 2 is 660.05 m3, with the maximum construction pressure ranging from 12.40 MPa to 33.6 MPa and the maximum pressure drop ranging from 0.7 MPa to 16 MPa. The accumulated fracturing fluid in borehole No. 3 is 581.75 m3, with the maximum construction pressure ranging from 16.1 MPa to 29.8 MPa and the maximum pressure drop ranging from 2.2 MPa to 12.9 MPa. Figure 11 is the pressure–flow curve during the water injection period of the fracturing stage. It can be reflected from the curve that the pressure of the fracturing stage rises to a critical value and then a sudden pressure drop occurs. This is due to the fact that the fracturing fluid breaks through and produces fractures. This critical value is the fracture initiation pressure, and then with the continuous injection of the fracturing fluid, the fractures are expanded, and the pressure and flow are in a serrated stable state.
The initiation pressure and the ratio of injection pressure to injection flow in each fracturing stage of the borehole are shown in Figure 12. The pressure changes of different fracturing stages are not regular in the region, indicating that the strength, integrity, fracture development, permeability, and other factors of the rock layer are not stable and vary greatly. In general, the initiation pressure (as shown in Figure 12a) can serve as an indicator of the maximum pressure exerted by the fracturing fluid on the packer during its passage through the rock layer. A higher value signifies a stronger rock layer in this particular area. Furthermore, the ratio of injection pressure to injection flow (as shown in Figure 12b) provides insights into the extent of rock fracture expansion within the scope of hydraulic fracturing [31,32]. However, it is difficult to accurately verify and evaluate the effect of fracturing through simple construction parameters. Therefore, a reliable geophysical detection method is needed more to obtain the data changes during and before fracturing, and to effectively evaluate the effect of fracturing in real time [33].

4.2. Analysis of Microseismic Monitoring Results

During the fracturing process of drilling field 1, the joint well-ground microseismic monitoring system provided real-time monitoring and successfully conducted an effective waveform analysis. In this regard, the joint well-ground monitoring system of borehole No. 1 recorded a total of 115 microseismic events in working face 1802, with a maximum microseismic energy of 19,522.5 J. The focal distribution is shown in Figure 13a.
A total of 311 microseismic events were recorded by the joint monitoring system of borehole No. 2 at the working face of 1802, with the maximum microseismic energy of 7755.6 J. The focal distribution is shown in Figure 13b.
A total of 270 microseismic events were recorded by the joint monitoring system of borehole No. 3 at the working face of 1802, with the maximum microseismic energy of 38,211.8 J. The focal distribution is shown in Figure 13c.
After cross-referencing the daily production operation of the mine, it was determined that the vibration events in Area A during the fracturing of the three boreholes were attributed to excavation activities related to the 1802 special drainage roadway. Similarly, in Area B, where the No. 2 borehole was being fractured, pipeline installation operations were underway. Hence, it can be inferred that both the vibration events in Area A and No. 2 borehole’s Area B are unrelated to hydraulic fracturing activities. Further examination revealed that the microseismic event in No. 2 borehole’s Area C resulted from blasting operations and not hydraulic fracturing-induced phenomena.
Figure 14 shows all frac microseismic events generated by borehole No. 1 and the corresponding frac stages. Among them, five microseismic events were induced during the fifth stage of fracturing, with very small energy levels, and the maximum energy recorded was only 54 J. Based on the coordinates and spatial location of these events, it can be observed that the vibrations occurred above the roof of the 1802 return air lane but at a distance greater than 60 m from the fracturing site. Furthermore, these five vibrations took place during a stable pressure and flow stage rather than when there was a significant reduction in water pressure. Therefore, it is concluded that these microseismic events were not directly caused by fracture propagation during water injection. Instead, two possibilities exist.
Firstly, after hydraulic fracturing altered the long-distance formation stress field, weak surfaces and natural fractures experienced breakage due to changes or disturbances in this stress field. The fact that these vibrations occurred at locations where roof leakage was detected suggests communication between fractures induced by vibration and those formed due to water pressure, thus forming a pathway for water flow [34,35].
Secondly, it is speculated that reduced cohesion resulting from water entering pre-existing fractures (which are relatively well developed at leakage sites) led to decreased effective normal stress (increased pore pressure), thereby inducing the fracture slip instability [36,37].
Subsequently, as shown in Figure 14, the occurrence of a mine earthquake event with an energy of 1.9 × 104 J outside the return air lane shortly after the seventh stage fracturing has been fully demonstrated by monitoring results, showing that it is challenging to induce high-energy microseismic events at the location of rock strata where hydraulic fracturing takes place. Therefore, this event is also attributed to far-field rock strata instability caused by hydraulic pressure. Additionally, it indicates significant stress conduction between hydraulic pressure and fractures during the hydraulic fracturing process, which cannot be achieved solely through small-scale fractures [38].
Furthermore, it has been confirmed that no operation was conducted during the fracturing of the borehole No. 3 B area. However, a microseismic event was induced during the fifth stage fracturing with very low energy (only 78 J), as shown in Figure 15. Considering the coordinates and spatial location of this microseismic event, vibrations occurred over a distance exceeding 60 m from the fracturing interval and took place during flow reduction and pressure reduction stages. Hence, it can be determined that this microseismic event is not directly caused by crack propagation in the water injection process but it is rather due to weak surface ruptures and natural fractures being disturbed or altered by changes in the long-distance formation stress field induced by hydraulic fracturing activities [39,40].
As shown in Figure 16, the microseismic waveforms before and after fracturing were marked to determine the source location, and the analysis focused on the 220 m mark level with the highest ray coverage density. The marked area in the figure is the tunneling operation of the outer drainage roadway of the transportation lane. These operations induce redistribution of the original stress field in the roof, leading to stress concentration. Apart from this tunneling area, a comparison of imaging data reveals a stress concentration zone with abnormally high wave velocity values extending approximately 400 m along the working face from the middle of drilling field 1 to the side of the transport lane prior to the fracturing construction. However, after the fracturing construction, both coverage area and abnormal wave velocity values within this zone are significantly reduced (as shown in Figure 16c). Furthermore, there is an increased extent of low abnormal wave velocity values observed through No. 1 and No. 2 boreholes, indicating that hydraulic fracturing has expanded fracture fissures in alignment with the maximum principal stress direction resulting in fractures within hard rock layers of roof strata [41,42].
For areas where abnormal wave velocity values remain higher than 0.25 even after fracturing (especially within intermediate regions between borehole No. 3 and the transport lane), additional pressure relief measures such as coal blasting should be reinforced to prevent potential rock burst incidents during the mining operation.

4.3. Analysis of Transient Electromagnetic Detection Results

According to the transient electromagnetic detection scheme of drilling field 1 formulated in Section 2.2, 53 points (530 m) were detected from the cutting hole of the working face 1802, each of which included three angles of 15°, 45°, and 60°. Two probes were conducted before and after fracturing. Figure 17 shows the apparent resistivity distribution of drilling field 1 at three detection angles. The apparent resistivity before fracturing is subtracted from the apparent resistivity after fracturing, which can reflect the fracture expansion after fracturing the high-pressure water injection, if the difference between the two is greater than 0, that is, the apparent resistivity decreases, indicating that the fracture expands to this area after fracturing [43].
The three-dimensional distribution of resistivity before and after fracturing along the horizontal and vertical directions of the borehole is illustrated in Figure 18 [44]. It can be observed that, apart from the influence of the advance bracket within a 150 m range near the cut hole, there has been a significant reduction in resistivity throughout the entire drilling field 1 in the horizontal direction. In terms of vertical direction, borehole No. 1 exhibits an overall superior fracturing effect with reduced resistivity near all stages of fracturing; particularly notable is the substantial decrease in resistivity at stages 5–12. On the other hand, while not evident during stages 1–4, borehole No. 3 demonstrates a remarkable reduction in resistivity during stages 5–14, indicating a significant fracturing effect [45,46].
By comparing the injection pressure–injection flow ratio distribution cloud chart and the CT inversion results of the vibration wave after fracturing, specifically examining the distribution of abnormal wave velocity values in Figure 16b, it can be observed that these values correspond to the spatial distribution area of apparent resistivity. Moreover, a positive correlation between them is evident; requisite fracturing is insufficient in areas with high wave velocity abnormal values and high apparent resistivity, thereby indicating that the integration of these two detection methods can effectively verify and analyze the effectiveness of fracturing.

4.4. Fracture Radius Analysis

According to the combined seismic-electromagnetic detection results, the fracturing radius is determined through analysis, specifically by calculating the distance of microseismic sources induced by fracturing and assessing the range of apparent resistivity reduction post-fracturing. This enables an evaluation of the efficacy of hydraulic fracturing at each stage [47,48].
As shown in Figure 19, the fracturing radius r 2 is satisfied as follows:
r 2 > L − r 1
In the formula, r 1 is the focal rupture radius, and L is the vertical distance between the fracturing location and the 1802 return air lane, which is 66 m.
Therefore, according to Formula (3), the fracturing radius can be estimated as long as the rupture radius of the source is calculated [49].
Generally, the larger the energy, the larger the rupture radius. Therefore, a microseismic signal with the largest energy induced by disturbance during the fracturing of the fifth section of borehole No. 1 is calculated to calculate the rupture radius of the source. Brune’s model is selected as the calculation model, and the calculation formula is as follows [50]:
r = K β 2 π f c
K is the constant selected based on Brune’ s model, P wave is 2.01, S wave is 1.32; β is the vibration wave velocity, P wave is 4670 m/s, S wave is 2696 m/s; f c is the corner frequency, expressed in Hz. As shown in the signal spectrum in Figure 20, the calculated average corner frequency is 65 Hz, and the calculated corner frequency, P-wave velocity and S-wave velocity of each channel are inserted into Formula (4), and the calculated average source radius is 16.1 m.
Therefore, by substituting r 1 = 16.1 and L = 66 into Formula (3), the calculation shows that the fracturing radius r 2 is at least greater than 49.9 m.
r 2 > L − r 1 = 66 − 16.1 = 49.9   m
The overall distribution of apparent resistivity detected by transient electromagnetic detection along the horizontal section of the borehole is illustrated in Figure 21. In proximity to the cut-off hole, the apparent resistivity detection value is influenced by the advance bracket and deviates from its actual value. It can be observed that fracturing has led to a certain reduction in apparent resistivity near the borehole, with some areas experiencing a reduction range exceeding 2 Ω·m, indicating effective diffusion of fracturing fluid.
Utilizing software (Voxler 4) analysis, we have determined that the radius encompassing areas with the rate of change of apparent resistivity greater than 0.5 (as shown in Figure 21c) ranges from 29.8–54.7 m, representing an effective zone for hydraulic fracturing operations. By integrating this information with seismic source rupture calculations, our comprehensive analysis suggests that the effective fracturing radius falls within a range of 29.8–54.7 m [51].

4.5. Fracturing Effect Evaluation

The pressure relief effect of roof hydraulic fracturing can be analyzed and evaluated based on microseismic data collected in the fractured area during face mining operations. Currently, there are two operational working faces at Hetaoyu Coal Mine: the 2804 fully mechanized caving working face and the 1802 fully mechanized caving working face. Among these, only the 1802 working face has implemented segmented hydraulic fracturing using directional long boreholes. By analyzing statistical data on pushing and mining advancement, it is observed that both working faces maintain a relatively stable mining speed during the initial stage of mining. To comparatively analyze daily frequency and total energy of microseismic events, all monitored mining microseismic events from the initial stage (first 200 m of advancement and first 50 days of stable pushing and mining) for both working faces were selected for analysis (as shown in Figure 22). A total of 1864 microseismic events with a total energy of 4,354,002.878 J were recorded in the working face 2804 at 200 m before the initial mining. A total of 622 microseismic events were recorded at the 1802 face, with a total energy of 1,495,330.004 J. A total of 995 microseismic events were recorded in the first 50 days of mining at the 2804 working face, with a total energy of 2,126,645.502 J. A total of 611 microseismic events were recorded at the 1802 face, with a total energy of 1,502,743.076 J.
By comparing the daily microseismic frequency and energy statistics of the two working faces depicted in Figure 22, it can be observed that the 1802 working face, which implements roof hydraulic fracturing, exhibits a general reduction in the number of microseismic events and total energy during the initial mining stage. Moreover, only a few high-energy events are induced. The total energy associated with each day’s microseismic activity is directly proportional to its vibration frequency. In Figure 22a,c, as the 1802 working face advances inward from the cut hole, there is an increase in microseismic events after passing through the hydraulic fracturing area; however, these events exhibit relatively low average energy levels, indicating a distribution characterized by “high frequency and low energy”. This observation reflects the effective pressure relief achieved through directional long borehole segmented hydraulic fracturing [52].
The spatial evolution of microseismic frequency and energy in the initial mining stage of 1802 and 2804 working faces is illustrated in Figure 23. For the 1802 working face with fracturing, microseismic events generated during mining are primarily concentrated on both sides of the working face. Moreover, compared to the non-fractured working face, there is a significant reduction in vibration frequency and energy, indicating that the segmented hydraulic fracturing construction of directional long boreholes has effectively alleviated pressure.
By considering the position of hydraulic fracturing boreholes in drilling field 1, it can be observed that the stress concentration phenomena occur near the fifth fracturing segment of borehole No. 1 and the fifth fracturing segment of borehole No. 3. These areas experience high-frequency vibrations with a substantial total energy release during mining operations, which aligns with well-ground joint microseismic monitoring results during fracturing activities as well as CT inversion results for seismic wave analysis (Section 4.2) and resistivity distribution along the borehole after fracturing (Section 4.4). This suggests that these specific regions exhibit high-strength overburden rock layers where extensive fracture expansion occurs during hydraulic fracturing processes, leading to far-field stress disturbances induced by such activities. Notably, unlike other areas where post-fracture resistivity reduction is more pronounced, this particular region demonstrates less obvious changes in resistivity after fracturing while still maintaining an elevated stress intensity within its overburden rock layer, consequently resulting in concentrated energy release phenomena during subsequent mining operations, thus validating the effectiveness of employing combined seismic-electromagnetic detection methodologies adopted within this study.

5. Conclusions

This paper adopts a combined seismic-electromagnetic detection method that combines an underground–ground microseismic monitoring with mining transient electromagnetic detection to quantify and analyze the spatial distribution of microseismic events and apparent resistivity during hydraulic fracturing in the hard roof of the 1802 working face at the Hetaoyu Coal Mine. The subsequent findings are as follows:
  • Hydraulic fracturing poses challenges in directly inducing microseismic events and generates low levels of vibration energy. However, the pressure generated by water injection and stress conduction resulting from fracture expansion can indirectly induce far-field vibrations. This is particularly evident in areas with naturally weak surfaces, highly fractured rock masses, or concentrated stress zones featuring faults, anticlines, and other unique geological structures. After fracturing, these rock layers may experience instability slips accompanied by elevated abnormal wave velocity values.
  • The transient electromagnetic detection of the apparent resistivity in different layers of coal and rock mass reveals that the roof rock layer exhibits a lower apparent resistivity after fracturing, indicating successful penetration of fracturing fluid into generated cracks and overall good fracturing effectiveness. However, certain areas still exhibit high apparent resistivity post-fracturing, aligning with the distribution of hydraulic flow and microseismic sources, suggesting suboptimal fracturing outcomes.
  • According to the microseismic monitoring data during the mining period in the fractured area, a significant reduction in both frequency and number of high-energy mining vibration events was observed across the entire working face after fracturing, as compared to the non-fractured working face. Furthermore, these findings align with those obtained from combined seismic-electromagnetic detection, confirming that the concentrated vibration areas are consistent. Additionally, it was found that microseismic events occurring in the regional mining post-fracturing exhibit a distribution characterized by “high frequency and low energy”. These results indicate that staged fracturing technology effectively alleviates stress concentration phenomena associated with hard roof strata while reducing the potential for rock burst disasters on the working face.
  • The combined seismic-electromagnetic detection accurately reflects the impact of hydraulic fracturing on the overlying strata. The results obtained from the combined approach of joint well-ground microseismic monitoring and transient electromagnetic detection are in correspondence. This combination strengthens the purpose of effect verification and facilitates analysis of potential shortcomings in fracturing operations. In future studies, it is imperative to precisely determine an optimal fracturing layer, ensure effective fracturing within critical rock layers, promote fracture development, and implement robust measures to mitigate erosion and relieve excessive pressure in areas where satisfactory fracturing outcomes have not been achieved.

Author Contributions

Conceptualization, J.B.; methodology, J.B. and S.Y.; software, A.L. and B.J.; validation, Q.L., A.L. and J.B.; formal analysis, S.Y and J.B.; investigation, S.Y., Q.L. and F.L.; resources, F.L.; data curation, A.L. and J.B.; writing—original draft preparation, J.B.; writing—review and editing, W.C. and S.G.; visualization, A.L. and B.J.; supervision, X.M.; project administration, S.G. and W.C.; funding acquisition, X.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (52374101), the Jiangsu Basic Research Program (Natural Science Foundation) (BK20231499), and the technology projects of the Huaneng Group headquarters (MTJSZX23-H05 and HNKJ21-H56-10).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The seismic data and apparent resistivity data used in this work are from Hetaoyu Coal Mine, and they are confidential.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tai, Y.; Yu, B.; Kuang, T.; Shi, B.; Meng, X.; Fu, M. An innovative technology of fracturing hard strata from the ground for precontrol of rock burst in a coal mine. Energy Sci. Eng. 2022, 11, 565–584. [Google Scholar] [CrossRef]
  2. He, H.; Cheng, R.; Zhao, J.; Men, Z.; Mu, Z. Research on the initiation pressure criterion of directional hydraulic fracturing in coal mine. Heliyon 2023, 9, e17638. [Google Scholar] [CrossRef]
  3. Yu, B.; Duan, H. Study of roof control by hydraulic fracturing in full-mechanized caving mining with high strength in extra-thick coal layer. J. Rock Mech. Eng. 2014, 778–785. [Google Scholar] [CrossRef]
  4. Guanhua, N.; Hongchao, X.; Zhao, L.; Lingxun, Z.; Yunyun, N. Improving the permeability of coal seam with pulsating hydraulic fracturing technique: A case study in Chang** coal mine, China. Process Saf. Environ. Prot. 2018, 117, 565–572. [Google Scholar] [CrossRef]
  5. Dunlop, E.C.; Salmachi, A.; McCabe, P.J. Investigation of increasing hydraulic fracture conductivity within producing ultra-deep coal seams using time-lapse rate transient analysis: A long-term pilot experiment in the Cooper Basin, Australia. Int. J. Coal Geol. 2019, 220, 103363. [Google Scholar] [CrossRef]
  6. Mou, Q.; Yan, Z.; Zhang, J. High efficiency gas drainage technology of hydraulic fracturing with di-rectional long drilling in underground coal mine. J. Coal Sci. Technol. 2020, 13, 296–303. [Google Scholar] [CrossRef]
  7. Li, H.; Huang, B.; Zheng, W.; Zhao, X.; Tannant, D. Effect of proppant distribution in hydraulic fractures on coalbed methane extraction. Results Eng. 2023, 20, 101550. [Google Scholar] [CrossRef]
  8. Huang, B.; Wang, Y.; Cao, S. Cavability control by hydraulic fracturing for top coal caving in hard thick coal seams. Int. J. Rock Mech. Min. Sci. 2015, 74, 45–57. [Google Scholar] [CrossRef]
  9. Mejia, C.; Roehl, D. Induced hydraulic fractures in underground block caving mines using an extended finite element method. Int. J. Rock Mech. Min. Sci. 2023, 170, 105475. [Google Scholar] [CrossRef]
  10. Liu, J.; Liu, C.; Yao, Q.; Si, G. The position of hydraulic fracturing to initiate vertical fractures in hard hanging roof for stress relief. Int. J. Rock Mech. Min. Sci. 2020, 132, 104328. [Google Scholar] [CrossRef]
  11. Zhu, S.; Liu, J.; Jiang, F.; Shang, X.; Sun, X.; Zhang, X.; Song, D.; Zhang, M.; Wang, A.; Xie, H. Classification, prediction, prevention and control of roof movement-type mineearthquakes and induced disasters in China’s coal mines. J. Coal 2022, 47, 807–816. [Google Scholar]
  12. Cao, A.; Bai, X.; Cai, W.; Wen, Y.; Li, X.; Ma, X.; Huang, R. Mechanism for stress abnormality and rock burst in variation zone of roof-stratum thickness. Chin. J. Geotech. Eng. 2023, 45, 512–520. [Google Scholar]
  13. Su, B. Application of hydraulic fracturing technology in hard and stable roof in coal mine. Coal Eng. 2019, 51, 54–57. [Google Scholar]
  14. Liu, W.; Li, G.; Liang, S.; Yang, C.; He, F. Control technology of hard roof hydraulic fracturing in extra-thick coal seam. Coal Eng. 2021, 53, 67–72. [Google Scholar]
  15. Pang, L.; Hu, Q.; Jing, J.; Jia, C.; He, W.; Hou, S.; Feng, T. Regional pressure relief technology with directional hydraulic fracturing of deep buried thick hard roof working face. Coal Eng. 2023, 55, 67–73. [Google Scholar]
  16. Zhao, R.; Fan, T.; Li, Y.; Wang, J.; Ma, Y.; Wang, B.; Liu, L.; Fang, Z. Application of borehole transient electromagnetic detection in the test of hydraulic fracturing effect. Coal Geol. Explor. 2020, 48, 41–45. [Google Scholar]
  17. Zhong, K.; Chen, W.; Zhao, W.; Qin, C.; Cao, H.; Xie, H. Monitoring and evaluation of segmented hydraulic fracturing effect in rock burst prevention on hard roof of coal mine. J. Cent. South Univ. (Nat. Sci. Ed.) 2022, 53, 2582–2593. [Google Scholar]
  18. Wang, J.; Zhang, H.; Zhao, L.; Zha, H.; Cheng, T.; Shi, X.; Qian, J.; Liu, Y.; Malinowski, M.; Linghu, J. Evaluation of coalbed methane hydraulic fracturing effect based on surface microseismic monitoring location and imaging. Geophys. Prospect. Pet. (Chin.) 2019, 62, 31–42+55. [Google Scholar]
  19. Men, H.; Zhao, H.; Dou, G.; Jia, Z.; Gao, Y.; Yan, B.; Xie, F.; Wu, L. Application analysis of hydraulic fracturing technology in extra thick coal seam roof with ultra-long boreholes. China Min. Ind. 2022, 31, 111–118. [Google Scholar]
  20. Kim, J.; Um, E.S.; Moridis, G.J. Integrated simulation of vertical fracture propagation induced by water injection and its borehole electromagnetic responses in shale gas systems. J. Pet. Sci. Eng. 2018, 165, 13–27. [Google Scholar] [CrossRef]
  21. Kang, H.; Feng, Y.; Zhang, Z.; Zhao, K.; Wang, P. Hydraulic fracturing technology with directional boreholes for strata control in underground coal mines and its application. J. Coal Sci. Technol. 2023, 31–44. [Google Scholar] [CrossRef]
  22. Zhang, D.; Yang, H.; Rao, Z.; Ou, Z.; Tang, P. Research on Application of Transient Electromagnetic Method in Hydraulic Fracturing. Geotech. Geol. Eng. Int. J. 2020, 38, 507–516. [Google Scholar] [CrossRef]
  23. Peng, Y.; Qiu, L.; Zhu, Y.; Liu, Q.; Song, D.; Cheng, X.; Wang, C.; Liu, Y.; Sun, Q. Multi-scale multivariate detection method for the effective impact range of hydraulic fracturing in coal seam. J. Appl. Geophys. 2023, 215. [Google Scholar] [CrossRef]
  24. Zhang, M.; Deng, H. Coal mining methods and safety evaluation of mining under deep water rich roof, rock burst and extra thick coal seam. J. Coal Technol. 2022, 9, 31–34. [Google Scholar]
  25. Qin, Z.; Chen, C.; Li, F.; Zhang, Y.; Du, T. Seismic mechanism and prevention technology of solid coal roadway in deep buried structural area. J. Coal Sci. Technol. 2021, 49, 87–92. [Google Scholar]
  26. Zhu, Y.; Wang, J.; Sun, F.; Lv, H.; Lin, J.; Chen, Z. Micro-seismic monitoring and instrument for hydraulic fracturing in the low-permeability oilfield. Chin. J. Geophys. 2017, 60, 4282–4293. [Google Scholar] [CrossRef]
  27. Wang, Y.; Liu, N.; Chen, F.; Li, Y.; Sun, X. Analysis of horizontal well staged fracturing effect based on the joint monitoring of surface and underground micro-seismic monitoring technology. Min. Strat. Control Eng. J. 2023, 5, 87–97. [Google Scholar]
  28. Li, Y.; Xu, D.; Ma, Z.; Zhou, H.; Guo, W. Application of CT inversion monitoring and early warning technology in microseismic anomaly area. J. Ind. Autom. 2021, 47, 39–45. [Google Scholar] [CrossRef]
  29. Li, X. Study on Response Characteristics and Application of Mine Transient Electromagnetic Monitoring; China Mining University: Xuzhou, China, 2023. [Google Scholar] [CrossRef]
  30. Han, D.; Li, D.; Shi, X. Effect of application of transient electromagnetic method in detection of water-inrushing structures in coal mines. Procedia Earth Planet. Sci. 2011, 3, 455–462. [Google Scholar] [CrossRef]
  31. Ma, S.; Huang, Y.; Zhang, L. Experimental Study on the Effect of Water Injection on Fault Stability and Its Enlightenment to the Mechanism of Induced Earthquake. 2020 China Joint Academic Annual Meeting of Geosciences. (eds.) Proceedings of the 2020 China Joint Academic Annual Meeting of Geosciences (6)-Project 16: Deep and Shallow Tectonic Characteristics and Dynamic Mechanisms of the North-South Seismic Belt and Northeast China, Project 17: Comprehensive Study of Focal Physical Processes and Seismic Hazard, Project 18: Strong Earthquake Mechanism and Prediction of Active Block Boundary Zone (pp.33); Beijing Botong Electronic Publishing House: Beijing, China, 2020. [Google Scholar]
  32. Huang, Y.; Ma, S.; Li, X. Research progress on water injection-induced earthquakes. Chin. J. Earthq. Eng. 2023, 45, 387–400. [Google Scholar] [CrossRef]
  33. Eyinla, D.; Henderson, S.K.; Emadi, H.; Thiyagarajan, S.R.; Arora, A. Optimization of hydraulic fracture monitoring approach: A perspective on integrated fiber optics and sonic tools. Geoenergy Sci. Eng. 2023, 231, 212441. [Google Scholar] [CrossRef]
  34. Naoi, M.; Chen, Y.; Nishihara, K.; Yamamoto, K.; Yano, S.; Watanabe, S.; Morishige, Y.; Kawakata, H.; Akai, T.; Kurosawa, I.; et al. Monitoring hydraulically-induced fractures in the laboratory using acoustic emissions and the fluorescent method. Int. J. Rock Mech. Min. Sci. 2018, 104, 53–63. [Google Scholar] [CrossRef]
  35. Sherratt, J.; Haddad, A.S.; Rafati, R. Modifying the orientation of hydraulically fractured wells in tight reservoirs: The effect of in-situ stresses and natural fracture toughness. Geomech. Energy Environ. 2023, 36, 100507. [Google Scholar] [CrossRef]
  36. Lei, P.; Wei, X.; Ru, L. Application of transient electromagnetic method technology in thedetection of shallow buried multilayer goaf. Shaanxi Coal 2023, 42, 94–99+112. [Google Scholar]
  37. Yu, B.; Zhu, W.; Li, Z.; Gao, R.; Liu, J. Mechanism of the instability of strata structure in far field for super-thick coal seam mining. J. Coal 2018, 43, 2398–2407. [Google Scholar] [CrossRef]
  38. Gao, R.; Kuang, T.; Meng, X.; Huo, B. Effects of Ground Fracturing with Horizontal Fracture Plane on Rock Breakage Characteristics and Mine Pressure Control. Rock Mech. Rock Eng. 2020, 54, 3229–3243. [Google Scholar] [CrossRef]
  39. Igonin, N.; Verdon, J.P.; Eaton, D.W. Seismic anisotropy reveals stress changes around a fault as it is activated by hydraulic fracturing. Seismol. Soc. Am. 2022, 93, 1737–1752. [Google Scholar] [CrossRef]
  40. Yu, B.; Gao, R.; Xia, B.; Kuang, T. Ground fracturing technology and application of hard roof in large space. J. Coal 2021, 46–48, 800–811. [Google Scholar] [CrossRef]
  41. Stoeckhert, F.; Molenda, M.; Brenne, S.; Alber, M. Fracture propagation in sandstone and slate–Laboratory experiments, acoustic emissions and fracture mechanics. J. Rock Mech. Geotech. Eng. 2015, 7, 237–249. [Google Scholar] [CrossRef]
  42. Moghaddam, R.H.; Golshani, A. Experimental study on fracture propagation in anisotropy rock under cyclic hydraulic fracturing. Eng. Fract. Mech. 2024, 295, 109775. [Google Scholar] [CrossRef]
  43. Zhang, Q.; Li, D.; Liu, G. An electromagnetic monitoring method based on underground charging conductor for hydraulic fracture diagnostics. Geoenergy Sci. Eng. 2023, 223, 211551. [Google Scholar] [CrossRef]
  44. Niu, Y.; Du, W.; Li, C. Research on 3D visualization method of transient electromagnetic detection data of mine working face. Shanxi Coal 2022, 42, 87–94. [Google Scholar]
  45. Li, Y.; Fan, T.; Zhao, R.; Liu, L.; Zhao, J. Application of three-component TEM in hydraulic fracturing effect detection of hard roof. J. Coal Sci. Technol. 2022, 50, 101–107. [Google Scholar]
  46. Cao, W.; Yildirim, B.; Durucan, S.; Wolf, K.-H.; Cai, W.; Agrawal, H.; Korre, A. Fracture behaviour and seismic response of naturally fractured coal subjected to true triaxial stresses and hydraulic fracturing. Fuel 2021, 288, 119618. [Google Scholar] [CrossRef]
  47. Jia, J.; Wang, D.; Li, B. Study on influencing factors of effective fracturing radius of hydraulic fracturing. China Work Saf. Sci. Technol. 2022, 18, 58–64. [Google Scholar]
  48. Cao, A.; Dou, L.; Jiang, H.; Lv, C.; Guo, X.; Wang, Y. Characteristics of energy radiation and stress drop in different failure modes of mining-induced coal-rock mass. J. Min. Saf. Eng. 2011, 28, 350–355. [Google Scholar]
  49. Li, Y. Microseismic analysis of hydraulic fracturing process. Acta Seismologica Sinica 1996, 18, 15–23+135. [Google Scholar]
  50. Jiang, P.; Dai, F.; Xu, N.; Li, T.; Li, B. Analysis of correlation between fracture scale and frequency characteristic of rock mass and its engineering verification. Rock Soil Mech. 2016, 5, 483–492. [Google Scholar]
  51. Xia, Y.; Deng, Y.; Jin, Y. Advances in Numerical Simulation of Fluid Flow in Fractured Reservoirs. Sci. Found. China 2021, 35, 964–972. [Google Scholar]
  52. Sun, R.; Xia, Y.; Gao, J. Study on anti-impact effect of hydraulic fracturing for long holes in medium and high thick hard roof. J. Coal Mine Saf. 2023, 69–77. [Google Scholar]
Sensors 24 01771 g001
Figure 1. Three dimensional geological map of Hetaoyu mine.
Figure 1. Three dimensional geological map of Hetaoyu mine.
Sensors 24 01771 g001
Sensors 24 01771 g002
Figure 2. Sectional fracturing construction effect of directional long drilling.
Figure 2. Sectional fracturing construction effect of directional long drilling.
Sensors 24 01771 g002
Sensors 24 01771 g003
Figure 3. Distribution of fracturing stages of each borehole in drilling field 1.
Figure 3. Distribution of fracturing stages of each borehole in drilling field 1.
Sensors 24 01771 g003
Sensors 24 01771 g004
Figure 4. Schematic diagram of underground–ground integrated microseismic monitoring system.
Figure 4. Schematic diagram of underground–ground integrated microseismic monitoring system.
Sensors 24 01771 g004
Sensors 24 01771 g005
Figure 5. Schematic diagram of source P wave marking location.
Figure 5. Schematic diagram of source P wave marking location.
Sensors 24 01771 g005
Sensors 24 01771 g006
Figure 6. Layout scheme of the underground–ground integrated microseismic sensors.
Figure 6. Layout scheme of the underground–ground integrated microseismic sensors.
Sensors 24 01771 g006
Sensors 24 01771 g007
Figure 7. Host of the transient electromagnetic meter.
Figure 7. Host of the transient electromagnetic meter.
Sensors 24 01771 g007
Sensors 24 01771 g008
Figure 8. Transient electromagnetic detection direction.
Figure 8. Transient electromagnetic detection direction.
Sensors 24 01771 g008
Sensors 24 01771 g009
Figure 9. Transient electromagnetic detection angle. (a) Return air lane; (b) transport lane.
Figure 9. Transient electromagnetic detection angle. (a) Return air lane; (b) transport lane.
Sensors 24 01771 g009
Sensors 24 01771 g010
Figure 10. Transient electromagnetic detection site of drilling field 1.
Figure 10. Transient electromagnetic detection site of drilling field 1.
Sensors 24 01771 g010
Sensors 24 01771 g011
Figure 11. Fracturing pressure–injection water flow curve diagram. (a) The fifth stage of borehole No. 1; (b) the fifth stage of borehole No. 3.
Figure 11. Fracturing pressure–injection water flow curve diagram. (a) The fifth stage of borehole No. 1; (b) the fifth stage of borehole No. 3.
Sensors 24 01771 g011
Sensors 24 01771 g012
Figure 12. Fracturing pressure cloud map of each fracturing section of boreholes in drilling field 1. (a) Fracture initiation pressure (MPa); (b) ratio of injection pressure to injection flow in the stable stage (MPa/m3).
Figure 12. Fracturing pressure cloud map of each fracturing section of boreholes in drilling field 1. (a) Fracture initiation pressure (MPa); (b) ratio of injection pressure to injection flow in the stable stage (MPa/m3).
Sensors 24 01771 g012
Sensors 24 01771 g013
Figure 13. Microseismic distribution map monitored during fracturing at drilling field 1.
Figure 13. Microseismic distribution map monitored during fracturing at drilling field 1.
Sensors 24 01771 g013
Sensors 24 01771 g014
Figure 14. Distribution of hydraulic fracturing microseismic signals in hole No. 1.
Figure 14. Distribution of hydraulic fracturing microseismic signals in hole No. 1.
Sensors 24 01771 g014
Sensors 24 01771 g015
Figure 15. Distribution of hydraulic fracturing microseismic signals in hole No.3.
Figure 15. Distribution of hydraulic fracturing microseismic signals in hole No.3.
Sensors 24 01771 g015
Sensors 24 01771 g016
Figure 16. Cloud chart of abnormal wave velocity value in drilling field 1.
Figure 16. Cloud chart of abnormal wave velocity value in drilling field 1.
Sensors 24 01771 g016
Sensors 24 01771 g017
Figure 17. Working face 0 m~530 m apparent resistivity contour profile.
Figure 17. Working face 0 m~530 m apparent resistivity contour profile.
Sensors 24 01771 g017
Sensors 24 01771 g018
Figure 18. Apparent resistivity distribution along borehole plane and vertical profile.
Figure 18. Apparent resistivity distribution along borehole plane and vertical profile.
Sensors 24 01771 g018
Sensors 24 01771 g019
Figure 19. Schematic diagram of fracturing radius calculation.
Figure 19. Schematic diagram of fracturing radius calculation.
Sensors 24 01771 g019
Sensors 24 01771 g020
Figure 20. Fracturing disturbance microseismic signal spectrum.
Figure 20. Fracturing disturbance microseismic signal spectrum.
Sensors 24 01771 g020
Sensors 24 01771 g021
Figure 21. Apparent resistivity distribution along the horizontal section of the boreholes.
Figure 21. Apparent resistivity distribution along the horizontal section of the boreholes.
Sensors 24 01771 g021
Sensors 24 01771 g022
Figure 22. Statistics of initial mining microseismic data of working face.
Figure 22. Statistics of initial mining microseismic data of working face.
Sensors 24 01771 g022
Sensors 24 01771 g023
Figure 23. Cloud map of spatial evolution of microseismic frequency energy in 1802 and 2804 working face.
Figure 23. Cloud map of spatial evolution of microseismic frequency energy in 1802 and 2804 working face.
Sensors 24 01771 g023
Table 1. Drilling design fracturing parameters of drilling field 1.
Table 1. Drilling design fracturing parameters of drilling field 1.
BoreholeDesign Hole Depth/mAperture/mmNumber of Stages
1#52812012
2#51912012
3#56012012
Table 2. Length of fracturing section of each borehole in drilling field 1.
Table 2. Length of fracturing section of each borehole in drilling field 1.
Stage of FractureFracturing Range of Hole No. 1/mFracturing Length of Hole No. 1/mFracturing Range of Hole No. 2/mFracturing Length of Hole No. 2/mFracturing Range of Hole No. 3/mFracturing Length of Hole No. 3/m
1512.27–519.857.58503.22–510.87.58555.99–562.576.58
2452.27–459.857.58476.22–483.87.58525.99–532.576.58
3422.27–429.857.58461.22–468.87.58495.99–502.576.58
4389.27–396.857.58440.22–447.87.58453.99–472.576.58
5359.27–366.857.58422.22–429.87.58410.99–417.576.58
6329.27–336.857.58398.22–405.87.58385.99–392.576.58
7305.27–312.857.58377.22–384.87.58334.99–341.576.58
8278.27–285.857.58278.22–285.87.58278.99–285.576.58
9251.27–258.857.58257.22–264.87.58221.99–228.576.58
10233.27–240.857.58251.22–270.819.58356.99–370.5713.58
11209.27–216.857.58230.22–246.816.58300.99–314.5713.58
12191.27–198.857.58209.22–225.816.58230.99–250.5719.58
13//191.22–207.816.58197.99–205.577.58
14//164.22–180.816.58161.99–169.577.58
Table 3. Basic parameters of microseismic sensors.
Table 3. Basic parameters of microseismic sensors.
Underground Sensors (SOS)Ground Sensorss (ET-GSY)
Range0.625 mm/s (small range); 0.5 m/s (Medium rang); 1 m/s (wide range)±1 g/±2 g/±4 g
Type of transmissionCurrent modeBuilt-in three-way EpiSensor force-balanced accelerometer
Bandwidth0.1~600 HzDC~200 Hz
Working temperature−5 °C~50 °C−30~+70 °C
Supply district18~42 V9~28 VDC
Table 4. Relationship between wave velocity outliers and stress concentration probability.
Table 4. Relationship between wave velocity outliers and stress concentration probability.
Stress ConcentrationWave Velocity Outlier A n Probability of Stress Concentration
Weak grade0~0.15<0.6
Medium grade0.15~0.250.6~1.4
Strong grade>0.25>1.4
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

Bian, J.; Liu, A.; Yang, S.; Lu, Q.; Jia, B.; Li, F.; Ma, X.; Gong, S.; Cai, W. A Combined Method of Seismic Monitoring and Transient Electromagnetic Detection for the Evaluation of Hydraulic Fracturing Effect in Coal Burst Prevention. Sensors 2024, 24, 1771. https://doi.org/10.3390/s24061771

AMA Style

Bian J, Liu A, Yang S, Lu Q, Jia B, Li F, Ma X, Gong S, Cai W. A Combined Method of Seismic Monitoring and Transient Electromagnetic Detection for the Evaluation of Hydraulic Fracturing Effect in Coal Burst Prevention. Sensors. 2024; 24(6):1771. https://doi.org/10.3390/s24061771

Chicago/Turabian Style

Bian, Jiang, Aixin Liu, Shuo Yang, Qiang Lu, Bo Jia, Fuhong Li, Xingen Ma, Siyuan Gong, and Wu Cai. 2024. "A Combined Method of Seismic Monitoring and Transient Electromagnetic Detection for the Evaluation of Hydraulic Fracturing Effect in Coal Burst Prevention" Sensors 24, no. 6: 1771. https://doi.org/10.3390/s24061771

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