Next Article in Journal
Effects of Mooring Compliancy on the Mooring Forces, Power Production, and Dynamics of a Floating Wave Activated Body Energy Converter
Next Article in Special Issue
A Thermo-Hydro-Mechanical-Chemical Coupling Model and Its Application in Acid Fracturing Enhanced Coalbed Methane Recovery Simulation
Previous Article in Journal
AC Breakdown Characteristics of c-C4F8/N2 Gas Mixtures in an Extremely Non-Uniform Electric Field
Previous Article in Special Issue
Field Investigation of Hydraulic Fracturing in Coal Seams and Its Enhancement for Methane Extraction in the Southeast Sichuan Basin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influences of Fracturing Fluid Injection on Mechanical Integrity of Cement Sheath under Four Failure Modes

1
School of Petroleum and Natural Gas Engineering, Chongqing University of Science and Technology, Chongqing 401331, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
3
State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences (CAS), Wuhan 430071, China
*
Author to whom correspondence should be addressed.
Energies 2018, 11(12), 3534; https://doi.org/10.3390/en11123534
Submission received: 29 October 2018 / Revised: 30 November 2018 / Accepted: 17 December 2018 / Published: 19 December 2018

Abstract

:
The significant decreased wellbore temperature and increased casing pressure during fracturing fluid injection present a big challenge for the mechanical integrity of cement sheath in fracturing wells. Based on the theories of elastic mechanics, thermodynamics, and a multi-layer composed thick-wall cylinder, this paper proposed a new mechanical model of cement sheath for fracturing wells, coupling pressure, and thermal loads, which consider the failure modes of de-bonding, radial cracking, disking, and shear failure. The radial nonuniform temperature change and the continuous radial stress and radial displacement at two interfaces have been considered. With the proposed model, the radial distributions of failure stress and the corresponding safety factor for cement sheath during fracturing fluid injection have been analyzed and compared under four failure modes. Results show that the decreased wellbore temperature will produce significant tri-axial tensile stress and induce cement failure of de-bonding, radial cracking, and disking. The increased casing pressure will significantly lower the risk of de-bonding but also aggravate radial cracking and shear failure. For integrity protection of cement sheath, increasing the injected fluid temperature, maintaining higher circulation pumping pressures, and adopting cement sheath with a low elasticity modulus have been suggested for fracturing wells.

1. Introduction

For maintaining gas well’s long-term and safe production, mechanical integrity of cement sheath has been given more emphasis in recent years. Field experiences and laboratory studies have both proven that cement sheath is very likely to fail after downhole operations, such as pressure testing, hydraulic fracturing, acidizing, steam injection, and more [1,2,3]. The cement sheath fails mainly because of the loads from variations of wellbore temperature and pressure. For fracturing wells, because of the high displacement and pump pressure fracturing fluid injection into wellbore, the wellbore temperature decrease (by −70 °C) and casing pressure increase (by 70 MPa) are both very serious. In addition, this leads to cement sheath failure and continuous annular casing pressure. Presently, Portland cement is still frequently-used for wellbore purposes. However, because of its high elasticity and lower tensile strength, it is more likely to fail. Consequently, an analysis of the failure stress and the corresponding failure mode in the cement sheath during the fracturing process is significant for the wellbore safety of fracturing wells.
Currently, some theoretical mechanical models of cement sheath coupling pressure and temperature loads have been established for high pressure and high temperature (HPHT) gas wells and steam injecting wells. Thiercelin et al. [4] first proposed the mechanical model for the cement sheath in which the influence of thermo-elastic properties of materials on cement failure was verified. Li et al. [5] proposed a thermal stress model for a coupled casing-cement-formation system in a thermal well and investigated the casing failure behaviors. After that Li et al. [6] further established a mechanical model of cement sheath coupling the effects of temperature and pressure loads in a non-uniform in-situ stress field. Bois et al. [7] also proposed a cement sheath model. The failure modes were analyzed by investigating the casing deformation induced by wellbore temperature variation. Haider et al. [8] developed a composite axisymmetric casing-cement-formation model for a CO2 storage well. The results indicated that the decrease in wellbore temperature could induce radial tensile stress in the cement sheath. Bui et al. [9] also presented a mechanical model for cement failure prediction, where the non-uniform stress field and the thermal stress were involved. Xi et al. [10] established a transient thermo-mechanical model of casing-cement sheath-formation assembly by using the analytical method and the numerical method comprehensively. Our previous study [11] proposed a mechanical model of casing-cement-formation system and investigated the influence of wellhead casing pressure on cement sheath integrity for HPHT gas wells. After that, we further proposed a model of thermal stress for cement sheath during hydraulic fracturing [12], where radial non-uniform temperature change was involved, and the main influencing factors were also investigated. All the models mentioned above could be used to simulate the stress and failure of cement sheath in some ways. However, because of the complexity of wellbore geometry and the diversity of cement sheath failure modes, the present models cannot provide a relatively accurate and comprehensive failure analysis for cement sheath.
In this paper, a new mechanical model of cement sheath, coupling pressure, and thermal loads and considering the failure modes of de-bonding, radial cracking, disking, and shear failure, has been proposed for fracturing wells, based on the elastic mechanics, thermodynamics theory, and the theory of multi-layer composed thick-wall cylinder. The radial non-uniform temperature change in the combined system has been considered and both the radial stress and radial displacement at the casing-cement sheath interface and the cement sheath-formation interface are supposed to be continuous. Based on the proposed model, the radial distributions of failure stress and corresponding safety factor for cement sheath during fracturing fluid injection have been analyzed and compared under four failure modes. In addition, some measures have been proposed for cement sheath protection, which can be of great significance for decreasing the risk of sustained casing pressure after fracturing operations.

2. Mechanical Model Development and Solution

2.1. Basic Assumptions

To maintain the wellbore mechanical and hydraulic integrity during fracturing and long-term production, the open hole is usually cemented with a steel casing. When the high pump pressure fluid with continuous large displacement is injected into the wellbore, the coupled loads from increased casing pressure pi and decreased wellbore temperature T(r) will simultaneously exert on the casing-cement sheath-formation combined system, which induces cement sheath stresses. The schematic diagram of the casing-cement-formation combined system during fracturing fluid injection is shown in Figure 1. Before developing the mechanical model of cement sheath, to make the present model tractable, the following assumptions are made [11,12].
(1)
The casing, cement sheath, and formation are continuous, homogeneous, and isotropic.
(2)
The casing-cement sheath-formation system is completely cemented and axisymmetric.
(3)
Both casing-cement sheath and cement sheath-formation interface satisfy the continuity conditions for radial stress and radial displacement.
(4)
The radial non-uniform temperature varies along the combined casing-cement sheath -formation system except for casing as consideration of its thin wall and well heat conduction performance.
(5)
The temperature at the inner wall of casing is equal to the wellbore fracturing fluid temperature and that at the outer wall of formation maintains an original formation temperature during fracturing.

2.2. Cement Sheath Stress Induced by the Decrease of Wellbore Temperature

The continuous fracturing fluid injection will significantly decrease wellbore temperature and constrict the combined system, which will lead to thermal stresses in the cement sheath. Our previous study investigated thermal stresses for cement sheath during hydraulic fracturing based on the elastic mechanics and thermodynamics theory [13], and the general solutions of radial thermal displacement and thermal stresses for thick wall cylinder were given as follows [12,13].
{ u T = 1 + μ 1 μ α r a r T ( r ) r d r + C 1 r + C 2 r σ r T = α E 1 μ 1 r 2 a r T ( r ) r d r + E 1 + μ ( C 1 1 2 μ C 2 r 2 ) σ θ T = α E 1 μ 1 r 2 a r T ( r ) r d r α E T b ( r ) 1 μ + E 1 + μ ( C 1 1 2 μ + C 2 r 2 ) σ z T = α E T ( r ) 1 μ + 2 μ E C 1 ( 1 + μ ) ( 1 2 μ )
In Equation (1), T(r) stands for the radial temperature change at radius of r in the cylinder after and before fracturing fluid injection. Considering wellbore temperature varies radially constant in casing and logarithmically in cement sheath and formation, if the temperature at the inner boundary of casing is Ti before fracturing and it decreased to Tt after fracturing fluid injection, while the temperature at the outer boundary of formation is constantly Te, the radial wellbore temperature change T(r) can be determined by the following [12].
T ( r ) = { T t T i ( r i r r 1 ) T t T i + ( T i T t ) ln ( r / r 1 ) ln ( r o / r 1 ) ( r 1 r r o )
For the integral expression a r T ( r ) r d r in Equation (1), it can be deduced as:
a r T ( r ) r d r = { 1 2 ( T t T i ) ( r 2 r i 2 ) ( r i r r 1 ) 1 2 ( T t T i ) ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r r 1 ) r 2 2 + r 1 2 2 ] ( r 1 r r o )
Submitting Equations (2) and (3) into Equation (1), the thermal radial displacement and thermal stresses for combined system can be determined, and the detailed derivative results are shown in Appendix A.

2.3. Cement Sheath Stresses Induced by the Increase of Casing Pressure

Taking wellbore temperature as no change during fracturing fluid injection and, since only the increased casing pressure acts on the casing-cement sheath-formation combined system, casing expansion and induced stresses in the cement sheath will occur. Assuming the plane strain of the combined system (εz ≈ 0), the following can be obtained [13].
σ z P = μ ( σ r P + σ θ P )
Furthermore, the radial displacement of a cylinder induced by the increased inner pressure can be expressed below.
u P = r E [ ( 1 μ 2 ) σ θ P ( μ + μ 2 ) σ r P ]
If the contact pressure induced by increased casing pressure is pc1 at the casing-cement sheath interface and pc2 at the cement sheath-formation interface respectively, the induced stress in the cement sheath can be calculated with the Lame Formula [13]. Now, we will solve pc1 and pc2 with the interfacial continuity conditions.
For the casing, because it is subjected to the increased casing pressure pi at its inner wall and the contact pressure pc1 at its outer wall, its radial and tangential stresses can be obtained with the Lame Formula as follows:
{ σ r s P = p i r i 2 r 1 2 r i 2 ( 1 r 1 2 r 2 ) p c 1 r 1 2 r 1 2 r i 2 ( 1 r i 2 r 2 ) σ θ s P = p i r i 2 r 1 2 r i 2 ( 1 + r 1 2 r 2 ) p c 1 r 1 2 r 1 2 r i 2 ( 1 + r i 2 r 2 )
By submitting Equation (6) into Equation (5), the radial displacement of the casing at its outer wall (r = r1) can be obtained.
u s P | r = r 1 = r 1 E s [ ( 1 μ s 2 ) ( 2 r i 2 r 1 2 r i 2 p i r 1 2 + r i 2 r 1 2 r i 2 p c 1 ) + ( μ s + μ s 2 ) p c 1 ]
For the cement sheath, because it is subjected to the contact pressure pc1 at its inner wall and the contact pressure pc2 at its outer wall, its radial and tangential stresses can also be obtained with the Lame Formula below.
{ σ r c P = p c 1 r 1 2 r 2 2 r 1 2 ( 1 r 2 2 r 2 ) p c 2 r 2 2 r 2 2 r 1 2 ( 1 r 1 2 r 2 ) σ θ c P = p c 1 r 1 2 r 2 2 r 1 2 ( 1 + r 2 2 r 2 ) p c 2 r 2 2 r 2 2 r 1 2 ( 1 + r 1 2 r 2 )
By submitting Equation (8) into Equation (5), the radial displacement of the cement sheath at its inner wall (r = r1) can be obtained.
u c P | r = r 1 = r 1 E c [ ( 1 μ c 2 ) ( r 2 2 + r 1 2 r 2 2 r 1 2 p c 1 2 r 2 2 r 2 2 r 1 2 p c 2 ) + ( μ c + μ c 2 ) p c 1 ]
According to the continuity condition at the casing-cement sheath interface, there exists u s P | r = r 1 = u c P | r = r 1 . Thus, we can get the following.
[ E c ( μ s + μ s 2 ) r 1 2 + r i 2 r 1 2 r i 2 E c ( 1 μ s 2 ) r 2 2 + r 1 2 r 2 2 r 1 2 E s ( 1 μ c 2 ) E s ( μ c + μ c 2 ) ] p c 1 + 2 r 2 2 r 2 2 r 1 2 E s ( 1 μ c 2 ) p c 2 = 2 r i 2 r 1 2 r i 2 E c ( 1 μ s 2 ) p i
Similarly, the radial displacement of the cement sheath at its outer wall (r = r2) can also be obtained as follows.
u c P | r = r 2 = r 2 E c [ ( 1 μ c 2 ) ( 2 r 1 2 r 2 2 r 1 2 p c 1 r 2 2 + r 1 2 r 2 2 r 1 2 p c 2 ) + ( μ c + μ c 2 ) p c 2 ]
For the formation, because it is subjected to the contact pressure pc2 at its inner wall and the changed formation pressure pf at its outer wall, its radial and tangential stresses can also be obtained with the Lame Formula below.
{ σ r f P = p c 2 r 2 2 r o 2 r 2 2 ( 1 r o 2 r 2 ) p f r o 2 r o 2 r 2 2 ( 1 r 2 2 r 2 ) σ θ f P = p c 2 r 2 2 r o 2 r 2 2 ( 1 + r o 2 r 2 ) p f r o 2 r o 2 r 2 2 ( 1 + r 2 2 r 2 )
Then, by submitting Equation (12) into Equation (5), the radial displacement of the formation at its inner wall (r = r2) can be obtained.
u f P | r = r 2 = r 2 E f [ ( 1 μ f 2 ) ( r o 2 + r 2 2 r o 2 r 2 2 p c 2 2 r o 2 r o 2 r 2 2 p f ) + p c 2 ( μ f + μ f 2 ) ]
According to the continuity condition at the cement sheath-formation interface, there exists u c P | r = r 2 = u f P | r = r 2 . Thus, we can get the following.
[ E f ( μ c + μ c 2 ) r 2 2 + r 1 2 r 2 2 r 1 2 E f ( 1 μ c 2 ) r o 2 + r 2 2 r o 2 r 2 2 E c ( 1 μ f 2 ) E c ( μ f + μ f 2 ) ] p c 2 + 2 r 1 2 r 2 2 r 1 2 E f ( 1 μ c 2 ) p c 1 = 2 r o 2 r o 2 r 2 2 E c ( 1 μ f 2 ) p f
Solving Equations (10) and (14) simultaneously, the contact pressure pc1 and pc2 can be obtained. Lastly, from Equation (8) we can easily get the radial stress σ r c P and tangential stress σ θ c P respectively. Furthermore, the axial stress σ z c P in the cement sheath caused by the increased casing pressure can be calculated by Equation (4) as:
σ z c P = μ c ( σ r c P + σ θ c P )

2.4. Combined Stresses of Cement Sheath during Fracturing Fluid Injection

In fact, the cement sheath will simultaneously suffer from induced stresses caused by decreased wellbore temperature and increased casing pressure during the fracturing fluid injection. Therefore, the combined radial, tangential, and axial stress of cement sheath during fracturing fluid injection can be expressed as follows.
{ σ r c = σ r c T + σ r c P σ θ c = σ θ c T + σ θ c P σ z c = σ z c T + σ z c P

2.5. Failure Criterion and Safety Factor for Cement Sheath

Four kinds of failure modes may lead to the loss of zonal isolation for cement sheath including inner debonding, outer debonding, radial cracking, shear damage, and disking. Bois et al. [14] illustrated that in Figure 2. Debonding, radial cracking, and disking can be determined with the maximum tensile stress criterion, while shear failure can be determined with the Mohr-Coulomb Criterion or the Mogi-Coulomb Criterion. Compared to the Mogi-Coulomb Criterion, because the Mohr-Coulomb Criterion does not consider the effect of intermediate principal stress on the strength of the material, it cannot reflect the real failure under a true tri-axial stress state in the rock [15,16]. The Mogi-Coulomb Criterion was selected in the paper.
If the radial tensile stress at the casing-cement sheath interface or the cement sheath-formation interface exceeds the tensile strength of cement sheath, the de-bonding failure can be described as follows [17].
σ r c | r = r 1 o r r 2 > σ t c
If the tangential tensile stress in the cement sheath surpasses its tensile strength, cracking will begin to form and propagate along the radial direction. When multiple radial cracking are linked, the leakage path can be generated. The radial cracking failure can be expressed by Equation (18) below.
σ θ c > σ t c
Similarly, if the axial tensile stress in the cement sheath exceeds its tensile strength, disking will begin to form and run across the cement wall. The disking failure can be described below.
σ z c > σ t c
Shear failure is a critical failure mode, which typically causes an entire failure for cement sheath. The Mogi-Coulomb Criterion assumes that shear failure will occur when the octahedral shear stress τoct in the cement sheath exceeds the maximum allowable shear stress τcmax, which can be expressed as the following equation [15,16].
1 3 ( σ r c σ θ c ) 2 + ( σ r c σ z c ) 2 + ( σ θ c σ z c ) 2 = τ o c t > τ c max
If a linear form is considered for the Mogi-Coulomb criterion and, noting that the compressive stress is negative in this paper, the maximum allowable shear stress τcmax in Equation (20) can be determined by the equation below [18].
τ c max = c c cos ϕ c σ 1 + σ 3 2 sin ϕ c
Considering the following relationship between the uniaxial compressive strength of cement sheath σcc and its cohesion cc, the internal friction angle ϕc is shown below.
σ c c = 2 c c cos ϕ c 1 sin ϕ c
If σcc and ϕc are both known, by combining Equations (20)–(22), we can express the Mogi-Coulomb Criterion below.
1 3 ( σ r c σ θ c ) 2 + ( σ r c σ z c ) 2 + ( σ θ c σ z c ) 2 > 1 - sin ϕ c 2 σ c c sin ϕ c σ 1 + σ 3 2
According to the failure criteria above, the safety factors for four failure modes can be calculated respectively, as shown in Table 1.

3. Mechanical Integrity Analysis of Cement Sheath during Fracturing Fluid Injection

3.1. Basic Calculation Parameters

A shale gas fracturing well in the Sichuan Basin in China was selected for the analysis of mechanical integrity of cement sheath during fracturing fluid injection. The 215.9 mm open hole was cemented with 139.7 mm × 9.17 mm P110 casing. The basic calculation parameters, including wellbore geometry parameters, material property parameters, and operation parameters, are listed in Table 2.

3.2. Results and Discussion

To make clear the failure laws of cement shear during fracturing fluid injection, the proposed cement sheath stress model has been validated by the FEA method and then the corresponding failure stresses and safety factors of cement sheath under a different failure mode have been calculated according to the basic parameters in Table 2. Simultaneously, when seeking the measures of protecting cement sheath, the safety factors have also been investigated under a different elasticity modulus of the cement sheath Ec.

3.2.1. Cement Sheath Stresses Validated by the FEA Method

Up to now, as the complexity of the casing-cement sheath-formation combined system and the difficulty of applying coupled loads, fewer experimental studies were conducted about the in situ cement sheath integrity during fracturing. To validate the cement sheath stress model proposed in this paper, a 2D axisymmetric casing-cement sheath-formation model has been established with Ansys Workbench 14.5. Figure 3 shows the finite element analysis model of casing-cement sheath-formation combined system and the element type of Plane 183 was selected for this analysis. For simulating cement sheath stresses induced by increased casing pressure and decreased wellbore temperature during fracturing, a Steady-State Thermal analysis followed by a Static Structural analysis has been conducted. The axial displacement on the upper and lower bounds of the whole model and the radial displacement on the outer boundary of formation are all set to be zero. Furthermore, 70 MPa net increased casing pressure is exerted on the inner boundary of casing. The initial temperature is uniform at 100 °C on the whole model and it decreased to 30 °C on the casing inner wall while still maintaining 100 °C on the formation outer wall during fracturing.
Figure 4 shows the comparison of cement sheath stresses between analytical model results and finite element analysis ones. It can be seen that the radial stress, tangential stress, and axial stress in the cement sheath calculated with two methods are all very close. The average relative errors of three kinds of stress are 7%, 5%, and 2% respectively, which indicated that the proposed model can be applicable for cement sheath stress calculation during fracturing to some extent.

3.2.2. De-Bonding

De-bonding occurs when radial tensile stress at the cement interfaces exceeds its own tensile strength, which may mainly result from decreased wellbore temperature during fracturing fluid injection [19]. Figure 5 shows the radial distribution of the radial stress in the cement sheath considering increased casing pressure pi of 0–70 MPa and decreased temperature at the casing inner wall of 70 °C. It can be seen from Figure 5 that the decreased wellbore temperature will induce significant radial tensile stress, which is very harmful to cement sheath and may cause de-bonding. However, the increased casing pressure will produce radial compressed stress in the cement sheath at the same time and offset the radial tensile stress. In addition, with the increase of casing pressure pi, the radial stress in the cement sheath gradually changes from tensile to compressed stress, which indicates that the risk of de-bonding decreases. When the increased casing pressure is 50 MPa, the radial stress is nearly compressive. Consequently, maintaining higher circulation pumping pressures, avoiding pumping off during injection, and lowering the pump pressure slowly during the operation completion will all be beneficial for cement sheath protection.
Figure 6 shows the safety factors of cement sheath for de-bonding at two interfaces under different increased casing pressures pi. It can be seen form Figure 6 that the safety factors at the casing-cement interface and the cement-formation interface both increase exponentially with the increase of casing pressure pi. When pi is less than 32 MPa, both interfaces will separate. When pi is between 32 MPa and 42 MPa, only the cement-formation interface will separate. When pi is larger than 42 MPa, both interfaces will not separate. Generally, the outer de-bonding is more relatively easy compared to inner de-bonding. This is because increased casing pressure pi has a larger effect on the casing-cement interface. For the integrity of cement sheath, the value of pi must be maintained above 42 MPa during the whole fracturing fluid injection stage.

3.2.3. Radial Cracking

Radial cracking appears when tangential tensile stress in the cement surpasses its own tensile strength, Figure 7 and Figure 8 shows the radial distribution of the tangential stress and safety factors for radial cracking in the cement sheath considering increased casing pressure pi of 0–70 MPa and decreased temperature at the casing inner wall of 70 °C. It can be seen from Figure 7 that the decreased wellbore temperature and increased casing pressure will both induce significant tangential tensile stress, which leads to the risk of radial cracking. Figure 8 shows that the safety factors are all less than one under the present operation parameters, which means that radial cracking occurs completely. To some extent, this answers for the sustaining annular pressure after gas production. Meanwhile, we can also know easily that the maximum tangential tensile stress is nearly located at the casing-cement sheath interface and is caused by a decreased temperature of 70 °C, which is nearly 10 MPa less than that caused by increased casing pressure of 70 MPa. This indicates that radial cracking is mainly induced by increased casing pressure pi and it occurs from the inner wall to the outer wall. However, pi cannot be reduced during fracturing fluid injection because the need of fracture generation. Consequently, some other measures must be taken for protecting cement sheath, such as selecting cement with low elasticity and high tensile strength.

3.2.4. Disking

Disking occurs when axial tensile stress in the cement sheath exceeds its tensile strength when the gas leaking will run across the cement wall. Figure 9 and Figure 10 show, respectively, the radial distribution of the axial stress and safety factors for disking in the cement sheath considering increased casing pressure pi of 0–70 MPa and decreased temperature at the casing inner wall of 70 °C. It can be seen from Figure 9 that the decreased wellbore temperature will induce significant axial tensile stress, which may cause disking and is very threatening for cement sheath. Yet, the increased casing pressure will produce axial compressed stress and offset a portion of the axial tensile stress but is very limited. Figure 10 shows that the safety factors for disking are all less than one under the present operation parameters, which indicates that disking occurs through the cement sheath. This is also the reason for sustaining annular pressure. Taking Figure 9 and Figure 10 into consideration simultaneously, it can be found that disking mainly results from the decreased wellbore temperature and it will generate from the inner wall to the outer wall. Consequently, increasing the temperature of injected fracturing fluid can reduce the risk of disking.

3.2.5. Shear Failure

Shear failure generally causes the entire failure of cement sheath. Figure 11 and Figure 12 show the radial distribution of the octahedral shear stress and safety factors for shear failure in the cement sheath considering increased casing pressure pi of 0–70 MPa and decreased temperature at the casing inner wall of 70 °C. It can be seen from Figure 11 that the decreased wellbore temperature and increased casing pressure induce significant octahedral shear stress. When the decreased temperature at the casing inner wall is 70 °C with the increase of pi, the octahedral shear stress τoct decreases first and then increases. The maximum τoct is located at the casing interface. From Figure 12, we can conclude that the shear failure of cement sheath nearly does not occur if the increased casing pressure pi is less than 50 MPa while the shear failure will initiate from the inner wall to the outer wall. Considering the present operation parameters, some measures must be taken, such as selecting cement with low elasticity, for protecting cement sheath.

3.3. Comparison and Discussion

Figure 13 and Figure 14 show the radial distribution of the failure stress and corresponding strength in the cement sheath considering increased casing pressure pi of 0 or 70 MPa and decreased temperature at the casing inner wall of 70 °C, respectively. It can be seen from Figure 13 that the decreased wellbore temperature has produced significant tri-axial tensile stress in the cement sheath and three of them exceed its own tensile strength, which will induce cement failure of de-bonding, radial cracking, and disking, especially de-bonding and disking. While the octahedral shear stress is less than the maximum allowable shear stress for cement sheath, it shows that shear failure will not appear. Consequently, increasing properly the injected fluid temperature can lower the failure risks of cement sheath for fracturing wells.
Figure 14 shows that, when considering the increased casing pressure pi of 70 MPa simultaneously, the radial stress has turned into compressive but tangential and axial stresses, which are still tensile and still clearly larger than the tensile strength of cement sheath, which will lead to cement failure of radial cracking and disking. Meanwhile, the octahedral shear stress has also exceeded the maximum allowable shear stress near the casing interface, which indicates that shear failure occurs at this location. Combining the analysis above, we know that the increased casing pressure will significantly lower the de-bonding risk but aggravate radial cracking and shear failure. For protecting cement sheath integrity, higher circulation pumping pressures should be maintained in order to avoid de-bonding during whole fracturing operations.
In order to lower the risks of radial cracking, disking, and shear failure for cement sheath under a high pumping pressure from fracturing, the common method is cemented with low elasticity modulus cement sheath. Figure 15, Figure 16 and Figure 17 show, respectively, the radial distribution of the cement sheath safety factors for three kinds of failure modes under a different elasticity modulus of cement sheath Ec. The considered increased casing pressure pi is 60 MPa and the decreased temperature at the casing inner wall is 50 °C. It can be seen from Figure 15, Figure 16 and Figure 17 that, with a decrease of Ec, three kinds of safety factors increase gradually and the increasing range is larger and larger. Therefore, lowering the elasticity modulus of cement sheath is very effective for cement protection. When Ec is about lower than 4 GPa, cement sheath will not occur to radial cracking and disking. Meanwhile, when Ec is about lower than 8 GPa, the shear failure may not occur. In conclusion, cement sheath with a low elasticity modulus must be adopted in fracturing wells.

4. Conclusions

(1)
The radial distributions of failure stresses and corresponding safety factors for cement sheath have been analyzed under four failure modes, including de-bonding, radial cracking, disking, and shear failure.
(2)
The decreased wellbore temperature will produced significant tri-axial tensile stress and induce cement failure of de-bonding, radial cracking, and disking.
(3)
The increased casing pressure will significantly lower the risk of de-bonding but also aggravate radial cracking and shear failure.
(4)
Increasing the injected fluid temperature, maintaining a higher circulation pumping pressures and adopting cement sheath with a low elasticity modulus can protect cement sheath in fracturing wells.
(5)
Considering some assumptions before model development, some conclusions may not be entirely appropriate. Further model improvement including casing centralization, anisotropy of the formation stresses, and a transient wellbore temperature change should be investigated.

Author Contributions

Conceptualization, H.X., B.Y., and T.M. Investigation, H.X. and N.P. Methodology, H.X. and N.P. Writing—original draft, H.X., T.M., and N.P. Writing—review and editing, B.Y. and T.M.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51804060), the Basic and Frontier Research Programs of Chongqing Science & Technology Commission of China (Grant No. cstc2018jcyjAX0614), the Science and Technology Research Program of Chongqing Municipal Education Commission (Grant No. KJ1713321), the Research Foundation of Chongqing University of Science & Technology (Grant No. CK2016B13), and the Program of Introducing Talents of Discipline to Chinese Universities (111 Plan) (Grant No. D18016).

Conflicts of Interest

The authors declared that they do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

List of Symbols

σ r T , σ θ T , σ z T Radial, tangential, and axial thermal stress of cylinder, respectively, MPa
σ r P , σ θ P , σ z P Radial, tangential, and axial stress of cylinder caused by increased inner pressure, respectively, MPa
σ r s T , σ θ s T , σ z s T Radial, tangential, and axial thermal stress of casing, respectively, MPa
σ r c T , σ θ c T , σ z c T Radial, tangential, and axial thermal stress of cement sheath, respectively, MPa
σ r f T , σ θ f T , σ z f T Radial, tangential, and axial thermal stress of formation, respectively, MPa
σ r s P , σ θ s P , σ z s P Radial, tangential, and axial stress of casing caused by increased casing pressure, respectively, MPa
σ r c P , σ θ c P , σ z c P Radial, tangential, and axial stress of cement sheath caused by increased casing pressure, respectively, MPa
σ r f P , σ θ f P , σ z f P Radial, tangential, and axial stress of formation caused by increased casing pressure, respectively, MPa
σ r c , σθc, σ z c Combined radial, tangential, and axial stress of cement sheath caused by decreased wellbore temperature and increased casing pressure during fracturing fluid injection, respectively, MPa
σ1, σ2, σ3Maximum, intermediate, and minimum principal stresses in the cement sheath, respectively, MPa
uTRadial thermal displacement in the cylinder, mm
uPRadial displacement in the cylinder caused by increased inner pressure, mm
u s T , u c T , u f T Radial thermal displacement for casing, cement sheath, and formation, respectively, mm
u s P , u c P , u f P Radial displacement caused by increased casing pressure for casing, cement sheath, and formation, respectively, mm
E, Es, Ec, EfElasticity modulus of ordinary cylinder, casing, cement sheath, and formation, respectively, MPa
μ, μs, μc, μfPoisson’s ratio of ordinary cylinder, casing, cement sheath, and formation, respectively, dimensionless
α, αs, αc, αfLinear thermal expansion coefficient of ordinary cylinder, casing, cement sheath, and formation, respectively, 1/°C
aInternal radius of thick wall cylinder, mm
rRadial distance from the axis of wellbore, mm
riInside radius of casing, mm
r1Outside radius of casing or inside radius of cement sheath, mm
r2Outside radius of cement sheath or inside radius of formation, mm
roOutside radius of formation, mm
T(r)Temperature change value at the radius of r after fracturing fluid injection, °C
TiWellbore temperature before fracturing fluid injection, °C
TtWellbore temperature after fracturing fluid injection, °C
TeFormation temperature, °C
C1, C1s, C1c, C1fUndetermined coefficients, dimensionless
C2, C2s, C2c, C2fUndetermined coefficients, m2
[A], {B}Coefficient matrix and constant vector for Equation (A5), respectively
σtcTensile strength of cement sheath, MPa
σccUniaxial compressive strength of cement sheath, MPa
τcmaxThe maximum allowable shear stress for cement sheath, MPa
τoctOctahedral shear stress, MPa
ccCohesion of cement sheath, MPa
ϕcAngle of internal friction of cement sheath, (°)
piIncreased casing pressure, MPa
pc1Contact pressure at the casing-cement sheath interface induced by increased casing pressure, MPa
pc2Contact pressure at the cement sheath-formation interface induced by increased casing pressure, MPa
pfChanged formation pressure, MPa
SFSafety factor of cement sheath, dimensionless

Appendix A. Cement Sheath Stress Induced by a Decrease of Wellbore Temperature

Submitting Equations (2) and (3) into Equation (1), the thermal radial displacement and thermal stresses for the combined system can be determined and the detailed derivative results are given below [12].
For the casing (rirr1):
{ u s T = 1 + μ s 1 μ s α s r 1 2 ( T t T i ) ( r 2 r i 2 ) + C 1 s r + C 2 s r σ r s T = α s E s 1 μ s 1 r 2 1 2 ( T t T i ) ( r 2 r i 2 ) + E s 1 + μ s ( C 1 s 1 2 μ s C 2 s r 2 ) σ θ s T = α s E s 1 μ s 1 r 2 1 2 ( T t T i ) ( r 2 r i 2 ) α s E s ( T t T i ) 1 μ s + E s 1 + μ s ( C 1 s 1 2 μ s + C 2 s r 2 ) σ z s T = α s E s ( T t T i ) 1 μ s + 2 μ s E s C 1 s ( 1 + μ s ) ( 1 2 μ s )
For the cement sheath (r1rr2):
{ u c T = 1 + μ c 1 μ c α c r { T t T i 2 ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r r 1 ) r 2 2 + r 1 2 2 ] } + C 1 c r + C 2 c r σ r c T = α c E c 1 μ c 1 r 2 { T t T i 2 ( r 2 r 1 2 ) + ( T i - T t ) 2 ln ( r o / r 1 ) [ r 1 ln ( r r 1 ) r 2 2 + r 1 2 2 ] } + E c 1 + μ c ( C 1 c 1 2 μ c C 2 c r 2 ) σ θ c T = α c E c 1 μ c 1 r 2 { T t T i 2 ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r r 1 ) r 2 2 + r 1 2 2 ] }   α c E c 1 μ c [ T t T i + ( T i T t ) ln ( r / r 1 ) ln ( r o / r 1 ) ] + E c 1 + μ c ( C 1 c 1 2 μ c + C 2 c r 2 ) σ z c T = α c E c 1 μ c [ T t T i + ( T i T t ) ln ( r / r 1 ) ln ( r o / r 1 ) ] + 2 μ c E c C 1 c ( 1 + μ c ) ( 1 2 μ c )
For the formation (r2rro):
{ u f T = 1 + μ f 1 μ f α f r { T t T i 2 ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r / r 1 ) r 2 2 + r 1 2 2 ] } + C 1 f r + C 2 f r σ r f T = α f E f 1 μ f 1 r 2 { T t T i 2 ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r / r 1 ) r 2 2 + r 1 2 2 ] } + E f 1 + μ f ( C 1 f 1 2 μ f C 2 f r 2 ) σ θ f T = α f E f 1 μ f 1 r 2 { T t T i 2 ( r 2 r 1 2 ) + ( T i T t ) 2 ln ( r o / r 1 ) [ r 2 ln ( r / r 1 ) r 2 2 + r 1 2 2 ] }   α f E f 1 μ f [ T t T i + ( T i T t ) ln ( r / r 1 ) ln ( r o / r 1 ) ] + E f 1 + μ f ( C 1 f 1 2 μ f + C 2 f r 2 ) σ z f T = α f E f 1 μ f [ T t T i + ( T i T t ) ln ( r / r 1 ) ln ( r o / r 1 ) ] + 2 μ f E f C 1 f ( 1 + μ f ) ( 1 2 μ f )
Now, let us solve the six undetermined coefficients of C1s, C2s, C1c, C2c, C1f, and C2f in Equation (A1) to Equation (A3). Suppose the combined system is completely cemented at two interfaces. The boundary conditions and the continuity conditions of the radial thermal stress and thermal radial displacement at the two interfaces are as follows.
{ { σ r s T | r = r i = 0 σ r f T | r = r o = 0 { σ r s T | r = r 1 = σ r c T | r = r 1 u s T | r = r 1 = u c T | r = r 1 { σ r c T | r = r 21 = σ r f T | r = r 21 u c T | r = r 2 = u f T | r = r 2
Then submitting Equations (A1) to (A3) into Equation (A4), we get the system of linear equations with unknown number of C1s, C2s, C1c, C2c, C1f, and C2f as follows.
[ A ] × { C 1 s ,   C 2 s , C 1 c , C 2 c , C 1 f , C 2 f } T = { B }
In Equation (A5), the coefficient matrix [ A ] and constant vector { B } can both be obtained from the final actual boundary and continuity conditions. The Gaussian main elimination method was selected for solving Equation (A5) [20]. Then, submitting C1c and C2c into Equation (A2), we can get the radial thermal displacement and thermal stress induced by a decrease of wellbore temperature during fracturing at any radius of r in the cement sheath.

References

  1. Shadravan, A.; Schubert, J.; Amani, M.; Teodoriu, C. HPHT cement sheath integrity evaluation method for unconventional wells. In Proceedings of the SPE International Conference on Health, Safety, and Environment, Long Beach, CA, USA, 17–19 March 2014. [Google Scholar]
  2. Wang, W.; Taleghani, A.D. Cement sheath integrity during hydraulic fracturing: An integrated modeling approach. In Proceedings of the SPE Hydraulic Fracturing Technology Conference, The Woodlands, TX, USA, 4–6 February 2014. [Google Scholar]
  3. Vrålstad, T.; Skorpa, R.; Opedal, N.; Andrade, J.D. Effect of thermal cycling on cement sheath integrity: Realistic experimental tests and simulation of resulting leakages. In Proceedings of the SPE Thermal Well Integrity and Design Symposium, Banff, AB, Canada, 23–25 November 2015. [Google Scholar]
  4. Thiercelin, M.J.; Dargaud, B.; Baret, J.F.; Rodriquez, W.J. Cement design based on cement mechanical response. SPE Drill. Complet. 1998, 13, 266–273. [Google Scholar] [CrossRef]
  5. Jing, L.I.; Lin, C.Y.; Yang, S.C.; Zhi, Y.; Chen, S. Theoretical solution of thermal stress for casing-cement-formation coupling system. J. China Univ. Pet. 2009, 33, 63–69. [Google Scholar]
  6. Li, Y.; Liu, S.; Wang, Z.; Yuan, J.; Qi, F. Analysis of cement sheath coupling effects of temperature and pressure in non-uniform in-situ stress field. In Proceedings of the International Oil and Gas Conference and Exhibition in China, Beijing, China, 8–10 June 2010. [Google Scholar]
  7. Bois, A.P.; Garnier, A.; Galdiolo, G.; Laudet, J.B. Use of a mechanistic model to forecast cement-sheath integrity for CO2 storage. In Proceedings of the SPE International Conference on CO2 Capture, Storage, and Utilization, New Orleans, LA, USA, 10–12 November 2010. [Google Scholar]
  8. Haider, M.G.; Sanjayan, J.; Ranjith, P.G. Modeling of a wellbore composite cylinder system for cement sheath stress analysis in geological sequestration of CO2. In Proceedings of the 46th U.S. Rock Mechanics/Geomechanics Symposium, Chicago, IL, USA, 24–27 June 2012. [Google Scholar]
  9. Bui, B.T.; Tutuncu, A.N. Modeling the failure of cement sheath in anisotropic stress field. In Proceedings of the SPE Unconventional Resources Conference Canada, Calgary, AB, Canada, 5–7 November 2013. [Google Scholar]
  10. Xi, Y.; Li, J.; Liu, G.H.; Zha, C.Q.; Yan, P. Analysis on cement sheath integrity under transient thermo-mechanical coupling effect. Oil Drill. Prod. Technol. 2017, 39, 417–423. [Google Scholar]
  11. Xu, H.; Zhang, Z.; Shi, T.; Xiong, J.Y. Influence of the WHCP on cement sheath stress and integrity in HTHP gas well. J. Pet. Sci. Eng. 2015, 126, 174–180. [Google Scholar]
  12. Xu, H.; Peng, N.; Ma, T.; Yang, B. Investigation of thermal stress of cement sheath for geothermal wells during fracturing. Energies 2018, 11, 2581. [Google Scholar] [CrossRef]
  13. Xu, B.; Wang, J. Theory of Elastic Mechanics; Tsinghua University Press: Beijing, China, 2007. [Google Scholar]
  14. Bois, A.P.; Garnier, A.; Galdiolo, G.; Laudet, J.B. Use of a mechanistic model to forecast cement-sheath integrity. SPE Drill. Complet. 2012, 27, 303–314. [Google Scholar] [CrossRef]
  15. Al-Ajmi, A.M.; Zimmerman, R.W. Relation between the Mogi and the Coulomb failure criteria. Int. J. Rock Mech. Min. Sci. 2005, 42, 431–439. [Google Scholar] [CrossRef]
  16. Ma, T.S.; Chen, P.; Yang, C.H.; Zhao, J. Wellbore stability analysis and well path optimization based on the breakout width model and Mogi-Coulomb criterion. J. Pet. Sci. Eng. 2015, 135, 678–701. [Google Scholar] [CrossRef]
  17. Nelson, E.B.; Guillot, D. Well Cementing, 2nd ed.; Schlumberger: Amsterdam, The Netherlands, 2006. [Google Scholar]
  18. Jaeger, J.C.; Cook, N.G.W.; Zimmerman, R.W. Fundamentals of Rock Mechanics, 4th ed.; Blackwell Publishing: Hoboken, NJ, USA, 2007. [Google Scholar]
  19. Andrade, J.D. Cement Sheath Integrity during Thermal Cycling. Ph.D. Thesis, Norwegian University of Science and Technology, Trondheim, Norway, 2015. [Google Scholar]
  20. Li, Q.; Wang, N. Numerical Analysis; Tsinghua University Press: Beijing, China, 2008. [Google Scholar]
Figure 1. The schematic diagram of casing-cement sheath-formation combined system during fracturing fluid injection.
Figure 1. The schematic diagram of casing-cement sheath-formation combined system during fracturing fluid injection.
Energies 11 03534 g001
Figure 2. Failure modes in the cement sheath. (a) Inner de-bonding; (b) Outer de-bonding; (c) Radial cracking; (d) Disking; (e) Shear damage.
Figure 2. Failure modes in the cement sheath. (a) Inner de-bonding; (b) Outer de-bonding; (c) Radial cracking; (d) Disking; (e) Shear damage.
Energies 11 03534 g002
Figure 3. The 2D axisymmetric finite element analysis model of casing-cement sheath-formation system.
Figure 3. The 2D axisymmetric finite element analysis model of casing-cement sheath-formation system.
Energies 11 03534 g003
Figure 4. Comparison of stresses calculated with the analytical model and the finite element analysis model.
Figure 4. Comparison of stresses calculated with the analytical model and the finite element analysis model.
Energies 11 03534 g004
Figure 5. Radial distribution of the radial stress in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Figure 5. Radial distribution of the radial stress in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g005
Figure 6. Safety factor of cement sheath for de-bonding at two interfaces under different increased casing pressures pi and decreased temperature of 70 °C at the casing inner wall.
Figure 6. Safety factor of cement sheath for de-bonding at two interfaces under different increased casing pressures pi and decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g006
Figure 7. Radial distribution of the tangential stress in the cement sheath under different increased casing pressure pi and decreased temperature of 70 °C at the casing inner wall.
Figure 7. Radial distribution of the tangential stress in the cement sheath under different increased casing pressure pi and decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g007
Figure 8. Radial distribution of the cement sheath safety factors for radial cracking under different increased casing pressures pi and decreased temperature of 70 °C at the casing inner wall.
Figure 8. Radial distribution of the cement sheath safety factors for radial cracking under different increased casing pressures pi and decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g008
Figure 9. Radial distribution of the axial stress in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Figure 9. Radial distribution of the axial stress in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g009
Figure 10. Radial distribution of the cement sheath safety factors for disking under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Figure 10. Radial distribution of the cement sheath safety factors for disking under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g010
Figure 11. Radial distribution of the octahedral shear stress τoct in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Figure 11. Radial distribution of the octahedral shear stress τoct in the cement sheath under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g011
Figure 12. Radial distribution of the cement sheath safety factors for shear failure under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Figure 12. Radial distribution of the cement sheath safety factors for shear failure under different increased casing pressures pi and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g012
Figure 13. Radial distribution of the failure stress and corresponding strength of the cement sheath under only a decreased temperature of 70 °C at the casing inner wall.
Figure 13. Radial distribution of the failure stress and corresponding strength of the cement sheath under only a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g013
Figure 14. Radial distribution of the failure stress and corresponding strength of the cement sheath under an increased casing pressure pi of 70 MPa and a decreased temperature of 70 °C at the casing inner wall.
Figure 14. Radial distribution of the failure stress and corresponding strength of the cement sheath under an increased casing pressure pi of 70 MPa and a decreased temperature of 70 °C at the casing inner wall.
Energies 11 03534 g014
Figure 15. Radial distribution of the cement sheath safety factors for radial cracking under a different elasticity modulus of cement sheath.
Figure 15. Radial distribution of the cement sheath safety factors for radial cracking under a different elasticity modulus of cement sheath.
Energies 11 03534 g015
Figure 16. Radial distribution of the cement sheath safety factors for disking under a different elasticity modulus of cement sheath.
Figure 16. Radial distribution of the cement sheath safety factors for disking under a different elasticity modulus of cement sheath.
Energies 11 03534 g016
Figure 17. Radial distribution of the cement sheath safety factors for shear failure under a different elasticity modulus of cement sheath.
Figure 17. Radial distribution of the cement sheath safety factors for shear failure under a different elasticity modulus of cement sheath.
Energies 11 03534 g017
Table 1. Safety factor for cement sheath.
Table 1. Safety factor for cement sheath.
Failure ModeFailure CriterionSafety factor (SF)
Inner de-bonding σ r c | r = r 1 > σ t c S F = σ t c / σ r c | r = r 1
Outer de-bonding σ r c | r = r 2 > σ t c S F = σ t c / σ r c | r = r 2
Radial cracking σ θ c > σ t c S F = σ t c / σ θ c
Disking σ z c > σ t c S F = σ t c / σ z c
Shear failure τ o c t > τ c max S F = τ c max / τ o c t
Table 2. Basic calculation parameters.
Table 2. Basic calculation parameters.
No.SymbolValueUniteNo.SymbolValueUnit
1ri60.68mm11μs0.30dimensionless
2r169.85mm12μc0.19dimensionless
3r2107.95mm13μf0.21dimensionless
4ro1079.5mm14Ti100°C
5αs1.15 × 10−51/°C15Tt30–70°C
6αc1.03 × 10−51/°C16Te100°C
7αf1.03 × 10−51/°C17pi0–70MPa
8Es206GPa18σtc2.3MPa
9Ec10GPa19σcc32MPa
10Ef30GPa20φc30°

Share and Cite

MDPI and ACS Style

Xu, H.; Ma, T.; Peng, N.; Yang, B. Influences of Fracturing Fluid Injection on Mechanical Integrity of Cement Sheath under Four Failure Modes. Energies 2018, 11, 3534. https://doi.org/10.3390/en11123534

AMA Style

Xu H, Ma T, Peng N, Yang B. Influences of Fracturing Fluid Injection on Mechanical Integrity of Cement Sheath under Four Failure Modes. Energies. 2018; 11(12):3534. https://doi.org/10.3390/en11123534

Chicago/Turabian Style

Xu, Honglin, Tianshou Ma, Nian Peng, and Bin Yang. 2018. "Influences of Fracturing Fluid Injection on Mechanical Integrity of Cement Sheath under Four Failure Modes" Energies 11, no. 12: 3534. https://doi.org/10.3390/en11123534

APA Style

Xu, H., Ma, T., Peng, N., & Yang, B. (2018). Influences of Fracturing Fluid Injection on Mechanical Integrity of Cement Sheath under Four Failure Modes. Energies, 11(12), 3534. https://doi.org/10.3390/en11123534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop