Strength Reduction of Coal Pillar after CO 2 Sequestration in Abandoned Coal Mines

CO2 geosequestration is currently considered to be the most effective and economical method to dispose of artificial greenhouse gases. There are a large number of coal mines that will be scrapped, and some of them are located in deep formations in China. CO2 storage in abandoned coal mines will be a potential option for greenhouse gas disposal. However, CO2 trapping in deep coal pillars would induce swelling effects of coal matrix. Adsorption-induced swelling not only modifies the volume and permeability of coal mass, but also causes the basic physical and mechanical properties changing, such as elastic modulus and Poisson ratio. It eventually results in some reduction in pillar strength. Based on the fractional swelling as a function of time and different loading pressure steps, the relationship between volumetric stress and adsorption pressure increment is acquired. Eventually, this paper presents a theory model to analyze the pillar strength reduction after CO2 adsorption. The model provides a method to quantitatively describe the interrelation of volumetric strain, swelling stress, and mechanical strength reduction after gas adsorption under the condition of step-by-step pressure loading and the non-Langmuir isothermal model. The model might have a significantly important implication for predicting the swelling stress and mechanical behaviors of coal pillars during CO2 sequestration in abandoned coal mines.


Introduction
Greenhouse gas emissions are the most important contributor to global climate change.Among all kinds of greenhouse gass, the contribution rate of CO 2 to greenhouse efficiency was 63% [1].According to statistics, fossil fuel combustion and industrial emissions of CO 2 accounted for about 78% of total CO 2 emission in the ten years from 2000 to 2010 [2].Currently, CO 2 storage in oil and gas fields, brine water, deep unmined coal seams and deep sea are considered to be effective and practical ways to reduce atmospheric CO 2 level, helping to slow global climate change and temperature rise trends, in the event that fossil fuels remain in use as a primary energy source.In the 2013 technology roadmap, the international energy agency (IEA) proposed an integrated approach to drop greenhouse gas emission by reducing the use of fossil fuels, improving energy efficiency, implementing new energy sources and carbon capture and sequestration (CCS) technology [3].In addition, with the successful exploitation of coalbed methane and shale gas, its considerable economic benefit prompted many countries to begin in order to regard shale gas, coalbed methane, as alternative unconventional energy.Many scholars have proposed using CO 2 to enhanced coalbed methane (ECBM) production, and the presence of CO 2 will mechanically weaken the coal and thus create fractures, helping to increase the permeability, improve the coalbed methane production yield and simultaneously sequestrate CO 2 [4][5][6][7][8].In addition, for the residual space volume constituted by goaf areas and principle infrastructures in abandoned coal mines, some researchers proposed CO 2 sequestration in abandoned coal mines following the example of natural gas storage in it [9][10][11][12], which will be a potential option for CO 2 disposal because China has a large number of scrapped coal mines, and 541 key coal mines will be gradually closed by 2020.CO 2 can be stored in abandoned coal mines in three states: adsorbed on the remaining coal, free in empty space or dissolved in mine water.
The adsorption of CO 2 in coal can result in coal matrix swelling due to the fact that it has the highest adsorption potential compared with other fluids such as CH 4 and N 2 .Currently, two mechanisms are applied to explain the adsorption-induced swelling in coal.On the one hand, several authors have widely researched the polymer structure, degree of cross-linking, three-dimensional polymeric network structures, as well as flexibility characteristics of coal macromolecules from the perspective of chemistry and molecules, and consider that the lower molecular-weight solvent, such as CO 2 and CH 4 , can enter the macromolecule cross-linked polymer mesh, causing the coal matrix macromolecular structure rearrangement, resulting in swelling [13][14][15][16][17].One the other hand, some scholars attribute the swelling being due to the formation of microfractures as the result of different pore systems, maceral components and mineral stiffness [18][19][20][21].Hol et al. [22] considered CO 2 induced both reversible (i.e., adsorption-induced swelling and elastic compression) and irreversible (i.e., adsorption-induced microfracturing) strains under unconfined conditions.
Both ECBM and CO 2 sequestration in coal seams are concerned with the coal reservoir permeability behavior, the adsorption-induced swelling of coal matrix can compress the pore space and cleat system to result in the distinct decrease of the permeability of coal mass, and several models have been proposed from the consideration of effective stress, cleat volume compressibility, gas sorption-induced strain effect as well as pressure pulse decay [23][24][25][26][27][28][29].Ranjith et al. [30] developed a triaxial equipment to study the gas fluid flow and found that coal mass permeability for CO 2 decreased largely with the increase of effective stress than that of N 2 , due to the matrix swelling by CO 2 adsorption in coal.In addition, Verma and Sirvaiya [31] utilized the artificial neural network (ANN) to predict the Langmuir volume and pressure constants during CO 2 adsorption in coal, and that the ANN method was more accurate than other models in their study.
The effects of CO 2 adsorption on the strength of coal have been studied widely [32][33][34].According to Gibbs' theory, when a more reactive, higher-chemical-energy adsorbate is used to displace the original adsorbate in solid adorbent, the surface energy of rock mass would reduce, which can lead to some reduction in initiation tension stress for fracturing and eventually the coal becomes more prone to damage.In addition, considering the thermodynamics of adsorption of gases in porous solids, the changes in surface energy at the interface between the gas adsorbate and solid adsorbent result in swelling through the conversion between surface free energy and elastic strain energy [35,36].Based on the theory and experiment proposed by Meyers, a theoretical model was derived by Pan and Connell [37] through the energy balance approach.Hol et al. [38] developed a thermodynamic model based on statistical mechanics, and the model combined adsorption in a stress-supporting solid with the poroelastic to derive the relationship of stress-strain-sorption of coal under unconfined swelling condition.Liu et al. [39] revised the model derived by Hol et al. [38], and a corrected expression was obtained based on both statistical mechanics and kinetic approaches.Furthermore, Hol et al. [40] found the apparent bulk modulus determined for CO 2 -equilibrated state was approximately 25% lower compared to the evacuated state through experiment data analysis.Ranjith et al. [32] studied the crack closure, crack initiation and crack damage of coal subjected to saturation with CO 2 .Ranjith and Perera [41] considered the effects of the cheat system on strength reduction of coal after CO 2 adsorption.Perera et al. [42] experimented with adsorption of gaseous and super-critical CO 2 on bituminous coal from the Southern Sydney Basin, Australia, and studied the mechanical properties of coal sample before and after adsorption.The results showed that, compared with the natural state uniaxial compressive strength (UCS), the gaseous CO 2 saturation reduced UCS by 53% and elastic modulus by 36% using gas saturation pressure of 6 MPa.However, the supercritical CO 2 saturation reduced UCS by 79% and elastic modulus by 74% using super-critical saturation pressure of 8 MPa.It is shown that the phase of CO 2 has a significant effect on physical properties of coal, the adsorption capacity, swelling effect and strength parameters.
Coal is a discontinuous structure comprised by many natural fractures and cleats (Figure 1), and coal seams are normally conceptualized by a matchstick model (Figure 2).After the colliery was scrapped, leaving pillars with a large number of cracks by mining action on both sides of the goaf.Meanwhile, CO 2 adsorption causes the coal mass to break down along the cleat system easily due to the fact that the adsorption mainly affects the cohesion of coal, and the reduction of cohesive force leads to the apparent plastic deformation areas.However, for the internal friction angle, it decreases to a certain extent and no longer keeps changes.Pillars with a large number of cleats and fractures act as sealing walls when CO 2 is stored in goafs and drifts.It means that the stability of the pillar decides the safety and sealing of CO 2 sequestration.Based on the above discussions and conclusions, it is significantly important to discuss the strength reduction of the coal pillar after CO 2 injection in the abandoned coal mines.In this paper, we focus our attention on strength reduction of the pillar when CO 2 sequestration in abandoned coal mines, through the relationship of sorption-strain based on unconfined conditions, and swelling-stress under uniaxial conditions.Finally, a strength reduction model is proposed to qualitatively understand the effects of CO 2 adsorption on the strength and failure mechanics of coal pillar in abandoned coal mines.reduced UCS by 53% and elastic modulus by 36% using gas saturation pressure of 6 MPa.However, the supercritical CO2 saturation reduced UCS by 79% and elastic modulus by 74% using super-critical saturation pressure of 8 MPa.It is shown that the phase of CO2 has a significant effect on physical properties of coal, the adsorption capacity, swelling effect and strength parameters.
Coal is a discontinuous structure comprised by many natural fractures and cleats (Figure 1), and coal seams are normally conceptualized by a matchstick model (Figure 2).After the colliery was scrapped, leaving pillars with a large number of cracks by mining action on both sides of the goaf.Meanwhile, CO2 adsorption causes the coal mass to break down along the cleat system easily due to the fact that the adsorption mainly affects the cohesion of coal, and the reduction of cohesive force leads to the apparent plastic deformation areas.However, for the internal friction angle, it decreases to a certain extent and no longer keeps changes.Pillars with a large number of cleats and fractures act as sealing walls when CO2 is stored in goafs and drifts.It means that the stability of the pillar decides the safety and sealing of CO2 sequestration.Based on the above discussions and conclusions, it is significantly important to discuss the strength reduction of the coal pillar after CO2 injection in the abandoned coal mines.In this paper, we focus our attention on strength reduction of the pillar when CO2 sequestration in abandoned coal mines, through the relationship of sorption-strain based on unconfined conditions, and swelling-stress under uniaxial conditions.Finally, a strength reduction model is proposed to qualitatively understand the effects of CO2 adsorption on the strength and failure mechanics of coal pillar in abandoned coal mines.

Adsorption-Induced Swelling Strain
CO2 adsorption-induced volume strain has been studied by many researchers under unconfined conditions or uniaxial strain conditions [28,38,42].In this section, we consider the adsorption-induced strain from the perspective of unconfined conditions.It is easy to measure the volumetric strain value of test coal samples when carrying out a gas adsorption experiment in the reduced UCS by 53% and elastic modulus by 36% using gas saturation pressure of 6 MPa.However, the supercritical CO2 saturation reduced UCS by 79% and elastic modulus by 74% using super-critical saturation pressure of 8 MPa.It is shown that the phase of CO2 has a significant effect on physical properties of coal, the adsorption capacity, swelling effect and strength parameters.
Coal is a discontinuous structure comprised by many natural fractures and cleats (Figure 1), and coal seams are normally conceptualized by a matchstick model (Figure 2).After the colliery was scrapped, leaving pillars with a large number of cracks by mining action on both sides of the goaf.Meanwhile, CO2 adsorption causes the coal mass to break down along the cleat system easily due to the fact that the adsorption mainly affects the cohesion of coal, and the reduction of cohesive force leads to the apparent plastic deformation areas.However, for the internal friction angle, it decreases to a certain extent and no longer keeps changes.Pillars with a large number of cleats and fractures act as sealing walls when CO2 is stored in goafs and drifts.It means that the stability of the pillar decides the safety and sealing of CO2 sequestration.Based on the above discussions and conclusions, it is significantly important to discuss the strength reduction of the coal pillar after CO2 injection in the abandoned coal mines.In this paper, we focus our attention on strength reduction of the pillar when CO2 sequestration in abandoned coal mines, through the relationship of sorption-strain based on unconfined conditions, and swelling-stress under uniaxial conditions.Finally, a strength reduction model is proposed to qualitatively understand the effects of CO2 adsorption on the strength and failure mechanics of coal pillar in abandoned coal mines.

Adsorption-Induced Swelling Strain
CO2 adsorption-induced volume strain has been studied by many researchers under unconfined conditions or uniaxial strain conditions [28,38,42].In this section, we consider the adsorption-induced strain from the perspective of unconfined conditions.It is easy to measure the volumetric strain value of test coal samples when carrying out a gas adsorption experiment in the

Adsorption-Induced Swelling Strain
CO 2 adsorption-induced volume strain has been studied by many researchers under unconfined conditions or uniaxial strain conditions [28,38,42].In this section, we consider the adsorption-induced strain from the perspective of unconfined conditions.It is easy to measure the volumetric strain value of test coal samples when carrying out a gas adsorption experiment in the laboratory under unconfined conditions.As shown in Figure 3, the sample is a cube with the side length l.Two lengths are measured as the base values for volume calculation.One located at parallel (l pa ) to its bedding plane, and the other is perpendicular (l pe ).The parallel and perpendicular displacements were ∆l pa and ∆l pe , respectively.
Minerals 2017, 7, 26 4 of 15 laboratory under unconfined conditions.As shown in Figure 3, the sample is a cube with the side length l.Two lengths are measured as the base values for volume calculation.One located at parallel (lpa) to its bedding plane, and the other is perpendicular (lpe).The parallel and perpendicular displacements were Δlpa and Δlpe, respectively.For the sake of convenience in the computation, it was assumed that parallel lengths and displacements of each block were equal and the perpendicular lengths and displacements were also equal.In this paper, the hypothesis that the coal mass satisfies the characteristics of isotropic and homogeneous is assumed.The reference value of volume was measured at the vacuum (initial volume V0), and the volume increment was calculated based on the reference value with the pressure increasing step-by-step.
Here, according to the hypothesis and illustration above, the volume, as a function of time, can be written as: pa pe At adsorption time t, the swelling of coal as a function of time is: The swelling before ith adsorption is Q(t)i−1 (the swelling at the end of time exposure to Pi−1).At the ith adsorption pressure step, Q(t)i is the swelling at the end of time exposure to Pi.The fractional swelling increment is qi(t) during the coal sample is exposed to Pi [45]: where qi(t) is the volume strain change of the ith pressure step, which ranges from 0 to 1. qi(t) = 0 means that adsorption just recently occurs at the Pi pressure step, and the swelling increment instantaneously changes with the adsorption amount.qi(t) = 1 indicates that the swelling has reached equilibrium at the given gas pressure Pi.In order to facilitate the use of the elastic mechanic theory, formula ( 4) is converted to the form as follows: v ( ) Here, the Δεvi(t) means the fractional volumetric swelling increment from Pi−1 to Pi.
Except for swelling, the given adsorption pressure generates volume compression due to closure of fractures and cleats, and coal matrix solid is also compressed during the loading process.For the sake of convenience in the computation, it was assumed that parallel lengths and displacements of each block were equal and the perpendicular lengths and displacements were also equal.In this paper, the hypothesis that the coal mass satisfies the characteristics of isotropic and homogeneous is assumed.The reference value of volume was measured at the vacuum (initial volume V 0 ), and the volume increment was calculated based on the reference value with the pressure increasing step-by-step.
Here, according to the hypothesis and illustration above, the volume, as a function of time, can be written as: At adsorption time t, the swelling of coal as a function of time is: The swelling before ith adsorption is Q(t) i−1 (the swelling at the end of time exposure to P i−1 ).At the ith adsorption pressure step, Q(t) i is the swelling at the end of time exposure to P i .The fractional swelling increment is q i (t) during the coal sample is exposed to P i [45]: where q i (t) is the volume strain change of the ith pressure step, which ranges from 0 to 1. q i (t) = 0 means that adsorption just recently occurs at the P i pressure step, and the swelling increment instantaneously changes with the adsorption amount.q i (t) = 1 indicates that the swelling has reached equilibrium at the given gas pressure P i .In order to facilitate the use of the elastic mechanic theory, Equation ( 4) is converted to the form as follows: Here, the ∆ε vi (t) means the fractional volumetric swelling increment from P i−1 to P i .

of 15
Except for swelling, the given adsorption pressure generates volume compression due to closure of fractures and cleats, and coal matrix solid is also compressed during the loading process.In this process, the elastic modulus and Poisson's ratio are not fixed values.Goodman [46] suggested that the strain variation by pressure independent action was: Exposure to the given P i pressure step, the fractional strain variation of the sample is ∆ε Pi from P i−1 to P i where ∆P i is the pressure increment from P i−1 to P i .E s is the elastic modulus of coal matrix solid, which is not equivalent to the Young's modulus (E P ) that takes the elasticity of micropores into account.
The relationship between the elastic modulus E s and the Young's modulus E P can be expressed as [47]: ρ s and ρ are the density of the solid phase (skeletal density) and apparent density, respectively.According to Bentz et al. [48], Poisson' ratio change can be expressed in Equation ( 9) after sorption where ν s is the Poisson ratio of the solid frame, ranging from −1 to 0.5, ν is the effective Poisson ratio, and φ is the porosity.In the process of deduction, it was assumed that the pore shape was cylindrical and the pores were randomly distributed.The relationships of E P /E s − ρ/ρ s as well as ν − ν s based on the research data of Bentz et al. [48] are shown in Figure 4.
Minerals 2017, 7, 26 5 of 15 In this process, the elastic modulus and Poisson's ratio are not fixed values.Goodman [46] suggested that the strain variation by pressure independent action was: Exposure to the given Pi pressure step, the fractional strain variation of the sample is ΔεPi from Pi−1 to Pi where ΔPi is the pressure increment from Pi−1 to Pi. Es is the elastic modulus of coal matrix solid, which is not equivalent to the Young's modulus (EP) that takes the elasticity of micropores into account.The relationship between the elastic modulus Es and the Young's modulus EP can be expressed as [47]: ρs and ρ are the density of the solid phase (skeletal density) and apparent density, respectively.According to Bentz et al. [48], Poisson' ratio change can be expressed in Equation ( 9) after sorption where νs is the Poisson ratio of the solid frame, ranging from −1 to 0.5, ν is the effective Poisson ratio, and ϕ is the porosity.In the process of deduction, it was assumed that the pore shape was cylindrical and the pores were randomly distributed.The relationships of EP/Es − ρ/ρs as well as ν − νs based on the research data of Bentz et al. [48] are shown in Figure 4.The fractional volumetric strain increment and the total volumetric strain is Δεi and εi, respectively.Two parts are included, gas adsorption induced matrix swelling strain Δεvi, taken as positive, and adsorption pressure loading induced strain ΔεPi, taken as negative: Equation ( 11) is an expression of the total volumetric strain after equilibrium by pressure increasing step-by-step under unconfined conditions.This expression of adsorption-induced The fractional volumetric strain increment and the total volumetric strain is ∆ε i and ε i , respectively.Two parts are included, gas adsorption induced matrix swelling strain ∆ε vi , taken as positive, and adsorption pressure loading induced strain ∆ε Pi , taken as negative: Minerals 2017, 7, 26 6 of 15 Equation ( 11) is an expression of the total volumetric strain after equilibrium by pressure increasing step-by-step under unconfined conditions.This expression of adsorption-induced volumetric swelling strain considers the initial volume V 0 and the volume at the time of t.In addition, the changes of Young's modulus and Poisson ratio are also contained.
Under the hypothesis of isotropic and homogeneous, based on the relationship between linear strain and volumetric strain in elastic mechanics, and the linear swelling strain in the vertical direction can be expressed as follows:

Adsorption-Induced Stress under Uniaxial Conditions
After the mining work finished, only two parts of pillars and goafs are left in the original work face (Figure 5).Due to the overburden pressure and stress redistribution by mining activities, the pillar bears the pressure of overlying strata as well as mining stress and preforms uniaxial condition (Figure 6).If the CO 2 gas (fluid) is injected in goaf, this adsorption behavior will occur in the pillar.Under the uniaxial condition, adsorption-induced swelling strain occurs on the two sides adjacent to goafs and drifts, but the swelling strain of the direction perpendicular to bedding is inhibited, resulting in the swelling stress occurring in the vertical direction (z axial, Figure 3).The swelling stress can decrease crack initiation stress, resulting in the damage of the coal pillar.In this section, we consider the derivation of swelling stress in the direction perpendicular to bedding.volumetric swelling strain considers the initial volume V0 and the volume at the time of t.In addition, the changes of Young's modulus and Poisson ratio are also contained.Under the hypothesis of isotropic and homogeneous, based on the relationship between linear strain and volumetric strain in elastic mechanics, and the linear swelling strain in the vertical direction can be expressed as follows:

Adsorption-Induced Stress under Uniaxial Conditions
After the mining work finished, only two parts of pillars and goafs are left in the original work face (Figure 5).Due to the overburden pressure and stress redistribution by mining activities, the pillar bears the pressure of overlying strata as well as mining stress and preforms uniaxial condition (Figure 6).If the CO2 gas (fluid) is injected in goaf, this adsorption behavior will occur in the pillar.Under the uniaxial condition, adsorption-induced swelling strain occurs on the two sides adjacent to goafs and drifts, but the swelling strain of the direction perpendicular to bedding is inhibited, resulting in the swelling stress occurring in the vertical direction (z axial, Figure 3).The swelling stress can decrease crack initiation stress, resulting in the damage of the coal pillar.In this section, we consider the derivation of swelling stress in the direction perpendicular to bedding.Based on the analysis of the section above, under uniaxial conditions, the constitutive relation of swelling stress and swelling strain in the vertical direction (z axial) follows Equation ( 13): Combining Equations ( 12)-( 14), stress increment in the vertical direction in the process of exposure from Pi−1 to Pi can be calculated by Formula (15).The swelling stress includes two parts.The first part is the adsorption-induced swelling stress, and the second is volume-compressed stress: volumetric swelling strain considers the initial volume V0 and the volume at the time of t.In addition, the changes of Young's modulus and Poisson ratio are also contained.Under the hypothesis of isotropic and homogeneous, based on the relationship between linear strain and volumetric strain in elastic mechanics, and the linear swelling strain in the vertical direction can be expressed as follows:

Adsorption-Induced Stress under Uniaxial Conditions
After the mining work finished, only two parts of pillars and goafs are left in the original work face (Figure 5).Due to the overburden pressure and stress redistribution by mining activities, the pillar bears the pressure of overlying strata as well as mining stress and preforms uniaxial condition (Figure 6).If the CO2 gas (fluid) is injected in goaf, this adsorption behavior will occur in the pillar.Under the uniaxial condition, adsorption-induced swelling strain occurs on the two sides adjacent to goafs and drifts, but the swelling strain of the direction perpendicular to bedding is inhibited, resulting in the swelling stress occurring in the vertical direction (z axial, Figure 3).The swelling stress can decrease crack initiation stress, resulting in the damage of the coal pillar.In this section, we consider the derivation of swelling stress in the direction perpendicular to bedding.Based on the analysis of the section above, under uniaxial conditions, the constitutive relation of swelling stress and swelling strain in the vertical direction (z axial) follows Equation ( 13): Combining Equations ( 12)-( 14), stress increment in the vertical direction in the process of exposure from Pi−1 to Pi can be calculated by Formula (15).The swelling stress includes two parts.The first part is the adsorption-induced swelling stress, and the second is volume-compressed stress: Based on the analysis of the section above, under uniaxial conditions, the constitutive relation of swelling stress and swelling strain in the vertical direction (z axial) follows Equation ( 13): Combining Equations ( 12)-( 14), stress increment in the vertical direction in the process of exposure from P i−1 to P i can be calculated by Formula (15).The swelling stress includes two parts.The first part is the adsorption-induced swelling stress, and the second is volume-compressed stress: If, at the given nth pressure step, P n , the pillar perfectly reaches adsorption equilibrium, and the fractional swelling variation does not change and the linear strain also remains constant.In this case, the swelling stress is the cumulative value of stress increments at overall adsorption pressure steps: where σ iz (t) is the total swelling stress considering the adsorption occurring in the pillar in the vertical direction where swelling strain is inhibited under the uniaxial condition.The swelling stress is a cumulative value by step-by-step pressures loading from 1st to nth steps.

Mechanical Strength Change after Gas Adsorption
The pillar coal matrix swells after CO 2 injection in abandoned coal mines, while it shrinks when gas is expelled from coal mass.Adsorption expansion will inevitably lead to some changes in the mechanical properties of the pillar, such as the above mentioned parameters, elastic modulus, Poisson's ratio, bulk modulus as well as shear modulus, due to these parameters being highly effected by adsorption-induced fracturing.All of these variations result in the decrease of pillar strength during CO 2 sequestration in abandoned coal mines.Hu et al. [49] demonstrated that the adsorption of gas exhibits dual effects on the physical properties of coal, mechanical (swelling stress) and non-mechanical (erosion effect) effects.The swelling effect leads to the decrease of the interaction between coal particles, while the erosion effect reduces the surface energy and lowers the surface tension of the coal mass.Gas adsorption mainly affects the cohesion, and the decrease of cohesion leads to obvious plastic deformation of coal, but, for the internal friction angle, it decreases to a certain extent and the hold over no longer changes.Hol et al. [40] verified that gas sorption can lead to a decrease in bulk modulus, while an increase in swelling caused the strain hysteresis to be oversized during the process of loading-unloading.Hagin and Zoback [50] compared the adsorption characteristics of CO 2 with that of helium, and found that the Young's modulus decreased after the CO 2 saturation adsorption.Simultaneously, the static bulk modulus reduced by an order of magnitude.Yang and Zoback [51] observed that CO 2 injection into coal samples resulted in the volume increase, and the coal sample became more viscous and less elastic.Ranjith and Perera [41] considered the effects of cleat density and direction on CO 2 adsorption-induced strength reduction.
On the basis of Griffith line elastic fracture theory, when cracks propagate, a portion of the elastic energy is converted into a solid surface energy.The fracture is confined in the extended critical state when the release rate of the elastic energy is equal to the increase rate of the surface energy.In the plane stress state, the critical stress is: where γ is the surface energy per unit area, N/m; π is the surface tension change value from the vacuum to the given adsorption condition, N/m.In the case of the plane strain state, E is replaced by E − (1 − υ 2 ).Furthermore, the relationship between the solid linear expansion strain and the change in surface tension is as follows: where γ 0 and γ s are the vacuum state and surface tension after gas adsorption, respectively, and λ can be described as Equation ( 19) [52]: where S is the specific surface area of the coal sample, and K s is the apparent modulus.Equation ( 18) is rewritten as: By Equation ( 17), the critical propagation stress at the end of the given pressure P n for the nth adsorption step is expressed as function of time and swelling stress: It can be written as: The square of the critical strength in vacuum conditions is σ 2 c0 , as the following At the end of exposure to the nth adsorption pressure step, the strength reduction rate of the pillar is the ratio of σ 2 ci (t) to σ 2 c0 That is, Equation ( 26) is the calculation formula of the strength reduction rate considering the swelling stress increment under the condition of gas adsorption at partial pressure loading step-by-step.The formula can be used to calculate the strength reduction value of gas adsorption conveniently.The curve of ratio σ 2 ci (t) to σ 2 c0 is showed in Figure 7, and the strength reduction is nearly 17% after adsorption equilibrated based on the swelling strain data at pressure 2 MPa from [45].Compared with other formulas of linear pressure adsorption models, Equation ( 14) is a progressive accumulation form that could be applied to calculate the swelling stress of nonlinear adsorption pressure loading-unloading.Moreover, in the case of re-adsorption after the coal sample having already swelled to a certain degree, Equations ( 16) and ( 26) are also helpful for determining the swelling stress fractional increment and strength reduction value.

Swelling Strain and Swelling Stress
In this paper, we consider the volumetric swelling strain εi(t) and isotropic linear strain (εiz(t)) in the z axial direction (Figure 3) as function of time after equilibrium by adsorption pressure increasing step-by-step under unconfined conditions.In connection with swelling strain, Hol et al. [38] and Liu et al. [39] considered thermodynamic models of gas adsorption and studied the effect of stress on the adsorption concentration of gas as well as sorption behavior.Their starting point is different from this paper based on directly volumetric changes.In the paper of Liu et al. [39], they established the relationship of internal energy, chemical potential, entropy change as well as stressstrain work on a single molecule of gas absorbed by the coal matrix cube.The strain was divided into mean extensional strain and deviatoric strain.In the deviation process of volumetric strain in this paper, it is easy to measure the strain increment without considering the thermodynamic process of adsorption.However, the coupling based on thermodynamic between stress-strain-sorption is significant important to understand the effect of pressure and temperature on adsorption in coal matrix.The swelling strain was divided into a reversible part and irreversible part under unconfined conditions [53].Similarly, Wang et al. [54] divided strain into two parts at an isothermal condition.One is the mechanical deformation meeting the Hooke law stress-strain relationship and is calculated by the effective stress.The other is the deformation induced by gas adsorption or desorption.In Section 2, we do not mention it because not only an elastic strain but also an irreversible strain are converted to swelling stress considered from a macro perspective in the vertical direction under uniaxial loading.
In addition, based on the energy conservation law, and from the viewpoint of surface tension change, Wu et al. [55] and Bai et al. [56] utilized the principle that the expansion of the strain energy equals the surface tension work in order to establish the constitutive equation of swelling stress and swelling strain.Some researchers held the view that gas adsorption satisfies the Langmuir isothermal adsorption model (Figure 8), and considers the linear relationship between adsorption capacity and swelling strain.Although a large number of studies have shown that gas adsorption presents a single molecule arrangement structure that meets the Langmuir model, Lin [57] and Yu et al. [58] suggested that CO2 adsorption on the coal mass can be described in a multi adsorbed layer model, i.e., the Brunauer-Emmett-Teller (BET) adsorption type model.

Swelling Strain and Swelling Stress
In this paper, we consider the volumetric swelling strain ε i (t) and isotropic linear strain (ε iz (t)) in the z axial direction (Figure 3) as function of time after equilibrium by adsorption pressure increasing step-by-step under unconfined conditions.In connection with swelling strain, Hol et al. [38] and Liu et al. [39] considered thermodynamic models of gas adsorption and studied the effect of stress on the adsorption concentration of gas as well as sorption behavior.Their starting point is different from this paper based on directly volumetric changes.In the paper of Liu et al. [39], they established the relationship of internal energy, chemical potential, entropy change as well as stress-strain work on a single molecule of gas absorbed by the coal matrix cube.The strain was divided into mean extensional strain and deviatoric strain.In the deviation process of volumetric strain in this paper, it is easy to measure the strain increment without considering the thermodynamic process of adsorption.However, the coupling based on thermodynamic between stress-strain-sorption is significant important to understand the effect of pressure and temperature on adsorption in coal matrix.The swelling strain was divided into a reversible part and irreversible part under unconfined conditions [53].Similarly, Wang et al. [54] divided strain into two parts at an isothermal condition.One is the mechanical deformation meeting the Hooke law stress-strain relationship and is calculated by the effective stress.The other is the deformation induced by gas adsorption or desorption.In Section 2, we do not mention it because not only an elastic strain but also an irreversible strain are converted to swelling stress considered from a macro perspective in the vertical direction under uniaxial loading.
In addition, based on the energy conservation law, and from the viewpoint of surface tension change, Wu et al. [55] and Bai et al. [56] utilized the principle that the expansion of the strain energy equals the surface tension work in order to establish the constitutive equation of swelling stress and swelling strain.Some researchers held the view that gas adsorption satisfies the Langmuir isothermal adsorption model (Figure 8), and considers the linear relationship between adsorption capacity and swelling strain.Although a large number of studies have shown that gas adsorption presents a single molecule arrangement structure that meets the Langmuir model, Lin [57] and Yu et al. [58] suggested that CO 2 adsorption on the coal mass can be described in a multi adsorbed layer model, i.e., the Brunauer-Emmett-Teller (BET) adsorption type model.
Langmuir isotherm adsorption model equation is given by: where q is the adsorbed amount of gas during adsorption reaching equilibrium, a is the gas limit adsorption capacity at the reference pressure, b is the adsorption equilibrium constant, and p p is the pore pressure.Langmuir isotherm adsorption model equation is given by: where q is the adsorbed amount of gas during adsorption reaching equilibrium, a is the gas limit adsorption capacity at the reference pressure, b is the adsorption equilibrium constant, and pp is the pore pressure.
In the study of the acid gas CO2 and H2S storage in coal seam, Chikatamarla et al. [59] also calculated the frictional expansion strain by the Langmuir model: where V, P, VL, PL are gas adsorption volume, adsorption pressure, Langmuir volume and Langmuir pressure, respectively.In addition, εV and εL are, respectively, the strain at given pressure and strain at infinite pressure, and pe is the Langmuir pressure constant, which is equal to the pressure value when the strain at the pressure is half of the Langmuir maximum strain.For the cylindrical specimen, according to the isotropic assumption of swelling/collapse in the process of adsorption/desorption, the volume strain is expressed as: where V and ΔV are the initial volume and the volume change, respectively.εa and εr are, respectively, the axial strain and the radial strain.
Robertson and Christiansen [60] introduced the strain factor function to modify the Langmuir strain constant.The strain factor function is given by ob where pob denotes the overburden pressure, a and b are empirical constants, Vr denotes the vitrinite reflectance of coal, and γ denotes the bulk density of the gas.p, εL and PL are consistent with Equations ( 28) and (29).
The BET adsorption model is a generalization of the Langmuir model considering multi adsorbed molecule layers, which is expressed by [57,61] where P0 is the saturation pressure of gas, P is the adsorption pressure, Vm is the maximum adsorption volume at the time the entire adsorbent surface is covered with a complete single In the study of the acid gas CO 2 and H 2 S storage in coal seam, Chikatamarla et al. [59] also calculated the frictional expansion strain by the Langmuir model: where V, P, V L , P L are gas adsorption volume, adsorption pressure, Langmuir volume and Langmuir pressure, respectively.In addition, ε V and ε L are, respectively, the strain at given pressure and strain at infinite pressure, and p e is the Langmuir pressure constant, which is equal to the pressure value when the strain at the pressure is half of the Langmuir maximum strain.For the cylindrical specimen, according to the isotropic assumption of swelling/collapse in the process of adsorption/desorption, the volume strain is expressed as: where V and ∆V are the initial volume and the volume change, respectively.ε a and ε r are, respectively, the axial strain and the radial strain.Robertson and Christiansen [60] introduced the strain factor function to modify the Langmuir strain constant.The strain factor function is given by where p ob denotes the overburden pressure, a and b are empirical constants, V r denotes the vitrinite reflectance of coal, and γ denotes the bulk density of the gas.p, ε L and P L are consistent with Equations ( 28) and (29).
The BET adsorption model is a generalization of the Langmuir model considering multi adsorbed molecule layers, which is expressed by [57,61] where P 0 is the saturation pressure of gas, P is the adsorption pressure, V m is the maximum adsorption volume at the time the entire adsorbent surface is covered with a complete single molecular layer, and C is a constant related to the net heat of adsorption.On the basis of the isotropic assumption and the Langmuir adsorption model, the above mentioned literature obtained the stress-strain relationship of swelling.All of those models follow a certain degree of representation in accordance with adsorption induced-swelling.However, they are not perfect because the Langmuir model has a very high accuracy and precision in the case of describing low-pressure gas adsorption.However, it is not suitable for high temperature gas adsorption, particularly, when CO 2 is sequestrated in deep unmined coal seams, where it is confined in states of high temperature and high pressure.Under the situation that CO 2 may exhibit a supercritical state that includes dual characteristics of gas and liquid, and the density of supercritical CO 2 fluid is close to that of the liquid, but the viscosity of which is similar to that of the gas.The diffusion coefficient of supercritical CO 2 is nearly one hundred times of that of the liquid.
Although this paper is also based on the isotropic hypothesis to establish the stress-strain relationship as a function of time and each adsorption pressure step, we do not consider the Langmuir isothermal adsorption model as the reference model but directly analyze it based on volume increment.Furthermore, other models only consider the initial state and the final equilibrium state and have exclusively been expressed in single integral steps, whilst the swelling increment of the prior swelled sample under re-loading is not taken into account.

Stress and Strength Reduction
Coupling relationship of swelling strains, swelling stressed and total gas uptake are affected by coal properties.Fractured coal possesses dual porosity system: (a) the cleat macroporosity system, and (b) the microporosity of coal matrix.Espinoza et al. [62,63] proposed a double porosity poromechanical model and considered the strains as a function of stresses, fracture pore pressure, and the pressure-dependent adsorption stress developed by the coal matrix.
Liu et al. [64] considered the seepage model of fracture-matrix interaction during coal deformation, and put forward the concept of internal swelling stress, σ I (Figure 9).The strain is considered including two parts: the first part is the coal matrix strain induced by the internal swelling stress, and the second part is the volume strain of the fracture.Their work modified the form of effective stress, so that the corrected expression of effective stress can directly reflect its impact on the permeability.Nevertheless, the elastic modulus, Poisson's ratio as well as other parameters are variable rather than fixed values during the occurrence of coal matrix adsorption swelling.These variates are not taken into account in the derivation of theory and formula, which leads to some errors in the accuracy of the model.
On the basis of the isotropic assumption and the Langmuir adsorption model, the above mentioned literature obtained the stress-strain relationship of swelling.All of those models follow a certain degree of representation in accordance with adsorption induced-swelling.However, they are not perfect because the Langmuir model has a very high accuracy and precision in the case of describing low-pressure gas adsorption.However, it is not suitable for high temperature gas adsorption, particularly, when CO2 is sequestrated in deep unmined coal seams, where it is confined in states of high temperature and high pressure.Under the situation that CO2 may exhibit a supercritical state that includes dual characteristics of gas and liquid, and the density of supercritical CO2 fluid is close to that of the liquid, but the viscosity of which is similar to that of the gas.The diffusion coefficient of supercritical CO2 is nearly one hundred times of that of the liquid.
Although this paper is also based on the isotropic hypothesis to establish the stress-strain relationship as a function of time and each adsorption pressure step, we do not consider the Langmuir isothermal adsorption model as the reference model but directly analyze it based on volume increment.Furthermore, other models only consider the initial state and the final equilibrium state and have exclusively been expressed in single integral steps, whilst the swelling increment of the prior swelled sample under re-loading is not taken into account.

Stress and Strength Reduction
Coupling relationship of swelling strains, swelling stressed and total gas uptake are affected by coal properties.Fractured coal possesses dual porosity system: (a) the cleat macroporosity system, and (b) the microporosity of coal matrix.Espinoza et al. [62,63] proposed a double porosity poromechanical model and considered the strains as a function of stresses, fracture pore pressure, and the pressure-dependent adsorption stress developed by the coal matrix.
Liu et al. [64] considered the seepage model of fracture-matrix interaction during coal deformation, and put forward the concept of internal swelling stress, σI (Figure 9).The strain is considered including two parts: the first part is the coal matrix strain induced by the internal swelling stress, and the second part is the volume strain of the fracture.Their work modified the form of effective stress, so that the corrected expression of effective stress can directly reflect its impact on the permeability.Nevertheless, the elastic modulus, Poisson's ratio as well as other parameters are variable rather than fixed values during the occurrence of coal matrix adsorption swelling.These variates are not taken into account in the derivation of theory and formula, which leads to some errors in the accuracy of the model.However, desorption differs from adsorption.Gas desorption can induce coal matrix shrinkage and stress relaxation.Espinoza et al. [65] studied the desorption-induced shear failure under zero-lateral strain condition and simulated the stress path far from the wellbore.Pore pressure in the coal cleats and desorption-induced variations of stresses jointly affect the stress path, i.e., poroelastic coefficient and variation changes of fluid pressure, as well as lateral stress decrease by desorption.The two part change effective stress of coal results in the stress path reaching the failure envelope.Wang et al. [66] considered gas desorption weakens coal by reducing effective stress.Desorption-induced shrinkages can relax lateral stress and aggravate shear failure of coal.Ranjith et al. [32] considered CO 2 adsorption and mechanical behavior of coal, and held that peak strength could decrease when coal is saturated with CO 2 .However, they did not comprehensively interpret the strength reduction mechanism.In this research, pillars are assumed to reach equilibrium on both sides of goaf after the mining working finished, and pillars are under the uniaxial condition.The coal pillar strength reduction is theoretically analyzed, and adsorption can increase the swelling stress on the faces of microfractures and cleats, resulting in effective stress decreasing.Under the condition of uniaxial conditions, if the loading pressure remains constant, and no swelling occurs, cracks could maintain the equilibrium state.However, when adsorption behaviors take place, the crack surface tension changes, and swelling stress emerges by swelling strain being inhibited.The stress needed to cause crack initiation is decreased, leading to coal pillars potentially being more prone to failure.The model proposed in this paper can describe the effect of swelling stress on strength clearly.
Coal is a solid that contains a large number of slit-like pores and fractures inter-connected by narrow capillary constriction and connected to the surface [67], which presents absolutely anisotropy and heterogeneity.Simultaneously, the adsorption induced-swelling also has anisotropic characteristics for the reason that most microfractures developed parallel to the bedding plate, approximately following maceral-maceral and bedding interfaces and swelling anisotropic [53].In this paper, it is assumed that both coal and swelling have isotropic elastic properties for conveniently facilitating the application of existing theories and simplifying the process of calculation.From the perspective of CO 2 sequestration in abandoned coal mines, the strength change of coal pillars must be emphasized to prevent CO 2 leakage from goafs.To accurately predict the swelling stress in the vertical direction and strength reduction of pillars, much more research and many more approaches are still required.

Conclusions
This study has focused on volumetric swelling strain and strength reduction of pillars when CO 2 is stored in abandoned coal mines.The volumetric swelling strain is theoretically derived as a function of time by adsorption pressure increasing step-by-step under unconfined conditions.In connection with the conditions of coal pillars in abandoned coal mines, and a uniaxial loading model is proposed by simplifying the actual condition.Swelling strain in a direction perpendicular to bedding is inhibited when CO 2 adsorption is in the pillar.The effect of adsorption on pillar strength reduction is theoretically analyzed and deduced.Our findings can be summarized as follows: 1.
There are a large number of coal mines that will be closed and some of them are located in deep formations in China.CO 2 storage in abandoned coal mines could be a potential option for greenhouse gas disposal.

2.
The volume strain and swelling stress, as a function of time, and different loading pressure steps are deduced.Equation ( 15) is used to describe swelling stress considering coal has already had prior swelling deformation under the condition of step-by-step non-linear loading and a non-Langmuir isothermal model.The model presented in this paper is different from other models, in which only the initial state and the final equilibrium state are considered, and the incremental swelling process is neglected.

3.
A theoretical model based on linear swelling stress-strain work is proposed to calculate the reduction ratio of coal pillar strength under uniaxial conditions.This theoretical model can be used to describe strength reduction during adsorption under adsorption pressure loading step-by-step.

Figure 2 .
Figure 2. A matchstick model of a coal seam.

Figure 2 .
Figure 2. A matchstick model of a coal seam.

Figure 2 .
Figure 2. A matchstick model of a coal seam.

Figure 3 .
Figure 3. Sketch of original coal sample and swelling.

Figure 3 .
Figure 3. Sketch of original coal sample and swelling.

Figure 5 .Figure 6 .
Figure 5. Schematic diagram of pillar and goaf in the abandoned coal mine.P

Figure 5 .
Figure 5. Schematic diagram of pillar and goaf in the abandoned coal mine.

Figure 5 .Figure 6 .
Figure 5. Schematic diagram of pillar and goaf in the abandoned coal mine.P

Figure 6 .
Figure 6.Schematic diagram of the uniaxial condition of the pillar.