FR 2853101 A1 20041001 FR 0303908 A 20030328 La présente invention concerne une méthode dite de pseudoïsation, groupement ou lumping pour décrire un fluide formé d'un mélange de multiples constituants par un nombre réduit de constituants (au moins trois), et une méthode permettant l'opération inverse, méthode dite de dégroupement ou delumping , c'est-à-dire une méthode permettant l'obtention de la composition détaillée du fluide à partir de la description par un nombre réduit de constituants. La présente invention concerne donc une méthode de modélisation: - permettant d'estimer les propriétés ou le comportement des phases liquide et/ou vapeur de mélanges de multiples composants, tels que rencontrés in situ dans les réservoirs de pétrole ou de gaz ou en surface lors de l'exploitation de ces mêmes gisements, à l'aide d'une représentation à nombre réduit de pseudo-composants, - et offrant la possibilité de prédire en fonction du temps la composition détaillée de fluides produits au cours de la production. L'invention est particulièrement utile pour accélérer les calculs lors de la simulation de l'exploitation de gisements souterrains d'hydrocarbures. Un tel modèle permet à l'ingénieur de gisement de réduire les temps de calcul des simulations de comportement de gisement soumis à une exploitation, tout en conservant une bonne modélisation des mécanismes d'échanges entre phases hydrocarbure, et d'établir des profils compositionnels détaillés, ceux-ci étant nécessaires en particulier au dimensionnement prévisionnel et à la gestion d'installations de surface, telles que séparateurs, unités de traitement, lignes de transport, etc., donc utiles en ingénierie de surface et de procédés. ETAT DE LA TECHNIQUE La modélisation des écoulements dans un réservoir pétrolier ou dans un stockage sousterrain est fondée essentiellement sur l'application au réservoir préalablement maillé (ou à une portion de celui-ci) de la loi de Darcy bien connue décrivant l'écoulement de fluides en milieux poreux, de lois de bilans de matière dans chaque unité de volume, des relations thermodynamiques régissant l'évolution des propriétés de phase des fluides tels que la viscosité, la masse volumique, sur des conditions initiales, sur des conditions aux limites de fermeture de la structure, et sur des conditions aux puits. Le modèle dit Black Oil , désigné ci-après en abrégé par B.O., est l'un des plus utilisés en simulation pétrolière. Il permet de décrire un écoulement tridimensionnel et triphasique (eau-huile-gaz) compressible. Les effluents pétroliers intervenant dans ce modèle sont décrits généralement par un constituant eau, et deux constituants pour le fluide de gisement, le mot de constituant recouvrant ici la notion de composant (comme H20 pour l'eau) et la notion de pseudo-composant (regroupement de composants). Les constituants intervenant dans ce modèle sont au nombre de trois : un constituant eau (W), un constituant hydrocarbure lourd (H), et un constituant hydrocarbure léger (V). Dans un modèle de type B.O. dit strict , le constituant (W) n'est présent que dans la phase eau, le constituant (H) n'est présent que dans la phase liquide hydrocarbure (appelée huile ou condensat), le constituant (V) se partage entre les phases hydrocarbure liquide et vapeur (phase appelée gaz). Un modèle B. O. dit étendu diffère d'un modèle B. O. strict en ce que le constituant (H) se partage entre les phases hydrocarbure liquide et vapeur. Mais si l'usage des modèles B. O. est applicable à un grand nombre de cas industriels, il n'en reste pas moins à déconseiller dans certains cas, en particulier des cas de gisements de gaz à condensats soumis à une injection de gaz sec. Un autre modèle bien connu de simulation dit compositionnel est aussi utilisé, où les fluides hydrocarbure sont représentés par un nombre plus élevé de constituants, au moins trois, souvent davantage, le constituant eau n'étant présent que dans la phase aqueuse. La modélisation des flux de ces fluides mieux détaillés se traduit par des temps de calcul d'autant plus importants que le nombre de constituants est élevé. Pour que les calculs de modélisation puissent être menés dans un délai raisonnable, les fluides en place sont décrits comme constitués d'un nombre de composants ou pseudocomposants beaucoup plus réduit que le nombre réel de composants. Le passage d'une représentation détaillée des fluides à une représentation à nombre moindre de constituants s'effectue par des opérations dites de lumping (groupement ou regroupement) ou de pseudoïsation . Dans ce qui suit, sauf indication contraire, le terme de pseudoïsation sera employé pour toute méthode permettant de réduire le nombre de constituants. Différentes méthodes de pseudoïsation ont déjà été proposées pour sélectionner, définir les pseudo-composants, l'ingénieur se trouvant souvent à devoir trouver un compromis entre la précision, et le temps (et le coût) des calculs. Pour la simulation de l'exploitation de gisements de gaz à condensat soumis à une injection de gaz sec, ce sont des représentations à quelque 6 ou 8 constituants qui sont en général utilisées, ce qui se traduit par des temps de calcul d'autant plus importants qu'il est souvent nécessaire de réduire la taille des mailles pour limiter l'erreur numérique et par conséquent d'augmenter le nombre de mailles. Des efforts considérables sont donc encore aujourd'hui consacrés à l'élaboration de méthodes de pseudoïsation à l'usage de l'industrie, qui permettraient de réduire encore le nombre de constituants, tout en modélisant le comportement des fluides avec précision, et en offrant la possibilité de recouvrer une information compositionnelle détaillée. Les opérations qui permettent de prédire quels auraient été les résultats d'une simulation de réservoir si l'on avait utilisé une modélisation finement détaillée (où les fluides sont représentés par un nombre plus élevé de composants) sont connues des gens de l'art sous l'appellation de delumping , que l'on désignera par la suite aussi par éclatement . Par le brevet WO 00/37898, on connaît une méthode de pseudoïsation, applicable à des simulations compositionnelles, basée sur la sélection d'un nombre de composants de base 'dominants' égal au nombre de pseudo-composants que l'on veut obtenir à l'issue de la procédure. Dans cette méthode, les composants non dominants sont représentés dans tous les pseudo-composants, et un composant dominant particulier n'est représenté que dans un seul pseudo-composant. La transformation mathématique à la base de la méthode de pseudoïsation permet, par transformation inverse, de recouvrer l'information compositionnelle détaillée. La composition d'un pseudo-composant pris en particulier peut faire apparaître des fractions molaires négatives de constituants de base, comme dans l'exemple donné table E7 du brevet en référence. On conçoit qu'une telle représentation perd un certain sens physique dès lors que l'on considère isolément un pseudo-composant particulier. Cela semble pouvoir poser un problème de robustesse, lorsque, comme il arrive en pratique dans les cas d'injection de gaz, un résultat 'local' de la simulation, par exemple dans certaines mailles, indique la disparition de l'un ou plusieurs des constituants utilisés dans la simulation dynamique. Par ailleurs, la mise en #uvre de l'invention est décrite comme nécessitant de nombreux calculs itératifs et un large espace mémoire. Par des publications antérieures, on connaît des méthodes de pseudoïsation, applicables également à des simulations compositionnelles, où chaque pseudo-composant est formé en regroupant plusieurs constituants de base, un constituant de base particulier ne se trouvant représenté que dans un seul pseudo-composant. Le regroupement peut être effectué selon un choix de critères fixés a priori, tels que ceux cités dans le brevet mentionné ci-dessus, ou par une procédure d'optimisation, par exemple telle que proposée par K. Liu dans le papier Reduce the Number of Components for Compositional Réservoir Simulation , SPE 66363, présenté à la SPE Réservoir Simulation Symposium, Houston, Texas, 11-14 Février 2001. Un papier, bien connu des gens de l'art, écrit par D. E. Kenyon et G. Alda Behie Third SPE Comparative Solution Project : Cycling of Retrograde Condensate Reservoirs , SPE 12278, Journal of Petroleum Technology, Août 1987, illustre une situation non exceptionnelle : il arrive de constater, pour un même cas d'étude, une grande disparité dans les résultats de simulations compositionnelles lorsqu'elles sont réalisées avec des représentations de fluides différentes, et qui plus est, avec des logiciels de simulation différents. Dans ce papier, le nombre des constituants varie de 5 à 16, et la disparité des résultats est en partie imputable aux divers choix de représentations compositionnelles des fluides : la tendance générale observée est que plus le nombre de constituants est faible, plus la saturation en liquide hydrocarbure (appelé huile ou condensat) dans une maille particulière peut être sous-estimée, et plus la récupération d'huile en surface est alors surestimée. Pour recouvrer l'information compositionnelle détaillée au cours d'une simulation compositionnelle de réservoir, on peut utiliser des méthodes de delumping telles que celle exposée dans le brevet WO 99/42937 et aussi dans le papier de C. Leibovici et J. Barker A Method for Delumping the Results of a Compositional Réservoir Simulation SPE 49068, présenté à la SPE Annual Technical Conference and Exhibition New-Orleans, 27- 30 Septembre 1998. La méthode permet de prévoir l'évolution de la composition détaillée au cours du temps à partir des calculs, en particulier de coefficients d'équilibre, effectués dans une simulation de type compositionnel de fluides décrits par un certain nombre réduit de pseudo-composants, le nombre de composants étant au moins de trois. Une représentation de type Black-Oil peut être considérée comme issue d'une opération de pseudoïsation particulière fournissant deux pseudo-composants. La composition détaillée de chacun des pseudo-composants, qui n'est pas utile à la construction de la représentation, n'est pas connue a priori, sans que ceci soit un obstacle rédhibitoire au recouvrement d'une information compositionnelle détaillée, par une opération de delumping. Ainsi, par le brevet FR 00/09008, on connaît une méthode qui permet de prévoir l'évolution de la composition détaillée au cours du temps à partir des calculs effectués dans une simulation dynamique de type Black-Oil. Le principe de l'étape d'éclatement (ou delumping) des brevets WO 99/42937 et FR 00/09008 est de déterminer, - à partir des calculs effectués pendant la simulation, avec la représentation thermodynamique regroupée (compositionnelle ou BO), dans chaque maille et à chaque pas de temps-, le coefficient #D0 et les n coefficients #Dp (soit n+1 coefficients, n étant le nombre de paramètres de l'équation d'état) d'une équation générale connue, publiée antérieurement dans un papier de C.F. Leibovici, E.H. Stenby, K. Knudsen, A Consistent Procedure for Pseudo-Component Delumping , Fluid Phase Equilibria, 1996, 117, 225-232. : (1) Ln(Ki) = #D0 + # #Dp #pi p=1 où K; est la constante d'équilibre du constituant i et les #pi sont des paramètres fixes de caractérisation du constituant i dans l'équation d'état pour une représentation thermodynamique donnée. Une fois le coefficient #D0 et les n coefficients #Dp déterminés, on les utilise pour calculer les constantes d'équilibre des constituants de la représentation thermodynamique détaillée (Nrb composants) en appliquant l'équation (1) aux Nrb composants avec leurs propres paramètres fixes de caractérisation dans la représentation thermodynamique détaillée. L'un des intérêts de cette méthode est que, dans l'étape d'éclatement ( delumping ), il n'est pas nécessaire de résoudre les Nrb équations d'équilibre associées à l'équation d'état (équations qui traduisent l'égalité des fugacités de chaque constituant dans chaque phase) aux différents pas de temps de la simulation dynamique des écoulements, ce qui constitue une économie de temps de calcul. Le papier de W. H. Goldthorpe Simulation of Gas Injection Processes in Gas-Condensate Reservoirs Using a Binary Pseudo-Component Représentation , SPE 19470, présenté à la SPE Asia-Pacific Conference, Sydney, Australie, 13-15 Septembre 1989, illustre les résultats de simulation que l'on peut obtenir avec une modélisation Black-Oil avancée dans un cas d'exploitation d'un réservoir de gaz à condensat par un procédé d'injection de gaz. Le cas pris comme exemple est issu de la publication ci-dessus mentionnée de D.E. Kenyon et G. Alda Behie. Au vu de la disparité des résultats dans cette publication, les résultats de W. H. Goldthorpe avec une représentation Black-Oil, jugée par comparaison avec les résultats d'une simulation effectuée avec une représentation détaillée, apparaissent bien plus cohérents, mais on peut constater dans la figure 8 du papier que la solution n'est pas satisfaisante dans la phase de re vaporisation car la saturation en huile dans une maille particulière (la même que dans le papier de référence) est très différente de celle de la prédiction compositionnelle détaillée, et semble par ailleurs être tronquée de valeurs négatives pendant huit années d'exploitation simulée. L'état de la technique antérieure est donc que l'on ne dispose pas de méthode de pseudoïsation à la fois simple et robuste: permettant de réduire à trois le nombre de pseudo-composants dans des simulations compositionnelles, pour bénéficier de temps de calcul notablement réduits, notamment pour les simulations de cas d'injection de gaz avec effets de re vaporisation, difficilement traitables en Black-Oil, - garantissant un sens physique aux résultats des simulations compositionnelles et, par voie de conséquence, aux résultats d'opérations associées de delumping. LA METHODE SELON L'INVENTION La méthode de pseudoïsation selon l'invention permet d'estimer les propriétés ou le comportement des phases hydrocarbure liquide et/ou vapeur, à partir de données relatives à un ensemble de référence constitué de mélanges hydrocarbonés dans une série d'états thermodynamiques découlant de conditions d'exploitation rencontrées, ou prévues comme telles, pour des gisements souterrains d'hydrocarbures. Elle comporte les étapes suivantes : - on regroupe chacun des dits mélanges hydrocarbonés en au moins trois constituants (V, I, H), aucun de ces constituants ne correspondant à une sélection particulière de composants ou pseudo-composants de base qui serviraient à une description compositionnelle détaillée des fluides, en considérant que les phases gaz issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le troisième constituant (H) est exclu, et que les phases huile issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le premier constituant (V) est exclu ; - on détermine par bilan de matière les compositions des produits de séparation comprenant, pour les produits gazeux, au moins le premier et le deuxième constituant (V, I) en proportions variables, et pour les produits liquides, au moins le deuxième et le troisième constituant (I, H) en proportions variables ; et - on détermine la composition à au moins trois constituants de chaque mélange hydrocarboné de l'ensemble de référence par combinaison des produits de sa séparation au prorata des quantités de chaque produit de séparation. Suivant un mode de mise en #uvre, on regroupe chacun des mélanges hydrocarbonés en seulement trois constituants (V, I, H), les phases gaz issues de la dite séparation sont des mélanges en proportions variables du premier constituant (V) et du deuxième constituant (I), on regroupe les phases huile issues de la dite séparation sont des mélanges en proportions variables du deuxième constituant (I) et du troisième constituant (H), et on détermine la composition à trois constituants. Les conditions de surface sont celles rencontrées ou prévues lors de l'exploitation du gisement mais elles sont redéfinissables selon le contexte de l'étude menée, et peuvent donc être différentes des conditions de surface rencontrées ou prévues lors de l'exploitation du gisement. Suivant un détail de mise en #uvre, le bilan de matière est un bilan massique et on attribue à chacun des trois constituants (V, I, H) une masse molaire après une analyse quantitative des masses molaires des produits de séparation de l'ensemble de référence. Suivant une variante de mise en #uvre, on définit les données nécessaires aux calculs d'équilibre et à la modélisation des propriétés de phases dans la représentation groupée, en utilisant les compositions des phases dans la représentation groupée et des données connues, ou estimées, a priori portant au moins sur la densité et la viscosité de phases huile et gaz en équilibre appartenant à l'ensemble de référence. Lorsque les dites données incluent des données compositionnelles détaillées des phases préalablement représentées par une équation d'état détaillée , on définit les paramètres d'une première (le cas échéant unique) équation d'état de la représentation groupée, utile à la modélisation des propriétés de phases, en utilisant ces données compositionnelles. Suivant une autre variante préférée de mise en #uvre, on ajuste les paramètres d'une seconde équation d'état de la représentation groupée, utile aux calculs d'équilibre, pour reproduire les coefficients d'équilibre de la représentation groupée. Pour réaliser cet ajustement, on utilise par exemple les paramètres par constituant de la représentation groupée dans une équation d'état utile aux calculs des propriétés des phases. Suivant une autre variante de mise en #uvre utile en particulier à des calculs ultérieurs d'éclatement ou delumping , on détermine les coefficients d'équilibre des fluides dans une représentation compositionnelle détaillée, à partir de variables et/ou paramètres intervenant dans le calcul des propriétés de phases, dès lors que les paramètres utiles aux calculs des propriétés de phases dans la représentation groupée ont été estimés de façon à reproduire les paramètres des phases dans l'équation d'état de la description compositionnelle détaillée. La méthode peut également comporter des étapes d'éclatement ou delumping pour prédire en fonction du temps et dans au moins une zone thermodynamique, une composition détaillée d'un fluide contenu dans un gisement d'hydrocarbures ou produit par au moins un puits, ces étapes étant par exemple les suivantes : - on représente le réservoir sous la forme d'un réseau de mailles (m) dont chacune constitue un volume élémentaire rempli de fluides sous forme de une ou plusieurs phases, avec au moins une phase non aqueuse ; - on définit, pour chaque zone ou domaine thermodynamique, les fluides par une représentation détaillée de base, de façon à déterminer la quantité de chaque constituant de base (i) dans chaque phase hydrocarbure dans chaque maille (m) à l'instant défini comme initial pour le calcul de delumping ; par zone thermodynamique pour laquelle on adopte une représentation groupée des fluides, on détermine une équation d'état construite, préalablement à la simulation dynamique de réservoir avec la représentation groupée, pour reproduire les paramètres de phases, dans l'équation d'état de la représentation détaillée, des fluides hydrocarbonés au cours de chemins thermodynamiques tenus pour représentatifs de ceux que vont suivre les fluides hydrocarbonés pendant la simulation maillée ; - on réalise à un pas de temps t une simulation compositionnelle avec un nombre limité de constituants où les propriétés des phases sont calculées par une équation d'état, ladite simulation permettant de calculer au moins dans chaque maille (m) et à des pas de temps consécutifs, une pression pour une phase hydrocarbure, la température lorsque celle-ci varie, les débits des phases entre mailles, entre mailles perforées et puits, et les valeurs de paramètres et/ou propriétés de phase intervenant dans l'expression formelle des coefficients d'équilibre de la représentation détaillée, et on mémorise ces diverses quantités ; - on estime au pas de temps suivant (t+1) la fraction molaire de chaque constituant i dans la composition détaillée globale du fluide hydrocarboné dans la maille (m) par bilan matière sur la maille (m) ; on détermine en utilisant les quantités mémorisées, au même pas de temps (t+1) et dans chaque maille (m), les coefficients d'équilibre de chaque constituant (i) dans la représentation détaillée ; - on détermine, au même pas de temps (t+1), la fraction vaporisée dans chaque maille (m) ; et - on estime la composition détaillée de chaque phase hydrocarbonée, au même pas de temps (t+1) et dans chaque maille (m) en particulier sous la forme des nombres de moles de chaque constituant i dans la représentation détaillée dans chaque phase au même pas de temps (t+1). En d'autres termes, on peut aussi définir l'essentiel de la méthode de pseudoïsation comme suit : - Ayant à disposition des données, au moins de densités et viscosités, concernant des phases huile et gaz en équilibre in situ pour un certain nombre d'états thermodynamiques 'e' dans des conditions de pression et température Pe et Te, - connaissant directement ou par l'intermédiaire de calculs préliminaires, les masses molaires des phases 'huile' en provenant de la séparation de l'huile in situ, de la séparation du gaz in situ, notées respectivement (MMoo) et (MMGOe), les masses molaires des phases gaz en surface provenant de la séparation de l'huile in situ, de la e e séparation du gaz in situ, notées respectivement (MMOGe) et (MMGGe), et les fractions molaires vapeur (ou fraction molaire de phase gaz), #OOe pour la séparation de l'huile in situ, #GOe pour la séparation du gaz in situ, - on détermine les compositions des phases en surface , et par voie de conséquence les coefficients d'équilibre, dans la représentation groupée utilisant les constituants (V, I, H) avec l'hypothèse que les phases gaz en surface peuvent contenir les constituants (V) et (I), à l'exclusion du constituant (H), et que les phases huile en surface peuvent contenir les constituants (I) et (H), à l'exclusion du constituant (V). - on détermine les compositions des phases in situ dans la représentation groupée utilisant les constituants (V, I, H) par recombinaison des compositions des phases en surface en utilisant les fractions molaires vapeur #OOe, #GOe, ces déterminations fournissant les coefficients d'équilibre des constituants (V, I, H) dans les conditions de réservoir in-situ , - on utilise les données connues a priori et les compositions des phases dans la représentation groupée pour définir les entrées nécessaires à la modélisation des propriétés de phases et aux calculs d'équilibre dans la représentation groupée et, plus particulièrement lorsque les données a priori incluent des données compositionnelles détaillées des phases préalablement représentées par une équation d'état, on utilise ces informations additionnelles pour définir les paramètres d'une équation d'état de la représentation groupée utile en premier lieu à la modélisation des propriétés de phases, le cas échéant pour définir les paramètres d'une seconde équation d'état de la représentation groupée utile aux calculs d'équilibre. Les constituants V, I, H s'apparentent respectivement à un constituant volatile, un constituant intermédiaire et un constituant lourd, sans qu'aucun d'entre eux ne corresponde à une sélection particulière de composants ou pseudo-composants de base qui serviraient à une description compositionnelle détaillée des fluides. L'idée essentielle mise en #uvre pour le dégroupage ou delumping des simulations compositionnelles est de calculer les coefficients d'équilibre des fluides dans la représentation compositionnelle détaillée à partir de variables et/ou paramètres intervenant dans le calcul des propriétés de phases dès lors que les paramètres utiles aux calculs des propriétés de phases dans la représentation groupée ont été estimés de façon à reproduire les paramètres des phases dans l'équation d'état de la description compositionnelle détaillée PRESENTATION 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 présentés où : - la figure 1 schématise la séparation distincte d'une phase huile et d'une phase gaz en équilibre dans le réservoir ; - les figures 2 à 5 comparent les résultats d'une simulation compositionnelle avec une représentation à trois pseudo-composants (résultats désignés par 'ternaire'), obtenue par application de la méthode selon l'invention, à ceux d'une simulation compositionnelle détaillée à seize constituants de base (résultats désignés par 'référence ; - les figures 6 et 7 montrent, pour les représentations détaillée et ternaire, respectivement les variations du terme d'attraction et du co-volume (exprimés sous forme adimensionnelle) des phases huile et gaz en fonction de la pression au cours de la simulation thermodynamique d'une opération de dépressurisation à volume constant à la température de réservoir ; - les figures 8. 1 à 8. 16 d'une part, et les figures 9.1 à 9. 16 d'autre part, montrent, par constituant, respectivement la composition détaillée du gaz et de l'huile, au cours de la même opération de dépressurisation, telles qu'obtenues à partir de la représentation de référence, et à partir de l'éclatement de la représentation à trois pseudo-composants ; - la figure 10 récapitule les résultats illustrés aux figures 8. 1 à 8. 16 et 9. 1 à 9. 16, en montrant l'erreur absolue maximale obtenue sur la teneur en chaque constituant durant la dépressurisation de 237. 4 à 69 bars en fonction de la teneur du constituant dans la phase à la pression de rosée du fluide initial ; - les figures 11. 1 à 11. 16 d'une part, et les figures 12. 1 à 12. 16 d'autre part, montrent, par constituant, respectivement la composition détaillée du gaz et de l'huile, au cours de la simulation thermodynamique d'une opération de vaporisation d'un condensat par injection d'un gaz sec à la pression constante de 169 bars, telles qu'obtenues à partir de la représentation de référence, et à partir de l'éclatement de la représentation à trois pseudo-composants ; et - la figures 13 récapitule les résultats illustrés aux figures 11.1à 11.16 et 12.1à 12. 16, en montrant l'erreur absolue maximale obtenue sur la teneur en chaque constituant pour l'opération de vaporisation en fonction de la teneur du constituant dans la phase au démarrage de l'opération de vaporisation. DESCRIPTION DETAILLEE Préambule La modélisation du comportement de mélanges fluides nécessite un mode de définition des constituants utilisés pour décrire la composition des fluides et comporte : - des calculs d'équilibre thermodynamique, pour savoir si, dans des conditions données de pression et de température, l'état d'équilibre d'un mélange, caractérisé par une composition globale, est monophasique, ou si le mélange se sépare en différentes phases, le calcul d'équilibre permettant alors de connaître les proportions de chaque phase, et la composition de chaque phase. Un calcul d'équilibre sur un mélange, dont la composition globale est décrite par les fractions molaires zi de chaque constituant (i), conduit à la détermination de la fraction vapeur # et à la détermination de la composition de chaque phase. La fraction vapeur # est déterminée en général en résolvant l'équation de Rachford-Rice connue des spécialistes qui est appliquée par exemple dans le document suivant : - Rachford, H.H.Jr et Rice, J.D.; Procédure for Use of Electronic Digital Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium", J. Pet. Technol., 1952,14, 19, (2) E zi (Ki - 1)~ = (2) i 1+(Ki-1)0 les compositions des phases 'huile' et 'gaz', décrites respectivement par les fractions molaires xi et yi, étant déterminées par : Zi Ki Zi (3) Xi = 1 + (Ki - 1) e; Yi = 1 + (Ki - 1) e Une phase hydrocarbure est dite sous-saturée à une pression et une température T donnée lorsqu'elle ne se trouve pas en équilibre avec une autre phase hydrocarbure. Selon une convention souvent admise par les gens de l'art, une phase hydrocarbure sous-saturée est dite 'huile' (ou 'condensat') lorsque sa température critique est inférieure à la température T, et est dite 'gaz' lorsque sa température critique est supérieure à la température T. Lorsque sa température critique est égale à la température T, la phase peut être qualifiée, en principe indifféremment, de 'huile' ou 'gaz', et est donc considérée comme rattachée à l'une ou l'autre des qualités 'huile' ou 'gaz'. - des calculs de propriétés des phases telles que masse volumique, viscosité, - le cas échéant, des calculs, dits de delumping, permettant la restitution d'une information compositionnelle détaillée lorsque les fluides ont été modélisés de façon non détaillée. Les différentes étapes de la méthode de modélisation proposée sont présentées ci-après, la présentation étant structurée en cinq sous-chapitres: - mode de définition des constituants, - propriétés des phases, - équilibres thermodynamiques, delumping, - remarques concernant la généralisation de l'approche. Mode de définition des constituants a) On considère un ensemble E , dit de référence, de fluides non aqueux, dits mélanges hydrocarbure même s'ils peuvent contenir certains composants autres que des hydrocarbures tel que l'azote, le dioxyde de carbone, l'anhydride sulfureux, susceptibles de correspondre : - aux phases fluides hydrocarbure sous-saturées trouvées, supposées ou prédites dans une ou plusieurs zones thermodynamiques du réservoir considéré ; - aux phases fluides hydrocarbure en équilibre trouvées, supposées ou prédites dans une ou plusieurs zones thermodynamiques du réservoir considéré ; - aux phases fluides hydrocarbure correspondant aux fluides hydrocarbure injectés (ou supposés l'être) dans le réservoir. Dans les cas pratiques les plus fréquents, les fluides d'injection sont des gaz contenant principalement des hydrocarbures légers et intermédiaires. b) Lorsque l'on dispose des données nécessaires, on définit une représentation détaillée à Nrb composants et/ou pseudo-composants de base pour caractériser les mélanges fluides. Cette situation étant très courante, on place la suite de l'exposé dans une telle situation car cela permet d'introduire progressivement les conventions d'usage utiles à la description de la méthode. Le nombre Nrb étant typiquement supérieur à 10, chaque constituant est décrit par un certain nombre de paramètres, dont sa masse molaire MMi, l'indice i étant le numéro attribué au constituant ( i=l,..., Nrb). Les fractions molaires Xi et yi décrivent respectivement la composition des phases hydrocarbure liquide ('huile') et vapeur ('gaz') dans la représentation de base. Un élément 'e' de l'ensemble E est alors décrit par : pression P , une température e T, et : e e e - (2 x Nrb) fractions molaires (xi .. xNrb) et (y1 .. yNrb) lorsque deux phases hydrocarbure en équilibre coexistent, - ou (Nrb) fractions molaires (x1e .. xNrbe) ou (y1e .. yNrb) lorsque une seule phase hydrocarbure (huile ou gaz) sous-saturée est présente. Dans ce qui suit, on écrit les équations pour un élément 'e' où deux phases liquide et vapeur coexistent. Pour un élément 'e', où une seule phase sous-saturée est présente, il suffit de sélectionner le jeu d'équations relatives à la seule phase représentée. c) On considère que l'on dispose, pour chaque phase de l'ensemble E de référence, de données de fluides susceptibles de correspondre aux compositions des phases fluides hydrocarbure trouvées, supposées ou prédites correspondant aux mélanges fluides produits par une chaîne de séparation, la chaîne de séparation pouvant être réduite à une seule détente ou comporter un étage ou plusieurs étages intermédiaires de séparation. Une chaîne de séparation est ici, comme il est d'usage dans l'ingénierie de gisement, considérée comme définie par les valeurs de pressions et températures des étages de séparation. Il est pratique de considérer la chaîne de séparation existante ou prévue par laquelle le mélange fluide est ou est supposé être produit. On peut être amené, dans des cas pratiques, à considérer différentes chaînes de séparation, par exemple une chaîne de séparation différente par zone thermodynamique, ou un changement de chaîne de séparation dans une zone thermodynamique donnée. Dans ce qui suit, on simplifie d'abord la présentation en ne considérant qu'une chaîne de séparation identique pour tous les éléments de l'ensemble E , puis on considère le cas de plusieurs chaînes de séparation. Les produits de séparation d'une phase huile de composition (x1e .. xNrbe) sont en général une phase liquide (phase huile) et une phase vapeur (phase gaz), dont les masses molaires sont supposées connues et désignées respectivement par MMoo et MMOGe, dans les proportions molaires supposées connues de #Oe de phase vapeur et (1 - #Oe) de phase liquide. #Oe peut prendre toutes les valeurs réelles entre 0. 0 et 1.0, bornes incluses. L'indice O (en caractère gras) rappelle qu'il s'agit des produits de séparation de la phase huile. Les indices 0 et G (en caractère non gras) donnent la qualité huile ou gaz du produit de séparation. e e Les produits de séparation d'une phase gaz de composition y1e yNrbe) sont en général une phase liquide (phase huile) et une phase vapeur (phase gaz), dont les masses molaires sont e e supposées connues et désignées respectivement par MMGOe et MMGGe, dans les proportions e e e molaires supposées connues de #Ge de phase vapeur et (1 - #Ge) de phase liquide. #Ge peut prendre toutes les valeurs réelles entre 0.0 et 1.0, bornes incluses. L'indice G (en caractère gras) rappelle qu'il s'agit des produits de séparation de la phase gaz. La Fig. 1 schématise la séparation distincte d'une phase huile et d'une phase gaz en équilibre dans le réservoir dans la situation générale où la séparation (simple ou complexe) d'un mélange initialement monophasique, liquide (huile) ou gaz, produit deux phases, l'une gazeuse, l'autre liquide. Les relations suivantes traduisent la conservation de la masse lors des opérations de séparation respectivement de l'huile et/ou du gaz de l'élément 'e': (4) MMo = (1 - Oo) MMoo + 6o MMoG MMG = (1 - #Ge) MMGOe + #Ge MMGGe où MMo et MMGsont les masses molaires des phases huile et gaz d'un élément 'e'. d) La méthode de pseudoïsation proposée consiste à décrire chaque phase hydrocarbonée de l'ensemble E en au moins trois constituants (V, I, H). On définit les pseudo-composants en considérant une hypothèse désignée ci-après par M , selon laquelle les phases gaz issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le troisième constituant (H) est exclu, et que les phases huile issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le premier constituant (V) est exclu. Les constituants (V), (I), (H) s'apparentent respectivement à un constituant volatile, intermédiaire et lourd, sans qu'aucun d'entre eux ne corresponde à une sélection particulière de composants ou pseudo-composants de base qui serviraient à une description compositionnelle détaillée des fluides. La modélisation associée est désignée ci-après en abrégé par BRO, pour 'BROWN OIL', pour évocation de la couleur de l'huile qui passe du noir au marron avec un contenu croissant en constituants intermédiaires. En un premier temps, pour faciliter l'exposé, on considère selon une interprétation restrictive de l'hypothèse maîtresse M que les mélanges produits par la chaîne de séparation sont: - un mélange en proportions variables des pseudo constituants volatile (V) et intermédiaire (I) pour un produit gazeux (la masse molaire du constituant V étant inférieure à celle du constituant I) ; - et un mélange en proportions variables des pseudo constituants intermédiaire (I) et lourd (H) pour un produit liquide (la masse molaire du constituant I étant inférieure à celle du constituant H). La restriction sur l'hypothèse M sera levée dans le paragraphe tt). Les relations suivantes s'appliquent donc: ) MMoo = xoi MMI + xoH H XOI + XOH = 1 e MMOG = yov MMV + yOI MMI ; yov + yoi = 1 e MMGO = xGI MMI + XGH MMH ; xGI + xGH = 1 e MMco = yGV IVIMV + yGI MMi ; YGV + YGI = 1 Pour conserver un sens physique à la représentation ternaire, il est nécessaire que: - les masses molaires MMv, MMI, MMH soient définies positives ; e e e e e e - les fractions molaires xoI , xoH , Yov , yoI , xci, xGH yGV yogi soient comprises entre 0. 0 et 1.0, bornes incluses. Il en découle que le choix des masses molaires MMv, MMI, MMH doit respecter les inégalités suivantes: MMV # Min {MMGGe, MMOGe} (eE (eeE) Ge, MMOGe} # MMI # Min {MMGOe, MMOCe Max { MMGG, G 1(eeE) MMI Min {MMGo, MMoo } 1(eeE) MMH ? Max NIM(, 0 e oo )(e6E) (e#E) Dans les équations précédentes, les accolades désignent un ensemble de valeurs sur lesquelles on recherche une valeur minimum (Min) ou maximum (Max). Les inégalités étant non strictes, l'absence de l'un des constituants (V, I) dans un produit de séparation gazeux donné est possible mais ne constitue pas une règle générale et, de même, 5 l'absence de l'un des constituants (I,H) dans un produit de séparation liquide donné est possible mais ne constitue pas une règle générale. e) Les valeurs des masses molaires des constituants (V), (I), (H) doivent être choisies en respectant les conditions ci-dessus et être toutes différentes les unes des autres. Dans l'exemple fourni à titre d'illustration, on a posé, par analogie avec une démarche Black- 10 Oil, et sans que ce soit un exemple limitatif de réalisation: cet e e MMy Min { MMGG, MMoG 1(eeE) ; MMH --- Max IMM, o@ MMoo} 1(eEE) L'exemple d'illustration, largement utilisé dans la suite de la description sous la désignation de cas 'SPE3', correspond au premier cas test de la publication déjà citée de D. E. Kenyon et G. Alda Behie d'exploitation d'un gisement de gaz à condensat par 15 cyclage de gaz. f) Une fois les masses molaires des constituants (V), (I) et (H) définies, la composition des fluides issus de la chaîne de séparation considérée en c) se déduisent, par les hypothèses du paragraphe d) et les équations (5), par : (6) XOV = 0 ; e MMH - MMOO e Mmoo - MMi XOI MM H- MM ; I XOH MM H- MM e MMI - 1,0G e MMOG - MMy Yod Yov = MMi-MMy ' YOI = MMi-MMy XGV = 0 ; e MMH - MMGO e ~ MMGO - MI XGI MM H- MM ; r XGH= MM H- MM e ~ MMi - GG e NIMGC, - MMy GH = 0 YGV = MM MM ; YGI = I ~ MMV ; g) Pour obtenir la composition des phases de l'élément 'e' dans la représentation à trois constituants, on utilise les relations suivantes (dans lesquelles une notation en caractères gras est adoptée pour désigner les fractions molaires des fluides in situ) qui traduisent la conservation molaire lors des opérations de séparation respectivement de l'huile et/ou du gaz de l'élément 'e': (7) XK e = (1 - 80) e XOK e + 80 e YOK e } K - - V ,. 1 H YK = (1 - #Ge) XGK + #Ge yGKe En introduisant les relations (6) dans le jeu d'équations (7), on obtient dans le réservoir pour une phase huile de l'élément 'e' la composition suivante: (8) e e BAMI - IVIMOG Xv = 8 :MM 1- MM V XI = 1 - 80 MMH - MMOO fi MMOG - MMv XI= 1-uo :MM H- :MM 1 +O :MM 1- :MM V e e zizi XH = l - e0 MMH - MMI et pour une phase gaz: (9) e eMMi-MMe Yv = 8G MMi-MMv MMH MAMY e N4MH -IVIMGO,MMGG - MMV YI=1-eG MMH - MMI +0 G MMI - MMV e e GO - I YH = (1 - 0G) MMH-mmi h) A l'issue des étapes a) à g), on dispose donc, pour une utilisation ultérieure, des données minimales pour définir les constituants (V), (I), (H) à savoir leurs masses molaires. On note que, pour l'ensemble des éléments 'e' correspondant à des états d'équilibre entre phases hydrocarbure, on dispose, d'après les jeux d'équations (8) et (9), des coefficients d'équilibre de chacun des constituants (V), (I), (H) : (10) [ Kv e = e = e 8 MMI - MMG Xv 60 MMI - MMoG (1 - e MMH - MMGo 8 MMGG - MMv KI e = Xl e = MMH - MMI +Bg MMi-MMy K& yov~ e ypi MMqg - MMy MMH - MMi e ~YOH Kov KOI = zou - MM MM ; KOH xov XOI H - 00 I V xOH e e e Ke - YGV - . e - Ylli - GG - é H - ~ I . Ke - yen XGV XGI H - GO I V XGH i) Le jeu d'équations (7), écrit pour respectivement les constituants (V) et (H), donne les relations suivantes: (12) Xv e = 80 e Yov e xH = {1- eo) XoH yV = #G yGV ; yH@ = (1 - #G@) xGH@ dont il découle: (13) e XV e XH Yov = 0 0 e XOH = 1 - eô e e e Yv e YH YGV = 8G XGH = 1 - eG puis: (14) e 1 - Xv -e e XH XV e e yo.= 1- XOI=1-1-eô KOI = 1-ee0 e 1 - e e e s Yv YH , yv YGI - 1 - 6c e XGI = 1 - 1-OG e i- 8 é e KGI= e 1 - YH e 1- eG Le jeu d'équations (14) montre que, dès lors que les masses molaires des constituants (V), (I), (H) sont fixées, et que les compositions des phases d'un élément 'e' dans la représentation ternaire ont été établies, les coefficients d'équilibre du constituant (I) peuvent être obtenus à partir de la seule connaissance des fractions molaires des constituants (V) et (H) dans les phases huile et gaz de l'élément 'e' et des fractions vapeur #Oe et #Ge en sortie de chaîne de séparation pour chacune des phases. Une première implication est que, en cas de chaînes de séparation multiples, il est possible de générer les données de coefficients d'équilibre du constituant (I) relative à une autre chaîne de séparation (ci-après désignée par CS~B) que celle utilisée pour obtenir les compositions huile et gaz de l'élément 'e' (ci-après désignée par CS~A). Ceci est schématisé dans le tableau suivant : Mélange Opération Données utilisées Mélange Coefficients détaillé ternaire d'équilibre CS~B - 0 G 0 {xv, xl5 xh}cs~a 1 - XCS-A 1 - e -~~> ~-~ ~~--~~~> ka 01 = 80 e CS (Xi e nti J0!= e CES-A X1 .. Xrby XHCS~A CS~B {#Oe}CS~B # 1 - #O CS~B ~-- - ~> CS A e e e e e e e - {MMGo, MMGG, SG }cS-A {Yv, YI, YH}cS-A Y CS-A # # KGe SGCS-B (YI e .. YNrb) e KGI = 1 - 3'H éS CS~B {#Ge}CS~B # 1 - #G CS~B # # Les conditions nécessaires à une telle procédure sont le respect des inégalités suivantes: e e e e e yo, cs >- 0 et xo, CS~B 0 d' où xv cs 80 cs 1 - XH CS-A e e e e e yci CS-B 0 et xqi CS-B 0 d' où yv cs-A On peut, si besoin est, réviser le choix des masses molaires des constituants (V), (I), (H), en tenant compte des masses molaires des phases issues de la séparation CS~B dans l'écriture des inégalités du paragraphe d). j) A titre d'exemple, la démarche non limitative de réalisation suivie pour obtenir la représentation BRO pour le cas 'SPE3' présenté en illustration a consisté en différentes étapes : - générer avec un simulateur thermodynamique tel que par exemple le simulateur dit 'PVT PACKAGE') un ensemble d'états 'ed' par une opération de dépressurisation à volume constant à la température du réservoir (93.3 C), à partir de la pression de saturation du fluide (237.4 bar) jusqu'à 69 bar, en utilisant une équation d'état de Peng-Robinson, d'emploi très largement répandu dans le contexte pétrolier, pour le fluide détaillé à 16 constituants dont la composition initiale et les paramètres sont donnés dans le Tableau 1 ci-après : Tableau 1: Cas SPE3 - Description de référence du fluide de réservoir Nom Fraction MM Tc Pc Vc (J) 8.* N2 Cl C02 C2 C3 molaire F psi ft3/lbm N2 .0194 28. 016 -232. 33 492. 84 1. 4433 0. 0350 N2 0.000 # # # C1 .6599 16. 043 -115. 76 672.98 1. 5858 0. 0130 C1 0. 120 0. 000 C02 .0121 44.010 88.07 1069. 7 1. 5057 0. 2250 C02 0.000 0. 15000 0.000 # C2 .0869 30. 070 90. 34 709. 67 2. 3707 0. 1050 C2 0. 120 0. 00000 0.150 0.000 C3 .0591 44.097 206. 28 617. 28 3. 2037 0. 1520 C3 0. 120 0. 00000 0. 150 0. 000 0.000 IC4 . 0239 58. 124 274. 96 528. 95 4. 2129 0. 1918 IC4 0. 120 0.00000 0. 150 0.000 0.000 C4.0278 58. 124 305. 64 550. 56 4. 1007 0. 2010 C4 0. 120 0. 00000 0. 150 0.000 0.000 IC5 . 0157 72. 151 370. 11 483. 41 4. 9337 0. 2060 IC5 0. 120 0. 00000 0. 150 0. 000 0. 000 C5 .0112 72. 151 385. 92 489. 36 4. 9817 0. 2520 C5 0. 120 0. 00000 0. 150 0. 000 0.000 C6 . 0181 84. 000 452. 25 468. 47 5. 9268 0. 2809 C6 0. 120 0. 03000 0. 150 0. 010 0.010 C7.0144 100. 205 512. 64 396. 68 6. 9200 0. 3520 C7 0. 120 0. 03000 0. 150 0. 010 0.010 C8 .0150 114. 232 564. 64 362. 02 7.8811 0.3992 C8 0. 120 0. 03200 0. 150 0. 010 0.010 C9 .0105 128. 259 613. 31 331. 99 8. 7621 0. 4439 C9 0. 120 0. 03400 0.150 0.010 0. 010 C10 . 0073 142. 286 655. 07 304. 15 9. 6591 0. 4869 C10 0. 120 0. 03600 0. 150 0. 010 0.010 C11 . 0049 147.000 699. 91 349. 54 10. 172 0.4770 C11 0.120 0. 04300 0. 150 0. 010 0.010 C12P.0138 210. 000 797. 52 324. 79 14.000 0.5266 C12P 0. 120 0. 04408 0. 150 0. 010 0. 010 Dans ce tableau, les paramètres d'interactions binaires forment une matrice symétrique de trace nulle ; les cases où figurent un symbole # sont remplies par symétrie. Les colonnes manquantes doivent être remplies d'une part par la propriété de symétrie, d'autre part par des paramètres d'interactions binaires nuls. - A divers paliers de dépressurisation répartis entre la pression de saturation du fluide initial et 69 bar, on a effectué une opération de vaporisation à pression constante (la pression du palier considéré) de la phase huile de l'état 'ed' considéré par un gaz dont la composition est tenue pour représentative du gaz d'injection en cours de simulation (dans le cas SPE3, la composition du gaz injecté varie en cours de simulation puisque c'est le gaz produit qui est réinjecté), chaque opération de vaporisation comportant plusieurs paliers de quantités de gaz injecté pour aller jusqu'à la vaporisation complète de la phase liquide. A chaque palier de quantité de gaz injecté où les phases sont saturées correspond un état 'edv'. - Toujours avec le même simulateur, les opérations de séparation des phases huile et gaz de chacun des états 'ed' et 'edv' formant l'ensemble E ont été simulées, la chaîne de séparation CS~A utilisée comportant un premier étage à 26.7 C et 56. 2 bar, un second étage de détente de la phase liquide issue du premier étage de séparation à 26.7 C et 4. 5 bar, une dernière détente de la phase liquide issue du deuxième étage de séparation à 15.6 C et 1. 01325 bar, les phases gazeuses issues des différentes séparations étant mélangées. Comme le cas SPE3 présente une modification, au cours de la simulation, des conditions de pression du premier étage de séparation (la pression passant de 56.2 bar à 21. 7 bar), les opérations de séparation des phases huile et gaz de chacun des états 'e' ont également été simulées avec cette autre chaîne de séparation CS~B. Il convient de noter que l'enchaînement des diverses opérations de simulation ci-dessus est aisément réalisable avec la plupart des logiciels de simulations thermodynamiques disponibles pour l'industrie. Les conditions nécessaires décrites dans le paragraphe d) trouvées à partir des masses molaires des produits de séparation des phases huile et gaz de chacun des états de l'ensemble E par les chaînes CS~A et CS~B sont en grammes/mole: MMV # 21. 592 ; 27.865 # MMI # 72. 334 ; MMH # 162. 06 Les masses molaires attribuées aux trois constituants dans la représentation BRO sont données dans le Tableau 2. Les compositions à trois constituants des phases de l'ensemble E par les calculs décrits dans les paragraphes f) et g) ont été obtenues en considérant les produits de séparation de la chaîne CS~A. La composition du fluide initial dans la représentation BRO est donnée dans le Tableau 2 ci-après : Tableau 2 : Cas SPE3 - Description BRO du fluide de réservoir Nom du MM Fraction constituant molaire V 21.592 .856714 I 62. 334 .092358 H 162.06 .050928 Propriétés des phases La règle de Gibbs, qui donne le nombre de degrés de liberté d'un système thermodynamique à partir du nombre de constituants et du nombre de phases, appliquée à un mélange à trois constituants, donne trois degrés de liberté lorsque deux phases hydrocarbure sont à l'équilibre : les propriétés de phases dépendent de la pression, la température et d'une variable compositionnelle. Pour chaque phase fluide pour laquelle on connaît, pour des conditions de pression et température données, les valeurs de propriétés de phase utiles pour la modélisation des écoulements (facteur de compressibilité, viscosité,..), on peut rentrer ces propriétés de phase sous forme de tables fonctions de la pression, de la température, et d'un index de la composition dans la représentation à trois constituants. Ces données peuvent être utilisées ensuite comme données d'entrées dans un simulateur, le simulateur disposant en général de méthodes internes d'interpolation ou d'extrapolation pour estimer les propriétés de phases en des points intermédiaires. k) Alternativement à l'usage de tables, on peut rechercher, par exemple par régression, les paramètres de corrélations ou d'équations d'état permettant de calculer, par zone thermodynamique, les propriétés des phases en fonction de la composition et des conditions de pression et de température, ces paramètres étant ensuite introduits dans le modèle de simulation en lieu et place des tables précédemment citées. Une situation fréquente est de disposer d'une équation d'état fiable pour le fluide de réservoir décrit de façon détaillée dans la représentation de base. Ainsi, lorsqu'une équation d'état (et une corrélation associée pour le calcul des viscosités) définie pour la représentation détaillée à Nrb constituants, permet de modéliser le comportement thermodynamique des fluides au cours des étapes de pressions et températures d'un chemin thermodynamique suivi dans une ou plusieurs zones thermodynamiques, on peut utiliser cette équation d'état pour générer l'ensemble de référence de compositions de fluides et les propriétés déjà citées des phases correspondantes, en y adjoignant les paramètres des phases dans l'équation d'état, tels que le co-volume, le terme d'attraction, la pression critique,.. qui peuvent être considérés comme des propriétés spécifiques des phases. Pour chacune des compositions de fluides obtenues, cette même équation d'état, ou une autre, peut être utilisée pour générer les compositions et propriétés des mélanges fluides produits par une chaîne de séparation, et permettre ainsi l'application de la méthode de pseudoïsation décrite dans les paragraphes a) à i). Lorsque les paramètres des pseudo-composants dans une équation d'état, et dans une corrélation pour le calcul de viscosité, sont ajustés, par exemple par régression, pour reproduire les propriétés spécifiques des phases telles que le co-volume, le terme d'attraction, la température critique,.. obtenues avec la représentation détaillée, on obtient une équation d'état et une corrélation utilisable pour les calculs de propriétés de phase telles que le facteur de compressibilité, la viscosité, mais cette première équation d'état, désignée ci-après en abrégé par EOS~PRO, n'est pas forcément valide pour les calculs d'équilibre notamment si le respect d'une condition d'égalité des fugacités à l'équilibre n'a pas été pris en compte. Une démarche non limitative de réalisation pour obtenir les paramètres de l'équation d'état EOS~PRO de la représentation BRO est celle suivie pour le cas 'SPE3' présenté en illustration. L'équation d'état utilisée, celle de Peng-Robinson, est rappelée ci-après pour aider à la compréhension de la démarche, mais la démarche est aisément transposable à d'autres équations d'état. 1) L'équation de Peng-Robinson s'écrit pour la phase P : (15) p = R T - a(Tp (15) P=-##-###"t###? (15) Vp - bp (Vp+bp) -2bp où R est la constante des gaz parfaits, T la température, P la pression, V le volume molaire de la phase P, b son co-volume et a (T) terme d'attraction, ces dernières quantités étant exprimées par : (16) bP = #ci bi i (17) a(T)p = c [ t c aj(T)] avec a;(T) = ai(T) a(T) (1 - Si,) où ci désigne la fraction molaire du constituant i dans la phase P, fraction molaire dont les dépendances physiques sont omises pour alléger la présentation, bi et a;(T) respectivement le co-volume et le terme d'attraction du constituant i, #ij le paramètre d'interaction binaire entre les constituants i et j, les #ij formant une matrice symétrique de trace nulle. Sous forme adimensionnelle, l'équation de Peng-Robinson s'exprime comme: 3 3 (18) Zp - (1 - Bp) Zp +(Ap-3Bp -2Bp)Zp-(ApBp-Bp -Bp)=0 où Z est le facteur de compressibilité de la phase P, et B et A sont les formes adimensionnelles des co-volume et terme d'attraction de la phase P (les dépendances des variables avec la pression et température étant omises pour alléger l'écriture des équations): (19) Bp = b -p (20) Ap = ap-YZ R T m) Le facteur de compressibilité de la phase P est uniquement fonction des deux paramètres de phase adimensionnels A et B. Il s'ensuit qu'estimer les paramètres par constituant dans la représentation BRO (à savoir bv, bi, bH, av(T), ai(T), aH(T), #VI, #VH, #IH) permettant de reproduire les valeurs de a (T) b, obtenues pour chaque phase aux différents paliers des simulations PVT avec la représentation de base (détaillée à 16 constituants), permet de reproduire les valeurs de référence des facteurs de compressibilité des phases. La procédure de régression utilisée pour obtenir les différents paramètres de la représentation BRO est la minimisation de fonctions objectif de la forme: # e#E e # e 2 e#E e e # e (21) O(1) = e c [ Le minimum de O(p) est obtenu par la résolution du système d'équations, de même # # dimension que le vecteur p, formé par la nullité des dérivées partielles de la fonction O(p) par rapport aux différents paramètres p#: (22) #p# O(#) = 0 soit : (23) 1 % [ (24) 0(bv,bj,bH) ec- E M('; e E YK e bK - bcl 2 + o 1 x K e bK - bo K e K Le système à résoudre est simplement un système linéaire de trois équations à trois inconnues: (25) mvv bv + mvi bi + mVH bH = rv mIV bv + mu bi + mIH bH = rI # mHV bV + mHI bI + mHH bH = rH les termes mKK' et rK des seconds membres des équations étant donnés par : e#E e e e#E e e e e mKK' eue (mG YK' yK'e + M 0 e XKe xK0 ; rK eue ( û50 xKe bo + ûfG ege bG) e e o) Pour rechercher les paramètres av(T), ai(T), aH(T), 8vi, #VH, #IH par simplement des résolutions de systèmes linéaires de trois équations à trois inconnues, l'équation (17) peut être réécrite comme la somme de deux termes, le second terme comprenant les interactions binaires : (26) a(T)p = a1(T)P 2 + a2(T)p avec : (27) a1(T)P = # ci #ai(T) et a2(T)p = - Cj Ci cj où #j(T) = a(T) a(T) 8 j 1 Dans le terme a2(T)p, les paramètres Oij(T) forment un ensemble de trois paramètres indépendants [alpha]VI(T), [alpha]VH(T), am(T). En utilisant les termes al(T), obtenus lors des simulations PVT avec la représentation détaillée à 16 constituants, comme quantités observables o par phase et par état e, les inconnues étant les quantités #\/av(T), -\/aï(T), aH-(T), le problème à résoudre est bien analogue à celui décrit pour l'obtention des paramètres de co-volumes des constituants de la représentation BRO. De même, les paramètres [alpha]ij(T) sont obtenus en résolvant le système de trois équations à trois inconnues formé à partir des observables a2(T). p) Les paramètres d'interactions binaires peuvent être ensuite obtenus par: (28) - avi(T) 8VH #YH(T) 8IH = - #rn(T) av(T) 7a,(T) YH - ay(T) aH(T) ffl aI(T) aH(T) q) Certains simulateurs acceptent, comme paramètres des constituants pour l'équation d'état, l'introduction directe des co-volumes bi des constituants et de la matrice des paramètres aij(T) intervenant dans le calcul du terme d'attraction (cf. l'équation (17)), cette matrice devant être définie le cas échéant pour différentes températures. Au cas où le cas traité nécessiterait une dépendance en température (température variable dans le réservoir), la démarche précédente pourrait être adaptée, par exemple en générant l'ensemble E par des opérations PVT à diverses températures de réservoir, en considérant des états e à diverses températures pour obtenir les co-volumes par constituant (ceux-ci devant être les mêmes pour les différentes températures), en considérant les états e obtenus à une température donnée pour générer la matrice des paramètres aij(Te) correspondante et en répétant la démarche pour autant de températures Te que nécessaire. r) Lorsque le simulateur n'accepte pas l'introduction directe du co-volume par constituant b, et des matrices de paramètres aij(Te): - le co-volume par constituant est en général obtenu à partir de l'introduction des paramètres Tci, Pc;, optionnellement #bi (par défaut tous égaux à 0.077796) par : (29) bi = R #b1Tc1/Pci - les quantités ai(T) sont en général obtenues par: .TFM dp (30) a;(T) = a0; [ 1 + Si ( 1 - VkJJ avec i = C. ffi" n7--O La quantité a0i est en général obtenue à partir de l'introduction des paramètres Tci, Pc1, optionnellement #ai (par défaut tous égaux à 0. 457235) par: (31) aOi R2 SZa; pc (31) ni Pci Dans le polynôme #i de degré dp, les paramètres #1 sont les facteurs acentriques des constituants en général considérés comme constants, le degré du polynôme et les constantes Cn pouvant varier d'un simulateur à l'autre. s) Lorsque les pressions, températures et volumes critiques par constituant sont requis, on peut utiliser comme observables les pressions, températures et volumes critiques des phases de l'ensemble E , définies par des lois de mélanges identiques à celles utilisées pour les co-volumes : (32) Pcp = LCi PCi ; Tcp = LCi TCi ; VcP = ici ici i i i L'obtention respectivement des paramètres Pcv, Pci, Peu, Tcv, Tci, TcH et Vcv, Vci, VcH peut donc s'effectuer par résolution successive de trois systèmes de trois équations à trois inconnues totalement analogues à ceux générés pour le calcul des co-volumes par constituant, développé précédemment. t) Les paramètres #bk (K=V, I,H) sont alors obtenus, à partir de l'équation (29) et des paramètres bK, PcK et TcK tels que déterminés dans les paragraphes n) et s), par: r PACK BOK (33) bK = R TcK u) Dans l'équation (30), dans le cas le plus fréquent où la température de réservoir est considérée comme constante, on peut s'affranchir de la dépendance en température des quantités #aV(T), #aI(T), #aH(T), en adoptant une valeur nulle pour les facteurs acentriques des constituants V, I, H, d'où: (34) #aK(T) #a0K, VK, K=V, I,H et on obtient alors les paramètres #aK (K=V,I,H) , à partir des paramètres aK(T), PcK et TcK tels que déterminés dans les paragraphes o) et s), par: PCK aK(T) (35) .QaK = R2 Tel Une autre possibilité, utilisée dans le cas SPE3, est l'attribution aux paramètres #aK (K=V,I,H) de la valeur par défaut. Dans ce cas, les paramètres #a0K sont déterminés avec les PcK et TcK tels que déterminés dans le paragraphe s), et les paramètres K sont ensuite obtenus par réarrangement de l'équation (30): (3 ) &=### (36) SK= 1 - TCK T 1 Tcx puis, les facteurs acentriques #K ont été recherchés comme solutions des trois équations polynomiales (K=V,I,H): dp (37) # Cn #Kn - #K = 0 n=0 v) Lorsque pour le cas considéré, la température du réservoir est variable, on peut réécrire l'équation (30) comme : (38) #aK(Te)= ssK + [gamma]K #K(Te) avec les quantités : #k(Te) = 1 - #Te/TcK, calculées en utilisant les températures critiques des constituants telles que déterminées dans le paragraphe s). On peut ensuite rechercher les quantités ssK et [gamma]K (ssK = #a0K et [gamma]K = #a0K K) comme des paramètres indépendants par régression linéaire des couples de valeurs #K(Te), #aK(Te), ces dernières valeurs étant telles que déterminées dans les paragraphes o) et q). Les paramètres #aK et K sont alors obtenus par: (39) K et zak Les facteurs acentriques recherchés sont solutions des trois équations polynomiales (K=V,I,H): dp (40) Cn COKn ~ - 0 n=0 K w) Des facteurs de correction de volume sont souvent utilisés pour ajuster les densités des phases. Lorsque des facteurs de correction de volume par constituant sont requis, on peut utiliser comme observables les corrections de volume des phases de l'ensemble E , définies par des lois de mélanges identiques à celles utilisées pour les co-volumes: (41) CvP = #c1 Cvi i L'obtention des paramètres Cvv, CvI, CvH, peut s'effectuer alors par résolution d'un système de trois équations à trois inconnues totalement analogue à celui généré pour le calcul des co-volumes par constituant développé dans le paragraphe n). x) Dans le cas SPE3, les calculs de viscosité des phases hydrocarbure ont été effectués avec une corrélation largement utilisée dans le domaine de l'ingénierie pétrolière, celle de Lohrenz, John, Bray, B. G. et Clark, C.R.: Calculating Viscosities of Réservoir Fluids from Their Compositions , Journal of Petroleum Technology, 1964, pp. 1171-1176. Les paramètres de phases utilisés explicitement dans cette corrélation sont les pressions, températures, volumes critiques des phases. La détermination des paramètres critiques (pression, température, volume) des constituants V, I, H, telle qu'exposée dans le paragraphe s), a permis une modélisation appropriée des viscosités des phases avec la représentation BRO. La démarche est aisément transposable à d'autres corrélations. y) La ligne directrice suivie dans l'ensemble des paragraphes n) à w) est de rechercher les paramètres d'équation d'état EOS~PRO utiles aux calculs, dans les conditions de réservoir, de densités des phases (via les facteurs de compressibilité) et de viscosités par des résolutions simples (non itératives) de systèmes d'équations linéaires. Son efficacité pour reproduire les co-volumes et termes d'attraction des phases est illustrée, pour les fluides du cas SPE3, par les Fig. 6 et 7 pour une opération de dépressurisation à volume constant à la température de réservoir (93.33 C) à partir de la pression de rosée bars du fluide de réservoir initial (237. 4 bars) jusqu'à 69 bars. Sur chacune des figures, sont illustrées les quantités calculées à partir de l'équation d'état de la représentation détaillée de référence et les quantités calculées à partir de l'équation d'état EOS~PRO de la modélisation BRO, ces quantités étant pratiquement confondues. Cette même ligne a été utilisée avec le même succès pour traiter d'autres cas non présentés ici. On a pu constater qu'une telle ligne est largement adaptable à des contextes de température variable. Par ailleurs, les traitements nécessaires ne représentent aucune difficulté particulière de codage et ont été aisément automatisables. Les paramètres de l'EOS~PRO utilisés pour le cas SPE3 en illustration sont donnés dans le tableau 3 ci-dessous : Tableau 3 : Cas SPE3 - Paramètres de l'EOS~PRO de la représentation BRO Nom T P V w S2b Cv 8. v I l I H F psi ft3/lbm ft3/lbm V -61.212 666. 02 1. 869 -0. 02043 0. 079345 -0. 0686 V 0.0 0. 058447 0. 002065 1 368.292 581.21 4. 389 0. 31255 0. 082877 -0. 0621 i 0. 000000 0. 001604 H 663.138 263. 46 10. 99 0. 21233 0. 065003 0. 2664 H # 0. 000000 z) Dans le cas SPE3, on a choisi d'utiliser pour l'EOS~PRO une équation d'état du même type que celle utilisée pour la représentation de base. Ceci n'est pas obligatoire. Il est possible d'utiliser une équation d'état d'un autre type, dès lors que cette autre équation d'état utilise des paramètres de phase identifiables aux paramètres de phase utilisées par l'équation d'état de référence. Plus généralement encore, on peut obtenir les paramètres nécessaires à la représentation BRO par diverses techniques dites de régression, inversion, optimisation,.. que l'on peut trouver dans une littérature abondante. On peut par exemple les rechercher par une ou plusieurs résolutions itératives de systèmes d'équations formés par des dérivées partielles, par rapport à des sous-ensembles ou à l'ensemble des paramètres à la représentation BRO, d'une fonction objectif du type de l'équation (21) ou d'un autre type. On pourrait même chercher à utiliser un formalisme de ce type pour optimiser l'une ou plusieurs des masses molaires des constituants V, I, H, en la (les) faisant apparaître explicitement, - par remplacement dans la fonction objectif des compositions des phases par leurs expressions données par les équations (8) et (9)-, et en choisissant des observables sensibles à leur(s) valeur(s). Un telle démarche pourrait être envisageable si l'on se trouvait dans une situation, non rencontrée jusqu'à présent, où la ligne directrice suivie dans l'ensemble des paragraphes n) à w) conduisait à une détermination de paramètres de la représentation BRO ne permettant pas une reproduction fidèle des observables. Enfin, pour des raisons connues en soi, on peut être amené à donner a priori des valeurs à certains des paramètres de l'équation EOS~PRO et n'effectuer alors des calculs d'estimation que pour un nombre réduit de paramètres. aa) La première expression généralement attendue des résultats de simulation sont des résultats en volume. Les observables que l'on doit alors chercher à restituer sont les volumes molaires des phases en sortie de séparation. Pour le cas SPE3, on a utilisé deux corrélations distinctes, l'une pour le calcul des volumes molaires des phases gaz en sortie de séparation, l'autre pour le calcul molaire des phases huiles en sortie de séparation, les conditions de sortie de séparation étant les conditions standard usuelles de la profession à savoir une atmosphère et 15.56 C (60 F). De nombreuses corrélations ou équations d'état sont disponibles pour estimer le volume molaire d'une phase gaz notamment dans des conditions proches des conditions atmosphériques. A titre indicatif, c'est la corrélation de P. M. Dranchuk et J. H. AbouKassem (présentée dans le papier : Calculating Z Factors For Natural Gases Using Equations of State , JCPT, Juillet-Septembre 1975) qui a été utilisée pour le cas SPE3 où les paramètres utiles au calcul du volume molaire du gaz sont les pressions et températures critiques des constituants (V) et (I), paramètres qui ont été régressés de façon itérative de façon à reproduire les volumes molaires des phases gaz obtenus par les simulations PVT de séparation avec la représentation à 16 constituants. La variation, en fonction de la fraction molaire du constituant I, du volume molaire des phases huile en des conditions proches des conditions atmosphériques a été trouvée, dans le cas SPE3, pratiquement linéaire, plus généralement, pour des cas non présentés ici, représentable par un polynôme de degré peu élevé (degré deux à trois). Pour le cas SPE3, c'est une dépendance linéaire qui a été adoptée dont les paramètres ont été obtenus par simple régression linéaire. Les paramètres nécessaires aux calculs des volumes molaires en conditions de séparation ont été obtenus pour les deux chaînes de séparation CS~A et CS~B avec la même méthode. Equilibres thermodynamiques bb) Les calculs d'équilibre peuvent être effectués soit à partir de données de coefficients d'équilibre, soit à partir de l'écriture de l'égalité des fugacités de tous les constituants dans chaque phase, ces deux possibilités étant examinées respectivement ce qui suit. cc) Selon la règle de Gibbs, pour le calcul des équilibres dans le réservoir, les coefficients d'équilibre des pseudo-composants (V), (I), (H) tels que calculés par le jeu d'équations (10) peuvent être introduits dans des tables comme fonctions de la pression, de la température, et d'un index de la composition dans la pseudo représentation. Ces données peuvent être utilisées ensuite comme données d'entrées dans un simulateur, le simulateur disposant en général de méthodes internes d'interpolation ou d'extrapolation pour estimer les coefficients d'équilibre en des points intermédiaires. On peut aussi, avant introduction des données dans un simulateur, ajouter des points aux tables en utilisant des méthodes, connues en soi, d'interpolation ou d'extrapolation ou utiliser, à la place des tables, des corrélations bâties pour restituer les coefficients d'équilibre. Une dépendance des coefficients d'équilibre avec, tout à la fois la pression, la température et l'index de composition, n'est pas toujours nécessaire. Ainsi, dans nombre de cas d'application, la température du réservoir est considérée comme constante. Dans certains cas, on peut être amené à abandonner la dépendance par rapport à un index compositionnel. Ces cas où l'on réduirait le nombre des degrés de liberté des coefficients d'équilibre restent dans le domaine de l'utilisation de la méthode proposée. dd) Pour le calcul des équilibres à l'issue d'une chaîne de séparation, les coefficients d'équilibre des pseudo-composants (V) et (H) sont respectivement infini et nul quelques soient les conditions de pression et de température en sortie de chaîne de séparation et quelque soit la composition du mélange soumis à la séparation. Si besoin est, dans un simulateur qui n'accepterait pas le concept de coefficients d'équilibre infini ou nul (ce concept étant souvent implémenté en forçant l'absence du constituant respectivement dans la phase huile et dans la phase gaz), on peut utiliser une valeur finie très grande et une valeur finie très faible, toutes deux positives, en lieu et place de coefficients d'équilibre infini ou nul. Pour les données relatives aux coefficients d'équilibre du constituant (I) nécessaires à un calcul d'équilibre, le jeu d'équations (11) indique une définition différente des coefficients d'équilibre du constituant (I) selon que l'on considère la séparation d'une phase 'gaz' issue du réservoir, ou la séparation d'une phase 'huile' issue du réservoir. Les entrées, en principe nécessaires aux calculs d'équilibre, sont donc, par phase, une table de coefficients d'équilibre du constituant (I) en fonction d'un index compositionnel de la phase soumise à la séparation, de la pression et de la température en sortie de la chaîne de séparation considérée. Les plages de variation des compositions des phases huile et gaz étant en général distinctes, il peut être pratique de 'fusionner' en une seule table les données utiles pour les phases huile et gaz, notamment si le simulateur n'accepte pas l'introduction de données distinctes pour chaque phase. Ces données sont utilisables comme données d'entrées dans un simulateur, le simulateur disposant en général de méthodes internes d'interpolation ou d'extrapolation pour estimer les coefficients d'équilibre en des points intermédiaires. On peut aussi, avant introduction des données dans un simulateur, ajouter des points aux tables en utilisant des méthodes, connues en soi, d'interpolation ou d'extrapolation ou utiliser, à la place des tables, des corrélations bâties pour restituer les coefficients d'équilibre. ee) Si l'équation d'état EOS~PRO est bâtie pour restituer les propriétés des phases, elle n'est pas forcément valide pour les calculs d'équilibre, c'est-à-dire pour permettre de restituer les compositions des phases de l'ensemble E obtenues dans le paragraphe g). Comme précédemment, on utilise ci-après l'équation d'état de Peng-Robinson comme support de présentation. Pour cette équation d'état, les coefficients d'équilibre, d'après la condition d'égalité des fugacités de chaque constituant dans chaque phase à l'équilibre, s'obtiennent par : N (42) Ln(Ki) = Do + D1 b; + ~ r R2 T2 , . C B y - ô xJ al(T) j=l avec les paramètres adimensionnels de phase définis selon les équations (19) et (20) et: (43) + ' Bp (43) =Zp+(l-2)Bpj Do=J1 (ZO - Bo) Do=Ln Zo-Bo P Zp-1 1 ZG - 1 Ao lo AG LG Dl = R T Bo BG +2~ V r2Bo2 2 -J2 Bo 2 et les facteurs de compressibilité des phases calculés à partir de l'équation (18). Pour traiter un cas où la température de réservoir est constante, le problème est de déterminer les neuf paramètres indépendants formant un vecteur p: by, bI, bH, aVV, an, aHH, aVI, aVH, aIH, par exemple par minimisation d'une fonction objectif construite en utilisant comme observables OKles coefficients d'équilibre de l'ensemble E tels que calculés par les équations (10). La fonction objectif selon l'équation (21) peut alors être adaptée comme : # e#E e e # e (44) O(#) = # # #K[#K(#) - oK] e K=V,I,H où les fonctions #K(#) sont les exponentielles des seconds membres des équations (42). 3 Alternativement, dans l'équation (44), on pourrait utiliser comme observables oK les e # logarithmes des coefficients d'équilibre de l'ensemble E et comme fonctions #K(p) sont les seconds membres des équations (42). En suivant la démarche de K. E. Starling dans le papier: A New Approach for Determining Equation-of-State Parameters Using Phase Equilibria Data , SPE 1481, e SPEJ, Décembre 1966, on peut choisir les poids @ comme les inverses des valeurs des observables oK. La nullité des dérivées partielles de la fonction O(p) par rapport aux neuf paramètres d'intérêt fournit un système de neuf équations non linéaires à neuf inconnues, résoluble uniquement de façon itérative, de nombreuses méthodes adaptées à la résolution des systèmes non linéaires étant discutées dans la littérature. Un calcul itératif de ce type nécessite une première approximation du vecteur de paramètres. On a pu constater, notamment pour le cas SPE3, qu'il est opportun de tirer profit des paramètres obtenus pour l'EOS~PRO pour démarrer les calculs itératifs. Plus généralement, notamment pour reproduire des coefficients d'équilibre à diverses température, on peut considérer un autre ensemble ou sous-ensemble de paramètres pour le vecteur p et d'autres techniques de recherche de paramètres. ff) A titre d'illustration, pour le cas SPE3, ce sont les six paramètres Peu et #bK qui ont constitué le vecteur p, les valeurs des autres paramètres étant reprises de l'équation EOS~PRO. La procédure utilisée dans le cas SPE3 est une procédure connue en soi en cours d'amélioration. L'ensemble des paramètres de l'EOS~EQ utilisés pour le cas SPE3 est donné dans le tableau 4 suivant : Tableau 4 : Cas SPE3 - Paramètres de l'EOS~EQ de la représentation BRO Nom Tc Pc E #b #KK' V I H F psi V -61.212 848.47 -0.02043 0. 089345 V 0.0 0. 058447 0. 002065 1 368.292 274.23 0.31255 0. 070000 # # 0.000000 0. 001604 H 663.138 309.90 0.21233 0. 065003 H 0. 000000 gg) Dans le cas SPE3, les données introduites pour les calculs d'équilibre en surface pour chacune des chaînes de séparation CS~A et CS~B sont les coefficients d'équilibre du constituant (I) sous forme d'une table unique pour la séparation des phases huile et gaz, table obtenue comme décrit dans le paragraphe dd) avec pour seul abscisse l'index de la composition de la phase avant séparation. hh) Les Fig. 2 à 5 comparent les résultats, sur les quinze années d'exploitation simulées du cas SPE3, d'une simulation compositionnelle à trois pseudo-composants BRO (résultats désignés par 'ternaire'), où la représentation compositionnelle a été obtenue par application de la méthode selon l'invention, à ceux d'une simulation compositionnelle détaillée à seize constituants de base (résultats désignés par 'référence'). La simulation compositionnelle de référence utilise une même équation d'état de PengRobinson pour les calculs d'équilibre et pour les calculs de propriétés de phase et pour toutes les zones thermodynamiques. La simulation BRO à trois constituants utilise, pour les calculs dans les mailles de réservoir, l'équation d'état (EOS~EQ) pour les calculs d'équilibre et l'équation d'état (EOS~PRO) pour les calculs de propriétés de phases. Pour les calculs en conditions de surface, la simulation BRO utilise des coefficients d'équilibre introduits sous forme de tables et des corrélations simples pour le calcul des volumes des phases. Les résultats présentés Fig. 2, qui représentent l'évolution de pression en fond de puits producteur, et Fig. 5, qui représentent l'évolution de saturation en huile (ou condensat) dans la maille perforée la plus profonde à l'endroit du producteur, sont parmi les résultats les plus révélateurs du bon calcul des équilibres et des propriétés de phases dans les conditions in situ. On peut constater que la modélisation BRO fournit une solution en pression pratiquement confondue avec la solution de référence (à 16 constituants) et une solution d'évolution de saturation en huile en bon accord avec la solution de référence sur tout le temps simulé. La maille utilisée en Fig. 5 est la même que celle utilisée dans le papier de référence de D.E. Kenyon et G. Alda Behie et dans le papier de W. H. Goldthorpe cité dans l'état de la technique. On peut constater que les résultats de la Fig. 5, jugés par comparaison avec les résultats obtenus par W. H. Goldthorpe avec une modélisation BO, montrent une bien meilleure qualité de reproduction de l'évolution de la saturation en huile dans le réservoir alors même que la procédure utilisée pour l'obtention des paramètres de l'EOS~EQ est améliorable. ii) Les résultats présentés Fig. 3 et 4 de débit d'huile et de production cumulée d'huile sur toute la durée du temps simulé montre que la modélisation BRO fournit une solution très proche de la solution de référence. La réduction de temps de calcul par la modélisation BRO est considérable, le facteur d'accélération étant de 46. Delumping jj) On représente tout d'abord le réservoir sous la forme d'un réseau de mailles (m) dont chacune constitue un volume élémentaire rempli de fluide(s) sous forme de une ou plusieurs phases, avec au moins une phase non aqueuse. Les phases non aqueuses restent dites phases hydrocarbure même si elles peuvent contenir certains composants autres que des hydrocarbures tel que l'azote, le dioxyde de carbone, l'anhydride sulfureux. kk) On définit, pour chaque zone ou domaine thermodynamique, le fluide par une représentation détaillée à Nrb composants et/ou pseudo-composants. Il est à noter que l'on peut traiter les cas qui requièrent l'usage de plusieurs représentations thermodynamiques, par exemple si l'on peut distinguer plusieurs chemins thermodynamiques locaux durant la modélisation (on peut être amené à distinguer des zones produites uniquement par dépressurisation et des zones soumises à une injection de gaz). Plusieurs zones ou domaines de variations des grandeurs thermodynamiques ou de composition, désignés souvent par les spécialistes comme zones thermodynamiques, pourront donc être définis et utilisés. 11) Par zone thermodynamique pour laquelle on adopte une représentation groupée des fluides, on détermine une équation d'état EOS~PRO construite, préalablement à la simulation dynamique de réservoir avec la représentation groupée, pour reproduire les paramètres de phases, dans l'équation d'état de la représentation détaillée, des fluides hydrocarbure au cours de chemins thermodynamiques tenus pour représentatifs de ceux que vont suivre les fluides hydrocarbure pendant la simulation maillée. mm) On réalise, de façon connue en soi, une simulation compositionnelle avec un nombre limité de constituants où les propriétés des phases sont calculées par une équation d'état EOS~PRO, ladite simulation permettant de calculer au moins dans chaque maille (m) et à des pas de temps consécutifs (t, t+1, etc. ), une pression dans une phase hydrocarbure (pm), la température (Tm) (si celle-ci varie), les débits molaires par phase d'injection ou de production (Sgm) et (Som), pour chaque paire de mailles (m, h), les débits molaires des phases liquide (uomht) et vapeur (ugmht), la fraction vapeur (6m). On mémorise ces quantités, ainsi que les paramètres dans l'équation d'état EOS~PRO pour les phases hydrocarbure liquide et vapeur (typiquement les co-volumes et termes d'attraction) intervenant dans l'expression des coefficients d'équilibre. Pour éviter de les recalculer, on peut mémoriser également les facteurs de compressibilité des phases (Zom) et (ZGm) obtenus par résolution de l'équation d'état EOS~PRO. Pour une équation d'état de Peng-Robinson, en conservant les notations introduites par les équations (16), (17), (26) et (27), les paramètres de phase sont désignés par (bom) et (bGm) pour les co-volumes des phases, (aom) et (acm), - ou (alors), (algra), foom). fecm) -, pour les termes d'attraction des phases. Les différents paramètres utiles au calcul des coefficients d'équilibre, lorsque évoqués d'une façon générale, seront désignés ci-après par (E#Pmt), l'indice P étant pour chaque phase huile et gaz, et l'indice n un numéro de référence pour un paramètre particulier. nn) On estime au pas de temps t+1 la fraction molaire de chaque constituant i dans la composition détaillée globale (zim) du fluide hydrocarbure dans la maille (m) en connaissant les nombres de moles de chaque constituant i dans la représentation détaillée dans chaque phase au pas de temps t, respectivement (Moimt) pour la phase huile et Ngimt) pour la phase gaz à partir des équations ci-après: Nrb Nrb (45) NOm t = Noim t Ngm t Ngi m t Nm Nom t + Ngm i i t NOm t N im t NOim + Ngim (46) Xim NOM t im = Ngm Zim = Nm t t t t t t t Zim Nm - Ot (Yim' Sgm + xim Som) - dt 1 (Yim' ügmh + Xim' Uomh) t+1 hE J(m) (47) Z im = ~Fi Nm t+I (48) N m = Nm - Dt (Sgm + Som) - Ot L (Ugmh + Uomh) hE J(m) Les équations (45) et (46) permettent respectivement de calculer, au pas de temps (t) et dans la maille (m), les nombres de moles par phase hydrocarbure et global Nm à partir de la connaissance des nombres de moles par constituant et par phase, et les fractions molaires des constituants par phase et globales. L'équation (48) traduit le bilan molaire total hydrocarbure sur la maille (m), en tenant compte des échanges de matière, pendant l'intervalle de temps #t, avec toutes les mailles (h) voisines de (m) qui forment l'ensemble J(m). Dans l'équation (47), l'écriture des termes (yij'm) et xij'm), - dans lesquels m'=m pour un débit de la maille (m) vers la maille (h) ou dans le puits, et m'=h pour un débit de la maille (h) vers la maille (m), et m' correspondant au fluide injecté dans le cas des puits d'injection, S étant alors négatif-, suppose implicitement l'utilisation d'un schéma simple amont pour les flux compositionnels. Une écriture plus générale de ces termes est (yimh) et (ximht) où ximht et yimh décrivent les compositions des phases liquide et gaz, obtenues d'une façon connue en soi, s'écoulant entre les mailles (m) et (h). oo) On détermine, au pas de temps (t+1) et dans chaque maille (m), les coefficients d'équilibre (Kim) du constituant i dans la représentation détaillée en utilisant une formulation où apparaissent explicitement les paramètres E#P, en y introduisant les paramètres variables (E#P mt+1) issus de la simulation avec la représentation groupée, et, le cas échéant, en utilisant certaines des variables compositionnelles décrites dans le paragraphe précédent calculées au pas de temps (t). Ainsi, pour une équation d'état de Peng-Robinson, les coefficients d'équilibre donnés par l'équation (42) peuvent être réécrits sous la forme: ~~~ t+1 (49) t+1 = Hl + d/*1 + -FI t+1 + r- 2 't+i t+1 (49) Ln(K M) Dot+ml D1 m bi D2 ID a(T ID) + . r;; 2 Hl 2 iim ." 2 R T ID avec : (50) t+l (ZGt1 - BGtl] JZ'-B. Do m = Ln "in#+r Zorn -Born Lp m = Ln Zpt + ( 1 - - ) Bptm P=0,G ZPm + (1 -"\12) Bpm t+1 t+1 t+1 t+1 r t+l '7+1 t+1 t+1 t+1 t+1t+1 Pm Zorn -1 ZGm -1 Ap m Lom AG m Lam Dim T t+1 - t+1 +2 B 2 2B RI Om Gm 2 -\1 2 Hom 2 -\1 2 BG m t+t t+1 t+1 t+1 P m Loin t+1 Lo m t+1 D2m = r- 2 t+1 2 '~+raiGm -alOm '\/2K Tm BGm Bom Dans les équations (50), on reconnaîtra les paramètres de phases adimensionnels B et A définis précédemment par les équations (19) et (20). Dans l'équation (49), le terme Iimt+1représente la contribution des interactions binaires entre constituants et ce terme est nul lorsque les coefficients d'interactions binaires de la représentation détaillée sont tous nuls. Dans un tel cas, les coefficients d'équilibre peuvent être calculés directement à partir des paramètres (E3P mt+1) issus de l'équation d'état EOS~PRO de la simulation dynamique utilisant la représentation groupée. Dans le cas contraire, une expression rigoureuse du terme Iimt+1est: t+1 Nrb ~~~~~ ~~~~~ t+1 Nrb (si) Hl La rn -V Hl -V t+1 t+1 4 rn -V Hl -V Hl Hl ] BOrn j=i BG m j=1 Ce terme fait donc intervenir les fractions molaires, dans la représentation détaillée, des constituants dans les phases hydrocarbure au pas de temps (t+1) xi mt+1, yi mt+1, non encore estimées. Pour éviter une résolution implicite lourde, on peut, dans l'équation (51), t+1 t+1 substituer les fractions molaires xi m , yi m par les fractions molaires au pas de temps (t) issues des équations (46). C'est avec une telle approximation, non limitative de réalisation, que les résultats de delumping : - illustrés Fig. 8. 1 à 8. 16, Fig. 9. 1 à 9. 16 et récapitulés Fig.10 pour une opération de dépressurisation du fluide du cas SPE3, - illustrés Fig. 11.1à 11.16 et Fig. 12. 1 à 12. 16 et récapitulés Fig.13 pour une opération de vaporisation à une pression de 169 bar d'un condensat (le condensat obtenu à 169 bars lors de l'opération de dépressurisation précédente), ont été obtenus. Dans les Fig. 8. 1 à 8. 16 et 9. 1 à 9. 16 d'une part, Fig. 11. 1 à 11. 16 et 12. 1 à 12. 16 d'autre part, la teneur en chaque constituant dans la phase considérée est exprimée en pourcentage molaire de la phase, et, dans chaque figure, l'échelle est dilatée pour être adaptée à la variation de la teneur du constituant dans la phase. L'abscisse sur les Fig. 11.1 à 11.16 et Fig. 12. 1 à 12. 16 est la fraction molaire de vapeur au cours de l'opération de vaporisation, qui croît de 0 (à l'état initial, seule la phase condensat, ou huile, est présente) à 1, lorsque tout le condensat est vaporisé. La Fig. 10 et la Fig. 13 montrent que l'erreur absolue maximale obtenue sur la teneur en chaque constituant (toujours exprimée en pourcentage molaire) reste en moyenne inférieure à 0.1 % durant les opérations de dépressurisation et de vaporisation. On constate donc que les compositions issues de l'opération de delumping sont en bonne cohérence avec les compositions détaillées de référence des fluides, bien que l'approximation utilisée ne soit pas parmi les plus sophistiquées. Une autre approximation possible, toujours non limitative de réalisation, est d'estimer les fractions molaires xi mt+1, yi mt+1 à partir des coefficients d'équilibre Kimt calculables à partir des fractions molaires données par les équations (46) par: /PAS t+l (52) t+1 Zim Kim - 1 X im = t -1) 0 i t+4 m t t+l t t+1 Icim Z im Kim - 1 j' im = -I- (Kim - ) t 1 e t+ m En général, les coefficients d'interactions binaires sont petits devant 1, et un calcul itératif, par exemple par substitutions successives, ne s'avère pas nécessaire, mais reste possible. pp) On détermine, à chaque pas de temps (t+1), la fraction vaporisée (#mt+1) dans chaque maille (m), soit à partir des résultats de la simulation dynamique avec la représentation groupée soit, si une meilleure précision est recherchée, en résolvant l'équation de Rachford-Rice à partir des fractions molaires de chaque composant i dans la composition détaillée globale (zim ) du fluide hydrocarbure dans la maille (m) au pas de temps (t+1) Nrb t+l t+l /e-\ Nrb Zim (Kim - 1) (53) lm t,, t+l (S3) 1= (53) z 1 + (Kim -1) e m Les nombres de moles par phase dans chaque maille à chaque pas de temps (t+1) sont donnés par: (54) No m = N m (1 - e m ) Ngim = N m e qq) On estime la composition détaillée de chaque phase hydrocarbure, à chaque pas de temps (t+1) et dans chaque maille (m), en utilisant les relations suivantes : = t+1 zim (Kjm - 1) (55) x im = Z'm tK'm 1 t+ pour la phase huile 1 +(Kim -1) 0m t+1 t+1 t+1 /c t+1 Kim Zim (Kim - 1) (56) Yim t+1 t+1 pour la phase gaz + (Kim -1) bm rr) Enfin, on obtient les nombres de moles de chaque constituant i dans la t+1 représentation détaillée dans chaque phase au pas de temps t+1, respectivement (Noim) pour la phase huile et (Ngim) par les équations (46) remaniées en: ,#.-. t+i t+i t+l t+i c+1 t+l (57) NOim = No m Xim ; Ngm = Ng m jim ss) Les paragraphes nn) à rr) détaillent le contenu des opérations pour passer du pas de temps t au pas de temps t+1, en supposant connus les nombres de moles de chaque constituant i dans la représentation détaillée dans chaque phase au pas de temps t dans chaque maille de la simulation avec la représentation groupée, respectivement (Noim) pour la phase huile et (Ng;m) pour la phase gaz. Pour démarrer les calculs, il suffit donc de connaître ces quantités à un instant ou pas de temps donné, qui peut être en particulier l'instant initial. Ceci peut être obtenu par divers moyens connus en soi, un moyen consistant en un calcul d'initialisation d'une simulation compositionnelle (calcul de l'état initial uniquement) utilisant la représentation compositionnelle détaillée, les entrées autres que compositionnelles étant les mêmes que celle de la simulation compositionnelle avec la représentation groupée. Dès lors, on sait décrire l'évolution de la composition détaillée dans chaque maille au cours du procédé d'exploitation modélisé dans un modèle compositionnel utilisant une représentation groupée grâce à une mémoire des paramètres de l'équation d'état de la représentation détaillée conservée dans la méthode de calcul des propriétés de phases. Certaines étapes de calcul peuvent apparaître non nécessaires car redondantes. Certaines de ces redondances ont l'avantage de permettre de diminuer l'espace mémoire nécessaire au stockage des données, mais peuvent à l'évidence être évitées si les considérations d'espace mémoire ne sont pas prioritaires. Dans le cas où l'EOS~PRO se révèle valide pour le calcul des coefficients d'équilibre (Ki, m et t+1), le delumping des résultats d'une simulation de réservoir effectué avec la méthode selon l'invention, peut s'effectuer par application de la méthode décrite dans le brevet WO/42.937 précité. Remarques concernant la généralisation de l'approche tt) Les détails fournis dans la description de la méthode utilisent un formalisme molaire ; on aurait pu tout aussi bien utiliser un formalisme massique. uu) On peut remarquer que le raisonnement basé sur l'hypothèse maîtresse M décrite dans le paragraphe d), appliqué à l'ensemble des constituants de la représentation de base, ensemble ci-après désigné par B , dans les paragraphes précédents, aurait pu alternativement être appliqué à un sous-ensemble des constituants de la représentation de base, ci-après désigné par R . L'ensemble complémentaire de R dans B étant désigné par S , les équations (5) seraient alors réécrites (avec P=O,G) comme : 5g e iE S 'tE S MMpo = XPI I + XPH H + 1 Xpei Mmi ; XPI + XPH = 1 - xpi i e e e 'ES e e ES 1 MMPG = ypv MMv + YPI MMI + 1 ypi MMi ; YPV + YPI = 1 - 1 YPi i i Une telle démarche peut être intéressante lorsque, pour le cas traité, il est besoin de garder explicites certains constituants de base dans la représentation groupée. Ainsi, dans un cas pratique où les conditions de réservoir ne permettraient pas de négliger la dissolution de certains constituants (tels le dioxyde de carbone) dans l'eau, ces constituants pourraient former l'ensemble S . Les simples changements de variables indiqués ci-après permettent de ramener les équations (5) à la situation où l'on a considéré l'ensemble des constituants de la représentation de base: (59) ieS e e ic= se MMpo - xpi MMi MMpG - E ypi i Mmpo el = ives e MM pe GR 1 ieS e 1 - Lr XPi 1 - Gr yPi i i e e e e e XPI e XPH e Ypy e Ypl XPI e 1 R = 'ES XPH e 1 R = XPH ; Ypv e 1 R = ics YPI e 1 R = XPI ieS e xPH ieS e YPV ieS e yPI ieS e 1 - L. Xpi expo 1 - L. YPi 1 - L. YPi i i i i Les équations (6) étant résolues avec les nouvelles variables, une transformation inverse de variables permet alors d'enchaîner les calculs des paragraphes g) et suivants. La représentation groupée comporte alors les constituants (V), (I), (H) et les constituants de base de l'ensemble S . Dans chaque phase de l'ensemble E , les fractions molaires des constituants de base formant l'ensemble S sont les mêmes dans la représentation groupée que dans la représentation de base. vv) Les compositions des phases de l'ensemble E dans la représentation groupée étant ainsi obtenues, on peut suivre les grandes lignes des étapes décrites précédemment pour modéliser les propriétés de phases et les équilibres thermodynamiques. Dès lors qu'une équation d'état est recherchée pour la modélisation avec la représentation groupée des propriétés de phases et/ou des équilibres thermodynamiques, l'indication est de conserver pour les constituants de base formant l'ensemble S les paramètres de ces constituants dans l'équation d'état de la description détaillée. Les seuls paramètres à estimer restent alors les paramètres des constituants (V), (I), (H) dans l'équation d'état, et le cas échéant les paramètres d'interactions binaires de ces constituants avec les constituants de l'ensemble S . Lorsque ces derniers sont posés comme nuls, on peut suivre la ligne directrice suivie dans l'ensemble des paragraphes n) à w) qui est de rechercher les paramètres d'équation d'état EOS~PRO utiles aux calculs, dans les conditions de réservoir, de densités des phases (via les facteurs de compressibilité) et de viscosités par des résolutions simples (non itératives) de systèmes d'équations linéaires. REVENDICATIONS 1) Méthode de pseudoïsation pour estimer les propriétés ou le comportement des phases hydrocarbure liquide et/ou vapeur à partir de données relatives à un ensemble de référence constitué de mélanges hydrocarbonés dans une série d'états thermodynamiques découlant de conditions d'exploitation définies d'un gisement souterrain d'hydrocarbures, caractérisée en ce que : - on regroupe chacun des dits mélanges hydrocarbonés en au moins trois constituants (V, I, H), aucun de ces constituants ne correspondant à une sélection particulière de composants ou pseudo-composants de base qui serviraient à une description compositionnelle détaillée des fluides, en considérant que les phases gaz issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le troisième constituant (H) est exclu, et que les phases huile issues de la séparation en conditions de surface de chacun des mélanges hydrocarbonés sont des mélanges dont le premier constituant (V) est exclu ; - on détermine par bilan de matière les compositions des produits de séparation comprenant, pour les produits gazeux, au moins le premier et le deuxième constituant (V, I) en proportions variables, et pour les produits liquides, au moins le deuxième et le troisième constituant (I, H) en proportions variables ; et - on détermine la composition à au moins trois constituants de chaque mélange hydrocarboné de l'ensemble de référence par combinaison des produits de sa séparation au prorata des quantités de chaque produit de séparation. 2) Méthode selon la revendication 1, dans laquelle on regroupe chacun des mélanges hydrocarbonés en seulement trois constituants (V, I, H), les phases gaz issues de la dite séparation étant des mélanges en proportions variables du premier constituant (V) et du deuxième constituant (I), les phases huile issues de la dite séparation étant des mélanges en proportions variables du deuxième constituant (I) et du troisième constituant (H), et on détermine la composition à trois constituants. 3) Méthode selon la revendication 1 ou 2, dans laquelle les conditions de surface sont celles rencontrées ou prévues lors de l'exploitation du gisement. 4) Méthode selon la revendication 1 ou 2, dans laquelle les conditions de surface sont différentes de celles rencontrées ou prévues lors de l'exploitation du gisement 5) Méthode selon la revendication 1 ou 2, dans laquelle le bilan de matière est un bilan massique et dans laquelle on attribue à chacun des trois constituants (V, I, H) une masse molaire après une analyse quantitative des masses molaires des produits de séparation de l'ensemble de référence. 6) Méthode selon la revendication 1 ou 2, dans laquelle on définit les données nécessaires aux calculs d'équilibre et à la modélisation des propriétés de phases dans la représentation groupée, en utilisant les compositions des phases dans la représentation groupée et des données connues, ou estimées, a priori portant au moins sur la densité et la viscosité de phases huile et gaz en équilibre appartenant à l'ensemble de référence. 7) Méthode selon la revendication 6, dans laquelle, lorsque les dites données incluent des données compositionnelles détaillées des phases préalablement représentées par une équation d'état, on définit les paramètres d'une première équation d'état de la représentation groupée, utile à la modélisation des propriétés de phases, en utilisant ces données compositionnelles. 8) Méthode selon la revendication 6, dans laquelle on ajuste les paramètres d'une seconde équation d'état de la représentation groupée, utile aux calculs d'équilibre, pour reproduire les coefficients d'équilibre de la représentation groupée. 9) Méthode selon la revendication 8, dans laquelle on utilise, dans la procédure d'ajustement pour définir les paramètres de la seconde équation d'état de la représentation groupée utile aux calculs d'équilibre, les paramètres par constituant de la représentation groupée dans la première équation d'état. 10) Méthode selon l'une des revendications précédentes, dans laquelle on détermine les coefficients d'équilibre des fluides dans une représentation compositionnelle détaillée, à partir de variables et/ou paramètres intervenant dans le calcul des propriétés de phases, dès lors que les paramètres utiles aux calculs des propriétés de phases dans la représentation groupée ont été estimés de façon à reproduire les paramètres des phases dans l'équation d'état de la description compositionnelle détaillée. 11) Méthode selon l'une des revendications précédentes, comportant un éclatement ou lelumping pour prédire en fonction du temps et dans au moins une zone hermodynamique, une composition détaillée d'un fluide contenu dans un gisement i'hydrocarbures ou produit par au moins un puits. 12) Méthode selon la revendication 11 dans laquelle : on représente le réservoir sous la forme d'un réseau de mailles (m) dont chacune constitue un volume élémentaire rempli de fluides sous forme de une ou plusieurs phases, avec au moins une phase non aqueuse ; on définit, pour chaque zone ou domaine thermodynamique, les fluides par une représentation détaillée de base, de façon à déterminer la quantité de chaque constituant de base (i) dans chaque phase hydrocarbure dans chaque maille (m) à l'instant défini comme initial pour le calcul de delumping ; par zone thermodynamique pour laquelle on adopte une représentation groupée des fluides, on détermine une équation d'état construite, préalablement à la simulation dynamique de réservoir avec la représentation groupée, pour reproduire les paramètres de phases, dans l'équation d'état de la représentation détaillée, des fluides hydrocarbonés au cours de chemins thermodynamiques tenus pour représentatifs de ceux que vont suivre les fluides hydrocarbonés pendant la simulation maillée ; on réalise à un pas de temps t une simulation compositionnelle avec un nombre limité de constituants où les propriétés des phases sont calculées par une équation d'état, ladite simulation permettant de calculer au moins dans chaque maille (m) et à des pas de temps consécutifs, une pression pour une phase hydrocarbure, la température lorsque celle-ci varie, les débits des phases entre mailles et au niveau des perforations de production et d'injection, et les valeurs de paramètres et/ou propriétés de phase intervenant dans l'expression formelle des coefficients d'équilibre de la représentation détaillée, et on mémorise ces diverses quantités ; on estime au pas de temps suivant (t+1) la fraction molaire de chaque constituant i dans la composition détaillée globale du fluide hydrocarboné dans la maille (m) par bilan matière sur la maille (m) ; - on détermine en utilisant les quantités mémorisées, au même pas de temps (t+1) et dans chaque maille (m), les coefficients d'équilibre de chaque constituant (i) dans la représentation détaillée ; - on détermine, au même pas de temps (t+1), la fraction vaporisée dans chaque maille (m) ; et - on estime la composition détaillée de chaque phase hydrocarbonée, au même pas de temps (t+1) et dans chaque maille (m). 13) Méthode selon l'une des revendications précédentes, où différentes étapes sont traduites de façon à produire des résultats intermédiaires ou finaux utilisables dans des équations exprimées selon un formalisme molaire ou massique.