The Scientific World Journal

Volume 2013, Article ID 565891, 9 pages

http://dx.doi.org/10.1155/2013/565891

## Analyzing Axial Stress and Deformation of Tubular for Steam Injection Process in Deviated Wells Based on the Varied Fields

^{1}Uncertainty Decision-Making Laboratory, Sichuan University, Chengdu 610064, China^{2}College of Economics & Management, Sichuan Agricultural University, Chengdu 611130, China^{3}Research School of Engineering Technology, The Southwest Petroleum and Gas Corp, China Petroleum and Chemical Corp, Deyang 618000, China

Received 13 May 2013; Accepted 2 August 2013

Academic Editors: G. Carbone, S. Park, S. Torii, and Q. Yang

Copyright © 2013 Yunqiang Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The axial stress and deformation of high temperature high pressure deviated gas wells are studied. A new model is multiple nonlinear equation systems by comprehensive consideration of axial load of tubular string, internal and external fluid pressure, normal pressure between the tubular and well wall, and friction and viscous friction of fluid flowing. The varied temperature and pressure fields were researched by the coupled differential equations concerning mass, momentum, and energy equations instead of traditional methods. The axial load, the normal pressure, the friction, and four deformation lengths of tubular string are got ten by means of the dimensionless iterative interpolation algorithm. The basic data of the X Well, 1300 meters deep, are used for case history calculations. The results and some useful conclusions can provide technical reliability in the process of designing well testing in oil or gas wells.

#### 1. Introduction

The deviated wells had been wildly applicable for petroleum and natural gas industry. Deviated wells have their distinctive characteristics which are distinguished from that of other wells. (1) High temperature high pressure: the temperature distribution and pressure on the tubing are significantly different when outputs are varied (flow velocity) but neither has a simple linear relationship, because the fluid density is not constant. (2) Deep well: the sensibility of force and deformation influencing by the factors, such as the temperature, pressure, density of fluid, viscous friction and fluid velocity, and so forth, will become high with the increase of tubing length. The completion test of a deep well is a new problem. In the research of applied basic theory for deep well testing, tubular string mechanical analysis is very complex, but fluid temperature and tubing pressure affect the force of the tubular string heavily. Temperature, pressure, liquid density, and fluid velocity within tubing may change with of the hole depth, time, and operations, so that the axial force changes constantly. A large compression load at low end can induce the tubing plastic deformation and make the packer damaged. A large tension load at the top end may unpack the packer or cause the tubing to break. If the tubing failed, the whole borehole can hardly maintain its integrity and safety [1]. Therefore, it is very important for deviated wells to predict the axial forces for the safety.

Hammerlindl [2] had made a great contribution about tubular mechanics. He had put forth the four effects between the packer forces and length change of tubing: temperature effect, ballooning effect, axial load effect, and the helical buckling effect. There is a large amount of papers to research the effect of buckling behavior. Therefore it is considered that inflexion is caused on its axial force under certain conditions, by which colliding on parts of the drill string with well bore is induced. When buckled of tubular beyond wellhole’s control, the buckling configuration which will be transformed at the state of stabilization, sinusoidal buckling and helical buckling with the increase of load. The problem of buckling of the tube was first studied and put into practice by Lubinski et al. [3]. They had done the emulation experiment for the buckling behavior of tube in deviated wells and found the compute formula on critical buckling load of tube in deviated wells. Paslay and Bogy [4] found that the number of sinusoids in the buckling mode increases with the length of the tube. The buckling behavior by inner and outer fluid pressure of tubing was analyzed, and the mathematical relation between pitch and axial pressures was deduced based on the principle of minimum potential energy (see Hammerlindl [2]). The mptotic solution for sinusoidal buckling of an extremely long tube has been analyzed by Dawson and Paslay [5], based on a sinusoidal buckling mode of constant amplitude. Numerical solutions were also sought by Mitchell [6] using the basic mechanics equations. His solutions confirm the thought that, under a general loading, the deformed shape of the tube is a combination of helices and sinusoids while helical deformation occurs only under special values of the applied load. The formula about tubing forces had been put, however, which is too simple for shallow wells to accommodate the complicated states of deep wells. Up to now, many researches are centered on water injection tubular but not on steam injection. Among them, the values of temperature and pressure are considered as constant or lineal functions which will cause large errors on tubular deformation computing [7].

In fact, the tubular string deformation includes transverse deformation and longitudinal deformation. Because the transverse length (its order of magnitude is ) is much and much smaller than the longitudinal length (its order of magnitude is ), we mainly consider the axial (longitudinal) deformation for the tubular string deformation analysis in the paper. In the paper, the force states of tubular in the process of steam injection are analyzed. The varied fields are considered to compute the values of several deformations. The axial load and four deformation lengths of tubular string are obtained by the dimensionless iterative interpolation algorithm. The basic data of the X Well (deviated well), 1300 meters deep in China, are used for case history calculations. Some useful suggestions are drawn.

This paper is organized as follows. Section 2 gives a system model about tubular mechanics and deformation. And the varied fields were presented by model concerning mass, momentum, and energy balance. Section 3 gives the parameters, initial condition, and algorithm for solving model. In Section 4, we give an example from a deviated well at 1300 meters of depth in China, and the result analysis are made. Section 5 gives a conclusion.

#### 2. Model Building

##### 2.1. Basic Assumption

Before analyzing the force on the microelement, some assumptions are introduced as follows:(1)the curvature of the hole of the considered modular section is constant,(2)on the upper side or underside of the section which is point of contact of the pipe and tube wall, the curvature is the same with the hole curvature,(3)the radius of steam injection string, in contrast to curvature of borehole, is insignificant,(4)the string is at the state of linear elastic relationship.

##### 2.2. Forces Analysis of Tubular String

The forces of tubular string are shown in Figure 1. Consider the flow system depicted in Figure 1: a constant cross-sectional flow area , inner diameter , outer diameter , material density , packer fluid density , and a total length . Through this tubing gas flows from the bottom to the top with a mass flow rate . The distance coordinate in the flow direction along the tubing is denoted by . The cylindrical coordinate system , origin of which is in wellhead and axis is down as the borehole axis, is used.

As shown from Figure 1, the tubular string is mainly acted upon by the following forces at the process of steam injection.(1)*Initial Axial Force*. The initial axial force of tubular should include the deadweight, buoyant weight, and initial pull force.(2)*Thermal Stress*. On the process of steam injection, the temperature stress will act at the tubular with varied temperature.(3)*Axial Force by the Varied Internal and External Pressure*. Thanks to the varied pressure with internal and external pressure, the tubular will be acted by the bending force, piston force, and other axial forces.(4)*Friction Drag by Steam Injection*. On the process of steam injection, the flow in tubular will produce viscous flow which will cause the friction drag.

##### 2.3. The Axial Load and Axial Stress of the Tubular

###### 2.3.1. Initial Axial Load and Initial Axial Stress of Steam Injection Tubular

*Initial Axial Load*. The section to which the distance from the wellhead is was considered. The axial static load by the deadweight of tubular is as follows:
where is the deadweight of tubular, is the average unit length weight of tubing, is the length of tubular, is the density of tubular, and is the inclination angle.

The axial static load by the buoyant weight is as follows: where is the buoyant weight of tubular, is the density of packer fluid.

The axial load by the steam injection pressure where, represents the inner pressure at this section.

Therefore, summing (1), (2), and (3), the axial forces in the section are obtained as follows:
*Initial Axial Stress*. The axial stress can be derived from the following equation

###### 2.3.2. Axial Thermal Stress of Steam Injection Tubular

In the process of steam injection, the temperature of tubular will change with time and depth, which will make the tubular deform as follows: where, represents the steel elastic modulus of tubular, is the warm balloon coefficient of the tubular string, and is the temperature change with before and after steam injection.

###### 2.3.3. Axial Stress of Steam Injection Tubular by the Change with Pressure

The effect acting the tubular with pressure change which is called ballooning effect normally.

*Ballooning Stress Analysis*. The ballooning effect will be produced from pressure acted in inner and outer of the tube. Generally, there are two kinds of tubular in oil wells. One is the tubulars whose outer diameter is 88.9 mm, inner diameter is 76 mm, and thickness of tubes is 6.5 mm ; the other is the tubular whose outer diameter is 114.3 mm, inner diameter is 100.5 mm, and thickness of tubes is 6.9 mm . Neither is the thin-wall problem. Therefore, it should be solved by Lame’s formula [8].

The radial and tangential stresses in the thick-wall cylinder can be shown as Figure 2. The two can be calculated as follows: where is radial stress, is tangential stress, is radial coordinate, is tube internal pressure at point, and is tube external pressure at point.

###### 2.3.4. Axial Stress of Steam Injection Tubular by the Friction Loss

In fact, the flow in the tubular should be multiflow. On the process of steam injection, the flow will be run and it will give rise to friction effect to cause axial stress. In our paper, we consider the flow gas-liquid mix flow and the liquid head loss is gotten by the Darcy-Weisbach formula [9] as follows: where means heat loss of liquid flow, is frictional head losses coefficients, and is the velocity of liquid flow.

The friction drag in tubular is ( is density of liquid flow). The axial stress by fiction drag can be obtained as follows:

##### 2.4. Analysis of Axial Deformation

Based on the studies and analyses mentioned above, the axial deformation on the tubular is made up of the following parts.

###### 2.4.1. The Axial Deformation by the Axial Static Stress

For the microelement of the tubular , the unit deformation by the static stress can be computed by generalized Hooke law where represents Poisson’s ratios.

The axial deformation at an element can be obtained through integrating on the length of the element as follows:

Therefore, the total axial deformation by the static stress can be gotten accumulating each element as follows:

###### 2.4.2. The Axial Deformation with Temperature Changed

For the microelement of the tubular , the unit deformation by the temperature change is as follows:

The same principle is that the total axial deformation by the varied temperature fields can be gotten accumulating each element as follows:

###### 2.4.3. The Axial Deformation with the Friction Drag

For the microelement of the tubular , the unit deformation by the friction force is as follows:

###### 2.4.4. The Axial Deformation with the Tubular String Buckling

Researchers in general call the buckling a bending effect. The tubular is freely suspended in the absence of fluid inside as shown in Figure 3(a). Because the force applied at the end of the tubular which is large enough, the tubular will buckle as shown in Figure 3(b).

Lubinski et al. [3] had done many researches on the phenomenon. From their work, we can get the buckling effect. Define the virtual axial force of tubing as follows: where is the pressure inside the tubular at the packer length, is the pressure outside the tubular at the packer length, and is the area corresponding to packer bore.

By (16), whether the tubular will buckle or not can be judged. The string will buckle if is positive or remain straight if is negative or zero. The axial deformation of the tubular string buckling is where means tubing-to-casing radial clearance, is moment of inertia of tubing cross-section with respect to its diameter , denotes change with before and after injection, and is the unit weight of tubing, as In addition, the position of the neutral point is needed. The length () from the packer to the point can be computed as follows:

Generally, the neutral point should be in tubular (). However, at the multipackers, it will occur that the neutral point is outside the tubing between dual packers. In this paper, we leave the latter phenomenon.

To sum up, the whole deformation length can be represented as follows:

##### 2.5. The Analysis of the Varied Fields

In the course of dryness modeling, we can find that the numerical values of deformation ((10), (13), and (17)) were affected by the temperature and pressure. In fact, the two parameters varied according to the depth and time changing. So, the varied fields need to be researched. Under the China Sinopec Group Hi-Tech Project “Stress analysis and optimum design of well completion” in 2009 [6] undertaken by Sichuan University at early time. The varied fields had been deduced strictly based on the mass, momentum, and energy balance. The proof details can be shown in Xu et al. [11]. The varied fields is

#### 3. Numerical Implementation

##### 3.1. Calculation of Some Parameters

In this section, we will give the calculating method of some parameters.(1)Each point’s inclination: where represents segment point of calculation, represents measurement depth of inclination angle , and , is the step length of calculation. Transient heat transfer function [12]: (2)The density of wet steam. Since the flow of the water vapor in is the gas-liquid two-phase flow, there are many researches about this problem [13, 14]. In the paper, we adopt the M-B model to calculate the average density of the mixture.(3)The heat transfer coefficient from different positions of the axis of the wellbore to the second surface.

These resistances include the tubing wall, possible insulation around the tubing, annular space (possibly filled with a gas or liquid but is sometimes vacuum), casing wall, and cementing behind the casing as follows: and are the heat conductivity of the heat insulating material and the cement sheath, respectively. and are the coefficients of the convection heat transfer and the radiation heat transfer.

##### 3.2. Initial Condition

In order to solve model, some definite conditions and initial conditions should be added. The initial conditions comprise the distribution of the pressure and temperature at the well top. In this paper, we adopt the value at the initial time by actual measurement. Before steam injected, the temperature of tubular just is initial temperature of formation (, is geothermal gradient). At the same time, the pressure of inner tubular is assumed to be equal to the outer tubular before steam injected.

##### 3.3. Steps of Algorithm

To simplify the calculation, we divided the wells into several short segments of the same length. The length of a segment varies depending on variations in wall thickness, hole diameter, fluid density inside and outside the pipe, and wells geometry. The model begins with the calculation at one particular position in the wells: the top of the pipe.

*Step 1. *Set step length of depth. In addition, we denote the relatively tolerant error by . The smaller , is, the more accurate the results are. However, it will lead to rapid increasing calculating time. In our paper, we set (m), and .

*Step 2. *Give the initial conditions.

*Step 3. *Compute each point’s inclination.

*Step 4. *Compute the parameters under the initial conditions or the last depth variables.

*Step 5. *Let ; then we can get the by solving the following equation:

Let be the temperature at the injection time and radial at the depth . We apply the finite different method to discretize the equations as follows:
where is the interval of time and is the interval of radial, respectively. It can be transformed into the standard form as follows:

Then the different method is used to discretize the boundary condition. For , we have

For , we have
We can compute the symbolic solution of the temperature of the stratum. In this step, we will get the discrete distribution of as the following matrix:
where represents the injection time and represents the radial.

*Step 6. *Let the right parts of the coupled differential equations be functions , where . Then we can obtain a system of coupled functions as follows:
where at .

*Step 7. *Assume that are , respectively. Then we can obtain some basic parameters as follows:

*Step 8. *Calculate the pressure and temperature at point :

*Step 9. *Calculate the deformation , , and by previous equations.

*Step 10. *Repeat the third step to the tenth step until tubular length is calculated.

*Step 11. *Calculate the deformation and total deformation length as follows:

#### 4. Numerical Simulation

##### 4.1. Parameters

To demonstrate the application of our theory, we study a pipe in X well, which is in Sichuan Province, China. All the basic parameters are given as follows: depth of the well is 1300 m; ground thermal conductivity parameter is ; ground temperature is ; ground temperature gradient is ; roughness of the inner surface of the well is ; and parameters of pipes, inclined well, inclination, azimuth, and vertical depth are given in Tables 1, 2, and 3.

##### 4.2. Main Results and Results Analysis

After calculation, we obtain a series of results of this well as Table 4. The influence of outputs on the axial deformation of tubing was investigated as shown by Figure 4.

From the results as shown in Figure 4 and Table 4, some useful analysis can be drawn.(1)The amount of steam injected and inject pressure affected the stretching force with special severity.(2)The results were as follows: the length of tubular deformation was risen with increased injected pressure or injected velocity.(3)The length of tubular deformation increases with the increasing of outputs but more slowly.(4)The thermal stress is the main factor influencing the tubular deformation. Therefore, the temperature of steam injected should not be too high.(5)The lifting prestressed cementing technology has important meanings to reduce the deformation of tubular.(6)The creeping displacement of downhole stings will produce an upward contractility which causes packer depressed or lapsed. Therefore, the effective measures should be adopted to control the companding of tubular.

#### 5. Conclusion

In this paper, a total tubular deformation model about deviated wells was given. A coupled-system model of differential equations concerning pressure and temperature in high temperature-high pressure steam injection wells according to mass, momentum, and energy balances, which can reduce the error of axial stress and axial deformation, was given instead of the average value or simple linear relationship in traditional research. The basic data of the Well (high temperature and high pressure gas well), 1300 m deep in Sichuan, China, were used for case history calculations. The results can provide technical reliance for the process of designing well tests in deviated gas wells and dynamic analysis of production.

#### Nomenclature

: | Inner diameter (m) |

: | Microelement of the tubular |

: | Acceleration of gravity (m/s^{2}) |

: | Depth of top tubular located at the packer (m) |

: | Time of down stroke (s) |

: | Dimensionless time (dimensionless) |

: | Velocity of fluid in tubing (m/s) |

: | Velocity of down stroke (m/s) |

: | Distance coordinate in the flow direction (m) along the tubing |

: | Constant cross-sectional flow area (m^{2}) |

: | Effective area (m^{2}) |

: | Area corresponding to packer bore (m^{2}) |

: | The Joule-Thomson coefficient (dimensionless) |

: | Heat capacity of fluids (J/Kg · K) |

: | Outer diameter (m) |

: | Steel elastic modulus of tubular (Mpa) |

: | Axial forces in the section (N) |

: | Axial tensile strength (N) |

: | Friction force (N) |

: | Piston force for supporting packer’s pressure (N) |

: | Pumping force (N) |

: | Length of tubular (m) |

: | Buoyant weight of tubular (Kg) |

: | Dead weight of tubular (Kg) |

: | Pressure outside the tubular (Mpa) |

: | Pressure inside the tubular at the packer length (Mpa) |

: | Pressure in tubing (Mpa) |

: | Temperature in tubing (°C) |

: | Initial temperature of formation (°C) |

: | Mass flow rate (Kg/s) |

: | Total length (m) |

: | Material density (Kg/m^{3}) |

: | Packer fluid density (Kg/m^{3}) |

: | Inclination angle (∘) |

: | Warm balloon coefficient of the tubular string (dimensionless) |

: | Drop of any parameter |

: | Axial thermal stress () |

: | Density of fluid in the tubing (Kg/m^{3}) |

: | The tubular string buckling axial deformation (m) |

: | Total axial deformation by varied temperature fields (m) |

: | Total axial deformation by the varied pressure fields (m) |

: | Differential pressure at startup (Mpa) |

: | Change in tubing pressure at the length (Mpa) |

: | Change in annulus pressure at the length (Mpa) |

: | Differential pressure from top to bottom (Mpa) |

: | Temperature change with before and after well shut-in (°C) |

: | Change in density of liquid in the tubing at the length (Kg/m^{3}) |

: | Change in density of liquid in the casing at the length (Kg/m^{3}). |

#### Acknowledgments

This research was supported by the Key Program of NSFC (Grant no. 70831005) and the Key Project of China Petroleum and Chemical Corporation (Grant no. GJ-73-0706).

#### References

- D.-L. Gao and B.-K. Gao, “A method for calculating tubing behavior in HPHT wells,”
*Journal of Petroleum Science and Engineering*, vol. 41, no. 1–3, pp. 183–188, 2004. View at Publisher · View at Google Scholar · View at Scopus - D. J. Hammerlindl, “Movement, forces, and stresses associated with combination tubing strings sealed in packers,”
*Journal of Petroleum Technology*, vol. 29, pp. 195–208, 1977. View at Google Scholar · View at Scopus - A. Lubinski, W. S. Althouse, and J. L. Logan, “Helical buckling of tubular sealed in packers,”
*Journal of Petroleum Technology*, vol. 14, no. 6, pp. 655–670, 1962. View at Google Scholar - P. R. Paslay and D. B. Bogy, “The stability of a circular rod laterally constrained to be in contact with an inclined circular cylinder,”
*Journal of Applied Mechanics*, vol. 31, pp. 605–610, 1964. View at Google Scholar - R. Dawson and P. R. Paslay, “Drillpipe buckling in inclined holes,”
*Journal of Petroleum Technology*, vol. 36, no. 10, pp. 1734–1738, 1984. View at Google Scholar · View at Scopus - R. F. Mitchell, “Effects of well deviation on helical buckling,”
*SPE Drilling and Completion*, vol. 12, no. 1, pp. 63–68, 1997. View at Google Scholar · View at Scopus - P. Ding and X. Z. Yan, “Force analysis of high pressure water injection string,”
*Petroleum Dring Techiques*, vol. 36, no. 5, p. 23, 2005. View at Google Scholar - Z. F. Li, “Casing cementing with half warm-up for thermal recovery wells,”
*Journal of Petroleum Science and Engineering*, vol. 61, no. 2–4, pp. 94–98, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. M. Sun, “The analysis and computing of water injection tubular,”
*Drilling and Production Technology*, vol. 26, no. 3, pp. 55–57, 2003 (Chinese). View at Google Scholar - J. P. Xu, “Stress analysis and optimum design of well completion,”
*Technical Report of Sinopec*GJ-73-0706, 2009. View at Google Scholar - J. P. Xu, Y. Q. Liu, S. Z. Wang, and B. Qi, “Numerical modelling of steam quality in deviated wells with variable (T, P) fields,”
*Chemical Engineering Science*, vol. 84, pp. 242–254, 2012. View at Google Scholar - A. R. Hasan and C. S. Kabir, “Two-phase flow in vertical and inclined annuli,”
*International Journal of Multiphase Flow*, vol. 18, no. 2, pp. 279–293, 1992. View at Google Scholar · View at Scopus - H. D. Beggs and J. R. Brill, “A study of two-phase flow in inclined pipes,”
*Journal of Petroleum Technology*, vol. 25, no. 5, pp. 607–617, 1973, paper 4007-PA. View at Google Scholar · View at Scopus - H. Mukherjee and J. P. Brill, “Pressure drop correlations for inclined two-phase flow,”
*Journal of Energy Resources Technology*, vol. 107, no. 4, pp. 549–554, 1985. View at Google Scholar · View at Scopus