Geofluids

Geofluids / 2020 / Article
Special Issue

Multiscale and Multiphysical Approaches to Fluids Flow in Unconventional Reservoirs

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 8845292 | https://doi.org/10.1155/2020/8845292

Jingqiang Tan, Jun Xie, Lei Li, Qiao Lyu, Jianqiang Han, Zhengguang Zhao, "Multifractal Analysis of Acoustic Emissions during Hydraulic Fracturing Experiments under Uniaxial Loading Conditions: A Niutitang Shale Example", Geofluids, vol. 2020, Article ID 8845292, 19 pages, 2020. https://doi.org/10.1155/2020/8845292

Multifractal Analysis of Acoustic Emissions during Hydraulic Fracturing Experiments under Uniaxial Loading Conditions: A Niutitang Shale Example

Academic Editor: Wei Wei
Received02 Jun 2020
Revised13 Jul 2020
Accepted14 Jul 2020
Published10 Aug 2020

Abstract

Fracture characterization is essential for estimating the stimulated reservoir volume and guiding subsequent hydraulic fracturing stimulations in shale reservoirs. Laboratory fracturing experiments can help provide theoretical and technical guidance for field operations. In this study, hydraulic fracturing experiments on the shale samples from Niutitang Formation in Hunan Province (China) under a uniaxial loading condition are conducted. The multifractal method is used to analyze the acoustic emission (AE) signals and characterize fracture initiation and propagation. The hydraulic fracturing process can be divided into three stages based on the characteristics of AE signals: the initial stage, the quite stage, and the fracturing stage. The multifractal analysis results showed that: (1) the value of the spectrum width, , continues to increase as the energy accumulates until the fracturing stage starts; and (2) the difference in the multifractal spectrum values, , reflects the relationship between small and large signal frequencies and can quantify the fracture scale, i.e., the lower the , the larger the fracture scale and vice versa. The results were further verified using a time-frequency analysis of the AE signals and micro-CT scanning of the samples. This study demonstrates that the multifractal method is feasible for quantitatively characterizing hydraulic fractures and can aid field hydraulic fracturing operations.

1. Introduction

The shale gas revolution in North America has affected the world energy market and greatly promoted the process of global shale gas exploration and production [1, 2]. Shale gas is usually stored as adsorbed or free gas in organic-rich mud shale, which has (extremely) low porosity and permeability [3, 4]. Hydraulic fracturing (HF) improves reservoir permeability by injecting high-pressure fluid into the ground to form a fracture network that is effectively connected to the wellbore [57]. Microseismic monitoring is generally implemented as a routine technology to monitor the development and expansion of fracturing networks in real time [8, 9]. Additionally, microseismic monitoring technology has also been widely used for monitoring in other fields, including geothermal production, mining, carbon dioxide storage, and other industrial fields [10, 11].

However, in-field microseismic monitoring, underground geological conditions are quite complex, such as the unknown distribution of natural/preexisting fractures and heterogeneous stress conditions. To effectively control the variables and focus when studying the evolution mechanism of HF, a feasible and important strategy is to simulate field-scale HF monitoring using laboratory-scale experiments with acoustic emission (AE) monitoring [12, 13]. AE is primarily related to the initiation and expansion of cracks. The frequency of laboratory- and small-scale acoustic emissions is generally in a range between several kHz to MHz. By analyzing the AE characteristics of rock fractures, it is helpful to understand the rock fracturing mechanism and provide a theoretical and technical basis for field HF monitoring [14]. Numerous studies have reported laboratory-scale rock fracture analysis using AE monitoring, including the rupture process of brittle rock [15, 16], the mechanical and acoustic properties of shale rocks [17, 18], the characteristics of fracture propagation with bedding planes at different directions [19], and the focal mechanism of AEs during laboratory-scale HF experiments [20, 21].

The processes of rock cracking and propagation have self-similarity and can be described by the fractal method (e.g., fractal dimension) [22]. There are many studies that have characterized rock fracture evolution using the fractal or multifractal analyses of the accompanying AE signals [2326]. For example, Kong et al. [26] analyzed the fractal characteristics of AE parameters of methane-containing coal samples in triaxial compression experiments. The results of the correlation dimension revealed that the fractal dimension could describe the propagation of a crack. Compared with the conventional fractal method, the multifractal theory has the potential to describe the rock rupture evolution state comprehensively [27, 28]. Kong et al. [29] used the multifractal method to study the AE characteristics of sandstone subjected to thermal treatment. The AE time series had multifractal features, and the fractal dimension followed an increasing-decreasing trend against the temperature. Kong et al. [30] applied the multifractal theory to analyze the AE characteristics of coal samples subjected to uniaxial compression. The results revealed that the time-varying multifractal characteristics could describe the AE mechanisms and the difference in strong and weak AE signals. Tan et al. [31] studied the multifractal characteristics of shale acoustic emission energy in uniaxial experiments under different immersion conditions. The results showed that the multifractal spectrum of the AE parameters was influenced by soaking times, which was directly related to changes in the pores and cracks of the shale samples. However, most of the above studies have focused on conventional compression experiments with intact samples. Therefore, quantitative multifractal analysis of AE parameters that results from shale hydraulic fracturing experiments requires further study.

In this study, the characteristics of AE events in hydraulic fracturing experiments under a uniaxial loading condition are analyzed to evaluate shale fracturing at different stages. The AE signals are quantified using multifractal analysis to study the rupture scale and energy characteristics. Finally, the reliability of the results is further verified using a time-frequency analysis and X-ray computed tomography (CT) scanning. The multifractal analysis can help provide quantitative guidance for fracture characterization during field hydraulic fracturing operations.

2. Materials and Methods

2.1. Sample Preparation

The samples were collected from outcrop shale samples of Niutitang Formation in Yiyang, Hunan Province, China. The outcrop shales were cored into cylindrical samples (size of ) with errors of less than 0.5 mm, and small cylindrical holes (size of ) were drilled at the center of the samples. The size and other physical properties of the six samples are shown in Table 1. Three of the six samples had horizontal bedding characteristics (the direction of bedding plane is parallel to the axial direction), which are termed as H1, H2, and H3. The other three samples were characterized by vertical bedding (the direction of bedding plane is perpendicular to the axial direction), which are termed as T1, T2, and T3. There were no visible fractures at the surface of the samples. The mineral components of the shale samples were quartz, with an average value of 74.2%; feldspar, with an average value of 14.3%; and clay minerals, with an average value of 7.8%. The samples tested were brittle rock and were not influenced by weathering. During the experiment, six acoustic emission (AE) sensors were placed at heights of 15 mm, 50 mm, and 85 mm, as shown in Figure 1(a).


Sample IDBedding directionMass (g)Diameter (mm)Height (mm)Density (g/cm3)P-wave velocity (m/s)

T-1Vertical461.549.8199.462.38/
T-2459.149.9699.542.353805
T-3463.449.999.672.383780
H-1Horizontal482.549.84100.332.474457
H-2492.449.91100.212.514450
H-3500.349.9499.082.584620

2.2. Experimental System and Scheme

Compared with the true triaxial experiments in field applications, there is no confining pressure under uniaxial loading conditions. It is no doubt that the confining pressure affects the characteristics of both the propagation of the fractures and the accompanying AE events. [3234]. However, hydraulic fracturing experiments under uniaxial loading conditions can better control associated variables and focus on the mechanism of fracture evolution. Besides, it can also clarify the fracture formation characteristic under different conditions, and reveal the fracture microscopic morphology in guiding field hydraulic fracturing [12]. Hydraulic fracturing experiments under uniaxial loading conditions are used to analyze AE characteristics associated with fracture propagation in various fields including geothermal production, salt well mining, and shale gas exploitation [35, 36].

The uniaxial compressive loading and hydraulic fracturing experiments were conducted using a servo-hydraulically controlled deformation apparatus, the GCTS RTR-2000 rock mechanic system at the Institute of Acoustics, Chinese Academy of Sciences. The system can provide a maximum axial force of 4600 kN and a maximum confining pressure of 210 MPa. Figure 1 shows the experimental system and setup for the shale samples. The PAC AE monitoring system was used for the experiments, which is made by the Physical Acoustic Corporation Ltd., USA for the real-time monitoring of AE signals during hydraulic fracturing (HF) experiments. The components of the acoustic emission monitoring system include AE sensors, preamplifiers, and a central system with recording, processing, and display units (Figure 1(d)). The sampling rate of the AE acquisition was 1 MHz, and the bandpass filter was set to 125–750 kHz. The threshold value of the AE signals was set to 35 dB. The two primary parameters obtained using the AE monitoring system included the AE counts and energy, which can reflect the deformation and fracture process of rock samples. The number of times that the signal exceeded the threshold was the AE counts, and the total area of the envelope for the AE waveform was defined as the AE energy. Other AE parameters and a typical AE waveform along with the time-frequency characteristics are shown in Figure 2.

The experimental scheme included the following five main steps: (1)Sealing Wellbore. A 50 mm long steel pipe was inserted into the injection hole (Figure 1(a)). Epoxy resin was used to seal the injection section of the shale sample. The sealed sample was preserved for 24 hours to ensure a stable sealing.(2)Sensor Installation. The heat shrinkage tube was closely bonded to the shale sample, and the AE sensors were placed on the tube surface using honey (Figures 1(a) and 1(c)).(3)Axial Loading. The specimen was loaded in the axial direction, and the axial stress value was set to 35 MPa, referring to the measured minimum compressive strength. The axial loading rate was 2 MPa/min.(4)Fluid Injection. Water was injected into the borehole at a rate of 2 MPa/min. When a sudden drop in pore pressure (injection pressure) occurred, this would indicate a break or damage of the sample and the system would shut down.(5)CT Scanning. After the HF experiments, the fractured samples were removed for micro-CT scanning that produced both 3D images and 2D slices of the cross sections.

2.3. Multifractal Method

The fractal method is used to describe the self-similarity of different local features in irregular phenomena in nature, and it is generally represented by the fractal dimension [37, 38]. Fractal analysis has been successfully applied to many branches in petrology and rock physics, such as characterizing the pores and fracture networks in rocks [3941]. To better quantitatively describe the spatial distribution characteristics of complex systems, the multifractal method has been proposed and utilized [27]. For example, the observed AE signals from rock samples could be multifractal. Hence, there is a need for more than one scaling exponent to characterize AE signals and related fractures [42]. The multifractal method has been used to describe the structural characteristics of rock samples using the generalized fractal dimension and multifractal spectrum [5]. In addition, many studies have adopted multifractal theory to study the AE activity of rock fracturing [26, 30]. In this study, the box-counting method was used to analyze the multifractal characteristics of the time series of the AE parameters [31, 43]. The detailed calculation process of the multifractal spectrum of the AE time series is as follows [31]:

First, a box with a width of was used to divide the time series into subsets. The normalized probability of each subset is , and the singularity strength is defined as:

When the time series has multifractal regularity, the following relationship can be obtained: where is the multifractal spectrum function, which can be also roughly defined as the fractal dimension with the same singularity subset. For calculating the multifractal spectrum of a time series, the following spectrum function was generally used: where is the quality index, ; and the quality index, , is the slope of the relationship of and under double logarithmic coordinates. Also, can determine the fractal characteristics of the time series: if the value of changes with , the AE time series has multifractal characteristics. If it is a constant value for a certain , the time series has a single regularity. Chhabra and Jensen [43] proposed the normalized single parameter measurement cluster as:

The multifractal spectrum function, , and singularity strength, , can be computed by:

When the value of is determined, and can be calculated, and we can obtain the multifractal spectrum of the time series. The heterogeneity of the time series can be indicated by the spectrum width [30, 44]. The larger is, the greater the difference between the large and small signals is. In addition, the value of represents the frequency at which the singularity strength, , occurs in the time series. Therefore, indicates the frequency of occurrence of the large and small signals: means the large-amplitude signals dominate, and the larger-scale fractures dominate the fracture system. In contrast, means the amplitude/energy of most AE signals is small, and small-scale fractures dominate the fracture system.

3. Results and Discussion

3.1. Mechanical Behaviors

During the fluid injection process, the pore pressure and AE signals were recorded. Table 2 lists the vertical stress and pore pressure loading on the shale samples. Figure 3 shows the relation of the pore pressure and loading time. The curves of the vertical samples were shifted 400 s horizontally to present a clear comparison between the horizontal and vertical samples. The maximum pore pressures of the horizontal samples were 14.52 MPa, 16.78 MPa, and 18.27 MPa, while those of the vertical samples were 26.53 MPa, 27.98 MPa, and 29.04 MPa, which are approximately double compared with the horizontal samples. Since the mineral composition was basically the same and the samples had no visible preexisting fractures, the results suggest that the bedding structure played a controlling role on the process of rock deformation and fracturing. The results are consistent with Lin et al. [45], and it showed the effects of the anisotropy on fracture propagation in shale samples. As the pore pressure increased, the effective normal stress decreased continuously. Once the breakdown pressure was reached, the fracture initiated along the direction of the minimum principal stress and propagated along that of the maximum principal stress, which was in the vertical direction here. The bedding direction of the horizontal samples was parallel to the axial direction, and lower pore pressures were required to break the samples compared with the vertical samples. In addition, the bedding direction of which was nearly perpendicular to the direction of the fracture propagation [4649]. The results are also consistent with the principle of the extensively-used horizontal well drilling technology in shale gas exploitation in which wells are drilled along the bedding direction (horizontal) to achieve a smaller breakdown pressure and a larger contacting volume with the target layer [50]. It is worth mentioning that shale anisotropy also has important effects on the field fracturing operation [51]. In the current study, we only focus on testing and validating the effectiveness of the multifractal method, without considering the shale anisotropy. We believe that further studies are required to address the anisotropic issues when generalizing the method to the field scale.


SamplesVertical stress (MPa)Injection rate (MPa/min)Maximum pore pressure (MPa)

H135214.52
H235216.78
H335218.27
T135226.53
T235227.98
T335229.04

3.2. Multifractal Analysis of AE Signals

There are several studies that have classified the fracturing process based on fracturing experiments. Damani et al. [52] reported that the fracture pressure response can be divided into six broad regions according to the pump pressure and cumulative AE events. Jiang et al. [53] divided the shale supercritical CO2 fracturing process into the fracture initiation stage and the fracture propagation stage according to the AE parameter curves.

Figure 4 shows the relationship of the AE counts, the AE energy, and the loading pore pressure for the horizontal samples (H1, H2, H3). Figure 5 shows the corresponding results of the vertical samples (T1, T2, T3). It is obvious that the AE counts and the AE energy shared a similar trend in all of the samples. In addition, the accumulative counts and energy steadily increased with the loading pore pressure. There is a difference between the detailed plots of the different samples, which was largely caused by the heterogeneity of the samples. As the loading pore pressure increased, the effective normal stress inside the rock sample decreased continuously. When the Mohr stress circle reached the rupture envelope surface, the rock broke down and generated a large amount of AE events [5, 54, 55]. Therefore, the hydraulic fracturing process can be roughly divided into three stages based on these plots [52, 53, 56]. First is the initial stage, which corresponds to the beginning of the water injection, where a small number of AE events occurred as the pore pressure increased instantaneously. Then, the quiet stage in which the fracturing system tended to be stable and the AE activity was inactive (“quiet”). Finally, the fracturing stage indicated the breakdown of the rocks and a sudden increase in the AE counts and energy. An important difference between the horizontal and vertical samples was the duration of the quiet stage. According to the dashed green lines in Figures 4 and 5, the duration of the quiet stage (stage II) for the horizontal and vertical samples was approximately 400 s and 700 s, respectively. This is consistent with the larger pore pressure required for the fracturing of vertical sample, which was largely caused by the direction of the bedding planes and discussed in the above section.

The multifractal method was used to further quantify the AE characteristics associated with the hydraulic fracturing experiments. It has been demonstrated that the time series of the AE parameters has multifractal characteristics [31]. The multifractal spectrum of the AE counts and energy time series throughout all stages is shown in Figure 6. The complete results for both the individual and overall stages are listed in Table 3. The results are closely related with the characteristics of fracture evolution. The values of the spectrum width, , of the AE counts for the six rock samples vary between two and four, which means the dispersion of the AE signals was relatively obvious. The primary reason is that the fracturing process can produce both large and small AE events/signals. The multifractal spectrum of the AE energy showed that the values were between three and five, which shows a consistent result of the AE counts. The entire process of shale hydraulic fracturing can be divided into different stages, including the initiation stage, the quiet stage, and the fracturing stage, and these naturally produce AE events with varying characteristics (e.g., counts and energy).


SampleCountsEnergy
Stage IStage IIStage IIIAll stagesStage IStage IIStage IIIAll stages

H11.2962.4403.2913.2081.1133.2434.0243.986
0.0070.607-0.1470.714-0.0290.301-0.2620.135
H20.0061.1972.5742.5570.0061.5553.8003.780
-0.003-0.4300.2641.643-0.003-0.0900.2821.639
H30.8411.0762.2912.9462.0082.3743.4933.874
0.142-0.1800.1291.0600.0610.035-0.0010.336
T13.1511.6083.9273.9213.9661.7384.7644.756
0.993-0.3360.2550.6900.804-0.3570.2970.586
T20.0041.1773.1333.3290.5211.8394.4024.596
0.0020.5520.000-1.884-0.5370.2140.011-1.839
T32.1541.6311.8042.1403.2872.3793.1923.261
0.567-0.0731.1431.6310.6670.0461.1211.562

The different values of the samples with the same bedding direction partially revealed the heterogeneity of the AE parameters/signals and that of the shales. In sample T3 for example, the values of the AE count and energy for all of the stages were 2.140 and 3.261, respectively. The values were relatively smaller than those of the other samples, indicating a small difference in the AE activities and/or fracture scales during the different stages. The difference in the multifractal spectrum values, , reflects the amplitude/energy of the AE events and the scale of the fractures. Most of the values of the horizontal samples were positive values, which means the fracturing process was dominated by small-scale fractures. The values of the AE count and energy for sample H2 were 1.643 and 1.639, respectively. The values were larger than those of the other two horizontal samples, and subsequent CT scanning results also showed several small-scale fractures (see Figures 7 and 8). However, the AE count and AE energy values of the vertical sample, T2, for all the stages were negative, and these were -1.884 and -1.839 for the AE count and energy, respectively. The CT scanning results of sample T2 showed large-scale fractures (see Figures 9 and 10) that resulted from the relatively large injection pressure and the connection between new fractures and weak bedding planes (preexisting fractures). More detailed discussions of the time-frequency analysis and CT scanning are presented in later sections.

Figure 11 shows the variations in and values for the different samples at different stages. For the horizontal samples, increased along with loading times and stages. During the initial stage, the injection pressure was small and only excited small AE events. As the injection pressure increased, the strain energy continuously accumulated. At the fracturing stage, small fractures were extended and connected to form large-scale fractures. Therefore, both small and large AE events were generated, and the difference between small and large signals increased. Several previous studies also revealed that the value rose gradually along with hydraulic fracturing [30, 44]. However, for the vertical samples, the results were more complex. A possible reason was the inconsistency in the directions of the principal stress and the bedding planes. And there were few AE events at the quiet stage and led to low values. For samples H2 and T2, the values of during different stages varied greatly, revealing a difference in the AE activity and/or fracture scales throughout all the stages. In addition, most large-scale fractures and AE events occurred during the final fracturing stage.

3.3. Time-Frequency Analysis

AE waveforms and frequencies contain important information during hydraulic fracturing, including internal fracture size and fracture type [21, 57, 58]. The microscopic evolution characteristics can be inferred by statistically analyzing the time-frequency characteristics of the AE signals during the fracturing process [35, 59]. First, the time domain AE signal was converted into the time-frequency spectrum using a wavelet transform to obtain the main frequency of its signal. AE events were classified to four types: low amplitude in low frequency band (LL), high amplitude in low frequency band (HL), low amplitude in high frequency band (LH), and high amplitude in high frequency band (HH) (Table 4). The maximum amplitudes and main frequency contents are statistically summarized in Figure 12. It can be seen that a small number of LL-type and LH-type AE events were generated during the initial stage of hydraulic fracturing, and the amplitude was on the order of 10-2. As the injection pressure gradually increased, the effective principal stress decreased continuously, and strain energy was accumulated. LL-type AE events were continuously generated during this stage. Once the breakdown pressure was reached, the crack initiated and propagated to form small-scale fractures, producing more LH-type AE events. Simultaneously, small-scale cracks might have been connected to each other, forming large-scale fractures and producing HL type AE signals. Therefore, large-scale and small-scale fractures coexisted simultaneously during the fracturing stage. This was consistent with the variations of the multifractal spectrum values at different stages and further verified the reliability of the results. And the HL- and LH-type AE signals were used as precursors before samples cracking [35]. Though the mineral composition of the rock samples is almost the same, the maximum pore pressures of the vertical and horizontal rock samples are also in the same range; the energy and the count values of AE events show a large difference. It is obvious that the anisotropy of shale seriously affects its physical properties. For the two horizontal rock samples H2 and H3, the number of events in the first and second stages is relatively small. However, in the third stage, the HH-type and the HL-type AE events of rock sample H3 are much larger than those of H2.


Low amplitude (0-0.1)High amplitude (0.1-5)

Low frequency band (0-200 kHz)Low amplitude in low frequency band (LL)High amplitude in low frequency band (HL)
High frequency band (200-400 kHz)Low amplitude in high frequency band (LH)High amplitude in high frequency band (HH)

3.4. CT Scanning

The internal structural characteristics of the horizontal sample H2 and the vertical sample T2 were studied using micro-CT scanning. Figures 7 and 9 show the overview of the micro-CT scanning results of rock samples H2 and T2. Figures 8 and 10 show four horizontal slices of the scanning results. The fractures initiated in the vicinity of the injection hole and continuously propagated along the axial direction (the vertical direction, which is also the direction of the maximum principal stress). The fractures were finally connected and formed a fracture plane, which was also nearly parallel to the axial direction. The morphological characteristics of the fractures are closely related to the confining pressure. When the axial pressure is much greater than the confining pressure, the fracture will propagate perpendicular to the direction of the local minimum principal stress, along with the borehole wall [60, 61]. The different scales of the cracks of the two samples were also consistent with the results of the multifractal spectrum values.

For the horizontal sample H2, the scale of the fractures was small, and only a few fractures were visible. The fractures remained in the vicinity of the injection hole and did not propagate vertically far away from the injection hole (see Figure 8). The fracture widths corresponded to heights of 20 mm and 40 mm and were relatively large, being 0.1–0.15 mm and 0–0.1 mm, respectively [62]. There were no visible fractures in the horizontal slices corresponding to the heights of 60 mm and 80 mm. A possible reason is that the injection (pore) pressure was small, and only small-scale fractures were produced. To quantitatively study the morphological characteristics of the fractures, an image processing software (Image J) was used to calculate the length and curvature of the visible fractures in the horizontal slices. The detailed results are listed in Table 5. The tortuosity was defined as the ratio of the total fracture length to the two direct lengths of the ends of the fracture [63]. The tortuosity of the visible fractures was small, and the fractures penetrated the rocks nearly straightly.


Direct length (mm)Length (mm)Width (mm)Tortuosity

H2I47.5047.690-0.151.00
II47.6248.110-0.151.01
III////
IV////

T2I48.9049.680.10-0.401.02
II47.9048.800.20-0.301.02
III46.8848.470.10-0.201.03
IV47.3647.830.05-0.151.01

For the vertical sample T2, the fracture width and scale were larger than those of horizontal sample H2. This was a result of the relatively larger maximum injection pressure. Fracture width gradually decreased from the bottom (0 mm) to the top (100 mm). The reason was that the fractures initiated at the bottom portion of the sample, where the injection hole was located. Another important and distinct characteristic was the nearly horizontal fracture plane at a height of approximately 80 mm, which was potentially caused by weakness of the bedding planes and/or nonvisible preexisting fractures. The fracture length, width, and tortuosity of samples H2 and T2 are quite different. The horizontal sample H2 is dominated by small-scale fractures, while the vertical sample T2 is dominated by large-scale fractures.

4. Conclusions

Hydraulic fracturing experiments were conducted under a uniaxial loading condition on shale samples from the Niutitang Formation in Yiyang, Hunan Province (China). The AE signals were recorded during the fracturing process, and the multifractal method was used to characterize fracture initiation and propagation. The study indicated that the multifractal method can be used to quantitatively characterize hydraulic fracturing and can contribute to optimize in situ hydraulic fracturing operations. The results were further verified using a time-frequency analysis of the AE signals and micro-CT scanning of the samples. The primary conclusions are listed below. (1)The direction of the bedding planes largely affected the required injection pressure that fractured the rock. The required pressure needed for rocks with bedding planes parallel to the direction of maximum principal stress was less (nearly half) than those with vertical bedding planes under essentially the same conditions.(2)The hydraulic fracturing process could be divided into three stages based on the characteristics of AE signals: an initial stage, the quiet stage, and the fracturing stage. During the initial stage, a small number of AE events were excited due to the instantaneous change in pore pressure. Also, the strain energy accumulated and a few AE events occurred during the quiet stage. During the fracturing stage, the fractures initiated and propagated promptly and produced a large amount of AE events. A multifractal analysis of the time series of the AE counts and energy quantitatively characterized the heterogeneity of the AE signals during hydraulic fracturing. The value of spectrum width, , continued to increase as the energy accumulated until the fracturing stage began. The value reflects the relationship between small and large signal frequencies and quantifies the fracture scale. Hence, the lower the , the larger the fracture scale, and vice versa.(3)The time-frequency analysis and CT scanning results further demonstrated the different magnitudes and scales of the AE signals and the fractures. Most AE events and large fractures were produced during the fracturing stage. The required fracturing pressure for the vertical samples was larger, while more large-scale fractures and complex fracture planes might be generated due to the existence of bedding planes that vertically or obliquely cross the maximum principal stress. The time-frequency characteristics can verify the change of the multifractal spectrum parameter . And the LH- and HL-type AE signals can be used as precursors of rock failure. CT scanning results were consistent with the multifractal spectrum parameter . These results further verified the stage classification and multifractal results of the AE signals.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The work is sponsored by the National Natural Science Foundation of China (grant no. 41872151), the National Thousand Young Talents Program (Jingqiang Tan, grant no. 413170110), the Innovation-driven Project of Central South University (grant no. 502501005), Hunan Provincial Natural Science Foundation of China (grant no. 2019JJ50762), the China Postdoctoral Science Foundation (grant no. 2019M652803), the Graduate Research Innovation Project of Central South University (grant no. 2020zzts661), and the Open Research Fund Program of State Key Laboratory of Acoustics, Chinese Academy of Sciences (grant no. SKLA201911).

References

  1. R. W. Howarth, A. Ingraffea, and T. Engelder, “Should fracking stop?” Nature, vol. 477, no. 7364, pp. 271–275, 2011. View at: Publisher Site | Google Scholar
  2. J. D. Hughes, “A reality check on the shale revolution,” Nature, vol. 494, no. 7437, pp. 307-308, 2013. View at: Publisher Site | Google Scholar
  3. C. Wang, L. Hong, R. Gao et al., “Status-quo and challenges of enhanced oil recovery in low permeability reservoirs,” Unconventional Oil & Gas, vol. 5, no. 3, pp. 102–108, 2018. View at: Google Scholar
  4. B. Liu, Y. Yang, J. Li, Y. Chi, J. Li, and X. Fu, “Stress sensitivity of tight reservoirs and its effect on oil saturation: A case study of Lower Cretaceous tight clastic reservoirs in the Hailar Basin, Northeast China,” Journal of Petroleum Science and Engineering, vol. 184, article 106484, 2020. View at: Publisher Site | Google Scholar
  5. K. Liu, D. Gao, Y. Wang, and Y. Yang, “Effect of local loads on shale gas well integrity during hydraulic fracturing process,” Journal of Natural Gas Science and Engineering, vol. 37, pp. 291–302, 2017. View at: Publisher Site | Google Scholar
  6. L. Wang, Y. Tian, X. Yu et al., “Advances in improved/enhanced oil recovery technologies for tight and shale reservoirs,” Fuel, vol. 210, pp. 425–445, 2017. View at: Publisher Site | Google Scholar
  7. P. Jian, Z. Zhuang, X. Tang, W. Duan, C. Guo, and T. Lin, “Comprehensive analysis for hydrocarbon accumulation main control factors and enrichment of unconventional oil and gas,” Unconventional Oil & Gas, vol. 4, no. 3, pp. 110–118, 2017. View at: Google Scholar
  8. S. C. Maxwell and T. Urbancic, “The role of passive microseismic monitoring in the instrumented oil field,” The Leading Edge, vol. 20, no. 6, pp. 636–639, 2001. View at: Publisher Site | Google Scholar
  9. L. Li, J. Tan, D. A. Wood et al., “A review of the current status of induced seismicity monitoring for hydraulic fracturing in unconventional tight oil and gas reservoirs,” Fuel, vol. 242, pp. 195–210, 2019. View at: Publisher Site | Google Scholar
  10. M. Weingarten, S. Ge, J. W. Godt, B. A. Bekins, and J. L. Rubinstein, “High-rate injection is associated with the increase in U.S. mid-continent seismicity,” Science, vol. 348, no. 6241, pp. 1336–1340, 2015. View at: Publisher Site | Google Scholar
  11. F. Grigoli, S. Cesca, E. Priolo et al., “Current challenges in monitoring, discrimination, and management of induced seismicity related to underground industrial activities: a European perspective,” Reviews of Geophysics, vol. 55, no. 2, pp. 310–340, 2017. View at: Publisher Site | Google Scholar
  12. M. S. Paterson and T. Wong, Experimental Rock Deformation—The Brittle Field, Springer, Berlin; New York, 2nd edition, 2005.
  13. C. Grosse and M. Ohtsu, Acoustic Emission Testing, Springer Berlin Heidelberg, Berlin, Heidelberg, 2008. View at: Publisher Site
  14. X. Lei, K. Masuda, O. Nishizawa et al., “Detailed analysis of acoustic emission activity during catastrophic fracture of faults in rock,” Journal of Structural Geology, vol. 26, no. 2, pp. 247–258, 2004. View at: Publisher Site | Google Scholar
  15. S. Stanchits, A. Surdi, P. Gathogo, E. Edelman, and R. Suarez-Rivera, “Onset of hydraulic fracture initiation monitored by acoustic emission and volumetric deformation measurements,” Rock Mechanics and Rock Engineering, vol. 47, no. 5, pp. 1521–1532, 2014. View at: Publisher Site | Google Scholar
  16. Z. Moradian, H. H. Einstein, and G. Ballivy, “Detection of cracking levels in brittle rocks by parametric analysis of the acoustic emission signals,” Rock Mechanics and Rock Engineering, vol. 49, no. 3, pp. 785–800, 2016. View at: Publisher Site | Google Scholar
  17. Y. Liang, Q. Li, Y. Gu, and Q. Zou, “Mechanical and acoustic emission characteristics of rock: effect of loading and unloading confining pressure at the postpeak stage,” Journal of Natural Gas Science and Engineering, vol. 44, pp. 54–64, 2017. View at: Publisher Site | Google Scholar
  18. J. Clarke, L. Adam, J. Sarout, K. van Wijk, B. Kennedy, and J. Dautriat, “The relation between viscosity and acoustic emissions as a laboratory analogue for volcano seismicity,” Geology, vol. 47, no. 6, pp. 499–503, 2019. View at: Publisher Site | Google Scholar
  19. Y. Wang, C. H. Li, and Y. Z. Hu, “Experimental investigation on the fracture behaviour of black shale by acoustic emission monitoring and CT image analysis during uniaxial compression,” Geophysical Journal International, vol. 213, no. 1, pp. 660–675, 2018. View at: Publisher Site | Google Scholar
  20. E. Aker, D. Kühn, V. Vavryčuk, M. Soldal, and V. Oye, “Experimental investigation of acoustic emissions and their moment tensors in rock during failure,” International Journal of Rock Mechanics and Mining Sciences, vol. 70, pp. 286–295, 2014. View at: Publisher Site | Google Scholar
  21. B. Zhang, X. Tian, B. Ji, J. Zhao, Z. Zhu, and S. Yin, “Study on microseismic mechanism of hydro-fracture propagation in shale,” Journal of Petroleum Science and Engineering, vol. 178, pp. 711–722, 2019. View at: Publisher Site | Google Scholar
  22. H. P. Xie, An Introdction of Fractal Methods on Rock Mechanics, Science Press, Beijing, 1996.
  23. H. P. Xie, J. F. Liu, Y. Ju, J. Li, and L. Z. Xie, “Fractal property of spatial distribution of acoustic emissions during the failure process of bedded rock salt,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 8, pp. 1344–1351, 2011. View at: Publisher Site | Google Scholar
  24. X. Z. Wu, X. X. Liu, Z. Z. Liang, X. You, and M. Yu, “Experimental study of fractal dimension of AE serials of different rocks under uniaxial compression,” Rock & Soil Mechanics, vol. 33, no. 12, pp. 3561–3569, 2012. View at: Google Scholar
  25. R. Zhang, F. Dai, M. Z. Gao, N. W. Xu, and C. P. Zhang, “Fractal analysis of acoustic emission during uniaxial and triaxial loading of rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 79, pp. 241–249, 2015. View at: Publisher Site | Google Scholar
  26. X. Kong, E. Wang, S. Hu, R. Shen, X. Li, and T. Zhan, “Fractal characteristics and acoustic emission of coal containing methane in triaxial compression failure,” Journal of Applied Geophysics, vol. 124, pp. 139–147, 2016. View at: Publisher Site | Google Scholar
  27. R. Lopes and N. Betrouni, “Fractal and multifractal analysis: a review,” Medical Image Analysis, vol. 13, no. 4, pp. 634–649, 2009. View at: Publisher Site | Google Scholar
  28. A. N. Pavlov and V. S. Anishchenko, “Multifractal analysis of complex signals,” Physics-Uspekhi, vol. 50, no. 8, p. 819, 2007. View at: Publisher Site | Google Scholar
  29. B. Kong, E. Wang, Z. Li, X. Wang, L. Chen, and X. Kong, “Nonlinear characteristics of acoustic emissions during the deformation and fracture of sandstone subjected to thermal treatment,” International Journal of Rock Mechanics and Mining Sciences, vol. 90, pp. 43–52, 2016. View at: Publisher Site | Google Scholar
  30. X. Kong, E. Wang, X. He, D. Li, and Q. Liu, “Time-varying multifractal of acoustic emission about coal samples subjected to uniaxial compression,” Chaos, Solitons & Fractals, vol. 103, pp. 571–577, 2017. View at: Publisher Site | Google Scholar
  31. J. Tan, C. Hu, Q. Lyu et al., “Multi-fractal analysis for the AE energy dissipation of CO2 and CO2+ brine/water treated low-clay shales under uniaxial compressive tests,” Fuel, vol. 246, pp. 330–339, 2019. View at: Publisher Site | Google Scholar
  32. Z. Cheng, Y. Fu, and C. Tang, “Confining pressure effect on acoustic emissions during rock failure,” Chinese Journal of Rock Mechanics and Engineering, vol. 16, no. 1, pp. 65–70, 1997. View at: Google Scholar
  33. T. Xu, C. Tang, S. Wang, and Y. Zhang, “Numerical tests on confining pressure effect in rock failure process,” Journal of Central South University(Science and Technology)(China), vol. 35, no. 5, pp. 840–844, 2004. View at: Google Scholar
  34. Y. Li, B. Wang, L. Song, J. Liu, J. Zuo, and D. Liu, “Ultrasonic wave propagation characteristics for typical anisotropic failure modes of shale under uniaxial compression and real-time ultrasonic experiments,” Journal of Geophysics and Engineering, vol. 17, no. 2, pp. 258–276, 2020. View at: Publisher Site | Google Scholar
  35. G. Chen, X. Sun, J. Wang, D. Wang, and Z. Zhu, “Detection of cracking behaviors in granite with open precut cracks by acoustic emission frequency spectrum analysis,” Arabian Journal of Geosciences, vol. 13, no. 6, p. 258, 2020. View at: Publisher Site | Google Scholar
  36. J. Yang, Z.-L. Mu, and S.-Q. Yang, “Experimental study of acoustic emission multi-parameter information characterizing rock crack development,” Engineering Fracture Mechanics, vol. 232, article 107045, 2020. View at: Publisher Site | Google Scholar
  37. B. B. Mandelbrot, “Fractals: form, chance, and dimension,” Physics Today, vol. 32, no. 5, pp. 65-66, 1979. View at: Publisher Site | Google Scholar
  38. C. C. Barton and P. R. Pointe, Eds., Fractals in the Earth Sciences, Springer US, Boston, MA, 1995. View at: Publisher Site
  39. R. Liu, Y. Jiang, N. Huang, and S. Sugimoto, “Hydraulic properties of 3D crossed rock fractures by considering anisotropic aperture distributions,” Advances in Geo-Energy Research, vol. 2, no. 2, pp. 113–121, 2018. View at: Publisher Site | Google Scholar
  40. A. Sakhaee-Pour and W. Li, “Fractal dimensions of shale,” Journal of Natural Gas Science and Engineering, vol. 30, pp. 578–582, 2016. View at: Publisher Site | Google Scholar
  41. Y. Xia, J. Cai, E. Perfect, W. Wei, Q. Zhang, and Q. Meng, “Fractal dimension, lacunarity and succolarity analyses on CT images of reservoir rocks for permeability prediction,” Journal of Hydrology, vol. 579, article 124198, 2019. View at: Publisher Site | Google Scholar
  42. L. Telesca, V. Lapenna, and M. Macchiato, “Multifractal fluctuations in earthquake-related geoelectrical signals,” New Journal of Physics, vol. 7, pp. 214–214, 2005. View at: Publisher Site | Google Scholar
  43. A. Chhabra and R. V. Jensen, “Direct determination of the f(α) singularity spectrum,” Physical Review Letters, vol. 62, no. 12, pp. 1327–1330, 1989. View at: Publisher Site | Google Scholar
  44. S. Hu, E. Wang, Z. Li, R. Shen, and J. Liu, “Time-varying multifractal characteristics and formation mechanism of loaded coal electromagnetic radiation,” Rock Mechanics and Rock Engineering, vol. 47, no. 5, pp. 1821–1838, 2014. View at: Publisher Site | Google Scholar
  45. C. Lin, J. He, X. Li, X. Wan, and B. Zheng, “An experimental investigation into the effects of the anisotropy of shale on hydraulic fracture propagation,” Rock Mechanics and Rock Engineering, vol. 50, no. 3, pp. 543–554, 2017. View at: Publisher Site | Google Scholar
  46. E. Detournay, “Propagation regimes of fluid-driven fractures in impermeable rocks,” International Journal of Geomechanics, vol. 4, no. 1, pp. 35–45, 2004. View at: Publisher Site | Google Scholar
  47. A. P. Bunger, E. Detournay, and D. I. Garagash, “Toughness-dominated hydraulic fracture with leak-off,” International Journal of Fracture, vol. 134, no. 2, pp. 175–190, 2005. View at: Publisher Site | Google Scholar
  48. D. Xu, R. Hu, W. Gao, and J. Xia, “Effects of laminated structure on hydraulic fracture propagation in shale,” Petroleum Exploration and Development, vol. 42, no. 4, pp. 573–579, 2015. View at: Publisher Site | Google Scholar
  49. X. Zhang, Y. Lu, J. Tang, Z. Zhou, and Y. Liao, “Experimental study on fracture initiation and propagation in shale using supercritical carbon dioxide fracturing,” Fuel, vol. 190, pp. 370–378, 2017. View at: Publisher Site | Google Scholar
  50. M. J. Economides and K. G. Notte, Reservoir Stimulation, Wiley, Chichester, UK, 3rd edition, 2000.
  51. M. Naoi, Y. Chen, K. Nishihara et al., “Monitoring hydraulically-induced fractures in the laboratory using acoustic emissions and the fluorescent method,” International Journal of Rock Mechanics and Mining Sciences, vol. 104, pp. 53–63, 2018. View at: Publisher Site | Google Scholar
  52. A. Damani, C. H. Sondergeld, and C. S. Rai, “Experimental investigation of in situ and injection fluid effect on hydraulic fracture mechanism using acoustic emission in Tennessee sandstone,” Journal of Petroleum Science and Engineering, vol. 171, pp. 315–324, 2018. View at: Publisher Site | Google Scholar
  53. Y. Jiang, C. Qin, Z. Kang et al., “Experimental study of supercritical CO2 fracturing on initiation pressure and fracture propagation in shale under different triaxial stress conditions,” Journal of Natural Gas Science and Engineering, vol. 55, pp. 382–394, 2018. View at: Publisher Site | Google Scholar
  54. M. Cacace, G. Blöcher, N. Watanabe et al., “Modelling of fractured carbonate reservoirs: outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany,” Environmental Earth Sciences, vol. 70, no. 8, pp. 3585–3602, 2013. View at: Publisher Site | Google Scholar
  55. Y. Zheng, J. Liu, and Y. Lei, “The propagation behavior of hydraulic fracture in rock mass with cemented joints,” Geofluids, vol. 2019, 15 pages, 2019. View at: Publisher Site | Google Scholar
  56. S. D. Goodfellow, N. Tisato, M. Ghofranitabari, M. H. B. Nasseri, and R. P. Young, “Attenuation properties of Fontainebleau sandstone during true-Triaxial deformation using active and passive Ultrasonics,” Rock Mechanics and Rock Engineering, vol. 48, no. 6, pp. 2551–2566, 2015. View at: Publisher Site | Google Scholar
  57. P. J. de Groot, P. A. M. Wijnen, and R. B. F. Janssen, “Real-time frequency determination of acoustic emission for different fracture mechanisms in carbon/epoxy composites,” Composites Science and Technology, vol. 55, no. 4, pp. 405–412, 1995. View at: Publisher Site | Google Scholar
  58. K. Ohno and M. Ohtsu, “Crack classification in concrete based on acoustic emission,” Construction and Building Materials, vol. 24, no. 12, pp. 2339–2346, 2010. View at: Publisher Site | Google Scholar
  59. J. Zhang, “Investigation of relation between fracture scale and acoustic emission time-frequency parameters in rocks,” Shock and Vibration, vol. 2018, 14 pages, 2018. View at: Publisher Site | Google Scholar
  60. X. Li, Z. Feng, G. Han et al., “Breakdown pressure and fracture surface morphology of hydraulic fracturing in shale with H2O, CO2 and N2,” Geomechanics and Geophysics for Geo-Energy and Geo-Resources, vol. 2, no. 2, pp. 63–76, 2016. View at: Publisher Site | Google Scholar
  61. B. Zhang, B. Ji, and W. Liu, “The study on mechanics of hydraulic fracture propagation direction in shale and numerical simulation,” Geomechanics and Geophysics for Geo-Energy and Geo-Resources, vol. 4, no. 2, pp. 119–127, 2018. View at: Publisher Site | Google Scholar
  62. J. Tan, R. Hu, W. Wang, and J. Dick, “Palynological analysis of the late Ordovician - early Silurian black shales in South China provides new insights for the investigation of pore systems in shale gas reservoirs,” Marine and Petroleum Geology, vol. 116, article 104145, 2020. View at: Publisher Site | Google Scholar
  63. J. He, C. Lin, X. Li, Y. Zhang, and Y. Chen, “Initiation, propagation, closure and morphology of hydraulic fractures in sandstone cores,” Fuel, vol. 208, pp. 65–70, 2017. View at: Publisher Site | Google Scholar

Copyright © 2020 Jingqiang Tan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views288
Downloads326
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.