Next Article in Journal
Recording of Full-Color Snapshot Digital Holographic Portraits Using Neural Network Image Interpolation
Next Article in Special Issue
Analysis of the Effect of Ultra-Fine Cement on the Microscopic Pore Structure of Cement Soil in a Peat Soil Environment
Previous Article in Journal
Intelligent Jammer on Mobile Network LTE Technology: A Study Case in Bucharest
Previous Article in Special Issue
Analysis of the Influence of Initial Stress on the Bandgap Characteristics of Configuration-Controllable Metamaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Mechanical Properties of Heat-Treated Shale under Different Temperatures

1
Institute of Defense Engineering, Academy of Military Sciences (AMS), Peoples Liberation Army (PLA), Beijing 100850, China
2
School of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
3
Shock and Vibration of Engineering Materials and Structures Key Lab of Sichuan Province, Southwest University of Science and Technology, Mianyang 621010, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(22), 12288; https://doi.org/10.3390/app132212288
Submission received: 29 August 2023 / Revised: 7 November 2023 / Accepted: 8 November 2023 / Published: 13 November 2023

Abstract

:
As a typical rock, shale’s reservoir depth is about 1500–4000 m, and the temperature of the shale reservoir at this depth is 150 °C. Therefore, in order to study the dynamic strength of shale at this temperature, it is necessary to consider the effects of temperature and strain rate on the dynamic strength of shale, and then establish the damage constitutive model of shale. This paper took black shale from the Sichuan Basin as the research object, combined it with the separated Hopkinson bar experiment and temperature control system, and conducted the Hopkinson bar experiment on shale at room temperature, 60 °C, 90 °C, 120 °C, and 150 °C, and at three groups of air pressures of 0.2 MPa, 0.3 MPa, and 0.4 MPa. The stress–strain curves of shale at the same strain rate and different temperature and at the same temperature and different strain rate were obtained. In the temperature difference range of this experiment, the dynamic strength of the sample presented two opposite trends (increasing and decreasing) with the increase in temperature, which was determined via the direction of the bedding. The peak strength linearly increased with the increase in strain rate. Based on the Weibull statistical distribution and the D–P failure criterion, a statistical damage constitutive model of shale dynamic strength considering the effects of temperature and strain rate was obtained. By modifying the parameters F0 and m, the dynamic strength statistical damage constitutive model of shale was in good agreement with the experimental results.

1. Introduction

Deflagration fracturing technology [1], also known as high-energy gas fracturing (HEGF) technology, is a unique oil and gas field stimulation technology [2]. The advantages of HEGF technology applied to rocks are as follows: First, stress waves can lead to the formation of multiple radial fractures that communicate with natural fractures. High temperatures and high-pressure gases can effectively extend the fractures and increase the content of shale gas flowing into the wellbore to a large extent. Second, deflagration fracturing is suitable for fracturing stimulations of liquid-sensitive oil and gas reservoirs. Third, the fracture itself is cut, resulting in the dislocation between layers to form support, so the fracture is not easy to close; fourth, in the actual project, deflagrant compression cracking does not require too much large equipment, which can meet the requirements of economic practicality; fifth, as a waterless fracturing technology, detonation fracturing has achieved a better environmental protection effect.
HEGF generally forms a pre-fracture through perforation, and then injects a liquid drug into the perforation section to generate a large amount of high-temperature and high-pressure gas after ignition, which exerts an effect on the pre-fracture and further expands and extends it, as shown in Figure 1a. According to the physical model, the geometric model of fracture propagation in high-energy gas fracturing of horizontal wells is defined, as shown in Figure 1b. σH is the maximum principal stress, σh is the minimum principal stress, L is the reservoir length, W is the reservoir width, and Pt is the gas pressure distribution in the fracture [3].
In terms of shale mechanics research, Horsrud et al. [4,5,6] introduced a series of shale formations and analyzed their series mechanical parameters via heating. It was found that the shale has strong anisotropy and that the temperature has an obvious influence on the anisotropy, porosity, and permeability of the shale. The heating of the shale significantly reduces the strength, stirring speed, and sound velocity of the shale, and tests from other relevant scientific procedures have concluded that the shale may cause changes in the interlayer structure due to pressure and temperature unloading. Based on this, a physics-based micromechanical model was created in order to quantify the damaged state of shale under loading through macroscopic measurements. Niandou et al. [7,8] found that the volume strain of shale mainly comes from compression deformation through the experiments of compressibility and triaxial compression tests of shale under different confining pressures and loads. The transition to dilatancy only occurs in the region close to the peak stress. Increasing the magnitude of the isotropic stress field has little effect on the velocity anisotropy. Luo et al. [9] studied shale with different bedding (0°, 30°, 45°, 60°, and 90°) and conducted a study through the triaxial SHPB test with a confining pressure of 0~15 MPa. It was found that the dynamic mechanical properties of shale mainly depend on the interaction between bedding dip angle and various confining pressures. Huang et al. [10] used the split Hopkinson pressure bar test to study the dynamic compressive mechanical response of shale specimens with five different bed dip angles after heat treatment at 25–200 °C. It was found that with the increase in heat treatment temperature, the failure mode of samples with bedding dip angles of 45° and 60° were significantly affected, but the bedding dip angles of 0°, 30, and 90° were not fundamentally changed.
Since hydraulic fracturing technology is mostly used in the exploitation of shale gas abroad, research on the mechanical properties of shale is mainly focused on static research, and the exploration on the dynamic mechanical behavior of shale is relatively rare. Although our country started late in shale research, some substantial progress has been made. Dong et al. [11,12,13,14] conducted uniaxial and triaxial compression tests on shale, analyzed the mechanical properties, strength, and failure mode of shale, and revealed the anisotropic failure mechanism of shale. Kang et al. [15] divided the pyrolysis weight loss of oil shale into three stages based on the new technology of in situ thermal injection mining of oil shale. When the temperature exceeds 200 °C, the influence of temperature on bulk density, specific gravity, porosity, etc., increases. When the temperature exceeds 300 °C, there is a significant increase in the number, length, and width of shale fractures. At the same time, microcracks are formed in the direction perpendicular to the bedding plane. Deng [16] et al. studied the anisotropy of shale generated by layered waves and gave the P-wave and S-wave velocities of samples in different directions under different conditions. Wang Yi [17] found that fracturing fluid intrudes into rock mass along with fractures, which is the main reason for the decline in shale strength. Liu et al. [18,19] carried out uniaxial compression mechanical experiments of shale under four kinds of low strain rates and analyzed the influence of shale strain rates on the relevant mechanical parameters of shale. Zhang [20] conducted the separation Hopkinson bar compression experiment for the high strain rate of shale, and preliminarily explored the dynamic constitutive relationship of shale. Yang et al. [21] used the separated Hopkinson pressure bar system to carry out impact tests on shale. It was concluded that the damage degree of shale significantly increases with the increase in impact velocity under different strain rates, and that the ductility and failure resistance of rock materials are improved under passive confining pressure.
Crack initiation and its effects on the dynamic behavior of rock, particularly regarding ultrasonic wave velocity, have been studied extensively in recent years. Numerical manifold method (NMM) simulations have been used to investigate the propagation of stress waves through the fractured rock mass [22]. A study by Fan [23] demonstrated that the stress wave velocity is affected by the size and orientation of the fractures, as well as the properties of the surrounding rock’s matrix. Micro-fractures have also been found to impact the dynamic behavior of rocks significantly. A study by Zhou [24] on wave propagation through rock mass with micro-fractures found that the velocity and amplitude of the waves are affected by the orientation and distribution of the micro-fractures. Thermal shock-induced damage is another factor that can affect the dynamic behavior of rocks. The spatial gradient distribution of damage to granite induced via thermal shock has been studied. It was found that the thermal shock-induced damage significantly affects the ultrasonic wave velocity and attenuation [25]. Finally, experimental investigations of thermal effects on the dynamic behavior of granite have demonstrated that changes in temperature can cause significant changes in the ultrasonic wave velocity, attenuation, and dispersion [26]. These studies have contributed to a better understanding of the effects of crack initiation and thermal shock-induced damage on the dynamic behavior of rocks. Different material grains have different anisotropy and thermophysical properties, which leads to incongruous deformation across grain borders and thermal stress between or within grains [27]. When the thermal stress is greater than the material’s strength, cracking takes place, causing a lot of cracks and a corresponding drop in strength. However, the mechanical properties of various rocks vary with temperature as a result of variations in mineral composition, particle size, and microstructure. Research has shown that Yangbajing geothermal wells in China and larderello geothermal wells in Italy both have temperatures above 300 °C [28]. As a result, the issue of rock engineering under high-temperature environments has emerged as a key area for research in the field of rock mechanics [29].
The physical and mechanical properties of shale exhibit significant anisotropy, impacting the design and construction of high geothermal tunnels and deep resource mining. The research on the stress–strain relationship of shale under different conditions and different strain rates is not perfect domestically and abroad. Therefore, this paper took shale from the Sichuan Basin as the research object and conducted Hopkinson bar compression experiments at normal temperature, 60 °C, 90 °C, 120 °C, and 150 °C, and at three groups of air pressures of 0.2 MPa, 0.3 MPa, and 0.4 MPa. The differences and similarities of the stress–strain effects of shale under the same strain rate with different temperatures, and the same temperature with different strain rates, were analyzed in depth. Based on the experimental results, a constitutive model considering shale temperature was established. The research results can provide an important theoretical basis for the strength, fracture, and deformation theory of shale. The research process is shown in Figure 2.

2. Preparation of Shale Specimen

2.1. Sampling and Processing

The samples were taken from Liutang Township, Shizhu County, Chongqing, on the edge of the Sichuan Basin, which can better represent the rock profiles of black shale. Engineering bricks with an inner diameter of 56 mm were used for core sampling and a core length of 150 mm. The cores were quickly sealed with cling film and kraft paper and sampled and numbered after being taken out for use and subsequent processing later. The sampled shale is shown in Figure 3. The rough sampling method described above leads to various defects in the cores and therefore requires secondary machining to meet the strict requirements of the SHPB test that the specimen must be in close contact with the incident bar and transmission bar. The finishing process ensured that the final experimental specimens have a more uniform length, a flatter section, good parallelism, and no bending of the specimens. The processing apparatus used was a C6132 machine Ltd. (Tengzhou, China) [30] of Yunnan Yunji Group Import & Export Co., Ltd. (Yuxi City, China) to finely process the sample. The cores were machined into standard rock samples with a height of 50 mm and a diameter of 50 mm, and the tolerance of dimensions, parallelism, and diameter were ± 0.02 mm.

2.2. Shale CT Scanning and Material Composition

CT scanning technology has been widely used in industry, including X-ray tomography, Compton scattering tomography, and Mossbauer effect tomography, which are mainly used to detect the interior of industrial products. In this paper, X-ray scanning technology was employed in the CT scanning experiments, which has the characteristics of high spatial resolution, non-destructive detection, thin scanning layer, good image quality, and fast speed, which can better reconstruct the cross-section of the shale sample and detect the various defects, such as cracks, delamination, and impurities, in the samples.
Five groups of CT scans of the shale were analyzed. A comparison of the original shale samples and CT results is shown in Figure 4, showing the original and scanned images of samples 1# to 5#. The 1# original sample has obvious impurities and small radiating cracks, which could be seen clearly through the CT scan. No cracks were observed in the original sample of 2#, but the CT scanning results showed that the cracks were located to the right of the center of the shale circle and did not penetrate the sample. The original samples of 3# and 4# were observed via visual and CT scans to have distinct joints’ planes, and the joints’ planes were roughly parallel. The 5# sample was relatively dense overall, with little tramp mass. Based on the CT scan results of the shale, there was a high probability of defects, such as joints and fractures, in the shale itself, but the shale was relatively dense throughout. It can be inferred that the existence of defects in the shale itself has a certain impact on its mechanical properties.
The structure of natural shale is complex, and the composition and content of the material directly affect the mechanical properties of the shale. The shales were quantitatively and qualitatively analyzed via XRD diffraction tests. The composition of each substance is outlined in Table 1 [31].

3. Splitting and Stretching Test

The static splitting test of three kinds of loading angle specimens was carried out using an electronic universal testing machine and a gray-dot camera. The main crack propagation, secondary crack generation, and final rupture of the specimens were different during the bad form loading process.
Figure 5 shows the final failure diagram of the specimen with a loading angle of 0°. The test starts from the initial crack tip during the loading process, and gradually expands to form a “one” shape main crack with the increase in the load. The principal crack expansion morphology of the three specimens was similar. With the increase in load in the initial crack during loading, stress concentration occurs in the crack tip region. When the initial crack tip reaches the condition of cracking, the initial crack tip begins to crack and extends to the upper and lower ends of the loading along the diameter loading direction, forming a through principal crack that divides the specimen into two parts. Due to the Poisson effect in the loading process of the specimen, when the effective tensile stress exceeds the tensile strength of the specimen, the splitting tensile failure occurs along the diameter of the loading plane. Figure 5 shows the final failure diagram of the specimen with a loading angle of 0°. Table 2 shows the splitting tensile strength of shale.

4. Shale SHPB Experiment

At present, a split Hopkinson pressure bar is commonly used to test the dynamic characteristics of brittle materials, such as concrete and ceramics [32,33,34]. The impact experiment of the shale was carried out using a SHPB device with a diameter of 50 mm. The test was carried out on the split Hopkinson experimental system in the structural dynamics laboratory of Southwest University of Science and Technology. It is mainly composed of a bullet, an incident bar, a specimen, a transmission bar, an absorption bar, a measuring device, and a data processing system. The schematic diagram of the experimental device is shown in Figure 6, where the lengths of the bullet, incident bar, and transmission bar are 600 mm, 3000 mm, and 3000 mm, respectively. Each bar has a diameter of 50 mm and is made of a steel alloy with a modulus of elasticity (E) of 72 gPa, a yield strength of 300 MPa, and an ultimate strength of 460 MPa. The shale specimen was placed between the incident bar and the transfer bar, where the contact surfaces are parallel. Strain gauges were installed on the incident bar and the transmission bar to collect waveform data. The actual experiment is shown in Figure 7. The position of the strain gauge was two times the length of the bullet from the incident end face. The strain gauge of the incident bar and the strain gauge of the transmission bar were symmetrical with respect to the specimen.
The pulse-shaping technique was used in SHPB testing to promote dynamic stress equilibrium in the sample to reduce the effects of the heterogeneity of the shale material and the effects of high-frequency impacts, which can improve the stress homogeneity of the sample. In this paper, by comparing the effects of different materials as forming sheets, an appropriate material was selected so that the rise time of the incident wave was 300 µs or more, i.e., the stress of the specimen was much less than the peak stress when the stress was uniform. It can be assumed that the experimental data obtained achieve a constant strain rate loading process of shale materials and satisfy the assumption of stress uniformity.
In the process of the waveform shaping experiment, according to the experience of predecessors [35], butter, copper, brass, rubber, and other materials were added to the front end of the incident bar to shape the curves. The circular sheets with a diameter of 15 mm and a thickness of 2 mm were prepared, and the SHPB empty bar experiment was completed, and the incident waveform diagrams of the different materials are shown in Figure 8.
According to the experiment, the loading time reached 400 ms when rubber was added as the shaping piece, and the waveform was similar to a sine wave. After adding the rubber shaping sheet, high-frequency impacts were effectively avoided. After analyzing the incident pulses of shaping pieces of different materials, rubber was selected as the shaping piece in this experiment, and its average loading time reached more than 300 ms. By calculating the stress balance time required by the shale, it can be seen that the stress balance time of the test pieces is far less than the time required to reach the peak stress. Therefore, the selected rubber shaping sheet can sufficiently meet the requirements of stress balance in the later analysis.
At present, there are few SHPB compression experiments on shale with temperature. High-temperature dynamic mechanical experiments have been proposed much earlier. As research on the dynamic mechanical properties of brittle materials has increased, SHPB experiments on concrete and rock at high temperatures have gradually increased, but SHPB experiments on shale with temperature are very rare. Previous research has indicated that the shale reservoir in the Sichuan Basin is approximately 1500–4000 m, and that the temperature of the shale reservoir reaches 150 °C at this depth. Therefore, several groups of tests at different temperatures have been designed to investigate the dynamic mechanical properties of the shale. The shale heating device is equipped with the SHPB experimental procedure described above. The heating furnace and temperature controller are shown in Figure 9.
In the SHPB experiment on the shale with temperature, the temperature of the equipment was first adjusted to the specified temperature. When the temperature reached the specified temperature, the high-temperature device was held at the specified temperature for one hour before testing to ensure that the internal temperature of the shale also reached the specified temperature, which was closer to the actual environment of the shale.

5. Analysis of the Experimental Results

5.1. Stress Balance Check and Typical Failure Modes of Shale Specimens

In shale SHPB compression bar tests, it is an effective measure to ensure the validity of the experimental results by verifying that a state of stress equilibrium is reached before the shale is damaged. Formula (1) can be used to verify whether the stress of the SHPB test meets the balance requirement:
W i + W r   =   W t
where W i ,     W r , and W t represent the incident stress wave, reflected stress wave, and transmitted stress wave, respectively. When the data were processed, the corresponding incidence curve, reflection curve, and transmission curve were used to replace the calibration, and the results are shown in Figure 10.
As can be seen from Figure 10, the stress balance at both ends of the test specimen was well achieved at different impact pressures corresponding to the three shale failure modes, i.e., this test can meet the stress balance requirements at different strain rates. However, this experiment requires a high level of accuracy, and the stress equilibrium state of the specimen needs to be verified each time during the subsequent data processing stage. The results that reach the equilibrium state can be used, while the failed data can be abandoned, to ensure the effectiveness of the dynamic mechanical analysis of shale materials.
The failure mode of shale specimens at different strain rates were obtained from the tests by adjusting the experimental air pressure. When the pressure was changed, the impact velocity of the bullet changed, resulting in a change in the test strain rate and, thus, a significant difference in the final damage pattern of the shale specimens. In the process of loading, a large number of small cracks will appear in shale due to the increasing load, and small cracks will form large cracks. Shale is a kind of highly brittle rock mass, and three different failure modes are produced with different impact velocities. As shown in Figure 11, the damage degree of the shale under different impact velocities and different strain rates can be divided into a simple fraction. Based on the reference to the concrete damage modes, the shale can also be divided into three damage forms: “cracking and spalling”, “fracturing and crushing” [36]. Through a large number of experiments, waveforms that are representative of the shale damage patterns were selected for analysis to verify the validity of the experiments. The pressures were 0.2 MPa, 0.3 MPa, and 0.4 MPa, and the velocities were 4.025 m/s, 6.126 m/s, and 8.219 m/s, respectively. The corresponding average strain rates were 44.62/s, 73.6/s, and 104.07/s, respectively.
When the impact pressure of the test was 0.2 MPa, the bullet velocity was 4.025 m/s, and the average strain rate was 44.62/s. As shown in Figure 10a, the specimen is slightly spalled and cracked locally under the impact pressure, but the specimen is not damaged. This phenomenon is called “cracking and spalling”. When the impact pressure was 0.3 MPa, as shown in Figure 10b, the average strain rate was 73.6/s, and the whole specimen was divided into several large pieces, which can be considered damaged; we called this “rupture”. When the impact pressure was 0.4 MPa, the average strain rate reached 104.07/s, as shown in Figure 10c; the specimen was completely crushed under the impact pressure, and the main body of the specimen was completely divided into powder and blocks. This kind of phenomenon is called “crushing”. Therefore, it can be concluded that the damage degree of the shale intensifies with the increase in impact pressure.

5.2. Stress–Strain Curve of Shale Material at Room Temperature

The validity of the experimental data was verified through a large number of SHPB experimental tests, and the stress–strain curves of the shale at room temperature were obtained, as shown in Figure 12.
Due to the difference in strain rates, and the influence of internal voids, fractures, joints, and other defects, as well as their discontinuity, anisotropy, and heterogeneity, the stress–strain curves of shale are quite different, but the types of stress–strain curves are generally the same under the strain rates, with little difference in their values. The effects of internal defects, discontinuities, anisotropy, and heterogeneity of the shale on the shale SHPB bar compression experiments decrease with increasing strain rates, and the stress–strain curves of the shale are similar. The shale stress–strain curves begin to show strong characteristics of linear elastic. When the strain rate is low, the corresponding damage mode indicates that the shale is only locally spalled, which does not affect the overall damage. Therefore, the unloading section of the shale can still bear a degree of compressive pressure. However, the high strain rate shows two modes of damage, which are rupture and fragmentation. In this case, the shale has shown overall damage, and the compressive strength of the shale drops sharply to zero. Consistent with the findings of Abbas et al. [22], the destruction process can be divided into four stages. The first stage: when the stress gradually increases with the strain, the load seat and the natural micro-crack close. The second stage: with further compression, the shale gradually densifies, compacts, or consolidates. The third stage: the deformation of shale can be widely recovered in the elastic compression and micro-crack stage. The fourth stage: the accelerated micro-crack or macro-crack stage, whereby the stress–strain curve basically deviates from the linear elastic deformation until it reaches its peak stress and failure. The peak strength of shale increases significantly with an increase in strain rate, indicating that shale is a brittle material with a strain rate effect.

5.3. Dynamic Compressive Strength Analysis of Shale

The ability of a material to resist failure under impact loading is called dynamic compressive strength. The maximum stress (peak stress) in the stress–strain curve, i.e., when the material is subjected to an impact load, is usually taken as the strength value.
The relationship between the dynamic compressive strength of shale and the strain rate is shown in Figure 13 by integrating the relationship between the peak strength and the average strain rate of the shale at normal temperature.
Concrete materials are sensitive to strain rates, and it is an important task to study the dynamic properties of concrete. Current studies generally agree that the dynamic tensile and compressive strength of concrete increases with increasing strain rates [37,38,39]. According to the peak strength of the shale under different strain rates, it can be found that the dynamic compressive strength of the shale has a strong strain rate effect, especially when the strain rates are between 10/s and 100/s; moreover, the peak strength of shale presents a linearly increasing trend with the increase in strain rates, which is similar to that of concrete materials, indicating that shale is also a rate-sensitive material.

5.4. Analysis of Shale under the Same Temperature and Different Strain Rates

Figure 14a–d show the stress–strain curves of shale under the different strain rates at 60 °C, 90 °C, 120 °C, and 150 °C, respectively. It can be seen that the peak stress and ultimate strain of the specimens increase with the increase in strain rates at the same temperature condition, showing the correlation. This phenomenon is called the strengthening effect of strain rate on the impact compression mechanical properties of shale materials.
However, due to the limited number of specimens, the maximum strain rate of shale at these temperatures did not exceed 100/s, and the stress–strain curves for high strain rates of shale at these temperatures were not able to be obtained, meaning that the experimental results are somewhat limited.

5.5. Analysis of Shale under the Same Pressure and Different Temperatures

Figure 15 shows the stress–strain curves of shale specimens under the same pressure at different temperatures. It can be seen from Table 3 and Table 4 that under the same impact pressure condition, within the temperature difference range of this experiment, the dynamic strength of the specimen presents two opposite trends (increase and decrease) with the increase in temperature, which is determined by the direction of bedding [10]. As shown in Figure 14a, when the strain rates are close, the strain of shale will increase with the increase in temperature, especially when the temperature is 150 °C, and the corresponding strain value is much larger than that at 60 °C, 90 °C, and 120 °C. This indicates that shale will soften when the temperature exceeds 150 °C, resulting in an increase in strain, with the corresponding elastic modulus and hardness decreasing. Through calculation, it can be found that the elastic modulus of shale at 150 °C is 2–3 times smaller than that at 60 °C, 90 °C, and 120 °C, indicating that temperature will greatly affect the elastic modulus of shale and, thus, affect the hardness of shale.

6. Dynamic Damage Constitutive Model

6.1. Definition of Damage Variable

Krajcinovic and Silva [40] started with the random distribution of flaws within the rock material and used the rock’s microelement strength to follow the Weibull distribution to build the statistical damage constitutive equation of the rock fracture process. A granite thermal–mechanical coupling damage model was developed by Xu et al. [41] using the Drucker–Prager criterion and the Weibull distribution. Zhu et al. [42] developed a statistical damage model using the modified Mohr–Coulomb criterion as the failure condition for the microelement strength. An enhanced statistical damage constitutive model of rock thermal treatment was put out by Wang et al. [43] and is based on the longitudinal velocity correction of the thermal damage variable.
As a typical non-homogeneous material, shale contains a large number of randomly distributed microcracks, pores, and other defects. Assuming that shale consists of many microbodies, whose dimensions are large enough in a spatial sense to contain many microdefects but small enough in a mechanical sense to be considered as particles, the microbodies have the following properties: (1) The shale is isotropic, i.e., having an isotropic-damaged body. (2) The shale conforms to Hooke’s law before failure but has no bearing capacity after failure. (3) The strength of each microbody follows the Weibull distribution. The expression of its probability density function is as follows [44]:
P ( F ) = m F 0 ( F F 0 ) m 1 exp [ ( F F 0 ) m ]
where F is the distribution variable of micro-body strength, and F 0 and m are Weibull distribution parameters, reflecting the mechanical properties of shale materials.
When shale is subjected to impact loading, microelements are gradually damaged, and a statistical damage variable, D can be introduced to represent this process. Under a certain load, D can be expressed as:
D = N f N
where N is the total number of microelements, and N f is the number of destroyed microelements. When loading to a specific load level, F :
N f = 0 F N P x d x
By substituting Equations (2) and (4) into Equation (3), the expression of the damage variable, D can be obtained as:
D = 1 exp [ ( F F 0 ) m ]

6.2. Micro-Body Strength

According to Equation (5), the damage variable, D is affected by the micro-body strength, F which is related to the stress state of the rock. To express the effect of complex stress states on shale, the failure criterion of the microelement is considered to be consistent with the Drucker-Prager (D-P) failure criterion. The microbody strength, F based on the D-P failure criterion is as follows [40]:
F = α I 1 + J 2
α   =   sin φ 9 + 3 sin 2 φ I 1   =   ( σ 1 + 2 σ 3 ) E ε 1 σ 1 2 μ σ 3 J 2   =   ( σ 1 σ 3 ) E ε 1 3 ( σ 1 2 μ σ 3 )
where α is the strength parameter of the microelement, and φ is the internal friction angle of the rock. The internal friction angle of the shale selected in this paper was 28°, and the value of α was 0.1268. I 1 is the first invariant of the stress tensor, and J 2 is the second invariant of the stress deviation. Substituting Equation (7) into Equation (6), considering the uniaxial state of the shale in this paper, so σ 2 = σ 3 = 0 and ε 1 = ε , then the expression of microelement strength is:
F = ( α + 1 3 ) E ε

6.3. Establishment and Modification of the Constitutive Model

Using the strain equivalence assumption, the expressions of the statistical damage intrinsic model of the intensity of the uniaxial shale under impact loading are obtained as follows:
σ = E ε ( 1 D ) = E ε exp [ ( F F 0 ) m ]
According to Equation (9), the key to establish the constitutive model is to determine the two unknown parameters F 0 and m . Considering that the peak stress and corresponding strain of the stress–strain curves are easily obtained in the uniaxial impact test, this paper used the extremum method to determine the parameters F 0 and m . Since the derivative of the multivariate function at the extreme value point ( ε m , σ max ) is 0, both sides of Equations (8) and (9) are derived and simplified simultaneously as follows:
F 0 = ( α + 1 3 ) E ε m m 1 m
m = 1 ln [ σ max / ( E ε m ) ]
The constitutive model parameters calculated using Equations (10) and (11) are shown in Table 5:
The relationship between F 0 , m and the strain rate can be established, which can effectively modify the failure statistical constitutive model of shale. Scatter plots were drawn with F 0 and m as the ordinates and strain rate as the abscissa, respectively, and nonlinear fitting was carried out, as shown in Figure 16 and Figure 17.
The fit shows the following relationship between F 0 , m and the strain rate:
F 0 = 1.70563 ε 2 + 65.6271 ε 482.51175
m = 0.61848 ε 2 25.25049 ε + 262.5967
The fitting correlation coefficients were 0.960 and 0.989, respectively. Equations (12) and (13) were substituted into Equation (9) to obtain the modified equation of the statistical intensity constitutive damage model of the shale under different strain rates:
σ = E ε exp [ ( 0.728 E ε 1.70563 ε 2 + 65.6271 ε 482.51175 ) 0.61848 ε 2 25.25049 ε + 262.5967 ]

6.4. Verification of the Constitutive Model

The theoretical dynamic stress–strain curves with close strain rates of the shale were calculated via Equation (14) using the modified constitutive model, and the compared results with the experimental curves are shown in Figure 18.
From the model validation results, the theoretical curves are consistent with the experimental curves at the same strain rate. Based on the effect of the Weibull distribution parameters F 0 and m , a reasonable correction of the theoretical curve can better reflect the insensitivity of the peak stress of the shale to temperature. The modulus of elasticity decreases with increasing temperature. Notably, the shale modulus of elasticity significantly decreases at 150 °C. Of course, the model established in this paper cannot fully express the changing process of the dynamic stress–strain curves of the shale at each loading stage, and there are some limitations due to the limited data available. The errors caused via the determination of the model parameters and the choice of elastic modulus need to be further improved in future studies.

7. Conclusions

Through the shale SHPB compression bar test, it was found that the shale showed three typical damage forms of “cracking and spalling”, “breaking and crushing” at different impact velocities, and the degree of damage increases with the increase in impact pressure.
The Hopkinson bar compression tests of shale at room temperature, 60 °C, 90 °C, 120 °C, and 150 °C, and at 0.2 MPa, 0.3 MPa, and 0.4 MPa were carried out. The stress–strain curves of shale at the same strain rate with different temperatures and the same temperature with different strain rates were obtained. Under the same impact pressure condition, in the temperature difference range of this experiment, the dynamic strength of the specimen presents two opposite trends (increase and decrease) with the increase in temperature, which was determined via the direction of bedding.
Based on the Weibull statistical distribution and the D-P failure criterion, a statistical damage constitutive model of the shale dynamic strength considering temperature and strain rate effects was obtained. By modifying the parameters F 0 and m , the established dynamic statistical intensity damage constitutive model of the shale is in good agreement with the experimental results.

Author Contributions

Conceptualization, W.G.; Methodology, G.D.; Validation, G.S.; Formal analysis, Y.D.; Investigation, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, W.K. Application of high-energy gas fracturing technology in oil and gas resource development. J. Earth Sci. Environ. 2000, 22, 60–62. (In Chinese) [Google Scholar]
  2. Zhang, D.Y.; Zhang, X.Y. Research and application of high energy gas fracturing technology. Oil Drill. Prod. Technol. 1997, 19, 17–23. (In Chinese) [Google Scholar]
  3. Wang, H. Study on Evaluation of Compressibility of High Energy Gas Fracturing for Horizontal Wells in Shale Gas Reservoirs; Xi’an Shiyou University: Xi’an, China, 2016. [Google Scholar]
  4. Horsrud, P.; Sønstebø, E.; Bøe, R. Mechanical and petrophysical properties of North Sea shales. Int. J. Rock Mech. Min. Sci. 1998, 35, 1009–1020. [Google Scholar] [CrossRef]
  5. Hsu, S.-C.; Nelson, P.P. Characterization of Eagle Ford Shale. Eng. Geol. 2002, 67, 169–183. [Google Scholar] [CrossRef]
  6. Sarout, J.; Molez, L.; Guéguen, Y.; Hoteit, N. Shale dynamic properties and anisotropy under triaxial loading: Experimental and theoretical investigations. Phys. Chem. Earth Parts A/B/C 2007, 32, 896–906. [Google Scholar] [CrossRef]
  7. Niandou, H.; Shao, J.; Henry, J.; Fourmaintraux, D. Laboratory investigation of the mechanical behaviour of Tournemire shale. Int. J. Rock Mech. Min. Sci. 1997, 34, 3–16. [Google Scholar] [CrossRef]
  8. Kuila, U.; Dewhurst, D.; Siggins, A. Raven Stress anisotropy and velocity anisotropy in low porosity shale. Tectonophysics 2011, 503, 34–44. [Google Scholar] [CrossRef]
  9. Luo, N.; Fan, X.; Cao, X.; Zhai, C.; Han, T. Dynamic mechanical properties and constitutive model of shale with different bedding under triaxial impact test. J. Pet. Sci. Eng. 2022, 216, 110758. [Google Scholar] [CrossRef]
  10. Huang, L.; Guo, Y.; Li, X. Mechanical response to dynamic compressive load applied to shale after thermal treatment. J. Nat. Gas Sci. Eng. 2022, 102, 104565. [Google Scholar] [CrossRef]
  11. Dong, D.Z.; Zou, C.N.; Yang, H.; Wang, Y.M.; Li, X.J.; Chen, G.S.; Wang, S.Q.; Lv, Z.G.; Huang, Y.B. Progress and prospect of shale gas exploration and development in China. Acta Pet. Sin. 2012, 33, 107–114. (In Chinese) [Google Scholar]
  12. Heng, S.; Yang, C.H.; Zhang, B.P.; Guo, Y.T.; Wang, L.; Wei, Y.L. Experimental study on anisotropy characteristics of shale. Rock Soil Mech. 2015, 36, 609–616. (In Chinese) [Google Scholar]
  13. Peng, G.Z. Relationship between structural plane orientation and mechanical properties of shale rock under uniaxial compressive stress. Chin. J. Geotech. Eng. 1983, 2, 101–109. (In Chinese) [Google Scholar]
  14. Liu, G.J.; Xian, X.F.; Zhou, J.P.; Zhang, L.; Liu, Q.L.; Zhang, S.W. Experimental study on deformation characteristics of shale under load and influence of mineral components on rock brittleness. J. China Coal Soc. 2016, 41, 369–375. (In Chinese) [Google Scholar]
  15. Kang, Z.Q. Simulation Study on Pyrolysis Characteristics of Oil Shale and in-Situ Thermal Injection for Oil and Gas Production; Taiyuan University of Technology: Taiyuan, China, 2008. (In Chinese) [Google Scholar]
  16. Deng, J.X.; Tang, Z.Y.; Li, Y.; Xie, J.; Liu, H.L.; Guo, W. Influence of diagenetic process on seismic petrophysical characteristics of shale in Wufeng-Longmaxi Formation. Chin. J. Geophys. 2018, 61, 659–672. (In Chinese) [Google Scholar]
  17. Wang, Y.; Xu, J.; Mei, C.G.; Ou, B.; Li, L. Study on physicochemical and mechanical properties of hard brittle shale with fractures. J. Oil Gas Technol. 2011, 33, 104–108+9. (In Chinese) [Google Scholar]
  18. Liu, J.X.; Liu, W.; Yang, C.H.; Huo, L. Experimental study on mechanical properties of mud shale under different strain rates. Rock Soil Mech. 2014, 35, 3093–3100. (In Chinese) [Google Scholar]
  19. Liu, J.X.; Zhang, K.; Liu, W.; Shi, X.L. Experimental study on deformation and failure characteristics of shale under different confining pressures and strain rates. Rock Soil Mech. 2017, 38, 43–52. (In Chinese) [Google Scholar]
  20. Zhang, H.J. Shale Dynamic Constitutive Relationship and Its Application; Southwest Petroleum University: Sichuan, China, 2014. (In Chinese) [Google Scholar]
  21. Yang, G.L.; Bi, J.J.; Zhang, Z.F.; Li, X.G.; Liu, J.; Hong, P.X. Effect of bedding angle under passive confining pressure on dynamic strength and energy dissipation of shale. J. Min. Sci. Technol. 2021, 6, 188–195. (In Chinese) [Google Scholar]
  22. Abbas, H.A.; Mohamed, Z.; Kudus, S.A. Deformation behaviour, crack initiation and crack damage of weathered composite sandstone-shale by using the ultrasonic wave and the acoustic emission under uniaxial compressive stress. Int. J. Rock Mech. Min. Sci. 2023, 170, 105497. [Google Scholar] [CrossRef]
  23. Fan, L.; Yi, X.; Ma, G. Numerical manifold method (NMM) simulation of stress wave propagation through fractured rock mass. Int. J. Appl. Mech. 2013, 5, 1350022. [Google Scholar] [CrossRef]
  24. Zhou, X.; Fan, L.; Wu, Z. Effects of microfracture on wave propagation through rock mass. Int. J. Geomech. 2017, 17, 04017072. [Google Scholar] [CrossRef]
  25. Fan, L.; Gao, J.; Du, X.; Wu, Z. Spatial gradient distributions of thermal shock-induced damage to granite. J. Rock Mech. Geotech. Eng. 2020, 12, 917–926. [Google Scholar] [CrossRef]
  26. Fan, L.F.; Wu, Z.J.; Wan, Z.; Gao, J.W. Experimental investigation of thermal effects on dynamic behavior of granite. Appl. Therm. Eng. 2017, 125, 94–103. [Google Scholar] [CrossRef]
  27. Baldin, M.T.; Salamuni, E.; Lagoeiro, L.E.; Sanches, E. High-temperature deformation processes in migmatites from the Atuba Complex, southern Ribeira Belt: Constrains from microstructures and EBSD analyses. J. South Am. Earth Sci. 2022, 116, 103863. [Google Scholar] [CrossRef]
  28. SAnsari, S.A.; Kazim, M.; Khaliq, M.A.; Ratlamwala, T.A.H. Thermal analysis of multigeneration system using geothermal energy as its main power source. Int. J. Hydrog. Energy 2020, 46, 4724–4738. [Google Scholar] [CrossRef]
  29. Zheng, Y.; Jia, C.; Zhang, S.; Shi, C. Experimental and constitutive modeling of the anisotropic mechanical properties of shale subjected to thermal treatment. Geomech. Energy Environ. 2023, 35, 100485. [Google Scholar] [CrossRef]
  30. Zhang, K. Effect of Anisotropy on Mechanical Properties and Strain Rate Effect of Shale; Southwest University of Science and Technology: Mianyang, China, 2017. (In Chinese) [Google Scholar]
  31. Liao, J.S. Study on Mechanical Properties of Shale under Uniaxial Lateral Confinement; Southwest University of Science and Technology: Mianyang, China, 2018. (In Chinese) [Google Scholar]
  32. Kolsky, H. An Investigation of the Mechanical Properties of Materials at Very High Rates of Strain. Proc. Phys. Soc. Sect. B 1949, 62, 676–700. [Google Scholar] [CrossRef]
  33. ANSI/B0019QEADC; Metals Handbook—Ninth Edition. American Society for Metals: Detroit, MI, USA, 1985.
  34. Chen, W.; Ravichandran, G. Dynamic compressive failure of a glass ceramic under lateral confinement. J. Mech. Phys. Solids 1997, 45, 1303–1328. [Google Scholar] [CrossRef]
  35. Jia, B. Study on Static and Dynamic Characteristics of Concrete at High Temperature; Chongqing University: Chongqing, China, 2011. (In Chinese) [Google Scholar]
  36. Lv, T.H. Study on Mechanical Behavior of Concrete and Reinforced Concrete under Impact Compression Based on SHPB; University of Science and Technology of China: Hefei, China, 2018. (In Chinese) [Google Scholar]
  37. Li, Q.M.; Meng, H. About the dynamic strength enhancement of concrete-like materials in a split Hopkinson pressure bar test. Int. J. Solids Struct. 2003, 40, 343–360. [Google Scholar] [CrossRef]
  38. Tedesco, J.W.; Ross, C.A.; McGill, P.B.; O’Neil, B.P. Numerical analysis of high strain rate concrete direct tension tests. Comput. Struct. 1991, 40, 313–327. [Google Scholar] [CrossRef]
  39. Wu, H.; Zhang, Q.; Huang, F.; Jin, Q. Experimental and numerical investigation on the dynamic tensile strength of concrete. Int. J. Impact Eng. 2005, 32, 605–617. [Google Scholar] [CrossRef]
  40. Krajcinovic, D.; Rinaldi, A. Statistical Damage Mechanics—Part I: Theory. J. Appl. Mech. 2005, 72, 76–85. [Google Scholar] [CrossRef]
  41. Xu, X.; Karakus, M. A coupled thermo-mechanical damage model for granite. Int. J. Rock Mech. Min. Sci. 2018, 103, 195–204. [Google Scholar] [CrossRef]
  42. Zhu, Z.; Tian, H.; Wang, R.; Jiang, G.; Dou, B.; Mei, G. Statistical thermal damage constitutive model of rocks based on Weibull distribution. Arab. J. Geosci. 2021, 14, 495. [Google Scholar] [CrossRef]
  43. Wang, Z.; Feng, C.; Wang, J.; Song, W.; Wang, H. An improved statistical damage constitutive model for rock considering the temperature effect. Int. J. Geomech. 2022, 22, 04022203. [Google Scholar] [CrossRef]
  44. Zheng, Y.; Shi, H.R.; Liu, X.H.; Zhang, W.J. Failure characteristics and constitutive models of coal and rock under different strain rates. Explos. Shock Waves 2021, 41, 45–57. (In Chinese) [Google Scholar]
Figure 1. Physical model of fracture propagation under high-energy gas fracturing [3].
Figure 1. Physical model of fracture propagation under high-energy gas fracturing [3].
Applsci 13 12288 g001
Figure 2. The flow of research.
Figure 2. The flow of research.
Applsci 13 12288 g002
Figure 3. Sampled shale.
Figure 3. Sampled shale.
Applsci 13 12288 g003
Figure 4. 1# to 5 # original samples and CT scans.
Figure 4. 1# to 5 # original samples and CT scans.
Applsci 13 12288 g004aApplsci 13 12288 g004b
Figure 5. The ultimate failure mode of specimens with a loading angle of 0°.
Figure 5. The ultimate failure mode of specimens with a loading angle of 0°.
Applsci 13 12288 g005
Figure 6. SHPB experimental device schematic diagram.
Figure 6. SHPB experimental device schematic diagram.
Applsci 13 12288 g006
Figure 7. Physical picture of the SHPB experimental device.
Figure 7. Physical picture of the SHPB experimental device.
Applsci 13 12288 g007
Figure 8. The incident waveform of different shapers.
Figure 8. The incident waveform of different shapers.
Applsci 13 12288 g008
Figure 9. Shale heating device.
Figure 9. Shale heating device.
Applsci 13 12288 g009
Figure 10. Verification of stress balance.
Figure 10. Verification of stress balance.
Applsci 13 12288 g010
Figure 11. Typical failure modes of shale.
Figure 11. Typical failure modes of shale.
Applsci 13 12288 g011
Figure 12. Stress–strain curves of shale at different strain rates.
Figure 12. Stress–strain curves of shale at different strain rates.
Applsci 13 12288 g012
Figure 13. Compressive strength of shale under different strain rates.
Figure 13. Compressive strength of shale under different strain rates.
Applsci 13 12288 g013
Figure 14. Stress–strain curves of shale under the same temperature and different strain rates.
Figure 14. Stress–strain curves of shale under the same temperature and different strain rates.
Applsci 13 12288 g014aApplsci 13 12288 g014b
Figure 15. Stress–strain curves of shale under the same pressure and different temperatures. (a) Stress–strain curves of shale at different temperatures under 0.2 MPa. (b) Stress–strain curves of shale at different temperatures under 0.3 MPa.
Figure 15. Stress–strain curves of shale under the same pressure and different temperatures. (a) Stress–strain curves of shale at different temperatures under 0.2 MPa. (b) Stress–strain curves of shale at different temperatures under 0.3 MPa.
Applsci 13 12288 g015
Figure 16. Relationship between m and strain rate.
Figure 16. Relationship between m and strain rate.
Applsci 13 12288 g016
Figure 17. Relationship between F0 and strain rate.
Figure 17. Relationship between F0 and strain rate.
Applsci 13 12288 g017
Figure 18. Experimental and theoretical stress–strain curves of shale at the same strain rate and different temperatures.
Figure 18. Experimental and theoretical stress–strain curves of shale at the same strain rate and different temperatures.
Applsci 13 12288 g018
Table 1. Material composition of shale.
Table 1. Material composition of shale.
NumberGeological AgeClayQuartzPlagioclaseCalciteDolomitePyriteOtherClay MineralNon-Clay MineralsLithology
Dz01S1137461---13763Gray shale
Dz02S1s4337732314357Gray calcareous shale
Dz03S1s4322712114357Gray calcareous shale
Lt01S1132563262-3268Black shale
Lt02S113845171713862Black shale
Lt03O3w3346-116113367Linxiang Formation of Upper Ordovician
Ltc1S11315152-263169Black shale
Table 2. Splitting tensile strength.
Table 2. Splitting tensile strength.
Specimen NumberFailure Load Value (kN)Splitting Strength (MPa)Average Value (MPa)
H1039.9918.7618.36
H1139.9118.72
H1237.5617.62
Table 3. Values of each group at 0.2 MPa.
Table 3. Values of each group at 0.2 MPa.
Temperature60 °C90 °C120 °C150 °C
Peak stress112 MPa121 MPa123 MPa118 MPa
Ultimate strain0.00380.00330.00420.0075
Table 4. Values of each group at 0.3 MPa.
Table 4. Values of each group at 0.3 MPa.
Temperature60 °C90 °C120 °C150 °C
Peak stress151 MPa109 MPa146 MPa115 MPa
Ultimate strain0.00730.00680.00610.0056
Table 5. Results of the constitutive model parameter calculation.
Table 5. Results of the constitutive model parameter calculation.
TemperatureStrain Rate (s−1) F 0 (MPa) m
60 °C20.63141.7674.902
90 °C19.33153.2235.599
120 °C20.33144.4134.878
150 °C23.94112.60812.565
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, W.; Deng, G.; Sun, G.; Deng, Y.; Li, Y. Dynamic Mechanical Properties of Heat-Treated Shale under Different Temperatures. Appl. Sci. 2023, 13, 12288. https://doi.org/10.3390/app132212288

AMA Style

Gao W, Deng G, Sun G, Deng Y, Li Y. Dynamic Mechanical Properties of Heat-Treated Shale under Different Temperatures. Applied Sciences. 2023; 13(22):12288. https://doi.org/10.3390/app132212288

Chicago/Turabian Style

Gao, Weiliang, Guoqiang Deng, Guijuan Sun, Yongjun Deng, and Yin Li. 2023. "Dynamic Mechanical Properties of Heat-Treated Shale under Different Temperatures" Applied Sciences 13, no. 22: 12288. https://doi.org/10.3390/app132212288

APA Style

Gao, W., Deng, G., Sun, G., Deng, Y., & Li, Y. (2023). Dynamic Mechanical Properties of Heat-Treated Shale under Different Temperatures. Applied Sciences, 13(22), 12288. https://doi.org/10.3390/app132212288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop