Study on the Failure Mechanism of Lower Cambrian Shale under Different Bedding Dips with Thermosolid Coupling

To investigate the damage pattern and acoustic emission pattern of temperature on laminated shales, numerical experiments were carried out using the RFPA2D-Thermal numerical software under the effect of thermosolid coupling. During the tests, temperatures of 30°C, 60°C, and 90°C were controlled, and five sets of shales containing different laminar dips were numerically modeled at each temperature, with dips of 0°, 22.5°, 45°, 67.5°, and 90°. The test results show that (1) the increase in temperature reduced the linear elastic phase of the shale specimens in each group, with a significant reduction in the linear elastic phase of the shale at lamina dips of 22.5° and 45°. (2) The lamination effect decreased slightly as the temperature rose from 30°C to 60°C, and the most significant enhancement of the lamination effect on the shale occurred when the temperature reached 90°C. (3) The shale damage pattern is divided into five types (N, ʌ, v, slanted I-type, and cluttered-type), in which the lamina effect is stronger for high-angle lamina dips, and the lamina surface has a strong dominant effect on the entire shale crack expansion. At a temperature of 90°C, the lamina effect and temperature effect of the shale reached their maximum at the same time, and the thermal and load stresses inside the shale acted together causing the shale to show a complex damage mode. (4) The fractal dimension was used to analyze the damage pattern of the shale. The larger the fractal dimension was, the greater the crack rate of the specimen. The fractal dimension curve was flatter at a temperature of 60°C, while at 90°C, the fractal dimension rose rapidly, indicating the most favorable crack expansion in the shale at a temperature of 90°C.


Introduction
As the first country to develop shale gas, the United States has succeeded in reducing natural gas consumption and has seen a gradual increase in natural gas exports, contributing to the optimization of the US energy system [1][2][3][4][5]. In recent years, with the massive depletion of conventional energy sources, there has been an urgent need for a new qenergy source that can replace natural gas and oil, which has led to strong interest in shale gas from all corners of the world [6][7][8][9][10][11][12][13].
Currently, the commonly used technology in shale gas exploitation is hydraulic fracturing technology. In hydraulic fracturing technology, a series of rock mechanic problems are involved [14][15][16][17][18][19][20]. Hence, exploring the mechanical ability of shale is of outstanding importance for the evolution of reformed shale gas reservoirs. Shale gas is stockpiled in shale pores containing shale gas under high ambient stress and high-temperature conditions. Therefore, the exploitation of shale gas always includes the coupling of stress and temperature fields [21][22][23].
Inhomogeneous bedding has attracted the attention of scholars for a long time, but the temperature impact on the mechanical properties of shale rocks has only gradually received attention since the 1970s. Therefore, domestic and foreign researchers have acquired successful outcomes regarding the effect of temperature on the mechanical properties of shale rocks. Lou et al. [24] used numerical modeling experiments to probe the bedding effects of shale with different bedding dips and applied fractal dimensions to analyze the crack laws of shale. Wei et al. [25] created a thermal damage thermosolid coupling model, performed physical tests on granites, and found that the thermal effect caused changes in rock parameters thereby stimulating rock tensile-reduction damage. Zhou et al. [26] established compression, and the thermosolid coupling model of thermal expansion provided the thermosolid coupling equation for shale. Alm et al. [27] studied the effect of microcrack density on the elastic and mechanical properties of granite. After the granite was heated, the crack density of the rock was measured by scanning electron microscopy and other instruments, and finally, it was found that when the temperature exceeded 300°C, the crack density increased significantly. Lu et al. [28] carried out a large-scale triaxial compression test on shale. By controlling temperature variables, the relationship between temperature, peak strength of shale, Poisson's ratio, internal elastic modulus, and friction angle was found. Wang et al. [29] found that the fractality, connectivity and connectivity of parallel laminations were 1.31, 1.12, and 1.61 times higher than those of vertical ones by studying the coal rocks in the Fukang Mine area. Yao et al. [30] investigated the mechanical properties of outcrop shale under different laminar dip angles and found that the damage extension direction almost coincided with the laminar dip angle around 0~45°, and when the dip angle was 60~90°, only a small part of the fracture extension direction coincided with the laminar dip angle. Only a small part of the fracture extension direction coincides with the laminar dip angle.. Masri et al. [31] performed hydrostatic and conventional triaxial tests to explore the mechanical properties of the bedding angles of parallel and vertical bedding angles at different temperatures. As the temperature rose, the elastic modulus and compressive strength of the rock were weakened, but the overall deformation was strengthened. Gautam et al. [32] studied the deformation characteristics of sandstone at various temperatures. They observed that the stress-strain curve showed that sandstone underwent brittle deformation at 250°C, and plastic deformation occurred when the temperature exceeded 450°C. Ranjith et al. [33] performed uniaxial compression tests on Hawkesbury sandstone at different temperatures. The results showed that the maximum elastic modulus of Hawkesbury sandstone was a threshold value at 500°C, which explained the softening phenomenon of sandstone when the temperature exceeded 500°C. Lei et al. [34] used the MTS 815 MPa rock testing machine to monitor the microscopic failure process of sandstone during uniaxial compression testing. AE monitored the sandstone transition from brittleness to plasticity at approximately 600°C and found that the sandstone microscopic porosity decreased with increasing temperature. Mahanta et al. [35] conducted experiments on Manoharpur sandstone, Dholpur dolomite, and Bellary sandstone during these heating treatments. Their fracture strength increased by 40%, 25%, and 65%, respectively, within the range of 100°C. However, with temperature rising to 600°C, the fracture strength decreased by 59%, 36%, and 30%, respectively. Yang et al. [36] performed a dynamic compression test on shale rocks and controlled the temperature between 20 and 220°C. Their conclusion was that there was a temperature threshold between 20 and 220°C, before which the compressive strength increased with temperature and beyond which the compressive strength decreased with increasing temperature. Currently, many achievements have been made regarding the mechanical properties of rocks at diverse temperatures, dip angles, and confining pressures. They have mainly focused on the rupture modes of rocks at various temperatures. However, there have been few reports on the failure processes and acoustic emission characteristics under thermosolid coupling at different dip angles of shale bedding and thermosettings. Therefore, it is of great practical significance to carry out numerical simulation experiments to explore the mechanisms of thermosolid coupling on the propagation of shale fracturing and restructuring.
This paper utilizes the Lower Cambrian Niutitang Formation shale as the research object in the northern Guizhou area, uses statistical methods to describe the correlation between the bedding dip and temperature in the shale, and uses RFPA2D-Thermal to establish a numerical model of 5 groups of shale samples that have different dips. A constant confining pressure and several temperature conditions are set to simulate the natural environment of shale under high temperature and high pressure. Simultaneously, shale is analyzed in detail, in which the shale contains compressive strength and cracking processes at various temperatures, and the evolution of acoustic emission signals during the process is studied. Most areas of Guizhou are within the tectonic unit of the Yangtze block, which is divided into four blocks. The third block of shale gas in Feng'gang, Guizhou (referred to as the research area in the text), is an important part of the national shale gas resource experiment, with abundant resources and a good geological environment. The geographic location of the study area is located in the southeastern part of Zunyi city in the northern part of Guizhou Province (see Figure 1). The administrative divisions are located in the subregions of Meitan County, Feng'gang County, and Sinan County of Tongren city and parts of Meitan-Feng'gang-Tongren Sinan County from Zhengjiashan in Sinan County to Zhangjiapo and Tianjiawan Line in the east, Liujiawuji-Liangshuijing-Huangnipo Line in Meitan County in the west, Jiangjiagou-Ranjiayuanzi-Houtang Line in Tan County, Meitan County in the south, and North to Feng'gang Xianglu Mountain-Qinggang Park-Loquat Bay line. The area coverage in the study area is approximately 1167.49 km 2 . The study area has diverse terrain types that contain hills, midmountains, hills and dams interspersed with each other, and karst landforms with conventional landforms alternately spread in between. The terrain is generally mountainous and hilly. The western terrain is dominated by flat dams and hills. Most areas contain 2 Geofluids low hills and relatively flat valleys, with certain undulations. Some areas are cut by rivers (100~300 m), and the east is mountainous. The terrain is undulating, and the cut is deeper.

Influence of Folds and Fractures on Shale Temperature.
The occurrence of folds and fractures caused the strata to no longer lie on the same plane during the same period, presenting a spatially scattered structure with a large depth decrease. Affected by structural compression and deformation, the Niutitang Formation in the study area is generally shallower in the southeast and northwest (approximately 1500 m and 500 m, respectively) and relatively deep in the middle (mainly 2000~4000 m); among them, the overall burial depth is between 1000 m and 3500 m, and the Meitan compound anticline is relatively shallow, with a burial depth of 500 m to 1500 m, while the Yachuan compound syncline contains Triassic strata, and the burial depth is at a maximum of 5000 m (see Figure 2). Near the study area, detailed logging data record the burial depth of the Niutitang Formation shale. The Niutitang Formation, which is located in Songtao Wuluo well ZK408 in the eastern part of the survey area, is present at a depth of 616.90~640.80 m. In the Deye 1 well in the northern part of the study area, the floor depth of the Niutitang Formation is at 1661.66~1752.20 m. In the drilling record of Well Fengshen 1, the Niutitang Formation is approximately 2500 m underground in Liujiazhai Village, Dangwan Township, Feng'gang County. The temperature of the stratum increases as depth increases. According to valid statistics, every time the depth increases by 100 meters, the temperature increases by 3°. The temperature at the Niutitang Formation in the study area is approximately 45°C~120°C. The amount of methane adsorption is greater at low temperature and low pressure, and the amount of free methane at high temperature and high pressure is greater.
To approximate these real conditions, this study is conducted from 30°C to 90°C.

2.
3. Pore Characteristics of Shale. Nanoscale pores are prevalent in shales, but the connectivity between pores is negligible. Shale gas is mainly stored in the micropores of shale. The size and connectivity of pores have a large impact on the ability of shale to adsorb shale gas and its storage performance. The degree of pore development in shale determines the permeability of shale and the degree of gas reservoir storage. Using nuclear magnetic resonance (NMR) technology, the pores of 5 groups of samples in the rock reservoirs in the study area are quantitatively tested, and the pore throat and pore size distributions are obtained [37]. Figure 3 shows that the pore throat and pore size distribution of each group of samples are approximately similar. Three peaks appear in the pore size distribution, indicating that the pore size range of shale is primarily concentrated at 0.001~0.1 μm and 0.01~0.4 μm, and a very small amount is distributed at 1~10 μm.
Through NMR analysis of the Lower Cambrian shale in the Feng'gang No. 3 block, it can be seen that the development of pores in the Lower Cambrian shale is advantageous to the storage of shale gas. This indicates that the development of shale gas in this area is economically important. It can also be seen from Figure 3 that although the micronscale pores are well developed, the connectivity between the pores is not good, and increasing the temperature helps to increase the porosity and permeability of the shale. It is therefore important to study the mechanism of damage under thermosolid coupling for the extraction of shale gas.

Thermosolid Coupling Element Model
3.1. Element Thermal Damage Evolutionary Theory. At the microscopic scale of shale, an element is considered to be an elastic body at the initial failure stage, and the stressstrain curve of each element is linear. When the damage threshold is reached, the properties of the unit begin to change. When shale is compressed, stretched, or sheared, 3 Geofluids the tensile strength of shale is much smaller than the compressive strength, and tensile failure occurs, as shown in Equation (1). When shale reaches a critical state of stress, the shear stress exceeds the shear strength, and shear failure occurs, as shown in Equation (2) [38,39].
F 1 and F 2 are the two damage threshold functions, and σ 1 and σ 3 are the maximum principal stress and minimum principal stress, respectively. In terms of the sign, the compressive stress is positive, and the tensile stress is negative.
In actual situations, f t0 in the formula is expressed in terms of uniaxial tensile strength. f c0 is the uniaxial compressive strength of shale, and φ is the internal friction angle of the material. In the case of uniaxial compression and uniaxial tension, f c0/ f c0 = f c0 /f c0 = λ indicates that the formula is applicable.

Element Thermosolid Coupling Equation.
Assuming that the shale is an elastic body, the constitutive relationship satisfies the generalized Hooke's law. The shale in the study area is at high temperature and high pressure. Under this condition, Wee et al. proposed that the relationship between temperature variables and shale under combined loading is expressed as follows [40]: σ ij ′ = σ ij + ξpδ ij , δ represents the total stress tensor (symbol is positive for tension), δ ij is the Kronecker symbol, G is the shear modulus of the material, ν is Poisson's ratio, ξð≤1Þ is the compressibility of the material ½ξ = 1 − ðK ′/K S Þ, K S is the shale material modulus of elasticity, and α is the volume expansion coefficient (K −1 ). In this formula, the strain and stress caused by temperature changes are relative to the initial temperature; therefore, T is the temperature increase compared with the initial temperature.

Geofluids
During the research process, the interaction between mechanical energy and thermal energy is ignored. Based on the REV heat balance equation, it can be expressed as follows [26]: In this formula, T 0 is the initial temperature (K) in the unstressed state, ρ s is the density of the shale (kg/m 3 ), C s is the specific heat capacity (J/kg·K), and λ s is the heat transfer coefficient H of the shale (J/s·m·K). Equations (3) and (4) are a set of fully coupled nonlinear equations that control the thermosolid coupling of shale. This equation explains the evolution of heat transfer, compressibility, and thermal damage between shale materials under thermosolid coupling. The nonlinear differential equation is controlled by the second-order space partial derivative and the firstorder time partial derivative. The nonlinear differential equation in the space and time domains is unsolvable.
3.3. Numerical Simulation. The RFPA temperature version software is used to study crack propagation under the thermosolid coupling of shale. RFPA is software that is based on the finite element, statistical damage theory, and Coulomb failure criterion. The deformation and failure of the rock are analyzed to obtain the true failure mode. This calculation method considers fully the randomness of the material distribution, such as rock heterogeneity and pores, and combines it with the statistical distribution assumption of material properties into the finite element algorithm. The division unit meets the given strength criterion; therefore, the rock under real conditions and the destruction process can be presented through numerical simulation. In actual situations, many difficult-to-implement problems can be solved. Figure 4 shows the schematic diagram of controlling the displacement of the model with bedding inclinations of 0°~90°.
The discretization unit of RFPA obeys the Weibull distribution. On this premise, the relationship between the microscopic conditions and the mechanical properties of the macroscopic medium is established, and the Weibull statistical distribution function is introduced to describe the following [41]: where α is the mechanical property parameters of the primitive body, α 0 is the average value of the mechanical property parameters of the material, m is the parameter reflects the homogeneity of the material medium in a physical sense and is defined as the uniformity coefficient, and φðαÞ is the statistical distribution density of the primitive property α.
In the RFPA temperature version, the temperature distribution function T ðx, y, z, tÞ is used to express the functional expression of the temperature change with time at each certain point in the rock.
In the calculation, the heat q ðx, y, z, tÞ, q = dQ/dt/S, and Q per isothermal area are considered heat, t is time, and S is area, where q is proportional to the gradient of temperature function u, and the proportional coefficient H is the thermal conductivity.    Table 1) below. The experimental study is a numerical model of shale with different bedding angles under the condition of gradually increasing temperature. The bedding angles in the model are set to 0°, 22.5°, 45°, 67.5°, and 90°. A schematic diagram of the model's overall conditions is shown in Figure 5. The fixed temperatures for the shale numerical model are 30°C, 60°C, and 90°C. The temperature boundary condition is set as the first type boundary condition (temperature is constant at 30°C, 60°C, and 90°C), and there is no heat flow between the shale matrix and outside world. Figure 4 shows the stress-strain curves for shale at 30°C and dip angles of 0°, 22.5°, 45°, 67.5°, and 90°in a and 60°C and 90°C in (b) and (c). As seen from the graphs, the corresponding stress-strain curves for the shale change significantly with different inclinations, with temperature having less influence on the stress-strain curve than the effect of shale inclination. The damage of shale is broadly divided into four stages: the linear elastic stage, yielding stage, damage stage, and stability stage.

Analysis of Mechanical Characteristics of Shale Rupture.
(1) The first stage is the linear elastic stage. In this stage, stress and strain have a linear relationship and exhibit good elastic properties. As seen in Figure 4, the linear elastic phase of the shale decreases significantly with increasing temperature and laminar dip angles of 22.5°and 45°. This indicates that the increase in thermal stress within the shale reduces the cementation between the shale matrix and laminae, and the softening effect becomes more pronounced as the lamina dip angle increases (2) In the second stage of the yielding stage, influenced by the surrounding pressure and thermal stress, the temperature increase causes plastic deformation within the shale, and the shale stress-strain curve under the effect of thermosolid coupling shows a more obvious yielding stage (3) The third stage is the failure stage. When the strength reaches its peak, the stress decreases sharply, and shale damage occurs. The graph shows that the failure strength of the shale decreases to a certain extent at each laminar dip angle after the temperature rises, especially at laminar dips of 22.5°and 45°. When the dip of the shale is the same and is accompanied by a decrease in temperature, the reduction in peak strength is due to the difference in the coefficient of thermal expansion between the bedding and shale matrix within the shale. When the particles are heated, they cause uneven expansion and uneven deformation of the shale. The particles squeeze each other so that a temperature-induced thermal stress is developed in the shale. As the temperature increases, the thermal stresses increase, resulting in thermal damage to the shale. More microfractures or primary fractures appear in the shale. On a macroscopic scale, the temperature of the shale increases. In this case, the strength decreases, which is consistent with the results of Yan et al., who studied mudstones [42] (4) The fourth stage is the stabilization stage, where the tail end of the curve is smoothed. It can be seen from the graph that the stabilization phase at 30°C is shorter than that at 60°C and 90°C. At a temperature of 30°C, the damage mode of the shale mainly shows a brittle damage state. As the temperature increases, the properties of the shale change from brittle to plastic, showing an increase in the stability phase  Table 2 shows the compressive strength of shale with different bedding dip angles obtained by numerical simulation at 30°C, 60°C, and 90°C. The compressive strengths from Table 2 are plotted in Figure 6, from which it can be seen that the compressive strength of the shale reaches its maximum at 30°C, while the compressive strength decreases with increasing temperature. Additionally, temperature has a strong effect on the anisotropy of shale bedding, with increasing temperature leading to reduced cementation of the laminae to the shale matrix, and the effect is strongest at lamina dips of 22.5°and 45°. At temperatures of 30°C and 60°C, the compressive strength of the bedding inclination is greatest at 0°and least at 45°. As the azimuth of the bedding rises, the compressive strength of the shale shows a trend of first falling and then rising. When the temperature is 90°C, the compressive strength is greatest at a bedding azimuth of 90°. Alternatively, the compressive strength is lowest at 45°because of the weak cementation between the shale matrix and weak bedding surface. As the temperature increases, the thermal stress resulting from the incompatible expansion of the laminated particles and shale matrix gradually increases, and the compressive strength of the shale specimens decreases under the combined effect of the thermal stress and lamination effect.
As seen from the variation curve in Figure 6, the encrusted shale exhibits a strong lamina effect. To represent the effect of laminations on the mechanical properties of shale under thermosolid coupling, SðαÞ is used to characterize the structural effect of the laminations. The bedding effect is 1 minus the compressive strength at different inclination angles divided by the compressive strength measured at the inclination level. The expression is as follows [24]: In the formula, SðαÞ is the bedding effect coefficient, θðαÞ is the compressive strength of different bedding dips, and θ ð0°Þ is the compressive strength of the shale when the bedding dip is 0°. According to the above formula, Table 3 shows the bedding effect coefficients of different shale dip angles. To observe the changes in bedding effect coefficients more intuitively, bedding effect diagrams of different bedding dip angles at 30°C, 60°C, and 90°C are drawn.
As shown in Figure 7, the shale lamina effect rises more slowly at 60°C with increasing temperature, with a significant effect at 90°C. The laminar effect is not evident for dips of 0°versus 90°, and temperature changes within 30°C to    The temperature rises to 90°C, the thermal stress from the uneven expansion of the particles within the shale rises sharply, and the shale matrix turns to plastic deformation between the shale and laminae, leading to a weakening of cementation; therefore, the lamina effect increases, and the fold line in the diagram shows a rapid rise.

Failure and Deformation Characteristics of Shale under
Thermosolid Coupling. During the evolution of shale damage by thermosolid coupling, shale specimens reach peak strength followed by a sharp decrease in stress and sudden rupture. Figure 8 shows the damage process and the corresponding acoustic emission variation pattern for shales with different bedding dips at a fixed confining pressure. As seen from the diagram, the shale undergoes three stages: fracture initiation, fracture extension, and penetration. Figure 8(a) shows the fracture initiation phase of the shale damage process, Figure 8   (1) N-type (90°3 0°C , 67.5°6 0°C , 90°6 0°C , and 45°9 0°C ). During the evolution of fractures in these groups of shale specimens, cracks first appear in the laminae. As loading progresses, the cracks extend through the laminae and eventually penetrate the shale matrix along a 45°or 135°dip to form an N-fracture. When the laminar inclination is at a high angle, temperatures below 60°C have less effect on the damage pattern, and temperatures of 90°C strengthen the laminar effect. The lamina effect is stronger at high angular lamina dips, and the laminae have a strong dominant effect on the overall shale integrity. As seen from the acoustic emission diagrams in Figure 8(d), these shale specimens are mainly distributed in tensile stresses during the final damage phase, but there are a number of shear stresses in the laminae, which can be judged to be a gradual increase in tensile stresses in the shale matrix after shear damage has occurred in the laminae.
(2) ʌ-type (45°3 0°C and 0°9 0°C  (4) Slanted I-type (0°3 0°C , 67.5°3 0°C , 0°6 0°C , and 22.5°6 0°C ). Damage to shale specimens at low laminar angles is primarily subject to tensile and shear stresses, when the surrounding pressure plays a major factor in preventing multiple cracks from forming. The specimen with a high angle of 67.5°at a temperature of 30°C is then subject to the effects of perimeter pressure and laminae, with cracks expanding along the laminae and the perimeter pressure preventing the formation of multiple cracks, resulting in an oblique I-shaped damage pattern.
(5) Cluttered-type (22.5°9 0°C , 67.5°9 0°C , and 90°9 0°C ). At the higher temperatures of the 90°C complex, irregular disorganized damage is more likely to occur. At a temperature of 90°C, the laminar and temperature effects of the shale reach their maximum at the same time, and thermal and loading stresses within the shale act together to cause the shale to exhibit a complex damage pattern.
According to the above analysis, the surrounding pressure at low angles of the shale laminae is the dominant factor in the damage pattern with simpler forms of damage and more complex damage occurring as the temperature rises. The lamina effect is stronger at high angles in shale laminae, where the lamina effect is the main controlling factor, mainly manifesting itself as damage occurring first along the laminae, with the lamina effect acting in conjunction with the temperature effect to produce a more complex damage pattern as the temperature rises. During hydraulic fracturing, hydraulic fractures may extend along high angular laminar dips, inhibiting the expansion of the fracture network, while higher temperatures contribute to the formation of complex fracture networks, thereby increasing shale gas production.

Quantitative Fractal Characterization of Shale
Damage under Thermosolid Coupling. Analytical theory can be used to quantitatively describe irregular and complex matters in nature, including the evolution of damage to rocks. Shale converts its internal stored energy into elastic waves and releases them rapidly at the moment of destruction, a phenomenon known as acoustic emission [43,44]. During shale damage, acoustic emission occurs for each unit that breaks down; therefore, acoustic emission has a fractal character. In this paper, images of acoustic emissions at different stress levels are grayed out by ImageJ and imported into a MATLAB calculation program to find their fractal dimension. The fractal dimension is solved by the following [45]: where D s is the adaptive fractal dimension of the damaged 12 Geofluids region, r is the side length of the square box, and NðrÞ is the number of square boxes with side length r needed to cover the damaged region therein. Table 4 shows the fractal dimension values. Figures 9  and 10 show the trends between AE energy, fractal dimension, laminar dip, temperature, and stress level. As shown in Figure 9, the AE energy generally shows an increase with increasing stress levels. At stress levels less than 50%, the AE energy of each group is almost zero, showing little or no damage. When the stress level is between 50% and 80%, the AE energy rises slowly at this stage, with a more pronounced rise in AE energy in low-angle shale specimens compared with high-angle shale. When stress levels exceed 80%, AE energy rises rapidly, and shale damage breaks down. 22.5°, 45°, and 67.5°show significant increases in AE energy, corresponding to the macroscopic damage forms, with high release energy representing intense shale damage and injury. An increase in energy is observed as the temperature rises, at which point the macroscopic damage intensifies, with an insignificant increase in energy at 60°C and a significant increase at 90°C, when the corresponding damage pattern becomes more complex. Figure 10 shows the fractal dimension as a function of laminar dip, temperature, and stress level, with the x-axis coordinates representing stress level and the y-axis coordinates representing temperature. It is clear from the graphs that the fractal dimension generally increases with increasing temperature, and at α = 22:5°, 67.5°, and 90°, the fractal dimension increases more with increasing temperature, indicating that temperature enhances the laminar effect of the three groups of shale dip specimens, and the macroscopic phenomenon shifts from simple oblique I or N damage to more complex heterogeneous damage. Conversely, smaller fractal dimensions correspond to simpler forms of destruction.

Conclusion
In this paper, the effect of temperature on the damage process of shales containing laminated shales with different dip angles is studied by establishing a coupled thermosolid model for shales of the Niutitang Formation, and the following laws are summarized: (1) The temperature increase reduced the linear elastic phase of the shale specimens in each group, with a significant reduction in the linear elastic phase of shales with lamina dip angles of 22.5°and 45°. The elevated temperature causes partial plastic deformation within the shale, and the shale stress-strain curve under thermosolid coupling exhibits a more pronounced yielding phase (2) When the temperature rises from 30°C to 60°C, it makes the unit swell leading to an increase in pore space, which further compacts the shale matrix with the laminae under displacement-controlled loading, and the lamina effect decreases to a small extent. And after the temperature reaches 90°C, the thermal stress generated by the uneven expansion of particles inside the shale rises sharply. The shale matrix and laminae turn to plastic deformation between them leading to weakening of cementation, and the strengthening effect of laminae in shale is most significant (3) The shale damage patterns are grouped into five categories (N-Type, ʌ-Type, v-type, slanted I-type, and cluttered-type). Among them, the occurrence of Ntype damage is mainly due to the strong dominant effect of high-angle laminae on the integrity of the whole shale, and the warming reinforces the lamina effect. The ʌ and v shapes mainly occur when the dip angle of the laminae is 45°. The sloping Ishaped damage occurs in shale specimens with low laminar dip, which are mainly subjected to tensile and shear stresses, when the surrounding pressure plays a major role in preventing the formation of multiple cracks. Heterogeneous damage occurs at a temperature of 90°C. The lamina effect and temperature effect of the shale reach a maximum at the same time, and the thermal and loading stresses inside the shale act together to give the shale a complex damage pattern (4) The active role of temperature in shale damage is further quantified by the fractal dimension. The fractal dimension curve is relatively flat when the temperature is 60°C, while the fractal dimension rises rapidly by 90°C, indicating that the fractal extension of the shale is most favorable at a temperature of 90°C. At α = 22:5, 67.5°, and 90°, the fractal dimension increases with increasing temperature, indicating that the laminar effect of the three groups of shale dip specimens is strongly influenced by temperature, and the macroscopic phenomenon shifts from simple oblique I-or N-shaped damage to more complex heterogeneous damage

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.