Analysis of Thermo-Hydro-Mechanical Coupling of Coal and Gas with Absorption/Desorption Thermal Effect

Coal adsorption/desorption gas can cause temperature changes. Moreover, the pressure difference between the high gas pressure area and the low gas pressure area in the coal may affect the coal adsorption/desorption gas thermal effect. To study the effect of adsorption/desorption thermal effect on the evolution law of coal and gas outburst precursor information under different pressure differences, the adsorption/desorption thermal effect experiments were carried out. And the thermo-hydro-mechanical coupling model was established, which reflects the interaction between coal and gas. Under the action of ground stress, pore pressure, adsorption expansion stress, and thermal expansion stress, the model includes the dynamic evolution of coal porosity and permeability, gas adsorption-desorption-diffusion-seepage, and energy accumulation and dissipation. Using the coupling model, the mechanism of adsorption/desorption thermal effect on the evolution of outburst precursor information and the mutual coupling relation between coal and gas is analyzed. *e experimental results show that the relation between pressure difference and coal sample temperature variation is linear, and the pressure difference has a quadratic parabolic relation with the temperature accumulation of coal sample. As the pressure difference increases, the temperature variation of the coal sample gradually increases, but the temperature accumulation gradually decreases. Numerical simulation results show that as the pressure difference increases, the changes in coal temperature and gas pressure gradually decrease, and the thermal expansion strain caused by the adsorption/desorption thermal effect decreases, and the permeability increases. Comparing the calculation results with or without the adsorption/desorption thermal effect, the volumetric strain, temperature, and permeability of the coal have significant changes when the adsorption/desorption thermal effect is considered, and the gas drainage volume increases as the pressure difference increases. *e effect of thermal expansion strain caused by adsorption/desorption thermal effect on gas seepage is greater than the effect of temperature change on gas seepage.


Introduction
Coal and gas outburst is a complex dynamic disaster in the process of coal mining. In recent years, most of the major safety accidents in China's coal mines are caused by the abnormal gas outflow and the coal and gas outburst. Coal and gas outburst has always been a serious threat to coal mine safety production and has become a major issue that needs to be solved urgently in coal mining [1][2][3][4][5]. As the mining depth increases, the in situ stress, temperature, and gas pressure of the coal seam increase and the permeability decreases. e interaction mechanism of coal and gas under the coupling action of stress field, seepage field, and temperature field is more complicated, and the temperature field has become an important research target [6][7][8][9]. Meanwhile, the pressure difference between the high gas pressure area in front of the working face and the low gas pressure area near the coal wall is an important parameter to control the outburst strength of coal and gas. Moreover, the adsorption/desorption thermal effect may be different under different pressure differences [10][11][12][13].
To study the characteristics of coal adsorption/desorption gas thermal effect, the influence of adsorption/desorption thermal effect on the interaction of coal and gas, and the evolution law and mutual coupling relation of coal and gas outburst precursor information under multifield coupling, scholars have carried out a lot of research work and established multifield coupling model suitable for a variety of complex situations [14][15][16][17]. Coal releases heat during the process of gas adsorption and absorbs heat during the process of gas desorption [18][19][20]. Jacek [21] gave a calculation method of the maximum volume of work which can be performed by a gas while being liberated from a coal bed. Rahman [22,23] reported the theoretical frameworks for the thermodynamic quantities for the adsorption of methane onto various carbonaceous materials. Temperature variation causes deformation, and deformation generates heat, which causes temperature variation. Liu [24] proposed the thermal deformation effect of adsorption-desorption on the basis of the molecular migration and energy conversion during the adsorption-desorption. At present, multifield coupling is a hot research target [25][26][27][28]; the fields of coalbed methane extraction, shale gas extraction, coal and gas outburst prevention, and nuclear waste treatment all involve thermohydro-mechanical coupling issues [29,30]. araroop [31] believed that existing coal bed methane (CBM) simulators usually treat the coal seam as a dual pore and single permeability system, ignoring the influence of water in the coal matrix. e model includes the influence of moisture in the coal matrix and the shrinkage and expansion of coal. Considering the competitive effects of thermal expansion, nonisothermal gas adsorption, and thermal fracturing, Teng [32] established a fully coupled thermo-hydro-mechanical model of coal deformation, gas flow, and heat transfer. Zhu [33] uses the finite element method to establish and solve a fully coupled model of coal deformation, gas transport, and heat transport. A general model was established to describe the evolution of coal porosity under the combined influence of gas pressure, thermally induced solid deformation, thermally induced gas adsorption changes, and solid deformation caused by gas desorption. Zhou [34] proposed a coupled thermal pore elastic model, which considers the compressibility and thermal expansion of the convective heat flow of the components, as well as the changes in the porosity and related characteristics of the saturated soil. e model also considers the thermodynamically coupled water and heat flow. Fall [35] gave a detailed formula for coupling moisture and gas transmission in deformable porous media is given.
is model takes into account the fluid flow under damage control and the coupling of hydraulic and mechanical processes. e model also considers the coupling of diffusion coefficient and mechanical deformation and considers the changes in capillary pressure caused by changes in permeability and porosity. To analyze the influence of gas desorption on gas outburst, An [36] established a mechanical process model in the process of gas migration and tunneling. e study understands the influence of absorption shrinkage and mechanical changes on the protrusion. Tao [37] carried out physical simulation experiments on coal and gas outbursts and established a thermo-hydro-mechanical coupling model of gassy coal.
In summary, the existing research did not use experiments to analyze the influence of coal adsorption/desorption gas thermal effect on the evolution law of outburst precursor information. In addition, there are few quantitative studies on the influence of adsorption/desorption thermal effect on the interaction between coal and gas. Using the self-designed coal and gas thermo-hydro-mechanical coupling experiment system and the cyclic-step adsorption/desorption experiment method, the temperature change of coal adsorption and desorption gas under different pressure differences was discussed. Based on the experimental results, a thermohydro-mechanical coupling model of gassy coal with the adsorption/desorption thermal effect was established, and the model is solved by COMSOL Multiphysics. Comparing the calculation results with or without the adsorption/desorption thermal effect, the effect of the adsorption/desorption thermal effect on the evolution law of precursor information and the effect of it on the mutual coupling relation between coal and gas were studied.

Experimental System.
To study the temperature evolution law of coal adsorption/desorption gas under different pressure differences, an experiment system of coal adsorption and desorption gas was designed. e device used in this experiment is to add a thermostat and a standard tank to the original device. e experimental device is shown in Figure 1.

Coal Samples and Experimental Methods.
e coal sample used in the experiment was taken from the Xinjiang Coal Mine of Yangquan Coal Group, which belongs to a high gas outburst mine. e 3# coal seam is the outburst coal seam. Obtain coal sample with the size of 40 cm × 40 cm × 60 cm from the tunneling head and send it to the laboratory after being sealed and packed underground. Use a cutting machine to cut the coal sample into a square specimen with a size of 5 cm × 5 cm × 10 cm, and then smooth it with a grinder. e basic parameters of the coal sample are shown in Table 1, and the prepared coal sample is shown in Figure 2. e author designed a coal cycle-step adsorption/desorption gas experiment. e experimental method is shown in Table 2. e initial charging pressure is 0.4 MPa, and the pressure difference is 0.2 MPa, 0.4 MPa, and 0.6 MPa. e initial pressure of the desorption experiment is the maximum adsorption equilibrium pressure of the adsorption experiment, and the ambient temperature is maintained at 30°C during the experiment.

Experimental Procedures.
e experimental process is shown in Figure 3. e initial gas pressure is the same. After each adsorption equilibrium, the pressure is increased by the same pressure difference until the target pressure is reached. And then the pressure is reduced by the same pressure difference to the minimum pressure.
Step 1. Check the air tightness of the device, install the sample, connect the temperature measurement system and the gas pressure measurement system, turn on the thermostat, adjust the temperature of the thermostat to 30°C, and start the data logger.
Step 2. Debug the triaxial stress loading system to adjust the axial pressure and confining pressure to predetermined values.
Step 3. Connect the high-pressure gas cylinder to the experimental system.
Step 4. Turn on the vacuum pump to vacuum degas the gas pressure loading system.
Step 5. Carry out coal adsorption gas experiment. Adjust the pressure reducing valve at the outlet of the high-pressure gas cylinder to the required pressure value so that the pressure of the standard tank reaches the predetermined value. After the gas pressure in the standard tank is stable, the gas in the standard tank is sent to the adsorption/desorption device for continuous adsorption for 12 hours. Record temperature and pressure data.
Step 6. Carry out coal desorption gas experiment. After the adsorption experiment is cycled to the target pressure, the gas pressure in the adsorption/desorption device is reduced with an equal pressure difference, and the desorption is continued for 12 hours. Record temperature and pressure data.
Step 7. Repeat steps 5 and 6 to adjust the gas pressure difference to 0.2 MPa, 0.4 MPa, and 0.6 MPa.

Temperature Variation during the Process of Gas
Adsorption. Figure 4 shows the comparison of adsorption experiment results under different pressure differences. Under the same pressure difference, the temperature variation of coal sample decreases during the step adsorption process. As the gas pressure difference increases, the temperature variation of coal sample caused by each pressurization gradually increases, and the temperature accumulation decreases. During the adsorption, the temperature change of the coal sample increases with the increase of the pressure difference.

Temperature Variation during the Process of Gas
Desorption. Figure 5 shows the comparison of desorption experiment results under different pressure differences. Under the same pressure difference, the temperature variation of coal sample decreases during the desorption. As the gas pressure difference increases, the temperature variation of coal sample caused by each pressure drop gradually increases, and the temperature accumulation decreases. During the desorption, the temperature change of the coal sample increases with the increase of the pressure difference.

Change Law of Temperature Accumulation of Coal
Sample. To further analyze the relation between the pressure difference and the temperature accumulation of coal sample, the curve in Figure 6 was fitted. e fitting equation is shown in Table 3. It can be seen from the figure that the temperature accumulation first increases with the increase of gas pressure and then tends to be flat, showing a trend of quadratic parabolic change. As the pressure difference increases, the temperature accumulation gradually decreases. When the pressure difference is constant, the temperature accumulation during the adsorption is slightly greater than that during the desorption.

Multifield Coupling Model and Governing Equations
Coal adsorption/desorption gas is an important part of the mutual coupling relation between coal and gas. e coal adsorption/desorption gas will be accompanied by the absorption and release of heat. erefore, when scholars establish the thermo-hydro-mechanical coupling model of coal and gas, the adsorption/desorption thermal effect has been studied as an important coupling relation [38][39][40][41][42]. Compared with previous research results, scholars have not    introduced the adsorption/desorption thermal effect into the coupled model from the experimental perspective, and there is no new research progress on the influence of adsorption/ desorption thermal effect on the interaction between coal and gas. Based on the experiments of coal adsorption/ desorption gas thermal effects with different pressure differences, this paper explores the influence of adsorption/ desorption thermal effects on the temporal and spatial evolution of outburst precursor information and the mutual coupling relation between coal and gas. Using the theory of continuum mechanics, the laws of static force balance, the mass conservation, and the energy conservation, considering the mutual influence of each physical fields, the thermohydro-mechanical coupling model of coal and gas is derived.

Model Hypothesis.
e thermo-hydro-mechanical coupling of coal and gas is an extremely complex physical process. e establishment of mathematical equations and the convenience of numerical calculations require certain simplified assumptions.
is article uses the reasonable assumptions and laws established by the predecessors to propose the following assumptions: e deformation of the coal is small, and the coal is in the linear elastic deformation stage, which obeys the generalized Hooke law. e effective stress change of the gassy coal skeleton follows the modified Terzaghi effective stress law. e seepage law of gas in coal seams conforms to Darcy's law.

Adsorption test Adsorption balance
Reduce pressure Step adsorption to target pressure and then desorption to minimum pressure

Curve
Fitting equation Shock and Vibration e adsorbed gas in the coal obeys the Langmuir adsorption equilibrium equation, and the free gas obeys the real gas state equation. e coal is full of single-phase gas. At the same point, the temperature of gas and coal are equal.
e inertial force and volume force in the process of gas seepage and coal deformation are negligible.

Dynamic Porosity and Permeability Models for Gassy Coal.
Compression deformation and adsorption expansion deformation caused by ground stress and gas pressure and thermal expansion deformation caused by temperature changes will cause the coal skeleton to deform. As the buried depth of the coal seam increases, the temperature, ground stress, and gas pressure will change significantly, and coal porosity is also in a process of dynamic change [43][44][45][46][47]. According to the basic assumptions, the porosity of coal can be defined as where V S is the coal skeleton volume, ΔV S is the change of coal skeleton volume, V P is the pore volume, ΔV P is the change of pore volume, ϕ 0 is the initial porosity, e is the volume strain, ΔV is the change of total coal volume, and V is the total coal volume. e volumetric strain increment of coal particles is mainly composed of three parts. ΔV S /V S0 is the strain increment caused by the pore gas pressure compressing coal particles. ΔV SF /V S0 is the strain increment caused by the expansion of coal particles during gas adsorption. ΔV ST /V S0 is the strain increment caused by thermoelastic expansion.
We obtain the porosity equation expressed in terms of effective stress and adsorption thermodynamic parameters: In addition, it is assumed that the coal permeability and porosity satisfy the following relation [48]: erefore, the permeability equation can be expressed as where k 0 is the initial permeability.

Deformation of Gassy
Coal. e continuous solid medium will produce stress and strain under the action of external force. According to the basic assumption, the stress-strain relation can be expressed as [49,50] where ε ij is components of strain tensor, σ ij is components of stress tensor, G is shear modulus, and K is bulk modulus. e temperature change will cause the expansion or contraction of the coal and produce thermal strain. Under the assumption of isotropy, the linear thermal expansion strain can be expressed as where ΔT is temperature change and ε T is the coefficient of volume thermal expansion of coal.
For coal with free gas, according to the porous linear elastic mechanics and the principle of effective stress, the linear compression strain caused by gas pressure is expressed as where K Y is volume compressibility and Δp is gas pressure change.
After coal adsorbs gas, the surface tension of coal particles decreases, and the coal will expand. On the contrary, after coal desorbs gas, the coal will shrink [51]. erefore, the linear adsorption expansion strain caused by coal particles adsorbing gas is expressed as where ρ is the density of coal, R is the universal gas constant, a is the limit adsorption capacity per unit mass of coal at reference pressure, b is the adsorption equilibrium constant of coal, and V m is the gas molar volume. Total strain includes strain caused by external force, pore compression strain caused by free gas, adsorption expansion strain, and thermal expansion strain under the coupling action of stress, gas pressure, and temperature, which can be expressed as According to the elasticity theory and the revised effective stress formula, the effective stress can be expressed by the balanced differential equation: Strain component and displacement component accord with Cauchy's equation:

Shock and Vibration
According to the strain-stress relation, an equation expressing stress with strain can be derived: Substituting equations (11) and (12) into equation (10), we obtain the improved form of Navier's equation expressed in displacement:

Gas Flow in Coal.
e ideal gas equation only applies to gases at low pressure and high temperature. e compressibility of real gas and ideal gas is different. erefore, real gas cannot simply be regarded as ideal gas. Gas can be regarded as an incompressible fluid under the conditions of low pressure difference, low flow rate, and small temperature difference. e state equation of gas in coal seam is usually expressed as where ρ g is the gas density, M is the gas molecular weight, T is the absolute temperature of the coal seam, M/R is the gas constant of specific gas, and Z is the compression factor. Assuming that the gas seepage in the coal seam conforms to Darcy's law, and the gravity effect is ignored, the gas seepage velocity can be expressed as where q g is Darcy's flow rate, k is permeability, ∇P is the gas pressure gradient, and μ is the dynamic viscosity coefficient of gas.
Coal is composed of solid framework and pores. e pores are filled with freely diffusing gas. Changes in gas pressure, temperature, or pores will cause gas desorption or adsorption and then cause gas seepage in the cracks. e mass conservation equation of gas seepage in coal can be expressed as where Q s is the source and sink and m is the gas content. Gas content includes free gas and adsorbed gas [52]. Assuming that free gas is an ideal gas, the gas content of coal can be expressed as where ρ c is the density of coal, M g is the molar mass of gas, p a is the gas pressure under standard conditions, and T a is the temperature under standard conditions. e first derivative of the gas content equation can be expressed as Substituting equations (14), (15), and (18) into equation (16), we have e derivative of porosity with respect to time is expressed as Substituting equation (20) into equation (19), we have 6 Shock and Vibration To simplify equation (21), we define some parameters as erefore, the equation (21) is simplified to Equation (23) is the governing equation of gas seepage in coal seams.

Heat Transfer Equation.
If the heat filtering effect is ignored, but the heat convection effect is considered, the constitutive equation of heat conduction can be expressed by Fourier's law as where q T is heat flow through fluid and solid media, λ M is thermal conductivity of coal, λ s is thermal conductivity of coal particles, λ g is thermal conductivity of gas, and C s is constant volume heat capacity of gas.
Assuming that the heat balance between the solid and the fluid is maintained, the heat balance equation on the characteristic unit body can be expressed as [53] z (ρC) M T zt where (ρC) M is the overall heat capacity of the porous medium, ρ s is density of coal skeleton, C s is constant volume heat capacity of coal matrix, K g is bulk modulus of gas, α g is the volumetric thermal expansion coefficient of gas under constant pore pressure and stress, and Q ht is the heat caused by coal adsorption/desorption gas. e first term on the left-hand side of equation (26) is the change rate of internal energy per unit volume caused by temperature change. e second term represents the heat dissipation due to the thermal expansion of the fluid. e third term represents the heat dissipation due to the thermal expansion of the coal matrix. e right-hand side represents the overall heat flow of fluid and solid medium [34].
Substituting equations (24), (25), and (27) into equation (26), according to the assumption of small deformation and the law of conservation of mass, it is further simplified as Equations (13), (23), and (28) constitute an equations' set describing the interaction between coal and gas, which considers the coupling among the stress field, seepage field, and temperature field. e solid mechanics module, PDE module, and solid heat transfer module of COMSOL Multiphysics can be used to solve the coupled model. Figure 7 shows the crosscoupling relation between coal deformation, gas seepage, and temperature changes. e stress changes the porosity and permeability of the coal and then affects the seepage process of free gas. e change of free gas pressure causes gas adsorption or desorption and causes the coal to expand or contract. e increase of temperature promotes the desorption of adsorbed gas and accelerates the seepage rate of free gas. e adsorption and desorption of gas is accompanied by the absorption and release of heat. Heat convection and internal energy dissipation in the process of gas seepage speed up the heat transfer. e temperature change produces thermal strain, which causes the coal to expand or contract. e deformation of coal is also accompanied by the storage or release of deformation energy [54,55]. e temperature change during coal and gas outburst includes two aspects: coal temperature rise in process of coal rupture and coal adsorption. Gas desorption and gas expansion lead to lower coal temperature. According to the assumption of linear elasticity, the temperature change during coal and gas outburst can be expressed as

Cross-Coupling Relations.
where ΔT 1 is temperature drop caused by gas expansion, ΔT 2 is temperature drop caused by gas desorption, and ΔT 3 is temperature rise caused by gas adsorption. According to the experimental results of coal adsorption/desorption gas thermal effect under different pressure differences, the energy change caused by coal adsorption/ desorption gas can be expressed as As the increase of temperature, the amount of adsorbed gas in the coal seam gradually decreases. According to the experiment results, the relation between adsorption constant and temperature is fitted.

Model Validation
To verify the rationality of the thermo-hydro-mechanical coupling model established in this paper, based on the physical simulation experiment of coal and gas outburst carried out by scholars, the numerical simulation study on the change law of gas pressure and temperature of coal in the aeration stage was carried out [37]. e geometric model is shown in Figure 8.  (1) (3)  Shock and Vibration e numerical calculation results and experiment results are shown in Figure 9. e experiment results showed that the gas pressure at point A remains unchanged for the first 8 seconds in the initial inflation stage, and then the gas pressure increases rapidly. After inflation to 110 s, the gas pressure reaches 1 MPa. e temperature at point B increases rapidly after the start of inflation. After inflating to 100 s, the temperature change rate decreases. e maximum temperature measured in the experiment is 295.5 K. Numerical calculation results under different pressure differences show that the gas pressure at point A remains unchanged for the first 10 seconds in the initial inflation stage, and then the gas pressure increases rapidly. After inflation to 110 s, the gas pressure reaches 1 MPa. e gas pressure changes under different pressure differences are the same, and the change variation is small. After the inflation starts, the temperature at point B increases rapidly. With the increase of the pressure difference, the maximum temperature is 297.3 K, 296.0 K, 295.5 K. e temperature change decreases as the pressure difference increases. e calculated result with a pressure difference of 0.6 MPa is close to the experimental result. e temperature change rate measured in the experiment after the start of inflation is greater than the numerical calculation result. According to the analysis, under the assumption of small linear elastic deformation, the energy released by coal deformation is small and the frictional heat of coal particles is ignored in numerical calculation. e difference in temperature changes under different pressure differences will be analyzed below.
Comparing the experimental results and the numerical calculation results, it is found that the numerical solution and the experimental solution are in good agreement. In summary, the thermo-hydro-mechanical coupling model of coal and gas can be used to simulate the evolution of precursor information during the preparation stage of coal and gas outburst.

Numerical Simulation and Discussion
Numerical calculations are carried out based on the thermohydro-mechanical coupling model of coal and gas, which considers the adsorption/desorption thermal effect. Comparing the calculation results with or without the adsorption/desorption thermal effect, the influence of the adsorption/desorption thermal effect on the evolution of coal and gas outburst precursor information is analyzed, and the mechanism of the adsorption/desorption thermal effect on the mutual coupling relation between coal and gas is studied.

Geometric Model.
e geometric model is shown in Figure 10. e working face surface is 5 m from the left boundary. e gas drainage borehole is 5 m away from the coal wall of the working face, the borehole diameter is 0.094 m, and the borehole depth is 20 m. e working face and the outside of the coal seam are free spaces, which can be regarded as exposed coal walls. Stress field boundary conditions: the upper boundary of the coal seam roof is subjected to a vertical uniform load of 9.8 MPa, the working face and the coal wall are free boundaries, the other boundaries around the model are supported by rollers, and the lower boundary of the coal seam floor is a fixed displacement. Boundary conditions of the seepage field: there is no gas flow between the coal seam and the rock, the initial gas pressure of the coal seam is 1.3 MPa, the gas pressure of the coal wall is 0.1 MPa, the suction pressure of the borehole is 14 kPa, and other boundaries have no flow. Boundary conditions of temperature field: there is no heat transfer between coal seam and rock layer, the initial temperature of coal seam is 292.4 K, the temperature of coal wall is 291.15 K, and other boundaries are zero flux. Along the advancing direction of the working face, a monitoring line is arranged in the coal seam. e parameters required for numerical calculation are shown in Table 4.

Evolution Law of Precursor Information without Adsorption/Desorption
ermal Effect. Figure 11 shows the distribution characteristics of the stress, gas pressure, and temperature of the coal in front of the working face and near the coal wall in the roadway without the adsorption/desorption thermal effect. e stress on the bottom of the working face is the smallest, and the coal walls on both sides and roof are stress concentration areas. As shown in Figure 12, the stress concentration area in the coal in front of the working face is 0-28.41 m, the gas discharge zone width is 10 m, the maximum volumetric strain in the stress concentration area is −0.00628, and the maximum coal temperature is 292.41 K. e coal temperature rises by 0.01 K, and the minimum permeability is 8.83 × 10 −19 m 2 . e increase in coal temperature without the adsorption/desorption thermal effect is caused by deformation energy. Under the assumption of small linear elastic deformation, the volumetric strain is small, and the coal deformation generates less heat. e location where the volumetric strain is Shock and Vibration the largest has the smallest permeability. Coal deformation is the result of the combined action of ground stress, adsorption expansion stress, and pore compression stress.

Evolution Law of Precursor Information with Adsorption/Desorption
ermal Effect. Figure 13 shows the distribution characteristics of the stress, gas pressure, and temperature of the coal in front of the working face and near the coal wall in the roadway considering the adsorption/ desorption thermal effect under different pressure difference. As shown in Figure 12, when the pressure difference is 0.2 MPa, the stress concentration area in the coal in front of the working face is 0-32.05 m, the width of the gas discharge zone is 13 m, the maximum volume strain is −0.00637, and the maximum coal temperature is 296.05 K, the coal temperature rises by 3.65 K, and the minimum coal permeability is 8.41 × 10 −19 m 2 ; when the pressure difference is 0.4 MPa, the stress concentration area in the coal in front of the working face is 0-29.47 m, the width of the gas discharge zone is 12 m, the maximum volume strain is −0.00632, and the maximum coal temperature is 294.39 K, the coal temperature rises by 1.99 K, and the minimum coal permeability is 8.60 × 10 −19 m 2 ; when the pressure difference is 0.6 MPa, the stress concentration area in the coal in front of the working face is 0-28.52 m, the width of the gas discharge zone is 11 m, the maximum volume strain is −0.00630, and the maximum coal temperature is 293.24 K, the coal temperature rises by 0.85 K, and the minimum coal permeability is 8.73 × 10 −19 m 2 . e gas pressure difference increased from 0.2 MPa to 0.4 MPa, the maximum difference in gas pressure was 10023 Pa, the increase in coal temperature decreased by 1.66 K, and the permeability increased by 1.94 × 10 −20 m 2 . e gas pressure difference increased from 0.4 MPa to 0.6 MPa, the maximum difference in gas pressure was 5002 Pa, the increase in coal temperature decreased by 1.14 K, and the permeability increased by 1.26 × 10 −20 m 2 . As the increase of the gas pressure difference, the width of the gas discharge zone decreases, and the temperature elevated amount and change rate of the coal temperature gradually  decrease. e deformation energy and the adsorption/desorption thermal effect cause the coal temperature to change. e higher the temperature, the greater the volumetric strain and the lower the permeability. e extreme point of volumetric strain and permeability overlaps with the position where the extreme point of temperature appears. Within the range of 0-10 m in front of the working face, near the extreme point, temperature has a more obvious effect on volumetric strain and permeability. Within 10 m from the working front, the variation range of volumetric strain is mainly controlled by ground stress, and the influence of temperature is relatively small. Figure 12 shows the influence of the adsorption/desorption thermal effect on the evolution of precursor information in the preparation stage of coal and gas outburst. Comparing the research results with or without the adsorption/desorption thermal effect, it is found that the adsorption/desorption thermal effect have obvious effects on coal gas pressure, temperature, volumetric strain, permeability, and other parameters. Comparing the evolution law of gas pressure, temperature, volumetric strain, and permeability of coal in front of the working face without the adsorption/desorption thermal effect, the temperature change caused by the adsorption/desorption thermal effect is larger, and the coal temperature increase is 3.64 K, 1.98 K, and 0.84 K; the adsorption/desorption thermal effect has little effect on coal gas pressure, and the maximum difference of gas pressure is 20285 Pa, 10596 Pa, and 5600 Pa; the adsorption/desorption thermal effect has a great influence on the volumetric strain of coal, and the increase in volumetric strain of coal is 9 × 10 −5 , 4 × 10 −5 , and 2 × 10 −5 ; the adsorption/desorption thermal effect has an obvious influence on permeability, and the decrease in permeability is 4.21 × 10 −20 m 2 , 2.27 × 10 −20 m 2 , and 1.01 × 10 −20 m 2 . e location where the gas pressure difference is the largest is about 1.5 m in front of the working face, where the slope of the gas pressure curve is the largest, and the temperature change is the largest. erefore, the effect of thermal expansion strain caused by adsorption/desorption thermal effect on coal permeability cannot be ignored. e thermal expansion strain caused by the temperature change in  the preparation stage of coal and gas outburst leads to the decrease of coal permeability and increases the outburst risk. Figure 14 shows the change law of gas pressure, permeability, and temperature in the coal and gas drainage volume around the borehole with or without the adsorption/ desorption thermal effect. Figure 14(b) shows that the temperature change is obvious within the influence range of the borehole. e temperature change near the borehole wall is large, and the farther away from the borehole, the smaller the temperature change. As shown in Figure 14(d), the maximum gas drainage volume is 964 m 3 without the adsorption/desorption thermal effect. As the pressure difference increases, the gas drainage volume increases gradually, and the maximum gas drainage volume of a single hole is 888 m 3 , 915 m 3 , and 940 m 3 . Comparing the calculation results without the adsorption/desorption thermal effect, the maximum gas drainage volume of a single hole is reduced by 76 m 3 , 49 m 3 , and 24 m 3 . Assuming that 1000 boreholes are arranged in the working face, the maximum gas drainage volume of a single hole is reduced by 7.6 × 10 4 m 3 , 4.9 × 10 4 m 3 , and 2.4 × 10 4 m 3 . With the increase of gas pressure difference, the volumetric strain of coal decreases and the permeability increases, which is beneficial to the desorption of gas. On the contrary, the coal temperature decreases with the increase of gas pressure difference, and the decrease of temperature is not conducive to the desorption and flow of gas. Analyzing the change law of gas drainage volume in boreholes, it can be seen that the effect of thermal expansion strain caused by adsorption/desorption thermal effect on gas seepage is greater than the effect of temperature change on gas seepage. In summary, the adsorption/desorption thermal effect not only causes the coal  temperature change, but also has a greater effect on coal seam gas seepage.

Conclusions
Based on the coal adsorption/desorption gas thermal effect experiment, this paper established a thermo-hydro-mechanical coupling model of coal and gas considering the adsorption/desorption thermal effect by analyzing the mutual coupling relation among the stress, gas, and temperature. e evolution of precursor information in the preparation stage of coal and gas outburst is simulated with or without the adsorption/desorption thermal effect, and the effect of the adsorption/desorption thermal effect on the mutual coupling relation between coal and gas is discussed. e main conclusions are as follows: (1) e relation between gas pressure difference and coal temperature variation is linear. Under the same pressure difference, the increasing range of the temperature variation of the coal sample gradually decreases with the increase of the gas pressure difference. e gas pressure difference has a quadratic parabolic relation with the temperature accumulation of coal sample. When the pressure difference is constant, the decreasing range of temperature accumulation decreases gradually with the increase of gas pressure difference. As the pressure difference increases, the temperature variation gradually increases, but the temperature accumulation gradually decreases.
(2) Based on the interaction between coal and gas, the mechanical model reflecting the dynamic evolution of porosity and permeability, gas adsorption-desorption-diffusion, and energy accumulation and dissipation was established with the adsorption/desorption thermal effect.
(3) As the gas pressure difference increases, the changes in coal temperature and gas pressure gradually decrease, the thermal expansion strain caused by the adsorption/desorption thermal effect decreases, and the permeability gradually increases. e evolution law of coal temperature depends on the rate and variation of gas pressure. (4) Considering the adsorption/desorption thermal effect, the volume strain, temperature, and permeability have significant changes. e gas drainage volume increases significantly with the pressure difference increases. e thermal expansion strain caused by adsorption/desorption thermal effect has a greater effect on gas seepage than temperature change.

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