Brevet FR2821675A1 — publié le 20020906

dans une raffinerie de pétrole. La présente invention concerne une méthode pour détecter et contrôler la formation d'hydrates en tout point d'une conduite o circulent des fluides pétroliers polyphasiques. Etat de la technique o Produire des hydrocarbures depuis des gisements en eaux profondes (1500 3000 m) soulève un grand nombre de difficultés, essentiellement à cause des pressions élevoes et des températures basses qui y règnent. Ces difficultés se rencontrent dans tous les domaines du savoir-faire pétrolier: le forage et les interventions sur puits, les installations de traitement (FPSO, FWHP, etc.) et leurs ancrages, les technologies sous- marines (pompage polyphasique, séparation sous marine des phases des fluides pétroliers), les risers et ombilicaux ou le "flow assurance ". Dans le domaine de la production proprement dit, l'exploitant est tenu de garantir le bon écoulement malgré de nombreux obstacles comme la formation de bouchons de fluide ou "slugging", les dépôts minéraux et organiques ou bien encore la formation d'hydrates. En offshore profond les pressions plus élevées et les températures plus basses favorisent la formation d'hydrates. Les interruptions de production durant lesquelles les fluides refroidissent accroissent encore les risques de formation d'hydrates. Une détection tardive de la présence d'hydrates notamment par méconnaissance des conditions de leur formation dans les conduites, peut entraîner des problèmes de production coûteux: blocage complet des conduites de production et/ou de transport. Il est donc important pour l'exploitant, de disposer de moyens s d'évaluation des risques de toute nature, de façon à mettre en place des techniques de régulation et d' intervention et ainsi de pouvoir garantir un bon écoulement des fluides en circulation. Pour protoger et isoler thermiquement les conduites sous-marines, on les enferme dans un tube extérieur résistant à la pression hydrostatique (" souvent 0 désigné par " bundles " par les gens de l'art). Le plus souvent, on regroupe plusieurs conduites pour former des faisceaux. Ces conduites regroupées ont souvent des fonctions différentes. Elles servent à la production d' hydrocarbures, à l ' injection de fluides dans le gisement ou réservoir, à l'injection de gaz dans les conduites (" gas lift "), à la circulation de fluide de réchanffage, etc. Dans l'espace entre elles, on interpose par exemple un calorifugeage à faible conductivité thermique laissé à la pression atmosphérique ou mis sous vide, avec des cloisonnements placés à intervalles réguliers pour des raisons de sécurité. Avec un outil efficace d'évaluation en continu des risques de formation d' hydrates et autres dépôts, l ' exploitant peut intervenir soit par le biais d' un réchauffage des conduites soit par l'injection d'additifs. Cependant un tel outil est difficile à mettre en place et coûteux en temps de calcul surtout si on doit tenir compte de la composition détaillée des fluides pétroliers. Pour étudier plus facilement leur comportement, il est connu de les décrire comme une combinaison d'un nombre de composants ou pseudo composants beaucoup plus réduit que le nombre réel de constituants. On parle alors d'une composition regroupée par opposition à une composition détaillée. Par le brevet FR 2 753 535 et la demande de brevet FR 0003867 du demandeur, on connaît des méthodes pour prédire la température d' apparition de dépôts tels que des cires ou paraffines dans des bruts pétroliers comportant le regroupement ou " lumping " de tous leurs constituants en un nombre plus réduit de pseudo composants représentant chacun plusieurs classes d'hydrocarbures et dont les paramètres physicochimiques sont déterminés par combinaison des paramètres correspondants d'un certain nombre d'hydrocarDures purs regroupés dans une base s de données. Les pseuJo composants de cette formulation regroupée, sont appliqués à un module thermodynamique permettant de déterminer différents paramètres indicatifs des conditions de formation des dépôts. Différents logiciels de simulation existent sur le marché permettant de modéliser le comportement de fluides polyphasiques en circulation dans des 0 pipelines. Par les brevets ou demandes de brevet suivants: US 5 550 761, FR 2.756. 044 (US 6 028 992) et FR 2 756 045 (US 5 960 187), FR 00/08 200 et FR 00/09 889 du demandeur, et aussi par les publications suivantes: - Faille I. et Heintzé E., " A rough finite volume scheme for modeling two-phase flow in a pipeline", Computers & Fluids 28 (1999), et Pauchon C. et al, "TACITE: a comprehensive mechanistic model for two-phase flow", 6th BHRG Multiphase International Conference, Cannes, France, June (1993), on connaît notamment le code de simulation TACITE (marque déposée) qui simule le comportement transitoire des fluides polyphasiques en circulation en prenant en compte à tout instant leur composition. Avec un tel code prévisionnel, l'ingénieur de production est en mesure de définir les caractéristiques des réseaux d'acheminement de fluides, notamment d'hydrocarbures: conduites, séparateurs, vannes, systèmes de contrôle, etc. Un des objectifs importants dans ces opérations de simulation, c'est de 2s prédire avec précision les caractéristiques des écoulements transitoires lors de situations telles que: variation du débit d'entrée, dépressurisation à la sortie, arrêts de production et redémarrage, nettoyage des conduites (" pigging "), apparition et propagation de bouchons liquides, formation de dépôts susceptibles de se former dans certaines conditions thermodynamiques, etc. Les modes d'écoulement de fluides polyphasiques dans des tubes sont extrêmement variés et complexes. Les écoulements diphasiques, par exemple, peuvent être stratifiés, la phase liquide s'écoulant dans la partie inférieure de la conduite, ou intermittents avec une succession de bouchons liquides et gazeux, ou bien encore dispersés, le liquide étant entrâîné sous forme de fines gouttelettes. Le mode d'écoulement et le glissement entre les phases varient notamment avec l'inclinaison des conduites par rapport à l'horizontale et dépendent du débit, de la lO température etc. Le code TACITE, on le rappelle, est adapté à suivre les variations de la composition des mélanges d'hydrocarDures. Comme déjà décrit dans les brevets précités du demandeur, il y a une équation de conservation de masse pour chaque pseuJo-composant. Il y a aussi une équation de quantité de mouvement et une équation pour l'énergie du mélange. Les équations correspondantes sont rappelées ci- après, après une définition des symboles et notations: À V vitesse absolue de phase; À U vitesse superficielle; À P fraction volumétrique par phase; À p densitéde phase; À H enthalpie de phase; À E énergie interne; À P pression du mélange; À Tw frottement pariétal; À Qw flux de chaleur pariétal entre le pipeline et le milieu environnant; À O inclinaison du pipeline; À g gravité; À S section du pipeline; À x fraction massique du composant dans une phase donnée; À p nombre de phases (de l à 3); À N nombre de composants; À k indice de phase; À i indice de composant; À m indice de mélange. À L'équation d'équilibre de masse pour chaque composant est: a TI S(pkRkXi) |+ a | s(pkRkX, Vk)| = 0 ( = 1 N). (1) À L'équation d'équilibre de la quantité de mouvement du mélange est: at TI S(pkRkvk)+ a (i S(PkRkVk) + p |= S(T,V -pmGsin0) (2) À L'équation d'équilibre d'énergie du mélange est: a 1I S(Pk Ek) | + a |E S (pkRkVk Hk) = S(Qw Pm gum sin 0). (3) Schéma numérique Le schéma numérique est conservatif et non dissipatif. Il fournit un bon bilan de masse et d'énergie en tout point de la conduite et à tout instant. Un schéma mixte implicite/explicite est utilisé pour optimiser le temps de calcul et la capacité de suivi des fronts des ondes de taux de vide, ce qui est particulièrement important dans le cas o la configuration du terrain suivi par le pipeline ou sa configuration propre, favorise la formation de bouchons de liquide avec propagation d'ondes de vide dans les deux directions opposées du pipeline, phénomènes que les spécialistes désignent couramment par "terrain slagging" ou "severe slugging". Module thermodynamique Le code de simulation TACITE comporte un "flash" thermodynamique intégré i.e. un sous-programme intégré de calcul des propriétés thermodynamiques (équilibre liquide-vapeur, composition de chacune des phases) à l 'ai de d'une équati on d'état. Ce flash réalise des calculs d'équilibre thermodynamique di et tri-phasiques pour des mélanges d'hydrocarbures incluant de l'eau. Les équations d'état cubiques de Peng-Robinson (1) et de Soave-Redlich-Kwong (2) sont utilisces pour modéliser les propriétés thermodynamiques à l'équilibre des phases. Ces équations sont définies dans les publications suivantes: 0 - D.-Y. Peng D. Y. et al, "A new two-constant equation of state", Ind. Eng. Chem. Fund. 15, 59(1976); Soave G., "Equilibrium constants from a modified RedlichKwong equation of state", Chem. Eng. Sci. 27,1197(1972). Dans les deux cas, les volumes moléculaires peuvent être corrigés par la méthode de Péneloux décrite dans la publication suivante: Péneloux A. et al, "A consistent correction for Redlich-Kwong-Soave volumes", Fluid Phase Equilibria 8, 7(1982). Le module dit flash assure un suivi précis de la composition des fluides aussi bien dans l'espace que dans le temps, durant toute la simulation. Ce suivi compositionnel rend le code TACITE particulièrement apte à prédire avec précision le risque de la formation d'hydrates, comme on le verra par la suite. Module hydrodynamique Le module hydrodynamique calcule le régime d'écoulement, la vitesse de glissement entre phases et les termes de frottement. La méthode selon l'invention La méthode selon l'invention permet de détecter en continu en tout point d'une conduite o circule un mélange polyphasique de fluides pétroliers, les conditions thermodynamiques de formation d'hydrates. Elle comporte l'utilisation s d'un module hydrodynamique mécanistique et d'un module thermodynamique compositionnel intogré pour définir les propriétés des phases et la résolution d'équations de conservation de la masse, de conservation de la quantité de mouvement et de transfert d'énergie au sein du mélange, en considérant que le mélange de fluides est sensiblement à l'équilibre à chaque instant, que la composition 0 du mélange polyphasique est variable tout le long de la conduite et que la masse de chaque constituant du mélange est définie globalement par une équation de conservation de la masse sans considération de son état de phase. Les fluides pétroliers sont regroupés en un nombre limité de pseudo-composants regroupant chacun certaines fractions d'hydrocarbures. La méthode est caractérisée en ce que l' on détecte les conditions thermodynamiques de formation d' hydrates après avoir opéré un regroupement particulier des fluides pétroliers en pseuJo-composants de façon à isoler les composants formateurs d'hydrates, avec définition pour chacun d'eux d'une fraction massique et d'un certain nombre de grandeurs physiques caractéristiques, et on applique aux modules thermodynamique et hydrodynamique, des donnces relatives à ces pseudo-composants sélectionnés de manière à déterminer en tout point la température (Td) de dissociation des hydrates. L' invention concerne aussi une méthode pour contr81er en continu et en tout point d'une conduite o circulent des fluides pétroliers polyphasiques, la forrnation 2s d' hydrates. Elle comporte la détection des conditions thermodynamiques de formation des hydrates comme défini cidessus, l' utilisation d' un organe de contr81e pour comparer la température des fluides pétroliers à cette température de dissociation (Td) et l' application commandée par cet organe de contrôle, de mesures destinées à combattre la formation des hydrates. Suivant un premier mode de mise en _uvre, on utilise des moyens de chauffage associés à la conduite pour élever la température des fluides pétroliers au dessus de la température de dissociation. Dans le cas o la conduite est incluse avec au moins un deuxième conduit dans un tube isolé du milieu extérieur, on utilise le deuxième conduit pour y faire circuler un fluide chaud. Suivant un deuxième mode de mise en _uvre, sur commande de l'organe de contrôle (C), on injecte des inhibiteurs d' hydrates dans la conduite. Présentation sommaire des figures D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés o: - la Fig.1 montre un organigramme comparé des opérations de simulation dans le cadre de l' application du code TACITE selon que l' on étudie ou non le risque de formation d'hydrates; la Fig.2 montre un organigramme de déroulement des opérations de simulation qui se répètent en tous les points d'un maillage de pipeline; - la Fig.3 montre schématiquement un mode de contrôle que l' on peut appliquer pour mettre en _uvre toute intervention visant à combattre la formation d'hydrates; la Fig.4 montre la section transversale d'un tube le long duquel courent plusieurs ensembles de tubes ou lignes de circulation de fluides; la Fig.S montre une comparaison des températures de dissociation d' hydrates calculées aux températures expérimentales dans le cas du méthane. L'histogramme est établi à partir de 158 points expérimentaux; - la Fig. 6 montre une comparaison des températures de dissociation d' hydrates s calculées aux températures expérimentales dans le cas de l'éthane. L'histogramme est établi à partir de 126 points expérimentaux; la Fig.7 montre une comparaison des températures de dissociation d' hydrates calculées aux températures expérimentales dans le cas du propane. L'histogramme est établi à partir de 93 points expérimentaux; 0 - la Fig. 8 montre une comparaison des températures de dissociation d' hydrates calculées aux températures expérimentales dans le cas de l'isobutane. L'histogramme est établi à partir de 43 points expérimentaux; - la Fig.9 montre une comparaison des températures de dissociation d' hydrates calculées aux températures expérimentales dans le cas du dioxyde de carbone. L'histogramme est établi à partir de 102 points expérimentaux; la Fig.10 montre une comparaison des températures de dissociation d'hydrates calculées aux températures expérimentales dans le cas de l'azote. L'histogramme est établi à partir de 57 points expérimentaux; - la Fig.11 montre une comparaison des températures de dissociation d'hydrates calculées aux températures expérimentales dans le cas du sulfure d'hydrogène. L'histogramme est établi à partir de 64 points expérimentaux; - la Fig. 1 2 montre une comparaison des températures de dissociation d' hydrates calculées aux températures expérimentales dans le cas de mélanges d' hydrocarbures. L' histogramme est établi à partir de 55 points expérimentaux provenant de 8 mélanges différents; la Fig.13 montre schématiquement un exemple de raccordement d'un pipeline avec un riser avec un point d'injection de gaz (pour gas lift); - la Fig.14 montre dans un diagramme bidimensionnel longueur de conduite temps, un exemple de domaine que l'on peut calculer par la méthode d' apparition s d'hydrates dans l'agencement de la Fig. 13;et la Fig.15 montre dans un diagramme bidimensionnel pression-température deux courbes locales d'apparition d'hydrates dans les conduites, selon que l'on tient compte ou non d'une modification de la composition du mélange en circulation. 0 DESCRIPTION DETAILLEE Présentation du module thermodynamique Structure des hydrates Les hydrates sont des composés solides constitués de molécules d'hydrocarbures (méthane, éthane, dioxyde de carbone etc.) piégées à haute pression et basse température à l'intérieur de cages approximativement sphériques formées par un réscau cristallin tridimensionnel de molécules d'eau. Cette structure est métastable en l' absence d'hydrocarbures et nécessite un taux de remplissage minimal pour assurer sa stabilisation. A la différence du réscau cristallin de la glace, celui des hydrates peut demeurer stable au-dessus de 0 C. Des mesures par diffraction aux rayons X réalisées dans les années 50 ont permis de mettre en évidence deux structures cristallines des hydrates selon la nature des molécules d'insertion présentes. Dans les deux cas, il s'agit de structures cubiques organisoes en réscaux de petites et grandes cavités. Le tableau ci-après montre l' occupation possible des différentes cavités selon la nature des molécules 2s formatrices d'hydrates mises en jeu (+: occupation possible de la cavité par la molécule; -: taille de la cavité insuffisante pour contenir la molécule). Structure I Struct Jre II Molécules Petite cavité Grande cavité Petite cavité Grande cavité N2 + + + + CO2 + + + + H2S + + + + nC4 + Le réseau de la structure I est de type cubique face centrce. Chaque maille contient 46 molécules d'eau organisces en: - 2 petites cavités de 3,95 À de rayon, s - 6 grandes cavités de 4,33 À de rayon. Le cristal de la structure II a le réscau du diamant. Une maille contient 136 molécules d'eau et 24 cavités: - 16 petites cavités de rayon 3,91 A, - 8 grandes cavités de rayon 4,73 . o Une troisième structure, la structure H. a été également identifiée. Son réseau cristallin, de type hexagonal, est constitué de 34 molécules d'eau et de 6 cavités: - 2 petites cavités de 3,91 À de rayon, - 3 cavités intermédiaires d'un rayon de 4,06 A, - 1 grande cavité de 5,71 À de rayon. Cette structure est beaucoup plus rare que les deux précédentes, elle se rencontre principalement dans le cas d'alcanes ou d'alcènes cycliques. Thermodvnamique de l'équilibre de phases en présence d'hydrates Tout calcul de prédiction d'hydrates est basé sur le critère d'équilibre de l'eau en phase hydrate et en phase aqueuse, c'est-à-dire sur l'égalité du potentiel chimique de l'eau dans la phase hydrate et du potentiel chimique de l'eau dans la phase aqueuse: aqueuse hydrate Ueau (P,T) IUeau (P,T). (1) Selon la théorie statistique de Van Der Waals et Platteeuw décrite dans la publication suivante:Van der Waals, J.H., et al., "Clathrate Solutions", Adv. Chem. Phys., 2,1, (1959), le potentiel chimique de l'eau dans la phase hydrate peut s'écrire: U hydrate (p T) = P hydrate (P. T) + RT2V, ln(1 ik ' (2) avec: - dUehayudrate-vide: potentiel chimique de l'eau dans le réseau cristallin en l'absence de molécule d'hydrocarbure, vi: rapport du nombre de cavités de type i (petite ou grande) sur le nombre de molécules d'eau dans un monocristal, - ik: probabilité pour qu'une cavité de type i soit occupée par une molécule hydrocarbonée de type k. D' après la théorie de l' adsorption de Langmuir, cette probabilité peut se calculer par: ik + C f ' (3) OUfk est la fugacité du composant k et Cik est le coefficient d' adsorption de Langmuir d'une molécule k dans une cavité de type i. Les coefficients de Langmuir peuvent être exprimés à partir du potentiel d'interaction entre les molécules d'eau et les molécules de gaz incluses dans les cavités de la structure cristalline d'hydrate: Cik = k J. exp:): r2dr, (4) o k est la constante de Boltzmann et W(r) l'énergie potentielle d' interaction entre les molécules d'eau formant la cavité et la molécule de gaz se trouvant à une distance r du centre de la cavité. I1 existe plusieurs formes de potentiel, les plus utilisées sont le potentiel de Lennard-Jones, le potentiel de Kihara et le potentiel à puits carré. Le potentiel de Lennard-Jones entre deux molécules s'écrit: jI2:6 avec: c: profondeur du puits, r: distance entre les centres des molécules, T: diamètre d'exclusion, r tel que W() = 0. Le potentiel de Kihara est identique à celui de Lennard-Jones, mais les molécules sont considérées comme des sphères de rayon non nul. Son expression analytique est la suivante: W(r) = oo Si r 2a (6) W(r) =4úL(râ2al2-(rï2a6 sir>2a avec: ú: profondeur du puits, - cr: tel que W(+ 2a) = 0, a: rayon moyen des 2 molécules, r: distance entre les centres des molécules. Enfin, le potentiel à puits carré est donné par: W(x)=oo si O Cc 0 Les différents paramètres figurant dans les expressions de ces potentiels sont ajustés sur des donnces expérimentales d'équilibre d'hydrates. Le potentiel chimique de l'eau dans la phase aqueuse s'écrit: Ueau (P,T) = eaU(p'T)+ RTln(x), (8) o y est le coefficient d'activité de l'eau, x la fraction molaire de l'eau dans la phase aqueuse et pe aU(P,T) le potentiel chimique de l'eau pure. C'est dans ce calcul de l'activité de l'eau en phase aqueuse que sont entre autres pris en compte les effets des sels et des alcools. En égalisant les équations (2) et (8), la condition d'équilibre recherchée (1) s'écrit: u au(P,T) -,UehaVudra'e-Vide(P,T)=RT>,vi ln1-1, ik-RTln(). (9) En fixant la pression, cette équation nous permet d'accéder à la température d'équilibre ou inversement, en fixant la température nous obtenons la pression d'équilibre. s Les différents modules de prédiction d'hydrates sont tous basés sur le formalisme décrit ci-dessus. Ces modules vont essentiellement se différencier dans le choix du potentiel d'interaction eau - hydrocarbure, dans le choix de l'équation d'état utilisée pour le calcul des fugacités ainsi que dans le module de calcul de l'activité de l'eau en phase liquide. Dans ce qui suit, nous présentons plus précisément le module que nous avons adapté pour les besoins de TACITE. Module implémenté dans TACITE Le module de prédiction d'hydrates implémenté dans TACITE s'inspire du module de Munck tel que décrit par Munck, J. et al, "Computations of the formation of gas hydrates", Chem. Eng. Sci 43, 2661 (1988). Il permet de calculer la température de dissociation des hydrates à pression donnée dans le cas d'un gaz, d'un liquide ou d'un mélange gaz - liquide diphasique en contact avec de l'eau liquide. Constantes de Langmuir: Ce module utilise un potentiel W(r) de type puits carré, ce qui conduit à une expression de la constante de Langmuir du composé k dans une cavité de type i sous la forme: Cik = T exp( T) (10) Les valeurs des paramètres A et B des différents composés formant des hydrates sont données dans le tableau ci-après: Petite cavité Grande cavité Molécule Structure A x 103 A x 103 B (K/atm) (K) (K/atm) (K) N2 1,617 2905 6,078 2431 II 0,1742 3082 18,00 1728 CO2 0,2474 3410 42,46 2813 0,0845 3615 851,0 2025 H2S I 0,0250 4568 16,34 3737 II 0,0298 4878 87,2 2633 C1 0,7228 3187 23,35 2653 II 0,2207 3453 100,0 1916 C2 0,0 0,0 3,039 3861 II 0,0 0,0 240,0 2967 C3 II 0,0 0,0 5,455 4638 iC4 II 0,0 0,0 189,3 3800 nC4 II 0,0 0,0 30,51 3699 Calcul des fuoacités: s Les fugacités des différents constituants présents dans le mélange, formateurs d' hydrates ou non, sont calculées à l' aide de l'équation d'état de Soave-Redlich Kwong déjà citée. Un tel calcul se fait par l'appel au module de flash à pression et température fixces intogré dans TACITE. Activité de l'eau en phase aqueuse: 0 Le module de Munck néglige la solubilité dans la phase aqueuse des hydrocarbures, de l'azote et du sulfure d'hydrogène. En revanche, la solubilité du dioxyde de carbone dans l'eau ne peut être négligée et est calculée à l'aide de la constante de Henry du CO2 dans l'eau: XCO2 = H (1 1) co2 o xcO2 est la fraction molaire de CO2 dans la phase aqueuse, fCO2 la fugacité du CO2 s et HCO2 la constante de Henry du CO2 dans l'eau donnce par: ln HCO2 = HA + TB + HC ln T + HDT. (12) Les valeurs des coefficients utilisés pour le calcul de la constante de Henry du CO2 figurent dans le tableau ci-après. HA HB HC HD (atm) (atm.K) (atm/lnK) (atm.K) ,27 -8764,5 -21,726 1,1055.10-4 o Constantes de référence: La différence de potentiel chimique entre l'eau liquide pure et l'eau dans la phase hydrate vide (premier membre de l'équation 9), peut s'écrire, moyennant quelques approximations: hydraievide O T P /Ueau (P,T)eau(P,T) _ /uo AHO +l\Cp(T-To) AV 3 - - dT+ _dP. (1) RT RTo To RT po RT 1S Dans cette équation, 0 représente la différence de potentiel chimique de l'eau dans le réseau hydrate vide et dans l'eau liquide à To = 273,15 K. IsHO est la différence d'enthalpie correspondante, ACP la différence de capacité calorifique et AV la différence de volume. La pression PO est la pression de vapeur à To; étant donné que PO est très petit comparé à P. on utilise PO = 0. La dépendance en température du terme PV est moyennce par: T = (T + 273,15)/2 (14) s AV et AHo sont tous deux considérés comme étant indépendants de la pression. AV est considéré aussi comme indépendant de la température, alors que la dépendance en température du terme enthalpique est prise en compte par le biais d'une différence de capacité calorifique constante /\Cp. Les différentes constantes figurant dans l'équation (13) sont présentées dans le tableau suivant: Propriétés Unités Structure I Structure II Ho J.moli 1264 883 M10 J.moli -485 8 -5201 / V cm3. moli 1151 808 ACp J.moli.K 39,16 39,16 La combinaison des équations (9) et (13) donne l'équation finale (15): 0 _ =0 P2( 0)dT + _ dP = ln()- Év; ln1E0ik À (15) Dans cette équation les inconnues sont la pression et la température. Nous avons arbitrairement choisi de fixer la pression et de calculer la température d'équilibre des hydrates correspondante. La recherche de cette température, c' est-à-dire la résolution de l' équation (15), s 'effectue de manière itérative par une méthode de Newton- Raphson partant d'une température initiale arUitrairement fixée à 0 C. Nous avons vérifié que la valeur de cette initialisation n'influe pas sur le résultat obtenu. Validation du module Nous avons validé le programme sur un ensemble de 700 points d'équilibre expérimentaux issus de 36 références bibliographiques différentes. Les résultats sont présentés sous forme d'histogrammes o figurent en abscisse, l'écart entre la température de dissociation calculée et celle mesurée, et en ordonnce, le pourcentage de points pour lesquels ces écarts de température sont compris entre les limites 0 indiquées en abscisse. Les points expérimentaux utilisés pour réaliser ces histogrammes correspondent à différents couples pression - température. Sur les figures 5 à 11, nous comparons ainsi les résultats du calcul avec les données expérimentales dans le cas de divers corps purs: méthane, éthane, propane, isobutane, dioxyde de carbone, azote et sulfure d'hydrogène. Nous constatons que sur l'ensemble des corps purs étudiés, le module utilisé donne d'excellentes prédictions des conditions d'équilibre des hydrates. Pour 99% des points, la température est prédite à +1 C. Dans tous les cas, l'écart maximal entre les températures de dissociation calculées et mesurées ne dépasse pas 1, 6 C. L'histogramme concernant les données de mélanges d'hydrocarbures, gaz naturels ou huiles, est présenté sur la Fig.12. Les données expérimentales utilisces pour calculer cet histogramme proviennent de huit mélanges différents, six de type gaz naturels et deux huiles. On observe à nouveau une très bonne description du comportement de ces systèmes par le module de prédiction utilisé. Simulation lobale avec ou sans calcul de risque d' apparition d'hvdrates L' impact se situe au niveau du regroupement des composants de base (lumping). Dans le cas standard, on cherche à en optimiser le nombre de façon à réduire le temps de calcul tout en gardant une bonne représentation des propriétés s d'équilibre et des propriétés des phases. Dans le cas o l'on cherche à prédire les hydrates, il faut pouvoir suivre la composition des composants appelés formateurs d'hydrates et donc, ces composants ne doivent pas étre regroupés avec d'autres. Outre les données se situant dans les fichiers de type ".PVT ", ".MAS " et ".TYP " décrits précisément ci-après, les différentes données d'entrée de TACITE sont regroupées dans 6 fichiers (cf. Fig. l): ".SCE ": description du scénario, c'est à dire de l'évolution dans le temps de la valeur des différentes conditions aux limites (débit d'entrée, pression de sortie, ouverture de vanne,...); - ".TOP ": description de la topographie de la conduite (élévation de la conduite lS en fonction de la longueur); - ".GEO ": description de la géométrie de la conduite (diamètre, longueur, porosité); ".THE ": description des isolations thermiques de la conduite; ".MSH ": description du maillage; - ".STO ": description du type de fichier de résultats et de la fréquence de stockage des résultats. Fichier ".PVT" Le fichier ".PVT" contient autant de lignes que de données thermodynamiques du fluide à indiquer, respectivement: - un entier qui définit le type d'équation d'état utilisée: 3 pour Peng-Robinson ou 4 pour Soave-Redlich-Kwong. - un entier qui définit le nombre de phases; il est mis à 1 et c'est le flash qui détermine leur nombre. s - le nombre de composants. - Boucle sur les composants avec indication, pour chacun d'eux, de leur: - masse molaire [kg/mol]. - température critique [K]. - pression critique [Pa]. lO - volume critique [m3/mol]. - facteur acentrique [-]. - parachor [((Nlm)/4m31mol)]. - volume molaire standard (inutilisé). - facteur de correction volumique [m3/mol]. lS - coefficients d'enthalpie idéale (tableau de 7 coefficients décrits ci-dessous). - Fin de la boucle. - la matrice des coefficients d'interaction binaire [-]. - la température de calage l [K], Température de calage 2 [K]. - la viscosité cinématique à la température l [m2/s], Viscosité cinématique à la température 2 [m2/s]. Fichier ".MAS" Ce fichier comporte également un certain nombre de lignes pour indiquer. - le nombre de composants initiaux. - la liste du numéro du groupe dans lequel le composant initial correspondant a été 2s mis. L'ordre est celui dans lequel se trouvaient les composants dans le fichier ".INP ". - le nombre de pseudo-composants. - les différentes fractions massiques des composants, ordonnées comme dans le fichier ".PVT". Fichier ".TYP" Ce fichier contient autant de lignes que de pseudocomposants définis dans le fichier ".PVT" et les informations sur les composants tiennent compte de l'ordre dans lequel les composants sont rangés dans le fichier ".PVT". La ième ligne concerne l'identificateur de la nature du ième composant. La convention suivante est retenue: Si le composant est ni l'eau, ni un formateur d' hydrates, son type est mis à 0. Si le composant est le sulfure d'hydrogène, son type est mis à 1. Si le composant est l'azote, son type est mis à 2. o Si le composant est le dioxyde de carbone, son type est mis 3. Si le composant est le méthane, son type est mis à 4. Si le composant est l'éthane, son type est mis à 5. Si le composant est le propane, son type est mis à 6. Si le composant est l'iso-butane, son type est mis à 7. Si le composant est le n-butane, son type est mis à 8. Si le composant est l'eau, son type est mis à 9. Le code de simulation TACll; utilise une méthode de résolution numérique de type à volumes finis pour simuler les écoulements dans les installations de production, les conduites étant discrétisces par exemple selon la méthode de maillage décrite dans la demande de brevet FR EN 00/08200 du demandeur. A tout instant, en tout point du maillage, les calculs de modélisation sont conduits selon l'organigramme de la Fig.2. / Le schéma numérique est alimenté par la donnée des termes de flux et les termes source, la modification des conditions aux limites (scénario). 2s Le schéma numérique fournit les masses des constituants, la quantité de mouvement du mélange, l'énergie interne du mélange (non-isotherme). La connaissance des masses et de l'énergie fournit, grâce au flash thermodynamique, la pression, la température et la composition massique des phases à l'équilibre. / Connaissant la pression et la composition massique, il est possible, si s l'utilisateur le souhaite (voir partie entourée d'un trait en pointillés), de: À calculer précisément grâce au module de prédiction d' hydrates, la température Td de dissociation des hydrates (température audessous de laquelle si des hydrates se sont formés, ils ne peuvent se dissocier); À d'utiliser un contr81eur C de type PID (Fig.3) pour comparer la 0 température calculée du fluide (Tfluide) à cette température Td de dissociation; et À d'agir si cette température de fluide est inférieure à la température de dissociation majorée d' un petit écart (pour anticiper le problème éventuel) sur l' un des moyens de contrôle prévu. Indépendamment de ce contrBle éventuel, le calcul continue, la connaissance de l'équilibre du mélange et des propriétés des phases, associée à la connaissance de la quantité de mouvement du mélange permet, grâce au module hydrodynamique, de fournir le glissement entre les phases et les frottements pariétaux et interfaciaux. I1 est alors possible de calculer les termes sources et les flux pour le prochain calcul. Contrôle en line Les opérations de lutte contre la formation d' hydrates sont généralement très coûteuses. La modélisation fine permise par la méthode selon l' invention, permet à l'exploitant de conna^tre à tout instant les conditions thermodynamiques régnant tout 2s le long du pipeline et donc de déterminer les endroits o les hydrates sont susceptibles de se former. Il est en mesure quand le besoin s'en fait sentir et seulement à ces moments-là, de mettre en _uvre des actions pour combattre les dépôts d'hydrates. Les nécessités d' intervention sont plus particulièrement critiques en cas d'arrêt momentané de la production ou quand on procède à des injections de gaz dans le mélange en circulation (aux emplacements de " gas lift ") par exemple. Un premier mode d' action consiste à réchauffer les fluides en circulation de manière à les amener au-dessus de la température de dissociation d' hydrates que l ' on a prédite, au moins aux endroits de la conduite o le risque est apparu. Les conduites pétrolières pour la production marine par grands fonds sont généralement incluses 0 avec d'autres lignes dans une conduite métallique T revêtue extérieurement d'une couche isolante. Cet ensemble de lignes (désigné par " bundles " par les gens de l' art) est réparti souvent en plusieurs faisceaux tels que F1, F2 (cf. Fig.4). Chaque faisccau comporte généralement une ou plusieurs lignes de transport 1, 2, 3 pour fluides monophasiques ou polyphasiques qui peuvent être rassemblées dans un tube commun 4 éventuellement recouvert lui-même extérieurement d'une couche isolante. L' intérieur de chaque tube 4 est rempli le plus souvent d' un fluide monophasique. L'espace intérieur de la conduite T. entre les différents faisceaux F1, F2, est rempli par exemple d'une mousse thermiquement isolante 5. L'une des lignes de transport 2, 3 peut être utilisée par exemple pour faire circuler si besoin est. un fluide de réchanffage. Pour modéliser les transferts thermiques entre les conduites pétrolières de l'ensemble de lignes et leur environnement (autres lignes, milieu extérieur) en tenant compte des moyens d' isolation thermique intercalés, on peut se reporter par exemple à la publication suivante: 2s - Duret, E. et al, "Pipeline bundles model implemented into a multiphase flow model", SPE 62948, SPE Annual Technical Conference and Exhibition, Dallas, Texas, 1-4 October (2000). Un deuxième mode d'action consiste à injecter dans les fluides des additifs inhibiteurs de type connu s ' opposant à la formation d' hydrates. Il s'agit d'un système de contrôle basé, non pas sur une mesure fournie par un capteur récl, mais sur une mesure fournie par un capteur dit " logiciel ". s Intérêt d'un code compositionnel pour un calcul d' hydrates L' utilisation d' un code compositionnel tel que le code TACITE permet en tout point du pipeline, un calcul exact, à la pression calculée du fluide et pour la composition locale calculée en cours de simulation, de la température de dissociation des hydrates, et par conséquent une bien meilleure simulation des risques 0 d' apparition d'hydrates que des codes o la composition est supposée fixe, correspondant à celle que l'on injecte à l'entrée de la conduite. On a pu vérifier par exemple que l'injection de gaz au pied d'un " riser " (technique bien connue des gens de l'art), modifiait les conditions thermodynamiques locales d'un mélange donné de fluides pétroliers, au point que l' on pouvait prédire la formation d'hydrates, alors qu'en supposant une composition fixe du mélange, ce risque n'était aucunement détecté. Nous illustrons ce point avec la simulation effectuée d'une injection de gaz au pied d'un riser se raccordant à un pipeline tel que celui schématisé sur la Fig.13, alimenté par un mélange dont la composition est donnée dans le tableau ci-dessous. Components Composition (% mol.) C, 10.53 C2 2.64 C3 2.63 IC4 2.63 NC4 2.63 IC5 2.63 NC5 2.63 NClo 13.16 NC20 7.90 H2O 52.62 La simulation est effectuce en considérant que le mélange initial a été modifié par une injection en bas du riser, de méthane favorisant la formationd'hydrates. Dans un diagramme longueur de conduite-temps, le domaine d' apparition d'hydrates que l'on obtient est schématisé à la Fig.14. s On a calculé et représenté sur le diagramme bidimensionnel pressiontemps de la Fig.15, deux courbes locales de formation d'hydrates, I'une (en pointillés) en considérant que la composition n'a pas été altérce par l'injection de méthane, comme le ferait un code de simulation non compositionnel, I'autre (en traits pleins) en tenant compte de la modification de composition intervenue. 0 On vérifie aisément sur la Fig.15, que le point qui y est figuré représentant le couple (pression, température) à l'état final serait considéré à tort comme sans risque puisqu'à droite de la courbe d'hydrates calculée à la composition initiale. Le domaine de risque d' apparition d'hydrates, présenté sur la Fig.14, serait également modifié en raison des variations de composition. REVENDICATIONS 1) Méthode pour détecter en continu en tout point d'une conduite o circule un mélange polyphasique de fluides pétroliers, les conditions thermodynamiques de formation d' hydrates dans laquelle on utilise un module hydrodynamique mécanistique et un module thermodynamique compositionnel intégré pour définir les propriétés des phases et l'on applique des équations de conservation de la masse et de la quantité de mouvement, et de transfert d'énergie au sein du mélange, en considérant que le mélange de fluides est sensiblement à l'équilibre à chaque instant, que la composition du mélange polyphasique est variable tout le long de la conduite 0 et que la masse de chaque constituant du mélange est définie globalement par une équation de conservation de la masse sans considération de son état de phase, et on regroupe les fluides pétroliers en un nombre limité de pseuJo-composants, caractérisce en ce que l' on détecte les conditions thermodynamiques de formation d'hydrates en: opérant un regroupement particulier des fluides pétroliers en pseudocomposants de façon à isoler les composants formateurs d'hydrates avec définition pour chacun d'eux d'une fraction massique et d'un certain nombre de grandeurs physiques caractéristiques; et - en appliquant aux dits modules, les donnces relatives à ces fractions particulières, de manière à déterminer en tout point la température (Td) de dissociation des hydrates. 2) Méthode pour contrôler en continu la formation d'hydrates, et en tout point d'une conduite o circule un mélange de fluides pétroliers polyphasiques, dans laquelle on utilise un module hydrodynamique mécanistique et un module 2s thermodynamique compositionnel intégré pour définir les propriétés des phases et l'on applique des équations de conservation de la masse et de la quantité de mouvement, et de transfert d'énergie au sein du mélange, en considérant que le mélange de fluides est sensiblement à l'équilibre à chaque instant, la composition du mélange polyphasique est variable tout le long de la conduite et la masse de chaque constituant du mélange est définie globalement par une équation de conservation de la masse sans considération de son état de phase, et on regroupe les fluides pétroliers s en un nombre limité de pseudo-composants, caractérisée en ce que: a) l'on détecte les conditions thermodynamiques de formation d'hydrates en: - opérant un regroupement particulier des fluides pétroliers en pseuJo- composants de façon à isoler les composants formateurs d'hydrates, avec définition pour chacun d'eux d'une fraction massique et d'un certain nombre de grandeurs 0 physiques caractéristiques; et en appliquant aux dits modules, les données relatives à ces fractions particulières, de manière à déterminer en tout point la température (Td) de dissociation des hydrates; b) on utilise d'un organe de contrôle (C) pour comparer la température des fluides pétroliers à cette température de dissociation (Td); et c) on applique sous le contrôle de cet organe de contrôle (C), des mesures destinées à combattre la formation des hydrates. 3) Méthode selon la revendication 2, caractérisoe en ce que l'on utilise des moyens de chauffage associés à la conduite (1) pour élever la température des fluides pétroliers au-dessus de la température de dissociation. 4) Méthode selon la revendication 3, caractérisée en ce que, la conduite (1) étant incluse avec au moins un deuxième conduit (2, 3) dans un tube (T) isolé du milieu extérieur, on utilise le deuxième conduit (2, 3) pour y faire circuler un fluide chaud. 2s s) Méthode selon la revendication 2, caractérisce en ce que, sur commande de