Numerical Analysis of Heat and Gas Transfer Characteristics during Heat Injection Processes Based on a Thermo-Hydro-Mechanical Model

Heat injection is an important artificial technique, which can significantly enhance the extraction efficiency of coal seam gas (CSG) and reduce the outburst risk caused by CSG. Although heat injection has been comprehensively investigated, the effect of temperature on the coal–gas interactions in CSG extraction is still not clear. In this study, a thermo-hydro-mechanical model was developed considering the expansion of coal mass and the change of adsorption capacity induced by heat injection. Subsequently, the reliability of the model was verified through a comparison with other theoretical models and field data. Finally, a numerical simulation and parameter analysis of the heat injection process were performed and compared with the traditional gas extraction method. The simulation results show that heat injection can significantly increase the gas production rate and cumulative gas production through the gas desorption and the permeability increase. The gas content in the coal seam dramatically decreases in the vicinity of the production and heat injection wells under the condition of heat injection, which greatly accelerates the gas drainage. The coal deformation caused by thermal-induced gas desorption has a more dominant effect on the porosity than other factors, i.e., pore pressure, thermal strain and compressive strain.


Introduction
Coal seam gas (CSG) is recognized as a fuel and chemical raw material for clean energy recovery.Simultaneously, it is also a hazardous gas in coal mines that can cause explosions, worker suffocation, and gas outbursts.Gas extraction technology is adopted in coal mines to reduce the gas pressure and eliminate the risk caused by coal mine gas.Underground wells and ground drilling are usually used in the pre-extraction process, and underground drainage is the main system for realizing actual reservoir conditions [1][2][3][4][5][6].However, various coal seams in the southwest of China are low-permeability coal seams [7], and methane stored in these seams cannot be effectively extracted, which leads to difficulty in methane extraction during disaster control [8][9][10].
The permeability of a coal seam is related to the interactions between coal and gas, such as gas adsorption, effective stress and temperature change, and other factors caused by the pore/crack evolution of coal [11][12][13].Numerous models have been developed to study the evolution of permeability of coal during the gas extraction process, including Warren and Root [14], Robertson and Christiansen [15], Cui and Bustin [16], Shi et al. [17], etc., where the coupled effect of effective stress, reservoir pressure and adsorption deformation on permeability are considered.From the view of effective strain, Liu et al. [18] studied the isotropic changes in coal permeability and used an index, i.e., modulus reduction ratio, to represent the anisotropy in permeability and deformability based on the relationship between volumetric strain and permeability.The experimental results of gas adsorption show that temperature has an important effect on the methane adsorption capacity of coal [19].Therefore, scholars also focused on the influence of temperature on the gas seepage behavior in coal seam extraction process.For example, a multi-physical field model of porous media considering solid deformation, fluid flow, and heat transfer was discussed by Tong et al. [20].Based on the work of Liu et al. [18], Xia et al. [21] investigated the complex interaction between the gas transport and energy transport in porous coal media and proposed a coal spontaneous combustion model of the gas, air-mixed flow, and energy transfer.
Improving the quality of gas drainage is an important objective of scientific researchers in coal mining.Various technological methods are adopted to increase the permeability and achieve efficient gas extraction.For example, the hydraulic fracturing technology, carbon dioxide flooding, thermal stimulation, and other engineering technologies have been successfully applied to enhance the extraction of CBM [22][23][24][25][26][27].Thermal stimulation techniques have received tremendous attention in recent years because of their superiority [28][29][30][31][32][33].The principle of thermal stimulation is to heat the coal seam by injecting heated gas or heated water or to adopt the electromagnetic technology to heat the coal seam, and then to achieve efficient gas extraction from the production well [19,34].It can effectively increase the recovery rate by increasing the reservoir temperature.When the temperature rises, it results in gas desorption from the coal mass.
Some thermal stimulation methods have been proposed and applied in engineering practices.For example, cobalt ferrite nanoparticles could be used as magnetic feed antennas to enhance the magnetic field strength and enhance oil recovery [35].Moreover, hot water injection is also adopted to increase the gas extraction rate.Compared with traditional production, this heat treatment can significantly improve the recovery rate [36].High-temperature gas injection is also adopted to achieve enhanced CBM recovery.The effect of high-temperature gas injection is normally better than hot water injection because water may hinder the flow of the gas.For example, high-temperature nitrogen injection can accelerate the rate of methane drainage [37].Besides N 2 injection, CO 2 injection is also an effective method to develop the recovery rate [38].In addition, microwave and radio frequency energy can also be used to increase the shale gas recovery rate by increasing the shale temperature [39,40].Temperature causes complex interactions between coal and gas, including gas adsorption and pore structure change.On the aspect of experiments, Yin et al. [41] conducted triaxial compression tests of coal samples under different temperature conditions and examined the effect of temperature on the deformation and geo-mechanical characteristics of gas-containing coal.The experimental results showed that temperature was an important physical parameter affecting the adsorption characteristics of coal.Wang et al. [27] studied the effect of pre-injection of nitrogen on coal structure and permeability through experimental and numerical simulations.
However, fewer studies have been conducted on gas extraction with heat injection using a thermo-hydro-mechanical model [40,42].The effect of temperature on the coal-gas interactions in coal seam gas (CSG) extraction for heat injection is still not clear.In this study, a thermo-hydro-mechanical model was developed to study the coal-gas interactions during heat injection.The coal mass expands and its adsorption capacity decreases with the increase in temperature.These factors are embodied in the expression of permeability based on the concept of effective strain.In addition, non-Darcy effect and Klinkenberg effect of the fluid are also considered in the model.In previous thermo-hydro-mechanical models [40,42], non-Darcy effect and Klinkenberg effect of the fluid were rarely considered, which may cause underestimates in the analyses of heat and gas transfer characteristics during heat injection process.After validating the model with other models and field data, numerical simulations were Energies 2018, 11, 1722 3 of 20 conducted to investigate the effects of temperature on gas extraction.Finally, evolutions of the permeability, gas content, and gas production with and without heating injection were compared.The effect of heat injection on coal seam gas recovery was analyzed quantitatively.

Governing Equations of Model
The gas transport occurs in three stages: flow in the fractures, gas diffusion and sorption in the matrix.The governing equations involving these physical processes are given in this section.

Gas Flow in Fractures
A non-Darcy flow is usually adopted in CSG extraction, and it can be expressed as [43]: where → µ is the velocity (m/s), ρ g is the density (kg/m 3 ), υ is the dynamic viscosity (Pa•s), is the non-Darcy coefficient, β = 1.75 , φ f is the porosity, and k g is the permeability (m 2 ).
The gas velocity can be expressed as: where → µ is the Forchheimer number.For porous media, the flow equilibrium equation can be expressed as: where Q s is the gas source (kg•m −3 •s −1 ), and mass content m is calculated as [44]: where m is the gas mass content (kg/m 3 ), ρ ga is the gas density at standard temperature and pressure (STP) (kg/m 3 ), ρ c is the coal density (kg/m 3 ), and V sg is the content of the absorbed gas (m 3 •kg −1 ).The gas absorption volume can be expressed as [45,46]: where V L and p L are the Langmuir volume constant (m 3 •kg −1 ) and pressure constant (Pa) at the temperature T t , T ar is the absolute reference temperature (K), T t is the reference temperature (K), c 1 and c 2 are the pressure coefficient (Pa −1 ) and temperature coefficient (K −1 ).The sorption-induced volumetric shrinkage strain ε s is assumed as: where V sg is the content of the absorbed gas (m 3 •kg −1 ) and α sg is the sorption-induced strain coefficient (kg•m −3 ).The ideal gas law can be expressed as: where ρ g is the gas density (m 3 •kg −1 ), M g is the molecular weight of the gas (kg•mol −1 ), T is the gas temperature (K), and p a is the standard atmospheric pressure (MPa).Therefore, the gas flow equation can be rewritten as:

Gas Diffusion Process
In the gas diffusion process, the source term can be written as: where Q s is the exchange between the matrix and the fractures (kg•m −3 •s −1 ), D is the gas diffusion coefficient (m 2 /s), c m is the gas concentration in the matrix (kg/m 3 ), c f is the gas concentration in the fractures (kg/m 3 ), τ is the sorption time (s), and σ c is the shape factor of coal matrix (m −2 ).The gas concentration in the matrix and fractures can be expressed as: The following equation is introduced to calculate the exchange rate: where m e is the equilibrium gas content (kg/m 3 ) at pressure p f .Then the diffusion equation can be written as: where Mc is the molar mass of methane (kg/mol), and τ is the sorption time (s).Effective permeability k g can be written as [47]: where k ∞ is the intrinsic permeability (m 2 ) and b is a coefficient that increases with the reduction in the permeability according to [48] where α k is the Klinkenberg effect coefficient and α k = 0.251.

Mechanical Equilibrium Equation
For dual porosity media, the effective stress can be expressed as: where σ eij is the effective stress (MPa), σ ij is the total stress (MPa), δ ij is the Kronecker delta tensor, p f is the gas pressure in the fractures (MPa), and α is the effective stress coefficient for the coal fractures.
The strain-displacement relation for coal is expressed as: A Navier-type equation is expressed as: where G is the shear modulus (MPa), µ is the Poisson's ratio, and f i is the body stress (MPa).

Coal Permeability
The porosity model of coal seam is expressed as: Then the porosity is expressed as: Then the porosity can be rewritten as: where p 0 is the initial pressure (MPa), and φ 0 is the initial porosity.
The permeability can be expressed as: The apparent permeability in the fracture system is obtained as:

Energy Evolution
Total heat flux q T is given by: where q T is the thermal flux (J•s −1 •m −2 ), ρ s is the mass density of the gas (kg•m −3 ), C g is the gas-specific heat constant at constant volume (J•kg , and λ M , λ s , and λ g are the thermal conductivities of coal, the solid components, and the gas components The thermal balance is obtained as [49]: where (ρC) M is the specific heat capacity of the solid, (ρC , ρ s is the mass density of the gas (kg•m −3 ), and C g and C s are the gas heat constant and solid heat constant at constant volume, respectively (J•kg −1 •K −1 ), i.e., the specific heat capacity of gas and the specific heat capacity of coal.The conservation of the mass of the two phases can be written as: The fully coupled model can be used to analyze the heat injection process.The cross-coupling variable (porosity and permeability) connects different physical fields.The coupled model with different physics is shown in Figure 1.For the mechanical field, the change of pore pressure induces the coal deformation.The change of porosity and permeability will affect the flow velocity and further affect the distribution of pore pressure in seepage field.At the same time, the change of gas velocity and the coal volume also affect the heat transfer.It can be seen from Figure 1 that the effective strain is composed of 4 different parts.The variables (i.e., body strain, pore pressure, adsorption, temperature) contained in the effective strain exist in different physical fields at the same time, and cannot be solved by a single physical field.The cross-coupling variables connect different physical fields.

Analytical Validation
To verify the effectiveness of the aforementioned coupled model in the calculation of CSG extraction, a finite element model is applied, and the simulation results are verified with a simplified analytical solution.A one-dimensional (1D) linear steady gas flow model has a length of 20 m.The gas is injected into the rock with a constant rate of gas at the inlet, and the gas pressure is maintained constant at the outlet.Assuming that the porosity is constant, the gas will eventually reach a steady state.The equation for a 1D linear flow can be reduced from Equation (8) to:

Analytical Validation
To verify the effectiveness of the aforementioned coupled model in the calculation of CSG extraction, a finite element model is applied, and the simulation results are verified with a simplified analytical solution.A one-dimensional (1D) linear steady gas flow model has a length of 20 m.The gas is injected into the rock with a constant rate of gas at the inlet, and the gas pressure is maintained constant at the outlet.Assuming that the porosity is constant, the gas will eventually reach a steady state.The equation for a 1D linear flow can be reduced from Equation (8) to: The boundary conditions are: At the given boundary conditions, the 1D steady-state analytic solution was given by Wu et al. [48]: The calculation parameters are listed in Table 1, and these parameters are also adopted in the analytical solutions.The evolution of the gas pressure calculated by the numerical model and analytical model is compared in Figure 2a, which verifies the validity of the numerical model.

Numerical Validation
To verify the validity of the model established in this paper, the thermal-hydraulic-mechanical model proposed by Liu et al. [45] is analyzed, which describes the interaction between coal and gas under the effect of the gas pressure, thermal strain, and heat-induced gas desorption.Compared with the porosity model represented by Equation ( 22), the porosity model of Liu et al. [45] is expressed as: A 50 m × 100 m rectangular area was considered as the model area, with an initial pressure of 4 MPa.The boundary conditions are shown in Figure 2b. Figure 2b shows that the variation in the gas pressure and permeability in this study agree well with the results obtained using the model of Liu et al.; this further verifies the validity of the proposed model.

Verification by Field Data
In order to further verify the model in the simulation of coal seam gas extraction, the calculation model is established based on a simplified form of a gas production process at constant temperature condition in the United States.The CSG recovery area is 3.2 × 10 5 m 2 and the thickness of coal seam is 3.8 m.All the actual parameters used in the calculation are derived from Mora [50], (see in Appendix A). Figure 3 indicates the comparison between the numerical result and the field production data [50], which shows that the proposed model is applicable in simulating and predicting the CSG recovery.
A 50 m × 100 m rectangular area was considered as the model area, with an initial pressure of 4 MPa.The boundary conditions are shown in Figure 2b. Figure 2b shows that the variation in the gas pressure and permeability in this study agree well with the results obtained using the model of Liu et al.; this further verifies the validity of the proposed model.

Verification by Field Data
In order to further verify the model in the simulation of coal seam gas extraction, the calculation model is established based on a simplified form of a gas production process at constant temperature condition in the United States.The CSG recovery area is 3.2 × 10 5 m 2 and the thickness of coal seam is 3.8 m.All the actual parameters used in the calculation are derived from Mora [50], (see in Appendix A). Figure 3 indicates the comparison between the numerical result and the field production data [50], which shows that the proposed model is applicable in simulating and predicting the CSG recovery.

Model Establishment
To analyze the effect of heat injection on gas extraction, a calculation model (OPQR) was established, as shown in Figure 4.The heat is injected into the well through a hot water recycling system with watertight surface.This system only injects heat into the coal seam without introducing water.For coal deformation, four boundaries are restrained by the normal displacement.For the gas flow, zero fluxes are applied to these boundaries.For the heat transfer, the injection borehole is given a constant temperature of 362.4 K (89.4 • C).
The initial pressure is 3.0 MPa, and the parameters in the calculation are listed in Table 2.The boundary conditions are listed in Table 3.A monitoring line is selected in the diagonal line of the coal mass, and three monitoring points A (30 m, 30 m), B (55 m, 55 m), and C (80 m, 55 m) are used to study the change law of the coal permeability, gas pressure, and production rate.
displacement.For the gas flow, zero fluxes are applied to these boundaries.For the heat transfer, the injection borehole is given a constant temperature of 362.4 K (89.4 °C).
The initial pressure is 3.0 MPa, and the parameters in the calculation are listed in Table 2.The boundary conditions are listed in Table 3.A monitoring line is selected in the diagonal line of the coal mass, and three monitoring points A (30 m, 30 m), B (55 m, 55 m), and C (80 m, 55 m) are used to study the change law of the coal permeability, gas pressure, and production rate.

Thermal Evolution during Heat Injection
The temperature distribution at different production times is shown in Figure 5. Here, it was considered that 1 month was equal to 30 days.The temperature of the coal seam increases with

Thermal Evolution during Heat Injection
The temperature distribution at different production times is shown in Figure 5. Here, it was considered that 1 month was equal to 30 days.The temperature of the coal seam increases with time.Because the injection holes are not located in the center of the model, the temperature contours do not exhibit a central symmetry pattern.The pattern is symmetrical along the axis of the two wells.

Thermal Evolution during Heat Injection
The temperature distribution at different production times is shown in Figure 5. Here, it was considered that 1 month was equal to 30 days.The temperature of the coal seam increases with time.Because the injection holes are not located in the center of the model, the temperature contours do not exhibit a central symmetry pattern.The pattern is symmetrical along the axis of the two wells.The temperature of point C gradually increases to 300 K at a production time of 1.34 × 10 7 s (approximately 155 days).The increase in the temperature also accelerates the flow of the gas, which is beneficial for gas extraction.In the initial short period (10 5 s), the temperatures of the three points do not show a clear increase.Subsequently, because the distance between point B and the heat injection well is the closest, the temperature increases first at point B, then at C point, and finally at point A. Simultaneously, it can be observed that point B exhibits a fluctuation in the temperature increasing process, without showing a rising trend similar to points A and C. Finally, with the increase in time, the temperatures at the three points eventually attain the same value of 363 K, which is the temperature of the heat injection.
Figure 6 depicts the temperature variation of three specific points at different distances from the heat injected well.The initial seam temperature is 293 K.The temperature of point A gradually increases to 300 K at a production time of 9.38 × 10 5 s (approximately 11 days).The temperature of point B gradually increases to 300 K at a production time of 4.25 × 10 5 s (approximately 49 days).The temperature of point C gradually increases to 300 K at a production time of 1.34 × 10 7 s (approximately 155 days).The increase in the temperature also accelerates the flow of the gas, which is beneficial for gas extraction.In the initial short period (10 5 s), the temperatures of the three points do not show a clear increase.Subsequently, because the distance between point B and the heat injection well is the closest, the temperature increases first at point B, then at C point, and finally at point A. Simultaneously, it can be observed that point B exhibits a fluctuation in the temperature increasing process, without showing a rising trend similar to points A and C. Finally, with the increase in time, the temperatures at the three points eventually attain the same value of 363 K, which is the temperature of the heat injection.

Gas Pressure Distribution and Permeability Evolution during Heat Injection
The gas pressure evolution characteristic in the coal seam is presented in Figure 7.The gas pressure reduces with the increase in the extraction time.The gas pressure of coal seam is affected by the coupled influence of gas composition, temperature, and porosity.The pore pressure in the coal seam increases with the distance from the well.The velocity of the pressure differential near the borehole is high, and the gas pressure gradually increases with distance from the well.
The evolution law of permeability under the condition of heat injection can be seen from Figure 8.The initial permeability is 1.09 × 10 −18 m 2 before drilling.Owing to the pressure drop, the permeability around the borehole increases, showing a funnel-shaped trend.The permeability of the boreholes above 20 m is increased by approximately 30% (i.e., k/k0 = 1.3) in two months, compared with the initial condition of gas production.The desorption and drainage are caused by the decrease in the pore pressure in the coal seam.With the development of heat injection, thermal expansion and new adsorption equilibrium of coal occur.In addition, the permeability of coal is increasing.The injection heat can improve the permeability of the coal seam and is conducive for gas drainage.The permeability ratio at point B increased from 1.26 at a production time of 1 month to 2.26 at a production time of 20 months.

Gas Pressure Distribution and Permeability Evolution during Heat Injection
The gas pressure evolution characteristic in the coal seam is presented in Figure 7.The gas pressure reduces with the increase in the extraction time.The gas pressure of coal seam is affected by the coupled influence of gas composition, temperature, and porosity.The pore pressure in the coal seam increases with the distance from the well.The velocity of the pressure differential near the borehole is high, and the gas pressure gradually increases with distance from the well.The evolution law of permeability under the condition of heat injection can be seen from Figure 8.The initial permeability is 1.09 × 10 −18 m 2 before drilling.Owing to the pressure drop, the permeability around the borehole increases, showing a funnel-shaped trend.The permeability of the boreholes above 20 m is increased by approximately 30% (i.e., k/k 0 = 1.3) in two months, compared with the initial condition of gas production.The desorption and drainage are caused by the decrease in the pore pressure in the coal seam.With the development of heat injection, thermal expansion and new adsorption equilibrium of coal occur.In addition, the permeability of coal is increasing.The injection heat can improve the permeability of the coal seam and is conducive for gas drainage.The permeability ratio at point B increased from 1.26 at a production time of 1 month to 2.26 at a production time of 20 months.The porosity of coal seam is controlled by the volume of the coal matrix blocks, which consists of the effective volumetric strain, compressive strain, adsorption-induced strain, and thermal-induced strain.Figure 9 depicts the variation in the contributions of these components to porosity.A porosity exceeding 1 implies that the porosity is increasing under the effect of this factor.The contributions of the various factors to the process of heat injection are obvious.Once exposed to the atmosphere, the CBM will be discharged into the borehole by the differential pressure.Thus, the gas pressure in the coal seam is reduced.This pressure drop leads to a change in the effective strain and coal mass desorption, Energies 2018, 11, 1722 14 of 20 which in turn change the porosity [37].In the early stage of gas drainage, thermal-induced strain has little effect on the porosity.With time, the temperature of point A increases gradually.Although the temperature of coal typically increases, the deformation is restrained under the restriction of an external boundary, so that the compressive stress slightly reduces the porosity of coal.When the strain caused by adsorption is considered, the porosity of coal increases rapidly, which indicates that the effect of the adsorption-induced strain on porosity is the main influencing factor compared with the compaction and thermal expansion of coal.The contribution of sorption induced volumetric strain is 1.48 while the contribution of thermal strain is 0.81 at the time of 300 days.These comprehensive factors lead to a monotonous increase in the porosity.The permeability of coal seam is directly proportional to the porosity.The increase in the temperature caused by the injection of heat promotes the massive desorption of gas, and finally shrinks the matrix, which improves the permeability.
of the effective volumetric strain, compressive strain, adsorption-induced strain, and thermal-induced strain.Figure 9 depicts the variation in the contributions of these components to porosity.A porosity exceeding 1 implies that the porosity is increasing under the effect of this factor.The contributions of the various factors to the process of heat injection are obvious.Once exposed to the atmosphere, the CBM will be discharged into the borehole by the differential pressure.Thus, the gas pressure in the coal seam is reduced.This pressure drop leads to a change in the effective strain and coal mass desorption, which in turn change the porosity [37].In the early stage of gas drainage, thermal-induced strain has little effect on the porosity.With time, the temperature of point A increases gradually.Although the temperature of coal typically increases, the deformation is restrained under the restriction of an external boundary, so that the compressive stress slightly reduces the porosity of coal.When the strain caused by adsorption is considered, the porosity of coal increases rapidly, which indicates that the effect of the adsorption-induced strain on porosity is the main influencing factor compared with the compaction and thermal expansion of coal.The contribution of sorption induced volumetric strain is 1.48 while the contribution of thermal strain is 0.81 at the time of 300 days.These comprehensive factors lead to a monotonous increase in the porosity.The permeability of coal seam is directly proportional to the porosity.The increase in the temperature caused by the injection of heat promotes the massive desorption of gas, and finally shrinks the matrix, which improves the permeability.

Gas Content and Gas Production during Heat Injection
The evolution law of the gas content in coal seam is analyzed in this section.The gas content of coal seam can be obtained from Equation (4). Figure 10a displays a method to increase the CBM recovery by heat injection, whereas Figure 10b exhibits a traditional CBM recovery method.The gas content in coal changes more rapidly in the case of heat injection.In the absence of a heat source, the gas content in the coal seam changes very slowly.After heat injection, the gas content varies not only in the vicinity of the drainage borehole but also in the thermally affected zone.For example, the gas content of point C is 18.51 m 3 /t for 1 month.However, the gas content is 15.59 m 3 /t for 2 months, while the gas content is 21.14 m 3 /t without heat injection for 2 months.Overall, the gas

Gas Content and Gas Production during Heat Injection
The evolution law of the gas content in coal seam is analyzed in this section.The gas content of coal seam can be obtained from Equation (4). Figure 10a displays a method to increase the CBM recovery by heat injection, whereas Figure 10b exhibits a traditional CBM recovery method.The gas content in coal changes more rapidly in the case of heat injection.In the absence of a heat source, the gas content in the coal seam changes very slowly.After heat injection, the gas content varies not only in the vicinity of the drainage borehole but also in the thermally affected zone.For example, the gas content of point C is 18.51 m 3 /t for 1 month.However, the gas content is 15.59 m 3 /t for 2 months, while the gas content is 21.14 m 3 /t without heat injection for 2 months.Overall, the gas content is consistent of with the absorbed gas and free gas.Temperature usually has an obvious influence on the absorbed gas owing to the thermal desorption effect.content is consistent of with the absorbed gas and free gas.Temperature usually has an obvious influence on the absorbed gas owing to the thermal desorption effect.The evolution of the absorbed gas content in the coal seam is shown in Figure 11.The absorbed gas content is related to the air pressure and reservoir temperature and can be obtained from Langmuir adsorption (Equation ( 5)).The initial coal seam absorbed gas content is 16.6 m 3 /t.The change in the absorbed gas content caused by the gas drainage is very slow.Even after the gas is drained for six months, the absorbed gas content far from the borehole is high.For example, for the drainage time 1 month, 2 months and 6 months, the absorbed gas content of point B are 11.05 m 3 /t, 8.08 m 3 /t Energies 2018, 11, 1722 16 of 20 and 6.11 m 3 /t, respectively, under the heat injection condition.However, without heat injection, the absorbed gas content of point B are 16.57m 3 /t, 14.51 m 3 /t and 13.14 m 3 /t, respectively.The gas drainage in low-permeability coal seam is restrained.After heat injection, the absorbed gas content is reduced not only near the borehole but also in the heat affected zone, owing to the heat-promoted gas desorption.A high temperature implies a decrease in the absorbed gas content.After gas drainage for six months, the absorbed gas content decreases to 3.8 m 3 /t, which is much lower than that without heat injection (13.1 m 3 /t).
The gas production rates of different extraction methods can be seen from Figure 12.An obvious difference can be seen from the figure.The gas injection rate of the gas production is relatively high with the increase in time.Simultaneously, the gas production rate continues to decline, and finally reaches a stable state.In the last stationary stage, the gas production rate of injection heating fluctuates, whereas ordinary extraction maintains a steady rate.The gas production rate of the injection heat is 1.52 × 10 4 m 3 /d and 0.94 × 10 4 m 3 /d at 10 2 s and 10 6 s, respectively, and the gas production rate of common extraction is 1.31 × 10 4 m 3 /d and 0.81 × 10 4 m 3 /d, which is 85.0% and 86.9% of that in case of the injection heat.Figure 12 shows the cumulative gas productions obtained by different extraction methods.With the two types of gas drainage methods, in the initial stage, the two types of drainage quantities are relatively low before 10 5 s, whereas the cumulative gas production of the two drainage methods continues to increase.With the increase in time, the cumulative gas production exhibits a significant difference.The heat injection drainage quantity is 4.61 × 10 7 m 3 , whereas the ordinary drainage accumulated gas production is 1.78 × 10 7 m 3 at 10 9 s.The gas production rates of different extraction methods can be seen from Figure 12.An obvious difference can be seen from the figure.The gas injection rate of the gas production is relatively high with the increase in time.Simultaneously, the gas production rate continues to decline, and finally reaches a stable state.In the last stationary stage, the gas production rate of injection heating fluctuates, whereas ordinary extraction maintains a steady rate.The gas production rate of the injection heat is 1.52 × 10 4 m 3 /d and 0.94 × 10 4 m 3 /d at 10 2 s and 10 6 s, respectively, and the gas production rate of common extraction is 1.31 × 10 4 m 3 /d and 0.81 × 10 4 m 3 /d, which is 85.0% and 86.9% of that in case of the injection heat.Figure 12 shows the cumulative gas productions obtained by different extraction methods.With the two types of gas drainage methods, in the initial stage, the two types of drainage quantities are relatively low before 10 5 s, whereas the cumulative gas production of the two drainage methods continues to increase.With the increase in time, the cumulative gas production exhibits a significant difference.The heat injection drainage quantity is 4.61 × 10 7 m 3 , whereas the ordinary drainage accumulated gas production is 1.78 × 10 7 m 3 at 10 9 s.

Conclusions
(1) In this work, to analyze the evolution of porosity during heat injection, a fully coupled thermo-hydro-mechanical model was proposed.This model is suitable for coal compaction, thermal expansion, heat-induced gas desorption, and adsorption of coal deformation.The coal deformation caused by thermal-induced gas desorption has a more prominent effect on the porosity than other factors.
(2) The increase in the temperature also accelerates the flow of the gas, which is beneficial for gas extraction.The permeability of the boreholes above 20 m is increased by approximately 30% in two months, compared with the initial condition of gas production.As the injection heat continues, the average temperature of the entire coal seam increases.At high temperatures, owing to the gas desorption, the permeability of the heat affected zone can be significantly improved by heat injection.The maximum of permeability ratio increases to 3.34 at 20 months.
(3) The gas content in coal changes more rapidly with heat injection.After gas drainage for six months, the absorbed gas content decreases to 3.8 m 3 /t, which is much lower than that without heat injection (13.1 m 3 /t).In the absence of a heat source, the gas content in the coal seam changes very slowly while the gas content changes significantly in the thermally affected zone after heat injection.The gas content of point C is 15.59 m 3 /t for 2 months, while the gas content is 21.14 m 3 /t without heat injection.
(4) Heat injection can significantly increase the gas production rate and cumulative gas production.After heat injection, the gas composition varies not only in the vicinity of the production well but also in the thermally affected zone.The gas production rate of common extraction is 1.31 × 10 4 m 3 /d and 0.81 × 10 4 m 3 /d, which is only 85.0% and 86.9% of that in case of the injection heat.
In this paper, the heat and gas transfer characteristics during heat injection process is analyzed through numerical simulation.The next step is to conduct laboratory tests and field tests to further enhance the practicability of the model.

Figure 2 .
Figure 2. Model validation through comparisons with previous models.(a) Distribution of gas pressure calculated by two models; (b) Evolutions of the coal permeability and gas pressure calculated by two models.

Figure 2 .
Figure 2. Model validation through comparisons with previous models.(a) Distribution of gas pressure calculated by two models; (b) Evolutions of the coal permeability and gas pressure calculated by two models.

Figure 3 .
Figure 3.Comparison of the gas production rate between the modeling result and historic data.Figure 3. Comparison of the gas production rate between the modeling result and historic data.

Figure 3 .
Figure 3.Comparison of the gas production rate between the modeling result and historic data.Figure 3. Comparison of the gas production rate between the modeling result and historic data.

Figure 6
Figure 6 depicts the temperature variation of three specific points at different distances from the heat injected well.The initial seam temperature is 293 K.The temperature of point A gradually increases to 300 K at a production time of 9.38 × 10 5 s (approximately 11 days).The temperature of point B gradually increases to 300 K at a production time of 4.25 × 10 5 s (approximately 49 days).The temperature of point C gradually increases to 300 K at a production time of 1.34 × 10 7 s (approximately 155 days).The increase in the temperature also accelerates the flow of the gas, which is beneficial for gas extraction.In the initial short period (10 5 s), the temperatures of the three points do

Figure 6 .
Figure 6.Temperature evolution at three specific points.

Figure 6 .
Figure 6.Temperature evolution at three specific points.

Figure 7 .
Figure 7. Evolution of gas pressure along the link-line of two wells at different times.

Figure 7 .
Figure 7. Evolution of gas pressure along the link-line of two wells at different times.

Figure 7 .
Figure 7. Evolution of gas pressure along the link-line of two wells at different times.

Figure 9 .
Figure 9. Contributions of different coal-gas interactions to coal porosity at point A.

Figure 9 .
Figure 9. Contributions of different coal-gas interactions to coal porosity at point A.

Figure 10 .
Figure 10.Distribution of the gas content in the coal seam by two extraction methods.(a) heat injection; (b) gas extraction.Figure 10.Distribution of the gas content in the coal seam by two extraction methods.(a) heat injection; (b) gas extraction.

Figure 10 .
Figure 10.Distribution of the gas content in the coal seam by two extraction methods.(a) heat injection; (b) gas extraction.Figure 10.Distribution of the gas content in the coal seam by two extraction methods.(a) heat injection; (b) gas extraction.

Figure 11 .
Figure 11.Evolution of the absorbed gas content in the coal matrix after different production times.

Figure 11 .
Figure 11.Evolution of the absorbed gas content in the coal matrix after different production times.

3 )Figure 12 .
Figure 12.Gas production rate and cumulative production with two extraction methods.

Table 1 .
Parameters used for one-dimensional linear steady gas flow.

Table 3 .
Initial and boundary conditions in the reservoir simulation.

Table 3 .
Initial and boundary conditions in the reservoir simulation.