Brevet FR2753535A1 — publié le 19980320

La présente invention concerne une méthode pour réaliser un modèle de la cristallisation des paraffines dans un fluide pétrolier en fonction de sa composition. La majorité des bruts pétroliers comportent des hydrocarbures lourds qui présentent une phase cireuse solide aux températures suffisamment basses. La cire qui se forme se dépose sur les parois des puits et des pipelines avec pour résultat notamment une augmentation de la perte de charge. C'est vrai particulièrement pour les pipelines immergés dans des eaux froides dont la température est souvent inférieure à celle du point de trouble (cloud point) du fluide. Dans des situations extrêmes, les dépôts de cire peuvent aller jusqu'au blocage complet des conduites de production et de transport. I1 est connu d'ajouter aux fluides à acheminer des additifs chimiques ou des solvants susceptibles de contrarier les dépôts cireux mais ces procédés sont coûteux. Pour lutter contre ces phénomènes qui peuvent avoir des conséquences économiques considérables, il importe de bien maitriser les conditions de formation des dépôts de façon à mieux les prévenir et notamment de savoir calculer la température d'apparition des paraffines ainsi que I'évolution des dépôts en fonction de la température. La méthode selon l'invention a pour objet de former un modèle thermodynamique compositionnel permettant de déterminer pour un mélange d'hydrocarbures de composition donnée, - les conditions thermodynamiques (température T et pression P) à partir desquelles la cristallisation se produit, et notamment, pour une pression donnée, la température de cristallisation commençante, i.e. celle à partir de laquelle les premiers cristaux solides se forment lors du refroidissement du mélange; et - la quantité et la composition des dépôts formés pour T et P fixées. L'élaboration de ce modèle thermodynamique compositionnel destiné à décrire le comportement d'un brut paraffinique passe par: 1) la représentation analytique du brut par des pseudo-constituants; 2) la détermination des grandeurs physico-chimiques intervenant dans les différentes équations du modèle et qui sont caractéristiques de chacun des pseudoconstituants du fluide étudié; 3) le choix des modèles thermodynamiques représentant le comportement non idéal des phases liquide, gazeuse et solide. I1 est alors possible, si l'on possède des données mesurées au laboratoire, d'ajuster sur celles-ci, le paramètre d'interaction binaire en phase solide A. La connaissance de ce paramètre permet, ensuite, de calculer les courbes de dépôts solides lorsque la composition du mélange étudié varie, par addition d'un fluxant en particulier. ÉTAT DE LA TECHNIQUE H existe de nombreuses publications portant sur la formation des dépôts de paraffines dans les mélanges d'hydrocarbures. Parmi elles, on peut citer par exemple: - Narayanan, L. et al, "Model for Predicting Wax Deposition from Crude Oils," Paper 55a, 1993, AIChE Spring National Convention, March 28 - April 1; - Pedersen, K. S., "Prediction of Cloud Point temperatures and Amount of Wax Precipitation." paper SPE 27629. Soc. of Pet. Engrs., Richardson, TX (1993); - Won, K. W., "Thermodynamics for Solid Solution-Liquid-Vapor Equilibria: Wax Phase Formation from Heavy Hydrocarbon Mixtures", in Fluid Phas Eq. (1986) 30:265. Nombre de méthodes de modélisation connues comportent l'utilisation d'une base de données spécifique à partir de laquelle les propriétés thermodynamiques sont calculées. Les méthodes décrites dans ces publications antérieures surestiment toutes l'ordre de grandeur des dépôts solides mesurées expérimentalement avec des bruts. La méthode selon l'invention permet de former un modèle rendant compte précisément de la cristallisation des paraffines dans un fluide pétrolier en fonction de sa composition, dans laquelle on analyse le fluide pétrolier pour y définir les fractions massiques respectives d'un certain nombre de pseudo-constituants en nombre limité regroupant chacun certaines classes d'hydrocarbures. Elle est caractérisée en ce que l'on constitue une base de données regroupant les valeurs de paramètres physico-chimiques déterminés concernant un certain nombre N d'hydrocarbures constituants purs, et on détermine les paramètres physicochimiques de chacun de ces pseudo-composants par une combinaison des paramètres physico-chimiques correspondants de ces hydrocarbures. On détermine par exemple les paramètres physico-chimiques d'une partie au moins des pseudo-composants par pondération linéaire des valeurs d'un certain nombre k de paramètres physico-chimiques correspondants Gd(i,j) des hydrocarbures purs en fonction des compositions molaires xi en accord avec la relation: - T désignant le nombre d'hydrocarbures représentés par le pseudo-constituant P; - xi, le pour-cent molaire dans le brut du ième hydrocarbure (i = 1 ..., ) représenté par le pseudo-constituant P, - Gd(i,j), le jème paramètre parmi les k paramètres du ième hydrocarbure; - Gd (P, j), le jème paramètre parmi k paramètres du pseudo-constituant P. La méthode comporte par exemple la description du fluide pétrolier par 9 pseudo-constituants. On choisit de définir deux pseudo-constituants regroupant les fractions les plus lourdes qui sont représentées par deux molécules fictives définies chacune parune répartition molaire entre différents groupes qui les constituent, et on utilise une méthode de contribution de groupe pour le calcul de leurs propriétés thermodynamiques. Les molécules fictives sont définies par exemple pour les pseudo-constituants incluant tous les hydrocarbures dans la fraction C20+ du brut, à l'exception des nalcanes de cette même fraction. La méthode comporte par exemple l'assimilation des dépôts de paraffines formés à une phase solide unique non idéale modélisable par un paramètre qui représente les interactions binaires en phase solide. Comme on le précisera dans le cours de la description, les choix effectués dans la méthode selon l'invention, conduisent au calcul de quantités de dépôts paraffiniques solides en excellent accord avec celles mesurées expérimentalement. 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 est un schéma des opérations analytiques à réaliser sur un brut; - la Fig.2 est un tableau montrant un exemple de calcul des pour-cent massiques de pseudà-constituants P1 à P9 dans un brut connaissant leur pour-cent massique dans les fractions C20- ou C20+, et les résultats issus de la distillation du brut. - la Fig.3 montre une comparaison des courbes de dépôts solides calculées et expérimentales obtenues par une technique RMN; - la Fig.4 montre une comparaison des courbes de dépôts solides calculées et expérimentales obtenues par une technique dite ATD; - la Fig.5 montre une comparaison des courbes de dépôts solides calculées et expérimentales obtenues par la technique ATD dans le cas d'un brut fluxé par addition d'un condensat léger; - la Fig.6 est un tableau qui montre la nature et le nombre de groupes constitutifs de la molécule fictive du pseudo-constituant P8; et - la Fig.7 est un tableau qui montre la nature et le nombre de groupes constitutifs de la molécule fictive du pseudo-constituant P9. La modélisation de la cristallisation des paraffines selon l'invention est conduite en accord avec un ensemble de lois physiques bien connues des spécialistes et que l'on va rappeler ci-après. Pour calculer les équilibres gaz, liquide, solide, il faut choisir un modèle thermodynamique pour la représentation de chaque phase. A température et pression fixées, on cherche, pour un brut de composition donnée, le rapport molaire entre les différentes phases solide (S), liquide (L) et gazeuse (G). On écrit les conditions d'équilibre liquide (L)/gaz(G) et solide(S)/liquide(L) d'un système composé de n constituants en équilibre, qui sont à la base des calculs de dépôt solide des bruts paraffiniques. A l'équilibre thermodynamique, la fugacité fi du composant i doit être égale dans les trois phases: j;s= gL = 1,9 Les fugacités du liquide et du gaz du composant i s'expriment comme suit: LL L L L fiL(T,PXX, ) *(p (%P,x) (A3) fiG(T,P,xiG)=P*xiG*#iG(T,P,xiG) (A4) où xL, xG et XS représentent respectivement les fractions molaires de "i" dans les phases liquide, gazeuse et solide. Les coefficients de fugacité en phases gazeuse et liquide sont calculés par la méthode de contribution de groupes d'Abdoul et Péneloux décrite par exemple par: - Abdoul, W., Rauzy, E., Peneloux, A. (1991); "Group-Contribution Equation of State for Correlating and Predicting Thermodynamic Properties of Weakly Polar and Non-Associating Mixtures. Binary and Multi-component Systems". Fluid Phase Equilibria, 69, pp. 47-102, L'expression des fugacités en phase solide s'écrit, quant-à-elle, en fonction des coefficients d'activité ys Préf est une pression de référence. vi est le volume molaire partiel, défini comme suit (avec V, le volume du système considéré et ni, le nombre de moles du constituant "i"): On en déduit l'expression des constantes de partage liquide/gaz KLG et solide/liquide KiSL: Si l'on considère qu'il existe une transition de phase solide S1 -solide S2 on écrit l'équation A8 sous la forme connue: s-s fLs (T, Préf) - g1S.Ls Mifusion T Ahj ~~~~~ sion *(l- s'.on+ T fjs(T,Préf ) RT RT CLs - C51 Tfusion TS-S ~ Tfusion (A9), et Pi Pi,(Is * ln + R *tln Tj + T CLs ( CS2 rFS-S T Ts-s Pi Pi ~~~~~~~~~ R t T T ) T T fLsTP L i,réf)=P*i5(T,Préf) = P*cpLS (T P,f) (A10) OÙ Tfusion TiS-S, #hifusion, #hiS-S sont les températures et enthalpies de fusion et de transition solide-solide. CLS ,CS1 ,CS2 sont les capacités calorifiques d'un constituant CPi Pi Pi "i" respectivement en phase "liquide surfondu" et en phase solide sous ses formes cristallines Sl ou S2. On obtient ainsi, pour la constante de partage solide/liquide KSjL, l'expression suivante: Pseudo-constituants Le modèle réalisé selon l'invention étant de type compositionnel moléculaire, il nécessite la connaissance de la composition molaire du brut et donc d'une quantité considérable de données puisque dans un brut, on peut identifier jusqu'à environ 500 hydrocarbures différents. Il s'avère impossible dans la pratique de considérer l'ensemble de ces données car: 1) les grandeurs physico-chimiques (Tc, Pc, o, températures et enthalpie de fusion, etc.) concernant notamment les composés non "n-paraffiniques", ne sont pas connues dès que leur nombre d'atomes de carbone est supérieur à 10 ; et 2) le temps de calcul nécessaire serait trop grand pour les applications industrielles du modèle réalisé. Aussi l'ensemble des hydrocarbures, qu'il s'agisse de l'ensemble des constituants parfaitement identifiés par chromatographie en phase gazeuse, ou de tous ceux dont on ne connait que le nombre d'atomes de carbone et la famille chimique, a été regroupé en un certain nombre de groupes ou pseudo-constituants. On a choisi de préférence de représenter les fluides par 9 pseudo-constituants, 5 pour la fraction C20- et 4 pour la fraction C20+. Afin de déterminer la température de cristallisation commençante et la quantité de dépôt du brut, il est nécessaire de connaître un certain nombre k de grandeurs ou paramètres physico-chimiques propres à chacun des pseudoconstituants représentatifs du fluide étudié, au nombre de neuf par exemple: 1. masse molaire, MM; 2. température de fusion, Tfusion; 3. température de transition solide-solide, T55; 4. enthalpie de fusion, fusion; 5. enthalpie des transitions solide-solide, Ah55; 6. température critique, Tc 7. pression critique, Pc; 8. facteur acentrique, oe; 9. paramètre de conformation moléculaire li. La méthode selon l'invention comporte la constitution d'une base de données où ont été regroupés les paramètres ou grandeurs physico-chimiques définis cidessus, relatives à un nombre N d'hydrocarbures purs(N=75 par exemple). C'est à partir de ces hydrocarbures purs que sont calculés les paramètres physico-chimiques correspondants que l'on affecte à au moins une partie des pseudo-constituants. Pour chacun des pseudo-constituants de la fraction C20-, on calcule les k grandeurs Gd (P, j) selon la relation suivante: - T, le nombre d'hydrocarbures représentés par le pseudo-constituant P; - xi, le pour-cent molaire dans le brut du ième hydrocarbure (i = 1, ..., l) représenté par le pseudo-constituant P, - Gd(i,j), la dème propriété parmi les 9 citées ci-dessus du ième hydrocarbure; - Gd (P, j), la jème propriété parmi les 9 citées ci-dessus du pseudo-constituant P. Lorsque l'on a calculé les masses molaires de chacun des pseudo-constituants, on obtient leur pour cent molaire dans le brut à partir des pour-cent massiques. Fractionnement d'un brut à modéliser: Comme le montre le schéma de la Fig.l, on fractionne d'abord un brut à modéliser par une distillation au C20 permettant ainsi de séparer les deux fractions C20+ et C20- et en plus de connaître leur rapport massique. Chacune des fractions ainsi obtenue est analysée différemment. Analyse de la fraction C20 La fraction C20- est analysée par chromatographie en phase gazeuse (CPG), de préférence selon la méthode dite "Carburane" décrite par: - DURAND J.P, et al, "Direct and Automatic Capillary GC Analysis for Molecular Weight Determination and Distribution in Crude Oils and Condensates up to C20", J. of Hîgh Resolution Chromatography, 12, 1989, p. 230. Ce type d'analyse est partie intégrante du protocole d'analyse des fluides pétroliers développé dans le cadre de 1'ARTEP, et il a conduit à l'identification complète de tous les hydrocarbures jusqu'au nC1 1, de tous les n-alcanes jusqu'au ni20, à la répartition entre les iso-alcanes, les naphtènes et les aromatiques par fraction d'un atome de carbone entre le nu 1 1 et le nC15, et à une répartition entre les saturés et les aromatiques par fraction d'un atome de carbone entre le nC15 et le nu20. Cette fraction C20- est représentée par 5 pseudo-constituants P1, P2, P3, P4 et P5 : - P1 et P2 représentent les n-alcanes dont le nombre d'atomes de carbone est inférieur ou égal à 20; - P3 est le pseudo-constituant représentatif des iso-alcanes de la fraction C20-; - P4 est le pseudo-constituant représentatif des naphtènes de la fraction C20-; - P5 est le pseudo-constituant représentatif des aromatiques de la fraction C20-. Analyse de la fraction C20+ Les grandeurs analytiques à acquérir sur la fraction C20+ sont: la répartition des n-alcanes jusqu'au nC40 ou plus, la masse molaire de la fraction, et la répartition entre les saturés, les aromatiques, les résines et les asphaltènes. On a considéré comme négligeable l'influence des asphaltènes sur les dépôts paraffiniques. Pour la mise au point du modèle thermodynamique, on a choisi de mesurer: 1) la répartition des n-paraffines par chromatographie en phase gazeuse (CPG) jusqu'au nC40; 2) la masse molaire de la fraction C20+ par tonométrie; 3) la distribution entre les saturés, les aromatiques et les résines par chromatographie en phase liquide de la fraction C20+. On obtient alors, pour cette fraction, la S.A.R., c'est-à-dire, la répartition en pour-cent poids entre chacune de ces familles chimiques. Cette fraction C20+ est représentée par 4 pseudo-constituants P6, P7, P8 et P9: - P6 et P7 représentent les n-alcanes dont le nombre d'atomes de carbone est supérieur à 20; - P8 est le pseudo-constituant représentatif des saturés non n-paraffiniques (l'ensemble des iso-alcanes et des naphtènes) de la fraction C20+; - P9 est le pseudo-constituant représentatif des molécules non saturées (aromatiques, résines, asphaltènes) de la fraction C20+. Pour les pseudo-constituants P8 et P9, on dispose des informations analytiques suivantes: - répartition massique entre les différentes familles chimiques : saturés et aromatiques + résines; - motifs structuraux par RMN du 13C de chacune de ces familles. La méthode selon l'invention comporte aussi la création pour chacun des pseudo-constituants P8, P9, d'une molécule fictive, caractérisée par sa répartition molaire entre les différents groupes qui la constituent (cf. Fig. 6 et 7). Cette répartition permet de définir les pseudo-constituants P8 et P9 et de calculer leurs propriétés physico-chimiques (Tc, Pc, co, etc.), à condition de connaitre leur masse molaire (respectivement MMp8 et MMP 9) ainsi que leur nombre d'atomes de carbone (Cp8 et Cpg). Les résultats d'analyses des bruts rassemblés dans la base de données ont permis, outre la détermination des deux molécules fictives, d'établir une relation entre la masse molaire de P8 et celle de P9: MMp8 = l,33*MMp9 (5) Connaissant la masse molaire de la fraction C20+ : MMC20+, le pour-cent poids dans le brut du pseudo-constituant P8 "saturés non n-paraffiniques" et celui du pseudo-constituant P9 non saturés en C20+, on détermine les masses molaires respectivement MMpg et MMp9 par les équations (5) et (6): Pi : pseudo-constituant n"i (i= 1,.., 9); p(Pi) : pour-cent poids dans le brut du pseudo-constituant Pi; MMPi : masse molaire du pseudo-constituant Pi. Les données concernant les pseudo-constituants représentatifs des n-alcanes de la fraction C20+ étant connues (P6 et P7), les équations (5) et (6) permettent le calcul des inconnues MMp8 et MMP9. A partir de la formule brute des molécules fictives de P8 et P9, on peut accéder au rapport: nombre d'atomes d'hydrogène du pseudo - constituant Pi RPi = nombre d'atomes de carbone du pseudo - constituant Pi On obtient ainsi : RP8 = 1,88 et RP9 = 1,2332 On détermine les nombres d'atomes de carbone CP 8 et Cp9 selon les équations (7) et (8): Cp5 = MMc+MlMP8S*MMH (7), et MMp9 Cp9 = MMC + RP9*MMH (8) avec: MMC, la masse atomique du carbone (voisine de 12 g) et; MMH, la masse atomique de l'hydrogène (voisine de 1 g). Détermination des grandeurs physico-chimiques des pseudo-constituants P8 et P9 Pour les pseudo-constituants P8 et P9, qui ne sont pas générés à partir d'hydrocarbures identifiés, il est nécessaire de calculer les neuf propriétés physicochimiques et de les rentrer dans le fichier de données. Afin de calculer l'ensemble des grandeurs physico-chimiques de ces pseudo-constituants par des méthodes de contributions de groupes, il est nécessaire de calculer en premier la grandeur n"9, à savoir les paramètres du modèle d'Abdoul et Péneloux (1987) précédemment défini. A partir des nombres de groupes des molécules fictives, des nombres d'atomes de carbone CPi ( i = 8 ou 9) respectifs des pseudo-constituants Pi, le nombre N(Pi) des groupes constitutifs de chacun des pseudo-constituants est calculé par la relation: N(Pi) = Cpi/lOO*Nombre de groupes de la molécule fictive de Pi. Répartition des pseudo-constituants dans le brut On calcule pour les pseudo-constituants P1 à P5 leur pour-cent poids dans la fraction C20- et pour les pseudo-constituants P6 à P9, leur pour-cent poids dans la fraction C20+. Le pour-cent poids de l'ensemble des 9 pseudo-constituants dans le brut est calculé par pondération à l'aide des résultats de la distillation du brut. Le tableau de la Fig.2 montre un exemple où les pour-cent poids dans le brut des pseudo-constituants P 1 à P9 ont été calculés à partir des résultats suivants: P, 45% - Pc =55% 20 Ajustement du paramètre d'interaction binaire en phase solide Pour la mise en oeuvre de la méthode selon l'invention, on a choisi d'assimiler les dépôts formés à une phase solide unique non idéale et d'exprimer l'enthalpie libre d'excès de cette phase à l'aide d'une relation empirique dont le ou les paramètres peuvent être ajustés sur les résultats expérimentaux (en particulier, la courbe de dépôt solide). On a choisi pour cela l'équation à un paramètre proposée par Margules bien connu des spécialistes. Celle-ci présente l'avantage d'avoir l'expression la plus simple possible parmi celles restituant l'enthalpie libre molaire d'excès: où Ajj est le paramètre d'interaction binaire en phase solide entre les constituants i et gE j. La connaissance de la fonction g permet le calcul des coefficients d'activité en phase solide: On a considéré que les paramètres d'interaction binaire en phase solide sont tous identiques et égaux au paramètre A, soit: I1 en résulte l'expression suivante pour les constantes de partage solide/liquide: gSLs - fusion T AhS.-S T ~~~~~ hfusionTi) *(1 RT RT T.fusion + RT *(l- T5-5 (Al6) et où les ys sont donnés par les équations (A12) à (A14). KiLG Les constantes de partage KLG (équation A7) et KSL (équation A15) permettent d'accéder aux compositions molaires xi de chacun des constituants dans chacune des phases en équilibre et donc au rapport voulu (c'est-à-dire, au rapport en nombre de moles de la quantité de dépôts solides S au nombre total de moles du brut F). L'équation de modélisation à un paramètre A utilisée, ne nécessite le calage que d'un seul paramètre A qui représente les interactions binaires en phase solide. On a vérifié que les résultats calculés avec cette équation pour représenter la phase solide, sont en accord avec les résultats expérimentaux (Fig.3, 4). Ce paramètre est identique pour tous les binaires mis en jeu dans le mélange solide. La connaissance du paramètre A permet de calculer les courbes de dépôts solides des fluides dont la composition a légèrement varié (par addition d'un fluxant par exemple, comme l'illustre la Fig.5). L'ajustement de A est réalisé à partir de la courbe de dépôt solide mesurée au laboratoire avec le brut d'origine. H est effectué par une méthode des moindres carrés sur l'ensemble des points de la courbe Validation expérimentale La représentation analytique en neuf pseudo-constituants de la méthode selon l'invention, a été appliquée à un brut et, pour chacun d'entre eux, on a déterminé les grandeurs physico-chimiques (triplet caractéristique, température et enthalpie de fusion, etc.). Les résultats calculés ont été comparés avec ceux issus de deux techniques expérimentales: la R.M.N et 1'A.T.D (Analyse Thermique Différentielle), et les résultats comparatifs sont présentés dans les figures 3 et 4. On a une bonne restitution de la température de cristallisation commençante (Tcc) lorsque les calculs sont effectués avec une phase solide idéale (c'est-à-dire en prédictif). Les résultats calculés sur les quantités de dépôts solides sont très proches des les résultats expérimentaux. Les quantités calculées sont comprises entre les résultats issus des deux techniques expérimentales. On a préparé un condensat léger de 5 n-alcanes compris entre le nCg et le nCg. Puis, on a réalisé le mélange suivant: 80% poids du brut +20% poids de condensat léger. On a ensuite mesuré la courbe de dépôt solide du brut ainsi fluxé par ATD et comparé les mesures avec les résultats calculés lorsque: 1) a = 0, c'est-à-dire en mode prédictif (phase solide idéale), 2) a = paramètre de la phase solide non idéale, déterminé précédemment, lorsque celui-ci a été calé sur les résultats expérimentaux obtenus par ATD sur le brut analysé. L'ensemble de ces courbes est présenté dans la Fig. 5. Compte-tenu d'une erreur expérimentale relative de l'ordre de 15 %, on vérifie que les résultats calculés sur les quantités de dépôts avec une phase solide non idéale sont compris dans le domaine d'incertitude des mesures expérimentales. REVENDICATIONS 1) Méthode pour modéliser la cristallisation des paraffines dans un fluide pétrolier en fonction de sa composition, dans laquelle on analyse le fluide pétrolier pour y définir les fractions massiques respectives d'un certain nombre de pseudoconstituants (P1-P9) en nombre limité regroupant chacun certaines classes d'hydrocarbures, caractérisée en ce que l'on constitue une base de données regroupant les valeurs de paramètres physico-chimiques déterminés concernant un certain nombre N d'hydrocarbures constituants purs, et on détermine les paramètres physico-chimiques d'au moins une partie de ces pseudo-composants par une combinaison des paramètres physico-chimiques correspondants de ces hydrocarbures constituants. 2) Méthode selon la revendication 1, caratérisée en ce que l'on détermine les paramètres physico-chimiques d'une partie au moins des pseudo-composants par pondération linéaire des valeurs d'un certain nombre k de paramètres physicochimiques correspondants Gd (i, j) des hydrocarbures purs en fonction des compositions molaires xi, en accord avec la relation: - T, est le nombre d'hydrocarbures représentés par le pseudo-constituant P; - xi, le pour-cent molaire dans le fluide du ième hydrocarbure (i = 1 ..., représenté par le pseudo-constituant P, - Gd(i,j), lajème paramètre parmi les k paramètres du ième hydrocarbure; et - Gd (P, j), la jème paramètre parmi k paramètres du pseudo-constituant P. 3) Méthode selon l'une des revendications 1 ou 2, caractérisée en ce qu'elle comporte la description du fluide pétrolier par 9 pseudo-constituants. 4) Méthode selon la revendication précédente, caractérisée en ce que l'on définit deux pseudo-constituants (P8, P9) regroupant les fractions les plus lourdes, qui sont représentés par deux molécules fictives définies chacune par une répartition molaire entre les différents groupes qui les constituent, et on utilise une méthode de contribution de groupes pour le calcul de leurs propriétés thermodynamiques. 5) Méthode selon la revendication précédente, caractérisée en ce que l'on définit des molécules fictives pour les pseudo-constituants incluant tous les hydrocarbures dans la fraction C20+ du brut, à l'exception des n-alcanes de cette même fraction. 6) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on assimile les dépôts de paraffines formés à une phase solide unique non idéale modélisable par un paramètre (A) qui représente les interactions binaires en phase solide.