Axial Loads
In this document we describe the calculation machinery behind axial loads calculations in Oliasoft WellDesign™
Last updated
Was this helpful?
In this document we describe the calculation machinery behind axial loads calculations in Oliasoft WellDesign™
Last updated
Was this helpful?
In this section we describe the calculation machinery behind axial loads calculations in Oliasoft WellDesign™. There are essentially two inputs that control the calculations, i.e. internal & external pressure and a temperature profile. The calculations are relative to some condition, normally the initial conditions which usually consist of the hydrostatic pressure and the geothermal temperature. The calculations can be divided into five separate calculations, ballooning- & piston- effects due to change in internal/external pressure, thermal effects, bending forces from the trajectory, and buckling effects. Axial loads calculations can be applied to arbitrary string configurations, and both packers, cement, and combinations are included.
Note: In this documentation we denote any tubular as tubing. All calculations however encompass any tubular, such as tubings, casings, liners, tie-backs etc.
Axial loads calculations are done on a per string basis, meaning that the string under consideration is not affected by other casing strings in the design. Necessary input, is a complete description of the casing string. That is, a complete description of the wellbore trajectory where the casing string is defined, including measured depth, true vertical depth, inclination, azimuth, and dogleg severity, and also where, if any, there are packers or cement. For packers, information about allowable movement and seal bore area is also necessary. Also, the dimensions of the casing are needed, i.e. inner- and outer- diameters, weight per meter, yield, anisotropy, and asymmetry. In addition, initial internal- & external- pressure profile, initial temperature profile, and initial axial loads are needed, which usually stem from initial conditions. When all these are given, the axial loads calculations are based on three inputs, internal- & external- load case pressure profiles, and load case temperature profile. All the calculations are relative to the initial pressure and temperature profile.
For a detailed description of the calculations in this document, we refer to [2].
Consider a tubular with inside- and outside- diameter and , respectively, implying an inside- and outside- area of
Let , denote the cross sectional area. If an axial force, , is applied to the tubing, then the axial stress is given by
Such a force also elongates or deforms the tubing, and the corresponding axial strain, , is defined as
where is the original length of the tubing, and is the length change due to .
Please note, in this document axial tensile forces are defined positive, while axial compressive forces are defined negative.
Elasticity is the property of a material to retain its original shape once the load is relieved. Hooke’s law states that the stress in a material is proportional to the strain which produced it, up to a limit called yield stress. The relation is expressed as
here the proportionality constant, , is called Young’s modulus.
An axial force, , applied to a tubing, not only generates an axial strain, , but also generates a radial strain, . In the elastic regime, where the tubing is not permanently deformed if the load is relieved, these quantities are proportional, with proportionality constant, , called Poisson's ratio, i.e.
We consider a tubing of length , and we assume all the dimensions of the the tubing to be known. Specifically, let denote the inside area of the tubing, denote the outside area of the tubing, and ,
Assume the material in the tubing has Young's modulus, , Poisson's ratio, , and coefficient of (linear) thermal expansion, .
Initially, the tubing is exposed to internal- and external- pressure profiles, and , respectively. All pressure related calculations below are relative to these pressures.
Initially, the tubing is exposed to a temperature profile, . All temperature related calculations are relative to this temperature profile.
Initially, the tubing has a known axial profile (usually from initial conditions). All axial loads calculations are relative to this initial profile.
For the pressure related calculations, we assume the tubing is exposed to an internal- and external- pressure profile, and , respectively. Let and .
For the temperature related calculations, we assume the tubing is exposed to a temperature profile, . Let .
Radial length change of a tubing, caused by a change in pressure, also results in a length change of the tubing, given by
In this calculation, we assume in addition that the tubing has a constant outer diameter, and a uniform wall thickness. If the tubing consists of sections with varying outer diameter and wall thickness, the calculation above are done per section, and then added together.
If the tubing is free to move, there is no axial force associated with ballooning. If, on the other hand, the tubing is fixed, this length change is converted to force through Hooke's law.
Piston effects are pressure-area effects that both results in change of axial forces and movement (if allowed). This effect, occurs where there is a change in geometry of the tubing, i.e. at any crossovers, at the bottom, at plugs inside the tubing, and at expansion devices. In other words, if there is a pressure, , across an area, , the resulting force, normal to the surface defined by , is given by
If movement is allowed, the length change due to this force, is given by Hooke’s law, equation 4,
A crossover occurs whenever there is a change in the diameters of the string, e.g. from a larger tubing above to a smaller tubing underneath, or vice versa. Also, if the outside diameter remains the same, and the inner diameter changes, internal pressure generates a point load. Denote by and , the differential inner- and outer- area, respectively, across such a crossover. Then, the crossover force is given by
This force act at the crossover, hence creates a discontinuity in axial load at the crossover.
A tubing plug occupies the internal area, , of the tubing, and if there is a differential pressure across the plug, a piston force occurs. This force is simply given by
where , and the and denote above and below the plug, respectively. The force act at the plug and creates a discontinuity in axial load at the plug.
At the base of the tubing, assumed open ended, we get a compressive axial force due to pressure acting on the cross sectional area, , given by
where denotes the differential external pressure at the base of the tubing.
Expansion devices come in different flavors, however, the piston force due to such a device is device independent, and is given by
where is the seal bore area.
If the temperature profile around a tubing is changed, then the length of the tubing changes according to
As for ballooning, if the tubing is free to move, there is no axial force associated with a change in temperature. Also, if the tubing is fixed, this length change is converted to force through Hooke's law, equation 4,
Cement is assumed to prevent axial movement, hence, in cemented areas of a string, ballooning, piston effects, and thermal effects directly create stress.
With the framework described above, a new axial load can be calculated based on the change in pressures and temperature. This is done using the following algorithm
a) Divide the casing string into subsections, e.g. in 100 m intervals, and use the end points of these intervals as calculation points. In addition, any change in the casing configuration needs to be included, including change in wellbore trajectory, i.e. calculation points at packers, cement, crossovers, plugs, change in dogleg severity, etc.
b) Calculate the length change of the tubing, if free to move, from ballooning and thermal effects.
c) Calculate axial forces and associated length changes of the tubing, if free to move, from piston effects.
d) Identify all points in the casing string which are fixed or have limited allowable movement, typically packers or cemented areas. Then convert the cumulative movement of the string, from ballooning, thermal-, and piston- effects, to an axial force at these points, using Hooke's law, equation 4. Note, if there is limited movement, and the cumulative movement is greater than the allowable movement, the calculation of a restoring force has to take the allowable movement into account.
e) Calculate a final axial load, not including bending from dogleg, by summing up all the contributions and add it to the initial axial load. At this point, it is important to acknowledge the direction of the forces, and also which parts of the tubing string to be affected by the different restoring forces.
In deviated wellbores, the effect of wellbore curvature must be considered when the axial load on a tubing, inside the wellbore, is calculated. When a tubing is forced to bend, a tensile/positive stress occurs on the outside of the bend, and a compressive/negative stress occurs on the inside of the bend. Assume the pipe takes the form of a circular arc, with radius of curvature . The radius of curvature is the inverse of the dogleg severity, , i.e . Then, the bending stress at the outside of the pipe is
where and sign indicates outside and inside of the bend, respectively.
Replace with to get the bending stress at the inside of the pipe.
The associated bending force is given by
Let denote axial load no bending, as calculated above. Axial load with bending, , from dogleg is then given by
We observe that bending stress is a purely geometrical effect, and is different from zero if the dogleg severity is different from zero, making it also a local effect.
Buckling is a more complex phenomenon than the effects described above. In the following, a short description of the mathematics involved is given, together with an algorithm for calculating the effects of buckling. For a more comprehensive description, reference is made to [1] and [2].
When a tubing is in compression, it will shorten according to Hooke's law. However, this is not the only effect of compression, if the compression force is big enough. At some point a critical compressional force is reached where the tubing is in an unstable condition, and will tend to buckle. It should be remarked that buckling is also dependent on internal- and external- pressure, and in theory, a tubing under tension can buckle if the internal pressure is big enough.
The following quantities are involved in buckling calculations:
In the presence of internal- and external- pressures, and , respectively, the effective axial load on the tubing is given by
where denotes the axial load without bending from dogleg.
When a tubing is displaced in a fluid, the effective buoyed pipe distributed weight is given by
where is the weight of the tubing in and are the internal- and external fluid density, respectively, and is the gravitational constant.
From this we calculate the casing contact load, given by
where and denote the inclination and azimuth of the wellbore trajectory, respectively, and the derivations are with respect to curve length, i.e. measured depth.
The critical compression force is given by
where is the second moment of area, given by , and is the radial clearance, i.e. the difference in radius between the inside of the wellbore wall and the outside of the tubing.
It is the ratio between and that determines whether there is buckling, or not, and also what type of buckling. In axial load calculations, we distinguish between no buckling/neutral , lateral buckling , and helical buckling , and the conditions are
Buckling related quantities we calculate are, the helix angle and the resulting dogleg severity, and , strain, , bending moment, , and bending stress, . The formulas for these quantities are related to the buckling mode. In the neutral mode, there is no buckling, hence all of these are zero.
The bend of the helix, aka the helix angle, is not constant for lateral buckling, hence the maximum is approximated. For helical buckling, on the other hand, this angle is constant.
The approximate solution for the maximum helix angle for lateral buckling is (ref 1)
The helix angle for helical buckling is given by
The resulting dogleg severity is given by
where is either or , depending on the buckling mode.
Buckling reduces the length of the tubing, and the related buckling strain, , is useful, which is a function of the helix angle and radial clearance.
For lateral buckling, the helix angle is not constant, and an average is used for the buckling strain
For helical buckling, the buckling strain is
The bending moment is directly related to the dogleg severity, and is given by
where is either or , depending on the buckling mode.
The related bending stress on the outside of the pipe is given by
Replace with to get the bending stress at the inside of the pipe.
Since bending stress from buckling is dependent on the axial load, and the quantity we want to compute is axial load including bending stress from buckling, the effect of buckling needs to be calculated in an iterative manner. One way to proceed, is to iterate on the length change due to buckling, and stop when the difference in length change between two iterations is less than some tolerance. Explicitly
Calculate axial load no bending, .
Calculate the length change from buckling, , due to .
Calculate a revised axial load, due to (if restricted movement).
Calculate the length change from buckling, , due to .
Calculate a revised axial load, due to (if restricted movement).
Iterate until convergence, i.e. , from .
Finally, add the bending stress from buckling to .
[1] R. F. Mitchell. Buckling analysis in deviated wells: A practical method. page 13, 1996.
[2] Jonathan Bellarby. Well Completion Design. Elsevier Science, 1 edition, February 2009.