Variation Law of Infrared Radiation Temperature of Unloading Fracture of Composite Coal-Rock

For composite mining coal-rock dynamic disaster, combining the theory of thermodynamics, damage mechanics, and other disciplines, the thermodynamic coupling mathematical model of composite coal-rock under an unloading condition is deduced, and the simulation model of composite coal-rock is established for numerical simulation. And the variation law of thermal infrared radiation under triaxial loading and unloading of composite coal-rock is analyzed and veri ﬁ ed by experiments. The results show the following ﬁ ndings: (1) The distribution of thermal infrared radiation temperature of composite coal-rock is di ﬀ erent in di ﬀ erent stages of stress. The overall temperature of the temperature ﬁ eld is lower than the initial temperature ﬁ eld in the three-dimensional stress loading stage and the stage of stress-keeping pressure, but the internal temperature of the coal body is the highest. In the ﬁ rst stage of “ loading and unloading, ” the temperature of a coal seam increases slightly, and the temperature of other parts of the rock layer increases except for the circular low-temperature zone. In the second stage of “ loading and unloading, ” an alternating zone of high and low temperatures appears in the rock mass, and the temperature ﬁ eld is enhanced, among which the temperature ﬁ eld reaches the strongest after unloading the con ﬁ ning pressure. After jumping over the maximum stress, the temperature ﬁ eld decreases as a whole at the instability and rupture stage. (2) The variation of surface average thermal infrared radiation temperature ( T ave ) of composite coal-rock can be divided into the initial ﬂ uctuation stage, the linear heating stage, the local decline stage, the temperature sudden increase stage, and the fracture decline stage. At three di ﬀ erent unloading rates of 0.003MPa/s, 0.03MPa/s, and 0.05MPa/s, the T ave of coal body, ﬂ oor rock, and roof rock reach the maximum before composite coal-rock instability and fracture, and the temperature change of the coal body is the most obvious. (3) Under di ﬀ erent con ﬁ ning pressure unloading rates, the T ave of roof rock, coal body, and ﬂ oor rock shows a strong linear relationship with stress after linear ﬁ tting. And the correlation between simulation and experimental results after ﬁ tting is above 0.89. The larger the con ﬁ ning pressure unloading rate is, the shorter the peak time of T ave arrives, and the larger the peak value. The comparison between the experimental results and the simulation analysis shows that the two results are consistent, and the research results can provide a theoretical basis for the prevention and control of dynamic disasters in coal and rock mining.


Introduction
With the increasing depth of coal mining, coal and rock dynamic disasters such as coal and gas outburst and rock burst are becoming more and more serious [1]. Reference [2] shows that nearly half of coal mine disasters in China have been found to occur in groups or accompanied by nearby earthquakes, and all types of disasters are involved. The stress disturbance appears to exist between mining areas. The process of coal mining is an unloading process essentially. This mechanical process will cause a lot of thermal infrared radiation and temperature change in coal and rock mass [3][4][5][6]. At the same time, the temperature change during the mining process will further affect the physical parameters and mechanical properties of the material, thus making the actual unloading process more complex. At present, most studies on the mechanical properties of coal and rock are carried out under uniaxial loading conditions. Reference [7] studied the brittle failure model of coal-rock composites under uniaxial compression through laboratory tests and numerical simulation. Reference [8] analyzed stress characteristics and mechanical properties and compared the similarities and differences of composite coal mass with soft and hard coal mass in deformation and failure characteristics. Moreover, the mechanical properties and deformation and failure laws of the composite coal mass under unidirectional loading were investigated. Reference [9] studied the mechanical properties of coal-like rocks and the law of damage evolution under different loading rates. The results show that the higher the loading rate, the earlier the specimen enters the stage of stable development of damage. In [10], the particle flow method was used to establish the composite model of coal and rock with different inclination angles, and the impact instability characteristics after unloading were analyzed.
Relevant results at home and abroad show that by studying the distribution and variation rule of thermal infrared radiation in composite coal-rock, the dynamic disaster of coal and rock can be effectively predicted [11][12][13][14]. Reference [15,16] used an infrared thermal imager to conduct a simulation test on the failure process of rock strata caused by excavation, and found that infrared radiation temperature can reflect the change of rock strata stress at different depths, while an infrared thermal image can reflect the failure process of rock mass. Reference [17] believes that infrared precursors are mainly manifested as infrared thermography and average infrared radiation temperature anomaly. Reference [18] carried out the infrared radiation observation test on the stress disaster of granite, studied the mechanism of the infrared anomaly phenomenon before the earthquake, and found that the infrared radiation temperature on the outer surface of the rock also changed with the constant change of the stress. Reference [19] carried out uniaxial compression loading tests on aqueous concrete samples, and tested and studied the variation characteristics of infrared radiation temperature, infrared thermography, and acoustic emission, as well as the precursor information of fracture. Reference [20] used infrared thermal imaging technology to study the changes in infrared radiation during coal-rock deformation. The research showed that the infrared radiation temperature on the coal and rock surface would be changed with the increase of the loadings. There was a stage differential in infrared radiation temperature variation in the cracking-developed zone, stress-concentrated area, and broken zone. Reference [21] carried out an observation test of infrared radiation characteristics in the process of loading failure of coal samples, discussed the infrared precursor information of instability failure of coal samples with different bursting liabilities, and obtained the infrared radiation temperature precursor characteristics of coal samples under the condition of no flushing, weak flushing, and strong flushing. Reference [22] monitored the change of an infrared radiation temperature field in surrounding rock of a roadway through infrared radiation technology, and verified that the radiation temperature field had a good corresponding relationship with roadway failure characteristics. Reference [23] found that infrared radiation changes would occur along with stress adjustment in the process of mechanical fracture of coal and rock. Reference [24] quantitatively analyzed the characteristic parameters of infrared radiation temperature and the evolution law of coal and rock damage in the process of coal and rock loading. Reference [25] studied the internal infrared radiation temperature change law of composite coalrock under load. The results show that the composite coalrock has thermoelastic effects and thermal friction effects during the loading process. In addition, the heat generated by the displacement and friction between the fissure surfaces of the coal is significantly higher than the heat generated by the thermoelastic effect, and the thermoelastic effect in the rock mass plays a leading role. Reference [26] found that the thermal infrared radiation temperature temporal evolution in the process of a fracture of coal and rock can be divided into four stages, namely, the T ave calm stage, the T ave slight fluctuation stage, the T ave rapid increase stage, and the T ave rapid decrease stage, which accurately described the infrared radiation law of coal and rock under uniaxial loading and unloading state. Reference [27] found in their study that the infrared radiation precursor of coal-rock fracture and instability showed a trend of step and type spurt upward. In summary, predecessors have studied the characteristics of thermal infrared radiation temperature and precursor laws of coal and rock under uniaxial loading, but there are few related studies on unloading failure of deep coal and rock mass.
Based on the above research results, related scholars have carried out related research on unloading fracture of coal and rock and multifield coupling research. However, there are few studies on the thermal-mechanical coupling mechanism of anisotropic coal and rock masses, and the research on thermal infrared radiation under the condition of triaxial unloading of composite coal-rock is not thorough enough. Reference [28] conducted a three-axis unloading compression test based on the pressure path of the underground engineering excavation. The research results show that the infrared radiation characteristics of rock specimens under different intermediate principal stresses are different, revealing that as the intermediate principal pressure increases, the failure mode of the rock is also different. Reference [29] conducted conventional triaxial compression (CTC) tests on the gas-bearing coal, gas-bearing coal-mudstone combination, and gas-bearing coal-sandstone combination using the RLW-500G triaxial experimental system. Reference [30] conducted laboratory tests under two stress paths, namely, conventional triaxial compression and prepeak unloading confining pressure, in order to explore the failure mechanism of the unloading side of coal roadway excavation, and analyzed its mechanical characteristics. According to the damage and failure law, the results show that under the unloading stress path, the coal body failure is a composite state in which shear failure is the main part, but there is still partial split failure. Based on the unloading mechanics path of the fixed axial displacement and the unloading surrounding pressure, the study on the damage features during the unloading process was conducted on the reconstituted coal with high plasticity with CT value, CT damage, and CT image of each scanning layer obtained from CT imaging experiments [31]. Reference [32] established a temperature-stress numerical coupling model of rock at high temperature, carried out numerical simulation on the evolution law of the temperature field and fracture field of overburden rock in a combustion cavity under the temperature-stress coupling condition, and studied the 2 Geofluids relationship between stress and temperature in the state of rock thermal damage. Reference [33] carried out physical experiments and numerical simulations on the system instability faults of gas coal-bearing rock under different circumstances, comprehensively analyzed the multiphysical coupling mechanism of gas coal-bearing rock damage and the occurrence conditions of outbreak-rock coupling dynamic disaster, and improved the understanding of the mechanism of outbreakrock coupling dynamic disaster. Reference [34] conducted an unloading coal and rock destruction experiment under the multiphase coupling action of water, force, and heat. By optimizing the noise reduction and the infrared radiation information of three-dimensional fusion, the three-dimensional infrared radiation field was reconstructed. And by considering the surface damage characteristics, the three-dimensional infrared radiation field and simulation software results reveal the spatiotemporal evolution of the infrared radiation information of the unloading damage of composite coal-rock, and realize the noncontact prediction of coal and rock dynamic disasters. Reference [35] proposed a thermal-mechanical coupling method for large-volume concrete on the inner wall of deep wells, considering various boundary conditions such as internal heat source, external cold source, plastic plate friction constraint, overlying shaft wall load, and formation lateral pressure. The early temperature-stress field evolution characteristics of the inner wall of the deep frozen wellbore were obtained, and the high-risk areas prone to early temperature cracks were determined. The abovementioned research results have made a very useful exploration of thermal infrared radiation precursors before unloading and rupture of coal and rock, and laid the foundation for subsequent research.
In summary, a large number of coal and rock mining dynamic disaster events show that before the occurrence of dynamic disaster events, there will be various degrees of precursor information, such as changes in stress and strain, temperature, and acoustic emission. By monitoring the changing trends and abnormal values of these signals, it is possible to provide effective early warning of shock disaster hazards, and timely effective control measures can be taken to reduce the casualties and equipment damage caused by rock shock disasters, which has engineering application prospects. However, most of the current researches on the prediction and prediction of coal and rock mining dynamic disasters are based on uniaxial conditions, which is far from the actual three-axis mining situation. The research on the infrared radiation temperature law based on the three-axis unloading condition needs to be in-depth. There is still much room for improvement in analyzing the unloading rupture mechanism from the perspective of multiphysics coupling. There are few studies on the thermodynamics of anisotropic coal and rock masses. Aiming at the problem of effectively predicting and forecasting coal and rock mining dynamic disasters, this paper intends to study not only the distribution of thermal infrared radiation temperature during the unloading and fracture process of anisotropic composite coal-rock, but also the relationship between thermal infrared radiation and stress on the surface of composite coal-rock and its changing law based on thermodynamics and damage mechanics. Joint early-warning simulations and experiments on coal and rock instability through multiphysics fields to explore the characteristics of infrared radiation temperature changes and coal and rock rupture precursor laws have important theoretical significance for the development of coal and rock mining dynamic disaster prediction and forecast.

Materials and Methods
2.1. Thermal-Mechanical Coupling Model of Composite Coal-Rock Unloading. According to the Stefan-Boltzman Law, the thermal infrared radiation energy of the coal and rock surface is proportional to the power of the physical temperature and emittance. When the material type and surface finish of coalrock are constant, the thermal infrared radiation intensity received by the surface is only related to the physical temperature of coal-rock [26]. Therefore, during the loading-unloading process, when the thermal infrared radiation intensity of the composite coal-rock under the action of stress changes, the reason is that the external stress causes the change of the thermal field of the coal-rock, which leads to the alienation of the temperature field of thermal infrared radiation on the surface of the coal-rock. This thermal change is actually a thermodynamic coupling effect, in which the thermal stress will produce thermal damage and other effects and accelerate the instability failure of the coal and rock mass [36].
According to the Elastic Damage Theory, the equivalent elastic modulus E of coal-rock after damage can be expressed as follows: where E 0 is the initial elastic modulus before damage; D is the damage variable, which is used to describe the damage degree and evolution of coal-rock under the action of external load. Taking the interior elements of coal and rock mass as the research object, the magnitude of the three principal stresses in the loading-unloading process is not zero [37], which conforms to the generalized Hooke law on the whole. The principal strain equation is as follows: where ε i (i = x, y, z) is the principal strain; E is the elastic modulus; σ i (i = x, y, z) is the principal stress; μ is the Poisson ratio.
The shear strain equation is as follows: where γ ii (ii = xy, yz, zx) is the shear stress; τ ii (ii = xy, yz, zx) is the shear strain; G = E/ð2ð1 + μÞÞ is the Lame constant.
Equation (2) and equation (3) jointly express the stressstrain relationship of composite coal-rock under unloading conditions. In a three-dimensional elastic system, strain can be characterized by principal strain and shear strain components. According to the principles of thermodynamics and reference [38], the relationship between temperature changes and strains of each component of the composite coal-rock is further deduced as follows: where T 0 is the initial temperature of different components of the composite coal-rock (°C); V is the volume of the different components of the composite coal-rock (m 3 ); C is the specific heat capacity of the different components of the composite coal-rock (J/(kg·K)); β = ð∂σ i /∂T 0 Þ ε , (i = x, y, z) or ð∂τ ii /∂T 0 Þ ε , (ii = xy, yz, zx) is the thermal stress coefficient. Then, the temperature field T of different components of composite coal-rock at any time is as follows: Substitute equation (1), equation (2), and equation (3) into equation (4). And according to equation (5), the specific relationship between temperature field and strain of each component of composite coal-rock can be obtained as follows: In order to establish the thermodynamic coupling model, the heat conduction equation should be added to determine the distribution of temperature field T in composite coalrock. According to the principle of heat balance in thermodynamics, it is assumed that the specific heat capacity C and thermal conductivity λ i (i = x, y, z) of composite coal-rock are constants that do not change with temperature. Introduce the solid three-dimensional anisotropic heat conduction differential equation [39] to describe the internal temperature field of composite coal-rock. The specific differential equation is as follows: where ρ is the density of material (kg/m 3 ); λ x , λ y , and λ z are the thermal conductivity coefficients of the material in three directions (W/(m·K)); f ðtÞ is the rate of heat release. However, when coal and rock are damaged under external force, the damage will lead to the increase of thermal conductivity and specific heat capacity of coal-rock. It is assumed that the influence of damage D on the specific heat capacity C and thermal conductivity λ i (i = x, y, z) of composite coalrock meets the following functional relationship [40]: where C 0 and C are the specific heat capacity of composite coal-rock of different components before and after damage (J/(kg·K)). λ 0 and λ i (i = x, y, z) are the thermal conductivity before and after damage of coal and rock of different components (W/(m·K)). α C and α λi are the influence coefficients of damage on the specific heat capacity and thermal conductivity of different components of composite coal-rock, respectively.
By substituting equation (8) into equation (6) and equation (7), the temperature of different components of composite coal-rock and the instantaneous temperature of each point in coal and rock can be calculated. The average infrared radiation temperature (T ave ) refers to the average value of the infrared radiation temperature corresponding to each point on the coal and rock surface at a certain moment, which reflects the variation characteristics of the infrared radiation temperature field on the coal and rock surface on the whole [25]. The specific equation is as follows: where T i is the radiation temperature of the ith point in the temperature field (°C). T ave is the average of T i (°C). n is the number of selected composite coal-rock points.
To sum up, the coupling mathematical model of the thermal infrared radiation temperature field and the stress field in the process of unloading and fracture of composite coal-rock is as follows:  [41], and the constitutive model in the thermal mode is the anisotropic heat conduction model. The material property parameters of the simulation model are shown in Table 1.
According to the "Coal Mine Safety Regulations," the working surface temperature does not exceed 26°C [42], and the temperature conditions are those that can be achieved by indoor experiments. Set the initial temperature and ambient temperature of the composite coal-rock simulation model in thermal mode to room temperature 20°C. The steps of the three-axial load-unloading simulation in the mechanical mode are as follows: (1) (1) and (2), triaxial loading is carried out again. When the stress reaches 80% of the peak intensity recorded in step (3), under the condition of ensuring that the simulation results are not distorted, three commonly used slow unloading confining pressure rates of 0.003 MPa/s, 0.03 MPa/s, and 0.05 MPa/s [43,44] are selected for unloading simulation until the composite coal and rock rupture, and the simulation ends.

Distribution Law of Thermal Infrared Radiation
Temperature. In order to better study the relationship between the unloading fracture mechanics stage of coal-rock and its infrared radiation temperature, the stress change process must first be clarified. The stress-time (σ-t) curve of triaxial loadingunloading under the unloading rate of 0.05 MPa/s confining pressure is shown in Figure 1. It is divided into five stages. OA is the three-dimensional stress loading stage. At this stage, due to the combined action of axial pressure and confining pressure, the primary voids and tiny cracks in the composite coal are gradually compacted, and the internal structure begins to change. This stage is similar to the compaction stage in uniaxial loading. AB is the stage of stress-keeping pressure. This stage is a process in which the axial pressure and the confining pressure resist each other, and the axial pressure presents a nonlinear change. BC is the first stage of "loadingunloading," which is elastic. The stress shows a linear increase trend. CD is the second stage of "loading-unloading," which is the process of converting elasticity to plasticity. The stress at this stage shows a transition from linear to nonlinear. DE is the instability and rupture stage. After unloading the confining pressure at a rate of 0.05 MPa/s at 142.5 s, until the peak stress at point D reaches the maximum value of 45.45 MPa, the composite coal-rock fractures, the curve drops rapidly, and the simulation ends.
In the simulation process, there is no need to slice the composite coal-rock simulation model, that is, the distribution of T ave can be observed through the overall infrared radiation temperature distribution cloud map. Figure 2 is a cloud map of the overall infrared radiation temperature distribution of the composite coal-rock under the confining pressure unloading at a rate of 0.05 MPa/s. At each stage of stress and after unloading the confining pressure, a time point is selected for analysis at 25 s, 60 s, 90 s, 130 s, 146 s, and 175 s, respectively. It can be seen from Figure 2 that the temperature distribution is different at different stages. However, due to the coal body's low density and low thermal conductivity, it has a large degree of force deformation and high-temperature accumulation. But the coal body has slow heat transfer and small temperature changes, and finally presents the characteristics of hightemperature distribution in the coal body.
From the perspective of different stages, the overall temperature of the temperature field during the three-dimensional stress loading (Figure 2(a)) and stress-keeping pressure stage (Figure 2(b)) is lower than the initial temperature. In the first stage of "loading-unloading" (Figure 2(c)), as the effect of axial stress becomes more obvious, the temperature of the coal body rises slightly. Except for the circular low-temperature zone, the temperature of the rest of the rock formation increases. In the second stage of "loading-unloading" (Figure 2(d)), the hightemperature area of the coal body spreads up and down on both sides of the rock mass. And the temperature field of part of the rock mass is significantly enhanced compared with the previous stage, and there is an alternating phenomenon of high-and lowtemperature regions in the lateral direction. After the confining pressure is relieved (Figure 2(e)), the temperature distribution of the temperature field expands according to the trend of the previous stage, and the temperature field reaches the strongest. After jumping over the highest point of stress, the coal-rock loses stability and ruptures, and the temperature field decreases during the instability and rupture stage (Figure 2(f)).

Changes of Thermal Infrared Radiation Temperature and Stress
(1) Study on Different Unloading Confining Pressure Rates. Figure 3 shows the stress-time (σ-t) curve of composite 1450 2400 Specific heat capacity (J/(kg·K)) 1200 920 Thermal expansion coefficient (°C) 6.435×10 -6 3.9×10 -5 Thermal conductivity (W/(m·K)) 0.258 3.081 5 Geofluids coal-rock triaxial loading-unloading simulation under different unloading confining pressure rates. By comparison, when the confining pressure and other conditions are the same and only the unloading rate is changed, the stress stages remain unchanged, and can still be divided into five stages. The overall stress trend is the same, and the previous loading curves are completely coincident. However, with the increase of unloading rate, the stress peak point of composite coal-rock fracture is reduced, the time from loading to fracture becomes shorter, and the overall stiffness and strength of composite coal-rock decreases, making it easier to fracture.
It can be seen from the foregoing that T ave can reflect the average radiation intensity of the infrared radiation released from the surface of the composite coal-rock, so this paper uses T ave as the analysis parameter for research. Figure 4 shows the average thermal infrared radiation temperature-time (T ave -t) curve of the coal surface under the condition of different rates of unloading confining pressure. It can be seen from Figure 4 that the three curves start to change from 142.5 s after unloading. And the peak arrival time of T ave at 0.05 MPa/s unloading confining pressure is earlier than that at 0.03 MPa/s and 0.003 MPa/s unloading confining pressure. The peak temperature is 22.49°C, which is greater than the 22.28°C and 22.03°C of the latter two. That is to say, under the same conditions, the greater the confining pressure unloading rate, the shorter peak time of T ave and the greater the peak value.
(2) The Law of Thermal Infrared Radiation Temperature Change. Figure 5 shows the T ave -σ-t curves of different components of composite coal-rock under the above three different confining pressure unloading rates, and the trends of the three are consistent. According to the average thermal infrared radiation temperature change curve of the coal surface, the temperature change process is divided into five stages: the initial fluctuation stage, the linear heating stage, the local decline stage, the temperature sudden increase stage, and the fracture decline stage. Take Figure 5(c) as an example for analysis. The five stages are as follows: (1) Initial fluctuation stage: T ave of different components of composite coal-rock shows a fluctuating trend of decreasing, increasing, and then decreasing to varying degrees. The difference is that the wave peak of the coal body appears later than the floor rock and the roof rock in sequence, and the T ave fluctuation degree is larger, and the maximum fluctuation is 0.38°C. At the same time, combined with Figure 1, it can be seen that, compared with the three-dimensional stress loading stage (OA), the internal stress field of the composite coal-rock gradually stabilizes after entering the stage of stress-keeping pressure (AB), and the degree of T ave fluctuation becomes smaller The instability and rupture stage The second stage of "loadingunloading" The first stage of "loadingunloading" The stage of stress keeping pressure   The first stage of "loadingunloading" The second stage of "loadingunloading"

The Correlation between Thermal Infrared Radiation
Temperature and Stress. The surface average infrared radiation temperature and stress of the composite coal body, roof rock, and floor rock are fitted to obtain the linear fitting relationship expression between the average thermal infrared radiation temperature and stress: where T ave is the surface average thermal infrared radiation temperature of each component of the composite coal-rock at any time (°C), K is the slope of the fitting curve, σ is the stress (MPa), T 1 is the starting temperature of the linear heating stage (°C). Figure 6 shows the fitting curves of the T ave and stress on the surface of the composite coal-rock with different components under different unloading rates. The expression of the linear relationship between the T ave and stress on the coal surface in Figure 6(c) is as follows: where the linear correlation coefficient obtained by fitting is 0.924.
A linear fitting is performed on the T ave and the axial stress of the roof rock surface, and the linear relationship between the T ave and the stress is obtained as follows: where the linear correlation coefficient obtained by fitting is 0.93553. A linear fitting is performed on the T ave and the axial stress of the floor rock surface, and the linear relationship between the T ave and the stress is obtained as follows: where the linear correlation coefficient obtained by fitting is 0.92295. According to the fitting parameters in Table 2 (unloading confining pressure at a rate of 0.003 MPa/s and 0.03 MPa/s), the linear correlation coefficients of the T ave and stress after fitting are all greater than 0.89. So, there is a strong linear correlation between the average thermal infrared radiation temperature and the stress on the surface of composite coal-rock under different working conditions. Figure 7 is a schematic diagram of the triaxial test equipment and field testing. The test system consists of a triaxial test machine, dynamic data acquisition, power control cabinet, and computer. In order to ensure that temperature data can be collected in a closed environment, an observation window is set at the confining cylinder of the test machine. This article uses the TAW-  10 Geofluids 2000 microcomputer-controlled high-temperature rock triaxial testing machine for testing. This equipment can collect, store, process, and display test data and test curves in real time during the test. The average infrared radiation tem-perature of the sample is collected by the ThermoView™ Pi20 infrared thermal imager; that is, the average infrared radiation temperature corresponds to all the pixels on the thermal infrared image of the sample surface at a certain moment. The spectral response of the infrared thermal imager is 8~14 μm, and the sensitivity is 0.03°C. The coal and rock samples were selected from the coal seams with strong impact at Feng Shuigou Coal Mine in Pingzhuang, Inner Mongolia. In the actual coal mine, the coal seam is mostly sandwiched in the rock layer, which basically presents the three-layer structure of "rock-coal-rock." According to this characteristic, the composite coal-rock samples used in the test are processed [45]. Through field exploration and comparison, materials were taken from the field, and fresh coal samples with both roof and floor siltstone were drilled. In the laboratory, the roof rock, coal sample, and floor rock were cut and bonded at a ratio of 1 : 1 : 1 into 20 cylindrical samples with a diameter of 50 mm and a height    14 Geofluids

Experimental Verification.
The experiment steps are as follows: (1) At room temperature, seal the composite coal-rock sample with a thermoplastic tube and put it in a confining cylinder (2) Connect the thermal imaging camera to the computer, aim the front of the thermal imaging camera at the observation window to the observation surface of the sample, turn on the power, open the software to connect the thermal imaging camera, and check the thermal imaging camera status (3) Open the nitrogen valve of the testing machine to inject nitrogen, and close the gas valve when it is full (4) Apply three-dimensional stress to the sample to a predetermined confining pressure value of 10 MPa (5) Pretest: keep the confining pressure unchanged, load the sample under axial compression at a displacement rate of 0.05 mm/min. Record the peak stress intensity, and obtain the average value of the peak stress intensity of the five samples (6) Triaxial loading-unloading test: divide samples 1 #~1 5 # into three groups equally, repeat steps (1), (2), (3), and (4), set the same triaxial loading conditions as the pretest. When the stress reaches about 80% of the average value of the peak strength, the confining pressure will be relieved. And the unloading rate is set to 0.003 MPa/s (samples 1 # , 2 # , 3 # , 4 # , and 5 # ), 0.03 MPa/s (samples 6 # , 7 # , 8 # , 9 # , and 10 # ), and 0.05 MPa/s (samples 11 # , 12 # , 13 # , 14 # , and 15 # ), keep the axial pressure to continue loading during the unloading process. When the axial stress exceeds the peak stress and the sample is damaged, immediately stop the confining pressure unloading operation. Keep the confining pressure unchanged, and continue to apply the axial stress until the stress curve no longer changes with the change of strain, and the test ends In the pretest, the average value of unloading confining pressure is 39 MPa. Since the samples are randomly grouped, the first samples, namely, sample 1 # , 6 # , and 11 # , are selected for comparison in the three groups of different confining pressure unloading rate tests. Figure 9 shows the test results of samples 1 # , 6 # , and 11 # . It can be seen from Figure 9 that the cracks of the sample mainly exist in the coal seam and show an inverted "Y-shaped" splitting fracture. And the rock mass has microcracks, but these are not obvious. This is due to the relatively high strength of the sample, the more prominent end effect, and the end is in a three-way stress state due to friction restraint.
The T ave -σ-t curves of samples 1 # , 6 # , and 11 # are shown in Figure 10, and the trends of the three curves are the same. Combining Figures 9 and 10, the greater the confining pressure relief rate, the greater the T ave peak, the more severe the mechanical changes caused by heat, the more obvious the thermal-mechanical coupling effect, and the shorter the sample rupture time. Therefore, the cracks in the sample are too late to fully expand, and the more the damage is concentrated on a certain macroscopic fracture surface, the lower the overall degree of fracture. During the loading-unloading process, the formation, expansion and fusion of microcracks within the coal sample occur. As the stress increases, the T ave shows an upward trend, and local peak value and global maximum value appear. Since the strength of the roof and floor rocks is much greater than the coal body, the coal body was unstable and fractured at the end of the test. Although cracks propagated inside the rock mass, the rock mass itself did not break obviously. The T ave of the roof and floor rocks showed an upward trend. After the confining pressure was relieved, the The symbol "#" constitutes the label of the "coal-rock-coal" sample, which is explained in Section 3.1.4.  16 Geofluids T ave reached the maximum value within a short time. Take sample 11 # in Figure 10(c) as an example for analysis. For the coal body, in the initial fluctuation stage, the composite coal-rock, which is an anisotropic and heterogeneous material, has microdefects in the internal framework particles, and the internal microcracks and micropores gradually rupture and close at the initial stage of three-dimensional stress loading. Under the combined action of slight dislocation friction (the apparent damage has not yet begun) and the gas discharged from the pores, coal sample damage begins to occur, and there is a certain degree of heat dissipation. T ave decreases by 0.19°C within 0~19.7 s. After entering stage of stress-keeping pressure, the internal stress field of the coal sample gradually stabilizes, and the damage is not severe at the beginning. At this time, the sample deforms in the elastic range, and the thermoelastic effect plays a leading role. After compaction, the particles in the coal body have a small range of friction, and the heat released causes the T ave fluctuation to rise again by 0.39°C. In the linear heating stage, with the increase of stress, new cracks in the coal sample begin to develop, and microfractures appear and continue to develop. These microcracks continue to rub and heat up during the formation process and are accompanied by friction between the coal and rock interface. For heat production, T ave showed a linear  Figure 12: The experimental fitting curve of thermal infrared radiation temperature and stress. Under the combined action of shear cracks leading to heat generation and local thermal stress-induced local damage, the internal microcracks of the coal sample increase, and the external microcracks appear, which increases the sample volume and expands and absorbs heat. T ave shows a trend that first decreases from the local maximum and then stabilizes, decreasing by 0.43°C within 121.5 s to 152.2 s. During the sudden temperature increase stage, since the confining pressure began to be unloaded at 152.2 s, the confining pressure constraints on the sample gradually decreased, and the development of cracks and frictional slippage in the coal body became more and more intense. The deformation of the sample accelerates, the frictional heat effect occurs, and a large amount of heat is generated. T ave increases suddenly in a short time, reaching a global maximum of 20.81°C at 156.9 s. At the same time, the expansibility of the mineral components in the coal gradually becomes prominent. The superimposed effect induces thermal stress inside the specimen, leading to aggravation of microfractures, the increase in volume is more significant, and some heat is absorbed, and T ave begins to show a downward trend again. In the fracture decline stage, after 180.7 s, due to the further development of original cracks and the production of a large number of new cracks, they clustered towards the main rupture. When the stress reaches the maximum value of 45.46 MPa, the sample reaches the loadbearing limit, the cracks penetrate macroscopically in the coal body, and the sample ruptures. There was violent dislocation friction between the fractured surfaces, resulting in a large number of high-temperature points, and T ave had a small step increase, but the overall trend was still decreasing.
For the roof rock, with the initial load gradually increasing, the cracks and pores are in the compaction stage, the fracture surface does not have frictional heating effect, and the temperature fluctuates up and down. During the linear heating stage, a small amount of microscopic cracks also appeared in the roof rock. The T ave curve showed a nearly linear upward trend, but there was a buffer between 80 and 100 s. T ave lags behind the coal body, reaching a local maximum of 20.28°C at 129.8 s, which increased by about 0.14°C. Due to the increase and expansion of internal microfractures, T ave starts to decrease from the local maximum, and reaches the global maximum at 20.70°C at 161.6 s later than the coal body after the confining pressure is unloaded and is lower than the maximum value of the T ave of the coal body.
The infrared radiation temperature change of the floor rock is similar to that of the roof rock. The difference is that the floor rock bears the gravity from the roof rock and the coal body. During the three-way stress loading stage and the stage of stress-keeping pressure, more gas is discharged from the internal voids, and the phenomenon of heat absorption phenomenon is obvious. The fluctuation trend of T ave is lower than that of the roof rock. In the linear heating stage, the temperature of the roof rock changes more obviously, the temperature is increased by 0.42°C. And the occurrence time of the local maximum is between the coal body and the roof rock, reaching 20.31°C at 125.03 s. Similarly, the global maximum reaches 20.73°C at 159.2 s, which is between the T ave maximum of the coal body and the roof rock. Figure 11 is the thermal imaging diagram of sample 11 # during the loading and unloading process. Six time points of 20 s, 40 s, 90 s, 130 s, 160 s, and 180 s are selected. It can be seen from the Figure 11 that the temperature change of the coal body during the loading and unloading process is more obvious. And the highest point of temperature is located in the middle axial area and the lower part of the coal body, and the temperature change range of the rock body is relatively small.
It can be seen from above that during the unloading and rupture process of the composite coal-rock, the T ave of the coal body appears earlier than the floor rock and the roof rock. In order to explore the precursory changes of infrared radiation temperature, only the T ave of the coal body can be compared with the stress peak time. Table 3 is a comparison of the peak arrival time and peak size of the average thermal infrared radiation temperature and stress on the coal surface under different unloading confining pressure rates for samples 1 #~1 5 # . The average thermal infrared radiation temperature peak arrival time of each sample is earlier than the stress peak arrival time. Combined with Figure 8, before the composite coal-rock is unstable and ruptured, the average thermal infrared radiation temperature of the roof rock, coal body, and floor rock exhibits a sudden increase and stepwise rise. Among them, the occurrence time of the precursors of the average thermal infrared radiation temperature of the coal body is ahead of the floor and roof rock, and the degree of change is the largest. Table 4 shows the fitting equations and related parameters of some samples at different confining pressure unloading rates. The fitting curves of samples 1 # , 6 # , and 11 # are shown in Figure 12. It can be seen that the T ave has a good linear correlation with the stress, and the correlation coefficients are all above 0.90. Compared with the fitting result in simulation, the correlation coefficient after experimental fitting is higher. This is the result of the large number of experimental sampling points and artificial points, which further proves that T ave has a good linear correlation with stress.
In summary, the experimental results show that the variation trend of the average thermal infrared radiation temperature on the surface of the sample is consistent with the simulation results. The average thermal infrared radiation temperature on the surface has a good correlation with the stress, which verifies the correctness of the simulation and analysis.

Conclusions
(1) Considering the combined effects of stress and thermal infrared radiation temperature on the unloading and rupture of anisotropic composite coal and rock, based on the basic theory of mechanics and the damage mechanics formula, it is deduced that the different components of the anisotropic composite coal and rock during the unloading and rupture process are at any arbitrary value. The relationship between 18 Geofluids thermal infrared radiation temperature and stress is averaged at all times, and a mathematical model of thermal-mechanical coupling under unloading and rupture conditions is established (2) The simulation results show that under unloading conditions, the thermal infrared radiation temperature of the composite coal is concentrated in the coal seam. After 142.5 s unloading the confining pressure, the overall temperature field reaches the strongest. The greater the rate of unloading confining pressure, the shorter the peak arrival time of T ave , and the larger the peak value. The T ave changes of roof rock, coal body, and floor rock are mainly divided into five stages, showing a trend of fluctuation, nearly linear rise, slow decline, and the sudden increase after unloading reaches the maximum and then slowly decreases The research results provide new ideas and new methods for monitoring and early warning of coal and rock dynamic disasters.

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

Conflicts of Interest
The authors declare no conflict of interest.