Research Article  Open Access
Numerical Investigation of Complex Thermal CoalGas Interactions in CoalGas Migration
Abstract
Understanding the influence of temperature on the gas seepage of coal seams is helpful to achieve the efficient extraction of underground coal seam gas. Thermal coalgas interactions involve a series of complex interactions between gas and solid coal. Although the interactions between coal and gas have been studied thoroughly, few studies have considered the temperature evolution characteristics of coal seam gas extraction under the condition of variable temperature because of the complexity of the temperature effect on gas drainage. In this study, the fully coupled transient model combines the relationship of gas flow, heat transfer, coal mass deformation, and gas migration under variable temperature conditions and represents an important nonlinear response to gas migration caused by the change of effective stress. Then, this complex model is implemented into a finite element (FE) model and solved through the numerical method. Its reliability was verified by comparing with historical data. Finally, the effect of temperature on coal permeability and gas pressure is studied. The results reveal that the gas pressure in coal fracture is generally higher than that in the matrix blocks. The higher temperature of the coal seam induces the faster increase of the gas pressure. Temperature has a great effect on the gas seepage behavior in the coal seams.
1. Introduction
Coal seam gas is an important natural source and an important part in the clean energy structure in China [1]. However, most Chinese coalbed methane cannot meet the requirements of effective recovery because of the low permeability characteristic [2, 3]. Therefore, many technological means are adopted to enhance the coal permeability and develop the gas extraction efficiency [4, 5], such as, hydraulic fracturing technology [6]. Besides, hot injection has been adopted to enhance the extraction of coalbed methane [7, 8]. Therefore, the thermal evolution characteristic evolution during the extraction of coalbed methane needs to be studied.
The thermal stimulation of gas reservoir triggers the complex interaction between coal, gas, and temperature. It changes the deformation behavior of coal, as well as the heat and gas flow. Temperature has an obvious effect on the adsorption capacity of coal seam gas, which is negatively correlated with the increase of temperature [9–11]. At the same time, the thermal expansion coefficient is heterogeneous and anisotropic due to the different response of different components to temperature changes. Noorishad et al. [12] evaluated the coupled thermalhydraulicmechanical characteristics of rock and investigated the effects of thermal stresses on permeability through the deformation change of the fractures. Mctigue [13] developed a fluidsaturated, porous, thermoelastic model and calculated the heating of a half space at constant temperature conditions. Harpalani and Schraufnagel [14] found that desorption induced the shrinking of the coal matrix and increased permeability through seepage tests. Zhou et al. [15] consider thermodynamically coupled water and heat flow and demonstrate the influence of thermodynamic coupling through numerical and analytical solutions. Zhu et al. [16] analyzed the influence of temperature on coal permeability based on a complex thermalhydrologicalmechanical model and found that sorptioninduced volumetric strain affected the seepage behavior of coal seam significantly. Cai et al. [17] found that thermal treatment can promote methane desorption and increase the permeability through the numerical simulation. Qu et al. [18] developed a new thermal model with gas flow and matrix adsorption and suggested that the increase in temperature leads to the decrease in coal swelling and larger cleat aperture and higher coal permeability.
The interaction mechanism of heat injection has not yet fully understood for the design of engineering. It is indispensable to study the coalgas interaction under variable temperature. In this study, a coupled thermalhydrologicalmechanical model is established, considering the sorption characteristic and permeability evolution. Then, this complex model is implemented into a finite element (FE) and solved through the numerical method. Finally, the effect of temperature on coal permeability and gas pressure is studied based on this coupled model.
2. Equations of Coupled Model
2.1. Gas Seepage Characteristic in Coal Fractures
For the coal seam reservoir, the nonDarcy flow can be expressed aswhere is velocity vector, is density, , is porosity, and is the permeability of coal.
The above equation can be expressed in the following form:where is the Forchheimer coefficient.
The mass conservation equation of gas migration can be expressed aswhere is the velocity vector and is the gas source, and the mass content can be expressed as [19]where is the porosity, is the gas density, is the coal density, and is the absorbed content.
In the coal seams, the gas absorption volume can be expressed as
The Sorption Strain can be expressed aswhere is the content of absorbed gas and is the strain coefficient.
The gas density can be obtained as
The mass conservation equation can be expressed as
2.2. Gas Diffusion in Coal Matrix
The gas migration process experiences three substeps: flow in fracture, gas diffusion, and sorption in matrix [20]. Figure 1 shows a conceptual model for gas transport. The source term can be expressed aswhere is the exchange from matrix to fractures, is the gas concentration in matrix, is the gas concentration, is the shape factor, and is the sorption time.
(a)
(b)
(c)
(d)
The relationship of gas concentration and gas pressure is written aswhere is the molar mass.
The diffusion equation can be written aswhere is the equilibrium gas content.
Then, the diffusion equation can be expressed as
The permeability can be written as [21]where is the Klinkenberg coefficient, , where is the Klinkenberg effect coefficient, = 0.251.
2.3. Solid Mechanics Equation
The stress equation of coal seam can be written as
The stressstrain equation can be written as
The overall stress balance equation can be written as
2.4. Coal Permeability
The relationship between porosity and effective stress can be expressed as [22]
Then, the porosity is expressed as
By substituting the porosity, it can be rewritten aswhere and . is the initial pressure, and is the initial porosity.
The evolution characteristics of permeability can be expressed as
2.5. Energy Conservation
Total heat flux is expressed aswhere is the thermal flux and is the mass density.
The thermal balance can be written as [23]where and is the mass density.
The conservation of mass yields
Considering and , it yields [16]
These equations (8), (12), (16), and (24) describe the fully coupling model of coal seam gas migration, including mechanical deformation of coal (gas desorption induces coal shrinkage and selfheat induces coal expansion), gas diffusion from the matrix, gas flow, and heat transfer in fractures. The coupled relationship is illustrated in Figure 2. The THM model indicates the nonlinear response of gas migration in underground coal seams. It is hard to get an analytic solution of these equations, so they are achieved through the numerical method with Comsol Multiphysics (a powerful partial differential equation solver) and appropriate boundary conditions.
3. Model Establishment and Analysis
3.1. Verification with Field Test
To verify the validity of the model established in this paper in the calculation of simulated extraction, the numerical simulation is carried out according to the practical example of coal seam extraction. The actual size of the coal seam is 568 × 568 m^{2}, and the parameters are obtained from Mora and Wattenbarger [24]. The numerical results and the field production data are compared in Figure 3. Good agreement was achieved between the results of numerical simulation and the field production data, which proves the validity of the model.
3.2. Model Establishment
A model is established to analyze the effect of adsorption on the distribution of coal permeability and gas pressure. The size of the model is 0.1 m × 0.1 m, as shown in Figure 4. The upper boundary is free. And the condition of the upper boundary is the constant gas pressure 2.5 MPa and the constant temperature 333 K. The initial pressure of the coal seam is 0.1 MPa and the parameters in the calculation are listed in Table 1.

3.3. Effect of Adsorption on Gas Pressure
There are two different pressures between matrix blocks and coal fractures. Because of the high permeability in coal fractures, gas can quickly flow through the fractures. The gas pressure in fractures is generally higher than that in the matrix blocks in the calculation model. The matrix also continues to adsorb gas, and the pressure rises gradually. As the gas is injected from the top of the model, the crack pressure reaches 2.5 MPa firstly at the upper boundary, and the matrix permeability is also increasing. Finally, pore pressure and matrix pressure reach the equilibrium state. Figure 5 shows the gas pressure distribution at different times. As the time goes on, the difference of pressure becomes smaller and smaller. The initial pressure in the coal seam is 0.1 MPa. The gas pressure maximums of coal matrix are 0.3 MPa, 1.8 MPa, and 2.3 MPa at 1 h, 4 h, and 6 h, respectively. Finally, two pressures reach 2.5 MPa at about 24 hours.
Figure 6 shows the gas pressure distribution at the monitoring line. The difference between the matrix pressure and the fracture pressure can be seen more easily from the detection line. The gas pressure at the bottom of the model is low, while the gas pressure at the entrance is high. Within the equilibrium time of 1h, 4h, 8h, 12h, and 24h, the minimum pressure of the coal fracture is 0.1 MPa, 1.42 MPa, 2.16 MPa, 2.39 MPa, and 2.5 MPa, respectively; meanwhile, the minimum pressure of the coal matrix is 0.1 MPa, 0.64 MPa, 1.66 MPa, 2.19 MPa, and 2.5 MPa. With the passage of time, the pressure of internal coal mass increases gradually. The ratios of minimum matrix pressure and the minimum fracture pressure are 45%, 76%, and 92% at 4 hours, 8 hours, and 12 hours, respectively.
3.4. Effect of Adsorption on Permeability Distribution
Based on the dualporous medium model, the fracture pressure and the matrix pressure are calculated through different equations. Therefore, we can further understand the dualporous medium model by analyzing the evolution rule of permeability. Figure 7 is the variation diagram of coal permeability and gas content of coal mass at different time. The permeability of coal fractures decreases gradually. It is caused by the adsorption swelling of the matrix and the reduction of the crack opening in the coal mass.
Coal mass expends and the fracture aperture decreases when it absorbs gas. Meanwhile, the coal mass is limited by the boundary conditions, so coal deforms under the combined influence of the change of gas pressure, expansion stress, and effective stress [25]. Figure 7 shows that the permeability decreases with time, and the volumetric strain caused by adsorption is the main reason. There is an inverse relationship between the permeability and the content of gas; that is, the gas content of coal increases ceaselessly and permeability decreases gradually. The maximum of permeability ratio decreased from 0.901 to 0.876, while the minimum of gas content increased from 0.00239 m^{3}/kg to 0.0481 m^{3}/kg when the time increased from 1 h to 8 h. As time goes, the permeability of coal gradually stabilizes and tends to be at a uniform state.
The changes of permeability and porosity of a point are analyzed. The point is at the center of the model, and its coordinates are (0.05 m, 0.05 m). Figure 8 is the figure of permeability, porosity, gas content, and temperature change at the monitoring point A. Figure 8 shows that the evolution of permeability and gas content at point A are consistent with Figure 7. The temperature and gas content increase with the increase of time. The permeability and porosity of coal decrease rapidly before 150 h, which is mainly caused by the large pressure difference. Pore pressure and adsorption expansion affect together, which causes the change of permeability and porosity. With the decrease of pressure difference in coal mass, the change speed of permeability and porosity decreases gradually and tends to be stable finally.
3.5. Effect of Temperature on Permeability Distribution
Figure 9 gives the temperature change law of the model. The temperature at the entrance is relatively high, and the temperature in the coal seam keeps rising. The temperature at point A increases from 303.8 K at 1 h to 308.6 K at 8 h, finally reaching 322.7 K at 40 h. The deformation characteristics of coal seam are given in Figure 10. After the coalbed is heated up, the adsorption increases and the volume expands. The bottom, left, and right boundaries of coal seam are restricted, so there is no normal displacement in this boundaries. When the gas pressure rises on the upper boundary of the coal mass, the coal mass absorbs gas constantly, which leads to expansive deformation of coal mass and expanding outwards of the upper boundary. These characteristics also show that the volumetric strain induced by gas adsorption plays a key role. If the effect of gas adsorption is not considered, the effective stress will decrease, and the permeability of coal may increase. If the influence of gas adsorption is ignored, pore pressure may be underestimated.
The distributions of gas pressure and permeability under different temperatures are given in Figures 11 and 12. It indicates that temperature has a greater impact on the gas migration in the coal seam. The influence of temperature is reduced as the coal seam gas pressure approaches the setting pressure. High temperature can promote the desorption of gas and then change the fracture aperture and porosity of coal mass. The influence of temperature on the porosity of coal seam is a complicated process. The temperature evolution is obtained under the coupling effect of multiple factors, such as restriction on the boundary condition and adsorptive expansion deformation of coal mass. At the initial stage, the gas pressure and porosity of coal mass differ greatly at different temperatures. Gas pressure and permeability change constantly with the increase of time, and finally attain states of equilibrium. The gas pressure of point A is 1.05 MPa, 1.27 MPa, and 1.41 MPa at 1h when the initial boundary temperature is 333 K, 353 K, and 373 K, respectively.
In order to quantify the influence of gas diffusion behavior, the fracture pressure distribution without Klinkenberg effect and different adsorption times was also plotted as shown in Figure 13. The adsorption time significantly affects the distribution of fracture pressure. When the adsorption time is large, the gas pressure of coal seam will be underestimated. The influence of diffusion effect is the same as that of Klinkenberg effect at the same order of magnitude. In the sorption stage, the total volume expansion effect and Klinkenberg effect of coal increase the permeability. When the Klinkenberg effect is not considered, the adsorption pressure will be underestimated.
4. Conclusions
A new model about heat transfer was developed to study the dynamic problem of gas adsorption in coal seam under variable temperature conditions. In the model, the coal mass was regarded as a dualporous medium, and each point of the coal mass has two different pressures, i.e., pore pressure and matrix pressure. The relationship between adsorption equilibrium and permeability evolution was analyzed through the numerical method. The main results are as follows:
The gas pressure in fractures is generally higher than that in the matrix blocks in the calculation model. The matrix also continues to adsorb gas, and the pressure rises gradually. Finally, pore pressure and matrix pressure reach the equilibrium state.
The adsorption of coal increases with the increase of gas pressure, and the gas content increases continuously. The porosity of coal mass decreases with the increase of time, while the temperature and gas content increase with the increase of time.
The higher temperature induces the faster increase of the gas pressure. Temperature has a greater impact on the gas migration. As the coal seam gas pressure approaches the setting pressure, the influence of temperature is reduced.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the China Postdoctoral Science Foundation (Grant no. 2018M633549), Fundamental Research Funds for the Central Universities (Grant no. 2015XKZD06), and National Natural Science Foundation of China (Grant nos. 51674247 and 51804302).
References
 Y. Xue, Z. Cao, F. Du, and L. Zhu, “The influence of the backfilling roadway driving sequence on the rockburst risk of a coal pillar based on an energy density criterion,” Sustainability, vol. 10, no. 8, p. 2609, 2018. View at: Publisher Site  Google Scholar
 Y. Zhang, J. Underschultz, L. Langhi, D. Mallants, and J. Strand, “Numerical modelling of coal seam depressurization during coal seam gas production and its effect on the geomechanical stability of faults and coal beds,” International Journal of Coal Geology, vol. 195, pp. 1–13, 2018. View at: Publisher Site  Google Scholar
 G. Hu, N. Yang, G. Xu, and J. Xu, “Experimental investigation on variation of physical properties of coal samples subjected to microwave irradiation,” Journal of Applied Geophysics, vol. 150, 2018. View at: Publisher Site  Google Scholar
 Z. Cao, Y. Zhou, P. Xu, and J. Li, “Mechanical response analysis and safety assessment of shallowburied pipeline under the influence of mining,” CMES: Computer Modeling in Engineering and Sciences, vol. 101, no. 5, pp. 351–364, 2014. View at: Google Scholar
 Y. Xue, T. Teng, L. Zhu et al., “Evaluation of the nonDarcy effect of water inrush from karst collapse columns by means of a nonlinear flow model,” Water, vol. 10, no. 9, p. 1234, 2018. View at: Publisher Site  Google Scholar
 A. Zhou, K. Wang, L. Wu, and Y. Xiao, “Influence of gas ventilation pressure on the stability of airways airflow,” International Journal of Mining Science and Technology, vol. 28, no. 2, pp. 297–301, 2018. View at: Publisher Site  Google Scholar
 G. Hu, J. Xu, T. Ren, Y. Dong, W. Qin, and Z. Shan, “Field investigation of using water injection through inseam gas drainage boreholes to control coal dust from the longwall face during the influence of abutment pressure,” International Journal of Surface Mining Reclamation and Environment, vol. 30, no. 1, pp. 48–63, 2016. View at: Publisher Site  Google Scholar
 Y. Xue, F. Dang, Z. Cao et al., “Numerical analysis of heat and gas transfer characteristics during heat injection processes based on a thermohydromechanical model,” Energies, vol. 11, no. 7, p. 1722, 2018. View at: Publisher Site  Google Scholar
 D. Wang, R. Lv, J. Wei, P. Zhang, C. Yu, and B. Yao, “An experimental study of the anisotropic permeability rule of coal containing gas,” Journal of Natural Gas Science and Engineering, vol. 53, pp. 67–73, 2018. View at: Publisher Site  Google Scholar
 G. Wang, L. Dou, C. Wu et al., “Unstable energy triggering mechanism of rock burst,” Journal of China University of Mining and Technology, vol. 47, pp. 190–196, 2018. View at: Google Scholar
 D. K. Wang, S. M. Liu, J. P. Wei, H. L. Wang, and B. H. Yao, “A research study of the intrananopore methane flow law,” International Journal of Hydrogen Energy, vol. 42, no. 29, pp. 18607–18613, 2017. View at: Publisher Site  Google Scholar
 J. Noorishad, C. F. Tsang, and P. A. Witherspoon, “Coupled thermalhydraulicmechanical phenomena in saturated fractured porous rocks: numerical approach,” Journal of Geophysical Research: Solid Earth, vol. 89, no. 12, pp. 10365–10373, 1984. View at: Publisher Site  Google Scholar
 D. F. Mctigue, “Thermoelastic response of fluidsaturated porous rock,” Journal of Geophysical Research Solid Earth, vol. 91, no. 9, pp. 9533–9542, 1986. View at: Publisher Site  Google Scholar
 S. Harpalani and R. A. Schraufnagel, “Shrinkage of coal matrix with release of gas and its impact on permeability of coal,” Fuel, vol. 69, no. 5, pp. 551–556, 1990. View at: Publisher Site  Google Scholar
 Y. Zhou, R. K. N. D. Rajapakse, and J. Graham, “A coupled thermoporoelastic model with thermoosmosis and thermalfiltration,” International Journal of Solids and Structures, vol. 35, no. 34, pp. 4659–4683, 1998. View at: Publisher Site  Google Scholar
 W. C. Zhu, C. H. Wei, J. Liu, H. Y. Qu, and D. Elsworth, “A model of coal–gas interaction under variable temperatures,” International Journal of Coal Geology, vol. 86, no. 2–3, pp. 213–221, 2011. View at: Publisher Site  Google Scholar
 Y. Cai, D. Liu, and Z. Pan, “Partial coal pyrolysis and its implication to enhance coalbed methane recovery: a simulation study,” Energy and Fuels, vol. 31, no. 5, pp. 4895–4903, 2017. View at: Publisher Site  Google Scholar
 H. Qu, J. Liu, Z. Pan, Y. Peng, and F. Zhou, “Simulation of coal permeability under nonisothermal CO2 injection,” International Journal of Oil, Gas and Coal Technology, vol. 15, no. 2, pp. 190–215, 2017. View at: Publisher Site  Google Scholar
 Z. Junxiang, L. Bo, and S. Yuning, “Dynamic leakage mechanism of gas drainage borehole and engineering application,” International Journal of Mining Science and Technology, vol. 28, no. 3, pp. 505–512, 2018. View at: Publisher Site  Google Scholar
 D. J. Remner, T. Ertekin, W. Sung, and G. R. King, “A parametric study of the effects of coal seam properties on gas drainage efficiency,” SPE Reservoir Engineering, vol. 1, no. 6, pp. 633–646, 1986. View at: Publisher Site  Google Scholar
 L. J. Klinkenberg, The Permeability of Porous Media to Liquids and Gases, Drilling and Production Practice, American Petroleum Institute, Washington, DC, USA, 1941.
 H. Zhang, J. Liu, and D. Elsworth, “How sorptioninduced matrix deformation affects gas flow in coal seams: a new FE model,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 8, pp. 1226–1236, 2008. View at: Publisher Site  Google Scholar
 F. Tong, L. Jing, and R. W. Zimmerman, “A fully coupled thermohydromechanical model for simulating multiphase flow, deformation and heat transfer in buffer material and rock masses,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 2, pp. 205–217, 2010. View at: Publisher Site  Google Scholar
 C. A. Mora and R. A. Wattenbarger, “Comparison of computation methods for CBM performance,” Journal of Canadian Petroleum Technology, vol. 48, no. 4, pp. 42–48, 2009. View at: Publisher Site  Google Scholar
 Z. Li, X. He, and L. Dou, “Control measures and practice for rock burst induced by overburden fracture in topcoal caving mining,” Journal of China University of Mining and Technology, vol. 47, pp. 162–171, 2018. View at: Google Scholar
Copyright
Copyright © 2018 Xiaoyan Ni et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.