ThermalHydraulicMechanical (THM) Coupling Behaviour of Fractured Rock Masses
View this Special IssueResearch Article  Open Access
Xiao Fukun, Meng Xin, Li Lianchong, Liu Jianfeng, Liu Gang, Liu Zhijun, Xu Lei, "ThermosSolidGas Coupling Dynamic Model and Numerical Simulation of Coal Containing Gas", Geofluids, vol. 2020, Article ID 8837425, 9 pages, 2020. https://doi.org/10.1155/2020/8837425
ThermosSolidGas Coupling Dynamic Model and Numerical Simulation of Coal Containing Gas
Abstract
Based on gas seepage characteristics and the basic thermosolidgas coupling theory, the porosity model and the dynamic permeability model of coal body containing gas were derived. Based on the relationship between gas pressure, principal stress and temperature, and gas seepage, the thermosolidgas coupling dynamic model was established. Initial values and boundary conditions for the model were determined. Numerical simulations using this model were done to predict the gas flow behavior of a gassy coal sample. By using the thermosolidgas coupling model, the gas pressure, temperature, and principal stress influence, the change law of the pressure field, displacement field, stress field, temperature field, and permeability were numerically simulated. Research results show the following: (1) Gas pressure and displacement from the top to the end of the model gradually reduce, and stress from the top to the end gradually increases. The average permeability of the section of the model tends to decrease with the rise of the gas pressure, and the decrease amplitude slows down from the top of the model to the bottom. (2) When the principal stress and temperature are constant, the permeability decreases first and then flattens with the gas pressure. The permeability increases with the decrease of temperature while the gas pressure and principal stress remain unchanged.
1. Introduction
Gas flow in coal seams is an extremely complex thermosolidgas coupling dynamic problems, which involve multiple factor and multiple coupled fields. In the process of coal mining, various gas dynamic disasters occurred, such as coal and gas outburst and abnormal methane discharge, which are closely related to the interaction of coal and gas under the thermosolidgas coupling effect. The permeability of coal is a physical parameter that reflects the characters of methane seepage and significant parameters in studying thermosolidgas coupling. When the gas flows in the coal seam, the permeability change is very complex, affected by the gas pressure, temperature, and principal stress. Therefore, establishing the thermosolidgas coupling model of gassy coal, studying the permeability under the thermosolidgas coupling effect has an important significance to the prevention and control of gas dynamic disasters and improving gas drainage efficiency [1]. The thermosolidgas coupling model of gassy coal and the change of permeability have been studied before [2–11]. Xiao [2] made a study of gas release through floors from lower adjacent seam; the influence of coal mining process on the working face floor and gas release is explored, which provides a basis for the study of underground permeability. According to the characteristic of gas flow in lowpermeability coal, Guozhong et al. [3] found the dynamically coupled deformation and gas flow model by the dynamically varied porosity, gas seepage equation, and controlling equations of deformation of the coal. The results show that gas flow predictions correspond to the test results from the sample. Lower permeability gives closer agreement to the test results: The minimum error is 0.9%. As Guangzhi et al. [4] consider the variation of the porosity and permeability of coal containing gas at differential deformation stages, a dynamic model for porosity and permeability is developed based on the previous research results. Furthermore, taking coal containing gas as a kind of isotropic elastoplastic material and taking into account the effect of gas adsorption, the solidgas coupling model for coal containing gas is constructed. In addition, the numerical simulation model is built by using the finite element method, and the research results may have a meaning for further enriching and improving solidgas coupling theories for coal containing gas. As Ping et al. [5] consider the Klinkenberg effect and gas adsorption swelling effect on gas transportation in coal, the porosity model, the dynamic permeability model, and gassolid coupling model were established finally, and the numerical simulation results in natural discharging conditions is obtained. Qixiang [6] studied the effect of gas pressure on permeability in the whole stressstrain process of coal or rock. The result shows that the permeability increases with the increase of gas pressure within a certain range of gas pressure. Xinle and Yongli [7] considered the effect of temperature on permeability and made relevant studies. The experimental results show that the permeability changes with the decrease of effective stress, and the permeability increases with the decrease of effective stress. At the beginning, the permeability degradation gradient at low temperature decreases faster than that at high temperature, and at the end, the permeability gradient at high temperature rises faster than that at low temperature. Jiawei et al. [8] studied the mechanics and permeability characteristics of raw coal and briquette. The results show that the increase of gas pressure leads to the decrease of coalrock strength; therefore, it leads to the decrease of the permeability. Yongjiang et al. [9] carried out the permeability test that was performed in ZYS1 equipment for the measurement of permeability to study the effects of confining stress, axial stress, and temperature on permeability of coal briquette specimens. Ruiduan et al. [10] found that the temperature has a close relationship with permeability. When the confining pressure, pore pressure, etc. are fixed values, in the high temperature area, the permeability increases with the increase of temperature. Shugang et al. [11] conducted experiments on gas seepage characteristics of coal outburst under variable gas pressure under different coaxial pressure and confining pressure. With the increase of gas pressure, it is proved that the seepage velocity of gas increases as a quadratic polynomial function. The air permeability of the outburst coal is in the form of “” which first decreases and then increases. Zhao et al. [12–15] studied the mechanical properties, crack evolution, and fracture precursors of rocks. Mu et al. [16] used the numerical simulation method to study the shear deformation characteristics of grouting rock.
The aforementioned was mainly through experiments to research the effect of gas pressure, temperature, and principal stress on the permeability, and then processing and analyzing experiment data, getting the change law; there is no systematic theoretical analysis, through the analysis of theoretical model which is rarely reported. Under uniaxial and pseudotriaxial conditions, the monotonic factors are considered. The combined effects of stress, temperature, and pressure are not taken into account. Therefore, based on the theory of gas flow and coal deformation, the thermoelastic expansion, absorption gas expansion, pore pressure, and compression deformation of coal are fully considered, together with the dynamic equations of porosity and permeability which are comprehensively considered by gas pressure, temperature, and principal stress that are established, and the thermosolidgas coupling model were established finally. In addition, the numerical simulation model is built by using the finite element method, and the numerical solution of the model for a special loading case is given in terms of the constraint conditions and corresponding parameters. The research results may have a meaning for further enriching and improving thermosolidgas coupling theories for coal containing gas.
2. The Basic Assumption
At a macro level, coal is a porous medium. The factors affecting the flow of gas in coal are extremely complex; therefore, in the establishment of the thermosolidgas coupling model of gassy coal to grasp the main factors, secondary factors were abnegated, and the relationship between the main influencing factors was found. Therefore, the following assumptions are made: (1) gassy coal is homogeneous isotropic, and the gas flow in the coal is laminar flow, obeying Darcy’s law. (2) Coal is the singlephase gas saturation, the gas adsorption and desorption were completed, and the gas follows the Langmuir equation in an instant. (3) The coal gas content is free gas and adsorbed gas. (4) The deformation of coal belongs to the small deformation.
3. Theoretical Model of Porosity and Permeability of Coal
In the process of coal mining, the condition of coal is very complex, which is affected by gas pressure, stress, temperature, etc., causing the changes of internal structure of coal, which will lead to the change of porosity and permeability of coal.
3.1. Theoretical Model of Gassy Coal Porosity
According to the definition of porosity of porous medium [17, 18], the expression of the gassy coal porosity is
In the formula, is the initial porosity of gassy coal, is the volumetric strain of gassy coal, is the initial frame volume of gassy coal, and is the change of frame volume of gassy coal.
Allowing for the effect of gas pressure, principal stress, and temperature on coal deformation, the increment of volumetric strain is mainly composed of the increment of thermal strains , the increment of gas pressure strain , and the increment of adsorption gas expansion strain , namely
In the formula, the increment of adsorption gas expansion strain of unit volume of coal [19] is
Plug formula (2) and (3) into formula (1), we get a theoretical model of gassy coal porosity:
3.2. Theoretical Model of Gassy Coal Permeability
Permeability, an important index of coal seam permeability, can reflect the gas migration in coal. Therefore, study on dynamic change law of gassy coal permeability has a guiding significance to the prevention and control of gas dynamic disasters.
According to the KozenyCarman [19] equation, the expression between permeability and porosity is
In the formula, is the permeability, is the pore surface area of the coal volume per unit volume, and is the dimensionless constant.
In the process of coal stress and strain, it can be considered that the perunit volume of coal particle surface area keeps invariant, and combining with the dynamic equation of porosity (4), there is
The theoretical model of the gassy coal permeability is
The formula (4)–(7) define the porosity and permeability equation of gassy coal under the thermosolidgas coupling condition, which reflects the effects of gas pressure and temperature on permeability. The influence of gas pressure and temperature on permeability is very difficult by using the finite element software COMSOL.
4. The ThermoSolidGas Coupling Model of Gassy Coal
4.1. The Stress Equation
Under the small deformation and elastic isotropic assumption of gassy coal, the stress field should obey the equilibrium equation, thermosetting gas constitutive equations, and geometric equations. (1)Equilibrium of equationIn the formula, is a volume force.
According to the principle of effective stress of gassy coal (not considering damage) (2)Thermosetting gas constitutive equations
In the formula, is Lame coefficient; is shear modulus; is volume deformation; is thermal stress coefficient; is gas pressure coefficient; and is gas adsorption coefficient. (3)Geometric equationJoint type (9), (10), and (11) obtained stress equation of gassy coal is as follows:
4.2. Seepage Equation
The paper only considers the gas flow in the coal body in one way, and to meet the Darcy’s law [19] (1)Gas flow equation (2)The state equation of gas flow(3)Continuous flow equation
Joint type (14), (15), and (16) obtained gas seepage equation of gassy coal:
4.3. Temperature Field Equation
Coal gas adsorption, desorption, and seepage are thermal effect and are a nonisothermal process; containing gas coal body under external force deformation also produces thermal effects, so with coal and gas coupling temperature field equation of coupled stress field and seepage field effect. The temperature change of coal containing gas is mainly caused by gas adsorption and desorption, adsorption and desorption of gas equivalent to a heat source. Therefore, to solve the problem of coal containing gas temperature field is the threedimensional unsteady heat conduction with internal heat source. The temperature field equations of coal containing gas seepage are
4.4. ThermoSolidGas Coupling Model
Combining the type (4), (7), (12), (17), and (18), we get the control equations of thermosolidgas coupling model of gassy coal. This equation reflects the interaction of porosity and permeability of gassy coal, gas flow characteristics, and dynamic changes of coal and gas. For the control equations of thermosolidgas coupling model of gassy coal, to the specific problems, giving the initial conditions and boundary conditions constitutes the thermosolidgas coupling model of gassy coal.
According to the process of gas migration in coal body, consider the following boundary conditions. (1)Initial Conditions. The initial condition refers to the time of , the original distribution of gas pressure in gassy coal (2)Boundary Conditions. Displacement boundary conditions are as follows:
Stress boundary condition:
Pressure boundary condition:
Flow boundary:
5. Numerical Simulation Analysis
According to the above analysis, the mathematical model of thermosolidgas coupling mathematical model is a very complicated nonlinear equation group, which needs to be solved by numerical method. In this paper, using the finite element software COMSOL [12] by partial differential equations calculation tool PDE in the system platform for secondary development, we fill in the need of porosity, permeability, stress field, seepage field, and temperature field equation and establish the containing coal gas heatsolidgas coupling mathematical model, suitable to determine the actual boundary conditions and initial conditions. The distribution of seepage field and deformation field of gas bearing coal body and the permeability change of gas bearing coal body under the action of different gas pressure, effective stress, and temperature is obtained by using the calculation tool, and the relationship between seepage field and deformation field and gas pressure, effective stress, and temperature is also obtained.
5.1. Geometric Model and Boundary Conditions
The built model is the actual size of gassy coal in this paper, and the model is 50 mm long, 50 mm wide, and 100 mm high. The material parameters of model are shown in Table 1. According to the test of raw coal seepage, the finite element model was set up. The upper end of the model is subjected to vertical principal stress of , the left and right end under vertical principal stress of MPa, the front end under vertical principal stress of MPa, and the bottom under the displacement constraints; at the same time, the model has the gravity load. The gas pressure at the upper boundary of coal has initial gas pressure of MPa. Initial temperature of coal body is , and gas flow in coal.

5.2. Numerical Results
Figure 1 shows the displacement field of gassy coal under the gas pressure of 0.6 (MPa) and the temperature of 30 (°C). As can be seen from the figure, with the principal stress and gas pressure being applied, coal is gradually compressed along the axis, axis, and axis, and strain increment of thermoelastic expansion, absorption gas expansion, and porepressure compression increased cause the change of porosity and permeability.
Figure 2 shows the pressure field and velocity field, displacement field, principal stress field coupling of gassy coal under the gas pressure of 0.6 (MPa) and the temperature of 30 (°C). In Figures 2 and 3, longitudinal section shows that the change of gas pressure field and cross section represent the change of the velocity field, displacement field, and stress field, respectively. In Figure 2, red indicates that the gas pressure is relatively large, and blue shows a smaller value. As can be seen from the analysis of Figure 2, the pressure of gas is 0.6 MPa on the top ends of model. With the gas migration in coal, the gas pressure decreases gradually from top to bottom. Figure 3 shows that the gas seepage speed gradually decreases from the top to bottom of coal; the reason is that when a large amount of gas enter into the coal, the internal cracks and pores in the coal will begin to absorb, adsorbing gas will increase, and the inner wall of the coal will become thicker. The Klinkenberg effect is more obvious, namely, the gas in the coal is easier to slide, resulting in the decrease of gas seepage speed. In Figure 1, with the principal stress and pressure of gas being applied, the top of coal will become the first and the maximum to withstand stress; therefore, the maximum deformation of compressed top coal, maximum relative compressive deformation of top coal, and the displacement are maximal. With the transfer of the stress, the top to bottom of coal gradually affected by stress; thus, the displacement decreases gradually from top to bottom. Figure 2 shows the change of stress of gassy coal, the color represents the stress value, the blue color represents the area of small stress value, and red represents the area of large stress value. The deeper the blue, the smaller the stress value; the deeper the red, the larger the stress value. As seen in the figure, the stress increases gradually from top to bottom under the action of gas pressure and principal stress.
Figure 4 shows the average permeability of the section of the thermosolidgas coupling model of the gassy coal with the change of gas pressure. As can be seen from the figure, the average permeability of he section of the model tends to decrease with the rise of the gas pressure, and the decrease amplitude slows down from the top of the model to the bottom. As far as the reason is concerned, with the increase of gas pressure, the adsorption capacity of coal to gas is enhanced, thus reducing the porosity of coal. Porosity affects the adsorption, analysis, and permeability of coal and thus reduces the permeability of coal. With the increase of gas content, coal gas has adsorption expansion deformation, which will not produce new fine holes, but expands inward, narrowing the internal cracks, thus reducing the permeability. With the increase of gas pressure, the suction capacity of coal gradually tends to balance, so the permeability decline rate gradually decreases and tends to be flat.
To study the variation of permeability in the thermosolidgas coupling, the dynamic model of gassy coal, the numerical value of different gas pressure, different principal stress, and different temperatures was calculated at the outlet end of the model, the numerical solution is introduced into MATLAB, the permeability curve of gassy coal is obtained by programming calculation, and the variation law of permeability is obtained. In Figure 5, when the temperature and the main stress are kept constant, the permeability of the gas pressure increases with the increase of the gas pressure at 0.30.9 MPa; and when the gas pressure is increased to 0.9 MPa, the trend of the decline of permeability is gentle. This is an increase of gas pressure, inevitably accompanied by an increase in gas flow, when a large number of gas penetration into the inside of the coal body, the existing cracks and pores in the coal body, will begin adsorption. When a large amount of gas is adsorbed, under the restraint of the three principal stresses, the coal body will not expand outward but will develop internally. Due to internal expansion effect, the original fractures and pores in the coal are compressed, or the gas molecules occupy; originally, supplied gas flow channel was diverted, resulting in gas more difficult from the coal permeability, and hence, the permeability decreases. The original offer of gas seepage channel was diverted, causing gas harder to penetrate from the coal body, so the permeability decreased. On the other hand, as the gas pressure and flow rate increase, the inner wall of the coal becomes thicker. The Klinkenberg effect is more obvious, namely, the gas in the coal is easier to slide; this leads to a decrease in permeability. With the increase of gas pressure, the suction capacity of coal gradually tends to be balanced. Swelling deformation of coal framework is also becoming smaller and smaller; therefore, the rate of permeability decline gradually decreases and tends to be flat. When the gas pressure and principal stress remain constant, the permeability decreases with the increase of temperature. This is because the viscosity coefficient of the gas increases with the increase of temperature, making it more difficult for the gas to flow in the coke. Expansion occurs inside the coal, and the pores inside the coal are compressed, which narrows the gas passage and thus reduces the permeability of the coal. When the gas temperature and pressure are fixed, the coal permeability reduces rapidly firstly and slowly where after with the principal stress increasing. With the increase of the principal stress, the deformation of coal body is gradually increasing, the fracture pore is compacted, the seepage channel becomes smaller, and so the permeability decreases. When the principal stress is higher, the coal body is gradually compressed, the compression effect is reduced, and the seepage channel is stable, so the permeability change tends to be gentle.
6. Conclusion
(1)Based on gas seepage characteristics, effects of gas pressure, principal stress, and temperature on gas coal are considered; the porosity model and the dynamic permeability model of coal body containing gas were derived. Furthermore, taking into account the basic theory of thermosolidgas coupling, the thermosolidgas coupling dynamic model was established(2)By using the thermosolidgas coupling model, the gas pressure, temperature, and the principal stress influence, the change law of the pressure field, displacement field, and other physical field and permeability were numerically simulated. Research results show that gas pressure and displacement from the top to the end of the model gradually reduce, and stress from the top to the end gradually increases. Average permeability of section of the model tend to decrease with the rise of the gas pressure, and the decrease amplitude slows down from the top of the model to the bottom(3)Applying the COMSOL Multiphysics finite element software to numerically solve the thermosolidgas coupling mode of gassy coal permeability, the effects of gas pressure, temperature, and the principal stress on relative permeability of gassy coal are analyzed. The results show that when temperature and principal stress are fixed, the permeability reduces rapidly firstly and slowly whereafter with gas pressure increasing. When gas pressure and principal stress are kept constant, the permeability decreases with the increase of temperature. If gas temperature and pressure are fixed, the relative permeability of coal reduces rapidly firstly and slowly whereafter with the principal stress increasing; the results have an important significance of the prevention and control of gas dynamic disasters and improving gas drainage efficiency
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflict of interest.
Acknowledgments
The research was supported by the National Natural Science Foundation of China (51774121), Joint Guidance Project of Natural Science Foundation of Heilongjiang Province (LH2019E087), the National Nature Science of China Foundation (51674107), and Innovative Scientific Research Project for Postgraduates of Heilongjiang University of Science and Technology (YJSCX2020103HKD).
References
 Z. Zhigang, T. Yunliang, and C. Guoqiang, “Catastrophe mechanism analysis of coal and gas outbursting coal heading face,” Rock and Soil Mechanics, vol. 29, no. 6, pp. 1644–1648, 2008. View at: Google Scholar
 X. Fukun, D. Liqun, and G. Zhihui, “Laws of floor breaking in coal mining face and gas extraction application,” Journal of China Coal Society, vol. 35, no. 3, pp. 417–419, 2010. View at: Google Scholar
 G. Hu, J. Xu, H. Wang, Z. Jin, and R. Wu, “Research on a dynamically coupled deformation and gas flow model applied to lowpermeability coal,” Journal of China University of Mining and Technology, vol. 40, no. 1, pp. 1–6, 2011. View at: Google Scholar
 G. Yin, D. Wang, D. Zhang, and G. Huang, “Solidgas coupling dynamic model and numerical simulation of coal containing gas,” Chinese Journal of Geotechnical Engineering, vol. 30, no. 10, pp. 1430–1436, 2008. View at: Google Scholar
 P. Guo, S. Cao, Z. Zhang, Y. Li, Y. Liu, and Y. Li, “Analysis of solidgas coupling model and simulation of coal containing gas,” Journal of China Coal Society, vol. 37, no. 2, pp. 330–335, 2012. View at: Google Scholar
 H. Qixiang, “Gas pressure of coal or rock materials should all be stress  strain process gas permeability characteristics,” Materials Review, vol. 24, no. 8, pp. 80–83, 2010. View at: Google Scholar
 Y. Xinle and Z. Yongli, “Experimental study on the effect of temperature on gas permeability of coal under the influence of gas and solid coupling,” Journal of Geomechanics, vol. 14, no. 4, pp. 374–379, 2008. View at: Google Scholar
 J. Li, J. Liu, Z. Zhang, J. Pei, X. Liu, and R. Zhang, “Investigation on mechanical properties and permeability of coal under gas pressure,” Journal of China University of Mining & Technology, vol. 42, no. 6, pp. 954–960, 2013. View at: Google Scholar
 Y. Yu, Z. Hua, and Z. Chunhui, “Effects of temperature and stress on the permeability of standard coal briquette specimen,” Journal of China Coal Society, vol. 38, no. 6, pp. 936–941, 2013. View at: Google Scholar
 C. Ruiduan, C. Haiyan, and X. Xuefu, “Experimental study on the influence of temperature on the permeability of coal samples,” Coal Engineer, vol. 1, pp. 13–16, 1998. View at: Google Scholar
 Y. Dai, W. Chen, G. Wu et al., “Study on elastoplastic damage model of unsaturated rock mass and its application,” Chinese Journal of Rock Mechanics and Engineering, vol. 27, no. 4, pp. 728–735, 2008. View at: Google Scholar
 Y. L. Zhao, Y. X. Wang, W. Wang, L. Tang, Q. Liu, and G. Cheng, “Modeling of rheological fracture behavior of rock cracks subjected to hydraulic pressure and far field stresses,” Theoretical and Applied Fracture Mechanics, vol. 101, pp. 59–66, 2019. View at: Publisher Site  Google Scholar
 Y. L. Zhao, L. Y. Zhang, J. Liao, W. J. Wang, Q. Liu, and L. Tang, “Experimental study of fracture toughness and subcritical crack growth of three rocks under different environments,” International Journal of Geomechanics, vol. 20, no. 8, article 04020128, 2020. View at: Publisher Site  Google Scholar
 Y. L. Zhao, Y. Wang, W. Wang, W. Wan, J. Tang, and J. Tang, “Modeling of nonlinear rheological behavior of hard rock using triaxial rheological experiment,” International Journal of Rock Mechanics and Mining Sciences, vol. 93, pp. 66–75, 2017. View at: Publisher Site  Google Scholar
 Y. Zhao, L. Zhang, W. Wang, J. Tang, H. Lin, and W. Wan, “Transient pulse test and morphological analysis of single rock fractures,” International Journal of Rock Mechanics and Mining Sciences, vol. 91, pp. 139–154, 2017. View at: Publisher Site  Google Scholar
 M. Wenqiang, L. Li, T. Yang, G. Yu, and Y. Han, “Numerical investigation on a grouting mechanism with slurryrock coupling and shear displacement in a single rough fracture,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 8, pp. 6159–6177, 2019. View at: Publisher Site  Google Scholar
 L. Ping, S. Zhaowu, and Z. Guiwang, “Characterization and experimental study on the permeability of rocksamples during complete stressstrain course,” Journal of University of Science and Technology of China, vol. 32, no. 6, pp. 678–684, 2002. View at: Google Scholar
 Z. Shining and B. Lin, “The theory of gas flow and storage in coal seams,” Tech. Rep., China Coal Industry, Beijing, 1999. View at: Google Scholar
 A. B. Comsol, COMSOL multiphysics version 3.2, COMSOL, Stockholm, 2005.
Copyright
Copyright © 2020 Xiao Fukun 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.