Fulltext
Fulltext
Fulltext
Trine Remmen
Preface
This thesis was carried out during the spring semester of 2018, and completes my Mas-
ter of Science in Petroleum Engineering with a specialization in Drilling Technology.
This thesis is carried out at the Norwegian University of Science and Technology at the
Department of Geoscience and Petroleum.
The thesis was developed with background from the fall project "Analysis of the Indus-
try Standard and Practice for Casing and Tubing Design". The fall project was written
together with Eilen Bauge and Kristian Wilhelmsen.
Buckling is a relevant subject within casing and tubing design, but the available models
are challenging to present. This thesis includes a comprehensive literature study and
two buckling models for installation and operational cycles respectively.
Trondheim, 2018-06-11
Trine Remmen
iii
Acknowledgement
First of all, I would like to thank my supervisor Bjørn Astor Brechan for his guidance
and good discussions. The technical support and feedback provided during this thesis
are gratefully appreciated. Your constructive feedback and enthusiasm have given me
energy throughout this spring semester.
I would also like to pay my gratitude to Bernt Sigve Aadnoy. Thank you for the good
discussions and for always being available when I needed it. Your motivation for the
subject has been truly inspiring.
I would also like to thank Sigbjørn Sangesland. Thank you for providing me with sup-
porting material and for guiding me towards the right path in the beginning. Your
knowledge of the buckling subject has been motivating.
My project partners from fall 2017, Eilen Bauge and Kristian Wilhelmsen, also deserves
my gratitude. I have learned so much from our good discussions, and the comprehen-
sive understanding I got from the fall project gave me a head start on the master thesis.
Finally, I want to thank my fellow students for all the good times through these years at
NTNU. For the friendships, for supporting each other through the tough subjects, and
for all the laughs. These years would never have been the same without you.
T.R.
v
Summary
The development of horizontal drilling technology have enabled the production of oil
and gas from thin layers. With the development of shallow, horizontal subsea wells,
new challenges arises. The occurrence of compressive forces affects horizontal reach
and tubing movement. These challenges can be seen in direct relation to buckling. The
consequences of buckling can be costly work-overs and even loss of integrity.
The integrity of a well construction is essential, and design factors are included in the
stress analysis to ensure safe and profitable hydrocarbon recovery. The highest design
factor is used for axial loads, as they are subject to the highest uncertainty. Exceeding
the buckling limit can introduce local bending forces. The tubing will permanently de-
form, if yield strength is exceeded by triaxial stress due to buckling. Deformation of the
tubing can lead to loss of drift, or even loss of integrity and discharge to the environ-
ment as a worst case scenario. A better understanding of buckling can reduce the risk
of severe axial failures.
Compressive forces are introduced early in the completion process. When the tubing
is lowered, the compressive forces are mainly induced by friction. The occurrence of
compression and geometrical imperfections can buckle the string, and affect the limit
of reach in horizontal wells. A model is therefore developed to predict buckling during
tubing installation. The model is applied on a generic well design for a shallow subsea
field to investigate the tubing integrity for producers and injectors.
Post installation, the tubing is anchored between the wellhead and the production packer.
Cyclic temperature and pressure loads induces tubing movement. Initiation of buck-
ling causes additional tubing movement. A second model is developed to evaluate load
cases and ensure that buckling does not affect the well integrity.
Literature review show that preexisting industry standard buckling models were devel-
oped at a time without widespread use of digital calculation aids. Consequently, the
models are developed by simplifying assumptions enabling use of only hand calcula-
tions. These models have further been implemented in industry software, and are still
in use today, despite these simplifying assumptions.
The digital aids are rapidly developing, and this opens the possibility of a new buckling
model with no simplifying assumptions. When the simulations matches a real scenario,
the model becomes more reliable. This thesis is the first step of achieving this long-term
objective.
vii
Sammendrag
Brønnteknologiske utviklinger har gjort det mulig å bore horisontalt for å produsere olje
og gass fra tynne lag. Med utviklingen av grunne, horisontale undervannsfelt kommer
nye utfordringer. Kompresjonskrefter påvirker både horisontal rekkevidde og rørbeveg-
else. Dette kan sees i direkte relasjon til bøying elle bukling av rør. Konsekvenser av
bukling kan være kostbar brønnservice og tap av brønnintegritet.
Digitale verktøy utvikles stadig, og dette åpner muligheten for å lage en ny bukling-
modell uten forenklede antagelser. Når simuleringene ligner et realistisk scenario, blir
modellen mer pålitelig. Denne avhandlingen er det første steget mot å oppnå dette
langsiktige målet.
Contents
Preface i
Acknowledgement iii
Summary v
Sammendrag vii
List of Figures xi
Nomenclature xv
Abbreviations xviii
1 Introduction 1
1.1 Industry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Literature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Objective and Approach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
3 Literature Study 9
3.1 True Axial Force vs. Effective Axial Force . . . . . . . . . . . . . . . . . . . . 10
3.1.1 Tubing string without packer . . . . . . . . . . . . . . . . . . . . . . . 10
3.1.2 Tubing string with packer . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.2 Critical Buckling Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.2.1 Euler’s equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.2.2 Pitch . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.2.3 Vertical well . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.2.4 Deviated well . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
ix
x CONTENTS
4 Results 33
4.1 Critical Buckling Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.1.1 Effect of relevant parameters . . . . . . . . . . . . . . . . . . . . . . . 34
4.1.2 D/t-ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2 Critical Buckling Ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2.1 Vertical well . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2.2 Deviated well . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2.3 Moment of inertia and radial clearance . . . . . . . . . . . . . . . . . 37
4.3 Buckling During Installation . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4.3.1 Construction of ideal well path . . . . . . . . . . . . . . . . . . . . . . 38
4.3.2 Construction of real well path . . . . . . . . . . . . . . . . . . . . . . . 40
4.3.3 Verification of static hook load . . . . . . . . . . . . . . . . . . . . . . 41
4.3.4 Results from drag analysis . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.3.5 Implementation of critical buckling limits . . . . . . . . . . . . . . . 43
4.3.6 Permanent corkscrewing . . . . . . . . . . . . . . . . . . . . . . . . . 44
4.4 Buckling During Operations . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.4.1 Packer force and neutral stability point . . . . . . . . . . . . . . . . . 46
4.4.2 Buckling in thermal environments . . . . . . . . . . . . . . . . . . . . 48
4.4.3 Set packer procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.4.4 Implementation of critical buckling limits . . . . . . . . . . . . . . . 52
4.4.5 Permanent corkscrewing . . . . . . . . . . . . . . . . . . . . . . . . . 53
6 Conclusion 71
7 Further Work 73
7.1 Buckling During Installation . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
7.2 Buckling During Operations . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
7.3 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Bibliography 75
A Additional Information 77
A.1 Derivation of Pitch . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
A.2 Construction of Build Sections . . . . . . . . . . . . . . . . . . . . . . . . . . 80
A.3 Shut-in . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
A.4 Start bullheading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
3.1 Geometrical imperfections induces bends in the tubing. The inside of the
bend has a smaller area than the outside of the bend. . . . . . . . . . . . . . 11
3.2 True and effective axial force from initial conditions and the squeeze ce-
menting load case presented in Lubinski et al. (1962). The packer config-
uration enables free movement of the tubing. . . . . . . . . . . . . . . . . . 12
3.3 Visualization of helix geometry within a confining cylinder . . . . . . . . . 14
3.4 Illustration of how the capstan effect, buckling contact and tubing weight
contributes to the drag force Bellarby (2009). . . . . . . . . . . . . . . . . . . 18
3.5 Reproduction of the length change and packer force calculation in the
model proposed by Lubinski et al. (1962) . . . . . . . . . . . . . . . . . . . . 23
3.6 Comparison of buckling strain by Lubinski et al. (1962) and Mitchell et al.
(1999). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.7 Comparison of buckling length change by Lubinski et al. (1962) and Mitchell
et al. (1999). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.8 The effective axial force is compared when buckling length change by Mitchell
et al. (1999) is used instead of Lubinski et al. (1962). . . . . . . . . . . . . . . 29
3.9 Triaxial stress plotted towards initial slack-off or pick-up force. This is a
reproduction of the results presented in Lubinski et al. (1962). . . . . . . . 32
xiii
xiv LIST OF FIGURES
5.1 Sinusoidal buckling limit is defined by the Paslay Dawson equation. The
transition between sinusoidal and helical buckling occurs between 1.4
and 2.8 times the Paslay Dawson limit. . . . . . . . . . . . . . . . . . . . . . 57
5.2 Buckling during installation when µ=0.45. The effects of curvature is in-
cluded in the sketch. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.3 Hypothetical scenario of buckling during installation when µ=0.6 . . . . . 59
5.4 Comparison of critical buckling limits for a tubing inside a 6 5/8" casing
and 9 5/8" casing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
5.5 The effective axial force exceeded the critical buckling limits in the vertical
section. No buckling is shown below the vertical section . . . . . . . . . . . 63
5.6 Location of neutral stability point in a shallow subsea well with wet tree . 66
3.1 Critical buckling limits for vertical, straight deviated and curved well sec-
tions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
xv
xvi NOMENCLATURE
Nomenclature
α = Well inclination
β = Coefficient that determines buckling regime
∆ρ i = Change in density inside tubing
∆ρ o = Change in density outside tubing
∆L = Section length in Adnoy’s model
∆L = Total length change by Lubinski
∆L 1 = Length change due to Hooke’s law
∆L 2 = Length change due to helical buckling
∆L 3 = Length change due to radial pressure
∆L 4 = Length change due to thermal effects
∆P i = Change in inner pressure at packer level
∆p i = Change in inner pressure at surface
∆P o = Change in outer pressure at packer level
∆p o = Change in outer pressure at surface
∆T = Change in temperature
δ = Drop of pressure per unit length due to flow
λ = Helix angle
λ0 = Helix angle as a function of length
λ0max = Maximum helix angle as a function of length for sinusoidal buckling
µ = Friction factor
µ = Poisson’s ratio
φ = Azimuth
ρi = Fluid density inside tubing
ρo = Fluid density outside tubing
θ = Dog leg angle
aφ = Azimuth build rate
Ai = Inner area
ai = Inclination build rate
Ao = Outer area
Ap = Packer bore area
As = Cross-sectional area
C = Curvature of the helix
Ct = Coefficient of thermal expansion
E = Young’s modulus
Ef = Final energy in the system
Ei = Initial energy in the system
xvii
xviii NOMENCLATURE
CT Coiled Tubing
MD Measured Depth
TD Total Depth
xxi
Chapter 1
Introduction
The development of horizontal drilling technology have enabled the production of oil
and gas from thin layers. This technology has been practiced on the Norwegian conti-
nental shelf since the 1990’s and especially the Troll field. Some recent fields are shallow
and deviated, which means that the drill pipe during drilling and casing/tubing during
installation will be in compression. Knowing the operational limits for buckling to un-
derstand where failure may occur and to know the limit of reach is important in life
cycle integrity well design.
Frictional drag is a concern when the tubing is lowered in a long-reach horizontal well.
If these drag forces becomes high, buckling can be initiated, and in the worst case cause
lock-up before reaching the desired depth. Consequences of not reaching the desired
depths can be lack of zone isolation. It is evident that the completion phases are very
time-consuming amounting to nearly half of the total well construction time (Aadnoy,
2010, 241). An evaluation of the completion strategy can therefore be valuable for the
company.
Lack of knowledge have caused costly failures in the past (Lubinski et al., 1962, 659). It is
therefore important to understand the theoretical foundation behind design software
programs used in the industry. A better understanding can mitigate potential errors
and wrong simulations.
A frictional drag model have been proposed by Aadnoy et al. (2010). These equations
are applicable for running completion in a well with both inclination and azimuth. The
1
2 CHAPTER 1. INTRODUCTION
drag forces are calculated by dividing the well into a 3D well path, i.e. ideal sections. A
limitation of the model is that it is derived from a soft-string model, i.e. bending stresses
are neglected. The predicted contact force will therefore be underestimated. Another
limitation is that the well is divided into large sections, and the uncertainty arises when
a point of interest lies within a section. A consequence is that dog leg variations within
each section are not accounted for. Dog leg variations will give a higher frictional force
than expected.
Lock-up in coiled tubing is a topic frequently discussed. Lock-up occurs when the ap-
plied surface force exceeds the tubings capacity to move further into the horizontal well.
This analogy is used to explain how compressive forces adds up when a long-reach hor-
izontal well is developed.
1.1 Industry
WellCat is part of the landmark EDM, and it is the most used soft ware tool for advanced
tubing and casing design in the industry. The fall project (Remmen et al., 2017) revolved
around the industry standards and practices for casing and tubing design. Load cases
were evaluated in detail, and the software WellCat was used to analyze the importance
of all input parameters for the casing design of a subsea well in the Norwegian sector.
The remaining challenge is to clarify the theoretical foundation behind the buckling
calculations performed in WellCat. The WellCat manual (Halliburton, 2001) states that
Bellarby (2009) and Mitchell et al. (1999) are the sources for buckling analysis.
The main limitation of full understanding of how the software performs design, is the
lack of technical details in the engineering performed. It is therefore challenging to
understand the theoretical basis behind the calculations.
1.2 Literature
A comprehensive literature study on buckling was conducted, and the most acknowl-
edged theory and methods are presented.
Lubinski et al. (1962) provided a model for helical buckling for tubing sealed in pack-
ers. The model is mainly developed for a vertical well with no friction. Mitchell et al.
(1999) developed a theoretical understanding of the previous mentioned model, and
developed equations for application of the model in deviated wells. It is challenging to
understand and reproduce the Lubinski et al. (1962) model. Some of the assumptions
taken does not apply for modern wells, and modifications needs to be done.
The main limitations of the theoretical foundations is the basic assumptions. Uncer-
tainty arises when the assumptions does not match the realistic scenario.
1.3. OBJECTIVE AND APPROACH 3
Axial forces are subject to the greatest uncertainty. Since they are dependent on pres-
sure, temperature, bending and drag. The consequences associated with axial failure
can be serious, whether it is a tensile or compressive failure. As axial failure could cause
considerable damage on people and surface equipment, a relatively high design factor
is used (Bellarby, 2009, 521-522).
This thesis can contribute to better understanding of when and where buckling occurs,
and the importance of a thorough investigation of this matter throughout the entire
lifetime of a completion. By a better understanding of the load limits, costly work-overs
and time consuming interventions can be avoided. Numerical aids improves the drag
analysis by dividing the tubing into indefinite small cross-sections.
Two models were developed in MatLab. The first model is a representation of a scenario
where a tubing is installed in a shallow, horizontal well. Buckling induced by drag forces
are presented by implementing buckling limits. Lock-up scenarios are then predicted.
The second model is an implementation of the buckling analysis used in WellCat. Both
models are generalized to apply for both vertical and deviated wells, and results for a
shallow, horizontal well are presented.
Lubinski et al. (1950) and Wu et al. (1995) included weight in the derivation of the criti-
cal sinusoidal and helical buckling limit, respectively, in a vertical well.
5
6 CHAPTER 2. STATE OF THE ART
Mitchell et al. (1999) found that a helical buckling limit could be derived by multiplying
the limit developed by Dawson et al. (1984) with a constant. The region between these
limits was defined as a transition where either sinusoidal or helical buckling can occur.
Dawson et al. (1984) neglected well curvature in the derivation of the critical buckling
limit.
He et al. (1995) showed that the established equation by Dawson et al. (1984) was con-
servative in curved sections. It was therefore developed a new critical buckling limit, in-
cluding the curvature effects. Their model indicated that curvature stabilizes the pipe,
hence resulting in a higher critical buckling limit than previously expected.
Aasen et al. (2002) did a thorough investigation of the available models, to find that
many buckling limits only varied by the coefficients.
When Lubinski et al. (1962) presented the model for helical buckling in a vertical well,
Mitchell et al. (1988) investigated the foundation of the model. The solution was put in
a more theoretical context to make the model applicable for deviated wells.
The model presented by both Lubinski et al. (1962) and Mitchell et al. (1988) was based
on the assumption that the pipe instantly took a helical shape when reaching the crit-
ical buckling limit. It was however recognized by Mitchell et al. (1997) that sinusoidal
buckling could occur before helical buckling, hence a new theoretical model was devel-
oped.
A fourth order differential equation for the helix angle was derived by Mitchell et al.
(1997). The solution of the differential equation was found to be very useful when pre-
dicting the effect of well deviation on buckling shape, tubing length change, contact
force and bending stress.
Mitchell et al. (1999) presented a more practical applicable model based on the equa-
tions in Mitchell et al. (1997). Correlations were developed to match the complex mod-
els. The correlations was compared to the results (Lubinski et al., 1962) presented, and
it was shown that the buckling length change depended on the buckling shape. A con-
vincing discovery was that the result presented by Lubinski et al. (1962) predicted an
increase in tubing length for an increase in deviation angle, while Mitchell et al. (1999)
predicted the opposite. Lateral gravity forces is proven to prevent buckling, and it was
2.3. BUCKLING DURING INSTALLATION 7
concluded that Mitchell et al. (1999) predicted a more realistic scenario. Another find-
ing was that (Lubinski et al., 1962) over predicted tubing movement. For a fixed tubing,
the axial loads would by this method be underestimated, and this could lead to a non-
conservative design. The model presented by (Mitchell et al., 1999) was proven to be
better when designing a tubing in a deviated well.
Mitchell et al. (1999) gave a better understanding of the consequences of helical buck-
ling versus sinusoidal buckling in a deviated well. The bending stress, buckling strain
and length change is depending on whether the buckling is sinusoidal or helical shaped.
The contact force is the last main factor affected by the buckling shape. When the buck-
ling shape changes from sinusoidal to helical, the contact force increases significantly.
Aadnoy et al. (2010) presented a three-dimensional friction model for the drag analy-
sis. This model includes contact force, dogleg angle and the capstan effect. The string
was considered to be either tension- or weight-dominated, depending on whether the
wellbore section is straight or curved. This assumption simplified the drag equations
significantly, and a set of only two equations could be used to calculate drag in a well-
bore with both inclination and azimuth. One can further use this drag force to calculate
the hook load on surface and the axial load throughout the entire completion.
8 CHAPTER 2. STATE OF THE ART
Literature Study
A study of critical buckling limits for vertical and deviated well sections has been con-
ducted. The effective axial force is defined, and its importance in buckling analysis is
justified.
2. Buckling during operations (i.e. production, injection and other load cases)
The model of buckling during operations is the existing model used for buckling analy-
sis in the industry leading software WellCat. This model is developed by Lubinski et al.
(1962) and Mitchell et al. (1999).
The aim of this chapter is to obtain a better understanding of failures related to sinu-
soidal or helical buckling, and to prevent potential failures in the future.
9
10 CHAPTER 3. LITERATURE STUDY
The inner and outer pressure will affect the structural stability, so one can not rely on the
compressive forces alone when buckling is calculated. The axial force combined with
inner and outer pressure is summed up to make buckling calculation more convenient.
The resulting force is called the effective axial force. When this effective force exceeds
a certain limit, buckling will occur. The effective axial force is not a true force that can
be measured, but it is useful to predict buckling tendencies. The effective axial force is
defined as:
Fe f f = F a − p i A i + p o A o (3.3)
These pressure terms are evaluated along the string.
more unstable. At a certain point it can snap into a buckled shape, even though the
entire pipe is in tension.
Figure 3.1: Geometrical imperfections induces bends in the tubing. The inside of the
bend has a smaller area than the outside of the bend.
Compression and internal pressure will therefore promote buckling, while tension and
external pressure reduces the chance of buckling (Bellarby, 2009, 492). Also note that a
sufficiently high internal pressure can cause the effective axial load to be negative, al-
though the true axial force is positive.
The hydrostatic pressures acting on the bottom of the tubing will now act on cross-
sections of both the tubing and the packer. The true axial force are presented by Lubin-
ski et al. (1962) as:
F a = (A p − A o )P o − (A p − A i )P i (3.4)
To calculate the axial force at surface, the string weight is simply added to equation
(3.4).
Lubinski et al. (1962) introduced the term "fictitious force" as the "effective force" at the
bottom of the string. The equation for fictitious force is:
F f = A p (P o − P i ) (3.5)
Equation (3.5) is obtained by inserting equation (3.4) in to (3.3). If the inner and outer
area is carefully expressed as the cross-sectional area in which the pressures are acting,
12 CHAPTER 3. LITERATURE STUDY
equation (3.1) and (3.3) can be correctly used for any general case. The string weight
must be included to find the effective force at surface.
When the tubing is fixed, and loads affect tubing movement, the packer will act with an
opposing force. This is presented in detail in section 3.4.
The true and effective axial force is plotted towards depth in figure 3.2. Note that the
effective axial force equals zero at the bottom of the string. The reason for this is that
the pressure terms cancels each other out at this point. This can be easily seen if the
equation for true axial force is inserted in the equation for effective axial force.
Where the true axial force crosses the z-axis, the string is neither in compression nor
tension. This is the point of zero axial load, and it is frequently denoted as the neutral
point. The neutral point will not give us useful information on whether the string buck-
les or not. Buckling analysis requires knowledge of the effective axial force. This point
where the effective axial force crosses the z-axis is defined as the neutral stability point
(NSP).
Fa Initial Conditions
1000
Feff Initial Conditions
2000 Fa Load Case
Feff Load Case
3000 NSP
4000
Depth [ft]
5000
6000
7000
8000
9000
10000
-8 -6 -4 -2 0 2 4 6
Force [lb] ×10 4
Figure 3.2: True and effective axial force from initial conditions and the squeeze ce-
menting load case presented in Lubinski et al. (1962). The packer configuration enables
free movement of the tubing.
Bellarby (2009) defines the neutral stability point as the boundary where buckling can-
not occur and where it may occur, and one can therefore use the effective buckling
force to evaluate at which depth the string buckles. This method is convenient, but
utterly conservative. It is therefore developed critical buckling forces to make realistic
predictions of buckling initiation.
3.2. CRITICAL BUCKLING LIMIT 13
4. The length of the column is very large compared to the other dimensions
Equation (3.6) is a solution of this differential equation, and indicates the buckling limit,
F cr :
π2 E I
F cr = n (3.6)
L2
where E I is the bending stiffness, and L is the tubular length. The coefficient, n, is de-
pendent on the end conditions. For fixed-free end conditions, n = 41 . For fixed-fixed
end conditions, n = 4 (Hearn, 1997, 28).
Euler’s equation is fundamental for understanding elastic stability theory, but it does
not consider tubulars confined within another circular cylinder. Laterally constrained
pipes also presents more complex characteristics, and the equation has to be modified
for application in well configurations (Cunha et al., 2003).
3.2.2 Pitch
Lubinski et al. (1962) derived an expression for the pitch, P , of a helix within a confining
cylinder. The assumption of a vertical suspended weightless pipe was the basis of this
theory. A helix geometry is displayed in figure 3.3.
14 CHAPTER 3. LITERATURE STUDY
where the force, F f , is the fictitious force (Lubinski et al., 1962), also known as effective
axial force at the bottom of the string. Equation (3.7) describes the pitch at the very
bottom of the tubing string.
Buckling limits were also derived by considering tubular weight. Lubinski et al. (1950)
derived the following solution for sinusoidal buckling in a vertical well:
Wu et al. (1995) presented the following solution for helical buckling in a vertical well:
Aasen et al. (2002) revisited several buckling models, and came to the conclusion that
the different sources had presented similar coefficients to determine the critical buck-
ling limit for vertical wells. The practical relevance of these coefficients are limited, as
the buckling limit is generally low for vertical wells.
3.2. CRITICAL BUCKLING LIMIT 15
This equation is called the Paslay-Dawson equation. The equation is based on the as-
sumption that the radial clearance, r c , is small in comparison with the pitch length, P .
That is, 2πr c
P << 1. Another assumption is that the well curvature is negligibly small.
Hence, the equation applies for straight, inclined wellbores.
Mitchell et al. (1999) used the Paslay-Dawson equation to define the critical limits and
the transition between sinusoidal and helical buckling. The buckling force must exceed
the Paslay-Dawson force to initiate sinusoidal buckling. The following equation defines
a region where the buckling may be either sinusoidal or helical:
s
E I w b sinα
F cr = 1.4 ∼ 2.8 × 2 (3.11)
rc
The variation between 1.4 and 2.8 reflects uncertainty of when sinusoidal buckling switch
to helical buckling while loading, and visa versa while unloading. The buckling force
must exceed 1.4 to initiate helical buckling, and when it exceeds 2.8, only helical buck-
ling occurs.
βF 1nbc E I
F i cr = (3.12)
rc
where F i cr is the improved critical buckling force and F 1nbc is the normal contact force
including curvature effects. β is a coefficient that determines whether the buckling
regime becomes sinusoidal or helical. β=4 is used for sinusoidal buckling and β=8 for
helical buckling.
q
F 1nbc = (wsinα + F a a α )2 + (F a sinαa φ )2 (3.13)
where F a is the axial force, a α is the inclination build rate and a φ is the azimuth build
rate.
Table 3.1: Critical buckling limits for vertical, straight deviated and curved well sections.
Dawson et al. (1984) neglected the curvature effects when deriving the equation for si-
nusoidal buckling in deviated wells. As curvature effects will stabilize the tubing, the
critical buckling limit by Lubinski et al. (1962) is conservative in curved sections.
He et al. (1995) show that curvature effects stabilizes the tubing, but the results were
only conducted from one small-scale experiment.
3.3. BUCKLING DURING INSTALLATION 17
The hook load is the actual weight of a pipe measured at surface. Generally, this force
includes the weight of the string in air, reduced by any force that tend to reduce that
weight. The most commonly known opposing force to the hook load, is the buoyancy
force. Another opposing force is friction, also known as the drag force.
F st at i c is the hook load when no movement is induced and F d ynami c is the hook load
after induced movement. The drag force will be negative when lowering the string, and
positive when hoisting the string.
Bellarby (2009) defined tubing-to-casing drag as a contact force that derives from three
main sources:
1) Forces due to gravity. This will be the normal force, which is the component of the
tubing weight that acts perpendicular to the casing. The normal force will increase as
the inclination increases.
2) Frictional forces from buckling. When the tubing takes the form of a sinusoidal or
helical buckling shape, it will immediately contact the casing. If buckling is initiated, it
causes an increase in locally generated forces.
3) Frictional forces from the capstan effect. This effect occurs when the tubing passes
through doglegs. If the tubing is in tension, it will be pulled to the high side of the bend.
If the tubing is in compression, it will be pushed to the low side of the bend.
Figure 3.4 shows how all three contributions to the contact force can occur at the same
time. In this figure, the neutral stability point is located below the first build section.
18 CHAPTER 3. LITERATURE STUDY
Figure 3.4: Illustration of how the capstan effect, buckling contact and tubing weight
contributes to the drag force Bellarby (2009).
If the drag force becomes sufficiently high, it causes tubing lock-up when running on
hole. If the compressive forces exceeds the buckling limit, the contact force induced
by buckling will be added. The drag forces are especially important for extended-reach
wells with high inclination.
The drag analysis can be modeled by only two sets of equations, whereas they each rep-
resent the gravitational and tensional friction effects. The first equation is developed
for straight, inclined wellbore sections:
³ ´
F 2 = F 1 + ∆Lw cosα − µsinα (3.16)
where α is the inclination, φ is the azimuth and θ is the dog leg angle in radians. The
capstan effect is derived from the assumption that the pipe is weightless, hence weight
is added at the end of the bend.
When the string is located in a straight section, it is assumed that the friction contribu-
tion solely comes from the gravitational normal force. It is therefore said that equation
(3.16) is "weight-dominated". When the string is located in a curved section, it is as-
sumed that the tension is much higher than the weight of the pipe in the bend. The
string is assumed to be weightless. Equation (3.17) is therefore said to be "tension-
dominated".
The purpose of these assumptions was to simplify the equations. The equations are eas-
ier to use, and easier to understand. When tension or compression dominates, equation
(3.17) is correct. However, the "tension-dominated" equation can lead to errors at the
parts of the string where tension and compression does not dominate (Aadnoy et al.,
2010, 34). In a highly inclined build section, the weight can be the dominating effect.
The exact equation in App. E in Aadnoy et al. (2010) should therefore be applied.
The end-condition F 1 is assumed to be zero. This is correct if the hole is clean, and no
initial force pushes the string backwards when it is run in hole.
20 CHAPTER 3. LITERATURE STUDY
Lubinski et al. (1962) presented compression as a positive force and tension as a nega-
tive force. The opposite is more commonly used. In the following section, compressive
is considered negative, while tension is considered positive.
Lubinski et al. (1962) used the notation F a for the direct piston force acting on the tub-
ing, and F f for the fictitious piston force that promotes buckling. This implies that F a
is the true axial force and F f is the effective axial force at the bottom of the string. The
fictitious force F f is replaced by F e f f for consistency.
A freely hanging tubing only supported by pressure was imagined. Pressure acted on
the cross-sectional area of both the tubing and the packer. This opposing force caused
the tubing to shorten. This physical shortening in the tubing was expressed by Hooke’s
law:
Lp
∆L 1 = F (3.19)
E As
where L p is the tubing length, and F is the true axial force. When the outer and inner
pressure changed, the length change was obtained:
Lp h i
∆L 1 = (A p − A o )∆P o − (A p − A i )∆P i (3.20)
E As
It is assumed that helical buckling occurs when F e f f < 0. There will be a physical length
change due to the buckling itself, and this length change will always be negative, as the
3.4. BUCKLING DURING OPERATIONS 21
tubing shortens when it snaps into a helical shape. The length change is promoted by
the effective force acting on the bottom of the tubing. Any length change due to helical
buckling is calculated by a non-linear equation:
r c2
∆L 2 = − F e2f f (3.21)
8E I w
where the weight, w, was defined as:
w = wl + wi − wo
(3.22)
= wl + ρi Ai − ρo Ao
and ρ i is the fluid density in the tubing, and ρ o is the fluid density in the annulus. w s ,
w i and w o is the weight of tubing in air, weight of liquid inside tubing and weight of
liquid outside tubing respectively. The unit for w is [lb/in]. Equation (3.22) equals the
equation for buoyancy factor when two different fluids are present on the inside and
the outside of the pipe, times the weight of string in air (Aadnoy et al., 2010, 25).
The effective force, F e f f , was calculated with equation (3.5). When the outer and inner
pressures changed, the length change due to buckling was expressed as:
r c2 h i
∆L 2 = − A 2p (∆P o − ∆P i )2 (3.23)
8E I w
Equation (3.23) is valid under the assumption that the initial inner pressure equals the
initial outer pressures.
2 1+2µ
µ ∆ρ t − R ∆ρ a − 2µ δ 2 2µ (∆p i − R 2 ∆p o )L p
∆L 3 = − Lp − (3.24)
E R2 − 1 E R2 − 1
where µ is Poisson’s ratio. R is the ratio between OD and ID of the tubing, ∆ρ is the
change in density from initial conditions to load case conditions and δ is the drop of
pressure per unit length due to flow.
Fluid drag induced axial loads are often ignored in hand calculations and in many soft-
ware packages (Bellarby, 2009, 489-490). For ballooning in static conditions, the follow-
ing equation can be used:
−2µL
∆L 3 = (∆p i A i − ∆p o A o ) (3.25)
E (A o − A i )
In equation (3.25), ∆p is the average change in pressure between surface and packer.
22 CHAPTER 3. LITERATURE STUDY
A positive temperature change will cause thermal expansion, and a negative temper-
ature change will cause thermal contraction. The length change due to temperature
changes is expressed by (3.26) by both Lubinski et al. (1962) and Bellarby (2009).
∆L 4 = L p C t ∆T (3.26)
The total operational length change, ∆L will then be the sum of these:
∆L = ∆L 1 + ∆L 2 + ∆L 3 + ∆L 4 (3.27)
The freely suspended tubing will now have changed its physical length by ∆L.
The same tubing was then configured with a packer that permitted limited motion. The
tubing can be installed in compression on purpose. Lubinski et al. (1962) explained that
a shoulder on the tubing is removed, so the tubing can be elongated with the length
corresponding to this compressive force, F so . The length elongation is ∆L 5 . This expla-
nation is confusing because it does not apply for a practical well installation.
∆L = ∆L 1 + ∆L 2 + ∆L 5 = 0
³ ´
∆L 5 = − ∆L 1 + ∆L 2
(3.28)
³ Lp r c2 2
´
=− F so − F so
E As 8E I w
When the tubing is subjected to loading, the total operational length change will be:
∆L = ∆L 1 + ∆L 2 + ∆L 3 + ∆L 4 + ∆L 5 (3.29)
This length change will apply to the limited motion of the tubing, and the seal length
must be within this ∆L.
The final scenario was to imagine a tubing configuration where the packer permits no
motion. There will then be no physical length changes in the tubing, because it can not
physically move from its position. Instead of physically change its length, the packer
3.4. BUCKLING DURING OPERATIONS 23
The effective force was earlier defined as the force defining the neutral stability point.
The fictitious force is consecutive. The fictitious force was plotted towards length change,
where the negative part of F indicated buckling, while the positive part of F indicated
no buckling.
Lubinski et al. (1962) simulated a squeeze cementing load case. The effective axial force
at packer depth was calculated and plotted as if the packer permitted free motion. The
case considered however, was a packer permitting no tubular movement. The force
associated with the calculated length changes would therefore be transferred to the
packer. The packer force were shown to be tensile, thus it reduced the compressive
effective force at the packer depth.
F e∗f f = F e f f + F p (3.30)
The corresponding axial force for an anchored tubing is:
F a∗ = F a + F p (3.31)
50
0
dL [in]
-50
-100
-150
-200
-250
-8 -6 -4 -2 0 2 4 6 8
Axial Force [lb] ×10 4
Figure 3.5: Reproduction of the length change and packer force calculation in the model
proposed by Lubinski et al. (1962)
A new coordinate system was presented with F e f f as the origin. The shortening of the
tubing ∆L was plotted, and it was clear that the packer force was located at the positive
side of the x-axis in the new coordinate system. A flow chart is presented for further
understanding of this method:
24 CHAPTER 3. LITERATURE STUDY
dL'
Hooke's Law
Fso
Fa
Plot Ff and read of the
length on the y-axis
Sum up the total length
change, dL
Helical Buckling
dL'
Ff
F
Ff
3.4. BUCKLING DURING OPERATIONS 25
dL'
Ff Ff*
F
dL
Fp
Thermal Expansion
Add Fp to Fa and Ff to
obtain Fa* and Ff*
u 1 = r c cos(λ)
(3.32)
u 2 = r c sin(λ)
where the helix angle, λ, is given by a differential equation:
h i ³ ´
− E I λ0000 + 2E I (λ0 )3 − F λ0 + w l /r c sin(λ) = 0 (3.33)
F is the value of the axial force in the buckled state, also known as the effective axial
force or fictitious force.
According to Mitchell et al. (1999), general solutions to equation (3.33) are difficult to
obtain. The Galerkin method with cubic interpolation was used to calculate λ as a func-
tion of z. This method is too complex and slow for use in casing and tubing design, so
Mitchell et al. (1999) presented approximate correlations to make the model more prac-
tical.
Mitchell et al. (1999) presented correlations for the non-linear differential equations de-
veloped for buckling in deviated wells. By identifying that there can be both sinusoidal
and helical buckling regime, parameters like helix angle, bending stress and bending
strain were conditionally developed.
When a string is buckled like a sinusoidal, the bend of the helix, λ, will change through
the S-shape. Thus, the maximum helix angle was found for sinusoidal buckling. The
maximum helix angle, λmax , for sinusoidal buckling was given by Mitchell et al. (1999)
as:
1.1227 0.04
λ0max = p F e f f (F e f f − F cr )0.46 (3.34)
2E I
For helical buckling, the helix angle is assumed to be constant (Bellarby, 2009, 496). The
helix angle was given as:
s
Fe f f
λ0 = (3.35)
2E I
3.4. BUCKLING DURING OPERATIONS 27
The buckling strain, ²b , was defined as the buckling length change per unit length, thus
the elongation is unitless. Lubinski et al. (1962) followed the weightless string principle
and derived an expression for the ², using the geometry of a helix configuration (Figure
3.3). For the case of sinusoidal buckling, Mitchell et al. (1999) developed a correlation
of the average strain. It was assumed that sinusoidal buckling occurs up to 2.8 × F cr .
The buckling strain was conditionally determined:
0
Fe f f < F p
´0.92
r2
³
²b = −0.7285 4Ec I F e0.08
ff
F ef f − F p F p < F e f f < 2.8F p (3.36)
r c2
− 4E I F e f f F e f f > 2.8F p
In figure 3.6, it is clear that the amount of buckling strain is significantly reduced for
values of F below the helical buckling limit.
Lubinski
3.5 Mitchell Correlation
2.5
4EIe/r2
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
F/Fcr
Figure 3.6: Comparison of buckling strain by Lubinski et al. (1962) and Mitchell et al.
(1999).
To find the buckling length change, ∆L b (= ∆L 2 ), the buckling strain was integrated with
respect to z. z was integrated between the bottom of the string and the neutral stability
point. It was assumed that the neutral stability point was located within the string.
Z n
∆L z = ²z d z (3.37)
0
Further,
F z Fe f f
= (3.38)
z n
28 CHAPTER 3. LITERATURE STUDY
Inserting the equation for helical buckling strain in equation (3.37), and integrating
from 0 to n gives:
r2
Z n
∆L z = − c Fz d z
0 4E I
h r2 in
= − c F z2 (3.39)
8E I 0
r2
= − c nF e f f
8E I
The neutral stability point, n, can be calculated by:
Fe f f
n= (3.40)
w
where w has to be defined as in equation 3.22. Inserting equation (3.40) for (3.39), gives
the final equation for sinusoidal and helical buckling length changes. The same proce-
dure is performed with the correlations for buckling strain, and the conditional buck-
ling length changes are given as:
0 Fe f f < F p
r c2
∆L b = − 4E I (F e f f − F cr )[0.3771F e f f − 0.3668F cr ] F p < F e f f < 2.8F p (3.41)
r c2 2
− 8E I w F e f f F e f f > 2.8F p
The length change correlation for sinusoidal and helical buckling is presented in figure
3.7.
Lubinski
14 Mitchell Correlation
12
10
-8EIwdLb/r2
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
F/Fcr
Figure 3.7: Comparison of buckling length change by Lubinski et al. (1962) and Mitchell
et al. (1999).
Applying the buckling length change by Mitchell et al. (1999) into the Lubinski et al.
(1962) model, the packer force reduces slightly. The result is seen in figure 3.8.
3.4. BUCKLING DURING OPERATIONS 29
0
dL [in]
-50
-100
-150
-200
-250
-8 -6 -4 -2 0 2 4 6 8
Figure 3.8: The effective axial force is compared when buckling length change by
Mitchell et al. (1999) is used instead of Lubinski et al. (1962).
The buckling severity depends on whether the tubing is sinusoidal or helically buckled.
The pitch will also vary along the string. When buckling first is initiated, the pitch is
expected to be high. The more severely the tubing buckled, the smaller the pitch gets.
The dog leg curvature of the helix in [rad/in] was given as:
(
r c λ02
max Sinusoidal Buckling
κ= (3.42)
rc λ 02
Helical Buckling
The higher dog leg curvature, the higher bending moment. Bending moment, M , was
generally given as:
M = EIκ (3.43)
The bending moment was further given for the different buckling regimes:
0 F e f f < F cr
³ ´0.92
M = 0.6302r c F 0.08 F e f f − F cr F cr < F e f f < 2.8F cr (3.44)
ef f
0.5000r c F e f f F e f f > 2.8F cr
M × OD
σb = (3.45)
2I
30 CHAPTER 3. LITERATURE STUDY
Conditionally:
0
F e f f < F cr
³ ´0.92
σb = 0.3151 ODr
I
c
F e0.08 F e f f − F cr F cr < F e f f < 2.8F cr (3.46)
ff
0.2500 ODr
c
Fe f f F e f f > 2.8F cr
I
Lubinski et al. (1962) assumes that the tubing buckles when the tubing is subjected to
an effective axial force. This is a conservative buckling limit.
When applying the change in effective force in equation (3.23), it is assumed that the
initial inner pressure equals the initial outer pressure. Equation (3.23) has to be mod-
ified to apply for cases where the initial inner pressure does not equal the initial outer
pressure.
Another assumption is that the neutral stability point is within the string. This is of-
ten true for vertical wells, but for multiple completions or shallow, horizontal wells, the
neutral stability point can be shifted to a location above the string of interest (Lubinski
et al., 1962).
It was challenging to compare the results from Lubinski et al. (1962) and Mitchell et al.
(1999), as the entire example in Lubinski et al. (1962) is developed to work for a verti-
cal well. To compare the results, Mitchell et al. (1999) assumed a kick-off point (KOP)
at 2000 ft, and constructed a well path consisting of one vertical section and one sail
section. This is not a realistic scenario, and considerations regarding build-sections is
therefore not included.
Mitchell et al. (1999) assumed sinusoidal buckling for F cr < F e f f < 2.8F cr . It is an un-
certainty of whether helical or sinusoidal buckling occurs between 1.4F cr and 2.8F cr .
In theory, the maximum limit of stable sinusoidal buckling occurs at 2.8F cr . In prac-
tice however, helical buckling will occur before reaching that limit. This is because of
irregularities in the actual geometry. This assumption can lead to under-estimation of
buckling length changes.
3.5. PERMANENT CORKSCREWING 31
As an industry practice, it is established that casing and tubing strings are limited to the
elastic stress-strain regime. For thick-walled cylinders, Lame’s equations can be applied
for the radial and tangential stress components, σr and σθ respectively, as a function of
pressure (Bellarby, 2009, 515):
r i2 r o2 (p o − p i ) 1 p i r i2 − p o r o2
σr = +
r o2 − r i2 r2 r o2 − r i2
(3.48)
r 2 r o2 (p o − p i ) 1 p i r i2 − p o r o2
σθ = − i 2 +
r o − r i2 r2 r o2 − r i2
When Lame’s equations are inserted in the von Mises criterion, one can express one
relation for permanent corkscrewing at the outer wall, and another relation for perma-
nent corkscrewing and the inner wall.
σb
σz = σa ± rc (3.51)
ro
The axial stress, σa , is the total axial load divided on cross-sectional area, A s .
The bending stress, σb , is given depending on the buckling type and well deviation. The
sign in front of the bending stress, ±, is chosen as the one that gives the highest σo or σi .
Lubinski et al. (1962) used bending stress for helical buckling exclusively when calculat-
ing how permanent corkscrewing relates the the initial slack-off or pick-up force, and
figure 3.9 is a reproduction of this:
Sigmao Initial
12 Sigmao
Sigmai
10
Triaxial Stress [psi]
0
-8 -6 -4 -2 0 2 4 6 8
Initial Slack-Off and Pick-Up Force [lb] ×10 4
Figure 3.9: Triaxial stress plotted towards initial slack-off or pick-up force. This is a
reproduction of the results presented in Lubinski et al. (1962).
σo,i ni t i al is the initial yield after slack-off but before any pressures or temperatures have
changed. The dotted lines represents the available grades for this tubing. When the
triaxial stress exceeds a line, the line above should be chosen, and the available grade
can be read of at the y-axis.
Chapter 4
Results
Tubing design is determined by analyzing load cases and investigating critical parame-
ters. In this chapter, the critical parameters are clarified and evaluated, a drag analysis
is performed to predict lock-up situations, and operational load cases are run to calcu-
late the packer force and to locate the neutral stability point. Both models investigated
whether buckling affects the well integrity by permanent corkscrewing.
The theories are tested on a shallow, horizontal well path with a dry x-mas tree. Figure
4.1 presents a simple sketch. The purpose of using this design, is to observe how pres-
sure and temperature affects the axial force in a tubing installed in a horizontal well.
33
34 CHAPTER 4. RESULTS
A large radial clearance implies that the tubing is small relative to the confining casing.
The tubing will tend to lean towards the sides of the casing to reach the equilibrium po-
sition. When the tubing has significant radial space, the buckling severity will increase
due to reduced stability. The critical buckling limit for a vertical well is independent of
radial clearance. Radial clearance affects the critical buckling limit for deviated wells. A
well inclination of 45° is used. Figure 4.2a shows a correlation between radial clearance
and the critical buckling limit. The opposite effect will be seen for bending stiffness and
weight. A higher weight implies a larger tubing thickness. If the thickness increases,
bending stiffness increases accordingly because of the definition of moment of inertia.
It is evident that the trend will be positive for increasing values of bending stiffness and
weight.
Increased inclination will provide support and stability to the tubing. Figure 4.2b shows
how the critical buckling limit increases with inclination. A horizontal section will pro-
vide more stability than a section with a smaller inclination. Note that these are strictly
a simulation of sail sections with different inclinations. As equation 3.10 does not in-
clude curvature effects, the plot will not apply for curves sections. For a vertical well,
i.e. no inclinations, the limits presented in equation (3.8) and (3.9) should be used.
×10
4 Effect of radial clearance on critical buckling limit ×10
4 Effect of inclination on critical buckling limit
9 4.5
6 3
5 2.5
4 2
3 1.5
2 1
1 0.5
0 0
0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 0 10 20 30 40 50 60 70 80 90
Radial clearance [in] Inclination [deg]
Figure 4.2: Effect of radial clearance and inclination on critical buckling limit in a devi-
ated well section with 45°inclination.
4.1. CRITICAL BUCKLING LIMIT 35
4.1.2 D/t-ratio
The outer diameter to thickness ratio, D/t -ratio, will affect the geometric instability of
the tubing. If this ratio is high, i.e. a thin-walled pipe, the instability naturally increases.
One of the failure modes in casing and tubing design is collapse. Buckling is previously
stated to be suppressed by high external pressure. However, if the tubing is thin-walled,
it is showed that the tubing fails because of buckling before it fails because of collapse.
One can categorize pipes with different D/t -ratios, in three collapse regimes. The first
one is for low D/t -ratios. This is called yield collapse, and it is here Lame’s equations are
valid. The second region is plastic collapse, where the collapse limit is adapted to API
empirical data. The last one, but most relevant to the buckling scenario, is the elastic
collapse. Pipes within this D/t -category seems to collapse at a trend matching Euler’s
buckling equation. (Remmen et al., 2017)
Figure 4.3a presents a range of tubings with low D/t -ratio. The figure shows no correla-
tion between D/t -ratio and the critical buckling limit. Only L80 pipes were chosen, and
the D/t -ratios were calculated to be within yield or plastic collapse regime (Remmen
et al., 2017). In figure 4.3b it is clearly a correlation between high D/t -ratios and the
critical buckling limit. This result is expected for tubings in the elastic collapse region,
where the thin-walled piped experiences buckling before collapse.
×10
4 The effect of D/t-ratio on critical buckling limit ×10
6 The effect of D/t-ratio on critical buckling limit
10 2.2
1.8
7
1.6
6
5
1.4
4
1.2
1
2
1 0.8
9 10 11 12 13 14 15 16 17 18 20 25 30 35 40 45 50
D/t-ratio D/t-ratio
Figure 4.3: Effect of D/t-ratio on the critical buckling limit in a deviated well section
with 45°inclination and a horizontal well.
36 CHAPTER 4. RESULTS
where the notation i is the tubing initially chosen. The variables in equation (4.1) is
the moment of inertia, I , and the buoyed weight, w b . Generally, one can increase the
tubing weight to increase the critical buckling limit. One can also increase moment of
inertia by increasing tubing OD, and even the lowest corresponding weight will increase
the limit.
Note that the buckling ratio for a deviated well is independent of inclination. It can
therefore be applied for both inclined and horizontal sections. The radial clearance will
affect the buckling limit for deviated and horizontal wells. In an inclined or horizontal
well, the radial clearance affects the critical buckling limit.
For the 4" #10.7 tubing marked with green, it can be seen that this tubing is a better
choice than the 3 1/2" #12.7 tubing. The lower radial clearance in the deviated section
will stabilize the tubing. The critical buckling limit for vertical wells is not depending
on the radial clearance, hence other factors like bending stiffness and weight will be
dominating. The results are presented in table 4.1.
4.2. CRITICAL BUCKLING RATIO 37
Table 4.1: A set of tubing configurations from 2 7/8" to 4 1/2" is compared to the 3 1/2"
#9.2 tubing. The tubings are sorted from best to worst with regards to buckling.
VERTICAL INCLINED/HORIZONTAL
OD [in] w [lb/ft] rc [in] Fcr* OD [in] w [lb/ft] rc [in] Fcr*
4 1/2" 12.6 0.7105 1.64 4 1/2" 12.6 0.7105 2.34
3 1/2" 12.7 1.2105 1.36 4" 10.7 0.9605 1.52
4" 10.7 0.9605 1.29 3 1/2" 12.7 1.2105 1.35
4" 9.5 0.9605 1.14 4" 9.5 0.9605 1.35
3 1/2" 10.2 1.2105 1.11 3 1/2" 10.2 1.2105 1.11
3 1/2" 9.2 1.2105 1.00 3 1/2" 9.2 1.2105 1.00
3 1/2" 7.7 1.2105 0.85 3 1/2" 7.7 1.2105 0.86
2 7/8" 8.6 1.5230 0.81 2 7/8" 8.6 1.5230 0.67
2 7/8" 6.4 1.5230 0.61 2 7/8" 6.4 1.5230 0.51
Figure 4.4b presents how radial clearance affects the critical buckling ratio. Even though
the critical buckling ratio in the vertical section is not directly proportional to the radial
clearance, a correlation is seen. Radial clearance is calculated from the outer diameter,
and the outer diameter will strongly affect the moment of inertia. Hence, a correlation
will still be present here.
Effect of I on critical buckling ratio Effect of radial clearance on critical buckling ratio
2.4 2.4
1.8 1.8
1.6 1.6
Fcr*
Fcr*
1.4 1.4
1.2 1.2
1 1
0.8 0.8
0.6 0.6
0.4 0.4
1 2 3 4 5 6 7 8 9 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6
Moment of inertia [in4] Radial clearance [in]
Figure 4.4: Effect of moment of inertia and radial clearance. The markers above "1" in-
dicates the range of tubings that can withstand buckling better than the initially chosen
tubing.
38 CHAPTER 4. RESULTS
The installation of a completion string in a long horizontal well is often one of the most
critical operations. The purpose of investigating this subject is to predict lock-up situa-
tions and permanent corkscrewing due to helical buckling.
The friction model by Aadnoy et al. (2010) was developed for two well configurations.
One was compiled from the traditional ideal well path theory, while the other was com-
piled in the Landmark software Compass. Drag calculations in two dimensions are
presented for a range of friction coefficients, and the effective force is calculated and
compared to the critical buckling limit.
The traditional way of hand-calculating torque and drag is by dividing the well path into
ideal sections. The well path is presented in figure 4.5.
Figure 4.5: The ideal well path divided in sections. Each well section is indexed from
1 to 6. The subscripts between each section represents the order of drag-calculation.
Subscript 1 is before subscript 2 and so on.
The curved sections was calculated as if the build angle was constant throughout the
entire curve. The curved well length was then calculated by knowing the radius of cur-
vature. Figure 4.6 illustrates this technique.
4.3. BUCKLING DURING INSTALLATION 39
Figure 4.6: Illustration of an arbitrary build section to show the technique of calculating
section length in an ideal well path curve.
Consider the two vertical lines, H1 and H2 . These lines can be expressed by trigonome-
try:
H1 = Rsinα1 (4.3)
H2 = Rsinα2 (4.4)
The height of the section will now be the difference between these lines:
³ ´
h = H1 − H2 = R sinα1 − sinα2 (4.5)
h
R= (4.6)
sinα1 − sinα2
The measured length of the section is now calculated using the arc length equation:
∆L = R × α (4.7)
For each build section, the section length is calculated by equation (4.7). The technique
is illustrated for the first and second build section in Appendix C. The height and angle
of the sail section is known, and the section length is easily found by trigonometry. The
vertical and horizontal section lengths are given.
During drag-calculations, one begin with the lowermost section and calculate upwards.
It is convenient to use the ideal well path method for hand-calculations due to the re-
duced amount of iterations.
40 CHAPTER 4. RESULTS
Figure 4.7a presents the ideal well path, and figure 4.7b presents the real well path.
300 300
400 400
500 500
600 600
700 700
0 200 400 600 800 1000 1200 1400 1600 1800 2000 0 200 400 600 800 1000 1200 1400 1600 1800 2000
Horizontal departure [m] Horizontal departure [m]
Figure 4.7: Presentation of ideal well path and real well path.
Figure 4.8a shows the ideal and real well path in the same graph, while figure 4.8b shows
the deviations in the curved section.
Comparing ideal with real well path Comparing build up sections in ideal and real well path
0
450
100
460
Total vertical depth [m]
200
470
300
480
400
500 490
600
500
700
0 200 400 600 800 1000 1200 1400 1600 1800 2000 -2 0 2 4 6 8 10 12
Horizontal departure [m] Horizontal departure [m]
(a) Ideal and real well path plotted in the same (b) Ideal and real well path plotted in the same
graph graph. Zoomed in on the curved section.
Figure 4.8: Comparison of the ideal well path and real well path.
4.3. BUCKLING DURING INSTALLATION 41
The program is verified by comparing the static hook load to the dynamic hook load
with coefficient of friction equal to zero. The vertical projected height principle is used
to calculate the static hook load (Aadnoy, 2010).
Axial Force when lowering String Axial Force when lowering String
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
2000 2000
(a) µ=0 for the ideal well path (b) µ=0 for the real well path
Figure 4.9: Dynamic hook load and static hook load overlaps at µ=0
Figure 4.10, 4.11, 4.12 and 4.13 present the results from the drag analysis.
Notice how the hook load changes when the friction factor increases, and also notice
the difference between the ideal and the real well path.
Axial Force when lowering String Axial Force when lowering String
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
2000 2000
(a) µ=0.2 for the ideal well path (b) µ=0.2 for the real well path
Axial Force when lowering String Axial Force when lowering String
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
2000 2000
(a) µ=0.3 for the ideal well path (b) µ=0.3 for the real well path
Axial Force when lowering String Axial Force when lowering String
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
2000 2000
(a) µ=0.4 for the ideal well path (b) µ=0.4 for the real well path
Axial Force when lowering String Axial Force when lowering String
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
2000 2000
(a) µ=0.5 for the ideal well path (b) µ=0.5 for the real well path
The effective axial force is calculated. Well fluids will be present when running in hole,
and this will create an opposing force when lowering the string. This effect is included
when the true axial forces is calculated. Further, the effective axial force is calculated
and compared to the critical buckling limit.
Figure 4.14, 4.15 and 4.16 present the results after implementation of the critical buck-
ling limits. Sketches are included for a clear visualization of each buckling regime.
500
1000
MD [m]
1500
500
1000
MD [m]
1500
500
1000
MD [m]
1500
F [lb] ×10 4
When pressure and temperature is not changed in the well, i.e. during lowering of the
completion, permanent corkscrewing occurs when the bending stress at the outer wall
exceeds the yields strength. The derivation is conducted from the outer stress in equa-
tion (3.49). P i =P o gives:
P i (1 − R 2 )
σo = | + σa ± σb |
R2 − 1
= | − P i + σa ± σb |
−P i (A i − A o ) (4.8)
= | − Pi + ± σb |
As
= | ± σb |
= σb ≤ σ yi el d
Figure 4.17 shows a plot of triaxial stress versus the effective axial force. The triaxial
stress on the outer wall has been reduced to the bending stress. The relation between
bending stress and effective axial stress is linear.
If the bending stress exceeds the available grades, permanent corkscrewing will occur.
The effective axial force required to permanently deform the pipe is calculated and
marked in the plot. The minimum effective axial force required is ≈130000 psi, and
this force is not encountered in any of the drag scenarios. The lowest grade available is
therefore applicable for buckling scenarios during installation.
4.3. BUCKLING DURING INSTALLATION 45
×10
4 Permanent Corkscrewing
14
Triaxial stress (reduced to bending stress) [psi]
Sigmao Initial
12 Max Effective Axial Force
Grades
10
0
-1.5 -1 -0.5 0 0.5 1 1.5
Figure 4.17: The relation between effective axial stress and bending stress is linear.
When the bending stress exceeds the yield strength, i.e. grade, of the tubing, perma-
nent corkscrewing will occur.
46 CHAPTER 4. RESULTS
It is assumed that the completion fluid has the density equivalent to sea water, reservoir
pressure is 70bar and reservoir temperature is 17°C.
The neutral stability point is found by calculating the true and effective axial force for
an anchored tubing, thus including the packer force (B.S. Aadnoy, personal communi-
cation, May 27, 2018). The vertical projected height principle is used to calculate tubing
weight (Aadnoy, 2010). Inserting equation (3.30) and (3.31) in equation (3.40) results in
the following expression:
F f∗ Fe f f
n= = (4.9)
w w
The calculation procedure for axial forces at surface is presented:
F a∗ = TVD × w l + P o (A p − A o ) − P i (A p − A i ) + F p = F a + F p
F e∗f f = F a∗ − p i A i + p o A o = F e f f + F p
Further, the calculation procedure for axial forces at the packer is presented:
F a ∗ = P o (A p − A o ) − P i (A p − A i ) + F p = F a + F p
F e f f ∗ = F a ∗ −p i A i + p o A o = F e f f + F p
Capital P is the pressure at packer level, and lowercase p is the pressure at surface. Note
how the pressures are used in each equation. A flow chart for further understanding of
F e∗f f is presented.
4.4. BUCKLING DURING OPERATIONS 47
Packer permitting no
motion
Ballooning
dL
Ff
dL
Fp
1. Shut-in
2. Start bullheading
Figure 4.18 presents the temperature cycles for the load cases shut-in and start bull-
heading. It is assumed that this well is installed while circulating fluid with a tempera-
ture of 10 °C. The temperature cycle for start bullheading equals the temperature cycle
for production.
Figure 4.18: Temperature cycles for shut-in after installation and start bullheading after
production.
The initial condition is the temperature after installation. After shut-in, the temperature
will approach the geothermal gradient. After production, the temperature will increase
at wellhead. bullheading is initiated by applying a high surface pressure.
4.4. BUCKLING DURING OPERATIONS 49
4.4.2.1 Shut-in
The tubing will shorten due to ballooning and the negative thermal effect, and this
causes the packer to act with a tensile force. The packer force is now located within
the tensile region on the plot, and it means that the tensile force causes the tubing to
unbuckle. Figure 4.19 shows that the neutral stability point is located below the string.
Axial Force vs. Length Change True vs. Effective Axial Force
10 0 0
800
2
Depth [m]
Depth [ft]
dL [in]
1000 300
0
1200
-2 400
1400
-4
1600
500
-6
1800
-8 600
2000
-10
-6000 -4000 -2000 0 2000 4000 6000 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5
Axial Force [lb] Force [lb] ×10 4
(a) F f∗ > F f , thus tensile packer force (b) Neutral stability point due to shut-in
Figure 4.19: Packer force and neutral stability point during shut-in
Both scenarios were simulated to illustrate the effect. Figure 4.20a present the start
bullheading procedure in a gas filled well. Figure 4.20b show that the neutral stability
point was located within the build section.
Figure 4.20c show that the packer force was larger for the oil filled well than the gas filled
well. Figure 4.20d show that the neutral stability point was now located in the vertical
section.
50 CHAPTER 4. RESULTS
Axial Force vs. Length Change True vs. Effective Axial Force
15 0 0
5
NSP 200
800
Depth [m]
Depth [ft]
dL [in]
1000 300
0
1200
400
1400
-5
1600
500
-10 1800
600
2000
-15
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -5000 0 5000 10000 15000 20000
(a) Simulation of packer force in gas filled well (b) Neutral stability point in the gas filled well
Axial Force vs. Length Change True vs. Effective Axial Force
15 0 0
5
NSP 200
800
Depth [m]
Depth [ft]
dL [in]
1000 300
0
1200
400
1400
-5
1600
500
-10 1800
600
2000
-15
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.5 0 0.5 1 1.5 2 2.5
Axial Force [lb] ×10 4 Force [lb] ×10 4
(c) Simulation of packer force in oil filled well (d) Neutral stability point in the oil filled well
Figure 4.20: Packer force and neutral stability point for start bullheading.
4.4. BUCKLING DURING OPERATIONS 51
Axial Force vs. Length Change True vs. Effective Axial Force
40 0 0
10 800
Depth [m]
Depth [ft]
dL [in]
1000 300
0
1200
400
-10
1400
-20 1600
500
1800
-30
600
2000
-40
-3 -2 -1 0 1 2 3 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5
Axial Force [lb] ×10 4 ×10 4
Force [lb]
(a) F f∗ > F f , thus tensile packer force (b) Neutral stability point when packer is set
Figure 4.21: Packer force and neutral stability point during set packer
52 CHAPTER 4. RESULTS
Buckling induces bending stress, and this will be locally superimposed to the true axial
stress. A bend consists of both a tensile part and a compressive part. If the axial stress
is tensile, the tensile stress from the outer bend will be added. If the axial stress is com-
pressive, the compressive force from the inner bend will be added.
Figure 4.22 includes the bending stress when each buckling regime is clarified.
True vs. Effective Axial Force True vs. Effective Axial Force
0 0 0 0
800 800
Depth [m]
Depth [m]
Depth [ft]
Depth [ft]
1200 1200
400 400
1400 1400
1600 1600
500 500
1800 1800
600 600
2000 2000
(a) Set packer (b) Start bullheading with gas filled well
True vs. Effective Axial Force True vs. Effective Axial Force
0 0 0 0
800 800
Depth [m]
Depth [m]
Depth [ft]
Depth [ft]
1200 1200
400 400
1400 1400
1600 1600
500 500
1800 1800
600 600
2000 2000
-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5
Force [lb] ×10 4 Force [lb] ×10 4
Both sinusoidal and helical buckling will occur during the set packer load case. Both
buckling regimes will cause extra bending stresses, that locally superimpose to the true
axial force.
4.4. BUCKLING DURING OPERATIONS 53
Figure 4.23 plots initial slack-off and pick-up versus triaxial stress.
8000
2
7000
Triaxial Stress [psi]
1.5
5000
4000
1
3000
2000
0.5
Sigmao Initial Sigmao Initial
1000 Sigmao Sigmao
Sigmai Sigmai
0 0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Initial Slack-Off and Pick-Up Force [lb] ×10 4 Initial Slack-Off and Pick-Up Force [lb] ×10 4
Figure 4.23: Triaxial stress in tubing during start bullheading and set packer
It is clear that the triaxial stress does not exceed the minimum grade available for this
tubing, that is 40ksi. The tubing does not permanently corkscrew when slack-off/pick-
up forces are between -10000 and 10000 lb. However, it should be noticed that for slack-
off forces higher than 7500 lb, the bending stress induces higher triaxial stress at the
outer wall than the inner wall for start bullheading.
It is not expected that permanent corkscrewing will be an issue during any operations.
54 CHAPTER 4. RESULTS
Chapter 5
Table 5.1: A range of typical friction factors depending on lubrication and interfaces.
(Mitchell et al., 1996)
By evaluating the friction factors, one can determine which completion fluid that should
be used. It is evident that a completion fluid with lubrication is preferable.
When lowering the string, drag will occur due to the dynamic friction factor. If the com-
pressive forces exceeds the buckling limit and causes the string to get stuck or even lock
up, a static friction factor applies. The static friction factor is higher than the dynamic
friction factor. A tubing should be designed based on the worst-case-scenario, and it
can therefore a better choice to use the static friction factor to avoid potential lock-up.
To predict a realistic scenario, a variation of friction factors should be used.
The string is normally not rotated during tripping. When rotating the string, the friction
is reduced, thus it is common with lower friction factors in drilling operations. Dynamic
friction factors experienced while drilling are between 0.1 and 0.3.
55
56 CHAPTER 5. DISCUSSION AND EVALUATION
Some completions are designed for zone isolation. It is critical that the horizontal part
of the completion gets installed at the correct depth. A buckled liner or tubing will be
shorter than the initial length, and lack of predictions can cause tubing assemblies to
be at wrong depth in relation to the reservoir. An intervention to get a proper overview
of the tubing location can be costly, and is something that should be avoided.
It was expected that the both well paths would give similar curves. It was rather found
that the real well path was very much affected by the coefficient of friction. When the
coefficient of friction increases from 0.4 to 0.5, the tubing configuration locks-up for the
real well path, while the hook load remains in tension for the ideal well path.
In the ideal well path, the dog leg angle is constant throughout the curved sections. This
is not the case in the real well path. Since the real well path is adjusted to certain tar-
gets in Compass, small dog leg variations will occur. In Compass, this is seen as rapid
changes in inclination, creating an incremental curve instead of a smooth curve. In real
wells, dog leg variations will occur due to irregularities in the wellbore. This is called
tortuosity, and causes friction to build up while the string is lowered. Qualitatively, one
can assume that tortuosity is included in the friction factor (B.S. Aadnoy, personal com-
munication, May 30, 2018). One can therefore reach higher values than expected. One
should therefore carefully consider the dog leg variation and friction factors when in-
stalling a completion in a long-reach horizontal well.
The critical limit for helical buckling is higher than the compressive force in the hori-
zontal section. Drag in the straight section was calculated under the assumption that
weight dominates. One will not expect any helical buckling here, as the buckling limit
is high and the tubing is stable. However, the growth of sinusoidal buckling can be ex-
pected, and this causes associated increase in contact forces (Payne et al., 1996).
5.1. BUCKLING DURING INSTALLATION 57
For the deviated sections, equation (3.10) is used. The limit increases significantly in
the build up from 0 to 85 degrees. The limit in the tangential section is constant before
it increases slightly from 85 to 90 degrees.
According to Mitchell et al. (1999), there is a transition period between 1.4 and 2.8 times
the Paslay Dawson equation where helical buckling can occur. Figure 5.1 illustrates how
the limits are defined.
500
1000
MD [m]
1500
2000
Sinusoidal = Fcr
Minimum Helical = 1.4Fcr
Maximum Helical = 2.8Fcr
2500
-4 -3 -2 -1 0 1 2 3 4
F [lb] ×10 4
Figure 5.1: Sinusoidal buckling limit is defined by the Paslay Dawson equation. The
transition between sinusoidal and helical buckling occurs between 1.4 and 2.8 times
the Paslay Dawson limit.
For the highest friction considered, µ=0.45, one can not see any indications of helical
buckling in the straight, deviated sections. Using the Paslay-Dawson equation in the
curved sections, indicates a transitional region for helical buckling.
According to He et al. (1995), a curved section will provide additional support to a tub-
ing in compression. If one includes the contact forces initiated by curvature effects like
build-up/azimuth rate, the buckling limit can actually be higher than expected. If the
limit for sinusoidal buckling is higher than expected in the first build up section, the
pipe will remain straight at the low side of the bend.
58 CHAPTER 5. DISCUSSION AND EVALUATION
The tangential section will not experience any buckling as it is below the sinusoidal
buckling limit. Moving further down to the horizontal section in the well, the model
proposed by He et al. (1995) does not apply. Hence buckling can be present in this
straight section, even if it is suppressed in the build sections.
Figure 5.2: Buckling during installation when µ=0.45. The effects of curvature is in-
cluded in the sketch.
Consider a case where the effective axial force exceeds the helical buckling limit in the
horizontal section.
According to table 5.1, a friction factor of 0.6 matches the combination "clean surface,
metal on metal" or "unlubricated, metal on unmetal". It can therefore theoretically be
a scenario when running pipe in an open hole with water based mud.
If it is assumed that any curvature increases the buckling limit over the operational
force, the following scenario is expected: The tubing will helically buckle in the vertical
section. It will not be buckling in the curved section. Helical or sinusoidal buckling will
be present in the tangential section, while there is no buckling in the next curved sec-
tion. Finally, there will be either helical or sinusoidal buckling in the horizontal section.
Figure 5.3 presents the hypothetical scenario of installation of the completion with a
friction factor of 0.6.
5.1. BUCKLING DURING INSTALLATION 59
500
1000
MD [m]
1500
F [lb] ×10 4
(b) Sketch of buckling scenario when µ=0.6
Previous knowledge have shown that the tubing locks-up for a friction factor of 0.45.
Adding more weight on surface, will only cause the vertical section to buckle more
severely, without affecting the tubing in the horizontal section. Figure 5.3 is therefore
an unrealistic scenario for this well path.
Analyzing a real well path had the advantage that it included the tortuosity effect in the
build up sections. The dog leg variation caused friction to build up in the string, thus
pushing the neutral stability point higher up than expected.
A shortcoming of the model is that it was simulated with a constant friction factor. In
real well operations, the friction factor will vary due to fluid parameters like viscosity
and density. If the well is not clean, remains could cause friction to build up.
In this analysis it was assumed that the tubing string was lowered within a casing that
ensures a constant radial clearance. If the string was entering a section with a larger
radial clearance, it would decrease the buckling limit in the deviated section, hence in-
creasing the possibility of buckling. Figure 5.4 shows that the radial clearance affects
whether there is a possibility for helical buckling in the horizontal section or not. A
lower radial clearance will stabilize the tubing.
60 CHAPTER 5. DISCUSSION AND EVALUATION
Axial Force when lowering string Axial Force when lowering string
0 0
500 500
1000 1000
MD [m]
MD [m]
1500 1500
(a) µ=0.45 for tubing within 6 5/8" casing (b) µ=0.45 for tubing within 9 5/8" casing
Figure 5.4: Comparison of critical buckling limits for a tubing inside a 6 5/8" casing and
9 5/8" casing
Radial clearance is shown to be of importance. The largest casing increases the possi-
bility of buckling in the deviated sections. There are two situations that will cause the
radial clearance to increase: a smaller tubing outer diameter, or a larger casing inner
diameter.
If the tubing were to enter the open hole, both friction and radial clearance would in-
crease. This is not accounted for in the model.
When analyzing drag forces while lowering the string, it was convenient to include the
buckling limits in every section.
Contact force due to helical buckling was not included in the drag model. Weltzin et al.
(2009) did a full-scale experiment where no friction where experienced in the vertical
section, even though helical buckling was present. The theoretical model can therefore
be correct when not including friction in the vertical part. Weltzin et al. (2009) con-
cluded that friction from sinusoidal buckling could be of significance in deviated wells.
The contact force from sinusoidal buckling is calculated in the same manner as the nor-
mal force from the tubing weight (Mitchell et al., 1999). The contact force induced by
helical buckling would probably replace the normal force from the tubing weight.
An argument for simulating scenarios for high friction factors is to create a worst case
scenario, as is common practice in tubing design. The program should be able predict
the earliest possible time of a potential lock-up situation.
5.1. BUCKLING DURING INSTALLATION 61
A future improvement is to implement vectors with exact measurements for each part
of the well. A vector consisting of friction factors adjusted for the realistic scenario
would give a more reliable result. A vector could also include different radial clear-
ances. Variations could be due to wear, tubing connections or crossovers. An open hole
will also have varying diameter due to washouts and other irregularities.
The contact force due to helical buckling in the deviated sections should be included if
buckling is present. The induced contact force can make it more difficult to reach the
desired depth. This is particularly important for perforation and zone-isolation.
2. Well path configuration. Lubinski et al. (1962) uses a vertical well in the model.
The correct measurements for TVD and MD had to be implemented in the script
for accurate hydrostatic pressure calculation.
3. Packer configuration. Lubinski et al. (1962) uses a packer configuration that gives
A p 6= A i . A polished bore receptacle (PBR) is assumed to be installed, hence the
packer bore, A p , equals the inner area, A i .
4. The well design pressure (WDP) had to be considered. The wellhead on a shal-
low, horizontal well will have a lower tolerance for pressure than a deep, verti-
cal well. The load cases therefore needed to be adjusted to WDP. For this well,
WDP=1892psi.
5. The load cases involves fluids in static conditions. Equation (3.24) is therefore
replaced by equation (3.25) in the model to evaluate the ballooning effect in static
conditions.
Buckling was not associated with neither the shut-in or start bullheading load case.
These are considered to be of secondary importance, and the discussion of these is
therefore added to Appendix A.
5.2. BUCKLING DURING OPERATIONS 63
For the set packer situation, the neutral stability point is located in the vertical section.
Figure 5.5 shows results of the simulation. The effective force at the bottom of the string
is also higher than the effective force at the bottom of the vertical section.
in
Fa Load Case
200
Feff Load Case
Fa with Bending 100
400
NSP
Sinusoidal buckling
600 Minimum helical buckling 200
800
Depth [m]
Depth [ft]
1000 300
1200
400
1400
1600
500
1800
600
2000
-3 -2 -1 0 1 2 3 4
Figure 5.5: The effective axial force exceeded the critical buckling limits in the vertical
section. No buckling is shown below the vertical section
64 CHAPTER 5. DISCUSSION AND EVALUATION
The Lubinski et al. (1962) method is developed for a vertical well. It was assumed that a
fictitious force at the bottom of the tubing will induce helical buckling from the bottom
of the string up to the neutral stability point. The length change is in other words cal-
culated from the effective force at packer depth.
The buckling length change calculation is more complex when several buckling regimes
are evaluated in a deviated well.
The critical buckling limit for the vertical section is compared with the effective axial
force. The limits suggests that the vertical region below the neutral stability consists of
three scenarios: first no buckling, then sinusoidal buckling and finally helical buckling.
Figure 5.5 shows that buckling is suppressed in the build-up section. Due to the dis-
cussion in 5.1.3, the limits can be expected to be even higher. It is therefore confidently
assumed that no buckling occurs within this region.
Let the vertical section in figure 5.5 be considered as an isolated system. The effective
axial force acting on the bottom of the tubing in the vertical section will be the force
that induces the helical buckling. It is therefore this force that should be used to give
more accurate predictions for helical buckling length changes.
The vertical region consists of both sinusoidal buckling and helical buckling. Equation
(3.21) for helical buckling will overestimate the length change, while equation (3.41)
for sinusoidal buckling will underestimate the length change. The real length change
is found between these solutions. A weighted average should be calculated to find an
approximate length change over the vertical section.
The model is set by default to calculate the buckling length change associated with he-
lical buckling. One single iteration is required to calculate the length change associated
by sinusoidal buckling instead. This is assumed reasonable for vertical well sections
due to the uncertainty on the critical buckling limits in the vertical well. The transition
zone between helical and sinusoidal buckling can be relatively small, and the helical
buckling length change will give the worst-case scenario. Buckling length change is de-
viated sections should be modified, as the transition between sinusoidal buckling and
helical buckling would be of significance.
However, the effective axial force used in the calculation is the force at packer depth.
This is not correct, and a weakness in the model.
5.2. BUCKLING DURING OPERATIONS 65
2. Calculate the effective axial force as if the tubing was free to move
4. Calculate buckling length change with the effective axial force at the bottom of
the string
5. Find packer force and plot the effective axial force with the critical buckling limits
6. If buckling only occurs in the vertical section, find effective axial force at the bot-
tom of the helical buckling region and sinusoidal buckling region
7. Return to the buckling length change calculation, and calculate the a new weighted
average of the buckling length change using the force for helical and sinusoidal
buckling respectively.
8. Repeat procedure. Calculate packer force and plot the effective axial force with
the critical buckling limits.
The length change associated with buckling is often the smallest length change. The
consequences of over- or underestimating this length changes is not severe compared
to the other effects. However, this contributes to the uncertainty that revolves around
axial forces. A more accurate prediction of buckling length change will improve axial
analysis.
66 CHAPTER 5. DISCUSSION AND EVALUATION
Figure 5.6 shows where the neutral stability point was located during each load case.
True vs. Effective Axial Force True vs. Effective Axial Force
0 0 0 0
300 300
100 100
Depth [m]
Depth [m]
Depth [ft]
Depth [ft]
400 400
600 600
200 200
700 700
800 800
250 250
(a) Set packer (b) Start bullheading with gas filled well
True vs. Effective Axial Force True vs. Effective Axial Force
0 0 0 0
300 300
100 100
Depth [m]
Depth [m]
Depth [ft]
Depth [ft]
400 400
600 600
200 200
700 700
800 800
250 250
Figure 5.6: Location of neutral stability point in a shallow subsea well with wet tree
The set packer case causes the neutral point be located above the string. It means that
compression will be expected at the wellhead. At the wellhead, the main concern is
wellhead growth. This takes place when the tubing is under compression at the tubing
hanger. As the wellhead is an important element in the primary barrier envelope, this
should be carefully considered. To avoid this problem, a packer that require a lower set
pressure should be considered.
5.2. BUCKLING DURING OPERATIONS 67
Tubing and casing properties are defined as input. This makes it convenient to investi-
gate whether the impact of radial clearance on buckling. The model is further general-
ized to work for any load case. The load cases can easily be implemented in the model
by defining densities, pressures and temperatures.
The model developed by Lubinski et al. (1962) does not account for couplings or con-
nections. It is possible that high bending stresses causes the threads to deform, thus
making connections the weak points in the completion. Failures related to connections
might be both structural and leak failures (Bellarby, 2009, 544). A future improvement
would be to include the effects of bending stress on connections.
The buckling models presented in this thesis does not include friction. Vibrations will
most likely occur during production or injection, and cause the frictional forces to be
negligible with regards to buckling severity. However, will vibration effects not neces-
sarily mitigate friction in highly deviated wells. It is previously presented by Mitchell
et al. (1996) that friction can affect the axial stress during production or injection. A
future improvement would be to evaluate buckling by considering friction.
The entire model takes on average between 2 seconds to run. It is therefore a time
efficient method to evaluate length changes, packer forces, neutral stability point, per-
manent corkscrewing and critical buckling.
The theoretical expression for pitch is given by Lubinski et al. (1962) in equation (3.7).
The pitch is defined as the length between to spirals. The pitch will decrease with in-
creasing compression, and will the tubing buckle more severely at the bottom part of
the string. However, assumptions behind this equation limits the application. The pitch
in equation (3.7) is defined as the distance between spirals, just above the lower end of
the tubing. The pitch is not included in this model due to its limited practical applica-
tion.
A future improvement would be to estimate a variation of the pitch over an entire buck-
led interval. This would make it possible to determine the exact location of the string in
the well. This technology could also be applied in drilling engineering when hole clean-
ing is a problem. While circulating the well, it will be beneficial that the pipe is located
at high side. This will create a flow path below the pipe, and the accumulated mud or
dirt can be flushed out. (S. Sangesland, personal communication, March 13, 2018)
68 CHAPTER 5. DISCUSSION AND EVALUATION
Lubinski et al. (1962) developed the model from a vertical well, while Mitchell et al.
(1999) added a sail section directly to the vertical section. Neither of those assumptions
are valid when developing the model for a horizontal well. The main problems are the
temperature and pressure calculations. A future improvement is to include exact mea-
surements of pressure and temperatures in the entire well. In a deviated well, the hy-
drostatic pressure will be the same throughout the entire horizontal section, causing the
mean pressure to be underestimated by the traditional arithmetic average. A weighted
average of the pressure and temperature would give a better estimation.
The theory behind critical buckling limits suggests that the effective axial force needs
to exceed a critical limit before sinusoidal buckling occurs. I.e. the tubing is assumed
perfectly straight before the sinusoidal is developed. The critical buckling limits deter-
mined for a vertical well is debated (Cunha et al., 2003), and different coefficients to
determine sinusoidal and helical buckling is suggested. A long tubing run in a well can
be subjected to multiple defections. In a vertical well, the practical importance of the
different coefficients is therefore limited.
An experimental and analytic study of sinusoidal buckling in vertical wells was con-
ducted by Salies et al. (1994), and it was found that equation (3.8) by Lubinski et al.
E I 1/3
(1950) is a good approximation for a a tubing length equivalent to 7.94 × ( w ) . For
b
L
strings with lengths greater than this, > 7.94, the critical limit is less. A quick cal-
( wE I )1/3
b
culation shows that the length of the vertical section considered in this thesis is equiv-
E I 1/3
alent to 33 × ( w ) . Hence, the analytic limit for sinusoidal buckling can be less than
b
expected according to Salies et al. (1994). Further, the experiments conducted in Salies
et al. (1994) confirmed that the critical buckling limit was less. When the length was
further increased, the critical limit decreased towards a constant value.
The author assumes that geometrical imperfections will determine when the real criti-
cal limit is reached. The longer the pipe is, the more deflections is present. It is therefore
possible that sinusoidal buckling occurs between F e f f = 0 and F e f f = 1.94(E I w b2 )1/3 in
a vertical well. Full-scale experiments should be performed to give a better foundation
for the analytic models.
WellCat uses equation (3.10) to determine whether the tubing buckles or not in straight
deviated sections (Horgen, 2010). This critical buckling limit determines when stable
sinusoidal buckling occurs. By using the Paslay Dawson equation to predict the onset
of helical buckling, an implicit safety factor will be included. This is therefore a conser-
vative method to predict buckling. It can be useful for critical well applications (Aasen
et al., 2002), but it contributes to the uncertainty of axial force analysis.
The model presented in this thesis includes the critical buckling limits, thus opens the
possibility to investigate this further. The simulation of a horizontal well gives a more
complex result than a vertical well. The study of buckling during operations show that
WellCat is conservative in its buckling analysis for both vertical and straight deviated
wells.
70 CHAPTER 5. DISCUSSION AND EVALUATION
Chapter 6
Conclusion
• A critical buckling ratio is suggested. This method can be used to see if one tubing
is better than another with regards to buckling.
• A buckling model has been developed for installation of a completion. The model
accounts for dog leg variations in all well sections, and predicts lock-up in the
string. The drag forces were underestimated when dog leg variations was not
included.
• WellCat is found to perform conservative buckling analysis for both vertical and
straight, deviated wells.
• Buckling analysis is improved by supplementing the existing models with the lat-
est of theory. These models were modified to be applicable for all well designs.
• The model developed for well operations is an integrated model that calculates
the axial forces, the packer force and permanent corkscrewing. An initial installa-
tion strategy can also be chosen by implementing a variation of slack-off or pick-
up forces.
• Transforming the subsea well from a dry tree to a wet tree causes significant dif-
ference because of the reduction in vertical extent. Wellhead growth is expected
to be a problem for the subsea well during the set packer procedure.
• Buckling is not equivalent to failure. Permanent corkscrewing did not occur for
installation nor critical operations.
71
72 CHAPTER 6. CONCLUSION
Chapter 7
Further Work
• Include the critical buckling limit by He et al. (1995) for the curved sections.
• Implement vectors to account for a variation parameters like radial clearance and
friction factor in the installation phase.
• Develop an equation for the pitch over the entire buckled interval.
73
74 CHAPTER 7. FURTHER WORK
7.3 General
Further work is to investigate the tubing condition from installation to further opera-
tions. This model will be able to account for tubing wear from the installation process.
Bibliography
Aadnoy, B. S., Cooper, I., Miska, S. Z., Mitchell, R. F., and Payne, M. L. (2009). Advanced
drilling and well technology. SPE.
Aadnoy, B. S., Fazaelizadeh, M., Hareland, G., et al. (2010). A 3d analytical model for
wellbore friction. Journal of Canadian Petroleum Technology, 49(10):25–36.
Aasen, J. A., Aadnøy, B. S., et al. (2002). Buckling models revisited. In IADC/SPE Asia
Pacific Drilling Technology. Society of Petroleum Engineers.
Cunha, J. et al. (2003). Buckling of tubulars inside wellbores: a review on recent theoret-
ical and experimental works. In SPE Production and Operations Symposium. Society
of Petroleum Engineers.
Dawson, R. et al. (1984). Drill pipe buckling in inclined holes. Journal of Petroleum
Technology, 36(10):1–734.
He, X., Kyllingstad, A., et al. (1995). Helical buckling and lock-up conditions for coiled
tubing in curved wells. SPE Drilling & Completion, 10(01):10–15.
Hearn, E. J. (1997). Mechanics of Materials 2: The mechanics of elastic and plastic defor-
mation of solids and structural materials. Butterworth-Heinemann.
Horgen, T. (2010). Analysis of packer forces in vertical wells and frictionless sail sections.
Master’s thesis, University of Stavanger, Norway.
Lubinski, A., Althouse, W., et al. (1962). Helical buckling of tubing sealed in packers.
Journal of Petroleum Technology, 14(06):655–670.
Lubinski, A. et al. (1950). A study of the buckling of rotary drilling strings. In Drilling
and Production Practice. American Petroleum Institute.
75
76 BIBLIOGRAPHY
Mirhaj, S., Kaarstad, E., Aadnoy, B., et al. (2016). Torque and drag modeling; soft-string
versus stiff-string models. In SPE/IADC Middle East Drilling Technology Conference
and Exhibition. Society of Petroleum Engineers.
Mitchell, R. et al. (1988). New concepts for helical buckling. SPE Drilling Engineering,
3(03):303–310.
Mitchell, R. et al. (1996). Comprehensive analysis of buckling with friction. SPE Drilling
& Completion, 11(03):178–184.
Mitchell, R. et al. (1997). Effects of well deviation on helical buckling. SPE Drilling &
Completion, 12(01):63–70.
Mitchell, R. et al. (1999). Buckling analysis in deviated wells: A practical method. SPE
Drilling & Completion, 14(01):11–20.
Payne, M., Abbassian, F., et al. (1996). Advanced torque and drag considerations in
extended-reach wells. In SPE/IADC Drilling Conference. Society of Petroleum Engi-
neers.
Remmen, T., Bauge, E., and Wilhelmsen, K. (2017). Analysis of the industry standard
and practice for casing and tubing design. In Project Report. NTNU.
Salies, J., Cunha, J., Azar, J., Soren Jr, J., et al. (1994). Experimental and analytical study
of sinusoidal buckling in vertical wells. In SPE Eastern Regional Meeting. Society of
Petroleum Engineers.
Weltzin, T., Aas, B., Andreassen, E., Lindland, M., et al. (2009). Measuring drillpipe buck-
ling using continuous gyro challenges existing theories. SPE drilling & completion,
24(04):464–472.
Wu, J., Juvkam-Wold, H. C., et al. (1995). Coiled tubing buckling implication in drilling
and completing horizontal wells. SPE Drilling & Completion, 10(01):16–21.
Appendix A
Additional Information
Fa
sinθ = (A.1)
F
P
sinθ = p (A.2)
P 2 + 4π2 r 2
Combining equation (A.1) and (A.2) gives the expression for F a :
P
Fa = F p (A.3)
p 2 + 4π2 r 2
77
78 APPENDIX A. ADDITIONAL INFORMATION
Conservation of energy is used to find the force required to buckle a string. The princi-
ple goes as following:
Ei = E f (A.4)
where E i and E f is the initial and final energy in the system, respectively.
Energy can not just occur nor disappear, hence we say that energy is conserved. The
initial energy must equal the final energy. We consider the tubing string to be a static
system, hence the kinetic energy is zero. Before the tubing is subjected to deforma-
tion, potential energy will build up in the system. The potential energy consist of strain
energy of compression, strain energy of bending and the potential energy of a force sub-
jected to the material.
1 V 2 1 A s L F a2 1 L F a2
U st r ai n = σ = = (A.5)
2E 2 E A 2s 2 E A s
where V is volume, A s is the cross-sectional area and σ is the axial stress.
By combining equation (A.3) and (A.5), the strain energy of compression, Uc is ob-
tained:
F a2 L F 2P 2L
Uc = = (A.6)
2A s E 2A s E (P 2 + 4π2 r 2 )
When a string is bent, it will store energy from the bending. This is called the bending
strain energy, Ub , and it is given by:
LE IC 2
Ub = (A.7)
2
where C is the curvature of the helix.
When the curvature of the helix is given by
4π2 r
C= (A.8)
p 2 + 4π2 r 2
the strain energy of bending becomes:
8π4 r 2 E I L
Ub = (A.9)
P 2 + 4π2 r 2
From geometric similarities, the potential energy from the subjected force is obtained:
FPL F 2P 2L
U f = F Lh = p − 2 2 2
(A.10)
P 2 + 4π2 r 2 A s E (P + 4π r )
where U f is the potential energy from a subjected force, and L h is the length of the helix
measured along its axis.
The total potential energy, U , of the system is equal to the sum of strain energies Uc
A.1. DERIVATION OF PITCH 79
and Ub and the potential energy from the subjected force, U f . Equation (A.6), (A.9) and
(A.10) is combined to give the final expression for U :
F 2P 2L 8π4 r 2 E I L FPL
U =− + +p (A.11)
2A s E (P 2 + 4π2 r 2 ) (P 2 + 4π2 r 2 )2 p + 4π2 r 2
2
dU
=0
dP
(A.12)
P (P 2 + 4π2 r 2 ) 2
F − (P 2 + 4π2 r 2 )3/2 F + 8π2 E I P = 0
As E
The equation was the derivative of total potential energy with respect to the pitch with
constant r. Once a helical buckled pipe contacts the wall, r is constant, and P is there-
fore the only variable (Kwon et al., 1988).
This expression results in two roots. The smallest root corresponds to the smallest po-
tential energy, hence we solve the equation for the smallest root. The following expres-
sion is obtained:
p s
A s E P 2 + 4π2 r 2 h 32π2 I P 2 i
F= 1− 1− 2 2 2 2
(A.13)
2P A s (P + 4π r )
Two assumptions were made to simplify this expression. First, P 2 À 4π2 r 2 . This as-
sumption is found reasonable, as it is backed up by other authors (Mitchell et al., 1988).
The other assumption is that the second term in the square rootp will be small. For prac-
tical reasons, lets call this term a. Using Taylor expansion on 1 − a, the expression can
be simplified to 1 − a2 . Equation (A.13) can then be reduced to:
As E h 16π2 I P 2 i 8π2 E I
F= 1 − (1 − ) = (A.14)
2 As P 4 P2
By rearranging equation (A.14), equation (3.7) in section 3.2.2 is obtained.
80 APPENDIX A. ADDITIONAL INFORMATION
A.3 Shut-in
During installation, fluids will be circulated in the well. The circulating fluid is affected
by surface temperature, and it will be warmer than the sea water temperature. If shut-
in is required, the fluids in the well will become static, and the heat will transfer to the
surroundings. In practice, this means that the temperature will approach the geother-
mal gradient. The highest part of the tubing will experience the highest temperature
decrease. This will make the tubing contract, and consequently will more tension build
up in the tubing string.
It is assumed that gas fills up the well, and the inner pressure will increase. Figure A.4
shows the true and effective axial force when the well is shut in. The length changes
associated with this load will only be negative, thus the neutral stability point conse-
quently moves downwards. In this case, the point is located below the string. As the
tubing is entirely in tension through the build sections, it is expected that it is located
at the high side of the curvatures.
Fa Initial Conditions
200
Feff Initial Conditions
Fa Load Case 100
400
Feff Load Case
NSP
600 Sinusoidal buckling 200
Helical buckling
800 Depth [m]
Depth [ft]
1000 300
1200
400
1400
1600
500
1800
600
2000
Figure A.4: The negative length changes causes the entire tubing to be in tension during
the shut-in load case.
82 APPENDIX A. ADDITIONAL INFORMATION
The highest temperature difference in the well after production will be present at the
tubing head. This thermal effect will cause the tubing to expand. This positive length
change will contribute to a higher compressive load in the tubing.
Effects encountered due to increased inner pressure are ballooning and buckling. Bal-
looning will have a shortening effect on the tubing, thus reduce the compressive loads.
Increasing inner pressure also causes structural instability due to geometrical imper-
fections, and consequently decreases the effective axial force. Buckling will also occur,
and there will be a shortening associated with this effect.
Figure A.5 shows the true and effective axial force when start bullheading is initiated in
a gas-filled well. The negative length changes causes the tubing above the horizontal
section to be in tension. The inner pressure applied at surface will cause the effective
axial force to be shifted to a more compressive state than it was initially, but it is still not
enough to expect buckling to occur in any of the sections. The neutral stability point
will be located in the build section, and the neutral point is located at the bottom of the
build section. This indicates that most part of the tubing is located at the high-side of
the curve.
Fa Initial Conditions
200
Feff Initial Conditions
Fa Load Case 100
400
Feff Load Case
NSP
600 Sinusoidal buckling 200
Helical buckling
800
Depth [m]
Depth [ft]
1000 300
1200
400
1400
1600
500
1800
600
2000
Figure A.5: Critical buckling for start bullheading in a gas filled well.
A.4. START BULLHEADING 83
Figure A.6 shows the true and effective axial force when start bullheading is initiated in
a oil filled well. The negative length changes are dominating, indicating that the tubing
will be completely in tension. However, the high inner pressure will cause the effective
axial force to shift towards a more negative state than it was initially. The neutral stabil-
ity point will consequently move upwards, and be located in the vertical section. The
effective axial force is not high enough to cause buckling in any of the lower sections.
Fa Initial Conditions
200
Feff Initial Conditions
Fa Load Case 100
400
Feff Load Case
NSP
600 Sinusoidal buckling 200
Helical buckling
800
Depth [m]
Depth [ft]
1000 300
1200
400
1400
1600
500
1800
600
2000
Figure A.6: Critical buckling for start bullheading in a oil filled well.
The start bullheading load case covers the production load case with regards to buck-
ling. During production, the temperature increase will equal the start bullheading case.
However, the surface pressure is forced to be very low to be able to produce from the low
reservoir pressure. The inner pressure will thus be lower than in the start bullheading
scenario.
84 APPENDIX A. ADDITIONAL INFORMATION
Appendix B
This appendix includes the MatLab scrips used for compiling results from critical buck-
ling limit. The limit is compiled for sinusoidal and helical buckling in a vertical section,
straight deviated section of 45 ° and a horizontal section.
85
86 APPENDIX B. CRITICAL BUCKLING LIMIT
Contents
Critical Buckling Limit
Critical Buckling Force in a vertical section
Critical Buckling Force in an inclined section @ 45 deg
Critical Buckling Force in a horizontal section
Critical Buckling Force in an inclined section @ 0-90 deg
Plotting
clear
clc
% Tubing inputs
% Fluid properties
rhosw = 1.02; %sg sea water density
rhosteel = 7.89; %sg steel density
buoyedw = (1-rhosw/rhosteel).*weight/12; %lb/in buoyed weight
% D/t-ratio is a vector.
for i = 1:length(dtratio)
if dtratio(i) < 13.38
% Thick-walled tubing: Yield collapse
disp(['Yield Collapse for tubing ', num2str(OD(i)), ...
'in and #', (num2str(weight(i)))]);
elseif dtratio(i) < 22.47
% Plastic Collapse
disp(['Plastic Collapse for tubing ', num2str(OD(i)), ...
'in and #', (num2str(weight(i)))]);
else
% Thin-walled tubing: Elastic collapse
87
% Casing Inputs
%csgID = [32 32 32 32 32 32 32 32];
csgID = [5.921 5.921 5.921 5.921 5.921 5.921 5.921 5.921 5.921]; %in
% Radial Clearance
rc = csgID/2-OD/2; %in
varincl = linspace(0,90,10);
Plotting
figure(1)
plot(rc,Fcsinincl,'ko',rc,Fchelinclmax,'k*','markersize',10)
grid minor
legend('\fontsize{16}Sinusoidal buckling','\fontsize{16}Helical buckling')
title('\fontsize{16}Effect of radial clearance on critical buckling limit')
xlabel('\fontsize{16}Radial clearance [in]');
ylabel('\fontsize{16}Critical buckling limit [lb]');
figure(2)
plot(buoyedw,Fcsinincl,'ko',buoyedw,Fchelinclmax,'k*','markersize',10)
grid minor
legend('\fontsize{16}Sinusoidal buckling','\fontsize{16}Helical buckling' ...
, 'location', 'northwest')
title('\fontsize{16}Effect of weight on critical buckling limit')
xlabel('\fontsize{16}Weight [lb/in]');
ylabel('\fontsize{16}Critical buckling limit [lb]');
figure(4)
plot(I,ratio1,'ko',I,ratio2,'k*','markersize',10);
xL = xlim;
line(xL, [1 1],'color','black','linestyle','--'); %x-axis
grid minor
legend('\fontsize{16}Vertical well','\fontsize{16}Inclined/Horizontal well','location','nor
thwest');
title('\fontsize{16}Effect of I on critical buckling ratio');
xlabel('\fontsize{16}Moment of inertia [in4]');
89
ylabel('\fontsize{16}Fcr*');
figure(5)
plot(rc,ratio1,'ko',rc,ratio2,'k*','markersize',10);
xL = xlim;
line(xL, [1 1],'color','black','linestyle','--'); %x-axis
grid minor
legend('\fontsize{16}Vertical well','\fontsize{16}Inclined/Horizontal well','location','nor
theast');
title('\fontsize{16}Effect of radial clearance on critical buckling ratio');
xlabel('Radial clearance [in]');
ylabel('\fontsize{16}Fcr*');
figure(6)
plot(varincl,Fcsinvarincl,'ko',varincl,Fchelvarinclmin,'k*',varincl,Fchelvarinclmax,'kx','m
arkersize',10);
grid minor
legend('\fontsize{16}Sinusoidal buckling','\fontsize{16}Min helical buckling','\fontsize{16
}Max helical buckling','location','northwest');
title('\fontsize{16}Effect of inclination on critical buckling limit');
xlabel('Inclination [deg]');
ylabel('\fontsize{16}Critical buckling limit [lb]');
90 APPENDIX B. CRITICAL BUCKLING LIMIT
91
92 APPENDIX B. CRITICAL BUCKLING LIMIT
This appendix includes the MatLab scripts used for constructing well paths and per-
forming drag calculation with a variation of friction factors.
The first script constructs the ideal well path. The technique described in section 4.3.1
is used to calculate the build section length. Drag calculations are performed using the
traditional method: one section at the time, using one subscript for each section.
The second script imports data from Compass, and constructs the real well path. Drag
calculations are performed by creating 99 index points, thus including dog leg varia-
tions.
The third and last script uses the real well path, and implements the critical buckling
limits.
93
94 APPENDIX C. BUCKLING DURING INSTALLATION
Contents
INPUTS
Three-dimensional model by Bernt S. Aadnoy
Ideal Wellpath
clc
clear all
INPUTS
% Wellpath Calculations
TVD = [669.62 667 660 650 450 0]; % total vertical depth [m]
incl = [89.9 89.9 85 85 0 0]; % well inclination [deg]
incl = incl*pi/180; % well inclination [rad]
DL = zeros(1,length(L));
for i=1:(length(L)-1)
DL(i+1) = acos (sin(incl(i))*sin(incl(i+1))*cos(azi(i)-azi(i+1)) ...
+ cos(incl(i))*cos(incl(i+1)));
end
% If the static hook load equals the axial force with no friction,
% then the code is correct.
staticweight = zeros(1,length(TVD));
statichookload = zeros(1,length(TVD));
for i = 1:length(TVD)
staticweight(i) = TVD(i)*BF*wDP;
statichookload(i) = staticweight(1)-staticweight(i);
end
statichookload = statichookload*224.8089; % kN to lb
F = zeros(1,length(L));
hookload = zeros(1,length(my));
F(1) = WOB;
figure(j)
plot(F*224.8089,MD,'k',statichookload,MD,'k-.') % convert from [kN] to [lb]
xL = xlim;
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
set(gca,'YDir','reverse')
axis([-20000 20000 0 MD(1)])
grid minor
title('\fontsize{16}Axial Force when lowering String')
legend(['\fontsize{12}Dynamic Hook Load \mu =' num2str(my(j))] ...
96 APPENDIX C. BUCKLING DURING INSTALLATION
end
Ideal Wellpath
x1 = 0;
x2 = 0;
x3 = R(5)-sqrt(R(5)^2-(TVD(4)-TVD(5))^2);
x4 = x3 + (cos(pi/2-incl(3))*L(4));
x5 = x4 + sqrt(R(3)^2-(R(3)-(TVD(2)-TVD(3)))^2);
x6 = x5+cos(pi/2-incl(1))*L(2);
x = [x1 x2 x3 x4 x5 x6];
% Vertical section
verticalsection = linspace(TVD(6),TVD(5),10);
for i = 1:length(verticalsection)
xv(i) = 0;
yv(i) = verticalsection(i);
end
97
% Build section
build1 = linspace(TVD(5),TVD(4),10);
inclb = linspace(incl(5),incl(4),10);
for i = 1:length(build1)
yb(i) = TVD(5)-R(5)*sin(pi+inclb(i));
xb(i) = R(5)+R(5)*cos(pi+inclb(i));
end
% Tangent section
tangentsection = linspace(0,L(4),10);
for i = 1:length(tangentsection)
xt(i) = x(3) + tangentsection(i)*cos(-5*pi/180);
yt(i) = TVD(4) - tangentsection(i)*sin(-5*pi/180);
end
% Build section
build2 = linspace(TVD(3),TVD(2),10);
inclb2 = linspace(incl(3),incl(2),10);
for i = 1:length(build2)
yb2(i) = (TVD(2)-R(3))-R(3)*sin(pi+inclb2(i));
xb2(i) = x(5)+R(3)*cos(pi+inclb2(i));
end
% Horizontal section
horizontalsection = linspace(0,L(2),10);
inclh = acos((TVD(1)-TVD(2))/L(2));
for i = 1:length(horizontalsection)
xh(i) = x(5) + horizontalsection(i)*cos(-pi/2+inclh);
yh(i) = TVD(2) - horizontalsection(i)*sin(-pi/2+inclh);
end
figure(2)
plot(x,y,'ko',xv,yv,'k',xb,yb,'k',xt,yt,'k',xb2,yb2,'k',xh,yh,'k','markersize',10)
set(gca,'YDir','reverse')
legend('\fontsize{16}Section markers','\fontsize{16}Ideal well path')
axis([-100 2000 0 700])
grid minor
xlabel('\fontsize{16}Horizontal departure [m]')
ylabel('\fontsize{16}Total vertical depth [m]')
title('\fontsize{16}Ideal well path')
hold on
98 APPENDIX C. BUCKLING DURING INSTALLATION
Contents
INPUTS
Three-dimensional model by Bernt S. Aadnoy
Real well path
clc
clear all
INPUTS
load('TVD.txt')
TVD=TVD';
load('MD.txt')
MD=MD';
load('incl.txt')
incl=incl'*pi/180;
% Pressure Calculations
g = 9.81; % gravity [m/s2]
for i = 1:length(TVD)
Pi(i) = MW*g*TVD(i); % hydrostatic inner pressure [Pa]
Po(i) = Pi(i); % hydrostatic outer pressure [Pa]
end
% Wellpath Calculations
h = zeros(0,length(TVD));
R = zeros(0,length(TVD));
for i = 2:length(TVD)
h(i) = TVD(i-1)-TVD(i); % height of sections [m]
end
for i = 2:length(TVD)
if incl(i)==incl(i-1)
% the section is straight
R(i) = 0; % radius of curvature [m]
else
% the section is curved
R(i)=h(i)/(sin(incl(i-1))-sin(incl(i)));
end
end
100 APPENDIX C. BUCKLING DURING INSTALLATION
L = zeros(1,length(TVD));
for i = 2:length(TVD)
L(i) = MD(i-1)-MD(i) ; % total well length [m]
end
azi = zeros(1,length(TVD));
DL = zeros(1,length(L)); % rad
for i=1:(length(L)-1)
DL(i+1) = acos (sin(incl(i))*sin(incl(i+1))*cos(azi(i)-azi(i+1)) ...
+ cos(incl(i))*cos(incl(i+1)));
end
staticweight = zeros(1,length(TVD));
statichookload = zeros(1,length(TVD));
for i = 1:length(TVD)
staticweight(i) = TVD(i)*BF*wDP;
statichookload(i) = staticweight(1)-staticweight(i);
end
statichookload = statichookload*224.8089; % kN to lb
F = zeros(1,length(L));
hookload = zeros(1,length(my));
Feff = zeros(1,length(L));
F(1) = WOB;
my = 0.5;
for j=1:length(my)
for i=2:length(L)
if incl(i)-incl(i-1) == 0
% the section is straight
F(i) = F(i-1) + BF*wDP*L(i)*(cos(incl(i))-my(j)*sin(incl(i)));
else % the section is curved
F(i) = F(i-1)*exp(-my(j)*abs(DL(i)-DL(i-1)))+ BF*L(i)*wDP*( sin(incl(i))-sin(in
cl(i-1)) )/(incl(i)-incl(i-1));
end
end
figure(j)
plot(F*224.8089,MD,'k',statichookload,MD,'k-.') % convert from [kN] to [lb]
xL = xlim;
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
set(gca,'YDir','reverse')
axis([-20000 20000 0 MD(1)])
grid minor
title('\fontsize{16}Axial Force when lowering String')
legend(['\fontsize{12}Dynamic Hook Load \mu =' num2str(my(j))] ...
,'\fontsize{12}Static Hook Load', 'location','southeast');
xlabel('\fontsize{16}F [lb]')
ylabel('\fontsize{16}MD [m]')
end
load('HD.txt')
HD=HD';
figure(2)
plot(markersx,markersy,'r*',HD,TVD,'r','markersize',10)
legend('\fontsize{16}Section markers','\fontsize{16}Real well path')
set(gca,'YDir','reverse')
axis([-100 2000 0 700])
grid minor
102 APPENDIX C. BUCKLING DURING INSTALLATION
Contents
INPUTS
Three-dimensional model by Bernt S. Aadnoy
Buckling Limits
Permanent Corkscrewing
clc
clear all
INPUTS
load('TVD.txt')
TVD=TVD';
load('MD.txt')
MD=MD';
load('incl.txt')
incl=incl'*pi/180;
% Pressure Calculations in SI
g = 9.81; % gravity [m/s2]
Pi = zeros(1,length(TVD));
Po = zeros(1,length(TVD));
for i = 1:length(TVD)
Pi(i) = MW*g*TVD(i); % hydrostatic inner pressure [kPa]
Po(i) = Pi(i); % hydrostatic outer pressure [kPa]
end
% Wellpath Calculations
h = zeros(0,length(TVD));
R = zeros(0,length(TVD));
for i = 2:length(TVD)
h(i) = TVD(i-1)-TVD(i); % height of sections [m]
end
for i = 2:length(TVD)
if incl(i)==incl(i-1)
% the section is straight
R(i) = 0; % radius of curvature [m]
104 APPENDIX C. BUCKLING DURING INSTALLATION
else
% the section is curved
R(i)=h(i)/(sin(incl(i-1))-sin(incl(i)));
end
end
L = zeros(1,length(TVD));
for i = 2:length(TVD)
L(i) = MD(i-1)-MD(i); % total well length [m]
end
azi = zeros(1,length(TVD));
DL = zeros(1,length(L));
for i=1:(length(L)-1)
DL(i+1) = acos (sin(incl(i))*sin(incl(i+1))*cos(azi(i)-azi(i+1)) ...
+ cos(incl(i))*cos(incl(i+1)));
end
F = zeros(1,length(L));
hookload = zeros(1,length(my));
F(1) = WOB;
Fcsin = zeros(1,length(L));
Fcmin = Fcsin;
Fcmax = Fcsin;
Feff = zeros(1,length(L));
my = 0.45;
for j=1:length(my)
for i=2:length(L)
if incl(i)-incl(i-1) == 0
% the section is straight
F(i) = F(i-1) + BF*wDP*L(i)*(cos(incl(i))-my(j)*sin(incl(i)));
else % the section is curved
F(i) = F(i-1)*exp(-my(j)*abs(DL(i)-DL(i-1)))+ BF*L(i)*wDP*( sin(incl(i))-sin(in
cl(i-1)) )/(incl(i)-incl(i-1));
end
Fa(1) = F(1) + Pi(1)*Ai - Po(1)*Ao;
Fa(i) = F(i) + Pi(1)*Ai - Po(1)*Ao;
Feff(1) = Fa(1) - Pi(1)*Ai + Po(1)*Ao;
Feff(i) = Fa(i) - Pi(i)*Ai + Po(i)*Ao;
end
hookload(j)=F(6);
105
Buckling Limits
%csgID = 9.063;
% Radial clearance
rc = csgID/2-OD/2;
line(xL,[MD(53) MD(53)],'color','black','linestyle',':');
line(xL,[MD(64) MD(64)],'color','black','linestyle',':');
line(xL,[MD(84) MD(84)],'color','black','linestyle',':');
line(xL,[862.58 862.58],'color','black','linestyle',':'); %x-axis
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
set(gca,'YDir','reverse')
grid minor
title('\fontsize{16}Axial Force when lowering string')
legend(['\fontsize{12}Dynamic Hook Load \mu =' num2str(my(j))],...
'\fontsize{12}Effective Axial Load','\fontsize{12}Sinusoidal',...
'\fontsize{12}Minimum Helical', ...
'\fontsize{12}Maximum Helical','location','southeast');
xlabel('\fontsize{16}F [lb]')
ylabel('\fontsize{16}MD [m]')
end
Permanent Corkscrewing
Feff = linspace(-150000,150000,100);
for i = 1:length(Feff)
if Feff(i) < 0
% compression
so(i) = abs(OD*rc*(-Feff(i))./(4*I));
else
% tension
so(i) = 0;
end
end
figure(2)
plot(Feff(1:50),so(1:50),'b',Feffmax,grades(1),'ko','markersize',10)
hold on
xlim([-150000 150000])
ylim([0 140000])
xL = xlim;
yL = ylim;
line(xL, [grades(1) grades(1)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(2) grades(2)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(3) grades(3)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(4) grades(4)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(5) grades(5)],'color','black','linestyle','--'); %x-axis
line([0 0],yL,'color','black')
hold on
grid minor
This appendix includes the MatLab scripts used for compiling results from the method
proposed by Lubinski et al. (1962) and Mitchell et al. (1999).
The first script is a pure reproduction of the model presented in the paper "Helical
buckling of tubing sealed in packers" by Lubinski et al. (1962).
The second script is an implementation of the correlations for buckling strain and buck-
ling length change presented by Mitchell et al. (1999). The script compares the buckling
strain and buckling length change in the model presented by Lubinski et al. (1962) and
Mitchell et al. (1999).
The third script is used for compiling results. The script is modified to apply for the well
investigated, but the procedure of calculation is similar to the first script.
The fourth script is used for a subsea field with wet x-tree. The script is modified by
adding a sea level, and hydrostatic pressure due to seawater.
109
110 APPENDIX D. BUCKLING DURING OPERATIONS
Contents
INPUTS
Length Changes: Packers Permitting Free Motion
Length Changes: Packers Permitting Limited Motion
Length Changes: Packers Permitting No Motion
Plotting
Neutral stability point if the tubing is fixed
clc
clear all
tic
INPUTS
% Tubing Inputs
OD = 2.875; %in outer diameter
weight = 6.5; %lb/ft dry weight
ID = 2.441; %in inner diameter
grades = 10^3.*[40 55 75 80 105]; %psi yield strength
E = 30*10^6; %psi Young's modulus
poisson = 0.3; % Poisson's ratio
Ct = 6.9*10^(-6); %/F coeff. of thermal expansion
Lp = 10000*12; %in tubing length
Dp = 3.25; %in packer bore diameter
% Casing Inputs
csgOD = 7; %in outer diameter
csgweight = 32; %lb/ft dry weight
csgID = 6.094; %in inner diameter
rc = (csgID-OD)/2; %in radial clearance
% Fluid Properties
rhot = 1/231*[7.297 15]; %psi/in tubing initial-final
rhoa = 1/231*[7.297 7.297]; %psi/in annulus initial-final
drhot = rhot(2)-rhot(1); %psi/in tubing change in density
drhoa = rhoa(2)-rhoa(1); %psi/in annulus change in density
wi = rhot.*Ai; %lb/in tubing initial-final
wo = rhoa.*Ao; %lb/in annulus initial-final
w = ws + wi - wo; %lb/in total initial-final
% Pressures
pi = [0 5000]; %psi surface initial-final
po = [0 1000]; %psi surface initial-final
Pi = pi+rhot*Lp; %psi packer initial-final
Po = po+rhoa*Lp; %psi packer initial-final
111
% Temperature
dT = -20; %F avg. change in temperature
% The first length change is deformation due to the true axial force acting
% on the bottom of the tubing. Hooke's law is used:
dL1 = Lp/(E*As)*((Ap-Ao)*dP(2)-(Ap-Ai)*dP(1));
% The second length change occurs due to helical buckling. If the change in
% outer pressure is higher than the change in inner pressure, there will be
% no helical buckling (Lubinski et. al (1962)), hence:
if dP(2)>dP(1)
% No helical buckling
dL2 = 0;
else
% Helical buckling. Calculate length change with a non-linear equation.
dL2 = -rc^2*Ap^2*(dP(2)-dP(1)).^2/(8*E*I*w(2));
% Sinusoidal buckling length change by Mitchell:
% Ff(2) = -Ff(2); % Change signs for Mitchells equation
%dL2 = -(rc)^2/(4*E*I*w(2))*(Ff(2)-Fcr)*(0.3771*Ff(2)-0.3668*Fcr);
% Ff(2) = -Ff(2); % Change back for the remaining model
end
% The fourth length change occurs due to temperature change, also known as
% thermal expansion or thermal contraction.
dL4 = Lp*Ct*dT;
disp(['If the tubing is free to move, the total length change is ', ...
num2str(dL),' in'])
disp(['If the tubing is free to move, the neutral stability point is ', ...
num2str(nfree),' ft from the bottom'])
disp(['If the packer permits limited motion, the total length change is ',...
num2str(dL), ' in'])
If the packer permits limited motion, the total length change is -115.6379 in
dLp = dLf-dL;
if dLp(i) < 0
% Compression
% 'coeff' is the coefficients in the non-linear equation
coeff = [-rc^2/(8*E*I*w(2)) Lp/(E*As) -dLp];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff)';
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar = Ffstarroots(2);
% 'Fp' is the packer force
Fp = Ffstar - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar = Fa(2)+Fp;
if Fp(i) > 0
disp(['The tubing is shortening, and the packer' ...
' force is a tension: ',num2str(Fp(i)),'lb']);
else
disp(['The tubing is lengthening, and the packer' ...
' force is a compression: ',num2str(Fp(i)),'lb']);
end
else
% Tension
% 'coeff' is the coefficients in the linear eq: Hooke's law
coeff = [0 Lp/(E*As) -dLp];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff);
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar = Ffstarroots(2);
% 'Fp' is the packer force
Fp = Ffstar - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar = Fa(2)+Fp;
if Fp(i) > 0
disp(['The tubing is shortening, and the packer' ...
' force is a tension: ',num2str(Fp(i)),'lb']);
else
disp(['The tubing is lengthening, and the packer' ...
' force is a compression: ',num2str(Fp(i)),'lb']);
end
end
end
Plotting
figure(1)
plot(F,dLfnutt,'r',Ffstar,dLp,'r*',Ff(2),dLf,'ro','markersize',10)
hold on
xL = xlim;
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
line(xL, [0 0],'color','black'); %x-axis
hold on
line([Ff(2) Ff(2)],yL,'color','black','linestyle','-.')
114 APPENDIX D. BUCKLING DURING OPERATIONS
n = (Ffstar./w(2))/12; %ft
disp(['If the tubing is fixed, the neutral stability point is ', ...
num2str(n),' ft from the bottom'])
toc
If the tubing is fixed, the neutral stability point is -3965.0715 ft from the bottom
Elapsed time is 2.553992 seconds.
Contents
INPUTS
Lubinski Buckling Strain
Mitchell Buckling Strain
Plotting
Lubinski Buckling Length Change
Mitchell Buckling Length Change
Plotting
clc
clear all
INPUTS
% Tubing Inputs
OD = 2.875; %in outer diameter
weight = 6.5; %lb/ft dry weight
ID = 2.441; %in inner diameter
grades = 10^3.*[40 55 75 80 105]; %psi yield strength
E = 30*10^6; %psi Young's modulus
poisson = 0.3; % Poisson's ratio
Ct = 6.9*10^(-6); %/F coeff. of thermal expansion
Lp = 10000*12; %in tubing length
Dp = 3.25; %in packer bore diameter
% Casing Inputs
csgOD = 7; %in outer diameter
csgweight = 32; %lb/ft dry weight
csgID = 6.094; %in inner diameter
rc = (csgID-OD)/2; %in radial clearance
% Fluid Properties
rhot = 1/231*[7.297 15]; %psi/in tubing initial-final
rhoa = 1/231*[7.297 7.297]; %psi/in annulus initial-final
drhot = rhot(2)-rhot(1); %psi/in tubing change in density
drhoa = rhoa(2)-rhoa(1); %psi/in annulus change in density
wi = rhot.*Ai; %lb/in tubing initial-final
wo = rhoa.*Ao; %lb/in annulus initial-final
w = ws + wi - wo; %lb/in total initial-final
116 APPENDIX D. BUCKLING DURING OPERATIONS
pi = 3.1416;
% Fictitious Force
% For the general case of arbitrary variation of F (Mitchell et al)
Feff = Fcr*linspace(0,4.50,1000);
epsilonLUB = zeros(1,length(Feff));
for i = 1:length(Feff)
if Feff(i) > 0
% Helical buckling when the fictitious force is in compression (+)
epsilonLUB(i) = -(rc)^2/(4*E*I)*Feff(i);
else
% No buckling
epsilonLUB(i) = 0;
end
end
epsilonMIT = zeros(1,length(Feff));
for i = 1:length(Feff)
if Feff(i) < Fcr
% No buckling
epsilonMIT(i) = 0;
elseif Feff(i) < 2.8*Fcr
% Sinusoidal Buckling
epsilonMIT(i) = -0.7285*(rc)^2/(4*E*I)*(Feff(i)).^(0.08)*(Feff(i)-Fcr)^(.92);
else
% Helical Buckling
epsilonMIT(i) = -(rc)^2/(4*E*I)*Feff(i);
end
end
Plotting
figure(1)
plot(Feff/Fcr,-epsilonLUB*4*E*I/rc^2,'k--',...
Feff/Fcr,-epsilonMIT*4*E*I/rc^2,'k')
legend('\fontsize{16}Lubinski',...
'\fontsize{16}Mitchell Correlation','location','northwest')
grid minor
xlabel('\fontsize{16}F/Fcr')
ylabel('\fontsize{16}4EIe/r2')
title('\fontsize{16}Buckling Strain Comparison')
117
yL = ylim;
% Define PD-limit
line([1 1],yL,'color','black','linestyle',':');
% Define helical buckling limit (Mitchell)
line([2.8 2.8],yL,'color','black','linestyle',':');
dLbLUB = zeros(1,length(Feff));
for i = 1:length(Feff)
if Feff(i) > 0
% Helical buckling when the fictitious force is in compression (+)
dLbLUB(i) = -(rc)^2/(8*E*I*w(2))*Feff(i)^2;
else
% No buckling
dLbLUB(i) = 0;
end
end
dLbMIT = zeros(1,length(Feff));
for i = 1:length(Feff)
if Feff(i) < Fcr
% No buckling
dLbMIT(i) = 0;
118 APPENDIX D. BUCKLING DURING OPERATIONS
end
end
Plotting
figure(2)
plot(Feff/Fcr,-dLbLUB*8*E*I*w(2)/rc^2,'k--',...
Feff/Fcr,-dLbMIT*8*E*I*w(2)/rc^2,'k')
legend('\fontsize{16}Lubinski',...
'\fontsize{16}Mitchell Correlation','location','northwest')
grid minor
xlabel('\fontsize{16}F/Fcr')
ylabel('\fontsize{16}-8EIwdLb/r2')
title('\fontsize{16}Buckling Length Change Comparison')
yL = ylim;
% Define PD-limit
line([1 1],yL,'color','black','linestyle',':');
% Define helical buckling limit (Mitchell)
line([2.8 2.8],yL,'color','black','linestyle',':');
Contents
INPUTS
Length changes and packer forces
Critical buckling limits
Length Changes: Packers Permitting Free Motion
Length Changes: Packers Permitting Limited Motion
Length Changes: Packers Permitting No Motion
Plotting with constant slack-off/pick-up
Plotting with variating slack-off/pick-up
Neutral Stability Point
Plotting axial loads and neutral stability point
Permanent Corkscrewing
Plotting with constant slack-off/pick-up
Plotting with variating slack-off/pick-up
Plotting axial loads and neutral stability point with critical buckling
Including locally imposed bending loads on axial loads
Plotting bending stress
clc
clear all
INPUTS
tic
% Tubing Inputs
OD = 3.5; %in outer diameter
weight = 9.2; %lb/ft dry weight
ID = 2.992; %in inner diameter
grades = 10^3.*[40 55 75 80 105]; %psi yield strength
E = 30*10^6; %psi Young's modulus
poisson = 0.3; % Poisson's ratio
Ct = 6.9*10^(-6); %/F coeff. of thermal expansion
% Total length of well
Lp = 2523.1; %m value from Compass
Lp = Lp*39.37; %in converting from m to in
% Total vertical depth
TVD = 669.62; %m value from Compass
TVD = TVD*39.37; %in converting from m to in
% Tubing calculations
Ao = pi/4*OD.^2; %in2 outer area
Ai = pi/4*ID.^2; %in2 inner area
As = Ao-Ai; %in2 cross-sectional area
Ap = Ai; %in2 packer bore area PBR
R = OD./ID; %ratio OD/ID-ratio
ws = weight./12; %lb/in dry weight of tubing
I = (pi/64)*(OD.^4-ID.^4); %in4 moment of inertia
% Casing Inputs
csgID = 5.921; %in inner diameter
rc = (csgID-OD)/2; %in radial clearance
120 APPENDIX D. BUCKLING DURING OPERATIONS
% Fluid Properties
rhosw = 1/231*8.68; %psi/in sea water density
rhog = 1/231*1; %psi/in gas density
rhoo = 1/231*7.68; %psi/in oil density
% Pressures
Pr = 1015.3; %psi 70 bar reservoir pressure
Pg = rhog*TVD; %psi gas-filled well
SIWHP = Pr - Pg; %psi shut-in wellhead pressure
WDP = SIWHP + 508; %psi 35 bar safety margin
Ppacker = 4278.6; %psi 295 bar set packer pressure
% Density calculations
drhot = rhot(2)-rhot(1); %psi/in change in density
drhoa = rhoa(2)-rhoa(1); %psi/in change in density
wi = rhot.*Ai; %lb/in tubing initial-final
wo = rhoa.*Ao; %lb/in annulus initial-final
w = ws + wi - wo; %lb/in total initial-final
% Pressure calculations
Pi = pi+rhot*TVD; %psi packer initial-final
Po = po+rhoa*TVD; %psi packer initial-final
dp = [pi(2)-pi(1) po(2)-po(1)]; %psi change in surface pressure
dP = [Pi(2)-Pi(1) Po(2)-Po(1)]; %psi change in packer pressure
disp(['The length change related to this force is ', num2str(dLf), ' in'])
% The first length change is deformation due to the true axial force acting
% on the bottom of the tubing. Hooke's law is used:
dL1 = Lp/(E*As)*((Ap-Ao)*dP(2)-(Ap-Ai)*dP(1));
% The second length change occurs due to helical buckling. If the change in
% outer pressure is higher than the change in inner pressure, there will be
% no helical buckling (Lubinski et. al (1962)), hence:
if dP(2)>dP(1)
% No helical buckling
dL2 = 0;
else
% Helical buckling. Calculate length change with a non-linear equation.
dL2 = -rc^2*Ap^2*(dP(2)-dP(1)).^2/(8*E*I*w(2));
% Sinusoidal buckling length change by Mitchell:
% Ff(2) = -Ff(2); % Change signs for Mitchells equation
% dL2 = -(rc)^2/(4*E*I*w(2))*(Ff(2)-Fcr(2))*(0.3771*Ff(2)-0.3668*Fcr(2));
% Ff(2) = -Ff(2); % Change back for the remaining model
end
% Calculate radial pressure due to ballooning and fluid flow
%delta = 0; % Pressure drop per unit length
%dL3 = -(poisson/E)*((drhot-R^2*drhoa-(1+2*poisson)/(2*poisson)*delta)*Lp^2/(R^2-1)).
..
% -(2*poisson/E)*((dp(1)-R^2*dp(2))*Lp/(R^2-1));
122 APPENDIX D. BUCKLING DURING OPERATIONS
% The fourth length change occurs due to temperature change, also known as
% thermal expansion or thermal contraction.
dL4 = Lp*Ct*dT;
disp(['If the tubing is free to move, the total length change is ', ...
num2str(dL),' in'])
Fso = 0; %lb
disp(['If the packer permits limited motion, the total length change is ',...
num2str(dL), ' in'])
If the packer permits limited motion, the total length change is -19.1516 in
dLp = dLf-dL;
% Define vectors
Ffstar = zeros(1,length(Fso));
Fp = zeros(1,length(Fso));
Fastar = zeros(1,length(Fso));
for i = 1:length(dLp)
if dLp(i) < 0
% Compression
% 'coeff' is the coefficients in the non-linear equation
coeff = [-rc^2/(8*E*I*w(2)) Lp/(E*As) -dLp(i)];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff)';
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar(i) = Ffstarroots(2);
% 'Fp' is the packer force
Fp(i) = Ffstar(i) - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar(i) = Fp(i)+Fa(2);
if Fp(i) > 0
disp(['The tubing is shortening, and the packer' ...
' force is a tension: ',num2str(Fp(i)),'lb']);
else
disp(['The tubing is lengthening, and the packer' ...
' force is a compression: ',num2str(Fp(i)),'lb']);
end
else
% Tension
% 'coeff' is the coefficients in the linear eq: Hooke's law
coeff = [0 Lp/(E*As) -dLp(i)];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff);
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar(i) = Ffstarroots;
% 'Fp' is the packer force
Fp(i) = Ffstar(i) - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar(i) = Fp(i)+Fa(2);
if Fp(i) > 0
disp(['The tubing is shortening, and the packer' ...
' force is a tension: ',num2str(Fp(i)),'lb']);
else
disp(['The tubing is lengthening, and the packer' ...
' force is a compression: ',num2str(Fp(i)),'lb']);
end
end
end
124 APPENDIX D. BUCKLING DURING OPERATIONS
figure(1)
plot(F,dLfnutt,'r',Ffstar,dLp,'r*',Ff(2),dLf,'ro',Ffstar,dLf,'b>',...
'markersize',10)
hold on
xL = xlim;
yL = ylim;
line([Ff(2) Ffstar],[dLf dLf],'linestyle','-')
line([0 0], yL,'color','black'); %y-axis
line(xL, [0 0],'color','black'); %x-axis
hold on
line([Ff(2) Ff(2)],yL,'color','black','linestyle','-.')
line(xL, [dLf dLf],'color','black','linestyle','-.')
hold on
line([Ffstar Ffstar],yL,'color','black','linestyle','--')
line(xL, [dLp dLp],'color','black','linestyle','--')
grid minor
xlabel('\fontsize{16}Axial Force [lb]')
ylabel('\fontsize{16}dL [in]')
title('\fontsize{16}Axial Force vs. Length Change')
legend('\fontsize{16}Axial Force vs. Length Change' ...
, '\fontsize{16}Ff*','\fontsize{16}Ff','\fontsize{16}Fp', ...
'location','northwest')
125
% Initial Conditions
% Surface (z=0)
% True Axial Load
Fas(1) = TVD*ws + Po(1)*(Ap-Ao) - Pi(1)*(Ap-Ai);
% Effective Axial Load
Ffs(1) = Fas(1) - pi(1)*Ai + po(1)*Ao;
% Packer (z=zshoe)
% True Axial Load
Fap(1) = Po(1)*(Ap-Ao) - Pi(1)*(Ap-Ai);
% Effective Axial Load
% Assumes that the pressure is bled of after setting the packer, hence Fp=0
Ffp(1) = Fap(1) - Pi(1)*Ai + Po(1)*Ao;
126 APPENDIX D. BUCKLING DURING OPERATIONS
% Surface (z=0)
% The packer force will adjuct the true axial load, and this effect will
% also be included when calculating the effective axial load. Lubinski used
% a star "*" to describe these loads when the tubing is anchored.
% True Axial Load
Fas(2) = TVD*ws + Po(2)*(Ap-Ao) - Pi(2)*(Ap-Ai) + Fp;
% Effective Axial Load
Ffs(2) = Fas(2) - pi(2)*Ai + po(2)*Ao;
% Packer (z=zshoe)
% True Axial Load
Fap(2) = Po(2)*(Ap-Ao) - Pi(2)*(Ap-Ai) + Fp;
% Effective Axial Load
Ffp(2) = Fap(2) - Pi(2)*Ai + Po(2)*Ao;
figure(2)
plot([Fas(1) Fap(1)],[0 TVD]/12,'k', ...
[Ffs(1) Ffp(1)],[0 TVD]/12,'r', ...
[Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
0, TVD/12+n,'ro')
ylim([0 TVD/12]);
yL = [0 TVD/12];
xL = xlim;
line([0 0], yL,'color','black');
line(xL,[0 0],'color','black');
line(xL,[TVD/12 TVD/12],'color','black');
legend('\fontsize{16}Fa Initial Conditions', ...
'\fontsize{16}Feff Initial Conditions', ...
'\fontsize{16}Fa Load Case', ...
'\fontsize{16}Feff Load Case', ...
'\fontsize{16}NSP','location','northwest')
set(gca,'YDir','reverse')
grid minor
ylabel('\fontsize{16}Depth [ft]')
yyaxis right
ylim([0 670])
line(xL,[450 450],'color','black','linestyle',':');
line(xL,[650 650],'color','black','linestyle',':');
set(gca,'YDir','reverse')
xlabel('\fontsize{16}Force [lb]')
ylabel('\fontsize{16}Depth [m]')
title('\fontsize{16}True vs. Effective Axial Force')
127
Permanent Corkscrewing
for i = 1:length(Fso)
if Fso(i) < 0
soinitial(i) = abs(Fso(i)/As + OD*rc*Fso(i)./(4*I));
end
end
for i = 1:length(Ffstar)
if Ffstar(i) < 0
% No buckling
sigmab = 0;
% elseif Ffstar(i) < 5.55*Fcr && Ffstar(i) > 1.94*Fcr
% Sinusoidal Buckling
% sigmab = 0.3151*OD*rc/I*Ffstar(i).^(0.08)*(Ffstar(i)-Fcr).^(0.92);
% elseif Ffstar(i) > 5.55*Fcr
128 APPENDIX D. BUCKLING DURING OPERATIONS
else
% Helical Buckling
sigmab = OD*rc/(4*I).*Ffstar(i);
end
sigmaa = Fastar(i)/As;
end
figure(3)
plot(Fso,soinitial,'bo',Fso,so,'ro',Fso,si,'yo','markersize',10)
hold on
xL = xlim;
line(xL, [grades(1) grades(1)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(2) grades(2)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(3) grades(3)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(4) grades(4)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(5) grades(5)],'color','black','linestyle','--'); %x-axis
%yL = [0 120000];
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
hold on
grid minor
xlabel('\fontsize{16}Initial Slack-Off and Pick-Up Force [lb]')
ylabel('\fontsize{16}Triaxial Stress [psi]')
title('\fontsize{16}Permanent Corkscrewing')
legend('\fontsize{16}Sigmao Initial','\fontsize{16}Sigmao',...
'\fontsize{16}Sigmai','location','southeast')
129
% figure(3)
% plot(Fso(1:50),soinitial(1:50),'b',Fso,so,'r',Fso,si,'y','markersize',10)
% hold on
% xL = xlim;
% % line(xL, [grades(1) grades(1)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(2) grades(2)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(3) grades(3)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(4) grades(4)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(5) grades(5)],'color','black','linestyle','--'); %x-axis
% %yL = [0 120000];
% yL = ylim;
% line([0 0], yL,'color','black'); %y-axis
% hold on
% grid minor
% xlabel('\fontsize{16}Initial Slack-Off and Pick-Up Force [lb]')
% ylabel('\fontsize{16}Triaxial Stress [psi]')
% title('\fontsize{16}Permanent Corkscrewing')
% legend('\fontsize{16}Sigmao Initial','\fontsize{16}Sigmao',...
% '\fontsize{16}Sigmai','location','southeast')
Plotting axial loads and neutral stability point with critical buckling
figure(4)
plot([Fas(1) Fap(1)],[0 TVD]/12,'k', ...
[Ffs(1) Ffp(1)],[0 TVD]/12,'r', ...
[Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
0, TVD/12+n(2),'ro', ...
-[Fcsinv Fcsinv],[0 450*39.37]/12,'g', ...
130 APPENDIX D. BUCKLING DURING OPERATIONS
TVDvector = linspace(0,TVD/12,1000);
Favector = linspace(Fas(2),Fap(2),1000);
Feffvector = linspace(Ffs(2),Ffp(2),1000);
% Define compression as positive
Favector = -Favector;
Feffvector = -Feffvector;
Fbendvector = zeros(1,length(Favector));
for i = 1:length(Feffvector)
if TVDvector(i) < 450*39.37/12
if Feffvector(i) < Fcsinv
% No buckling
Fbendvector(i) = Favector(i);
elseif Feffvector(i) < Fchelv
% Sinusoidal buckling
if Favector(i) < 0
% Tubing is in tension
Fbendvector(i) = Favector(i) ...
- 0.3151*OD*rc/I*Feffvector(i).^(0.08)*(Feffvector(i)-Fcsinv).^(0.92);
else
% Tubing is in compression
Fbendvector(i) = Favector(i) ...
+ 0.3151*OD*rc/I*Feffvector(i).^(0.08)*(Feffvector(i)-Fcsinv).^(0.92);
end
else
% Helical buckling
if Favector(i) < 0
% Tubing is in tension
Fbendvector(i) = Favector(i) - OD*rc/(4*I).*Feffvector(i);
else
% Tubing is in compression
Fbendvector(i) = Favector(i) + OD*rc/(4*I).*Feffvector(i);
end
end
else
Fbendvector(i) = Favector(i);
end
end
figure(5)
plot([Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
Fbendvector,TVDvector,'k-.',...
0, TVD/12+n(2),'ro', ...
-[Fcsinv Fcsinv],[0 450*39.37]/12,'g', ...
-[Fchelv Fchelv],[0 450*39.37]/12,'b', ...
-Fcrsindev,TVDbuild*39.37/12,'g', ...
-Fcrheldev,TVDbuild*39.37/12,'b', ...
[-Fcrsinhor -Fcrsinhor],[650 669.62]*39.37/12,'g', ...
[-Fcrhelhor -Fcrhelhor],[650 669.62]*39.37/12,'b')
ylim([0 TVD/12]);
yL = [0 TVD/12];
132 APPENDIX D. BUCKLING DURING OPERATIONS
xL = xlim;
line([0 0], yL,'color','black');
line(xL,[0 0],'color','black');
line(xL,[TVD/12 TVD/12],'color','black');
legend('\fontsize{16}Fa Load Case', ...
'\fontsize{16}Feff Load Case', ...
'\fontsize{16}Fa with Bending', ...
'location','northwest') %, ...
% '\fontsize{12}Fa Load Case', ...
% '\fontsize{12}Feff Load Case', ...
% '\fontsize{12}NSP', ...
% '\fontsize{12}Sinusoidal buckling', ...
% '\fontsize{12}Helical buckling', ...
% 'location','northwest')
set(gca,'YDir','reverse')
grid minor
ylabel('\fontsize{16}Depth [ft]')
yyaxis right
ylim([0 670])
line(xL,[450 450],'color','black','linestyle',':');
line(xL,[650 650],'color','black','linestyle',':');
set(gca,'YDir','reverse')
xlabel('\fontsize{16}Force [lb]')
ylabel('\fontsize{16}Depth [m]')
title('\fontsize{16}True vs. Effective Axial Force')
toc
Contents
INPUTS
Length changes and packer forces
Critical buckling limits
Length Changes: Packers Permitting Free Motion
Length Changes: Packers Permitting Limited Motion
Length Changes: Packers Permitting No Motion
Plotting with constant slack-off/pick-up
Plotting with variating slack-off/pick-up
Neutral Stability Point
Plotting axial loads and neutral stability point
Permanent Corkscrewing
Plotting with constant slack-off/pick-up
Plotting with variating slack-off/pick-up
Plotting axial loads and neutral stability point with critical buckling
Including locally imposed bending loads on axial loads
Plotting bending stress
clc
clear all
INPUTS
tic
% Tubing Inputs
OD = 3.5; %in outer diameter
weight = 9.2; %lb/ft dry weight
ID = 2.992; %in inner diameter
grades = 10^3.*[40 55 75 80 105]; %psi yield strength
E = 30*10^6; %psi Young's modulus
poisson = 0.3; % Poisson's ratio
Ct = 6.9*10^(-6); %/F coeff. of thermal expansion
% Total vertical depth
TVDsw = 400; %m sea water depth
TVDsw = TVDsw*39.37; %in converting from m to in
TVD = 269.62; %m vertical well depth
TVD = TVD*39.37; %in converting from m to in
% Total length of well
Lp = 2523.1-TVDsw; %m value from Compass minus
% sea water depth
Lp = Lp*39.37; %in converting from m to in
% Tubing calculations
Ao = pi/4*OD.^2; %in2 outer area
Ai = pi/4*ID.^2; %in2 inner area
As = Ao-Ai; %in2 cross-sectional area
Ap = Ai; %in2 packer bore area PBR
R = OD./ID; %ratio OD/ID-ratio
ws = weight./12; %lb/in dry weight of tubing
I = (pi/64)*(OD.^4-ID.^4); %in4 moment of inertia
% Casing Inputs
135
% Fluid Properties
rhosw = 1/231*8.68; %psi/in sea water density
rhog = 1/231*1; %psi/in gas density
rhoo = 1/231*7.68; %psi/in oil density
% Pressures
Pr = 1015.3; %psi 70 bar reservoir pressure
Pg = rhog*TVD; %psi gas-filled well
SIWHP = Pr - Pg; %psi shut-in wellhead pressure
WDP = SIWHP + 508; %psi 35 bar safety margin
Ppacker = 4278.6; %psi 295 bar set packer pressure
Psw = rhosw*TVDsw; %psi sea water pressure
% Density calculations
drhot = rhot(2)-rhot(1); %psi/in change in density
drhoa = rhoa(2)-rhoa(1); %psi/in change in density
wi = rhot.*Ai; %lb/in tubing initial-final
wo = rhoa.*Ao; %lb/in annulus initial-final
w = ws + wi - wo; %lb/in total initial-final
% Pressure calculations
Pi = pi+rhot*TVD; %psi packer initial-final
Po = po+rhoa*TVD; %psi packer initial-final
136 APPENDIX D. BUCKLING DURING OPERATIONS
disp(['The length change related to this force is ', num2str(dLf), ' in'])
% The first length change is deformation due to the true axial force acting
% on the bottom of the tubing. Hooke's law is used:
dL1 = Lp/(E*As)*((Ap-Ao)*dP(2)-(Ap-Ai)*dP(1));
% The second length change occurs due to helical buckling. If the change in
% outer pressure is higher than the change in inner pressure, there will be
% no helical buckling (Lubinski et. al (1962)), hence:
if dP(2)>dP(1)
% No helical buckling
dL2 = 0;
else
% Helical buckling. Calculate length change with a non-linear equation.
dL2 = -rc^2*Ap^2*(dP(2)-dP(1)).^2/(8*E*I*w(2));
% Sinusoidal buckling length change by Mitchell:
% Ff(2) = -Ff(2); % Change signs for Mitchells equation
% dL2 = -(rc)^2/(4*E*I*w(2))*(Ff(2)-Fcr(2))*(0.3771*Ff(2)-0.3668*Fcr(2));
% Ff(2) = -Ff(2); % Change back for the remaining model
end
% Calculate radial pressure due to ballooning and fluid flow
137
% The fourth length change occurs due to temperature change, also known as
% thermal expansion or thermal contraction.
dL4 = Lp*Ct*dT;
disp(['If the tubing is free to move, the total length change is ', ...
num2str(dL),' in'])
Fso = 0; %lb
disp(['If the packer permits limited motion, the total length change is ',...
num2str(dL), ' in'])
If the packer permits limited motion, the total length change is 91.5191 in
dLp = dLf-dL;
138 APPENDIX D. BUCKLING DURING OPERATIONS
Ffstar = zeros(1,length(Fso));
Fp = zeros(1,length(Fso));
Fastar = zeros(1,length(Fso));
for i = 1:length(dLp)
if dLp(i) < 0
% Compression
% 'coeff' is the coefficients in the non-linear equation
coeff = [-rc^2/(8*E*I*w(2)) Lp/(E*As) -dLp(i)];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff)';
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar(i) = Ffstarroots(2);
% 'Fp' is the packer force
Fp(i) = Ffstar(i) - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar(i) = Fp(i)+Fa(2);
if Fp(i) > 0
disp(['The tubing is shortening, and the packer' ...
' force is a tension: ',num2str(Fp(i)),'lb']);
else
disp(['The tubing is lengthening, and the packer' ...
' force is a compression: ',num2str(Fp(i)),'lb']);
end
else
% Tension
% 'coeff' is the coefficients in the linear eq: Hooke's law
coeff = [0 Lp/(E*As) -dLp(i)];
% 'roots' calculates the F corresponding to dLp
Ffstarroots = roots(coeff);
% 'Ffstar' is the effective axial force when the tubing is fixed
Ffstar(i) = Ffstarroots;
% 'Fp' is the packer force
Fp(i) = Ffstar(i) - Ff(2);
% 'Fastar' is the true axial force when the tubing is fixed
Fastar(i) = Fp(i)+Fa(2);
if Fp(i) > 0
139
figure(1)
plot(F,dLfnutt,'r',Ffstar,dLp,'r*',Ff(2),dLf,'ro',Ffstar,dLf,'b>',...
'markersize',10)
hold on
xL = xlim;
yL = ylim;
line([Ff(2) Ffstar],[dLf dLf],'linestyle','-')
line([0 0], yL,'color','black'); %y-axis
line(xL, [0 0],'color','black'); %x-axis
hold on
line([Ff(2) Ff(2)],yL,'color','black','linestyle','-.')
line(xL, [dLf dLf],'color','black','linestyle','-.')
hold on
line([Ffstar Ffstar],yL,'color','black','linestyle','--')
line(xL, [dLp dLp],'color','black','linestyle','--')
grid minor
xlabel('\fontsize{16}Axial Force [lb]')
ylabel('\fontsize{16}dL [in]')
title('\fontsize{16}Axial Force vs. Length Change')
legend('\fontsize{16}Axial Force vs. Length Change' ...
, '\fontsize{16}Ff*','\fontsize{16}Ff','\fontsize{16}Fp', ...
'location','northwest')
140 APPENDIX D. BUCKLING DURING OPERATIONS
% Initial Conditions
% Surface (z=0)
% True Axial Load
Fas(1) = TVD*ws + Po(1)*(Ap-Ao) - Pi(1)*(Ap-Ai);
% Effective Axial Load
Ffs(1) = Fas(1) - pi(1)*Ai + po(1)*Ao;
% Packer (z=zshoe)
% True Axial Load
Fap(1) = Po(1)*(Ap-Ao) - Pi(1)*(Ap-Ai);
% Effective Axial Load
% Assumes that the pressure is bled of after setting the packer, hence Fp=0
Ffp(1) = Fap(1) - Pi(1)*Ai + Po(1)*Ao;
141
% Surface (z=0)
% The packer force will adjuct the true axial load, and this effect will
% also be included when calculating the effective axial load. Lubinski used
% a star "*" to describe these loads when the tubing is anchored.
% True Axial Load
Fas(2) = TVD*ws + Po(2)*(Ap-Ao) - Pi(2)*(Ap-Ai) + Fp;
% Effective Axial Load
Ffs(2) = Fas(2) - pi(2)*Ai + po(2)*Ao;
% Packer (z=zshoe)
% True Axial Load
Fap(2) = Po(2)*(Ap-Ao) - Pi(2)*(Ap-Ai) + Fp;
% Effective Axial Load
Ffp(2) = Fap(2) - Pi(2)*Ai + Po(2)*Ao;
figure(2)
plot([Fas(1) Fap(1)],[0 TVD]/12,'k', ...
[Ffs(1) Ffp(1)],[0 TVD]/12,'r', ...
[Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
0, TVD/12+n,'ro')
ylim([0 TVD/12]);
yL = [0 TVD/12];
xL = xlim;
line([0 0], yL,'color','black');
line(xL,[0 0],'color','black');
line(xL,[TVD/12 TVD/12],'color','black');
legend('\fontsize{16}Fa Initial Conditions', ...
'\fontsize{16}Feff Initial Conditions', ...
'\fontsize{16}Fa Load Case', ...
'\fontsize{16}Feff Load Case', ...
'\fontsize{16}NSP','location','northwest')
set(gca,'YDir','reverse')
grid minor
ylabel('\fontsize{16}Depth [ft]')
yyaxis right
ylim([0 270])
line(xL,[50 50],'color','black','linestyle',':');
line(xL,[250 250],'color','black','linestyle',':');
set(gca,'YDir','reverse')
xlabel('\fontsize{16}Force [lb]')
ylabel('\fontsize{16}Depth [m]')
title('\fontsize{16}True vs. Effective Axial Force')
142 APPENDIX D. BUCKLING DURING OPERATIONS
Permanent Corkscrewing
soinitial = zeros(1,length(Fso));
for i = 1:length(Fso)
if Fso(i) < 0
soinitial(i) = abs(Fso(i)/As + OD*rc*Fso(i)./(4*I));
end
end
so = zeros(1,length(Ffstar));
si = so;
Ffstar = -Ffstar;
Fastar = -Fastar;
for i = 1:length(Ffstar)
if Ffstar(i) < 0
% No buckling
sigmab = 0;
% elseif Ffstar(i) < 5.55*Fcr && Ffstar(i) > 1.94*Fcr
% Sinusoidal Buckling
% sigmab = 0.3151*OD*rc/I*Ffstar(i).^(0.08)*(Ffstar(i)-Fcr).^(0.92);
% elseif Ffstar(i) > 5.55*Fcr
else
% Helical Buckling
sigmab = OD*rc/(4*I).*Ffstar(i);
end
143
sigmaa = Fastar(i)/As;
end
figure(3)
plot(Fso,soinitial,'bo',Fso,so,'ro',Fso,si,'yo','markersize',10)
hold on
xL = xlim;
line(xL, [grades(1) grades(1)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(2) grades(2)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(3) grades(3)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(4) grades(4)],'color','black','linestyle','--'); %x-axis
line(xL, [grades(5) grades(5)],'color','black','linestyle','--'); %x-axis
%yL = [0 120000];
yL = ylim;
line([0 0], yL,'color','black'); %y-axis
hold on
grid minor
xlabel('\fontsize{16}Initial Slack-Off and Pick-Up Force [lb]')
ylabel('\fontsize{16}Triaxial Stress [psi]')
title('\fontsize{16}Permanent Corkscrewing')
legend('\fontsize{16}Sigmao Initial','\fontsize{16}Sigmao',...
'\fontsize{16}Sigmai','location','southeast')
144 APPENDIX D. BUCKLING DURING OPERATIONS
% figure(3)
% plot(Fso(1:50),soinitial(1:50),'b',Fso,so,'r',Fso,si,'y','markersize',10)
% hold on
% xL = xlim;
% % line(xL, [grades(1) grades(1)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(2) grades(2)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(3) grades(3)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(4) grades(4)],'color','black','linestyle','--'); %x-axis
% % line(xL, [grades(5) grades(5)],'color','black','linestyle','--'); %x-axis
% %yL = [0 120000];
% yL = ylim;
% line([0 0], yL,'color','black'); %y-axis
% hold on
% grid minor
% xlabel('\fontsize{16}Initial Slack-Off and Pick-Up Force [lb]')
% ylabel('\fontsize{16}Triaxial Stress [psi]')
% title('\fontsize{16}Permanent Corkscrewing')
% legend('\fontsize{16}Sigmao Initial','\fontsize{16}Sigmao',...
% '\fontsize{16}Sigmai','location','southeast')
Plotting axial loads and neutral stability point with critical buckling
figure(4)
plot([Fas(1) Fap(1)],[0 TVD]/12,'k', ...
[Ffs(1) Ffp(1)],[0 TVD]/12,'r', ...
[Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
0, TVD/12+n(2),'ro', ...
-[Fcsinv Fcsinv],[0 50*39.37]/12,'g', ...
145
TVDvector = linspace(0,TVD/12,1000);
Favector = linspace(Fas(2),Fap(2),1000);
Feffvector = linspace(Ffs(2),Ffp(2),1000);
% Define compression as positive
Favector = -Favector;
Feffvector = -Feffvector;
Fbendvector = zeros(1,length(Favector));
for i = 1:length(Feffvector)
if TVDvector(i) < 50*39.37/12
if Feffvector(i) < Fcsinv
% No buckling
Fbendvector(i) = Favector(i);
elseif Feffvector(i) < Fchelv
% Sinusoidal buckling
if Favector(i) < 0
% Tubing is in tension
Fbendvector(i) = Favector(i) ...
- 0.3151*OD*rc/I*Feffvector(i).^(0.08)*(Feffvector(i)-Fcsinv).^(0.92);
else
% Tubing is in compression
Fbendvector(i) = Favector(i) ...
+ 0.3151*OD*rc/I*Feffvector(i).^(0.08)*(Feffvector(i)-Fcsinv).^(0.92);
end
else
% Helical buckling
if Favector(i) < 0
% Tubing is in tension
Fbendvector(i) = Favector(i) - OD*rc/(4*I).*Feffvector(i);
else
% Tubing is in compression
Fbendvector(i) = Favector(i) + OD*rc/(4*I).*Feffvector(i);
end
end
else
Fbendvector(i) = Favector(i);
end
end
figure(5)
plot([Fas(2) Fap(2)],[0 TVD]/12,'k--', ...
[Ffs(2) Ffp(2)],[0 TVD]/12,'r--', ...
Fbendvector,TVDvector,'k-.',...
0, TVD/12+n(2),'ro', ...
-[Fcsinv Fcsinv],[0 50*39.37]/12,'g', ...
-[Fchelv Fchelv],[0 50*39.37]/12,'b', ...
-Fcrsindev,TVDbuild*39.37/12,'g', ...
-Fcrheldev,TVDbuild*39.37/12,'b', ...
[-Fcrsinhor -Fcrsinhor],[250 269.62]*39.37/12,'g', ...
[-Fcrhelhor -Fcrhelhor],[250 269.62]*39.37/12,'b')
ylim([0 TVD/12]);
yL = [0 TVD/12];
147
xL = xlim;
line([0 0], yL,'color','black');
line(xL,[0 0],'color','black');
line(xL,[TVD/12 TVD/12],'color','black');
legend('\fontsize{16}Fa Load Case', ...
'\fontsize{16}Feff Load Case', ...
'\fontsize{16}Fa with Bending', ...
'location','northwest') %, ...
% '\fontsize{12}Fa Load Case', ...
% '\fontsize{12}Feff Load Case', ...
% '\fontsize{12}NSP', ...
% '\fontsize{12}Sinusoidal buckling', ...
% '\fontsize{12}Helical buckling', ...
% 'location','northwest')
set(gca,'YDir','reverse')
grid minor
ylabel('\fontsize{16}Depth [ft]')
yyaxis right
ylim([0 270])
line(xL,[50 50],'color','black','linestyle',':');
line(xL,[250 250],'color','black','linestyle',':');
set(gca,'YDir','reverse')
xlabel('\fontsize{16}Force [lb]')
ylabel('\fontsize{16}Depth [m]')
title('\fontsize{16}True vs. Effective Axial Force')
toc