Abstract
Coal rock is a porous medium composed of organic matter and inorganic minerals, and it is very complex and highly heterogeneous. Coal bed methane (CBM) production is a thermo-hydro-mechanical (THM) coupling process in heterogeneous coal rock. THM coupling numerical simulation on the coal rock by considering the effect of mechanical heterogeneity is rarely reported. We use Weibull’s probability density distribution function to characterize the heterogeneity in elastic modulus of the coal rock, establish a THM coupling 3D finite element model of the coal rock by considering the variation in pore pressure caused by methane desorption, the linear thermal expansion effect, and coal rock skeleton shrinkage and deformation, and analyze variation in permeability, porosity, stress, temperature, and pore pressure within the coal rock representative elementary volume (REV) of variable mechanical heterogeneity with the cross-coupling correlation between permeability and porosity, and thermal field, stress field, and pressure field. The results show that the evolution of porosity and permeability in the coal rock is a THM coupling process related to mechanical heterogeneity, thermal expansion effect, pore pressure change caused by CBM desorption, and stressed deformation in the coal rock skeleton. The permeability and porosity fluctuate within the heterogeneous coal rock. The permeability and porosity fluctuate more frequently in the coal rock with stronger mechanical heterogeneity. The mechanical heterogeneity promotes local stress concentration. The time for variation in the stress through the whole the coal rock REV and the value of the first principal stress increase when the coal rock heterogeneity is enhanced. Under the THM coupling effect, the strong heterogeneity of the coal rock causes fluctuation in the thermal field. The evolution of coal porosity and permeability is a THM coupling process. This study provides theoretical guidance for CBM exploitation.
1. Introduction
China has an abundant CBM reserve, which is m3, and is mainly hosted in Qinshui Basin and Ordos Basin [1, 2]. CBM is highly efficient and clean energy. CBM, shale gas, and tight sandstone gas are all unconventional gas [3]. CBM is stored as an adsorbed state in the coal bed. CBM exploitation is completed through desorption, diffusion, and seepage and is controlled by temperature, pore pressure, and in situ stress. Thus, CBM exploitation is essentially a THM coupling process. The THM multifield coupling process in the coal rock is a hot topic in CBM exploitation. The THM coupling study is of great significance for accelerating gas desorption, improving coal rock permeability, and enhancing CBM recovery.
Previous experiments show that the coal rock permeability is closely related to confined pressure, temperature, and pore pressure. When the coal rock is heating up, the permeability generally shows a downward trend and is also possibly recovered. The variation in the coal rock permeability is the result of the coupling of matrix thermal expansion, gas desorption, and evaporation of residual water within cracks. Feng et al. [4] carried out experiments on mechanical deformation of anthracite in the process of heating from room temperature (i.e., 25°C) to 600°C and found the thermal effect on coal rock deformation in the thermal expansion stage from room temperature to 200°C, in the slow compression and deformation stage at 200-400°C, and in the strengthening compression stage at 400-600°C. The transition from brittle deformation to ductile deformation in anthracite occurs at 400-450°C. The thermal effect on the coal rock mechanical properties was tested under gas-free conditions, which is different from the underground gas-bearing conditions. In fact, the gas always exists dynamically within the coal rock. Due to gas desorption and adsorption in the coal rock, the effect of gas content on coal rock deformation is significant. The study shows that the adsorption capacity of CBM decays when the coal rock temperature rises. The temperature has a significant effect on coal rock’s mechanical properties which are closely related to coal rank, deformation stage, and pore pressure.
Yin et al. [5] developed the true triaxial servo-controlled THM coupling seepage device for experiments on gas-bearing coal rock and carried out a study on the gas permeability in the gas-bearing coal rock during full stress-strain. The results show that the coal rock at a higher temperature has a lower compressive strength and experiences smaller failure strain and is more prone to failure. The gas permeability in the coal rock mainly depends on the stress-strain stage. In the stages of initial compression and linear elastic deformation in the coal rock, the gas permeability drops with increase of the deviatoric stress. In the stages of yield, failure, and residual deformation in the coal rock, the gas permeability increases with increase of the deviatoric stress. The thermal effect on the gas permeability varies during coal rock deformation, and it is more obvious before yield deformation, where the high temperature causes deterioration of permeability. After yield deformation of coal rock, a large number of microcracks are generated in the coal rock and have a great impact on the gas permeability, and the thermal effect on the gas permeability is weakened. Liu et al. [6] carried out a study on the effect of initial adsorption equilibrium pressure and temperature on gas migration in the coal rock in experiments on CBM desorption and diffusion. The results show that the higher temperature is favorable for gas desorption and diffusion, and the gas desorption consumes the heat in the coal rock and causes a temperature drop in the coal rock, which is adverse to gas desorption and diffusion. Cyclic loading affects coal rock permeability due to reversible deformation and damage and crack propagation. Cai et al. [7] monitored the evolution of the gas permeability in semianthracite and anthracite under triaxial cyclic loading through CT scanning, acoustic emission, and P-wave velocity methods. The memory of cracks occurs when the axial stress is unloaded, which is similar to the acoustic emission Kaiser effect. The gas permeability in the coal rock is related to fracture connectivity, geometry, and stress. During loading of stress in the coal rock, the permeability decreases, and the P-wave velocity increases. Huang et al. [8] developed a THM coupling triaxial device with the function of X-ray scanning, and the confining pressure, axial force, injected fluid pressure, and temperature are up to 20 MPa, 400 kN, 10 MPa, and 100°C, respectively. They carried out the experiment on the gas permeability during triaxial stress loading and unloading of coal rock. The results show that the fractures are mainly caused by expansion of initial fractures during cyclic loading of triaxial stress. The fracture area ratio in the coal rock increases linearly with the stress cycle. As the fracture area ratio increases, the gas permeability in the coal rock increases. The experimental results are affected by the size of the coal rock samples. In fact, the small size of rock samples in the laboratory cannot be exactly equivalent to the real situation in the field application, even though some similarity criterion is usually used in the experiment to reduce the errors. Some importance characteristics such as formation heterogeneity and anisotropy cannot be completely included in rock samples with small size.
A lot of efforts have been put into numerical simulation of gas flow in the coal rock. Yang et al. [9] established a gas seepage-stress coupling model by considering gas desorption and adsorption and carried out a simulation. The results show that the confining pressure has a great effect on gas drainage, and the high confining pressure causes strong resistance in the gas flow, which leads to low gas production. Nevertheless, the thermal effect on the gas permeability was ignored. Ma [10] established a THM coupling model of the coal rock by considering the thermal effect during gas extraction. The simulation results show that rising temperature significantly enhances gas desorption which in turn causes deformation of the coal rock, and this is beneficial to improve the gas permeability. Howerver, this model does not consider the effect of mechanical heterogeneity on the porosity and permeability of coal rock. Based on this THM model, Wei [11] established a 2D THM damage coupling model and studied the variation in in situ stress and evolution of coal damage under different lateral pressure coefficients and hole wall warming conditions. Tao et al. [12] established a real three-field and two-way THM coupling mathematical model by applying theories of elasticity, seepage mechanics, and heat transfer. This model reveals the internal relationship between the gas flow field, deformation field, and thermal field in the gas-bearing coal rock. They verified the model by comparing the experimental results to analytical solutions. Hosking et al. [13] proposed a dual pore THM coupling model of nonisothermal gas flow during CO2 sequestration in the coal rock. The analytical solution and experimental results are compared to verify the reliability of the model. Xue et al. [14] developed a full THM coupling model which considers the nonlinear response of CBM production and verified the reliability of the model by matching it with the historic data. Zhao et al. [15] modified the THM coupling model by considering internal expansion and deformation, corrected the equations describing coal rock deformation and permeability, and used the dynamic diffusion model to describe the evolution of gas diffusion in the coal rock matrix. Moreover, a new model of matrix porosity was introduced based on the assumption that the change in the pore volume equals that in the matrix volume. Zhu et al. [16] extended the dual permeability model to the damage mechanics of coal rock caused by gas adsorption and dissolution and obtain the change in permeability when CO2, CH4, and N2 are injected into coal rock.
The heterogeneity of coal rock affects its mechanical properties. Wang et al. [17] carried out the uniaxial acoustic emission test on the black shale sampled at the depth of 1145 m from Well Youqian 1 in Chongqing and compared the results to those from coarse-grained granite samples. The shale samples show a significant difference in the wave velocity, strength, and brittleness, indicating strong microscopic heterogeneity. The total organic carbon content (TOC), mineral composition, and natural fractures are the main factors affecting the mechanical heterogeneity of shale. Wang et al. [18] applied Weibull’s probability density function to describe the mechanical heterogeneity of shale and carried out a numerical simulation on the effects of the in situ stress distribution, strain energy density, fluid pressure, and damage of shale under hydro-mechanical coupling during hydraulic fracturing. Nevertheless, the thermal effect on shale damage was not considered. Zhu et al. [19] proposed a finite element method based on digital images for the study of the rock microscopic heterogeneity, obtained the effect of microscopic heterogeneity on the macroscale mechanical properties, defined the attribute graph according to the pixel ratio of digital imaging, and used the variable elastic modulus and permeability to describe rock heterogeneity. Nevertheless, the model is two-dimensional and cannot reflect the real spatial process.
Coal rock is highly heterogeneous. The effects of mechanical heterogeneity on the temporal and spatial evolution of gas permeability and porosity have been rarely reported. We use Weibull’s probability density distribution function to characterize the elastic modulus heterogeneity of the coal rock, establish a 3D THM coupling finite element model by considering the pressure change caused by methane desorption, the linear thermal expansion effect, and the shrinkage and deformation of the coal rock skeleton, and analyze the variation in permeability and porosity in the coal rock REV with the cross-coupling correlation between the permeability and porosity and the thermal field, stress field, and pore pressure field.
The paper is organized as follows. Firstly, it provides the coal rock THM coupling physical model and assumptions. Secondly, it develops the governing equation, boundary conditions, finite element coupling algorithm of the THM coupling model, and random characterization methods of coal rock mechanical heterogeneity. Thirdly, it illustrates numerical simulation results, including case verification and sensitivity analysis. Finally, it describes the conclusions.
2. Finite Element Model
The solution to the 3D THM coupling problem needs a huge computational cost. To save the computation cost, we establish a physical model in a REV where a borehole with a radius of 0.2 m is set in the center. Here, the tetrahedral elements are used to mesh the RVE in the 3D space, and a total of 34,878 elements are generated. To simulate stress concentration around the borehole, local meshes are refined around the borehole (Figure 1). The discrete term of time derivative is expressed with an implicit format and an adaptive time step, and the Petrov-Galerkin finite element method ensures computational stability of the THM interaction model [20, 21]. The coal rock THM coupling finite element model is meshed (Figure 1). It is assumed that (1) the gas permeability is isotropic; (2) the fluid is slightly compressible, and the THM coupling effect on in situ stress, pore pressure, and temperature are considered; (3) the coal rock damage caused by variation in fluid pressure, rock deformation, and temperature is negligible; (4) the heterogeneity of thermal property is not considered in our simulation, and we only consider the mechanical heterogeneity of elastic modulus of coal rocks in this paper.
3. Mathematical Model
3.1. Mechanic Equation of Coal Rock Deformation
According to the theory of elastic mechanics, the mechanical equilibrium equation for coal rock deformation is expressed as follows:
Assuming that the coal rock is an isotropic medium and the coal rock is under linear elastic deformation, the stress-strain relation for coal rock deformation described by the Biot constitutive theory is expressed as follows:
The relationship between strain and displacement in the deformed coal rock is expressed as follows:
According to Equation (1), the volumetric strain is expressed as follows: where .
The volumetric contraction strain of coal rock caused by gas desorption is defined as
By substituting Equations (2)–(5) into Equation (1), we obtain the differential equation of coal rock deformation as follows:
The body force in Equation (6) is not considered in the next simulation because the THM coupling model is 2D, and thus, the value of the body force vector is equal to zero in the numerical simulation.
3.2. Equation of Gas Seepage in the Coal Rock
According to the law of conservation of mass, the differential equation of gas seepage in the coal rock is expressed as follows:
According to the Langmuir equation, the content of methane gas adsorbed in the coal matrix is expressed as follows:
The adsorbed gas content is expressed as follows:
According to the ideal gas law, we have
According to Equation (10), the gas density in the standard state is expressed as follows:
Considering that the gas in the coal rock is in the non-Darcy flow state, and according to the Forchheimer equation, the gas flow is expressed as follows: where and is the porosity of coal rocks.
When the gas is in the Darcy flow state, the gas flow is expressed as follows: where .
Generally, the relation of porosity between mean stress and pressure is expressed as follows:
Integrating both sides of Equation (13), we get
The subscript 0 indicates the initial state of the parameter.
Substituting Equation (4) into Equation (16), we get
According to the Klinkenberg slippage effect during gas flow, the effective gas permeability is related to pressure, and their relationship is expressed as follows:
is the absolute permeability under the very high pore pressure, and the Klinkenberg slippage is negligible; is the Klinkenberg factor related to the mean free path of gas molecules, and the mean free path depends on temperature, pore pressure, and gas molar mass. The Klinkenberg factor is expressed as follows:
According to the experimental estimation, .
The relationship between gas permeability and porosity is expressed as follows:
By substituting Equations (8) to (20) into Equation (7), we get [14]
3.3. Energy Conservation Equation
When the heat carried out by the gas flow is ignored, the total heat flux density is expressed as follows: where .
When the mutual conversion between thermal energy and mechanical energy is ignored, the thermal balance equation is expressed as follows: where , , and .
According to the law of mass conservation of matrix and gas in the coal rock, we have
and ; Equation (23) is transferred to [14]
It is the partial differential equation of the thermal field in the coal rock.
3.4. Mechanical Heterogeneity Characterization
According to the previous scholars’ experience, Weibull’s function can better describe the mechanical heterogeneity of porous media, which is in line with reality compared to other random functions. Weibull’s probability density distribution function is a method to characterize the rock heterogeneity, which is realized by constraint of the rock physical parameters with the Weibull distribution and random assignment of physical parameters with the Monte Carlo method [11, 18]. Weibull’s probability density function is expressed as follows:
The smaller value of indicates stronger heterogeneity. Here, the random distribution of the coal rock elastic modulus is considered. Figure 2 illustrates the 3D distribution of the elastic modulus in the RVE and the contour map at different planes. Heterogeneity degree , and the average elastic modulus .
(a)
(b)
3.5. Model Verification
According to the rock mechanics, there is an analytical solution to the stress field around the borehole in the homogeneous formation. The basic parameters are listed in Table 1. The analytical solutions and finite element numerical solutions to thecomponent in the stress field are shown when the bottom hole pressure is 30 MPa, 40 MPa, and 50 MPa [22–24], respectively (Figure 3). The solutions have a high degree of agreement, which verifies the reliability of the finite element numerical solution under the HM coupling conditions.
4. Results and Discussion
The parameters of the finite element model are listed in Table 1. The 3D stress field, thermal field, and pressure field are obtained by simulating the THM coupling process in the homogeneous coal rock and the coal rock with the heterogeneity degree , 5, 3, and 1. The variation in the porosity and permeability around the borehole with time and distance from the borehole under different conditions is analyzed.
4.1. Homogeneous Coal Rock
The 3D field of the first principal stress in the homogeneous coal rock REV is obtained (Figure 4). The principal stress is relatively higher around the borehole, indicating stress concentration. The stress concentration spreads from the borehole to the coal rock as the CBM is produced. Figure 5 shows the first principal stress in the coal rock during CBM production. When the production time s, shows a downward trend during CBM production. As the distance from the borehole increases, firstly increases and reaches a peak and then decreases to a stable value. When s, shows an upward trend during CBM production, and increases as increases. When m, is stabilized. varies significantly around the borehole, and its increment slows down when m, which is caused by stress concentration.
The 3D thermal field in the homogeneous coal rock REV is obtained (Figure 6). As the CBM is produced, the temperature decreases gradually around the borehole, which is caused by methane adsorption and desorption, and the change in the thermal field spreads out. The temperature of the coal rock during CBM production is shown in Figure 7. The temperature rises as increases, indicating that CBM production leads to cooling down around the borehole. When s, the thermal field in the coal rock is not changed.
The 3D pore pressure field in the homogeneous coal rock REV is obtained (Figure 8). The pore pressure field shows a circular distribution around the borehole. As the CBM is produced, the pressure wave gradually spreads out from the borehole. The pore pressure in the coal rock during CBM production is shown in Figure 9. As increases, the pressure gradually increases, indicating that production causes pressure depletion around the borehole. When s, the pore pressure shows a downward trend; when s, the pore pressure shows an upward trend.
The permeability and porosity in the coal rock with the homogeneous elastic modulus during CBM production are obtained (Figure 10). As the CBM is produced, the permeability and porosity decrease around the borehole. In the early production stage, the permeability and porosity decline faster. When s, the decrease slows down. The coal rock permeability and porosity decrease as increases. In the area within 0.15 m of the borehole, the permeability and porosity decrease slower when m and decrease faster when m.
(a)
(b)
4.2. Heterogeneous Strata
4.2.1. Heterogeneous Coal Rock
(1) Heterogeneity Degree . The 3D principal stress field (Figure 11) around the borehole in the coal rock shows the random distribution, which is caused by elastic modulus heterogeneity. As the CBM is produced, the variation in the principal stress spreads out. Initially, the principal stress variation area is concentrated around the borehole. When s, the variation in spreads over the entire heterogeneous coal rock REV. The heterogeneity of the stress field is transmitted from the top to the bottom. Figure 12 shows the during CBM production. fluctuates as the distance from the borehole increases. When s, decreases as the CBM is produced, which is similar to the phenomenon in the homogeneous coal rock (Figure 5). is relatively high around the borehole and is stabilized as increases. When s, increases during CBM production. increases as increases. When m, fluctuates around a stable value, which is caused by mechanical heterogeneity.
Figure 13 shows the 3D thermal field in the coal rock during CBM production. In the initial production stage, the temperature near the borehole shows a random distribution. When s, the temperature is uniform over the coal rock REV. As production goes on, the variation in the thermal field gradually spreads out. Compared with the stress field, the thermal field does not show obvious heterogeneity. This is because only the elastic modulus heterogeneity closely related to the stress field is considered in this study, and it has little effect on the thermal field. Figure 14 shows the temperature in the coal rock during CBM production, which is the same as that in the homogeneous coal rock. The weak mechanical heterogeneity has little effect on the coal rock thermal field.
Figure 15 shows the 3D pore pressure field in the coal rock during CBM production. The pore pressure field around the borehole shows the circular distribution. As CBM production goes on, the pressure wave gradually spreads out. When s, the pressure wave is stabilized and spreads from the top to the bottom of the REV, which is similar to the phenomenon in the homogeneous coal rock. Pressure heterogeneity propagation is not observed. This is because only the elastic modulus heterogeneity is considered and it has little effect on the pore pressure field. The pore pressure during CBM production is shown in Figure 16, which is the same as Figure 9. The weak mechanical heterogeneity has little effect on the coal rock pressure field. Another reason is that we do not consider the heterogeneous parameters of pore pressure equations in this paper. The mechanical heterogeneity of elastic modulus of coal rock has a negligible effect of pore pressure distribution despite the existence of multiple field coupling process.
The permeability and porosity in the coal rock with heterogeneous elastic modulus are obtained (Figure 17). The permeability and porosity in the heterogeneous coal rock fluctuate and vary as increases. As CBM production goes on, the permeability and porosity of the coal rock show a downward trend.
(a)
(b)
(2) Heterogeneity Degree . The principal stress field in the coal rock shows the random distribution (Figure 18). when is higher than that when . The reason is that stronger mechanical heterogeneity causes local stress concentration and an increase in principal stress. As CBM production goes on, the variation in the stress field gradually spreads out. Initially, the stress only varies around the borehole. When s, the wave in the stress field spreads over the coal rock REV. The stress field heterogeneity is transmitted from the top to the bottom of the REV. in the coal rock during CBM production is shown in Figure 19. fluctuates as increases, which is the same as the phenomenon when . Compared with that when , when fluctuates more frequently, which is caused by stronger mechanical heterogeneity.
When , the gas permeability and porosity fluctuate as increases (Figure 20). As CBM production goes on, both gas permeability and porosity show a downward trend. Moreover, they fluctuate more frequently, which is caused by stronger heterogeneity.
(a)
(b)
(3) Heterogeneity Degree . When , the first principal stress field shows the random distribution around the borehole (Figure 21). when is higher than that when . As CBM production goes on, the variation in the first principal stress gradually spreads out. Initially, only varies around the borehole. When s, the variation in the stress field spreads over the entire coal rock REV. The first principal stress in the coal rock during CBM production is obtained (Figure 22). As increases, fluctuates even more frequently, which is caused by stronger mechanical heterogeneity.
When , the gas permeability and porosity fluctuate more frequently (Figure 23). As CBM production goes on, both gas permeability and porosity in the coal rock show a downward trend.
(a)
(b)
(4) Heterogeneity Degree . When , shows the random distribution around the borehole (Figure 24). Stronger mechanical heterogeneity leads to a higher value of . As CBM production goes on, the variation in the stress area spreads out. When s, variation in spreads over the entire REV.
When , fluctuates most frequently (Figure 25). As increases, the thermal field also fluctuates, which is caused by the THM coupling effect (Figure 26). When the mechanical heterogeneity is strong (high value of ), the effect of the mechanical heterogeneity on the thermal field is negligible.
When , the gas permeability and porosity fluctuate more frequently. As the CBM is produced, both permeability and porosity in the coal rock show a downward trend (Figure 27).
(a)
(b)
5. Discussions and Analysis
The computation cost is increased with the increase of the degree of heterogeneity. This is because the stiffness matrix is not symmetric at different levels of heterogeneous parameter . With the increase of the degree of heterogeneity, the randomness of the stiffness matrix is enhanced, which leads to the extended computation time.
The 3D THM coupled model in this paper considers the complex process such as desorption and adsorption, thermal expansion, and pore pressure variation in heterogeneous coal rock. The mechanical heterogeneity of elastic modulus is incorporated into the THM model by Weibull’s random function. However, the thermal parameters such as thermal conductivity and specific heat capacity are not considered in the THM model, which can be done as the next work. Besides, the coal cleats widely distribute in the coal body. Coal rock body is simplified as a double-medium model, including a plenty of micropores and cleats. A more complex THM model can be extended by considering the cleat deformation and thermal expansion in future.
6. Conclusions
We characterize the mechanical heterogeneity of the coal rock based on Weibull’s probability density distribution function and establish a THM coupling 3D model in the finite element commercial software COMSOL by considering variation in the pore pressure caused by methane gas desorption and adsorption and rock deformation and thermal expansion effect. We carry out a numerical simulation on the temporal and spatial evolution of gas permeability and porosity through the correlation between gas permeability and porosity and the THM coupling effect. The conclusions are as follows: (1)The mechanical heterogeneity of coal rock has a significant effect on gas permeability and porosity. Compared with homogeneous coal rock, the heterogeneous coal rock shows the fluctuating distribution of permeability and porosity. The porosity and permeability vary significantly in different positions of the coal rock.(2)Similar to the homogeneous coal rock, the heterogeneous coal rock has decreasing permeability and porosity as CBM is produced. The time for variation in the first principal stress transmitted to whole coal rock REV is longer, and the first principal stress increases, because the mechanical heterogeneity enhances the local stress concentration in the coal rock.(3)When the coal rock has weak heterogeneity, the effect of mechanical heterogeneity on the thermal field and pressure field is negligible. When the coal rock has strong heterogeneity, the mechanical heterogeneity causes fluctuation in the thermal field due to the THM coupling effects.(4)The evolution of porosity and permeability in the heterogeneous coal rock is a complex THM coupling process, which is closely related to mechanical heterogeneity, thermal expansion effect, pressure change caused by methane gas desorption, and stressed deformation of coal rock skeleton.(5)In CBM exploitation, we should pay attention to the effects of mechanical heterogeneity on the porosity and gas permeability which are key parameters affecting CBM production. The mechanical heterogeneity varies within the coal bed and leads to multiple orders of magnitudes of gas permeability under the THM coupling effect. Thus, exploitation techniques should be optimized in the coal bed with different gas permeability and porosity.(6)The heterogeneity of thermal parameters such as thermal conductivity and specific heat capacity, rock porosity, and permeability on THM coupling process can be incorporated into the current model. In addition, the double-porosity THM model with consideration coal cleats can be developed on basis of the model in this paper.
Nomenclature
| : | Stress tensor (MPa) |
| : | Volumetric force vector (MPa) |
| : | Strain tensor (dimensionless) |
| : | Shear modulus (MPa) |
| : | Poisson’s ratio (dimensionless) |
| : | Gas pressure (MPa) |
| : | Biot coefficient (dimensionless) |
| : | Kronecker delta symbol (dimensionless) and when and when |
| : | Temperature (K) |
| : | Volumetric strain caused by adsorption (dimensionless) |
| : | Linear thermal expansion coefficient (1/K) |
| : | Displacement vector (m) |
| : | Volumetric strain, (dimensionless) |
| : | Average stress (MPa) |
| : | Effective stress (MPa) |
| : | Adsorbed gas content (m3) |
| : | Volumetric strain coefficient caused by gas adsorption (1/m3) |
| : | Gas content (m3) |
| : | Production time (s) |
| : | The first principal stress (MPa) |
| : | Distance from the borehole center (m) |
| : | Gas density (kg/m3) |
| : | Seepage velocity vector (m/s) |
| : | Source term (m3/s) |
| : | Porosity (dimensionless) |
| : | Gas density in the standard state (kg/m3) |
| : | Coal rock density (kg/m3) |
| : | Langmuir volume constant at a temperature of Tt (m3) |
| : | Pressure constant at a temperature of Tt (MPa) |
| Tar: | Absolute reference temperature in the unstressed state (K) |
| : | Referenced temperature in the gas desorption test (K) |
| : | Coal rock temperature (K) |
| : | Pressure coefficient (dimensionless) |
| : | Temperature coefficient (dimensionless) |
| : | Gas molecular mass (kg/mol) |
| : | Gas universal constant (J/(kg·K)) |
| : | Gas pressure in the standard state (MPa) |
| : | Temperature in the standard state (K) |
| : | Viscosity (mPa·s) |
| : | Gas permeability (m2) |
| : | Drag coefficient (dimensionless) |
| : | Darcy correction coefficient (dimensionless) |
| : | Thermal flux density (J/(m2·s)) |
| : | Specific heat capacity of gas at constant volume (J/(kg·K)) |
| : | Darcy flow rate (m/s) |
| : | Coal rock thermal conductivity (W/(m·K)) |
| : | Matrix thermal conductivity (W/(m·K)) |
| : | Gas thermal conductivity (W/(m·K)) |
| : | Variables such as elastic modulus and permeability |
| : | Shape parameter (dimensionless) |
| : | Average of a random variable |
| : | Exponential function (dimensionless) |
| : | Specific heat capacity of the coal rock medium filled by gas (J/(kg·K)) |
| : | Coal density (kg/m3) |
| : | Specific heat capacity of the matrix at constant volume (J/(kg·K)) |
| : | Gas thermal expansion coefficient at constant pressure (1/K). |
Data Availability
Datasets related to this article can be found by contacting the corresponding author.
Conflicts of Interest
The authors declare that they have no conflict of interests.
Acknowledgments
This work was supported financially by the Beijing Natural Science Foundation Project (3222030), the National Natural Science Foundation of China (51936001 and 52174045), the PetroChina Science and Technology Innovation Foundation Project (2021DQ02-0201), and the Science and Technology Project of PetroChina Company Limited “Research on Effective Reservoir Prediction and High-efficiency Well Optimization Technology for Deep/Ultra-Deep Oil and Gas Reservoirs” (2021DJ1003).