Gas (coalbed methane) extraction from low permeability coal seams has been always a difficult problem for prevention gas disaster in underground coal mine and utilization of coal mine gas in the world (Yuan 2015; Li and Fang 2014; Esterhuyse 2017). The variation of gas flow in low permeability coal seam is mastered under the impact of multiple factors such as gas, stress and temperature. This has a great significance to study the theory and technology for efficient coal seam gas extraction.

Many previous researches have proven that coal is a porous material with a large internal surface area, and it physically adsorbs gas. The adsorption and desorption processes for gas are basically reversible, and the higher the temperature is, the lower the adsorption quantity of gas (Zhang et al. 2013). At present, technical means were used to improve the gas extraction efficiency of low permeability coal seams, such as pressure relief extraction, hydraulic punching, hydraulic fracturing, etc. (Uth 2014). The principles of reducing the influence of ground stress and increasing the local permeability of coal seams are involved in improving the extraction efficiency in a short time. However, the attenuation of extraction rate is faster, especially in soft coal seams. Due to artificial fractures closing quickly under stress, the permeability decreases rapidly. Therefore, new technical ways should be found to improve the gas desorption and seepage velocity and delay the gas extraction rate attenuation in low permeability coal seams.

Drawing on the heat injection technology in the petroleum industry, the idea was putting forward the increasing of coal seam temperature by injecting heat to accelerate and increase gas production. The study found that raising the coal seam temperature can accelerate the desorption and releasing of adsorbed gas and increase the seepage rate of gas in coal seam fissures, which may become an effective method to promote gas extraction efficiency (Shahtalebi et al. 2016; Zhao et al. 2018; Teng et al. 2018). Previous research showed that some experimental and theoretical studies on the influence of temperature on coal seam gas adsorption and seepage have achieved many valuable results. Anderson et al. (2011) studied the influence of temperature on its a and b values using Langmuir adsorption theory and analyzed the changes in the a and b values below 200 °C. The coal desorption and permeability change under the coupling condition of temperature and stress fields was carried out by Li et al. (2009a, b) to establish the seepage control equation and obtain the gas desorption and seepage change mechanism with the increase in temperature. The influence of temperature change on the initial speed of methane emission was studied, and the result indicated that the relation between emission and temperature showed a proportional to quadratic function (Tailakov et al. 2015). The studies of the characteristics of methane desorption in coal under uniaxial stress and temperature concluded that temperature rise was a key reason for the large amount of desorption of adsorbed gas (Wang et al. 2013). Zhu et al. (2013) carried out research on the mechanical properties of coal samples for gas adsorption under the influence of temperature changes. It is believed that the adsorption capacity of coal samples was directly affected by the temperature changes, resulting in an increase in free gas content, thus reducing the strength of coal containing gas as a whole. Yin et al. (2013) and Zhang et al. (2018b) carried out permeability experiments to change the stress and temperature conditions of coal samples and discovered that permeability decreasing in the scope of 0–80 °C, while increasing sharply when the temperature exceed 80 °C. The COMSOL software was used to numerically simulate the solid–gas–heat coupling effect of gas-bearing coal & rock masses and obtained the gas seepage mechanism near the heading face but did not consider the temperature effect (Li et al. 2017). Ju et al. (2016) and Alam et al. (2015) tested the change mechanism of permeability under the change in temperature (below 80 °C) and confining pressure and believed that the increase in temperature will expand the coal matrix, resulting in pore closures and a decrease impermeability. Zhang et al. (2018a) introduced the theory of coal & gas two-state adsorption heat, constructed the control equation for the coal seam temperature field, improved the heat–fluid–solid coupling mathematical model of coal seam gas flow, and theoretically explained the interaction mechanism of adsorption and desorption, stress, temperature and seepage fields in the process of coal seam gas flow. Shan and Lai (2018) conducted seepage tests under different temperatures and pore pressures. It was believed that the increase in pore pressure can reduce permeability to a certain extent whereas the increase in temperature increases permeability. Wang et al. (2015) carried out experiments on the variation of coal sample permeability under different temperatures and axial pressures. It was concluded that the increase in axial strain will obviously reduce the permeability, and the reduce rate will gradually decrease with the increase in strain. The influence of temperature on permeability was related to the stress state inside the coal sample.

The above research results further prove that temperature has an important influence on gas adsorption and seepage flow. However, the coupling theory considering the heat–fluid–solid coupling and gas adsorption–desorption effects still needs improvement. In addition, there were few studies on the extraction of gas by injecting heat into coal seams close to the in-site conditions. Therefore, this paper aimed at the in-site conditions of coal heat injection and gas extraction, the heat–fluid–solid coupling model of gas-bearing coal was improved, and the variation state and extraction efficiency of coal seam gas were studied based on the actual conditions of low permeability coal seams after injecting heat by using numerical simulation.

Coal seams contain complex fissures and pores of different scales resulting in rich adsorption potential. A large number of gas molecules are adsorbed in the cracks and pores of coal. Coal seam temperature and stress interact with gas desorption and seepage coupling. The heat–fluid–solid coupling model of coal-containing gas needs to fully consider the changes in coal permeability and gas flow under the influence of multiple factors.

2.1 Fundamental assumptions

In order to make the three-field coupled seepage model of gas-bearing coal easy to calculate, the following assumptions were made:

1. (1)

   A coal seam is a uniform porous material with pores and fissures;

2. (2)

   Gas drainage has little influence on the overall deformation of a coal seam;

3. (3)

   The temperature change in the seepage process does not affect the gas dynamic viscosity;

4. (4)

   The free gas in the coal seam satisfies the ideal gas state;

5. (5)

   The coal matrix at any point in the coal seam is consistent with the gas temperature;

6. (6)

   The pore pressure is consistent with the gas pressure in the fissures;

7. (7)

   The influence of moisture on coal bed adsorption and seepage can be ignored.

2.2 Coal deformation control equation

The deformation field equation includes the solid balance, geometric and constitutive equations of a coal body under gas adsorption deformation. For the coal body that adsorbs gas, its variation control equation is:

1. (1)

   Solid equilibrium equation

According to the theory of elastic mechanics, the expression form of the equation was:

where σ_ij was the component of the stress tensor (N) and f_i was the component of volumetric force (N).

1. (2)

   Geometric equation

The strain and displacement components should satisfy the geometric equation:

where ε_ij was the strain tensor component and u_i was the displacement component (m).

1. (3)

   Constitutive equation of a coal body considering deformation of adsorbed gas

According to the elastic theory of porous media, the constitutive equation of coal considering gas adsorption deformation was (Wei 2012):

where G was the shear modulus(Pa); δ_ij was the Kronecker symbol (i, j = 1, 2, 3); Α was the Biot coefficient, α = 1 − k/k_s; k_s was the bulk modulus of the coal skeleton(Pa); K was the bulk modulus of coal(Pa); σ_kk = σ₁₁ + σ₂₂ + σ₃₃ was the normal stress component; α_s was the coefficient of thermal expansion(°C⁻¹); ΔT was the temperature increment; p was the gas pressure; ε_s was the coal body strain caused by gas adsorption, and the expression was:

where ε_L was the gas adsorption volume strain and b was the gas adsorption pressure constant.

Equations (1), (2), (3) and (4) were combined to obtain Navier-type equations, that was, a coal deformation control equation considering gas adsorption:

2.3 Gas seepage control equation

2.3.1 Coupling relationship among three fields

The literature established a three-field coupling relationship on the assumption that the coal matrix was a uniform block (Alam et al. 2015), and the specific process was as follows. The porosity in coal was closely related to its deformation, showing a positive or negative relationship. V_s was the volume of the coal skeleton (m³), and V_p was the coal body pore volume (m³), so the total coal body volume is V = V_s + V_p, and the porosity was:

Adsorption of gas makes the coal expansion strain consistent with the change in matrix block bulk strain. Homogeneous pressure s and gas pressure p borne by the coal body were used to express the constitutive equation of volume strain ΔV/V and pore volume strain ΔV_p/V_p of the coal matrix block:

Α was still the Biot coefficient of the coal matrix, α = 1 − K/K_s. β was the Biot coefficient of the coal body matrix pores, and β = 1 − K_p/K_s. K_p was the bulk modulus of pore gas (Pa). The average compressive stress \(\bar{\sigma }\) was

According to the Betti–Maxwell reciprocity theorem, it was calculated that:

where K_m represented the pore gas modulus.

Available from Eqs. (6)–(10):

Let

The above equation was simplified to:

The permeability of available coal seams was:

where k_m0 was the initial coal seam permeability.

2.3.2 Seepage control equation

The gas seepage mass balance equation in a coal body was:

where ρ_g was the density of gas in coal seam fissures, kg/m²; T was time, s; I was the gas source. m was the mass of gas in a unit volume of the coal seam, kg; q_g was the Darcy velocity vector, m/s. This equation ignored the gravity of gas, and its expression was:

where k was the gas permeability in a coal seam; μ was the viscosity coefficient of gas in a coal seam, N s/m².

The total mass of free gas in coal adsorption and fissures was m. The gas quality of the two parts was calculated separately, and the free gas quality in the crack was:

where m_f was the mass of free gas in the fractures, kg; and ϕ was the coal body porosity. According to the ideal gas state equation, the density of free gas could be expressed as:

where M_g was the gas molecule molar mass; R was Pussler’s gas constant of 8.31 J/(mol K), and T was the absolute temperature of the gas, °C. The mass of adsorbed gas was:

where m_a was the mass of gas adsorbed in the coal matrix, mg; ρ_ga was the density of gas under standard temperature and pressure, and ρ_c was the true density of the coal seam, mg. Thus, the gas quality of the whole coal seam was obtained by adding the gas quality of the two parts as follows:

where ρ was the apparent coal body density, mg/kg; A was the coal ash, %; and W was the coal seam moisture, %.

Since ϕ ≪ 1, the above formula could be simplified to:

Substituting formulas (15) and (21) into formula (14) could be obtained:

Formula (22) was sorted out, and the above parameters were substituted to obtain a fluid–solid coupling model considering the adsorption effect:

The above formula was the coal seam seepage equation of adsorbed gas containing the source and sink term I. The permeability k and porosity ϕ here varied with the strain, gas seepage and temperature fields, so the seepage field equation not only reflected gas seepage but also the coupling change in the coal seam stress and temperature fields. The solution could be obtained only when the stress field equation, seepage field equation and temperature field equation of coal body were obtained.

2.4 Temperature field control equation

The coal seam temperature field control equation could be derived from thermodynamic theory (Zhang et al. 2012; Sun et al. 2007). According to the first law of thermodynamics, the thermal conductivity of coal was as follows:

where ρ_s was the coal skeleton density; α_s was the coal coefficient of expansion and contraction affected by temperature; c_s was the thermal capacity of skeleton in the coal seam; J/(mol K); k_s was the thermal conductivity of the coal skeleton; Q_Ts was the strength of the heat source from the coal.

For the heat conduction of gas, on the basis of the energy conservation equation, the following results were obtained:

where ρ_g was the density of free gas, kg/m³. c_g was the heat capacity of free gas, J; k_g was the thermal conductivity of gas; ϕQ_Ts was the heat source intensity of the gas itself; h_g was the specific enthalpy of the gas; h_g = c_pΔT; c_p was generally the specific heat of gas at constant pressure, J.

When it was assumed that the thermally conductive solid and fluid were in thermal equilibrium, Eqs. (24) and (25) were added together to obtain the temperature field control equation of the coal seam:

Combined with multi-physical field coupling function of COMSOL Multiphysics, the coupling calculation of Eqs. (23) and (26) could be used to analyze the heat–fluid–solid coupling cases.

According to the above-mentioned heat–fluid–solid coupling model, a numerical model was established by using COMSOL Multiphysics to simulate the process of enhancing gas extraction by injecting heat into coal seam boreholes.

3.1 Numerical model and parameters

The size of the model was X × Y×Z = 50 m × 50 m × 10 m. The heat injection borehole was located in the center of the model (the coordinates of the hole on the y = 0 plane are x = 25, y = 0, z = 5). Drainage boreholes were arranged side by side on both sides of the heat injection borehole with a distance of 5 m. All boreholes had a radius of 0.1 m and a length of 30 m, to realize heat injection and gas extraction in coal seam. The periphery and bottom of the model were fixed displacement boundaries, and the upper boundary was a stress boundary. The buried depth of the coal seam was 500 m, and the upper boundary was loaded to simulate the overlying strata gravity. The surrounding, upper and lower boundaries of the model were seepage and thermal insulation boundaries. The initial coal seam pore pressure was 3 MPa, and the initial temperature was 30 °C. Considering that the model was symmetrical in X and Z directions, a monitoring line MN was set on the y = 0 plane to monitor the data, with the endpoint coordinates M(25,0,5) and N(50,0,5), as shown in Fig. 1. The property parameters of coal and gas in the numerical model were shown in Table 1. The parameters of coal were tested on the coal specimens sampled from Fangzhuang coal mine in Jiaozuo city, China.

Numerical model for a typical case of gas drainage in a coal seam

The heat injection temperature was various, and five numerical calculation schemes were designed, namely, no heat injection (keeping the original temperature at 30 °C) and heat injection temperature T_z at 60 °C, 90 °C, 120 °C and 150 °C. Other conditions of the model remained unchanged.

3.2 Initial conditions

Gas and stress initial conditions: when t = 0, the original coal seam gas pressure p₀(x, y, z) = 3 MPa. The stress field initial displacement u_i = 0(i = 1, 2, 3).

Boundary conditions: assuming that the roof and floor are impermeable strata, the boundary conditions of the seepage field were: p = p₀. Coal seam boundary gas flow q_s = 0. The drainage borehole boundary pressure p = P_a.

Stress boundary condition: the top boundary of the coal seam (z = 10 m) bears the gravity of the overlying strata of 10 MPa.

Displacement boundary conditions: lower boundary (z = 0) and surrounding boundaries (x = 0, x = 50, y = 0, y = 50) are displacement constraints.

Initial conditions of the temperature field: initial coal seam temperature was T₀ = 30 °C, and the heat injection hole temperature was T_z which according to the modelling case.

The unit of calculation time t was d, and the calculation time range was 0–1000 d for the above five schemes.

3.3 Numerical modelling results

3.3.1 Temperature and gas distribution in the coal seam after heat injection

Figures 2 and 3 reflected the different coal seam distribution gas pressure periods at 90 °C. The process of gas extraction in a coal seam was shown visually, and the gas pressure in the area closest to the pumping hole decreased rapidly. With the extraction process, the range of gas pressure decreased along the drainage hole, extending to the surrounding area and finally expanding to the whole model range.

Gas pressure distribution in a coal seam with heat injection temperature T_z = 90 °C

Gas pressure distribution on the monitoring line MN with heat injection temperature T_z = 90 °C

At 10 d after the extraction, the extracted range affected about 3 m on both sides of the bore hole but the pressure drop value was small, and the pressure was reduced only in the 0.5 m range of the drilling hole. At 100 d after the extraction, the extracted range affected about 10 m on both sides of the borehole. Due to the limited extraction time, the range of pressure reduction was still small. When extracting for 500 d, the pressure of the whole model range decreased and the drainage efficiency was obvious. The pressure dropped by 50% around 3 m range of the borehole. It was important to note that the pressure drop on both sides of the extraction bore was not symmetrical due to heat injection (Fig. 3). The pressure value on the left side was obviously lower than the right indicating that temperature affects gas desorption and seepage. High temperature on the left side means more gas desorption and more seepage into the drainage hole. The temperature on the right side was significantly lower than the left because it was farther away from the heat hole (Figs. 4, 5). Gas desorption and seepage flow were also lower than the left side. At 1000 d after the extraction, the gas pressure distribution form was similar to that of 500 d, but the overall pressure value decreased.

Temperature distribution in the coal seam with heat injection temperature T_z = 90 °C

Temperature distribution on the monitoring line MN with heat injection temperature T_z = 90 °C

Figures 4 and 5 showed the distribution of temperature values at different times in the coal seam when the heat injection temperature was 90 °C, showing the heat transfer process after injection in the coal seam borehole.

Heat was transferred outward after injection from the heat injection hole. When the heat was injected for 10 d, the rising temperature range was small, affecting only about 2.5 m around the heat injection hole. When the heat was injected for 100 d, the rising temperature range expands to 15 m, and the temperature of drainage hole that 5 m away from the heat injection hole increased to about 39.6 °C. When the heat was injected for 500 d, the temperature increased extended to the whole model range and the temperature of the drainage hole increased to 56.7 °C. At 1000 d of heat injection, the drainage hole temperature increased to 64.6 °C.

Figure 6 showed the influence of different heat injection temperatures T_z on the temperature distribution in the coal seam. The higher the heat injection temperature, the faster was the heat transfer. When heat was injected at 120 °C and 150 °C for 100 d, the temperatures at the extraction borehole were 44.4 °C and 49.2 °C, respectively. When heat was injected for 500 d, the temperatures at the extraction boreholes reached 70.1 °C and 83.5 °C, respectively.

Temperature distribution on the monitoring line MN with different injection temperatures after 100 d and 500 d injection times

Figure 7 showed the distribution of residual gas content in the coal seam after 100 d and 500 d of extraction at different heat injection temperatures. It could be seen that the higher the temperature and the longer the extraction time, the lower was the residual gas content in the coal seam.

Residual gas content distribution on the monitoring line MN with different injection temperatures after 100 d and 500 d injection times

When the drainage time reached 100 d, the influence range of the drainage borehole was little, limited to about 5 m, and the overall reduction of coal seam gas content was very low. Different heat injection temperatures had a certain influence on the reduction of coal seam gas content. The higher the injection temperature, the greater was the reduction of gas content. However, the influence of temperature was also limited to a small range near the heat injection hole and the drainage hole, and there was almost no difference in the influence of temperature on the gas content of the coal seam beyond the 5 m on the right side of the drainage hole.

At 500 d of extraction, the coal seam gas content decreased within the whole model range, and the largest decrease was near the extraction borehole, which was close to 50% of the original content. The higher the heat injection temperature, the lower was the residual gas content. According to the comparison on different temperature for gas extraction, the residual gas content decreased more significantly when temperature rise from 30 to 60 °C.

3.3.2 Effect of injection on the gas extraction rate and accumulative yield

Figures 8 and 9 showed the influence of heat injection temperatures on the gas extraction rate and cumulative gas extraction capacity. On the whole, the higher the injection temperature T_z, the higher was the rate of gas output. The drainage rate was largest at the beginning and then gradually decreased. The higher the temperature, the slower the drainage rate decreased. From 30 to 60 °C, the increase of extraction rate was most obvious. For 60 °C, 90 °C, and 120 °C cases, although the extraction rate was still proportional to temperature, the difference was not obvious, especially after the extraction of more than 500 d, the drainage rates were almost no change.

Output rate of gas versus drainage time

Cumulative output of gas versus drainage time

Heat injection had a significant influence on cumulative gas extraction. Within 0–1000 d of extraction time, the higher the temperature, the greater was the cumulative output and the increase in the output gap between different temperatures over time. Because the coal seam gas content was certain, it could be expected that the output would be consistent finally with time, but obviously raising the coal seam temperature greatly improved the gas extraction efficiency. At different temperatures, the output gap of between two cases in was greatest 30°C and 60 °C. For examples, when at 100 d and 500 d, the increasing percentage of gas were 55% and 47.8%, respectively. By contrast, the gap between 60 °C, 90 °C and 120 °C were smaller (Table 2).

3.3.3 Effect of heat injection temperature on gas extraction of coal seams with different permeability rates

Coal seam permeability is one of the main factors affecting the gas extraction efficiency and yield. The increase in the coal seam temperature will promote gas desorption, but the expansion deformation of the coal matrix caused by the temperature increasing will adversely affect permeability. Therefore, the initial coal seam permeability needs to be changed in the numerical model process.

According to the parameters in Table 1, the initial coal seam permeability k_f0 was 1.0 × 10⁻¹⁹ m². With this permeability as the benchmark to define the relative permeable rate Δk_f, that was, Δk_f = k_f/1.0 × 10⁻¹⁹ m², four cases Δk_f = 0.2, 1, 5, 10 were simulated. The drainage rate and cumulative drainage capacity were shown in Figs. 10 and 11.

Output rate and cumulative output of gas in coal seams with different permeability (coal seam temperature is 30 °C)

Output rate and cumulative output of gas in coal seams with different permeability (injection temperature is 90 °C)

When the coal seam temperature was 30 °C (that was, no heat being injected), the rate of extraction and the cumulative output were proportional to permeability. Because of 1.0 × 10⁻¹⁹ m² was a typical low permeability level for the coal seam, increasing the permeability would greatly increase the drainage efficiency. When the heat injection temperature was 90 °C, for all coal seams in this research with different permeability, the drainage rate increased. The corresponding drainage output were also significantly improved, but the enhancement of coal seam temperature on low permeability was more obvious considering the different permeability between the extraction rate and cumulative output for the relative difference between the reductions.

In order to analyze the influence of heat injection temperature and permeability on the extraction rate in different stages, based on the coal seam temperature of 30 °C (no heat injection), the relative drainage rate Δq was defined as the ratio of the drainage rate at different heat injection temperatures to the drainage rate at 30 °C at the same drainage time, and this could reflect the change in drainage rate caused by heat injection temperature.

Figure 12 illustrated the relative drainage rate of coal seams with different permeability in three typical time periods as a function of heat injection temperature. Through the simulation, it can be seen that, in a certain range, for coal seams with different permeability, the production rate can be improved by heat injection. This is due to the increase of temperature leading to the analysis of adsorbed gas and the increase of free coalbed methane concentration. At the same time, the higher the temperature, the smaller the adsorption deformation of the coal body, and the bigger the porosity and fracture opening in the coal body, and the more the permeability of the coal body, so the production rate. Obviously, the temperature increase had a more significant effect on the gas drainage rate in low permeability coal seams. Among the three extraction times selected, the relative drainage rate of the coal seam with the lowest permeability (Δk_f = 0.2) increased greatest with temperature, and that of the coal seam with the highest permeability (Δk_f = 10) increased smallest. Compared with different time periods, the slope of the relative drainage rate and temperature curve was biggest in the early extraction period (100 d), and the curve slope gradually decreased with time (500 d and 1000 d). The results showed that temperature increasing had the greatest effect on the drainage rate in the initial drainage time, and the effect decreased with time. According to the influence of different temperatures, in the range of 30–60 °C, the curve segment slope was significantly greater than other sections. This reflected that the same amount of energy was used; coal seams where the gas extraction raised from 30 to 60 °C had the highest efficiency. In the case of continuously increasing coal seam temperature (more than 60 °C), the increasing of extraction efficiency was reduced.

Relative drainage rate at three typical drainage times in coal seams with different injection temperatures and permeability

A variety of complex multi-field coupling problems were involved in the gas desorption seepage induced by heat injection in coal seams. The deformation of solids under stress, temperature, adsorption, desorption and the effect on permeability were considered and applied to the mathematical model. Although there was a certain difference with real coal seam conditions, the extraction rate and amount calculated numerically did not fully correspond to the real physical values. However, the simulated temperature, coal seam pressure and extraction rate may give some valuable results of in-site engineering.

The results showed that a temperature increase significantly affects the drainage rate but the effect was different in different temperature scope. From 30 to 60 °C, the ratio of drainage efficiency increasing was largest. The energy consumption for increasing drainage rate was smallest. It was no wonder that the result was consistent with previous research (Zhang 2011; Zhao 2012). When the coal temperature increased to 60 °C, the desorption rate increased most obviously. Although the continued increase in coal seam temperature could continuous to increase the drainage efficiency, with the excessive energy consumption and the impact of the mining environment, the difficulties of in-site implementation will be multiplied. Therefore, the coal seam was raised from the original temperature (20–30 °C) to 60 °C by injecting hot steam and other methods, which owned significant advantages, such as saving energy and easy operation while achieving a significant increase in the coal seam gas drainage rate.

The most common technique to extract gas for low permeability coal seams is bedding drilling, the simplest and most effective method to improve the drainage rate is reduce the spacing of drilling holes. Some coal mines’ drilling hole spacing, even for 1–2 m (Zhang 2011; Zhao 2012; Lu et al. 2015), increases the amount of drilling engineering and costs a lot of manpower and resources. However, the achievement was not as expected. Although coal seams are low thermal conductivity substances (thermal conductivity generally in the range of 0.2–0.5 W m⁻¹ K⁻¹), with steam injection heating, in about 3 months, the surrounding 5 m range of the coal seams can be heated to around 60 °C (Fig. 5). It doesn’t take long to heat the coal seam compared with the gas extraction time in coal mines. In addition, according to the simulation results, it is better to improve the drainage rate of low permeability coal seams by injecting hot heating steam. The low permeability coal seam drainage rate attenuation is fast, which is related to the difficulty of gas desorption from coal seam. The coal seams heating are required for a certain period of time, just coincident with the initial easy extraction time. The coal seam heating will slow down the drainage rate attenuation, increase the drainage rate and then enhance the double effect of the drainage output.

According to solid mechanics theory, the non-linear Darcy’s law of seepage and thermodynamics, the existing heat–fluid–solid mathematical model of coal containing gas was improved. The COMSOL Multiphysics software was used to establish a multi-field coupling model to perform the simulation. The process of promoting coal seam gas extraction by injecting heat into boreholes along the coal seam was simulated and analyzed. The following conclusions were obtained:

1. (1)

   The multi-physical field coupling model of heat transferring and adsorption in gas-bearing coal seams has been improved and solved by COMSOL Multiphysics. The desorption and seepage process of coal seam gas under the combined action of temperature, stress and adsorption affect can be well simulated to some extent.

2. (2)

   In the range of 30–150 °C, gases desorption and seepage flow will be promoted at higher temperatures. The gas drainage rate is proportional to temperature. During the extraction time of 0–1000 d, the increase in gas drainage rate was largest at the initial stage and reduces gradually with time.

3. (3)

   Different scope of temperature rise has different effects on increasing the gas extraction rate. From 30 to 60 °C, the ratio of increase is the largest. It will still increase with the temperature rise, but the increase ratio will reduce and there will be a remarkable increase in energy consumption, which will also worsen the in-site environment of underground coal mines.

4. (4)

   Comparing the gas extraction by injecting heat into coal seams with different permeability, the results show that the increase in temperature has the more obvious effect on the extraction efficiency of low permeability coal seams.

5. (5)

   Although the thermal conductivity of coal is low, injecting atmospheric pressure or pressurized steam into the coal seam will still raise the coal seam temperature within 5 m range of the hot injection borehole to about 60 °C in 3 months, which is much less than the gas extraction times in the coal mine. Therefore, coal seam heat injection was considered an efficient auxiliary means of gas extraction, with strong feasibility for low permeability coal seams in China.