FR 2775095 A1 19990820 FR 9808324 A 19980630 Méthode de simulation pour prédire en fonction du temps une composition détaillée d'un fluide produit par un réservoir La présente invention concerne une méthode de simulation pour prédire, en fonction du temps, une composition détaillée d'un fluide produit par un réservoir et plus particulièrement une composition détaillée d'un fluide contenu dans et produit par un gisement pétrolier dans lequel sont implantés un ou plusieurs puits de production. La simulation compositionnelle d'un gisement pétrolier est couramment utilisée pour fournir les profils prévisionnels de production du gisement qui permettent notamment de déterminer le schéma de production le mieux adapté à ce gisement. La simulation compositionnelle d'un gisement pétrolier est mise en oeuvre non pas en utilisant une description réelle du fluide mais en utilisant un fluide modélisé par un nombre de composants plus réduit que le nombre de composants du fluide réel. En effet, le nombre de composants du fluide de gisement étant relativement grand, la modélisation avec tous les composants conduirait à des temps de calcul trop importants. Ce temps de calcul prohibitif a amené les spécialistes à regrouper les composants purs en pseudo-composants, par exemple pseudo-composant regroupant l'azote et le méthane, pseudo-composant regroupant les hydrocarbures en C3 et C4, etc..., et à réaliser la simulation de gisement sur un nombre réduit N de composants purs et de pseudo-composants. Une telle composition réduite est désignée par l'expression "composition regroupée du fluide" ou "fluide regroupé". Généralement N représente 5 à 10 composants et pseudocomposants, ce qui est considéré comme étant suffisant pour bien représenter le comportement du fluide de gisement aux conditions de fond. Le fluide regroupé est décrit par une des équations d'état bien connues des spécialistes et adaptée audit fluide regroupé. Avant d'effectuer la simulation compositionnelle de gisement, on représente le réservoir ou gisement sous la forme d'un réseau de mailles dont chacune constitue un volume élémentaire dudit réservoir. Le nombre de mailles peut atteindre plusieurs milliers et chaque maille présente des propriétés qui lui sont propres, comme la géométrie, la porosité, la perméabilité, etc... De plus, à certaines mailles peut correspondre au moins un puits d'injection ou de production qui est implanté dans ledit réservoir. La simulation de gisement permet de calculer pour chaque maille un certain nombre de variables principales du fluide de gisement, lesdites variables pouvant être la quantité (nombre de moles) de fluide, la composition ou fraction molaire de chaque composant et pseudo-composant du fluide regroupé et la pression régnant dans chaque maille. Ces variables sont connues à l'instant t=O (début de l'exploitation du gisement) et sont ensuite calculées par la simulation de gisement à chaque instant t ou pas de temps m. Pour chaque maille, il est possible, à partir de ces variables principales, de calculer toute autre propriété du fluide présent dans la maille, comme le nombre de phases, la composition de chaque phase, etc. Pour chaque maille à laquelle est associé un puits de production, il est possible, à partir de ces variables principales et des contraintes de production imposées, de calculer également le débit de production et la composition du fluide produit par ledit puits. Les propriétés thermodynamiques du fluide regroupé peuvent être calculées en utilisant l'équation d'état qui lui est associée. La modélisation à N composants, si elle est suffisante pour représenter le comportement du fluide de gisement aux conditions de fond, n'est plus appropriée pour simuler le comportement du fluide dans les installations pour son exploitation en surface, laquelle exploitation nécessite de connaître une composition plus détaillée à Q composants et/ou pseudocomposants du fluide produit par le gisement, Q étant supérieur à N et par exemple de l'ordre de 16 à 30. Jusqu'ici, on obtenait cette composition détaillée à Q composants et/ou pseudo-composants à partir de la composition regroupée N en supposant que la composition des pseudo-composants restait constante dans le temps. Une telle manière de procéder est source d'erreurs car la composition des pseudo-composants varie dans le temps. notamment en fonction de la pression du gisement ou suite à l'injection d'un gaz dans le gisement. De ce fait, la composition détaillée obtenue reste très approximative et ne permet pas de prédire de manière satisfaisante le comportement du fluide dans les conditions d'exploitation en surface bien des années après le début de la production. La présente invention a pour but de proposer une méthode de simulation pour prédire une composition détaillée du fluide produit par un réservoir en tenant compte des paramètres du réservoir et qui soit beaucoup plus précise, une telle prédiction étant dénommée "delumping" par les spécialistes. Le "delumping" consiste à prédire quels auraient été les résultats d'une simulation de réservoir si la modélisation détaillée du fluide avait été utilisée et ce, à partir des résultats d'une simulation réalisée en utilisant la modélisation regroupée dudit fluide. L'invention a pour objet une méthode du type consistant à: a) représenter le réservoir sous la forme d'un réseau de mailles (j) dont chacune constitue un volume élémentaire rempli de fluide, b) définir le fluide par une modélisation regroupée à N composants et pseudo-composants (i) et déterminer une équation d'état décrivant le fluide dans cette modélisation regroupée, c) définir également le fluide par une modélisation détaillée à Q composants et/ou pseudo-composants, Q étant supérieur à N, et déterminer une autre équation d'état décrivant le fluide dans cette modélisation détaillée, d) réaliser, de façon connue en soi, une simulation compositionnelle du fluide regroupé à N composants et pseudo-composants (i), ladite simulation compositionnelle permettant de calculer au moins pour chaque maille (j) et à des pas de temps consécutifs (m, m+l, ...) la fraction vaporisée (OJ ), les constantes d'équilibre liquide-vapeur (km) de chaque composant (i), les débits d'injection ou de production (s,J ) et pour chaque paire de mailles (j, h) les débits des phases liquide (Uogh) et vapeur (ugmssh) du fluide à N composants et pseudo-composants, et elle est caractérisée en ce qu elle consiste en outre à: e) déterminer, à chaque pas de temps (m) et pour chaque maille (j) la composition des phases liquide et vapeur du fluide pour la modélisation détaillée à Q composants et/ou pseudo-composants, à partir des valeurs de la fraction vaporisée (ou ) et des constantes d'équilibre (km) du fluide regroupé, f) évaluer pour chaque maille, au pas de temps (m+ l ) la quantité molaire de chacun des Q composants et/ou pseudo-composants du fluide détaillé à partir des valeurs correspondantes au pas de temps (m) des débits des phases du fluide regroupé et des compositions des phases liquide et vapeur du fluide détaillé déterminées à l'étape e), et à g) évaluer, pour chaque puits de production, la composition détaillée du fluide produit entre des instants t et t' correspondant aux pas de temps m et m+l, à partir des débits des phases du fluide regroupé et des compositions des phases liquide et vapeur du fluide détaillé déterminées à l'étape e). Selon une autre caractéristique de la présente invention, les étapes e, f et g sont mises en oeuvre en meme temps que la réalisation de la simulation compositionnelle du fluide regroupé. Selon une autre caractéristique de la présente invention, les résultats de l'étape d sont mémorisés puis utilisés ultérieurement pour la mise en oeuvre des étapes e, f et g. La méthode selon la présente invention est mise en oeuvre pour un réservoir constitué, par exemple, par un gisement pétrolier. De manière connue, on réalise une simulation de gisement. A cet effet, on représente le gisement sous la forme d'un réseau maillé, certaines des mailles ou un groupe de mailles étant associés à un des puits de production implantés dans le gisement pétrolier à exploiter. Pour réduire le temps de calcul, la simulation compositionnelle du gisement est effectuée sur un nombre limité N de composants et pseudocomposants, par exemple de 5 à 10, définis de la manière indiquée précédemment, ces composants et pseudo-composants étant sélectionnés en fonction de la nature du gisement. Comme indiqué précédemment, le fluide regroupé à N composants et pseudo-composants est décrit par une équation d'état qui pourrait être par exemple celle de PENG-ROBfNSON adaptée à ladite composition regroupée. La simulation compositionnelle du gisement permet de calculer, à chaque pas de temps et pour chaque maille, un certain nombre de variables principales du fluide regroupé, lesquelles variables, dans ce cas, sont: la quantité de fluide: fm la fraction molaire du composant i : 7m avec i Par des mesures effectuées au préalable, on connaît la valeur desdites variables principales au début de l'exploitation du gisement, c'est-àdire à l'insant t=0 ou au pas de temps m=0. Pour passer du pas de temps m au pas de temps m+l, la simulation compositionnelle comprend plusieurs étapes détaillées ci-après. Dans une première étape et pour le pas de temps m, un calcul d'équilibre liquide-vapeur (flash) est effectué sur le fluide regroupé, pour chaque maille j afin de déterminer la fraction vaporisée #jm la fraction molaire de chaque composant i dans la phase liquide (x Jm avec i = 1. .N) et la fraction molaire dudit composant i dans la phase vapeur ou gaz (y" avec i = 1...N). Ces diverses fractions sont déterminées avec le système d'équations suivant: dans lesquelles 'Pi est la fugacité du composant i. On mémorise la fraction vaporisée #jm et les constantes yijm d'équilibre kijm = avec i variant entre 1 et N, afin de les utiliser dans une xijm autre étape dénommée "delumping" qui sera explicitée ultérieurement. Dans une deuxième étape et pour chaque maille j, on évalue certaines des propriétés de chaque phase liquide et gaz, à savoir masse volumique (#ojm et #gjm), viscosité, saturation et perméabilité relative. Dans une troisième étape, pour chaque maille j et au pas de temps m, on calcule les coefficients de l'équation pour la pression. Les équations pour la pression d'une maille j à l'instant t+At ou au pas de temps m+l, sont Pour résoudre ces équations, les coefficients à calculer sont - la compressibilité totale CTjm (fluide + roche), - les transmissivités généralisées - le terme source Dans ces équations: h sont les mailles voisines de la maille j, h E J(j), &alpha;ijm est le volume molaire partiel du composant i dans la maille j au pas de temps m, m et #gjm sont les mobilités des phases liquide (o) et gaz (g), t,h est la transmissivité entre les centres des mailles j et h, 1 étant la maille amont. Dans une quatrième étape, on résout pour toutes les mailles les équations pour la pression à l'instant t+At c'est-à-dire au pas de temps m+l. Ces équations forment un système linéaire du type A. où p = (p1m+1,. .pJm+1)T, J étant le nombre de mailles du réseau. Ces équations sont résolues par un procédé itératif standard. Dans une cinquième étape, on calcule au pas de temps m+l, les débits des phases entre chaque paire de mailles adjacentes et pour chaque puits associé à partir, par exemple, de la loi de DARCY. uojhm = -#ojm#olm.tjh(phm+1-pjm+1) ugjhm=-#gjm#glm.tjh(phm+1-pjm+1) p étant la masse volumique et 1 étant la maille amont avec l=j ou h. Selon une caractéristique de l'invention, on stocke ou mémorise les débits des phases qui seront utilisés dans le delumping. Dans une sixième étape, on évalue la quantité f et la composition z du fluide regroupé, dans chaque maille j et au pas de temps m+ l à partir des valeurs au pas de temps m, des débits des phases évalués à la cinquième étape et des compositions des phases calculées à la première étape. Pour ce faire, on applique les lois de conservation définie par les équations suivante s Par un processus itératif, on répète les étapes ci-dessus pour calculer les variables principales aux temps t+At, t+2At, t+3At, etc... ou encore aux pas de temps correspondants m+l, m+2, m+3, etc..., Incrément At pouvant être constant ou variable. Comme l'écoulement peut s'effectuer dans n'importe quelle direction entre les mailles voisines, il faut traiter toutes les mailles à chaque pas de temps. Selon la présente invention, on définit une modélisation détaillée à Q composants et/ou pseudo-composants ainsi qu'une autre équation d'état décrivant le fluide ainsi modélisé. Pour cela, on choisit préalablement les composants et/ou pseudocomposants dont on veut connaître les caractéristiques et qui sont nécessaires pour représenter le fluide de production dans les conditions de surface. Le nombre Q qui est supérieur au nombre N de la composition regroupée est généralement de l'ordre de 16 à 30. On connaît la composition du fluide détaillé à Q composants et/ou pseudo-composants au pas de temps m=0. Les données stockées aux première et cinquième étapes sont utilisées pour les étendre aux Q composants et/ou pseudo-composants de la composition détaillée, de manière à effectuer l'opération de "delumping". Pour l'opération de "delumping", on fait les hypothèses que le fluide du réservoir comporte seulement deux phases non aqueuses (liquide et gaz) et que l'eau ne joue aucun rôle dans l'équilibre liquide-gaz et qu'aucun des autres composants du fluide ne se dissout dans l'eau. Dans ce qui suit, les lettres en majuscule sont utilisées pour représenter le fluide détaillé et les lettres en minuscule sont utilisées pour représenter le fluide regroupé. Pour réaliser le "delumping", on fait, en outre, l'hypothèse qu'à chaque pas de temps m, - la fraction molaire de la phase vaporisée dans chaque maille est indépendante du nombre de constituants choisis pour modéliser le fluide, #jm=#jm - le nombre de moles de liquide et de vapeur s'écoulant entre les mailles (U, u) et dans ou hors des puits (S, s) est indépendant du nombre de constituants choisis pour modéliser le fluide. Uogh UOJ;1 et Ugh UgXh Sn; = Sn; et = = rn Pour effectuer le "delumping", il est nécessaire de connaître les constantes d'équilibre K pour le fluide détaillé en plus de la fraction molaire O- déjà déterminée par l'égalité #=#. Le calcul des constantes d'équilibre peut être effectué à partir de l'équation générale comme indiqué par C. LEIBOVICI dans son article "A Consistent Procedure for Pseudo-Component Delumping" paru dans FLUID PHASE EQUILIBRIA, 117 (1996), 225-232. Pour une équation d'état à deux paramètres, comme celle de PENG ROBINSON, l'équation généralisée de LEIBOVICI peut être écrite pour le fluide détaillé de la manière suivante, pour autant qu'il n'y ait pas de coefficient d'interaction binaire où aI et bI sont les paramètres de l'équation d'état pour le composant I et où les paramètres #jm,ssjm et &gamma;jm sont des paramètres calculés à partir des constantes d'équilibre (kjj) pour le fluide regroupé en minimisant la fonction: Les paramètres rljn,ssjm et &gamma;jm ainsi obtenus seront approximativement les mêmes pour les deux modélisations détaillée et regroupée du fluide, lesdits paramètres dépendant seulement des paramètres de phase de la pression et de la température. En présence de coefficient d'interaction binaire, les paramètres rl, ss et y peuvent etre déterminés par régression. De ce fait, lorsqu'on connaît la composition détaillée ZIjm dans la maille j au pas de temps m, on peut estimer les compositions détaillées des phases liquide et vapeur dans la même maille j et pour le même pas de temps m. I étant le composant et/ou pseudo-composant de la composition détaillée Ces compositions peuvent être normalisées afin que leur somme soit égale à 1 si nécessaire. Lorsqu'on connaît les compositions des phases détaillées pour toutes les mailles du réseau au pas de temps m, il est possible d'estimer les compositions détaillées Z1m+I au pas de temps consécutif m+ I à partir des équations ci-après avecl=1 Q dans lesquelles j'=j pour un débit de la maille j vers la maille h ou dans le puits et j'=h pour un débit de la maille h vers la maille j, et j' correspondant au fluide injecté dans le cas des puits d'injection, S étant alors négatif, et F est la quantité molaire de fluide dans chaque maille avec F = f. Le débit molaire du composant détaillé I produit par le puits W au pas de temps m est donné par: et la composition détaillée du fluide produit par le puits est obtenu en normalisant les débits ci-dessus Grâce à la présente invention, on obtient une composition détaillée à Q composants et/ou pseudo-composants qui permet de mieux définir les profils de production prévisionnels et ce, avec une précision accrue comparée aux méthodes antérieures. Cette plus grande précision des profils de production prévisionnels permet un choix de schéma de développement (type et taille d'installations de surface, nombre de puits, etc...) plus fiable, avec pour conséquence une augmentation de la rentabilité et une diminution du risque économique. Dans ce qui précède, certaines valeurs obtenues au cours de la simulation de gisement sont stockées ou mémorisées pour être utilisées ultérieurement dans le delumping. Il va de soi que ces mêmes valeurs pourraient être utilisées directement si on décide d'effectuer le delumping en même temps que ladite simulation de gisement. Un autre avantage de la présente invention est que la méthode peut être utilisée pour calculer la concentration d'hydrocarbures présents à l'état de trace dans le fluide produit sans avoir à effectuer une simulation de réservoir incorporant ces hydrocarbures, sous réserve que leurs paramètres d'équation d'état soient connus. Il est à noter qu'on peut traiter les cas où la composition initiale du fluide varie avec la profondeur ou latéralement à l'intérieur du réservoir, à condition de connaître la composition initiale détaillée en tout point du réservoir. Deux exemples sont donnés ci-dessous pour démontrer la rapidité et la précision de la méthode selon l'invention, les tests ayant été effectués sur un modèle de réservoir parallélépipédique comprenant quatre couches de différentes perméabilités, avec un pendage dans la direction x. Les prédictions ont été calculées pour déterminer le comportement du fluide dans les conditions en surface sur une période de 30 ans à partir du début d'exploitation, le fluide étant constitué d'une huile volatile. Dans le premier exemple, le réservoir comporte un contact eauhuile, alors que dans le deuxième exemple, le réservoir est dépourvu d'eau et ne contient que de l'huile et du gaz. Dans la modélisation du fluide détaillé, on définit le fluide par 16 composants et/ou pseudo-composants, comme indiqué dans le tableau 1, et on considère que les coefficients d'interaction binaire sont égaux à zéro. Dans la modélisation du fluide regroupé, on définit le fluide par 7 pseudo-composants, comme indiqué dans le tableau 2. A la température du réservoir, la pression du point de bulle de l'huile volatile est de 338,5 bars pour le fluide détaillé et de 336,6 bars pour le fluide regroupé. Exemple I La pression initiale du réservoir est de 400 bars et le réservoir est épuisé à travers un seul puits qui est foré près du sommet du réservoir. La production d'huile du puits est de 300 m3 par jour et la production maximum de gaz est 106 m3 normaux par jour. Une première simulation compositionnelle du réservoir fut d'abord effectuée en utilisant la modélisation détaillée du fluide à 16 composants et/ou pseudo-composants à l'aide des équations d'état adaptées à chacun des composants. Au début, le puits produit de l'huile au débit précité. Après environ 1250 jours, la pression du point de bulle est atteinte. Le débit de gaz augmente continuellement jusqu'à la valeur maximum de 106 N m3/jour qui est atteinte au bout d'environ 5000 jourS. Le puits débite à ce maximum alors que le débit d'huile diminue jusqu'à ce que la limite de pression au fond du puits soit atteinte au 6500ème jour environ. Après quoi, les deux débits d'huile et de gaz diminuent, la simulation se terminant à la fin de 10950 jours (30 ans). A ce moment là, la pression moyenne du réservoir est juste au-dessus de 50 bars. La production d'eau apparat au bout de 5000 jours. Une deuxième simulation compositionelle du réservoir est ensuite répétée avec la modélisation du fluide regroupé à 7 pseudocomposants. Les profils de production réalisés sont sensiblement les mêmes que ceux de la simulation avec le fluide détaillé. Cela suggère que la modélisation à 7 pseudo-composants est adéquate pour modéliser le comportement de phase du fluide dans les conditions du réservoir, tout au moins à des fins d'étude du gisement. Lorsqu'on applique le "delumping" aux résultats obtenus avec la deuxième simulation effectuée sur le fluide regroupé, on constate qu'il y a concordance entre les résultats obtenus après "delumping" suivant l'invention et ceux obtenus avec la première simulation qui est effectuée à l'aide des équations d'état adaptées à chacun des 16 composants et/ou pseudocomposants. La différence réside dans le fait que la première simulation du fluide détaillé nécessite un temps de calcul de 730 secondes alors que le temps de calcul global pour réaliser la deuxième simulation (232 s) et du delumping (27 s) n'est que de 259 secondes, c'est-à-dire que le temps de calcul est réduit de 65 % tout en ayant d'excellents résultats. Exemple 2 Dans cet exemple, le réservoir est au début totalement rempli d'huile. L'huile est produite par un puits débouchant à la partie inférieure du réservoir, avec un débit de 1000 m3/jour. Du gaz est injecté à travers un puits, dans la partie supérieure du réservoir, ledit gaz d'injection n'étant pas constitué par un mélange de composants regroupés, mais comprenant, en pourcentage molaire, 80 % C1, 15 % C2 et 5 % C3, cette composition étant utilisée pour la simulation avec le fluide détaillé ainsi que pour le "delumping". Pour la simulation de réservoir avec le fluide regroupé, la composition du gaz injecté est représenté en pourcentage molaire par 80 % N2 C1, 15 % CO2 C2 et 5 % C3C4 ; l'injection de gaz est réalisée à raison de 500 m3/jours pendant les 15 premières années et de 900 m3/jours les 15 années suivantes. Comme dans l'exemple 1, une première simulation avec un modèle de fluide détaillé à 16 composants est effectuée. Au début. le débit du gaz d'injection est insuffisant pour maintenir la pression du réservoir, de ce fait, la pression du point de bulle est vite atteinte et le rapport gaz-huile augmente. Après 15 ans, lorsque la pression du réservoir est de 120 bars, on augmente le débit du gaz d'injection de sorte qu'au bout de 30 ans, à la fin de la simulation, la pression du réservoir est d'environ 100 bars. Une deuxième simulation avec le modèle du fluide regroupé est effectuée. On constate les mêmes effets que dans l'exemple 1, c'est-à-dire que le modèle à 7 pseudo-composants est encore adéquat pour modéliser le comportement de phase du fluide dans les conditions du réservoir. Lorsqu'on applique le "delumping" selon l'invention aux résultats de la deuxième simulation, on constate que, grâce à l'invention, le temps de calcul pour effectuer la deuxième simulation (173 s) et le "delumping" (18 s) est de 191 secondes, ce qui représente environ 15 % du temps de calcul nécessaire pour effectuer la première simulation avec le modèle du fluide détaillé, les résultats étant excellents, comme dans l'exemple 1. Lorsqu'on introduit des coefficients d'interactions binaires, on constate que la concordance est encore satisfaisante, même si elle n'est pas du même ordre que celle des exemples 1 et 2. Pour les composants les plus lourds, la concordance est presque totale. Seuls les pseudo-composants CO2-C2 du fluide regroupé ont donné une concordance relativement faible. Par contre, le gain sur le temps de calcul est encore très appréciable puisqu'il ne représente que 35 % du temps de calcul de la simulation effectuée sur le modèle de fluide détaillé avec des équations d'état. TABLEAU I Composant Composition huile initiale (frac. mol.) N2 0,00080 CO2 0,02929 C1 0,64641 C2 0,04165 C3 0,04649 1C4 0,01569 NC4 0,02690 IC5 0,01024 NC5 0,00948 C8 0,01915 C7 0,02003 C8 0;01995 Cg 0,01391 C10 0,01120 CNG1* 0,06644 CNH?* 0,02309 TABLEAU 2 Composant Composition huile initiale (frac. mol.) N2C1 0,64721 C 2C2 0,07094 C3C4 0,08908 C5C7 0,05890 C8C10 0,04435 CNG1* 0,06644 CNH?* 0,02309 * CNG1 et CNH2 représentent la partie lourde du fluide (C1 1+) divisée en deux pseudo-composants, I'un (CNG1) trouvé principalement dans le gaz et l'autre (CNH2) dans l'huile. REVENDICATIONS 1. Méthode de simulation pour prédire, en fonction du temps, une composition détaillée d'un fluide produit par un réservoir, du type consistant à : a) représenter le réservoir sous la forme d'un réseau de mailles (j) dont chacune constitue un volume élémentaire rempli de fluide, b) définir le fluide par une modélisation regroupée à N composants et pseudo-composants (i) et déterminer une équation d'état décrivant le fluide dans cette modélisation regroupée, c) définir également le fluide par une modélisation détaillée à Q composants et/ou pseudo-composants, Q étant supérieur à N, et déterminer une autre équation d'état décrivant le fluide dans cette modélisation détaillée, d) réaliser, de façon connue en soi, une simulation compositionnelle du fluide regroupé à N composants et pseudo-composants (i), ladite simulation compositionnelle permettant de calculer pour chaque maille (j) et à des pas de temps consécutifs (m, m+l, ...) au moins la fraction vaporisé (OJ ), les constantes d'équilibre liquide-vapeur (kijm) de chaque composant (i), les débits d'injection ou de production (sjm) et pour chaque paire de mailles (j, h) les débits des phases liquide (UOJh) et vapeur (Um i du fluide à N composants et pseudo-composants, et elle est caractérisée en ce qu'elle consiste en outre à e) déterminer, à chaque pas le temps (m) et pour chaque maille (j) la composition des phases liquide et vapeur du fluide pour la modélisation détaillée à Q composants et pseudo-composants, à partir des valeurs de la fraction vaporisée (OJ ) et des constantes d'équilibre (km) du fluide regroupé, f) évaluer pour chaque maille, au pas de temps (m+ 1), la quantité molaire de chacun des Q composants et/ou pseudo-composants du fluide détaillé à partir des valeurs correspondantes au pas de temps (m), des débits des phases du fluide regroupé et des compositions des phases liquide et vapeur du fluide détaillé déterminées à l'étape e), et à g) évaluer, pour chaque puits de production, la composition détaillée du fluide produit entre des instants t et t' correspondant aux pas de temps m et m+l, à partir des débits des phases du fluide regroupé et des compositions des phases liquide et vapeur du fluide détaillé déterminées à l'étape e). 2. Méthode selon la revendication 1, caractérisée en ce que les étapes e, f et g sont mises en oeuvre en même temps que la réalisation de la simulation compositionnelle du fluide regroupé. 3. Méthode selon la revendication 1, caractérisée en ce que les résultats de l'étape d sont mémorisés puis utilisés ultérieurement pour la mise en oeuvre des étapes e, f et g.