Critical Conditions for Wellbore Failure during CO[sub.2]-ECBM Considering Sorption Stress.
Author(s): Hecheng Xiao [1,2]; Wenda Li (corresponding author) [1,*]; Zaiyong Wang [1,2]; Shuai Yang [3]; Peng Tan [4,5]1. Introduction
Since 95–98% of the total gas can be stored as absorbed gas in the microporous structure [1], deep, unmineable coalbeds with enormous CO[sub.2] storage capacity are considered to be an appealing alternative for carbon dioxide receptacles [2,3]. While large-scale carbon capture and storage (CCS) must evaluate all potential modes of failure, stress changes due to CO[sub.2] injection need to be limited [4]. In addition, significant stress changes caused by sorption-induced swelling also raise the coal failure potential [5], which directly impacts the safety and sustainability of CO[sub.2] -ECBM.
Coal failure during coalbed methane production has received a lot of attention from researchers, whereas coal failure during CO[sub.2] -ECBM has received less attention. Shi and Durucan [6] initially estimated the effective stress change during depletion using an analogy between thermal contraction and gas desorption-induced matrix shrinkage, and experimental data demonstrate excellent agreement [7]. Liu and Harpalani [8] provided an experimental investigation of the horizontal stress variations, and discovered that in addition to the poroelastic effect seen in conventional reservoirs, sorption-induced shrinkage also alters in-situ stress and increases coal failure potential. Additionally, experimental studies revealed that gas with a higher adsorption capacity corresponds to a greater horizontal stress loss in depletion [9,10]. This in turn suggests that CO[sub.2] injection would result in a greater horizontal stress change. The reservoir-scale failure of CO[sub.2] -ECBM was explored by Lu and Connell [11], who noted that stress changes are complex phenomena connected to reservoir pressure, gas mix, and adsorption parameters. Above all, however, these models account for constant vertical stress and a uni-axial strain condition. Due to the non-uniform pore pressure distribution in the reservoir, a numerical method was devised to discover that the common uni-axial strain condition has significant deviations in the coal failure analysis [12]. In order to properly analyze coal failure during CO[sub.2] -ECBM, non-uniform pore pressure distribution, adsorption characteristics, and gas composition must all be taken into account.
Furthermore, shear failure in the reservoir formation is promoted by near-wellbore stress concentration [13]. In addition, local wellbore failure occurred before reservoir-scale coal failure [14]. So, a basis for the safe CO[sub.2] -ECBM or sequestration is wellbore stability rather than coal stability, and it can provide an upper limit for CO[sub.2] storage capacity. Cui et al. [15] developed analytical solutions of the stress field under the uni-aixal strain, general stress, and plane strain conditions, shedding light on the stress distribution in a reservoir with an isotropic horizontal stress state. In the process of CBM depletion, Reisabadi et al. [16] examined the distinct coal failure index of the coal seams under various stress regimes. Masoudian et al. [17] simulated the stress change in coal with respect to the effect of the sorption-induced modification on mechanical behavior. Hu et al. [18] also developed a model for the cement-formation interface failure length during CO[sub.2] -ECBM, demonstrating the importance of the injection rate in determining the interface failure.
Furthermore, pure CO[sub.2] injection causes more sorption stress and permeability decline than N[sub.2] , so injecting a CO[sub.2] /N[sub.2] mixture gas is recommended. Li and Elsworth [19] discovered that a higher proportion of CO[sub.2] in the injected CO[sub.2] /N[sub.2] mixture results in lower shale gas recovery. Wen et al. [20] investigated the process of CO[sub.2] gas replacing CH[sub.4] using low-field NMR technology, and the effects of injection pressure and temperature were analyzed. Moreover, a micro pilot test of underground displacement using N[sub.2] injection demonstrated a significant stimulation effect when compared to a conventional emission test [21]. However, various adsorptive gases cause various sorption stresses, which change the formation and wellbore stability to a different extent. To the best of the authors’ knowledge, wellbore stability of coal seams with anisotropic stress state during mixture gas injection (including pure CO[sub.2] /N[sub.2] ) has not been thoroughly studied.
In this study, a semi-analytical method was applied to investigate the wellbore failure index during reservoir-scale CO[sub.2] -ECBM assuming zero displacement at the outer boundary condition and plane-strain condition near the wellbore. Furthermore, steady reservoir pressure distribution is assumed to simply simulate the mixture gas (or pure CO[sub.2] ) injection process. The analytical total stress field is derived by considering sorption-induced swelling and poroelastic effect, superposing the stress caused by in- situ stress, borehole pressure, reservoir pressure distribution, and gas adsorption. The critical borehole pressure and gas component (CBPGC) is calculated by combining the Mogi–Coulomb failure criterion with the total stress field. The CBPGC can serve as a benchmark for numerical simulation and quick assessment of the CO[sub.2] storage capacity of in- situ coal seams.
2. Model Development
The coal is a porous medium and adsorbent in which the stress state changes as the reservoir pressure and adsorption amount vary during CO[sub.2] -ECBM. In addition, the initial in- situ stress state, including the reservoir pressure, is essential to determine the ultimate stress state and failure potential. In this study, a wellbore with a radius of r[sub.w] covering a cylindrical domain with a radius of r[sub.b] is considered.
2.1. Stress–strain Constitutive Equation
Gases such as N[sub.2] , CH[sub.4] , and CO[sub.2] are absorbed primarily on the internal surface of microporous coal particles, causing significant matrix swelling. Using and extending Cui’s model [15], the volumetric strain associated with mixture gas sorption can be written as follows (1)e[sub.v] =eLbLp/1+bLp, where (2)b[sub.L] =?ci/pei,e[sub.L] =1/bL?eLici/pei,
Additionally, p[sub.ei] and e[sub.Li] are the Langmuir-type swelling constants of gas i, which represent the maximum swelling capacity and the pore pressure at which the volumetric strain is equal to 0.5e[sub.Li] , respectively [22,23]. Additionally, Langmuir-type swelling constants differ from Langmuir-type adsorption constants: V[sub.Li] and P[sub.Li] [9]. Additionally, c[sub.i] is the mole fraction of gas component i in free gas mixture, p is the reservoir (pore) pressure. Moreover, Cui et al. [15] established a relationship by using the experiment data: e[sub.Li] =e[sub.gi] V[sub.Li] ,p[sub.ei] =P[sub.Li] with e[sub.gi] volumetric strain coefficient of gas i, which gave similar values to those of Liu et al. [9]. In addition, e[sub.V0] =eLCH4p0/peCH4/1+p0/peCH4 is the initial volumetric strain caused by CH[sub.4] adsorption in the initial state, and p[sub.0] is the initial reservoir pressure.
The pressure variation during CO[sub.2] -ECBM has two impacts on effective horizontal stresses. In addition to the poroelastic (poromechanical) effect, the sorption-induced matrix shrinkage/swelling also modifies the effective horizontal stresses. Similar to the constitutive law of the non-isothermal poroelastic medium, the constitutive equations of coal seam considering sorption can read [9,15,22,23], (3)s[sub.ij] =E/1+?(e[sub.ij] +?/1-2?e[sub.b] d[sub.ij] )+?pd[sub.ij] +E/3(1-2?)e[sub.V] d[sub.ij] , where E is Young’s modulus of porous medium, ? is Poisson’s ratio, e[sub.b] is the bulk volumetric strain, d[sub.ij] is Kronecker’s delta and i or j is the directional index r,?andz in cylindrical coordinates. Additionally, ? is the Biot constants, which are assumed to be one due to the large amount of cleats and fractures in the coal.
2.2. Analytical Model of Stress Distribution around the Wellbore
The total stresses around the borehole during the CO[sub.2] -ECBM can be split into the two categories below: (1) Mode 1: the stresses induced by the in-situ principal stresses and the bottom-hole pressure, which is significant near the wellbore and referred to the “near wellbore effect”; (2) Mode 2: additional stresses due to fluid flow and gas sorption. Then, the analytical solution of the total stress around the wellbore can be obtained by superposition.
2.2.1. Simplifications and Assumptions
For the gas injection, the reservoir pressure, gas concentration, and stress will change transiently. While as the gas injection continues, the reservoir pressure profile approaches a steady state, which is considered here and allows us to obtain simple solutions and a perceptive understanding. Differing to Cui’s assumptions [15], the anisotropic in-situ stress state and depleted reservoir is considered here. The assumptions conclude:
(1) The tectonic stress affects the coal seams and the initial in-situ stress is anisotropic;
(2) A steady reservoir pressure profile changing logarithmically from a constant borehole pressure at the wellbore (r[sub.w] ) to a reservoir pressure p[sub.res] at the outer boundary (r[sub.b] );
(3) Gas injection can occur after depletion, in which case the reservoir pressure is unequal to the initial value before depletion (p[sub.res] ?p[sub.0] );
(4) A zero displacement at the outer boundary rather than each point (common uni-axial strain model) in the reservoir [15];
(5) A plane strain condition near the wellbore and a constant overburden s[sub.zz] far from the wellbore [14];
(6) At the final steady state, uniform gas composition of the injected gas is achieved.
2.2.2. Mode 1: Stresses Induced by the Initial In-Situ Stress and Borehole Pressure
According to the classical Kirsch solution, the stress distribution around the wellbore for a vertical wellbore is given as [14,24,25] (4)s[sub.rr] [sup.1] =p[sub.w] rw2/r2+sh0+sH0/2(1-rw2/r2)+sh0-sH0/2(1-4rw2/r2+3rw4/r4)cos2?,s[sub.??] [sup.1] =-p[sub.w] rw2/r2+sh0+sH0/2(1+rw2/r2)-sh0-sH0/2(1+3rw4/r4)cos2?,s[sub.zz] [sup.1] =s[sub.v0] -2?(s[sub.h0] -s[sub.H0] )rw2/r2cos2?,s[sub.r?] [sup.1] =-sh0-sH0/2(1+2rw2/r2-3rw4/r4)sin2?,s[sub.rz] [sup.1] =s[sub.?z] [sup.1] =0 where p[sub.w] is the borehole pressure, ? is the angular position around the wellbore. s[sub.H0] ,s[sub.h0] ,s[sub.v0] are the initial maximum, minimum horizontal and vertical stresses, which are determined by stress measurement in the initial in-situ state before depletion/depressurization. Additionally, the superscript “1, 2” denotes the stresses induced by the Mode 1 and 2, respectively.
2.2.3. Mode 2: Stresses Induced by Poroelastic Response and Gas Adsorption
When mixture gas (CO[sub.2] /N[sub.2] ) is gradually injected into the coal seams, the pore pressure eventually approaches a steady state, from which we can derive an analytical solution. Assuming that the wellbore has a radius r[sub.w] covering a cylindrical domain with a radius r[sub.b] and a constant pressure, p[sub.res] is fixed before injection, then the incremental steady reservoir pressure distribution relative to the initial reservoir pressure can be approximated as [15,26] (5)?p=p[sub.res] -p[sub.0] +pw-pi/ln(rw/rb)lnr/rb, where p[sub.w] is borehole pressure, p[sub.res] is the uniform reservoir pressure before mixture gas injection in the coal seams. It should be noted that appropriate injection timing can continue to the time when coal fails in depletion, so p[sub.res] =p[sub.0] can be assumed. The incremental volumetric strain due to displacement of CH[sub.4] with CO[sub.2] /N[sub.2] can be described by the extended Langmuir model as (6)?e[sub.V] =eLbLp/1+bLp-e[sub.V0] ,
Then, the isotropic elastic constitutive law in the incremental form of Equation (3) in the cylindrical coordinates reads (7)?s[sub.rr] =E/1+?(e[sub.rr] +?/1-2?e[sub.b] )+?p+E/31-2??e[sub.V] ,?s[sub.??] =E/1+?(e[sub.??] +?/1-2?e[sub.b] )+?p+E/31-2??e[sub.V] ,?s[sub.zz] =0=E/1+?(e[sub.zz] +?/1-2?e[sub.b] )+?p+E/31-2??e[sub.V] ,
With (8)e[sub.rr] =?u/?r,e[sub.??] =u/r,e[sub.zz] =?uz/?z,e[sub.b] =e[sub.rr] +e[sub.??] +e[sub.zz] ,
The stress equilibrium equation in Mode 2 is given as (9)?srr2/?r+srr2-s??2/r=0,
Combining Equations (7) and (8) with Equation (9), we can obtain the incremental radial displacement (10)?/?r[1/r?ru/?r]=-1+?1-2?/E??p/?r-1+?/3??eV/?r,
Then, integrating Equation (10) two times yields the solution (11)u=-1+?1-2?/EFp/r-1+?/3Fe/r+rC1/2+C2/r,
Then, induced radial, hoop stress due to fluid flow and gas sorption in the incremental form of Mode 2 is derived as follows (12)s[sub.rr] [sup.2] =(1-2?)/r2F[sub.p] +E/3r2F[sub.e] +EC1/2(1-?)-EC2/(1+?)r2, (13)s[sub.??] [sup.2] =-(1-2?)/r2F[sub.p] -E/3r2F[sub.e] +(1-2?)?p+E/3?e[sub.?] +EC1/2(1-?)+EC2/(1+?)r2,
With (14)F[sub.p] =?[sub.rw] [sup.r] ?prdr=r2/4[2(p[sub.i] -p[sub.0] )+q(2ln(r/rb)-1)], (15)F[sub.e] =?[sub.rw] [sup.r] ?e[sub.?] rdr=eL/bL[(b[sub.L] -e[sub.?0] bL/eL)r2/2-1/qr[sub.b] [sup.2] e[sup.-[2(1+bLpres)]/bLq] Ei(2(1+bLpres/bLq+ln(r/rb)))], where Ei(x)=-?[sub.-x] [sup.8] e[sup.-t] /tdt,q=(p[sub.w] -p[sub.res] )/ln(r[sub.w] /r[sub.b] ). Additionally, C[sub.1] ,C[sub.2] can be solved by the imposed boundary conditions. In addition, the term F[sub.p] ,F[sub.e] in Equations (14) and (15) denotes the poroelastic and sorption effect on the radial and hoop stresses in Equations (12) and (13), respectively. When a stress boundary condition s[sub.rr] [sup.2] (r=r[sub.w] )=0 at the borehole wall and zero displacement condition (u[sub.r] (r=r[sub.b] )=0) at the outer boundary (ZDBC) is specified, we can obtain from Equations (11) and (12) (16){C[sub.1] =-2(1-?)/(1+?)rw2+(1-?)rb2[(1+?)(1-2?)/EF[sub.p] (r[sub.b] )+(1+?)/3F[sub.e] (r[sub.b] )]C[sub.2] =(1+?)rw2/2(1-?)C[sub.1]
When a constant stress boundary condition (s[sub.rr] [sup.2] (r=r[sub.w] ,r[sub.b] )=0) at the borehole wall and outer boundary (ZSBC) is specified, i.e., no lateral constrain at the outer boundary, from Equations (11) and (12) we can obtain (17){C[sub.1] =-2(1-?)/rb2-rw2E[(1-2?)F[sub.p] (r[sub.b] )+E/3F[sub.e] (r[sub.b] )]C[sub.2] =(1+?)rw2/2(1-?)C[sub.1]
The incremental stress equations derived here are similar to those developed by Cui et al. [15] in absolute form of total stress. However, Cui’s model cannot study the wellbore stability of the coal seams with anisotropic in-situ stress state. In addition, CO[sub.2] is often injected into the depleted coal seams, i.e., the reservoir pressure before injection is not equal to the initial reservoir pressure (p[sub.0] ?p[sub.res] ). Additionally, the present model can deal with the situation.
2.2.4. The Total Stresses around the Wellbore with Anisotropic In-Situ Stress
Ultimately, the complete solution of stress field around a vertical wellbore during the CO[sub.2] -ECBM process can be obtained by superposition (18)s[sub.rr] [sup.] =s[sub.rr] [sup.1] +s[sub.rr] [sup.2] ,s[sub.??] [sup.] =s[sub.??] [sup.1] +s[sub.??] [sup.2] ,s[sub.r?] [sup.] =s[sub.r?] [sup.1] ,s[sub.zz] [sup.] =s[sub.zz] [sup.1] ,
Under the linear elastic theory, the maximum stress concentration and then failure index usually appears at the wellbore wall [27,28]. So stress distribution and the following failure index can be studied at the borehole wall. Then, the total stresses at the wellbore wall can be given as (19)s[sub.rr] [sup.w] =p[sub.w] , (20)s[sub.??] [sup.w] =s[sub.h0] +s[sub.H0] -2(s[sub.h0] -s[sub.H0] )cos2?-p[sub.w] +(1-2?)(p[sub.w] -p[sub.0] )+E/3(eLbLpw/1+bLpw-eLCH4p0/peCH4/1+p0/peCH4)+EC1/2(1-?)+EC2/(1+?)r2, (21)s[sub.zz] [sup.w] =s[sub.?0] -2?(s[sub.h0] -s[sub.H0] )cos2?, (22)s[sub.r?] [sup.w] =0,
2.3. Failure Criterion
Additionally, an appropriate shear failure criterion for wellbore stability should be applied for coal failure evaluation. The Mogi–Coulomb failure criterion, considering the effect of intermediate principal stress, can match well with field data [14,16,29]. In this study, the Mogi–Coulomb failure criterion is used to forecast coal failure during CO[sub.2] -ECBM [14]. (23)F=t[sub.oct] -a-(bS1+S3/2-p), (24)a=22Ccosf/3, (25)b=22Csinf/3, (26)t[sub.oct] =1/3[square root of (S2-S3)[sup.2] +(S2-S1)[sup.2] +(S1-S3)[sup.2] ], where C,f[sub.] are the rock cohesion and friction angle, and S[sub.1] =max{s[sub.rr] [sup.] ,s[sub.??] [sup.] ,s[sub.zz] [sup.] }, S[sub.2] =median{s[sub.rr] [sup.] ,s[sub.??] [sup.] ,s[sub.zz] [sup.] }, S[sub.3] =minimum{s[sub.rr] [sup.] ,s[sub.??] [sup.] ,s[sub.zz] [sup.] }, respectively. It should be noted that the s[sub.rr] [sup.] ,s[sub.??] [sup.] ,s[sub.zz] [sup.] in Equations (19)–(21) are the principal stresses at the borehole, since shear stresses at the borehole are s[sub.r?] [sup.] =s[sub.rz] [sup.] =s[sub.?z] [sup.] =0.
Then, the total stresses are substituted into the failure criterion to obtain the critical borehole pressure and gas component (CBPGC), defining the critical allowable p[sub.w] with specific gas components c[sub.CO2] in terms of wellbore stability during CO[sub.2] -ECBM.
3. Induced Stress during CO[sub.2] -ECBM and Parametric Analysis
Differing to the conventional reservoir, except for the poroelastic response, the swelling during adsorptive gas injection causes tremendous stress change, which is not well studied during CO[sub.2] -ECBM. So the stress change considering the poroelastic and adsorptive effect during mixture gas (CO[sub.2] /N[sub.2] ) injection must be carefully researched. Additionally, the influencing factors’ effect on the induced stress is also studied. Table 1 presents the typical input data in this study.
3.1. Model Verification
This paper aims to use the superposition method to provide an analytical solution considering the anisotropic in-situ stress state in the coal seams. Figure 1 illustrates the great agreement between the model results for the induced effective stress change and radial displacement due to CO[sub.2] injection and Cui’s model [15], verifying the accuracy of the present model. The induced effective horizontal stress change in Figure 1a can be given as ?s[sup.eff] =(s[sub.rr] [sup.2] +s[sub.??] [sup.2] )/2-?p because the induced effective stress change is caused by fluid flow and sorption. Furthermore, the ZDBC (zero displacement at the outer boundary) condition restricts the lateral deformation and results in an order-of-magnitude difference in radial displacement (seen in Figure 1b) compared to the ZSBC (zero stress change at the outer boundary). As a result, a more significant stress change is induced during CO[sub.2] -ECBM.
3.2. The Stress Caused by Poroelastic Response and Gas Adsorption
In contrast to conventional reservoirs, sorption-induced swelling results in a large shift in the reservoir’s stress, excluding the poroelastic effect. Figure 2 displays the variations between induced stress and radial displacement whether sorption is considered or not, emphasizing the significance of doing so. It should be noted that the steady reservoir pressure distribution results in negative (tensile) radial and hoop stresses, when only the poroelastic effect is considered, while compressive radial and hoop stresses are induced, including the sorption effect as shown in Figure 2a. The swelling/volumetric strains induced by the sorption of CO[sub.2] and CH[sub.4] are different, and then the displacement of CH[sub.4] by CO[sub.2] causes tremendous net swelling. However, the zero displacement at the outer boundary condition overwhelmingly restricts the coal to accommodate the net swelling. In contrast, the induced displacement, considering the sorption effect, is one order more than that merely considering the poroelastic effect in Figure 2b. Thus, considering the strong swelling and boundary condition, the stress state in the whole domain is markedly elevated in Figure 2a. Furthermore, the maximum radial and hoop stresses occur at the borehole; thus, reservoir instability can be analyzed using the wellbore stability issue. Additionally, the incremental reservoir pressure distribution ?p in Equation (6), rather than the real reservoir pressure distribution, causes the stresses change.
From Equations (12)–(15), with the exception of the reservoir pressure distribution (p[sub.0] ,p[sub.i] ,p[sub.w] ), the induced stress is clearly controlled by the mechanical properties (E,?), adsorption parameters (P[sub.e] ,e[sub.L] ), and gas component of CO[sub.2] (c[sub.CO2] ). Then, the parametric analysis of corresponding factors’ effect on the induced stress is made as follows.
3.2.1. The Influence of Langmuir-Type Constants (P[sub.e] ,e[sub.L] )
The influence of adsorption parameters on the distribution of the radial, hoop and effective horizontal stress variation s[sub.rr] [sup.2] ,s[sub.??] [sup.2] ,?s[sup.eff] with normalized radial distance in pure CO[sub.2] injection is displayed in Figure 3. Both adsorption parameters P[sub.e] ,e[sub.L] severely affect the stress distribution and subsequent wellbore instability. Additionally, we may deduce that the more e[sub.LCO2] and less P[sub.eCO2] , the higher the radial, hoop, and effective horizontal stresses can be generated. While the result differs from that under common uni-axial strain condition [15], the induced effective horizontal stress under ZDBC condition varies little throughout the radial distance in some circumstances.
3.2.2. The Influence of Young Modulus and Poisson Ratio of Coal
The influence of mechanical properties on the induced radial, hoop and effective horizontal stress distribution is depicted in Figure 4. The maximum values of induced radial, hoop and effective horizontal stresses grow as the Young modulus and Poisson ratio increase. In addition, the Poisson ratio has less substantial effect on the stress distribution than the Young modulus, as shown in Figure 4. Consequently, sorption-induced stress should be given more consideration in coal seams and shale with a high Young modulus and Poisson ratio.
3.2.3. The Influence of Gas Component
The coal seams’ permeability significantly decreases when pure CO[sub.2] is injected, hence mixture gas (CO[sub.2] /N[sub.2] ) injection is recommended, while wellbore stability during mixture gas injection has not been fully researched. Due to the distinctive sorption-induced swelling characteristics of N[sub.2] , CH[sub.4] , and CO[sub.2] , mixture gas injection with varied gas components can produce various stress changes and wellbore failure indexes. Figure 5 depicts the induced radial, hoop, and effective horizontal stresses with varied gas components. Since more volumetric fraction of N[sub.2] with weakest adsorptive property is injected into the coal seams as c[sub.CO2] decreases, maximum values of induced radial, hoop, and effective horizontal stresses become lower. Furthermore, effective horizontal stress near the wellbore becomes negative with c[sub.CO2] =0.2, suggesting that the displacement of CH[sub.4] with mixture gas (c[sub.CO2] =0.2) results in net coal shrinkage in this case.
4. Wellbore Stability Analysis during Mixture Gas Injection
The developed methodology is applied to investigate the total stress field around the wellbore and coal failure with the different borehole pressures and gas components of CO[sub.2] during CO[sub.2] -ECBM.
4.1. The Total Stress and Failure Index around the Borehole
The preceding investigation indicated that sorption-induced swelling has a significant impact on the stress field during mixture gas injection. In this section, the total stress distribution along with varying gas components c[sub.CO2] is examined. Comparisons of Figure 6 revealed that the gas component of CO[sub.2] significantly affects the total radial and hoop stress distribution along 0° and 90° directions, whose difference is caused by initial anisotropic in-situ horizontal stress. Additionally, when the gas component of CO[sub.2] c[sub.CO2] diminishes, a larger proportion of less adsorptive N[sub.2] is absorbed into the reservoir, which lowers radial stress inside the borehole and hoop stress at the borehole. In this instance, sorption-induced strain turns into net shrinkage in comparison to the initial CH[sub.4] adsorption, resulting in a reduction in sorption-induced stresses. Moreover, tensile hoop stress is generated when pure N[sub.2] (c[sub.CO2] =0,c[sub.N2] =1) is injected into the reservoir, which leads to tensile failure near the wellbore. Similar to but different from thermal fracturing caused by temperature change [30], tensile failure is induced by the displacement of the initial CH[sub.4] with weakly adsorptive N[sub.2] . Tensile failure, which was rarely noticed, should be carefully examined in the mixture gas injection with a large proportion of N[sub.2] (including pure N[sub.2] ) into the coal seams.
Next, failure index around the borehole in Equation (23) with various gas components of CO[sub.2] is presented in Figure 7. Given the large gas component of CO[sub.2] (c[sub.CO2] =0.8), wellbore failure due to significant sorption-induced stress (seen in Figure 6) emerges at the borehole wall. The maximum and minimum principal stresses at the borehole wall in this case are the hoop and radial stresses (i.e.,s[sub.??] [sup.] >s[sub.zz] [sup.] >s[sub.rr] [sup.] ), respectively. Furthermore, since the magnitude of isotropic stress induced by sorption is substantially greater than that of initial in-situ stress, the failure index distribution exhibits negligible anisotropy. Additionally, when c[sub.CO2] =0.2, the failure index around the wellbore lies in the range of -3.1~-6.7 MPa, indicating a stable wellbore state under the given conditions. Additionally, the highest failure index F[sub.max] is located at the borehole wall in a 0° direction. In contrast, when c[sub.CO2] drops further (c[sub.CO2] =0.1), wellbore failure occurs at the borehole wall in the 90° direction. In this instance, a large proportion of weakly adsorptive N[sub.2] induces the smallest hoop stress at the borehole in a 90° direction, as shown in Figure 6b, which is the minimum principal stress (i.e., s[sub.zz] [sup.] >s[sub.rr] [sup.] >s[sub.??] [sup.] ). So the wellbore failure index first declines and then climbs as c[sub.CO2] goes from 1 to 0.1. Then, we may deduce that for the provided parameters in this situation, the gas component of CO[sub.2] must fall within a rational scope in order to preserve wellbore stability.
Figure 8 depicts the evolution of the maximum failure index with borehole pressure p[sub.w] under the varying gas components of CO[sub.2] c[sub.CO2] to investigate wellbore stability in the mixture gas injection. It should be mentioned that the maximum borehole pressure p[sub.w] is constrained to the initial minimum horizontal stress s[sub.h0] to prevent the reservoir from hydraulic fracturing. In the absence of sorption, as borehole pressure p[sub.w] rises, the maximum failure index first sightly decreases and subsequently increases, as illustrated by the violet dashed line in Figure 8. The increase in pore pressure causes effective horizontal and vertical stresses to decline, which leads Mohr’s circle of stress to continuously shift to the left and approach the failure envelope, leading to an increase in wellbore failure index. In addition, the evolution of the maximum failure index is irrelevant with gas component c[sub.CO2] , when just the poroelastic effect is taken into account.
However, the maximum failure index’s evolution with mixture gas injection considering sorption-induced swelling differs from that when only the poroelastic effect is considered. Increased borehole pressure causes F[sub.max] to decline when the gas component of CO[sub.2] reduces from 0.8 to 0.2, indicating that the wellbore turns to be stable in the mixture gas injection. Additionally, for the given specific conditions, when c[sub.CO2] is higher than 0.8, mixture gas injection induces wellbore instability with arbitrary borehole pressure (certainly higher than the reservoir pressure before injection p[sub.i] ).
In contrast, F[sub.max] increases as borehole pressure p[sub.w] rises when the gas component of CO[sub.2] c[sub.CO2] continuously decreases form 0.2 to 0.05. As seen in Figure 6b, net shrinkage and correspondingly lower hoop stress are caused by a considerably smaller fraction of CO[sub.2] associated with a higher proportion of N[sub.2] . Thus, the anisotropy of the principal stress state (s[sub.zz] [sup.] >s[sub.rr] [sup.] >s[sub.??] [sup.] ) continues to be magnified, and the related failure index increases, as shown in Figure 8. Furthermore, it should be noted that a large proportion of N[sub.2] (c[sub.CO2] =0,c[sub.N2] =1) causes tensile hoop stress in Figure 6b, which also sets a lower limit of c[sub.CO2] for preventing wellbore tensile failure.
4.2. Critical Borehole Pressure and Gas Component (CBPGC)
A workflow with an iterative loop was developed to obtain the critical borehole pressure and gas component of CO[sub.2] (CBPGC), maintaining wellbore stability during CO[sub.2] -ECBM. CBPGC with pore pressure before injection p[sub.res] = 3 MPa is displayed in Figure 9. The upper limit of p[sub.w] is constrained by prohibiting hydraulic fracturing in the coal seams, which is approximately equal to p[sub.wmax] =s[sub.h0] . Meanwhile, the lower limit of p[sub.w] is constrained by p[sub.w] >p[sub.res] . Additionally, the wellbore failure index F[sub.max] =0 restricts the upper and lower limits of the gas component of CO[sub.2] .
Additionally, wellbore state moves into Failure Zone I in Figure 9, if a mixture gas with a large proportion of CO[sub.2] that exceeds the upper limit is injected. In this instance, as shown in Figure 6b, the displacement of initial CH[sub.4] with strongly adsorptive CO[sub.2] results in considerable compressive hoop stress, and then wellbore failure occurs near the borehole. By contrast, the wellbore state enters Failure Zone II if a mixture gas with a large proportion of N[sub.2] (the corresponding gas component of CO[sub.2] lower than the lower limit) is injected. Wellbore shear or tensile failure occurs when the initial CH[sub.4] is mainly displaced by weakly adsorptive N[sub.2] , which induces lower compressive and sometimes even negative hoop stress. Therefore, the borehole pressure and gas component of CO[sub.2] should be set within the CBPGC to prevent wellbore instability. Finally, the CBPGC can provide a benchmark for in-situ coal seams’ CO[sub.2] storage capacity in terms of wellbore stability.
5. Conclusions and Suggestions
Assuming steady reservoir pressure distribution and ZDBC (zero displacement at the outer boundary) condition, the stress field taking into account the sorption and poroelastic effect is derived in the coal seams with anisotropic in-situ stress state. Then, the critical borehole pressure and gas component (CBPGC) maintaining wellbore stability is obtained by combining with shear and tensile failure criteria. Due to the distinct sorption characteristics of CH[sub.4] , CO[sub.2,] and N[sub.2] , mixture gas injection with variable gas composition results in diverse stress change, as opposed to the case when only the poroelastic effect is considered. The following findings can be obtained from this study:
(1) The stress field is significantly influenced by the boundary condition and sorption-induced swelling characteristics. The ZDBC condition results in larger stress change in comparison to the constant stress condition at the outer boundary. Furthermore, the sorption-induced swelling of pure CO[sub.2] relative to CH[sub.4] induces compressive radial and hoop stresses, whereas tensile radial and hoop stresses are caused when only the poroelastic effect is considered;
(2) With the exception of the reservoir pressure distribution, mechanical properties and adsorption parameters both influence sorption-induced stresses. The larger sorption-induced stresses are caused by an increase in the Young modulus and Poisson ratio, which suggests that this effect should be taken into account more in the coal and shale with the high Young modulus and Poisson ratio. In addition, adsorption parameters also directly influence sorption-induced stress. More considerable sorption-induced stresses are induced by the larger e[sub.LCO2] and smaller P[sub.eCO2] ;
(3) The borehole pressure and gas component of CO[sub.2] should be restrained by the CBPGC to prevent wellbore from shear and tensile failure. Mixture gas injection with a large proportion of CO[sub.2] would result in considerable sorption-induced hoop stress and wellbore shear failure. By contrast, when mixture gas with a small proportion of CO[sub.2] is injected, the displacement of the initial CH[sub.4] with weakly adsorptive N[sub.2] would induce less compressive and even tensile hoop stress. In this instance, wellbore shear or tensile failure occurs, which is rarely noticed. Therefore, in light of wellbore stability, the CBPGC can provide a benchmark for in-situ coal seams’ CO[sub.2] storage capacity;
This paper derives a semi-analytical solution of stress field and failure index assuming the steady reservoir pressure distribution that occurs after the unsteady seepage stage. So the CBPGC can be treated as an ultimate result of mixture gas injection, and gives the upper and lower limits of CO[sub.2] storage capacity bound by wellbore stability. In addition, the injection scheme design should pay more attention to wellbore stability during unsteady seepage stage. Moreover, as cased, cemented wellbore are frequently constructed to maintain wellbore stability/integrity; interface failure and zone isolation should be thoroughly studied in the future.
Author Contributions
Conceptualization, H.X. and W.L.; methodology, W.L.; validation, Z.W. and S.Y.; formal analysis and investigation, Z.W.; resources, S.Y. and P.T.; writing—original draft preparation, H.X.; writing—review and editing, W.L.; supervision, W.L.; funding acquisition, W.L. All authors have read and agreed to the published version of the manuscript.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Data are contained in the article.
Conflicts of Interest
The authors declare no conflict of interest.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Glossary
Nomenclature
i Gas type e ?, e b Volumetric strain due to adsorption and deformation p e i, e L i Langmuir-type swelling constants V L i, P L i Langmuir-type adsorption constants E, ? Young modulus and Poisson ratio d i j Kronecker’s delta ? Biot coefficient p 0, p r e s Initial and depleted reservoir pressure, respectively p w Borehole pressure r, ?, z Directional index in cylindrical coordinates r w, r b Radius of wellbore and outer boundary s H 0, s h 0, s V 0 Initial maximum, minimum horizontal, and vertical stress, respectively s r r 1, s ? ? 1, s z z 1, s r ? 1, s r z 1, s ? z 1 Normal and shear stress components due to in-situ stress and borehole pressure s r r 2, s ? ? 2 Radial and hoop stress due to poroelastic response and gas adsorption, respectively s r r, s ? ?, s r ?, s z z Total stress components, respectively s r r w, s ? ? w, s z z w, s r ? w Total stress components at the wellbore ? p, ? e V Incremental reservoir pressure and volumetric strain, respectively C 1, C 2 Integration constants F p, F e Integration constants related to poroelastic response and sorption effect, respectively F Wellbore failure index C, f Rock cohesion and friction angle S 1, S 2, S 3 Maximum, median, and minimum stress, respectively t o c t Octahedral shear stress a, b Mogi–Coulomb coefficients c CO 2, c N 2 Gas component of CO2 and N2, respectively CBPGC Critical borehole pressure and gas component of CO2 ZDBC Zero displacement condition at outer boundary ZSBC Constant stress condition at outer boundary
References
1. P.N.K. De-Silva; P.G. Ranjith; S.K. Choi A study of methodologies for CO2 storage capacity estimation of coal., 2012, 91,pp. 1-15. DOI: https://doi.org/10.1016/j.fuel.2011.07.010.
2. M.S. Masoudian; D.W. Airey; A. El-Zein A chemo-poro-mechanical model for sequestration of carbon dioxide in coalbeds., 2013, 63,pp. 235-243. DOI: https://doi.org/10.1680/geot.SIP13.P.026.
3. Z. Fang; X.C. Li A preliminary evaluation of carbon dioxide storage capacity in unmineable coalbeds in China., 2014, 9,pp. 109-114. DOI: https://doi.org/10.1007/s11440-013-0228-5.
4. M.D. Zoback; S.M. Gorelick To prevent earthquake triggering, pressure changes due to CO2 injection need to be limited., 2015, 112,p. E4510. DOI: https://doi.org/10.1073/pnas.1508533112. PMID: https://www.ncbi.nlm.nih.gov/pubmed/26240342.
5. I.D. Palmer; Z.A. Moschovidis; J.R. Cameron Coal failure and consequences for coalbed methane wells.,
6. J.Q. Shi; S. Durucan Drawdown induced changes in permeability of coalbeds: A new interpretation of the reservoir response to primary recovery., 2004, 56,pp. 1-16. DOI: https://doi.org/10.1023/B:TIPM.0000018398.19928.5a.
7. J.Q. Shi; S. Durucan Variation in horizontal stress with pore pressure depletion in coal under uniaxial strain conditions: An update on the modelling of laboratory data., 2018, 187,pp. 94-97. DOI: https://doi.org/10.1016/j.coal.2018.01.009.
8. S. Liu; S. Harpalani Evaluation of in situ stress changes with gas depletion of coalbed methane reservoirs., 2014, 119,pp. 6263-6276. DOI: https://doi.org/10.1002/2014JB011228.
9. T. Liu; S. Liu; B. Lin; X. Fu; C. Zhu; W. Yang; Y. Zhao Stress response during in-situ gas depletion and its impact on permeability and stability of CBM reservoir., 2020, 266,p. 117083. DOI: https://doi.org/10.1016/j.fuel.2020.117083.
10. X. Hou; S. Liu; G. Li; Y. Zhu; A. Liu Quantifying and modeling of in situ stress evolutions of coal reservoirs for helium, methane, nitrogen and CO2 depletions., 2021, 54,pp. 3701-3719. DOI: https://doi.org/10.1007/s00603-021-02511-1.
11. M. Lu; L. Connell Coal failure during primary and enhanced coalbed methane production-theory and approximate analyses., 2016, 154,pp. 275-285. DOI: https://doi.org/10.1016/j.coal.2016.01.008.
12. M. Lu; L. Connell Coal seam failure during primary/enhanced gas production: How failure develops in fields., 2019, 221,p. 103432. DOI: https://doi.org/10.1016/j.coal.2020.103432.
13. D.N. Espinoza; J.M. Pereira; M. Vandamme; P. Dangla; S. Vidal-Gilbert Desorption-induced shear failure of coal bed seams during gas depletion., 2015, 137,pp. 142-151. DOI: https://doi.org/10.1016/j.coal.2014.10.016.
14. F.S. Huang; Y.L. Kang; H. Liu; L.J. You; X.C. Li Critical conditions for coal wellbore failure during primary coalbed methane recovery: A case study from the San Juan Basin., 2019, 52,pp. 4083-4099. DOI: https://doi.org/10.1007/s00603-019-01825-5.
15. X. Cui; R.M. Bustin; L. Chikatamarla Adsorption-induced coal swelling and stress: Implications for methane production and acid gas sequestration into coal seams., 2007, 112,p. B10202. DOI: https://doi.org/10.1029/2004JB003482.
16. M.Z. Reisabadi; M. Haghighi; A. Salmachi; M. Sayyafzadeh; A. Khaksar Analytical modelling of coal failure in coal seam gas reservoirs in different stress regimes., 2020, 128,p. 104259. DOI: https://doi.org/10.1016/j.ijrmms.2020.104259.
17. M.S. Masoudian; A. El-Zein; D.W. Airey Modelling stress and strain in coal seams during CO2 injection incorporating the rock–fluid interactions., 2016, 76,pp. 51-60. DOI: https://doi.org/10.1016/j.compgeo.2016.02.010.
18. C. Hu; F. Wang; C. Ai A model of cement-formation interface failure length in supercritical CO2-ECBM and storage injection well considering the coal swelling effect under high pressure., 2019, 9,pp. 2757-2767. DOI: https://doi.org/10.1007/s13202-019-0659-0.
19. Z. Li; D. Elsworth Controls of CO2-N2 gas flood ratios on enhanced shale gas recovery and ultimate CO2 sequestration., 2019, 179,pp. 1037-1045. DOI: https://doi.org/10.1016/j.petrol.2019.04.098.
20. H. Wen; J. Hao; L. Ma; X.H. Zheng Experimental Study on Replacing Coal Seam CH4 with CO2 Gas., 2022, 7,pp. 1395-1403. DOI: https://doi.org/10.1021/acsomega.1c06050.
21. H.R. Gong; K. Wang; G.D. Wang; X. Yang; F. Du Underground coal seam gas displacement by injecting nitrogen: Field test and effect prediction., 2020, 306,p. 121646. DOI: https://doi.org/10.1016/j.fuel.2021.121646.
22. T. Liu; B.Q. Lin; W. Yang; C. Zhai; T. Liu Coal permeability evolution and gas migration under non-equilibrium state., 2017, 118,pp. 393-416. DOI: https://doi.org/10.1007/s11242-017-0862-8.
23. Y. Wu; J.S. Liu; D. Elsworth Evolution of coal permeability: Contribution of heterogenous swelling processes., 2011, 88,pp. 152-162. DOI: https://doi.org/10.1016/j.coal.2011.09.002.
24. W.B. Bradley Failure of inclined boreholes., 1979, 101,pp. 232-239. DOI: https://doi.org/10.1115/1.3446925.
25. E. Fjær; R.M. Holt; P. Horsrud; A.M. Raaen; R. Risnes, 2nd ed. edition; Elsevier: Amsterdam, The Netherlands, 2008,
26. Y. Li; J. Wang; Z. Wang; Z. Pan Variation in Permeability during CO2–CH4 Displacement in Coal Seams. Part 2: Modeling and Simulation., 2020, 5,pp. 18432-18440. DOI: https://doi.org/10.1021/acsomega.0c02293. PMID: https://www.ncbi.nlm.nih.gov/pubmed/32743220.
27. S.K. Al-Shaaibi; A.M. Al-Ajmi; Y. Al-Wahaibi Three dimensional modeling for predicting sand production., 2013, 109,pp. 348-363. DOI: https://doi.org/10.1016/j.petrol.2013.04.015.
28. E. Detournay; H.D. Cheng Poroelastic response of a borehole in a non-hydrostatic stress field., 1988, 25,pp. 171-182. DOI: https://doi.org/10.1016/0148-9062(88)92299-1.
29. A.M. Al-Ajmi; R.W. Zimmerman Relation between the Mogi and the Coulomb failure criteria., 2005, 42,pp. 431-439. DOI: https://doi.org/10.1016/j.ijrmms.2004.11.004.
30. S. Yin Numerical analysis of thermal fracturing in subsurface cold water injection by finite element methods., 2013, 37,p. 15. DOI: https://doi.org/10.1002/nag.2148.
Figures and Table
Figure 1: Profiles of (a) effective horizontal stress change and (b) induced radial displacement during pure CO[sub.2] injection between the present model and results by Cui et al. [15]. ZDBC and ZSBC denote zero displacement and zero stress change at the outer boundary, respectively. [Please download the PDF to view the image]
Figure 2: (a) The induced stress and (b) radial displacement difference with normalized radial distance considering sorption or not. In this case, pure CO[sub.2] is injected by p[sub.w] = 12 MPa with reservoir pressure before injection p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 3: Profiles of (a) induced radial stress, (b) hoop stress and (c) effective horizontal stress distribution with varying adsorption parameters. In this case, pure CO[sub.2] is injected by p[sub.w] = 12 MPa with p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 4: Profiles of induced (a) radial stress, (b) hoop stress and (c) effective horizontal stress distribution with varying Young modulus and Poisson ratio. In this case, pure CO[sub.2] is injected by p[sub.w] = 12 MPa with p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 5: Profiles of induced (a) radial stress, (b) hoop stress and (c) effective horizontal stress distribution with varying gas components. In this case, mixture gas is injected by p[sub.w] = 12 MPa with p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 6: Profiles of (a) total radial stress and (b) total hoop stress with radius with c[sub.CO2] =0.8,0.2,0.1,0 along 0° and 90° direction. In this case, mixture gas is injected by p[sub.w] = 12 MPa with p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 7: Failure index distribution in the vicinity of the borehole with different gas component: (a) c[sub.CO2] =0.8; (b) c[sub.CO2] =0.2; (c) c[sub.CO2] =0.1. In this case, mixture gas is injected by p[sub.w] = 12 MPa with p[sub.res] = 3 MPa and ZDBC condition is imposed. [Please download the PDF to view the image]
Figure 8: Profile of maximum failure index’s evolution with varying borehole pressure and gas component. In this case, p[sub.res] = 3 MPa. [Please download the PDF to view the image]
Figure 9: CBPGC with pore pressure before injection p[sub.res] = 3 MPa. The critical p[sub.w] in the CBPGC is confined by s[sub.h0] and p[sub.w] >p[sub.res] , and the upper and lower limit of c[sub.CO2] in the CBPGC is obtained by wellbore failure index F[sub.max] =0. [Please download the PDF to view the image]
Table 1: Input parameters for simulation [9,15,16].
| Parameters | Variables | Values |
|---|---|---|
Poromechanical parameters | Young modulus E (MPa) | 4.35 |
Poisson ratio ? | 0.3 | |
Cohesion C(MPa) | 6.44 | |
Friction angle f | 45 | |
Biot coefficient ? | 1 | |
Adsorption parameters | Langmuir-type constants e[sub.L] | N[sub.2] : 0.0074; CH[sub.4] : 0.0106; CO[sub.2] : 0.0389 |
Langmuir-type constants P[sub.e] (MPa) | N[sub.2] : 24.71; CH[sub.4] : 6.02; CO[sub.2] : 4.31 | |
In-situ stress state | Vertical stress s[sub.v0] , MPa | 20.68 |
Maximum horizontal stress s[sub.H0] (MPa) | 16.55 | |
Minimum horizontal stress s[sub.h0] (MPa) | 14.48 | |
Initial reservoir pressure p[sub.0] (MPa) | 9.5148 | |
Injection parameters | Reservoir pressure before injection p[sub.res] (MPa) | 3 |
Borehole pressure p[sub.w] (MPa) | >3 | |
Pure CO[sub.2] injection c[sub.CO2] | 1 | |
Mixture gas injection c[sub.CO2] | 0.8, 0.2, 0.1, 0.05 | |
Mixture gas injection c[sub.N2] =1-c[sub.CO2] | 0.2, 0.8, 0.9, 0.95 | |
Wellbore radius r[sub.w] (m) | 0.1 | |
Studied cylindrical domain radius r[sub.b] (m) | 2000 |
Author Affiliation(s):
[1] Key Laboratory of In-Situ Property-Improving Mining of Ministry of Education, Taiyuan University of Technology, Taiyuan 030024, China
[2] College of Mining Engineering, Taiyuan University of Technology, Taiyuan 030024, China
[3] PetroChina Research Institute of Petroleum Exploration & Development, Beijing 100083, China
[4] CNPC Engineering Technology R&D Company Limited, Beijing 102206, China
[5] National Engineering Research Center for Oil & Gas Drilling and Completion Technology, Beijing 102206, China
Author Note(s):
[*] Correspondence: liwenda@tyut.edu.cn
DOI: 10.3390/su15043696
 Printer friendly
 Cite/link
 Email
 Feedback
| |
| Author: | Xiao, Hecheng; Li, Wenda; Wang, Zaiyong; Yang, Shuai; Tan, Peng |
|---|---|
| Publication: | Sustainability |
| Article Type: | Report |
| Geographic Code: | 9CHIN |
| Date: | Feb 1, 2023 |
| Words: | 7414 |
| Previous Article: | The Assessment of the Spatiotemporal Characteristics of the Eco-Environmental Quality in the Chishui River Basin from 2000 to 2020. |
| Next Article: | Influence of Complexity in Low-Fat Food Packaging on Chinese Consumers’ Purchase Intentions. |
| Topics: | |