18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

Modlisation tridimensionnelle du flambage des tiges de forage

Jol AKOWANOU, Hedi SELLAMI, Stphane MENAND et Michel TIJANI
Centre de Geosciences, Groupe Hydro-Go-Ingnierie,
Ecole des Mines de Paris
35 rue saint-honor, 77305 Fontainebleau - France
[email protected]

Rsum
Le flambage des structures lances est un thme qui fait toujours lobjet de nombreuses tudes. On traite ici du
cas particulier des longs tubes minces en acier utiliss dans les forages ptroliers. Un code de calcul unique
permettant daborder pour la premire fois notre connaissance, la problmatique du flambage des tiges de
forage dans un puits tridimensionnel, en tenant compte du frottement et de la rotation des tiges a t labor. Ce
code de calcul a t valid sur des donnes relles de forage de puits ptroliers ainsi qu partir dessais sur un
dispositif exprimental chelle rduite, reproduisant le phnomne du flambage des tiges dans un trou de
forage prsentant ou non de la tortuosit. Ce code de calcul devrait contribuer mieux prdire loccurrence du
flambage lors des oprations de forage, phnomne qui peut tre lorigine de nombreux dysfonctionnements,
comme la rupture par fatigue des tiges.

Abstract :
Buckling of long slender structures is still the subject of many studies. Here we deal with the long thin steel
tubes used in the oilfield drilling. For the first time to our knowledge a single computer code which allows to
tackle the problem of drillpipe buckling in three-dimensional wells, by taking into account the friction and the
rotation of the drillpipe has been developped. A comparison of the model is proposed with an experimental set
up that enables to reproduce the wellbore tortuosity. This computer code is then validated on real data of drilling,
as well as on lab tests with a reduced scale experimental facility which reproduces the buckling phenomena of
drillpipe in actual well conditions. This model should contribute to better predict the occurrence of buckling on
drilling operations, phenomena that can be in the beginning of many dysfunction, like the fatigue failure of
drillstring.

Mots-clefs :
Flambage ; tige de forage ; frottement ; rotation ; fatigue ;
1. Introduction
Aujourdhui, les puits devenant de plus en plus complexes, avoir un modle danalyse raliste
du flambage qui peut tre utilis pour la planification et la ralisation de ces puits ptroliers est un
atout.
Plusieurs mthodes de rsolution du flambage ont t labores ces dernires annes. Ces
approches, quelles soient analytiques, numriques ou exprimentales ignorent souvent leffet des
frottements tige et paroi du puits, ne tiennent pas compte de la rotation des tiges ou de la tortuosit du
puits de forage. Le milieu industriel aborde gnralement le problme du flambage en calculant
uniquement la charge critique (de compression) partir de laquelle les tiges flambent. Lambition de
ce papier est de contribuer une bonne comprhension et une meilleure analyse du flambage dans
les conditions relles de forage, cest dire en tenant compte du frottement, de la rotation et de la
tortuosit de la trajectoire. Le code de calcul dvelopp sera dcrit dans la premire partie de cet
article, puis le dispositif exprimental permettant de valider le code de calcul sera prsent. Enfin on
montrera dans la dernire partie de larticle, une confrontation entre ces simulations numriques et
lexprience dmontrant la robustesse du code de calcul.

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

2. Modle 3D de flambage des tiges de forage

2.1. Dfinition du problme
La distribution de la force axiale le long dune tige de forage met en vidence le risque de
flambement de la partie comprime (fig. 1-a). Ainsi lorsque la compression subie par la tige dpasse
une valeur critique, la tige flambe sous une forme sinusodale ou hlicodale (fig. 1-b).
En effet, on parlera du flambage sinusodal, quand la tige prend une configuration sinusodale,
lorsque la compression est plus grande quun premier effort critique (note Fcr _ sin ). Paslay et al
(1984) en se basant sur les travaux de Dawson (1964) ont trouv que la force critique de flambage
sinusodal dune tige dans un puits inclin vaut :

Fcr _ sin = 2

EIw sin(inc)
r

(1)

O EI est la rigidit de la tige, w le poids linique djaug de la tige, inc linclinaison du puits et r le
jeu entre la tige et le puits.
Dans cette configuration de tige flambe, les points de contact, voire les lignes de contact de la
tige avec la paroi font apparatre des efforts latraux dit de contact, qui gnrent une augmentation des
forces de frottement. Lorsque leffort axial crot et devient suprieur un deuxime effort critique
(note Fcr _ hel ), on parle de flambage hlicodal de la tige, ce qui correspond une forme hlicodale
prise par cette dernire. Si les diffrents auteurs arrivent la mme valeur de leffort critique de
flambement sinusodal, il existe des divergences vis vis de leffort critique de flambage hlicodal.
En effet, lquation de la force critique de flambage hlicodal dveloppe par Chen et al (1989) pour
des puits inclins est :

Fcr _ hel = * Fcr _ sin avec = 2 .

(2)

Wu (1993) propose de corriger cette quation en considrant une autre valeur de . Ce qui
donne :

Fcr _ Wu = (2 2 1) * Fcr _ sin .

(3)

Contrairement aux auteurs prcdents, Mitchell (1995) dans ses premiers travaux sur ltude
analytique du flambage hlicodal trouve que Fcr _ hel vaut plutt :

Fcr _ Mitc = (2 2 ) * Fcr _ sin .

(4)

Les diffrences entre les auteurs sont probablement dues aux diffrentes hypothses retenues dans leur
travail et notamment ltat de chargement ou de dchargement de la tige.
2.2. Description du modle
Lanalyse du flambage dune tige se fait tant sur le plan quantitatif (dtermination des forces
critiques) que qualitatif (sur le trajet dquilibre en post-flambage). Une dtermination analytique des
charges critiques de post-flambage nest possible que lorsque les structures tudies sont relativement
simples. Dans le cas contraire, lutilisation dun outil numrique devient indispensable.
Le code de calcul dcrit ici, appel ABIS, vient complter celui dcrit par Belaid et al (2005)
o lon peut retrouver les principales hypothses du modle et les tapes de rsolution. Ce modle fait
appel une rsolution numrique des quations globales et un algorithme de contact, permettant de
faire une analyse tridimensionnelle du comportement mcanique des garnitures de forage. Ce nouveau
code de calcul bas sur une approche semi-analytique permet aujourdhui de mener une analyse du
flambage de lensemble du train de tige de forage dans un puits tridimensionnel, et ceci sans avoir
recours la mthode des lments finis. Cette avance numrique est particulirement intressante
lorsque les longueurs des tiges atteignent plusieurs centaines de mtres, permettant un gain de temps
de calcul trs apprci.
Le principe de lanalyse de flambage dans ABIS consiste perturber une solution initiale
simple, puis de faire tourner un algorithme itratif jusqu stabilisation, et obtention de la solution
finale, quilibre du systme. La figure 2 est le rsultat dun calcul de flambage (avec ABIS) dune tige
de forage dans un parfait puits horizontal, avec la mise en exergue des vecteurs contacts 3D entre le
puits et la garniture.

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

La figure 3 montre lexemple du flambage dEuler. Il sagit dans ce cas de dterminer la

configuration prise par une poutre de longueur 20 m, encastre un bout et encastre glissant lautre
bout, non contrainte latralement et soumise une compression de 200 N. Comme on le remarque sur
cette figure, il a fallu seulement 10 itrations, pour obtenir la solution au problme pos.
3. Banc exprimental de flambage
3.1. Dispositif exprimental et matrice dessais
Le premier simulateur de flambage de tige a t initi par Sorenson (1984) pour une modlisation
du flambage des tiges conditionn par un cylindre circulaire. Par la suite, Wu (1993), puis Salis
(1994) se sont intresss la modlisation exprimentale du flambage des tiges dans des puits
verticaux, inclins et horizontaux. Ces travaux ont t tendus des puits courbes par Suryanarayana
(1995). Ces diverses tudes exprimentales ont t ralises pour des cas spcifiques et idaliss de
puits o la gomtrie est parfaitement droite. Les puits de forage tant naturellement tortueux, nous
pensons que la prise en compte de cet aspect dans notre analyse exprimentale est primordiale (fig. 4b). Cette tortuosit a t reprsente par des dog legs, terme utilis dans le milieu ptrolier pour dcrire
la courbure de la trajectoire du puits.
Le dispositif exprimental de flambage a t construit lchelle rduite via un facteur
dchelle (diamtre extrieure de la tige/diamtre intrieure du puits) qui permet de respecter la
similarit gomtrique.
Une tige en acier de diamtre extrieure 13.5 mm et dpaisseur 2 mm, de 16.75 m de long a
t introduite dans un tube transparent en PVC (de diamtre intrieur 42 mm) qui reprsente le puits.
Ce tube est maintenu par des manchons sur un banc parfaitement horizontal grce des supports
ajustables (fig. 4-a). De plus le tube transparent est assez flexible pour produire des dog legs le long du
chemin afin de reproduire la tortuosit dun puit de forage. Les extrmits de la tige sont fixes au bti
via des liaisons la fois de type rotule ou encastrement selon la nature de lessai effectuer et pour
connatre leffet des conditions aux limites. Deux capteurs de forces permettent la mesure des efforts
aux deux extrmits. La compression axiale est applique par un vrin hydraulique au niveau de lune
des extrmits. Un moteur asynchrone install cette extrmit permet de mettre la tige en rotation. A
cette extrmit un capteur de dplacement mesure le dplacement axial gnr. Les efforts (force et
couple) et le dplacement longitudinal de la tige sont enregistrs sur un ordinateur via un systme
dacquisition pour une analyse post exprimentale.
Ces essais sont effectus une vitesse de chargement axial contrle de 3.5 10 3 m/s pour
viter tout effet dynamique. Un essai type se droule comme suit. On commence par initialiser le
systme dacquisition et on dclenche le vrin hydraulique pour un chargement, durant lequel, on
enregistre les efforts et le dplacement longitudinal. Cette phase de chargement, est suivie par une
phase de dchargement.
Nous avons ralis les essais de flambage dans 4 types de trajectoire :
1- Trajectoire parfaitement horizontale : le tube est horizontal (inclinaison = 90, fig. 5a).
2- Trajectoire tortueuse avec variation dinclinaison : le tube reste dans le plan vertical, fig. 5b.
3- Trajectoire tortueuse avec variation dazimut : le tube reste dans le plan horizontal, fig. 5c.
4- Trajectoire tortueuse avec variation dinclinaison et dazimut : trajectoire 3D, fig. 5d.
Les tests sont rpts plusieurs fois pour vrifier la reproductibilit des essais. Par ailleurs la
dforme de la tige est mesure tous les 10 cm de tige en reprant sa position vis vis du tube.
4. Rsultats exprimentaux et prdiction du modle
Trajectoire parfaitement horizontale
La figure 6 montre la comparaison entre les rsultats dABIS et ceux dun essai ralis sur le
banc parfaitement horizontal. Lhistorique du chargement se dcompose en 4 phases repres par les
points O-A-B-C. Les observations suivantes ont t effectues.
En effet, le point A correspond au dclenchement du flambage sinusodal prdit par ABIS. Le
code de calcul permet galement de prdire lamplitude de cette sinusode en fonction de la

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

compression applique (Menand et al (2006)). La force critique Fcr _ sin partir de laquelle la tige se
met flamber sinusodalement est galement bien estime avec lquation (1). Ds lors que la tige a
flamb de faon sinusodale, une lgre augmentation de la compression entrane une grande
augmentation du dplacement, et la tige commence flamber de manire hlicodale. Cet tat
correspond au point B qui caractrise la formation de la premire hlice. Elle est obtenue
exprimentalement une compression de 1000 N, valeur assez proche de la force critique Fcr _ hel
prdite par lquation (2). Notons que la premire hlice thorique calcule par ABIS est obtenue pour
une force de compression gale 1250 N. Le trajet BC dtermine loccurrence des autres hlices.
Quant au point C, cest le point o lon observe la forme hlicodale complte pour une compression
de 2000 N avec la formation de 4/5 hlices. Cette valeur de la compression est proche de celle prdite
par ABIS (2100 N) pour avoir la formation complte de toutes les hlices. Au-del dune compression
de 2700 N, la tige est dcharge. Les hlices se dfont, la tige perd aussi sa forme sinusodale, et
retourne finalement son tat dorigine.
La figure 7-a montre une comparaison entre la configuration hlicodale de la tige mesure au
laboratoire et celle calcule avec ABIS (pour une force axiale fixe 2100 N). Comme on peut le voir,
except les deux extrmits de la tige, ABIS est capable de prdire correctement le pas de lhlice qui
est gale 2.90 m
Trajectoire tortueuse : variation dinclinaison
Le tableau 1 montre une synthse des rsultats exprimentaux obtenus aprs les tests de
variation dinclinaison (configuration de la fig. 5-b). On remarque que la force critique dobtention de
la premire hlice diminue mesure que les dog legs de la conduite augmente. Les simulations avec
ABIS des configurations testes au laboratoire confirment cette tendance. Il semblerait donc que la
tortuosit observe dans les puits de forage favorise lapparition du flambage hlicodal.
Trajectoire tortueuse : variation dazimut
Contrairement la trajectoire prcdente, les essais ont rvl que plus le gradient dazimut est
grand, plus la force critique de flambage hlicodal est grande (tableau 1, configuration de la fig. 5-c).
Cela montre que dans les puits horizontaux o lazimut uniquement varie, les tiges pourraient avoir
plus de rsistance flamber.
Trajectoire tortueuse et tridimensionnelle
La combinaison des deux configurations prcdentes donne lieu la configuration
tridimensionnelle de la figure 5d o on fait varier simultanment lazimut et linclinaison. Dans ce cas
de figure, cest leffet de la variation dazimut sur la force critique qui lemporte sur celui de
linclinaison. En effet la force critique de flambage est du mme ordre de grandeur que celle obtenu
pour le cas prcdent.
Les graphes de la figure 7 montrent la comparaison entre la dforme obtenue par ABIS et
celle mesure au laboratoire. On remarque clairement que le code de calcul permet de prdire
correctement la dforme de la tige dans un puits 3D lorsque celle-ci flambe, et ceci quelle que soit le
type de tajectoire.
5. Conclusion
Le dispositif exprimental prsent valide bien le nouveau modle ABIS pour des
configurations de puits complexe. Ce modle danalyse du flambage des tiges de forage a permis de
montrer limportance de la prise en compte des dog legs dans le calcul de flambage. Les effets de la
rotation et des connexions sur le flambage seront prochainement tudis pour se rapprocher au mieux
des conditions oprationnelles du milieu ptrolier.
References
[1] Dawson R. and Paslay P.R. (1984) Drillpipe Buckling in Inclined Holes, Journal of Petroleum
Technology 36(11), p 1734-1738.
[2] Paslay P.R. and Bogy D.B. (December 1964) The Stability of a Circular Rod Laterally
Constrained to Be in Contact With an Inclined Circular Cylinder, Journal of Applied Mechanics, p
605-610.

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

[3] Chen Y.C., Lin Y.H. and Cheatham J.B. (1989) An Analysis of Tubing and Casing Buckling in
Horizontal Wells, The 21st Annual OTC, Houston, TX.
[4] Wu J. and Juvkam-Wold H.C. and Lu R. (January-February 1993) Preventing Helical Buckling of
Pipes in Extended Reach and Horizontal Wells, Energy-Sources Technology Conferences &
Exhibition, Houston, TX.
[5] Mitchell R. (April 1995) Effects of Well Deviation on Helical Buckling, Production Operations
Symposium, Oklahoma City, OK.
[6] Belaid A., Sellami H., Menand S. et Tijani M. (Septembre 2005) Modlisation du comportement
mcanique dune structure lance rigide lintrieur dun trou trajectoire complexe : application
au dimensionnement du forage ptrolier, 17me Congrs Franais de Mcanique, Troyes.
[7] Sorenson K. (1984) Post-buckling behavior of a circular rod constrained within a circular
cylinder, Ph.D. thesis, Rice University.
[8] Wu J. and Juvkam-Wold H.C. (1993) Helical Buckling of Pipes in Extended Reach and Horizontal
Wells - part 2 : Frictional drag analysis, Transaction of the ASME, 115:196-201.
[9] Salis J.B. (1994) Experimental Study and Mathematical Modeling of helical Buckling of Tubulars
in Inclined Wellbores, Ph.D. Dissertation, University of Tulsa.
[10] Suryanarayana P.V.R. and McCann R.C. (1995) An experimental Study of Buckling and PostBuckling of Laterally Constrained Rods, Journal of Energy Resources Technology, vol. 117.
[11] Menand S., Sellami H., Tijani M. and Akowanou J. (September 2006) Buckling of tubulars in
Actual Field Conditions, papier spe 102850 presented at the Annual Technical Conference and
Exhibition, San Antonio, TX.

Gradient (/30 m)\ Force (N)

0
1
1.5
2
2.5

Pour linclinaison
1000
971
823
721
700

Pour lazimut
1000
1355
1760
2268
2634

Tableau 1 : Variation de la force critique de flambage hlicodale (exprimentale) en fonction du

gradient dinclinaison et dazimut.

Configuration droite

Tige de forage

Tension
(pas de flambage)

(a)

Point neutre

Compression
(risque de flambage)

Puits

Configuration sinusodale

(b)

Configuration sinusodale

Configuration hlicodale

Poids sur loutil

Vue de dessus

Section droite du puits

Fig. 1: (a) Profil de la force axiale dans un train de tige de forage (b) Configuration
successive dune tige comprime
5

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

Fig. 2 : Simulation ABIS : Visualisation

3D du flambage dans un puits horizontal
les vecteurs contacts sont reprsents par
les flches rouges.

Vrin
hydraulique

Manchons de fixation
paroi rigide bti

Fig. 3 : Flambage dEuler dune poutre de 20

mtres de long soumis la compression de
200 N.

capteur d'effort
bti

acquisition
capteur de
dplacement LVDT

Moteur

Pieds
ajustable

tige

Paroi
rigide
transparent

(a)

(b)

Fig. 4: Schma simplifi du simulateur ltat horizontal et avec des dog legs
a) horizontale

b) inclinaison

c) azimut

d) 3D

Fig. 5: Exemple de configurations de trajectoires adoptes lors des essais

6

18me Congrs Franais de Mcanique

Grenoble, 27-31 aot 2007

Exprimental

Numrique:ABIS

Compression (N)

2500

2000
Fcr_Mitc

Eq. (4)

Fcr_Wu

Eq. (3)

Fcr_hel

Eq. (2)

Fcr_sin

Eq. (1)

1500

1000

500

Dplacement (mm)

0
0

10

Fig. 6 : Courbes thorique et exprimentale de la compression en fonction du dplacement

longitudinal dans le cas du banc horizontal

Axe du tube
Numrique ABIS
Exprimental

Puits parfaitement horizontal

Puits tortueux verticalement

Axe du tube
Numrique ABIS
Exprimental

Puits tortueux horizontalement

Axe du tube
Numrique ABIS
Exprimental

Puits tortueux 3D

Fig. 7 : Essais de flambage dans le cas de trajectoires horizontale et tortueuses

7