Brevet FR2915298A1 — publié le 20081024

La présente invention concerne le domaine de l'extraction et du transport d'effluents pétroliers. En particulier, l'invention concerne une méthode pour dimensionner des installations industrielles où des effluents pétroliers, comportant une phase liquide et une phase gazeuse, 15 s'écoulent en régime intermittent. Présentation de l'art antérieur Dans l'industrie pétrolière, il est très important de dimensionner correctement les installations industrielles d'extraction et de transport d'effluents pétroliers. Pour ce faire, il est 20 nécessaire de prévoir les pertes de charge dues aux écoulements polyphasiques dans ce type de conduites. Généralement, l'écoulement d'effluents pétroliers est un écoulement comprenant une phase liquide et une phase gazeuse. Cet écoulement diphasique peut alors présenter différents régimes : 25 écoulement stratifié ( stratified flow ) : c'est un écoulement à phases séparées. écoulement dispersé ( dispersed flow ) : il s'agit par exemple d'un écoulement de liquide dans lequel se forme des bulles. 1 - écoulement annulaire ( annular flow ) : c'est un écoulement à phases séparées dans lequel le liquide mouille complètement la paroi en formant un film annulaire autour de l'écoulement de gaz. - écoulement intermittent ou à poches et bouchons ( slug flow ). L'invention concerne plus particulièrement ce dernier type d'écoulement. Un écoulement intermittent est observé pour des débits de gaz et de liquide "moyens". Sa structure présente une succession de poches de gaz, appelées poches de Taylor, et de bouchons de liquide qui peuvent contenir de petites bulles de gaz. C'est une configuration mixte entre un écoulement stratifié et un écoulement dispersé. Dans cette configuration il peut se produire un phénomène d'entraînement de gaz, le gaz étant entraîné depuis la. poche de gaz vers le bouchon liquide sous la forme de bulles millimétriques. La fraction de gaz dans le bouchon liquide dépend alors du flux de gaz entraîné en arrière de la poche. L'arrière de la poche est défini relativement au sens d'écoulement de l'effluent (écoulement de l'arrière vers l'avant). L'arrière correspond ainsi à l'amont (l'aval étant l'avant). Pour parvenir à estimer les pertes de charge de ce type d'écoulement (diphasique intermittent), il est nécessaire de calculer le taux de vide moyen, qui fait lui même intervenir le taux de vide dans les poches de gaz et le taux de vide moyen dans les bouchons de liquide. On appelle taux de vide d'un mélange gaz / liquide le rapport entre le volume occupé par le gaz et le volume total du mélange. Pour calculer le taux de vide dans le bouchon liquide il est connu d'utiliser une corrélation empirique. Parmi les publications récentes utilisant cette approche, on peut citer : - Felizola, H., Shoham, O. 1995. A unified model for slug flow in upward inclined pipes. Journal of Energy Ressources and Technology, 117, pp 7-12. - Gomez, L.E., Shoham, O., Taitel, Y. 2000. Prediction of slug liquid holdup: horizontal to upward vertical flow. International Journal of Multiphase Flow, 26, pp 517-521. On connaît également une méthode basée sur des modèles physiques pour calculer le taux de vide dans le bouchon de liquide et la poche de gaz. Cette méthode emploie un modèle physique pour le flux de gaz entraîné en arrière de la poche de gaz (dans un repère se déplaçant à la vitesse de la poche). Dans cette méthode, le modèle de flux de gaz entraîné est basé sur l'hypothèse qu'un jet turbulent de liquide, présent en arrière de la poche, est responsable de l'entraînement de gaz. Cependant, comme le montre le document suivant, les résultats associés à différentes propriétés des fluides en présence montrent qu'il peut conduire à des incohérences et ne donne pas une description fine du taux de vide. -Brauner, N, Ullmann, A. 2004. Modelling of gas entrainment from Taylor bubbles. Part A: slug flow. International Journal of Multiphase flow, 30, pp 239-272. En somme, aucune méthode ne permet d'estimer les pertes de charge d'un d'écoulement diphasique intermittent de façon précise, et en particulier aucune méthode ne tient compte des effets d'entraînement de gaz en arrière de la poche de gaz. Ces effets sont importants car l'aération du bouchon liquide influence de manière significative les pertes de charge. Or, la détermination des pertes de charge est indispensable pour déterminer les dimensions d'installations telles que des conduites d'effluents pétroliers. Ainsi, l'objet de l'invention est une méthode pour dimensionner des installations industrielles où des effluents pétroliers, comportant une phase liquide et une phase gazeuse, s'écoulent en régime intermittent. Au cours de cette méthode, on estime des pertes de charges au sein des installations (conduite d'extraction, conduite de transport d'effluents ...) en tenant compte des effets d'arrachement de gaz en arrière de la poche de gaz. La méthode selon l'invention L'invention concerne une méthode pour dimensionner des installations industrielles où un mélange diphasique comportant une phase liquide et une phase gazeuse s'écoule selon une configuration comportant une succession de bouchons de liquide et de poches de gaz en arrière desquelles du gaz est entraîné, dans laquelle on modélise l'écoulement de chaque phase au sein d'une poche de gaz et l'écoulement de chaque phase au sein d'un bouchon de liquide à l'aide d'un premier modèle physique, et on modélise le flux de gaz entraîné à l'aide d'un second modèle physique. La méthode comporte les étapes suivantes : -on définit le second modèle physique dans lequel ledit flux de gaz entraîné est proportionnel à une variation de pression entre une poche de gaz et un bouchon de liquide en arrière de cette poche, et dans lequel on prend en compte une condition critique de formation dudit flux de gaz entraîné, définie par une variation de pression telle que des forces de pression engendrées par cette variation de pression soient supérieures à des forces de tension superficielle entre le gaz et le liquide ; on initialise et calibre ledit second modèle à l'aide de mesures expérimentales ; - on détermine un gradient de pression dans la poche, un gradient de pression dans le 5 bouchon et le rapport d'une longueur de poche de gaz sur une longueur de bouchon de liquide, à l'aide d'une méthode itérative dans laquelle on ajuste un flux de gaz entraîné y'G,ent calculé à l'aide dudit second modèle, avec un flux de gaz moyen ytG obtenu par une équation de conservation des flux en écoulement intermittent établi à partir dudit premier modèle d'écoulement physique, 10 - on détermine des pertes de charges au sein des installations industrielles à l'aide desdits gradients de pression et dudit rapport, et - on détermine les dimensions des installations industrielles de façon à minimiser lesdites pertes de charges. 15 Selon l'invention, on peut déterminer le flux de gaz entraîné y'G,ent en considérant qu'une proportion KAp du travail effectué par les forces de pression est utilisée pour l'entraînement de gaz. On peut également, pour déterminer le flux de gaz entraîné Y' G,ent prendre en compte une valeur maximale du flux de gaz entraîné via une hypothèse d'écoulement homogène et sans glissement. 20 Le premier modèle physique peut comporter soit un modèle d'écoulement stratifié basé sur l'égalité des gradients de pression dans les deux phases, soit un modèle d'écoulement annulaire basé sur l'égalité des gradients de pression dans les deux phases. Il peut en plus comporter un modèle de type flux de dérive. Un critère d'arrêt pour la méthode itérative peut être défini par le critère de convergence 25 suivant : ~G,ent ~V G -V G La méthode comporte une modélisation physique de ce type d'écoulement, puis une estimation des pertes de charge au sein de la conduite. Enfin, on déduit le dimensionnement adéquat des installations à partir des pertes de charge. La méthode est basée sur l'égalité entre le flux de gaz entraîné en arrière de la poche de gaz et le flux de gaz obtenu par l'équation de conservation des flux en écoulement intermittent établi. La poche de gaz de l'écoulement est modélisée à l'aide d'un modèle d'écoulement stratifié et la zone de bouchon liquide est modélisée par une approche des flux de dérive. Le flux de gaz entraîné en arrière de poche est formulé à l'aide d'une approche physique faisant intervenir les vitesses et taux de vide dans les différentes zones de l'écoulement. L'égalité entre les flux de gaz permet de fermer le problème et d'obtenir le taux de vide moyen. Dans un premier temps, les modèles physiques permettant de modéliser les écoulements sont présentés. Puis la méthode itérative permettant de déterminer des pertes de charges est décrite. Modélisation physique d'écoulements intermittents gaz - liquide L'invention concerne les écoulements diphasique gaz û liquide de type intermittent. On rappelle qu'un écoulement intennittent est séparé en deux zones : une poche de gaz (P) et un bouchon de liquide (B). Pour modéliser de tels écoulements, on les caractérise par les 20 paramètres suivants : Concernant la poche de gaz : la vitesse de la poche de gaz VP le taux de vide moyenné sur la section poche RGP le taux de liquide dans la poche RLP 25 la vitesse moyecme du gaz dans la poche VGP la vitesse moyenne du liquide dans la poche VLP Concernant le bouchon de liquide : le taux de vide moyenné sur la section bouchon RGB15 le taux de liquide dans le bouchon RLB la vitesse moyenne du gaz dans le bouchon VGB la vitesse moyenne du liquide dans le bouchon VLB Ainsi que : le taux de vide moyenné sur la section totale RG le taux de liquide moyenné sur la section totale RL Selon l'invention, on réalise une modélisation physique des écoulements intermittents gaz -liquide, en définissant trois modèles physiques : un premier modèle physique (MP) décrivant l'écoulement de chaque phase au sein d'une poche de gaz (P), un second modèle physique (MB) décrivant l'écoulement de chaque phase au sein d'un bouchon de liquide (B), un troisième modèle physique (MF) décrivant le flux de gaz entrainé en arrière de la poche. Modèle physique (MP) décrivant l'écoulement dans la poche de gaz (P) On modélise l'écoulement dans la poche de gaz (P) à l'aide d'un modèle physique faisant intervenir le taux de vide dans la poche (RGP ), le taux de liquide dans la poche (RLP ), la vitesse moyenne du gaz dans la poche (VGP) et la vitesse moyenne du liquide dans la poche (yLP). Selon l'invention, on considère que l'écoulement de gaz dans la poche de gaz est de type stratifié en conduite inclinée ou annulaire en conduite verticale. Le modèle est alors basé sur l'égalité des gradients de pression dans les deux phases. En combinant les équations de conservation de la quantité de mouvement pour le gaz et pour le liquide afin d'éliminer le gradient de pression, il vient : S,e [RLPZGPxGP ù RGPV LP%LP + VGLPxGLP ù APRGPRLPg sin B = 0 . (8) Avec : pe : le périmètre de la conduite S : la section de la conduite xLp : le périmètre mouillé par le liquide (adimensionnel) xGP : le périmètre mouillé par le gaz (adimensionnel) xGLP : le périmètre interfacial gaz ù liquide (adimensionnel) rLp : la contrainte pariétale due au liquide. Elle est fonction de VLP ,RLP , xLP zcp : la contrainte pariétale due au gaz . Elle est fonction de Vcp , RLP , xcp VGLp la contrainte interfaciale gaz - liquide dans la zone d'écoulement stratifié. Elle est fonction de VG p , RLP , xcLP dp : différence des densités [kg/m3] B : angle de la conduite avec l'horizontale [rad] De plus, la vitesse de la poche de gaz, Vp , est calculée par la relation suivante : Vp =CpUm +Vdri,p, avec : D : diamètre de 1.a conduite U,,, : vitesse superficielle du mélange Vdro p vitesse relative de la poche par rapport au mélange et où le coefficient Cp et la vitesse Vdrtfr, p sont fonction de Um , e, D , et sont donnés par des modèles de la littérature, tels que ceux décrits dans les documents suivants : Fabre J., Liné A. 1996. Slug flow modelling. International encyclopaedia of heat and mass transfer. Innodata corp., 1015-1021. (7) [m]25 Modèle physique (MB) décrivant l'écoulement dans le bouchon de liquide (B) On modélise l'écoulement dans le bouchon de liquide (B) à l'aide d'un modèle physique faisant intervenir le taux de 'vide dans le bouchon (RGB) et la vitesse moyenne du gaz dans le bouchon ( VGB ). Dans le bouchon liquide, la vitesse de la dispersion de bulles, c'est-à-dire la vitesse moyenne du gaz dans le bouchon liquide (B), peut être modélisée par une approche des flux de dérive ( drift-flux ). Cette approche est décrite par exemple dans le document suivant : Zuber, N., Findlay, J.A. 1965. Average volumetric concentration in two-phase flow systems. J. Fleat Transfer Trans. ASME Ser., 87, pp 453-468. Suivant ce modèle, on écrit : VGB = COB Um + Vdrifi , (6) où CoB , fonction de RGB et B, et Vdr ft , fonction de 6', sont donnés par des modèles de la littérature. On peut se référer par exemple à : Guet, S., Ooms, G., Oliemans, R.V.A., Mudde, R.F. 2004. Bubble size effect on low liquid input drift flux parameters. Chemical Engineering Science, 59, pp 3315-3329. Modèle physique décrivant le flux de gaz entraîné en arrière de la poche Le flux de gaz entraîné en arrière de la poche de gaz (P) est modélisé à l'aide d'un modèle physique. On rappelle que l'arrière de la poche est défini relativement au sens d'écoulement (écoulement de l'arrière vers l'avant). 20 La figure 2 illustre le principe du modèle physique du flux de gaz entraîné vG ent selon l'invention : l'écoulement dans la conduite C, inclinée d'un angle 8 avec l'horizontale, comporte deux sections : une poche de gaz (P) avançant à la vitesse Vp, et un bouchon de liquide (B). On fait l'hypothèse que l'écoulement est développé et à l'état d'équilibre. Les flux de gaz 25 et de liquide sont alors équivalents à chaque frontière de l'écoulement. En écoulement établi, ce flux de gaz entraîné Y'G,en, est équivalent aux flux de gaz calculés par bilan sur chaque section (P et B) de l'écoulement (VGP et yiGB ), et au flux de gaz15 moyen VG . Le flux de gaz moyen est le flux de gaz obtenu par l'équation de conservation des flux en écoulement intermittent établi : yrG = RG (VP ù VG), VG étant la vitesse moyenne du gaz. YjG ùY" GP ù V1GB = V G,ent Les flux de gaz sont donnés par : Y~G = RG (VP ù VG ) V'f GB = RGB (v ùVGB) Y1GP ù RGP(VP VGP) Ces flux sont donnés dans un référentiel se déplaçant à la vitesse de poche VP . Le modèle selon l'invention emploie le même type d'égalité sur les flux pour la phase liquide : Y~L = 'LP =VLB (3) Les flux de liquide sont donnés par : (4) (Y~L =RL(VP -VL) =RLB(VùVLB) iY1LB =RLP(VP ùVLP) (5) Les taux de vide et de liquide sont par ailleurs liés par les égalités : 'RG +RL =1 RGP+RLP=1 RGB +RLB =1 Dans certaines conditions, du gaz est entraîné depuis la poche de gaz vers le bouchon liquide sous la forme de bulles millimétriques. Ceci conduit à une valeur non nulle du taux de vide dans le bouchon liquide ( RGB ≠ 0). Dans un repère mobile se déplaçant à la vitesse de la 15 poche (Vp), le flux de gaz entraîné est noté VG ,e . Les vitesses moyennes des phases étant très différentes dans la zone de poche de gaz et dans la zone de bouchon liquide, les effets d'entraînement de gaz peuvent être attribués à deux phénomènes : l'agitation turbulente dans le liquide dans la zone arrière de la poche, en raison du jet turbulent de liquide (Brauner et Ullmann, 2004), 20 le saut de pression hydrodynamique, nécessaire pour accélérer le liquide depuis la poche de gaz jusqu'au bouchon de liquide. (1) (2)10 Selon l'invention le flux VG,, de gaz entraîné en arrière de la poche de gaz (P) est modélisé à l'aide d'un modèle physique (MF). Dans ce modèle, on considère que ce flux est proportionnel à une source d'énergie et que le saut de pression en arrière de poche est la source d'énergie responsable de l'entraînement de gaz en arrière de poche. Le saut de pression nécessaire pour accélérer le liquide entre la poche de gaz et le bouchon est donné par : dP = PLy1L (VLB ù VLP) L'aération du bouchon liquide est gouvernée par une compétition entre cette force de pression dl' et la force due à la tension superficielle, z = cr : tension superficielle gaz - liquide [N/m] : taille maximale des bulles Par ailleurs à faibles vitesses débitantes il est fréquent que le bouchon liquide ne soit pas aéré, donc le flux de gaz entraîné est nul. La limite d'existence de l'entraînement est décrite par le cr d Mai dm. rapport de ces deux forces, 6 dPMa, Ic__dD rD c J L J, On définit une valeur critique AP, telle que seul l'excès d'énergie présent par rapport à cette valeur critique, c'est à dire dP > dP~, conduit à de l'entraînement de gaz. En considérant qu'une proportion KAp du travail effectué par les forces de pression est utilisée pour l'entraînement de gaz et en incluant de plus les conditions critiques d'entraînement de 20 gaz en arrière de poche, le flux de gaz entraîné s'exprime par : iUC.e = K6 d m (dP ù dP~) (9) Le paramètre K1 ainsi que la valeur critique du saut de pression dP~ nécessaire pour accélérer le liquide entre la poche de gaz et le bouchon sont détettuinés de manière spécifique. Détermination de la constante Al', 25 Le saut de pression critique Cc est spécifique à chaque calcul de flux de gaz entraîné. Il est déterminé par le modèle. On applique la considération suivante : s'il n'y a pas de gaz, la vitesse superficielle du gaz dans la conduite est nulle (UG =0) et le bouchon de liquide n'est pas aéré. Donc lorsque UG =0 le flux de gaz entraîné doit vérifier y'Ge = yrG =0 et dP = dPc . Comme UG =0 dans ces conditions la vitesse du mélange est égale à la vitesse superficielle du liquide (UL = Um ). Pour chaque condition d'écoulement considérée par le modèle une valeur critique dP, est en un premier temps calculée. Pour ce faire, on applique les modèles MB et MP en imposant yrce ='VG = 0 pour UL = U.. Les valeurs des vitesses du liquide dans le bouchon et dans la poche sont alors obtenues (VLB,c et VLP c ). Le flux de liquide critique associé est donné par y'L c = VP ù UL . Puis le saut de pression associé est calculé : ~c = PLVVL,c `VLB,c ù VLP,c) • On obtient ainsi une valeur de saut de pression dP correspondant à la valeur du saut de pression critique dPc. La valeur de la constante dl), est donc spécifique à chaque calcul, et doit être calculée en un premier temps. Détermination de la constante KaP La valeur de K11 est obtenue à partir de données expérimentales représentatives d'un écoulement d'un fluide de propriétés données dans une conduite de géométrie connue. En utilisant l'égalité des flux de gaz y'G = vG,e 66V'G Ka, = dMc (dPùdPe)ùf(VVG,VL,,dMaz,dPc,VLP,VLB) La constante KAp est ainsi déterminée en utilisant des mesures expérimentales de flux de gaz et de liquide, V1G et ytL . Pour ce faire une campagne de mesures expérimentales doit être conduite. Des données expérimentales sont collectées en conditions d'écoulement intermittent concernant : la vitesse de poche VP et le taux de vide R G . Les flux de gaz et de liquide sont alors déterminés expérimentalement : VrG = RGVP ù UG et V'L = RLVP ù UL . Les modèles de calcul MP et MB sont appliqués pour déterminer les valeurs de dPc, VLP et VLB associées aux points expérimentaux considérés. La valeur de la constante Kp ainsi obtenue est indépendante des conditions de vitesse débitante, et dépend essentiellement du diamètre de la conduite et des propriétés des fluides considérés (tension superficielle et niveau de contamination du liquide). Une valeur obtenue à l'aide de données expérimentales représentatives des conditions industrielles considérées (concernant le diamètre et :les propriétés des fluides) peut de ce fait être appliquée à des calculs de flux de gaz entraîné en conditions d'écoulement industriel. Enfin, selon l'invention, on considère que le flux de gaz entraîné ne peut pas dépasser une valeur maximale, notée vG Mnx . En effet, la poche de gaz présente une vitesse ascendante positive et supérieure à la vitesse moyenne du gaz dans le bouchon. Le taux de vide moyen RG est donc toujours inférieur au taux de vide obtenu avec une hypothèse d'écoulement homogène sans glissement RG,nos : RG R G VG = UG '' 0 . (10) Donc le flux de gaz doit vérifier : (11) ~Vc = RG(Vp ùVG) Estimation des pertes de charge au sein de la conduite La figure 1 illustre les étapes de la méthode. Les onze inconnues sont : (v ,RGB,RGP,RLB,RLP,VGB,VGP,VLB,VLP,RG,RL)• Ces onze inconnues sont déterminées avec les onze relations indépendantes qui constituent le modèle physique d'écoulement. Une méthode itérative est employée pour résoudre le problème. A cet effet une méthode de relaxation du flux de gaz moyen yJG est appliquée à chaque itération. Initialisation Il est nécessaire de définir une valeur initiale du flux de gaz moyen yrG . Cette valeur initiale n'a pas d'influence sur le résultat final et peut donc être choisie totalement arbitrairement. Boucle d'itération (ITE) La boucle d'itération comporte trois étapes successives : Tout d'abord on détermine le taux de vide moyenné sur la section poche (RGP ), le taux de liquide dans la poche (RL ), et la vitesse moyenne du gaz dans la poche (VGP) à l'aide d'un modèle physique décrivant la physique de l'écoulement de chaque phase au sein de la poche de gaz. Puis on détermine le 1:aux de vide moyenné sur la section bouchon (RGB ), le taux de liquide dans le bouchon (RLB ), la vitesse moyenne du gaz dans le bouchon (VGB ), et la vitesse moyenne du liquide dans le bouchon (VLB) à l'aide d'un modèle décrivant la physique de l'écoulement de chaque phase au sein des bouchons de liquide. Enfin, à l'aide des calculs précédents, on détermine le flux de gaz entraîné y'G ent à partir du troisième modèle. Ensuite, ce flux de gaz entraîné Y'G,ent est comparé au flux de gaz moyen y/G donné en entrée de la boucle. On calcule un nouveau flux de gaz moyen y/G à partir du flux de gaz entraîné yya ent 20 On peut par exemple prendre la moyenne des deux : Y~G + VI G,ent 2 Ce nouveau flux de gaz moyen y'G est donné en entrée de la boucle, et un nouveau calcul du flux de gaz entraîné yrG ent est réalisé. On rappelle que le flux de gaz moyen est le flux de gaz obtenu par l'équation de conservation des flux en écoulement intermittent établi : 25 YG=RG(vvùVG) Ainsi, en fournissant un nouveau G , on calcule un nouveau taux de vide moyenné sur la section totale, RG . En imposant Y1GB = Y1G et Y'GP =y/G, les modèles MB et MP permettent ainsi d'obtenir de nouvelles valeurs pour les vitesses et fraction de gaz et de liquide dans la poche et dans le bouchon. Les itérations sont arrêtées lorsque le flux de gaz entraîné vérifie un critère de convergence. On peut par exemple choisir le critère suivant : YYG,ent ù YIG G Ces équipements visent à amortir les fluctuations du débit de liquide générées par l'écoulement intermittent. Leur dimensionnement nécessite la connaissance des fractions de gaz et de liquide dans les différentes sections de l'écoulement ainsi que la longueur de la section poche et bouchon. Ainsi l'invention trouve une application industrielle dans l'exploitation de gisement pétrolier, que ce soit pour dimensionner les conduites de production et de transport d'hydrocarbures, ou pour la simulation du comportement hydrodynamique des conduites de production et de transport de, fluides pétroliers. REVENDICATIONS 1) Méthode pour dimensionner des installations industrielles où un mélange diphasique comportant une phase liquide et une phase gazeuse s'écoule selon une configuration comportant une succession de bouchons de liquide et de poches de gaz en arrière desquelles du gaz est entraîné, dans laquelle on modélise l'écoulement de chaque phase au sein d'une poche de gaz et l'écoulement de chaque phase au sein d'un bouchon de liquide à l'aide d'un premier modèle physique, et on modélise le flux de gaz entraîné à l'aide d'un second modèle physique, caractérisé en ce qu'elle comporte les étapes suivantes : - on définit le second modèle physique dans lequel ledit flux de gaz entraîné est proportionnel à une variation de pression entre une poche de gaz et un bouchon de liquide en arrière de cette poche, et dans lequel on prend en compte une condition critique de formation dudit flux de gaz entraîné, définie par une variation de pression telle que des forces de pression engendrées par cette variation de pression soient supérieures à des forces de tension superficielle entre le gaz et le liquide ; - on initialise et calibre ledit second modèle à l'aide de mesures expérimentales ; - on détermine un gradient de pression dans la poche, un gradient de pression dans le bouchon et le rapport d'une longueur de poche de gaz sur une longueur de bouchon de liquide, à l'aide d'une méthode itérative dans laquelle on ajuste un flux de gaz entraîné yrc,en,, 20 calculé à l'aide dudit second modèle, avec un flux de gaz moyen yrG obtenu par une équation de conservation des flux en écoulement intermittent établi à partir dudit premier modèle d'écoulement physique, -on détermine des pertes de charges au sein des installations industrielles à l'aide desdits gradients de pression et dudit rapport, et 25 - on détermine les dimensions des installations industrielles de façon à minimiser lesdites pertes de charges. 2) Méthode selon la revendication 1, dans laquelle on détermine le flux de gaz entraîné wG,ent en considérant qu'une proportion Kp du travail effectué par les forces de pression est utilisée pour l'entraînement de gaz. 3) Méthode selon l'une des revendications précédentes, dans laquelle on détermine le 5 flux de gaz entraîné vG en, en prenant en compte une valeur maximale du flux de gaz entraîné via une hypothèse d'écoulement homogène et sans glissement. 4) Méthode selon l'une des revendications 1 à 3, dans laquelle ledit premier modèle physique comporte un modèle d'écoulement stratifié basé sur l'égalité des gradients de pression dans les deux phases. 10 5) Méthode selon l'une des revendications 1 à 3, dans laquelle ledit premier modèle physique comporte un modèle d'écoulement annulaire basé sur l'égalité des gradients de pression dans les deux phases. 6) Méthode selon l'une des revendications précédentes, dans laquelle ledit premier modèle physique comporte un modèle de type flux de dérive. 15 7) Méthode selon l'une des revendications précédentes, dans laquelle la méthode itérative comporte un critère d'arrêt défini par le critère de convergence suivant : Y G,ent Y1G 'VG 8) Méthode selon l'urine des revendications précédentes, dans laquelle on détermine les dimensions des installations industrielles en déterminant des pertes de charges pour 20 différentes valeurs de propriétés géométriques données d'installations, et l'on sélectionne les installations ayant les propriétés géométriques minimisant les pertes de charges. 9) Méthode selon l'une des revendications précédentes, dans laquelle les installations industrielles sont des conduites d'effluents pétroliers. 10) Méthode selon l'une des revendications précédentes, dans laquelle les installations 25 industrielles sont des équipements de séparation pétrolière de type "Slug Catcher", et dans laquelle on détermine les dimensions de ces équipements de séparation en déterminant des fractions de gaz et de liquide et des longueurs relatives d'une poche et d'un bouchon liquide.