ThermalHydraulicMechanical (THM) Coupling Behaviour of Fractured Rock Masses
View this Special IssueResearch Article  Open Access
Zhonghu Wu, Huailei Song, Liping Li, Zongqing Zhou, Yujun Zuo, Wenjibin Sun, Hao Liu, Yili Lou, "Study on the Damage Evolution Process and Fractal of QuartzFilled Shale under ThermalMechanical Coupling", Geofluids, vol. 2021, Article ID 8843120, 14 pages, 2021. https://doi.org/10.1155/2021/8843120
Study on the Damage Evolution Process and Fractal of QuartzFilled Shale under ThermalMechanical Coupling
Abstract
Filling of brittle minerals such as quartz is one of the main factors affecting the initiation and propagation of reservoir fractures in shale fracturing, in order to explore the failure mode and thermal damage characteristics of quartzfilled shale under thermalmechanical coupling. Combining the theory of damage mechanics and thermoelasticity, RFPA^{2D}Thermal is used to establish a numerical model that can reflect the damage evolution of shale under thermalsolid coupling, and the compression test under thermalmechanical coupling is performed. The test results show that during the temperature loading process, there is a temperature critical value between 60°C and 75°C. When the temperature is less than the critical temperature, the test piece unit does not appear obvious damage. When the temperature is greater than the critical temperature, the specimen unit will experience obvious thermal damage, and the higher the temperature, the more serious the cracking. Under the thermalmechanical coupling of shale, the tensile strength and elastic modulus of shale show a decreasing trend with the increase of temperature. The failure modes of shale under thermalsolid coupling can be roughly divided into three categories: “V”shaped failure (30°C, 45°C, and 75°C), “M”shaped failure (60°C), and inverted “λ”shaped failure (90°C). The larger the fractal dimension, the more complex the failure mode of the specimen. The maximum fractal dimension is 1.262 when the temperature is 60°C, and the corresponding failure mode is the most complex “M” shape. The fractal dimension is between 1.071 and 1.189, and the corresponding failure mode is “V” shape. The fractal dimension is 1.231, and the corresponding failure mode is inverted “λ” shape.
1. Introduction
With the rapid consumption of conventional energy such as petroleum, unconventional clean energy such as shale gas has become a current research focus [1–5]. Shale gas refers to the natural gas resources that remain in shale reservoirs in the form of adsorption and dissociation during the evolution of geological structures. It has the characteristics of wide distribution, many developed layers, and large storage capacity [6–9]. It is predicted that my country’s recoverable shale gas can reach m^{3}, ranking first in the world [10, 11]. The shale matrix has the characteristics of small pore throats and extremely low permeability. More than 90% of shale gas wells need to undergo shale reservoir fracturing before they can be directly put into production [12–14]. The research on shale gas in my country started late. How to achieve largescale exploitation of shale gas has important strategic significance for the upgrading of my country’s energy structure and green and sustainable development.
The fracturing of shale reservoirs uses highpressure fluid to break the rock mass to form a complex fracture network with high permeability according to Jiang et al. [15]. Because shale gas reservoirs are in an environment of high confining pressure and high temperature, during the fracturing process, shale mining involves the mutual coupling of stress field and temperature field. Therefore, the problem of shale thermal damage is currently a hot topic for scholars at home and abroad. Meng et al. [16] conducted triaxial compression tests on shale samples whose temperature varied between 25 and 120°C to study temperature and shale peak strength, elastic modulus, Poisson’s ratio, cohesion, and the relationship between the internal friction angles. Masri et al. [17] conducted triaxial compression tests on shale at different temperatures and found that the elastic modulus and compressive failure strength of shale decreased significantly as the temperature increased, but the overall deformability of shale increased. Mohamadi et al. [18] conducted a series of triaxial tests on shale at different temperatures and discussed in detail the influence of temperature on shale strength. Jha et al. [19] conducted uniaxial compression tests and tensile strength tests at different temperatures from normal temperature to 900°C and studied in detail the influence of high temperature on the physical and mechanical properties and uniaxial compressive strength of shale. Rybacki et al. [20] conducted creep tests on shale with different temperatures (50~200°C) and different confining pressures (50~200 MPa) to study the influence of temperature and confining pressure on creep characteristics. Wang et al. [21] conducted multistage creep tests on shale samples under different temperatures and confining pressures to study the effects of temperature, confining pressure, and precrack on shale creep characteristics. Guo et al. [22] conducted uniaxial compression tests and acoustic emission tests on shale samples with different bedding surface dip angles at different temperatures and studied the effect of temperature on the physical and mechanical properties of shale. Yang et al. [23] conducted dynamic compression tests on shale specimens in the temperature range of 20~220°C, studied the mechanical properties of shale at different temperatures, and analyzed the dynamic deformation and failure process of shale.
Fractal theory has unique advantages in characterizing complex and irregular objects, and it provides a new theoretical basis for the evolution process of rock failure. Xie [24] combined fractal geometry with rock damage mechanics and proposed that the failure process of rock has fractal characteristics. Acoustic emission is a phenomenon of instantaneous elastic waves caused by the rapid release of local energy of the rock. Acoustic emission can reflect the change of the nature and state of the rock mass to a certain extent. Through the analysis of the change law of acoustic emission parameters, we can understand that the rock mass suffers using the characteristics of selfsimilarity of acoustic emission and further use fractal theory to analyze the failure process of rock mass. Yuan and Li [25] analyzed the fractal dimension of acoustic emission spatial distribution during rock damage and failure and established the relationship between stress and energy release and the fractal dimension of acoustic emission spatial distribution. Xie et al. [26] used a box dimension to calculate the fractal dimension of acoustic emission spatial distribution and studied the relationship between rock damage and acoustic emission fractal dimension. Zhang et al. [27] analyzed the fractal dimension of the spatial distribution of acoustic emission generated by the Brazilian split test and found that the fractal dimension can well reflect the initiation and expansion of rock microcracks. Lou et al. [28] analyzed the spatial distribution characteristics of acoustic emission of different bedding shale under fluidmechanical coupling through numerical simulation experiments. It is found that the fractal dimension of the spatial distribution of acoustic emission can well reflect the influence of bedding on the failure mode of shale.
Shale is a kind of brittle rock formed by transportation and sedimentation in the process of geological structure evolution. The shale is filled with a large number of mineral particles. Under the action of temperature, the shape and content of mineral particles have a significant impact on the macro and micromechanical properties and fracture modes of shale [29–32]. Although there are many studies on the mechanical properties of shale at different temperatures, they mainly focus on the macroscale aspects such as bedding effects, different confining pressures, and creep studies. There are few reports on the micromechanical properties and failure modes of quartzcontaining mineralfilled shale under temperature. Therefore, numerical simulation experiments are carried out on the microscale of shale to study the influence of thermalmechanical coupling on the compressive strength, failure process, and failure mode of shale. The test results have an important reference value for the fracture initiation and propagation mechanism of fracturing reformation in shale reservoirs.
This article uses statistical methods to describe the heterogeneity of quartz and shale matrix. Through the RFPA^{2D}Thermal software, the shale is subjected to numerical simulation tests under constant confining pressure and different temperature conditions, and the compressive strength, elastic modulus and failure process of shale under different temperature conditions are studied in detail. Calculate the fractal dimension of the acoustic emission distribution map and analyze the relationship between the fractal dimension and the failure mode under different temperature conditions. The research results will have an important reference value for fracture mechanism, secondary crack initiation and propagation prediction, and enhanced oil recovery in shale fracturing.
2. Regional Geological Characteristics
Guizhou is one of my country’s oil and gas storage bases. The regional tectonic unit is divided into the Yangtze quasiplatform. Between hilly basins, the terrain is high in the west and low in the east and descends from the center to the north, east, and south. According to the distribution of surface structures, the metamorphic rock forms and structural combination styles of the area and other structural features. The study area has mainly experienced four tectonic cycles during the evolution of geological structure: Wuling tectonic movement, CaledonianXuefeng tectonic movement, Yanshan tectonic movement, and Himalayan tectonic movement. Among them, the Yanshanian period is the most intense tectonic movement, which is the main cause of the current topography. The Himalayan period superimposed and transformed the structures formed in the Yanshan period [33]. Shale gas reservoirs are mainly located in the Lower Cambrian Niutitang Formation, which is also the most widely developed bottom layer in Guizhou, as shown in Figure 1.
The study of rock and mineral characteristics is an important factor that needs to be considered in the exploration and development of shale gas. The higher the content of brittle minerals in shale, the more likely it is to produce cracks under the action of tectonic stress or hydraulic fracturing, providing storage space and seepage channel for shale gas. It can be seen that the shale is filled with a large amount of quartz through the microsection identification. The identification process and results are shown in Figure 2. Generally, the lower the clay mineral content in shale, the higher the brittle mineral content such as quartz, feldspar, and calcite, and the stronger the rock brittleness, the easier it is to form natural cracks and induced cracks under tectonic stress field or artificial fracturing [35, 36]. The shale mineral composition in the study area is complex and diverse, of which clay minerals and brittle minerals are the main components in shale. This article mainly takes the black shale of the Lower Cambrian Niutitang Formation in Fenggang III District as the research object. Carrying out Xray diffraction quantitative analysis and clay mineral analysis, the instrument model is X’pert powder. Xray diffraction analysis is shown in Figure 3. The analysis results show that the quartz content is between 35.79% and 92.49%, with an average content of 62.09%. The content of plagioclase is 028.94%, with an average of 12.32%. The clay mineral content is 029.15%, with an average of 7.64%. It also contains a small amount of minerals such as pyrite, calcite, potash feldspar, and iron dolomite, with an average content of 7.50%, 3.99%, 3.24%, and 2.96%, respectively (Figure 4). As the buried depth increases, the quartz content gradually increases, while the plagioclase and pyrite content gradually decreases. Except for a few parts in the middle, the content of clay minerals is relatively high, and the overall decreases gradually (Figure 5).
Shale reservoirs have the characteristics of low permeability and micronano pore development. This lowpermeability and microporosity feature has an important influence on the content and storage of shale gas. As the main storage space of gas reservoirs, shale micropores determine the enrichment degree of shale gas reservoirs. Generally speaking, the more developed the microfractures and pores of shale, the higher the flow capacity and the richer the gas reservoir. In this paper, the relaxation of nuclear magnetic resonance (NMR) technology is used to characterize the pore structure of shale. The experimental instrument model is MesoMR23060HI. Five groups of samples in the study area are selected for detection and analysis. The results show that the pore throats of the Lower Cambrian shale in Fenggang III block are mainly nanopores, which are distributed between 0 and 0.1 μm. The pores of shale are mainly nanopores, and the pores are concentrated between 0.001 μm0.01 μm and 0.010.4 μm, and there are also a few micronsized pores. The pore size distribution is shown in Figure 6.
3. ThermalMechanical Coupling Model of Shale Microscopic Unit
3.1. Thermal Damage Theory of Shale Microscopic Units
Based on the theory of continuum mechanics, shale is considered to be an ideal linear elastic body, which satisfies the generalized Hooke’s law under external load [37]. A damage model of the microunit is established to characterize the evolution process of shale damage under the thermalmechanical coupling action and its influence on elastic modulus and compressive strength. Because the compressive strength of shale is much greater than its tensile strength, under the action of external load, the maximum tensile stress criterion and MohrCoulomb criterion are used as the basis for judging element damage, and the maximum tensile stress criterion has priority. When the stress state of the shale reaches the critical state of the MohrCoulomb criterion, the unit begins to undergo shear failure. The above two judgment failure criteria are expressed as Equation (1) and Equation (2) [38]:
where and are functions of the maximum tensile stress criterion and MohrCoulomb criterion, respectively; and are the first and third principal stresses, respectively; and are uniaxial tensile strength and uniaxial compression strength, respectively; is the internal friction angle of shale.
3.2. Shale ThermalMechanical Coupling Relationship
This article considers the high temperature and highpressure environment of shale, which involves the coupling relationship between temperature field and stress field. Through the control equations of stress field and temperature field, the damage evolution process of shale under external force and the influence of damage on stress field and temperature field are considered. Under the coupling effect of temperature field and stress field, microcracks in shale begin to sprout, and the damage of shale in turn affects the elastic modulus, compressive strength, and thermal conductivity of shale. Considering the effect of thermal stress on shale deformation caused by the temperature field, the stressstrain relationship of shale can be expressed as Equation (3) [39]:
where is the stress tensor; is the strain tensor; is the shear modulus of the shale, ; is the volumetric strain; is the Kronecker function; is the bulk modulus of shale; is the thermal stress coefficient of shale, , where and are Lame constant and thermal expansion coefficient, respectively; is the amount of temperature change.
The thermal conductivity of shale is closely related to temperature and reflects the unevenness of shale temperature field. Regarding the thermophysical properties of the shale microscopic unit, when the damaged unit is not initiated, the thermal conductivity of the shale is constant. The heat transfer coefficient of shale after element tensile damage is calculated as Equation (4) [40]. The specific heat capacity of shale is also affected by temperature. The specific heat capacity of the unit under tension of the shale is calculated according to Equation (5) [40]:
In the formula, is the thermal conductivity at a temperature of 0°C; is the specific heat capacity of the shale at a temperature of 0°C. is the impact factor, which is generally .
In the numerical model established in this paper, the mutual conversion between thermal energy and mechanical energy is ignored, and the energy conservation equation is calculated according to Equation (6) [41]:
where is the initial temperature of the shale; is the density of the shale; is the specific heat capacity of the shale; is the thermal conductivity of the shale.
4. Numerical Modeling
RFPA^{2D}Thermal is based on finite element theory and statistical damage theory. Considering the inhomogeneity of the rock, we simplified the complicated macrononlinear problem into a finetomicrolinear problem and combined the assumption of random distribution of the inhomogeneity with numerical calculation methods. The numerical simulation of the nonuniform rock failure process can be realized [42]. The shale model is discretized into a numerical model composed of mesoprimitives. It is assumed that the mechanical properties of the discretized microprimitives obey the Weibull statistical distribution law, thereby establishing the relationship between the mechanical properties of the micro and macromedium [43]:
In the formula, is the statistical distribution density of the mechanical property of shale primitives. is the mechanical property parameter of the shale element; is the average value of the mechanical property parameter of the shale element; is the property parameter of the distribution function, which reflects the uniformity of shale and is defined as the uniformity coefficient of shale.
In order to study the influence of temperature on the mechanical properties and failure modes of shale, this paper uses the RFPA^{2D}Thermal software to establish a numerical model and carries out simulation tests at different temperatures. In the numerical model, the larger the elastic modulus of the mineral particles, the brighter the color, which can be used to characterize the quartz mineral particles and the shale matrix. According to the physical experiments, the mechanical parameters and thermodynamic parameters of the quartz mineral and shale matrix in the model are shown in Table 1 [44].

This experiment establishes a numerical model of quartzfilled shale at different temperatures. The quartz content in the model is 62.09%. Shale pores are mostly nanosized pores and a small amount of microsized pores. This test does not consider the influence of shale primary pores. The loading model is shown in Figure 7. During the loading process, a fixed confining pressure of 10 MPa is applied first, and then, the temperature is loaded. When the temperature reaches 30°C, 45°C, 60°C, 75°C, and 90°C, the heating is stopped. The temperature is loaded step by step, and the increment is 5°C/step. After the temperature reaches a certain value, displacement loading is then carried out, and the displacement loading rate is 0.0005 mm/step.
5. Result Analysis
5.1. Thermal Damage Evolution Characteristics of Shale
Shale is often filled with brittle minerals such as quartz. Because the thermal expansion coefficients of brittle mineral particles and shale matrix are different, the thermal expansion at the boundary between quartz and shale matrix is inconsistent, resulting in tensile or compressive thermal stress in the junction between quartz particles and shale matrix. Figure 8 is the damage distribution diagram and the corresponding acoustic emission diagram of the shale specimen at different temperatures. In the acoustic emission diagram, the yellow represents the tensile damage caused by the current step, the red represents the shear damage caused by the current step, and the black represents the damage. During the temperature loading process, there is a critical value of unit thermal damage temperature between 60°C and 75°C. When the temperature is less than the critical value, the thermal stress does not cause obvious damage to the shale specimen. When the temperature is greater than the critical value, the local thermal stress is greater than the tensile strength of the shale specimen, and the specimen appears obvious microcracking. It can be seen from the acoustic emission diagram that the fracture is mainly tensile failure, and a small amount of shear failure occurs at 90°C. The higher the temperature, the more serious the damage, and the microfractures are mainly distributed at the cementation between quartz and shale matrix. This is because quartz particles are randomly distributed in shale reservoirs, and the thermal expansion at the boundary between quartz and shale matrix is uncoordinated.
Acoustic emission number is the elastic wave signal of element damage initiation and release of shale specimen under thermal stress, reflecting the damage evolution process of specimen. Figure 9 is the change rule of acoustic emission number with temperature. During the temperature loading process, as the temperature increases, the number of acoustic emissions increases rapidly, and the corresponding thermal damage becomes more severe. The thermal damage evolution process of shale is divided into not damage stage, microdamage stage, and damage stage. When the temperature is less than the critical temperature, the acoustic emission number is basically 0, and there is no obvious damage to the shale sample. When the temperature reaches the critical temperature, the number of acoustic emissions and the cumulative number of acoustic emissions begin to appear, and the shale sample begins to appear microdamage. When the temperature reaches 80°C, the thermal stress increases rapidly, and the thermal damage of the shale sample develops rapidly. This is consistent with the abovementioned damage evolution process. The higher the temperature, the greater the number of acoustic emissions and the faster the cumulative number of acoustic emissions increases, indicating that the higher the temperature, the more severe the thermal damage of the shale specimen.
5.2. The Change Law of Shale Strength under ThermalMechanical Coupling
During the displacement loading process, the thermalmechanical coupling effect of shale has a significant effect on the compressive strength and elastic modulus of shale filled with quartz minerals. Table 2 shows the strength and deformation parameters of shale obtained in the experiment, and Figure 10 shows the change trend graph of shale compressive strength and elastic modulus with temperature. With the increase of temperature, the tensile strength of shale shows a decreasing trend. The compressive strength is the largest at 30°C, which is 79.79 MPa, and at 90°C, the minimum is 10.15 MPa. The elastic modulus of shale shows a decreasing trend with the increase of temperature, and the most obvious change is between 60°C and 75°C. The maximum elastic modulus is 84.9 GPa at 30°C, and the smallest at 90°C is 54.88 GPa. The thermal expansion of the quartz mineral particles and the shale matrix is uneven. The higher the temperature, the greater the thermal stress generated, and the mechanical properties of shale will weaken. When the temperature is greater than the critical temperature, the thermal stress is greater than the maximum tensile strength of the shale microscopic unit, and thermal damage occurs in the sample, which leads to increased nonuniformity of the shale. The higher the temperature, the more serious the microfracture and the lower the bearing capacity of shale in the displacement loading stage.

5.3. The Influence of ThermalSolid Coupling on the Failure Process of Shale
Figure 11 shows the damage and failure process of shale under thermalmechanical coupling and the corresponding acoustic emission diagram. It can be seen from the figure that different temperature effects have significant effects on the damage evolution and failure modes of shale, and the failure modes can be roughly divided into three categories. It is a “V”shape failure when the temperature is 30°C, 45°C, and 75°C and “M”shape failure when the temperature is 60°C. When the temperature is 90°C, it is an inverted “λ”shape failure. When the temperature is 30°C, cracks will initiate at the lower left corner of the specimen first. After the cracks appear, stress concentration will easily occur at both ends of the cracks, causing the cracks to rapidly propagate to both ends, and finally, a “V”shaped crack is formed. When the temperature is 45°C, cracks appear in the upperright corner of the specimen under the combined action of thermal stress and loading stress. With continuous loading, an inverted “V”shaped crack is finally formed. When the temperature is 60°C, the crack network during the instability and failure of the specimen is the most complicated, and the two cracks expand downward at the same time, forming an “M”shaped fracture zone. When the temperature is 75°C, microcracks first appear in the upper right corner, and finally form an inverted “V”shaped fracture zone. When the temperature is 90°C, the microcracks first appeared at the left end of the specimen and expanded along a straight line; then, the rightside microcracks appeared, quickly expanded to the left, and finally converged to form an inverted “λ”shaped fracture zone. Under the action of thermalmechanical coupling, the higher the temperature, the rougher the cracks. The crack rupture zone at low temperature is smaller. The damage at high temperature mainly appears in the form of rupture zone, showing obvious plastic deformation. The main reasons for the above phenomenon are as follows: on the one hand, the uneven arrangement of quartz particles in the shale deposition process. The thermal expansion of quartz and shale matrix is quite different, resulting in unstable thermal stress under the action of temperature. As a result, thermal damage occurs at the boundary between the quartz and the shale matrix, and the bearing capacity of the shale specimen is reduced during displacement loading. On the other hand, the quartz mineral particles and crystals inside the rock undergo violent thermal movement under the action of temperature, which is more prone to microdamage. The strength of the quartz particles is greater than that of the shale matrix, the cracks are prone to bend when they pass through the quartz, and the damage range is larger under the action of external load.
In the acoustic emission diagram, yellow represents the tensile failure initiated at the current step, red represents the shear damage caused by the current step, and black represents the damaged unit. Looking at the acoustic emission diagram, it can be seen that the shale specimens are mainly tensileshear failure. The cracks when the specimen is instability and failure are connected by a tensile damage unit and a compressionshear damage unit. The cumulative damage development trend of microunits directly reflects the macrofailure mode of shale specimens. The joint action of thermal stress and external stress will cause tensile damage when the tensile strength of the specimen is first reached. As the loading process reaches the critical value of the MohrCoulomb strength criterion, the specimen unit appears compressionshear failure. The above results indicate that the heterogeneous microstructure of the quartz filled in shale under thermalmechanical coupling has a significant impact on its damage evolution process and failure mode.
5.4. Fractal Characteristics of Shale Damage Evolution Process under ThermalSolid Coupling
Figure 12 is a graph of the relationship between stress level and AE energy under thermalsolid coupling. It can be seen from the figure that temperature has a significant effect on the release of AE energy, and the law of energy release under the same stress level at different temperatures is significantly different. Due to the effect of thermal stress, at 90°C, the shale specimen begins to microcrack when the stress level is 10%, resulting in energy release. Except for the stress level at 100%, the released energy is slightly lower than 60°C, and the rest of the stress levels of the energy release is maximum. When the temperature is low (30°C and 45°C), the AE energy release increases with the increase of temperature, and it is the smallest at 30°C. When the temperature is high (60°C, 75°C, and 90°C), the release of AE energy at peak stress is greater. The AE energy release at 60°C is the largest, followed by 90°C, and the minimum is 75°C. When the temperature is lower than the critical temperature of damage, there is almost no microfracture in the shale specimen under the action of thermal stress. Therefore, energy release begins when the stress level reaches 60% at low temperature. When the temperature is greater than the critical temperature of damage, due to the different thermal damages at different temperatures, the temperature is higher, and there is an earlier release of AE energy. At 75°C, the stress level reaches 40% and the microcrack appears and starts to release energy.
Fractal theory can quantitatively describe complex objects in the world. It is widely used in the field of rock failure and helps to reveal the law of damage and fracture of rocks. In this paper, the fractal dimension of the image is calculated by the box dimension, and the fractal dimension program is written on the MATLAB software platform. Binarize the acoustic emission diagrams under different stress levels obtained in the experiment, and import the calculation program to obtain the fractal dimension of the acoustic emission diagrams. The formula is defined as follows:
where is the selfsimilar fractal dimension of the damage area; is the side length of the square box; is the number of boxes required to cover the damaged area in the entire graph with a square box with side length .
Table 3 shows the AE energy and fractal dimension values. Figure 13 is a graph showing the change trend of fractal dimension with stress level under thermalmechanical coupling. It can be concluded from the figure that under the effect of thermalsolid coupling, the fractal dimension increases with the increase of the stress level. Due to the incoordination of the thermal expansion coefficient between the quartz particles and the shale matrix during the temperature loading stage, thermal stress is generated, which leads to the initiation of thermal cracking of the specimen. The higher the temperature, the more serious the thermal cracking. When the temperature is less than the critical temperature of damage, the fractal dimension begins to appear after the stress level reaches 10%, and the lower the temperature, the later it appears. When the temperature is greater than the critical temperature of damage, the fractal dimension begins to appear at the beginning of displacement loading. The fractal dimension is the lowest at 30°C, which is 1.071. The fractal dimension is the largest at 60°C, which is 1.262, and the corresponding failure mode is the most complicated “M” shape. The fractal dimension is 1.071~1.189, and the corresponding failure mode is “V” shape. The fractal dimension is 1.231, and the corresponding failure mode is inverted “λ.” The failure mode of the specimen is quantitatively analyzed by the fractal dimension. The larger the fractal dimension at peak stress, the more complex the failure mode, the rougher the crack, and the more severe the internal damage. The abovementioned content is depicted in a threedimensional coordinate system. As shown in Figure 14, we can clearly see the changing trend of stress level, temperature, and fractal dimensions. From the perspective of fractal, the influence of thermalmechanical coupling on the failure mode of shale is confirmed.

6. Conclusion
Shale is often filled with quartz minerals, and the difference in thermal expansion between quartz particles and shale matrix is an important influencing factors of rock fracture under thermalmechanical coupling. In this paper, the thermalsolid coupling model is established to study the influence of temperature on the mechanical properties and damage evolution process of quartzfilled shale, and the following rules are summarized:
(1) The quartz content of the Niutitang Formation shale in the III District of Fenggang is 35.79%~92.49%, with an average content of 62.09%. As the depth of burial increases, the overall quartz content gradually increases. Shale pore throats are mainly nanopores distributed between 0 and 0.1 μm. Shale pores are concentrated between 0.001 μm~0.01 μm and 0.01~0.4 μm, mainly nanopores, with a few micronsized pores
(2) The effect of temperature has a significant effect on the thermal damage of shale filled with a large amount of quartz. The evolution process of thermal damage can be divided into nondamage stage, microdamage stage, and damage stage. There is a critical temperature value. When the temperature is greater than the critical temperature, the specimen unit will experience obvious thermal cracking. The higher the temperature, the more serious the cracking. The thermalmechanical coupling effect has a significant effect on the compressive strength and elastic modulus of shale filled with quartz minerals. With the increase of temperature, the tensile strength and elastic modulus of shale show a decreasing trend
(3) Temperature has a significant influence on the damage evolution and failure modes of shale under thermosolid coupling. The failure modes can be roughly divided into three categories. When the temperature is 30°C, 45°C, and 75°C, it is a “V”shape failure; when the temperature is 60°C, it is an “M”shape failure; when the temperature is 90°C, it is an inverted “λ”shape failure. The higher the temperature, the rougher the fracture will be. The damage at high temperature mainly appears in the form of rupture zone with obvious plastic deformation
(4) Using the selfsimilarity of the spatial distribution of acoustic emission points, we quantitatively analyze the influence of thermalmechanical coupling on the failure mode of shale specimens based on fractal theory. The larger the fractal dimension, the more complicated the failure mode of the specimen and the more serious the internal damage. When the temperature is 60°C, the fractal dimension is the largest, which is 1.262, and the corresponding failure mode is the most complicated “M” shape. The fractal dimension is between 1.071 and 1.189, and the corresponding failure mode is “V” shape. The fractal dimension is 1.231, and the corresponding failure mode is inverted “λ shape”.
Data Availability
The data used to support the study is available within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the Talent Introduction Project of Special Fund for Science and Technology of Water Resources Department of Guizhou Province (Project No. KT201804), Guizhou Postgraduate Innovation Fund (Project No. YJSCXJH[2020]087), Guizhou Science and Technology Fund (Project No. [2020]4Y046, Project No. [2019]1075, and Project No. [2018]1107), the National Natural Science Foundation of China (Project Nos. 51964007 and 51774101), and Project Scientific Research Project of Guiyang Rail Transit Line 2 Phase I Project (Project No. D2(I)FWYJ2019001WT). This study is also funded by the Teaching reform project of Guizhou University (Project No. JG201990) and Guizhou Province Mine Dynamic Disaster Early Warning and Control Technology Innovation Talent Team Project (Project No. [2019]5619).
References
 D. Zhang, P. G. Ranjith, and M. S. A. Perera, “The brittleness indices used in rock mechanics and their application in shale hydraulic fracturing: a review,” Journal of Petroleum Science and Engineering, vol. 143, pp. 158–170, 2016. View at: Publisher Site  Google Scholar
 C. Yan, J. Deng, Y. Cheng, M. Li, Y. Feng, and X. Li, “Mechanical properties of gas shale during drilling operations,” Rock Mechanics and Rock Engineering, vol. 50, no. 7, pp. 1753–1765, 2017. View at: Publisher Site  Google Scholar
 G. Li, G. Li, Y. Wang, S. Qi, and J. Yang, “A rock physics model for estimating elastic properties of upper Ordovician lower Silurian mudrocks in the Sichuan Basin, China,” Engineering Geology, vol. 266, article 105460, 2020. View at: Publisher Site  Google Scholar
 L. Li, B. Huang, Y. Li, R. Hu, and X. Li, “Multiscale modeling of shale laminas and fracture networks in the Yanchang formation, Southern Ordos Basin, China,” Engineering Geology, vol. 243, pp. 231–240, 2018. View at: Publisher Site  Google Scholar
 Z. Shen, L. Zhou, X. Su, H. Li, J. Tang, and X. Zheng, “An experimental investigation of the nonlinear gas flow and stressdependent permeability of shale fractures,” Energy Science and Engineering, vol. 8, no. 8, pp. 2808–2822, 2020. View at: Publisher Site  Google Scholar
 R. Wang, Y. Gu, W. Ding et al., “Characteristics and dominant controlling factors of organicrich marine shales with high thermal maturity: a case study of the Lower Cambrian Niutitang Formation in the Cen’gong block, southern China,” Journal of Natural Gas Science and Engineering, vol. 33, pp. 81–96, 2016. View at: Publisher Site  Google Scholar
 Z. H. Wu, Y. J. Zuo, S. Y. Wang et al., “Numerical study of multiperiod palaeotectonic stress fields in Lower Cambrian shale reservoirs and the prediction of fractures distribution: a case study of the Niutitang Formation in Feng’gang No. 3 block, South China,” Marine and Petroleum Geology, vol. 80, pp. 369–381, 2017. View at: Publisher Site  Google Scholar
 S. Yin, R. Xie, Z. Wu, J. Liu, and W. Ding, “In situ stress heterogeneity in a highly developed strikeslip fault zone and its effect on the distribution of tight gases: a 3D finite element simulation study,” Marine and Petroleum Geology, vol. 99, pp. 75–91, 2019. View at: Publisher Site  Google Scholar
 G. A. Yakaboylu, N. Gupta, E. M. Sabolsky, and B. Mishra, “Mineralogical characterization and strain analysis of the Marcellus shales,” International Journal of Rock Mechanics and Mining Sciences, vol. 130, article 104345, 2020. View at: Publisher Site  Google Scholar
 C. Zou, R. Zhu, Z.Q. Chen et al., “Organicmatterrich shales of China,” EarthScience Reviews, vol. 189, pp. 51–78, 2019. View at: Publisher Site  Google Scholar
 L. Wang, Y. Guo, C. Yang, J. Xiao, C. Lu, and Y. Song, “Mechanical characterization of continental shale in Sichuan Basin of China and its potential impact on reservoir stimulation,” Journal of Natural Gas Science and Engineering, vol. 79, article 103346, 2020. View at: Publisher Site  Google Scholar
 R. WANG, Z. HU, S. LONG et al., “Differential characteristics of the Upper OrdovicianLower Silurian WufengLongmaxi shale reservoir and its implications for exploration and development of shale gas in/around the Sichuan Basin,” Acta Geologica SinicaEnglish Edition, vol. 93, no. 3, pp. 520–535, 2019. View at: Publisher Site  Google Scholar
 Z. Wu, Y. Zuo, S. Wang, J. Sunwen, and L. Liu, “Experimental study on the stress sensitivity and influence factors of shale under varying stress,” Shock and Vibration, vol. 2018, Article ID 3616942, 9 pages, 2018. View at: Publisher Site  Google Scholar
 A. Michael, J. E. Olson, and M. T. Balhoff, “Orientation prediction of fracture initiation from perforated horizontal wells: application in shale reservoirs,” Journal of Petroleum Science and Engineering, vol. 193, article 107355, 2020. View at: Publisher Site  Google Scholar
 G. Jiang, J. Zuo, Y. Li, and X. Wei, “Experimental investigation on mechanical and acoustic parameters of different depth shale under the effect of confining pressure,” Rock Mechanics and Rock Engineering, vol. 52, no. 11, pp. 4273–4286, 2019. View at: Publisher Site  Google Scholar
 L. B. Meng, T. B. Li, L. W. Jiang, and H. M. Ma, “Experimental study on the influence of temperature on shale mechanical properties under conventional triaxial compression,” Advanced Materials Research, vol. 250253, pp. 1452–1455, 2011. View at: Publisher Site  Google Scholar
 M. Masri, M. Sibai, J. F. Shao, and M. Mainguy, “Experimental investigation of the effect of temperature on the mechanical behavior of Tournemire shale,” International Journal of Rock Mechanics and Mining Sciences, vol. 70, pp. 185–191, 2014. View at: Publisher Site  Google Scholar
 M. Mohamadi and R. G. Wan, “Strength and postpeak response of Colorado shale at high pressure and temperature,” International Journal of Rock Mechanics and Mining Sciences, vol. 84, pp. 34–46, 2016. View at: Publisher Site  Google Scholar
 M. K. Jha, A. K. Verma, P. K. Gautam, and A. Negi, “Study of mechanical properties of Vindhayan shaly rocks at elevated temperature,” Journal of the Geological Society of India, vol. 90, no. 3, pp. 267–272, 2017. View at: Publisher Site  Google Scholar
 E. Rybacki, J. Herrmann, R. Wirth, and G. Dresen, “Creep of Posidonia shale at elevated pressure and temperature,” Rock Mechanics and Rock Engineering, vol. 50, no. 12, pp. 3121–3140, 2017. View at: Publisher Site  Google Scholar
 Y. Wang, H. Wang, and X. Shi, “Creep investigation on shalelike material with preexisting fissure under coupling temperatures and confining pressures,” Advances in Civil Engineering, vol. 2019, Article ID 7861305, 10 pages, 2019. View at: Publisher Site  Google Scholar
 Y. Guo, C. Yang, L. Wang, and F. Xu, “Effects of cyclic loading on the mechanical properties of mature bedding shale,” Advances in Civil Engineering, vol. 2018, Article ID 8985973, 9 pages, 2018. View at: Publisher Site  Google Scholar
 G. Yang, J. Liu, X. Li, and J. Bi, “Effect of temperature on shale strength under dynamic impact loading,” Arabian Journal of Geosciences, vol. 13, no. 12, 2020. View at: Publisher Site  Google Scholar
 H. Xie, An Introduction to FractalRock Mechanics, Science Press, Beijing, 1996.
 R.f. Yuana and Y.h. Li, “Fractal analysis on the spatial distribution of acoustic emission in the failure process of rock specimens,” International Journal of Minerals, Metallurgy and Materials, vol. 16, no. 1, pp. 19–24, 2009. View at: Publisher Site  Google Scholar
 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
 S. W. Zhang, K. Shou, X. Xian, J. Zhou, and G. Liu, “Fractal characteristics and acoustic emission of anisotropic shale in Brazilian tests,” Tunnelling and Underground Space Technology, vol. 71, pp. 298–308, 2018. View at: Publisher Site  Google Scholar
 Y. Lou, Z. Wu, W. Sun et al., “Study on failure models and fractal characteristics of shale under seepagestress coupling,” Energy Science & Engineering, vol. 8, no. 5, pp. 1634–1649, 2020. View at: Publisher Site  Google Scholar
 Y. Zhao, L. Zhang, J. Liao, W. Wang, Q. Liu, and L. Tang, “Experimental study of fracture toughness and subcritical crack growth of three rocks under different environments,” International Journal of Geomechanics, vol. 20, no. 8, article 04020128, 2020. View at: Publisher Site  Google Scholar
 F. Sun, Y. Yao, G. Li, and W. Liu, “Simulation of real gas mixture transport through aqueous nanopores during the depressurization process considering stress sensitivity,” Journal of Petroleum Science and Engineering, vol. 178, pp. 829–837, 2019. View at: Publisher Site  Google Scholar
 W. Zhang, Q. Sun, Y. Zhang, L. Xue, and F. Kong, “Porosity and wave velocity evolution of granite after hightemperature treatment: a review,” Environmental Earth Sciences, vol. 77, no. 9, p. 350, 2018. View at: Publisher Site  Google Scholar
 Z. Baoyun, L. Dongyan, L. Ziyun, H. Wei, D. Qian, and L. Hang, “Mechanical behavior of shale rock under uniaxial cyclic loading and unloading condition,” Advances in Civil Engineering, vol. 2018, 2018. View at: Publisher Site  Google Scholar
 J. Liu, W. Ding, R. Wang et al., “Quartz types in shale and their effect on geomechanical properties: an example from the lower Cambrian Niutitang Formation in the Cen’gong block, South China,” Applied Clay Science, vol. 163, pp. 100–107, 2018. View at: Publisher Site  Google Scholar
 R. Wang, H. U. Zongquan, H. Nie et al., “Comparative analysis and discussion of shale reservoir characteristics in the WufengLongmaxi and Niutitang formations: a case study of the well JY1 in SE Sichuan Basin and well TX1 in SE Guizhou area,” Petroleum Geology & Experiment, vol. 40, no. 5, 2018. View at: Google Scholar
 Y. Zhao, Y. Wang, W. Wang, L. Tang, Q. Liu, and G. Cheng, “Modeling of rheological fracture behavior of rock cracks subjected to hydraulic pressure and far field stresses,” Theoretical and Applied Fracture Mechanics, vol. 101, 2019. View at: Publisher Site  Google Scholar
 Y. Zhao, C. L. Wang, and J. Bi, “Analysis of fractured rock permeability evolution under unloading conditions by the model of elastoplastic contact between rough surfaces,” Rock Mechanics and Rock Engineering, vol. 53, no. 12, pp. 5795–5808, 2020. View at: Publisher Site  Google Scholar
 Y. Zhao, L. Zhang, W. Wang, Q. Liu, L. Tang, and G. Cheng, “Experimental study on shear behavior and a revised shear strength model for infilled rock joints,” International Journal of Geomechanics, vol. 20, no. 9, article 04020141, 2020. View at: Publisher Site  Google Scholar
 C. H. Wei, W. C. Zhu, Q. L. Yu, T. Xu, and S. Jeon, “Numerical simulation of excavation damaged zone under coupled thermal–mechanical conditions with varying mechanical parameters,” International Journal of Rock Mechanics and Mining Sciences, vol. 75, pp. 169–181, 2015. View at: Publisher Site  Google Scholar
 F. Gao, Y. Xue, Y. Gao, Z. Zhang, T. Teng, and X. Liang, “Fully coupled thermohydromechanical model for extraction of coal seam gas with slotted boreholes,” Journal of Natural Gas Science and Engineering, vol. 31, pp. 226–235, 2016. View at: Publisher Site  Google Scholar
 Q. L. Yu, C. Zheng, T. H. Yang, S. B. Tang, P. T. Wang, and C. A. Tang, “Mesostructure characterization based on coupled thermalmechanical model for rock failure process and applications,” Chinese Journal of Rock Mechanics and Engineering, vol. 31, no. 1, pp. 42–51, 2012. View at: Google Scholar
 C. W. Zhu, C. H. Wei, J. Liu, H. Y. Qu, and D. Elsworth, “A model of coal–gas interaction under variable temperatures,” International Journal of Coal Geology, vol. 86, no. 23, pp. 213–221, 2011. View at: Publisher Site  Google Scholar
 H. Xin, T. S. Bin, and T. C. An, “Boundaryeffect and scaleheating rate equivalence effect of cracking behavior in the rock models subjected to heating from the central borehole,” Science China Technological Sciences, vol. 63, no. 5, pp. 809–818, 2020. View at: Publisher Site  Google Scholar
 C. A. Tang, H. Y. Liu, W. C. Zhu et al., “Numerical approach to particle breakage under different loading conditions,” Powder Technology, vol. 143144, pp. 130–143, 2004. View at: Publisher Site  Google Scholar
 Z. Wu, Y. Lou, S. Yin et al., “Acoustic and fractal analyses of the mechanical properties and fracture modes of beddingcontaining shale under different seepage pressures,” Energy Science & Engineering, vol. 8, no. 10, pp. 3638–3656, 2020. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 Zhonghu Wu 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.