Thermal-Hydraulic-Mechanical Coupling Simulation of CO 2 Enhanced Coalbed Methane Recovery with Regards to Low-Rank but Relatively Shallow Coal Seams

: CO 2 injection technology into coal seams to enhance CH 4 recovery (CO 2 -ECBM), therefore presenting the dual beneﬁt of greenhouse gas emission reduction and clean fossil energy development. In order to gaze into the features of CO 2 injection’s inﬂuence on reservoir pressure and permeability, the Thermal-Hydraulic-Mechanical coupling mechanism of CO 2 injection into the coal seam is considered for investigation. The competitive adsorption, diffusion, and seepage ﬂowing of CO 2 and CH 4 as well as the dynamic evolution of fracture porosity of coal seams are considered. Fluid physical parameters are obtained by the ﬁtting equation using MATLAB to call EOS software Refprop. Based on the Canadian CO 2 -ECBM project CSEMP, the numerical simulation targeting shallow low-rank coal is carried out, and the ﬁnite element method is used in the software COMSOL Multiphysics. Firstly, the direct recovery (CBM) and CO 2 -ECBM are compared, and it is conﬁrmed that the injection of CO 2 has a signiﬁcant improvement effect on methane production. Secondly, the inﬂuence of injection pressure and temperature is discussed. Increasing the injection pressure can increase the pressure difference in the reservoir in a short time, so as to improve the CH 4 production and CO 2 storage. However, the increase in gas injection pressure will also lead to the rapid attenuation of near-well reservoir permeability, resulting in the weakening of injection capacity. Also, when the injection temperature increases, the CO 2 concentration is relatively reduced, and the replacement effect on CH 4 is weakened, resulting in a slight decrease in CBM production and CO 2 storage.


Introduction
Since the industrial revolution, human activity has resulted in upraising the concentrations of greenhouse gases, causing an anthropogenic influence on global warming and climate change.If no reduction measures are taken, the global atmospheric temperature will rise by more than 1.5 °C and 2 °C by 2030 and 2050, respectively, according to the global average annual rate of anthropogenic CO 2 emissions from 2011 to 2020 (38.9 billion tons CO 2 /year) [1].Environmental repercussions of global climate change brought on by rising human-caused emissions of heat-trapping greenhouse gases are already ubiquitous, such as accelerated sea level rise, sea ice loss, melting glaciers and ice sheets, more intense heat waves, and plant and animal geographic ranges shifts.Human-caused global climate change is irreversible and will worsen in the coming decades [2].There is almost a universal agreement that the world cannot sit back and wait for definitive answers on this issue, and that preventive and mitigating measures must be implemented concurrently.CO 2 is recognized as by far the most significant greenhouse gas due to its relative abundance when compared to the other greenhouse gases, accounting for around 64% of the heightened "greenhouse effect".Therefore, reducing CO 2 emissions to the atmosphere is a significant challenge in mitigating anthropogenic effects on climate change.
Currently, heavy energy-consuming sectors are moving to developing nations [3], which results in the rising phase of CO 2 emissions in developing countries.The CO 2 emission reduction at the expense of reducing the amount of industrial activity harms economic development, generating two competing demands to keep both economic growth and substantial CO 2 emission reduction for developing countries, especially for China, the major developing country in the world [4].Carbon Capture, Utilization and Storage (CCUS), as one pathway to carbon neutrality goals, can meet these two competing needs.Nowadays, more than 100 countries have proposed corresponding carbon neutrality programs.According to the prediction of existing research, fossil energy will still play an important role by 2050, and Carbon Capture, Utilization, and Storage (CCUS) is the only technical choice to achieve large-scale zero-emission utilization of fossil energy [5].
CO 2 enhanced coalbed methane recovery (CO 2 -ECBM) is an important technical approach of CCUS, which has the dual benefits of greenhouse gas emission reduction and clean fossil energy recovery [6].Within the coal matrix, the adsorption capacity of CO 2 is higher than that of CH 4 so that CO 2 can replace most of the originally adsorbed CH 4 .Thus, the latter is desorbed and swept out.Simultaneously, along the run-through cleats and cracks in coal, CO 2 injection can maintain reservoir energy by improved pressure gradient so that the partial pressure of CH 4 is relatively compensated.CO 2 storage mechanism and storage capacity, and the effect of increasing production of CH 4 are the key technical problems of CO 2 -ECBM [7].
At present, CO 2 -ECBM technology is not very mature, and it is still in the pilot test or experimental mining stage in the world, and has not achieved large-scale commercial mining appraisal.Table 1 illustrates the existing CO 2 -ECBM field experiments and their details.It can be seen from Table 1 that the permeability of coal seams in the United States, Canada and other countries are better, which is one to two orders of magnitude higher than that of Chinese coal.In addition, due to the good adsorption of high-rank coal, most of the experiments are based on middle and high-rank coal seams, and there are few experiments on shallow low-rank coal.Among them, the CSEMP (CO 2 Storage and Enhanced Methane Production) project on which this study is based may be a relatively rare shallow low-rank coal CO 2 injection test project.From the Langmuir adsorption curve, we have known that the CO 2 adsorption capacity of low-rank coal is very small when the pressure is low, and when the pressure increases, the adsorption capacity increases slowly, but the adsorption capacity of high-rank coal is very large when the pressure is low, and the adsorption capacity increases greatly with the increase of pressure.However, it is worth mentioning that the CO 2 /CH 4 adsorption ratio of low-rank coal is much larger than that of high-rank coal, which highlights the superiority of CO 2 replacement of CH 4 in low-rank coal.
The gas migration in the CO 2 -ECBM process is a very complex process, including the competitive adsorption, diffusion, seepage, coal deformation and temperature exchange of coalbed methane and injected CO 2 gas.Based on the laboratory scale, previous scholars have carried out many studies on the adsorption characteristics of CO 2 and CH 4 , the mechanism of gas flow in coal as well as the replacement efficiency and utilization rate of CO 2 by coal.Barbara Dutka [8] carried out an experimental research into the sorption properties of coal briquettes under the various states of stress, proved that the existence of a relation between the stress state, sorption capacity and volumetric deformation of coal.Sun et al. [9] carried out an experimental study on the displacement of methane by CO 2 using eight bituminous and anthracite coals, to evaluate the displacement efficiency and utilization ratio of CO 2 .Zheng et al. [10] explored the characterizations of methane desorption and CO 2 displacement under different CO 2 injection pressures or different experimental temperatures using the self-designed NMR-based CO 2 displace methane equipment.Xu et al. [11] proposed a new multi-scale dynamic diffusion-percolation model to compare the old and new model analysis, as well as carefully studying the mechanism of gas flow in coal.Chen et al. [12] studied the the influence of coal water content and gas injection pressure on CH 4 replacement rate and CO 2 injection ratio based on the adsorption characteristics of coal seam gas injection.In addition to conducting physical experiments, many scholars have conducted indepth research on CO 2 -ECBM using numerical simulation.V. Vishal et al. [25] carried out research work to understand the technical feasibility of CO 2 driven enhanced CBM recovery in Indian coals using a commercial reservoir simulator, COMET.Hema J. Siriwardane et al. [26] performed a historical fit to the Pump Canyon ECBM field test in the United States using PSU-COALCOMP, where the anisotropy of the reservoir was considered in the model, and the study showed that a more accurate set of parameters could be determined using numerical simulations.The gas migration in the CO 2 -ECBM is a process of multi-phase fluid and multi-physical field coupling.
In recent years, scholars have also done some work on CO 2 -ECBM multi-physics coupling.The effect of the numerical case is different based on the difference in coal reservoir structure and physical properties.Wu et al. [27] established a gas-solid coupling model to study the migration of binary gases in coal reservoirs, and clarified the distribution of coalbed methane concentration in CO 2 -ECBM, but did not consider gas diffusion and competitive adsorption of binary gases.Wang et al. [28] carried out largescale laboratory tests and deduced the seepage diffusion control equation, but did not consider the deformation of coal.Fang et al. [29,30] established the coupling equations of the diffusion-adsorption-seepage-heat transfer fields of gas.The displacement processes under different pressures and temperatures are simulated by COMSOL.Subsequently, the stress field model and seepage field model were derived, and the numerical simulation of permeability evolution was carried out in COMSOL software.Hemant Kumar et al. [31] established a coupled finite element (FE) model of binary gas flow, diffusion, competitive sorption and permeability change to explore the effect of CO 2 injection on net recovery, permeability evolution and injectivity in uniform and homogeneously permeable reservoir.Wang et al. [32] built a dual porosity/permeability model through accurately expressing the volumetric strain of matrix and fracture which aims to reveal the reservoir permeability evolution during the process of CO 2 -enhanced coalbed methane (CO 2 -ECBM) recovery.Bai et al. [33] established a fluid-solid coupling model considering gas-water phases and Klinkenberg effect of CO 2 injection into coal seam for enhancement of gas drainage, parameters of gas pressure, gas content and gas extraction rate after coal seam CO 2 injection were analyzed by using COMSOL software and applied in engineering test.Erlei Su et al. [34] established an improved fully coupled gas migration model for CBM extraction.The permeability rebound and recovery times as well as rebound values are proposed to accurately quantify permeability evolution during CBM extraction.Also, the evolution of coal seam permeability under different CO 2 injection method is discussed.Most of the above scholars only consider a single or several physical fields in their research.In addition, many scholars [29,32,35] have studied the Thermal-Hydraulic-Mechanical coupling of CO 2 -ECBM, but most of these studies are based on deep unminable coal seams.At present, there are few numerical simulation studies on the Thermal-Hydraulic-Mechanical multi-physical coupling of shallow low-rank coal.The buried depth of the deep unminable coal seam is large.Under geological conditions of high temperature and high pressure, CO 2 is likely to exist in the supercritical state, so the change in its physical parameters may not be obvious.Therefore, previous studies regard fluid physical parameters as constants,which is barely acceptable.However, in the temperature and pressure range of shallow low-rank coal, the physical parameters of CO 2 are prone to mutation.Therefore, it is very inadvisable to regard the physical parameters as constants.
This study is based on the Canadian CSEMP shallow low-rank coal project to carry out CO 2 -ECBM numerical simulation research.The porosity-permeability evolution model is substituted into the fluid migration equation, combined with the reservoir deformation equation, the energy conservation equations and the fluid physical parameter evolution equation.Then, a set of coupled partial differential equations is derived.The COMSOL Multiphysics simulation software is used to solve the finite element method.Firstly, this study compared CBM (directly mining) with CO 2 -ECBM, affirms the advantages of CO 2 -ECBM.The evolution of reservoir pressure, permeability and adsorbed gas content in CO 2 -ECBM and CBM engineering are analyzed.Subsequently, the effects of different injection pressures and temperatures on CO 2 -ECBM were analyzed.The study has certain guiding significance for the preliminary design of the follow-up CO 2 -ECBM project based on shallow low-rank coal.

Fundamental Assumption
The CO 2 -ECBM process in CH 4 -bearing coal is the fully coupled process of the Thermal-Hydraulic-Mechanical fields, which involves the competitive adsorption of CH 4 and CO 2 , the deformation of the coal reservoir, and the heat exchange of the gas and coal skeleton.The following assumptions are required for establishing these fully coupled models: (1) Coal is a homogeneous isotropic body; (2) CH 4 is saturated in the reservoir; (3) The seepage and diffusion of CH 4 and CO 2 conform to Darcy's law and Fick's law, respectively; (4) The deformation of the coal reservoir is infinitesimal deformation.

Governing Model of the Hydraulic Field
The mass conservation equation of CO 2 and CH 4 gas seepage and diffusion in coal seam is [36]: where, m g is the gas mass in per unit volume of coal, kg/m 3 ; ρ i is the gas density, kg/m 3 ; q is seepage velocity, m/s The seepage of gas in coal seams obeys Darcy's law: Gas diffusion in coal seams obeys Fick's first law: Thus, the following can be obtained: where, D i is the diffusion coefficient, m 2 /s; c i is the diffusion mass concentration of gas, kg/m 3 ; Gas mass per unit volume of coal consists of free gas and adsorbed gas [37]: where, ρ c is the density of coal, kg/m 3 ; b i is the Langmuir adsorption constant of gas, MPa −1 ; a i is the langmuir volume constant of gas, m 3 /kg.Thus, available: where, k is permeability, m 2 ; µ i is the dynamic viscosity coefficient of gas, Pa•s; ϕ is porosity, zero dimension;

Governing Model of the Mechanical Field
The deformation of coal is formed under the joint actions of ground stress, gas pressure, gas adsorption/desorption and thermal stress, and the constitutive equation of the mechanical field can be deduced based on the elastic theory of porous media, which is as follows [38]: where, ε ij is the strain tensor component, m; G is the shear modulus, MPa; σ ij is a component of the stress tensor; σ hh is the normal stress component; δ ij is Kronecker symbol; α is Biot's coefficient.Equation ( 7) is the equilibrium equation characterizing the spatial equilibrium state of coal reservoir, and the strain component and displacement component satisfy Equation (8) [39] σ ij,j Based on the above equations, the Navier-Stokes equations which can characterize the stress field of coal reservoir can be derived: where, v, Poisson ratio; K s is the bulk modulus of coal matrix, MPa; K is the bulk modulus of coal, MPa; ε s is volume strain caused by gas adsorption desorption,

Governing Model of the Thermal Field
According to the law of conservation of energy and Fourier's law, the energy conservation equations of skeleton and gas can be obtained [40].
where, (ρC p ) e f f , is the effective heat capacity, J/m 3 •K, η e f f , is the effective convective heat transfer coefficient, J/m 2 •s, λ e f f , is the effective thermal conductivity, W/(m•K): where, C g1 , C g2 is the specific heat capacity of CH 4 , CO 2 , J/kg•K; λ g1 , λ g2 is the thermal conductivity of CH 4 and CO 2 , W/(m•K).

Permeability Evolution Equation
Based on gas pressure, temperature and volume change caused by adsorption and desorption of coal matrix, porosity evolution formula can be derived [41,42]: where, V p is the pore volume of coal, m 3 ; V is the total volume of coal, m 3 ; ϕ 0 is the initial porosity of coal seam; δV s is the volume change of coal skeleton, m 3 ; V s0 is the initial skeleton volume of coal body, m 3 The following can be obtained: A dynamic model of coal seam permeability can be derived by using the cubic theorem between permeability and porosity [43]:

Physical Parameters
The relevant thermal physical parameters change with temperature and pressure.This study is based on shallow low-rank coal with a buried depth of 450 m, which is relatively shallow.Therefore, the temperature range is 280-310 K and the pressure is 0-10 MPa.Refprop software has been used for physical parameters evaluation (Figures 1 and 2).It can be seen that all the thermal physical parameters change nonlinearly with temperature and pressure.Among them, the physical parameters of CO 2 show a sudden change trend under certain temperatures and pressure.Therefore, it is very unreasonable to regard the thermal physical parameters as constants as previous studies.This paper uses Matlab to call Refprop software, then uses polynomial function fitting in Matlab.It is worth mentioning that the Refprop software is based on the most accurate pure fluid and mixture models at present.It implements three models for the thermodynamic properties of pure fluids: equations of state explicit in Helmholtz energy, the modified Benedict-Webb-Rubin equation of state, and an extended corresponding states (ECS) model [44].For CO 2 and CH 4 , high accuracy Helmholtz energy equations of state is used.
The relationship between density, specific heat capacity, thermal conductivity, dynamic viscosity and pressure/temperature is fitted as follows.These equations are used for numerical solutions.
where, T is temperature, K; P i is pressure, Pa; µ i is the dynamic viscosity, Pa•s; C gi is the specific heat capacity, J/(kg•K); ρ i is the density, kg/m 3 ; λ gi is the thermal conductivity, W/(m•K); i = 1, CH 4 , i = 2, CO 2 .Tables 2 and 3 list the corresponding coefficients in the equation.
Cross coupling among CBM is expressed as: In summary, a fully coupled temperature-seepage-diffusion-deformation model for CO 2 -ECBM and CBM considering the variation of physical parameters is established in Figure 3.The thermal stress caused by the change of coal seam temperature affects the stress field of the coal seam skeleton.The heat generated by the energy dissipation in the coal seam skeleton has an impact on the coal seam temperature; The temperature change causes the gas pressure change to affect the gas; the thermal convection between the gas and coal seam skeleton affects the temperature of the coal seam.The porosity change caused by coal seam deformation affects gas flow.Gas pressure change causes coal seam deformation.
In this study, based on COMSOL Multiphysics software, the established mathematical model is solved by full physical field coupling of the finite element method.

Model Validation
The Canadian CSEMP project [45] is a rare approach CO 2 -ECBM micro-test project for shallow low-rank coal that has undergone multiple productions and multi-component injection processes.The test site is located near Buck Lake in Alder Plain, Canada.The test site consists of Wells 12-28, Well 7-28 and Enerplus Pembina 102/7-28-46-7 W5M (Well Injector).In this paper, based on the two-year initial production of well 12-28, a threedimensional model is established to simulate the 100*100 m area around well 12-28.The verification model is shown below Figure 4: Using the bottom hole pressure to fit the gas injection rate, the following results are obtained in Figure 5. Due to the expansion effect caused by the adsorption and desorption of CO 2 and CH 4 , the permeability gradually decreases, as well as the gas rate of injected CO 2 .The numerical simulation results are basically consistent with the measured data, and acceptable deviations are determined by the simplification of reservoir's initial anisotropy, and the ignore of CO 2 phase transition.However, the numerical simulation results and field data have the same trend, which can preliminarily illustrate the reliability of the numerical model.

Simulation Program
We set up three simulation programs.Firstly, this study will compare CBM (directly mining) with CO2-ECBM, in order to affirm the advantages of CO2-ECBM.The evolution of reservoir pressure, permeability and adsorbed gas content in CO2-ECBM and CBM engineering will be analyzed.Subsequently, the effects of different injection pressures and temperatures on CO2-ECBM will be analyzed (Table 4).

Basic Model
The traditional five-point well pattern layout is used as the research object, and a quarter of the extraction area near the bottom of the well is selected as the simulation area (Figure 6).The actual size is 100 × 100 m.A two-dimensional intercept line is set on the diagonal of the model, and observation points B and C are set near the injection well and the production well.The coordinates of point B and point C are (20 m, 20 m) and (80 m, 80 m) respectively.The initial formation pressure of the coal seam is 1.6 MPa, and the temperature of the coal seam is 289 K.The gas injection boundary condition is constant pressure boundary condition.The production well is connected to the surface, and the bottom hole flowing pressure is 0.1 MPa.Table 5 illustrates the numerical simulation core parameters.Figure 7 shows the distribution of reservoir pressure.During the CBM process, the reservoir pressure gradually decreases with time (Figure 7a).Affected by the pressure relief of the production well, the reservoir pressure begins to decrease from the production well and gradually diffuses into the reservoir.During the CO 2 -ECBM process (Figure 7b), reservoir pressure gradually increases with time.At the same production time, CO 2 injection has an enhancement effect on reservoir pressure.When CO 2 is injected, the reservoir pressure near the bottom of the well is greater.As the distance from the well increases, the pressure decreases in a step-wise manner, which indicates that CO 2 mainly accumulates at the bottom of the well and near the well.
According to the Langmuir adsorption theory, the pressure gradient in the reservoir will directly affect the adsorption and desorption of gas in the reservoir (Equation ( 5)).During the CO 2 injection process, a large amount of CO 2 is adsorbed in the near-well reservoir, which will lead to the rapid decrease of reservoir permeability around the bottom of the well.Finally, it will limit the flow capacity of CO 2 and cause a decrease in injection capacity. Figure 8 presents the distribution of gas pressure along the diagonal.After 1500 days of production, the gas pressure near the injection well is reduced to 1.06 MPa (Figure 8a).Taking the initial gas pressure as a reference, when the production is 300 days, 500 days, 1000 days, 1500 days, the gas pressure near the injection well is 87%, 81%, 75%, 66% of the initial pressure, respectively.During the CO 2 -ECBM process, the gas pressure gradually decreases from the injection well to the middle of the reservoir and then to the production well, and the gas pressure increases with time (Figure 8b).The pressure differential between the production well and the injection well is greatly increased compared with the direct exploitation.It is because of a higher pressure differential, gas will accelerate its migration to the production well.Figure 9 illustrates the distribution of gas pressure at reference point B. The reservoir pressure of point B near the injection well decreases with time.The rate of reservoir pressure reduction is relatively slow at the initial stage because point B is far away from the production well.However, as the influence range of the production radius continues to expand, the reservoir pressure of point B continues to decrease at a higher rate (Figure 9a).
During the CO 2 injection process, after an initial decrease, the reservoir pressure increases rapidly and then the increase rate begins to slow down and even slightly decreases (Figure 9b).The injection of CO 2 first increases the pressure of the reservoir near the injection well.With the progress of CO 2 injection, the pressure increase trend continues to spread to the internal part of the reservoir.After a period of time, it will inevitably affect the reservoir pressure at point B, and increase the pressure at this point.With the injection process, the decrease of the reservoir permeability inhibits the injection rate of CO 2 , and the increase of CO 2 pressure will slow down.However, the pressure reduction trend of CH 4 under the influence of the pressure relief of the production well is not be affected.Therefore, the reservoir pressure increases slowly and slightly decreases in the later period.

Variation of CH 4 Concentration
Figure 10 illustrates the distribution of CH 4 concentration.During the direct exploitation (Figure 10a), the CH 4 concentration first changes in the area near the production well and continues to extend into the reservoir.The CH 4 in the reservoir will continuously migrate to the production well because of the pressure differential.However, the pressure differential between the reservoir and the production well is relatively small, the gas migration speed is slow, and the CH 4 concentration changes gradually.During CO 2 injection mining (Figure 10b), due to the injection of CO 2 , the pressure differential from the injection well to the coal seam, to the production well is higher than that of the direct exploitation.Therefore, at the same time, the change in CH 4 concentration alone with CO 2 injection exploitation is higher than that of direct exploitation.Near the injection well, preferential adsorption of CO 2 to the coal seam after injection forces the desorption of adsorbed CH 4 .The desorbed CH 4 migrates to the production well under the effect of pressure differential, thus the CH 4 concentration near the injection well is also reduced.

Variation of Adsorbed Gas Concentration
The core of the CO 2 -ECBM is to achieve the purpose of stable CO 2 storage through the adsorption of CO 2 by coal matrix.The adsorption competition between CO 2 and CH 4 is one of the core problems in the process of CO 2 storage in the coal seam.Especially in the injection stage, with the injection of CO 2 , the bottom hole and near-well reservoir contain a large pressure range.According to the Langmuir adsorption theory, different pressures determine the adsorption capacity of different gases in the coal seam (Equation ( 5)).
Figure 11 illustrates the distribution of adsorbed CH 4 concentration in direct exploitation and CO 2 injection exploitation, as well as adsorbed CO 2 concentration in CO 2 injection exploitation.During direct exploitation, the adsorbed CH 4 concentration mainly changes near the production well, but the change is small.Affected by the pressure relief of the production well, the free CH 4 is first extracted, then the dynamic equilibrium of the gas in the pore is broken, and part of the adsorbed CH 4 is desorbed Figure 11a.However, the desorption of adsorbed CH 4 is not high during direct exploitation.When CO 2 is injected into the coal seam, it will preferentially adsorb in the coal matrix because the adsorption of CO 2 is greater than that of CH 4 , forcing the latter to desorb.It can be obtained from Figure 11b, With the continuous injection of CO 2 , the concentration of adsorbed CH 4 at the injection well and production well continues to decrease.It can also be seen from Figure 11c that as CO 2 is injected, the concentration of adsorbed CO 2 in the reservoir increased dramatically.
In brief, CH 4 desorption is slow and the amount of desorption is small during direct exploitation, which often leads to a small recovery rate of CH 4 .The CO 2 -ECBM technology has the dual benefits of greenhouse gas emission reduction and clean fossil energy development, and it's superiority to conventional exploration technologies.

Variation of Reservoir Permeability
The scientific connotation of the effectiveness of CO 2 -ECBM mainly includes the CO 2 storage mechanism and storage capacity, and the increase rate of CH 4 yield, which are also the basic theoretical and technical problems faced by CO 2 -ECBM.The effectiveness of CO 2 -ECBM also becomes the premise and foundation of the safety and economy of CO 2 -ECBM.The injectivity of CO 2 can be characterized by injection pressure, injection rate, cumulative injection volume and other parameters, which are mainly related to the permeability of the coal seam and its dynamic changes.Therefore, it is necessary to study the permeability evolution of the reservoir.With the injection of CO 2 , the permeability of coal seams will decrease, which often affects the injectivity of CO 2 [46].
The evolution of permeability in the CO 2 -ECBM process is a very complex phenomenon, which is determined by the stress effect and the adsorption swelling effect.On the one hand, effective stress increases and permeability decreases because of cleat compression [47].On the other hand, the matrix expands, the seepage channel is compressed, and the permeability will decrease, along with CO2 adsorption, as CH 4 desorption occurs, matrix shrinkage, seepage channel opens, permeability will increase [46].
During direct exploitation (Figure 12a), the reservoir permeability increases with time, and the closer to the production well is, the higher the permeability will be.This is because with the production, the pressure near the production well decreases, CH 4 desorbs, the matrix shrinks, then the permeability increases.In addition, the increase of effective stress leads to the decrease of permeability.However, the increase of permeability caused by matrix shrinkage deformation dominates.Therefore, the permeability increases compared with the initial permeability.
In the case of CO 2 injection (Figure 12b), for the area near injection well, the decrease of effective stress will lead to the increase of permeability near the gas injection well, while the matrix expansion caused by competitive adsorption will lead to the decrease of permeability.At this time, the expansion deformation caused by the competitive adsorption of CO 2 and CH 4 dominates, and the permeability near the gas injection well decreases rapidly.
For the area near the production well, when the production time is short, the injected CO 2 does not affect the production well.At this time, the stress effect reduces the permeability, and the matrix shrinkage effect caused by CH 4 desorption leads to the increase of permeability, but the matrix shrinkage effect caused by CH 4 desorption dominates.Therefore, the permeability of the production well is larger than the initial permeability.In the process of CO 2 -ECBM, the gas pressure continues to increase with time, and with the increase of injection pressure, the gas pressure in the reservoir changes more obviously (Figure 13).Especially in the early stage of gas injection, increasing the injection pressure can make the gas pressure in the reservoir increase rapidly in a short time, making CH 4 gas easier to produce.The pressure differential and concentration differential are the driving force of the seepage and diffusion of CO 2 and CH 4 in the coal seam.When CO 2 is injected, the reservoir pressure is greatly increased, and the pressure differential is also increased, which promotes the gas to flow and diffuse to the production well faster.This is one of the reasons why CO 2 injection can increase the production of coalbed methane.The permeability near the injection well is lower than the initial permeability, and the permeability decreases with the increase of injection pressure at the same production time.This is because the increase of gas injection pressure will lead to the increase of reservoir adsorption amount, which will lead to the decrease of permeability.The decrease of permeability will limit the flow capacity of CO 2 and weaken the injection capacity.In this study, point B near the injection well is taken to analyze the change of reservoir permeability near the well when the injection pressure changes.It can be seen that when the injection pressure increases, the permeability of a point near the injection well decreases rapidly.When the injection pressure is 6 MPa, the reservoir permeability decreases by 0.6 times of the initial permeability.At the same time, as the injection pressure increases, the pressure differential in the reservoir increases, prompting faster gas seepage and diffusion, and CH 4 production will increase.When the injection pressure is 6 MPa, the CH 4 production increases to 3.79 × 10 5 m 3 , which is 75% higher than that when the injection pressure is 1 MPa.In summary, increasing gas injection pressure can promote CH 4 production and CO 2 storage, but the increase of gas injection pressure will also lead to rapid attenuation of near-well reservoir permeability, resulting in the weakening of injection capacity.We analyzed the change of CH 4 production.The peak value of the change rate of CH 4 production was 0.153 at 3 MPa, indicating that the injection pressure continued to increase after the displacement pressure reached 3 MPa, and the recovery degree did not change much.Consequently, we can take 3 MPa as the optimal injection pressure in this case.Figure 15 illustrates the reservoir temperature distribution.When the injection temperature is the same as the reservoir temperature, the temperature increases near the injection well due to the adsorption of CO 2 , and the temperature near the production well decreases due to CH 4 desorption.After the gas injection, CO 2 adsorption induced an exothermic effect and the reservoir temperature near the injection well rose sharply [48].Meanwhile, the range affected by temperature gradually increased due to heat conduction.CH 4 desorption promoted a reduction in reservoir temperature near the production well.In addition to the thermodynamic effect of Hydraulic-Mechanical interaction and the temperature effect of non-isothermal adsorption, the internal energy change in the reservoir also has the Thermal-Hydraulic coupling of the fluid itself caused by the change of CO 2 phase state, that is, the phenomenon of liquefaction exothermic and gasification endothermic [49].In thermochemistry, it is the energy absorbed or released by a substance during a change of state (phase change) without a change in temperature that we call the latent heat of phase change.The heat which is required melt a certain quantity of a solid at the melting point into a liquid at the same temperature is called the latent heat of fusion.The conversion of the liquid into vapour requires an amount of 'latent heat' which is generally much greater than the latent heat of fusion of the same substance [50].In the process of CO 2 vaporization, the distance between molecules increases, the molecular force does negative work, and the molecular potential energy increases; where the molecular potential energy is the energy that molecules have with respect to their relative positions due to the mutual forces that exist between them.This process may release or absorb large amounts of heat, but this phenomenon has not been considered in this and previous studies.It is clear that increasing the injection temperature will reduce CH 4 production and CO 2 sequestration.Based on the adsorption potential theory, the increase in temperature has a negative impact on the adsorption of CO 2 and CH 4 in the coal seam, that is, the lower the temperature, the higher the amount of CO 2 and CH 4 adsorbed by coal [9].Meanwhile, with the increase in temperature, the main reason for the decrease in the ability of CO 2 to replace CH 4 is that the change rate of CO 2 potential energy and diffusion capacity is more significant than that of CH 4 , resulting in the decrease of the amount of CO 2 adsorbed by coal seam [10].

Conclusions
This paper combines the specific engineering parameters of the CSEMP shallow lowrank coal project and establishes a coupled Thermal-Hydraulic-Mechanical model, which achieves the dynamic evolution of fluid physical parameters and is of great interest as a new guide for the preliminary design of CO 2 -ECBM projects based on shallow low-rank coal.Based on COMSOL Multi-physics software, the finite element method is used to solve the multi-physical field coupling.The following conclusions are obtained: 1.
During direct mining, the reservoir pressure shows a downward trend, and the downward trend gradually slows down and tends to be stable.During CO 2 injection, the reservoir pressure increases with time, increases rapidly at the initial stage and gradually tends to be stable.The injection of CO 2 has a significant enhancement effect on coalbed methane extraction.2.
For direct mining, the increase in permeability caused by matrix shrinkage deformation dominates, and the reservoir permeability gradually increases with time, and the closer to the production well, the greater the permeability.For CO 2 injection, the volumetric strain effect caused by adsorption and desorption is dominant in both production wells and injection wells.The permeability near the injection well decreases rapidly and the permeability near the production well increases.

3.
Increasing the gas injection pressure can promote the production of CH 4 and the storage of CO 2 , but the increase of gas injection pressure will also lead to the rapid attenuation of near-well reservoir permeability and the weakening of injection capacity.When the gas injection pressure reaches 3MPa, the injection pressure continues to increase, and the recovery degree does not change much.4.
When the injection temperature increases, the CO 2 concentration is relatively reduced, and the replacement effect of CH 4 is weakened, resulting in a slight decrease in CBM production and CO 2 burial stock.

Future Work
For future works, the different data that have been yielded by the several adopted approaches can be integrated into a supervised learning process.This artificial intelligence technique has been extensively used to describe geological assets [51][52][53][54][55].It shows its strength by combining different correlated factors and seek their contribution towards the investigated phenomena.

Figure 2 .
Figure 2. CO 2 (a) density, (b) specific heat capacity, (c) thermal conductivity and (d) dynamic viscosity as a function of pressure/temperature.

Figure 5 .
Figure 5. History match of CH 4 production rate during primary CBM recovery.

Figure 8 .
Figure 8. Distribution curve of reservoir gas pressure along the model diagonal (a) CBM (b) CO 2 -ECBM.

Figure 9 .
Figure 9. Variation curve of reservoir pressure at reference point B. (a) Comparison of CO 2 -ECBM and CBM (b) CO 2 , CH 4 partial pressure and total reservoir pressure of CO 2 -ECBM.

Figure
Figure 14a,b illustrate the distribution of cumulative production and point B's permeability change during CO 2 injection exploitation under the pressure of 1 MPa to 6 MPa.The permeability near the injection well is lower than the initial permeability, and the permeability decreases with the increase of injection pressure at the same production time.This is because the increase of gas injection pressure will lead to the increase of reservoir adsorption amount, which will lead to the decrease of permeability.The decrease of permeability will limit the flow capacity of CO 2 and weaken the injection capacity.In this study, point B near the injection well is taken to analyze the change of reservoir permeability near the well when the injection pressure changes.It can be seen that when the injection pressure increases, the permeability of a point near the injection well decreases rapidly.When the injection pressure is 6 MPa, the reservoir permeability decreases by 0.6 times of the initial permeability.At the same time, as the injection pressure increases, the pressure differential in the reservoir increases, prompting faster gas seepage and diffusion, and CH 4 production will increase.When the injection pressure is 6 MPa, the CH 4 production increases to 3.79 × 10 5 m 3 , which is 75% higher than that when the injection pressure is 1 MPa.

Figure 13 .
Figure 13.300 days, 500 days, 1000 days, 1500 days Reservoir pressure cloud diagram (a) The injection pressure is 2 MPa (b) The injection pressure is 4 MPa (c) The injection pressure is 6 MPa.

4. 3 . 2 .
Figures 16 and 17 illustrate the CH 4 cumulative production and CO 2 cumulative storage of three different temperatures.It is clear that increasing the injection temperature will reduce CH 4 production and CO 2 sequestration.Based on the adsorption potential theory, the increase in temperature has a negative impact on the adsorption of CO 2 and CH 4 in the coal seam, that is, the lower the temperature, the higher the amount of CO 2 and CH 4 adsorbed by coal[9].Meanwhile, with the increase in temperature, the main reason for the decrease in the ability of CO 2 to replace CH 4 is that the change rate of CO 2 potential energy and diffusion capacity is more significant than that of CH 4 , resulting in the decrease of the amount of CO 2 adsorbed by coal seam[10].

Table 2 .
The correlation coefficient of the fitted equation for CH 4 .

Table 3 .
The correlation coefficient of the fitted equation for CO 2 .