FR 2811430 A1 20020111 FR 0009008 A 20000710 La présente invention concerne une méthode de modélisation permettant de prédire en fonction du temps la composition détaillée de fluides produits par un gisement souterrain en cours de production, et notamment d'un gisement d'hydrocarbures en cours de production. Un tel modèle permet à l'ingénieur de gisement d'établir des profils compositionnels détaillés qui, sont 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 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 la loi de conservation des masses dans chaque unité de volume, des relations thermodynamiques régissant l'évolution des paramètres 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 (E), un constituant hydrocarbure lourd (L), et un constituant hydrocarbure léger (V). Dans un modèle de type B.O. dit " strict ", le constituant (E) n'est présent que dans la phase eau, le constituant (L) n'est présent que dans la phase liquide hydrocarbure (appelée huile ou condensat), le constituant (V) se partage entre les phases hydrocarbures liquide et vapeur (phase appelée gaz). Même si l'usage des modèles B.O. est à déconseiller dans certains cas de gaz à condensats, il n'en reste pas moins applicable à un grand nombre de cas industriels. Un autre modèle de simulation bien connu dit " compositionnel " est aussi utilisé, o les fluides hydrocarbures sont représentés par un nombre plus élevé de composants, au moins trois, souvent davantage. La modélisation des flux de ces fluides mieux détaillés se traduit par des temps de calcul très importants (beaucoup plus longs que ceux requis pour une modélisation de type B.O.) du fait du plus grand nombre de constituants mais aussi parce 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. Pour des raisons pratiques, les fluides en place sont décrits comme constitués d'un nombre de composants ou pseudo-composants beaucoup plus réduit que le nombre réel de composants, de façon que les calculs de modélisation puissent être menés dans un délai raisonnable. Une composition réduite à quelque 5 ou 10 pseudo- composants s'avère en général suffisante pour représenter le comportement des fluides dans le gisement. Par la demande de brevet WO 99/42937 et aussi le papier de C. Leibovici et J. Barker "A Method for Delumping the Results of a Compositional Reservoir Simulation" SPE 49068, présenté à la SPE Annual Technical Conference and Exhibition New Orleans, 27-30 Septembre 1998, on connaît une méthode de prévision de l'évolution de la composition détaillée au cours du temps à partir des calculs effectués dans une simulation de type compositionnel de fluides décrits par un certain nombre réduit de pseudo- composants (principe de représentation 'groupée' ou 'regroupée' obtenue par une opération dite de 'lumping'), le nombre de composants étant au moins de trois. La méthode permet donc 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. C'est une opération connue des gens de l'art sous l'appellation de "delumping", que l'on désignera par la suite aussi par "éclatement". Le principe de l'étape d'éclatement antérieure est de calculer le coefficient ADo et les n coefficients ADp (soit n+l 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.: n (1) Ln(ki) = ADo + _ADp rlpi p=l o les IIpi sont des paramètres fixes de caractérisation du constituant i dans l'équation d'état pour une représentation thermodynamique donnée, à partir des constantes d'équilibre ki de chaque constituant de la représentation thermodynamique regroupée calculées pendant la simulation compositionnelle dans chaque maille et à chaque pas de temps. Si Nrg est le nombre de composants de la représentation thermodynamique regroupée, on dispose donc de Nrg équations pour déterminer n+1 coefficients. Une condition nécessaire est donc que Nrg soit au moins égal à n+l. Pour l'équation d'état de Peng-Robinson à deux paramètres, il faut donc une représentation thermodynamique regroupée à au moins trois composants. Une fois le coefficient ADo et les n coefficients ADp calculés, on les utilise pour calculer les constantes d'équilibre des composants de la représentation thermodynamique détaillée (Nrd composants) en appliquant l'équation (1) aux Nrd composants avec leurs propres paramètres fixes de caractérisation dans la représentation thermodynamique détaillée. En utilisant a) les constantes d'équilibre pour la représentation thermodynamique détaillée ainsi déterminées, b) les flux entre chaque maille et aux puits, c) la fraction vapeur dans chaque maille issus de la simulation compositionnelle regroupée, et, d) la composition détaillée globale dans chaque maille et aux puits injecteurs en début de pas de temps, on estime alors dans chaque maille la composition détaillée de chaque phase hydrocarbure au pas de temps t et la composition détaillée globale de chaque maille au pas de temps (t+1) suivant. L'un des intérêts de cette méthode est qu'il n'est pas nécessaire, dans l'étape d'éclatement (" delumping ") à chaque pas de temps, de résoudre l'équation d'état, que ce soit pour la représentation regroupée ou la représentation détaillée, ce qui constitue une économie de temps de calcul. Un inconvénient de la méthode est qu'elle n'est pas applicable aux simulations de type B.O, puisque les équations d'état pratiques ont au moins deux paramètres. Par la demande de brevet WO 98/5710726, on connaît également une méthode pour prédire l'évolution de la composition détaillée au cours du temps à partir des calculs de flux effectués dans une simulation de type B.O. o les phases hydrocarbures sont décrites par seulement deux composants (L) et (V). L'inconvénient de la méthode est qu'elle requiert l'utilisation de l'équation d'état de la représentation détaillée à chaque pas de temps, et donc d'être coûteuse en temps de calcul pendant l'étape de "delumping". La méthode selon l'invention La méthode selon l'invention combine une simulation prévisionnelle du comportement du gisement par une modélisation de type B.O. avantageuse puisqu'elle ne nécessite aucun recours à une équation d'état, et une étape d'éclatement ou de "delumping" qui ne nécessite à aucun pas de temps une quelconque résolution d'une équation d'état. Dans une modélisation de type B.O. on introduit, entre autres paramètres, les paramètres thermodynamiques des fluides (tels que viscosité, facteur de volume, densité, rapport de dissolution gaz-huile, etc.) dans un fichier des données d'entrée du simulateur (sous forme de tables, et ou par corrélation, comme fonction de la pression, de la température si celle-ci varie) et, le cas échéant (pour le modèle de type B.O. étendu), d'un paramètre additionnel conservant une mémoire de la composition du gaz tel par exemple que la densité du gaz). L'idée essentielle est ici d'introduire simplement dans les données d'entrée les (n+l) colonnes ou corrélations supplémentaires correspondant aux différences ADp, par exemple en les incluant dans un fichier annexe au fichier principal utilisé par la simulation. I est également possible d'introduire 2 x (n+1) colonnes ou corrélations supplémentaires correspondant aux paramètres Dp de chaque phase. Une condition pratique pour que des données puissent être entrées sous forme de table est que les dépendances fonctionnelles introduites soient suffisamment continues. Par là, on s'affranchit de la nécessité du calcul des constantes d'équilibre par le modèle. La méthode de simulation selon l'invention permet de prédire en fonction du temps et dans au moins une zone thermodynamique, 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. Elle est caractérisée en ce que: - on discrétise le gisement par un maillage, chaque maille contenant une ou plusieurs phases dont au moins une phase non aqueuse; - on détermine la variation de paramètres thermodynamiques des phases non aqueuses nécessaires à une modélisation de type " Black Oil " (B.O.) au cours d'étapes d'un chemin thermodynamique suivi par les fluides dans le gisement (tels par exemple que la viscosité des phases, le rapport de dissolution du gaz dans l'huile, les facteurs de volume de l'huile ou de compressibilité du gaz, etc.); - on définit les fluides par une représentation détaillée à Nrd composants et/ou pseudo-composants; - on détermine au moins une équation d'état à n paramètres permettant de simuler au moins durant la phase de préparation des données d'entrée (simulation qui est nécessaire strictement que durant cette phase de préparation des données d'entrée, et donc non nécessaire de façon répétitive au cours ou après la simulation de type B.O.), le comportement thermodynamique des fluides au cours des étapes d'un chemin thermodynamique suivi, permettant de générer, par zone ou domaine thermodynamique, n+1 relations fonctionnelles supplémentaires (sous forme de tables ou de corrélations); - on traduit le comportement thermodynamique de chaque phase non aqueuse en des données d'entrée convenant pour une représentation thermodynamique de type B.O., données d'entrée que l'on complète par les données d'entrée supplémentaires convenant à l'opération d'éclatement - on réalise une modélisation de type B.O., permettant de déterminer dans chaque maille et à des pas de temps consécutifs, des caractéristiques thermodynamiques de chaque phase non aqueuse et des données représentatives de déplacement des phases dans le gisement; et - on réalise une opération d'éclatement pour obtenir la composition détaillée des fluides dans chacune des mailles. Pour caractériser le chemin thermodynamrnique, on détermine par exemple les variations de la composition des fluides et des propriétés des différentes phases. Au cours de la simulation de type B.O., on détermine par exemple dans chaque maille des paramètres thermodynamiques tels que la pression (pj), la température (Tj) si celle-ci varie, les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire de mailles(j,h), les débits volumiques des phases liquide (uojh) et vapeur (Ugjh). L'opération d'éclatement comporte par exemple une détermination des constantes d'équilibre à partir de données d'entrées spécifiques à l'opération d'éclatement, et une conversion de résultats exprimés en volume par des résultats utilisables dans des équations de conservation des quantités molaires ou massiques. Présentation sommaire des figures D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés o: - les Fig. 1 à 3 montrent respectivement les variations en fonction de la pression de trois coefficients ADo, AD1, AD2, au cours de la vaporisation différentielle d'une huile de gisement simulée à l'aide de l'équation de Peng-Robinson à deux paramètres. Description détaillée Les différentes étapes de la méthode de modélisation permettant notamment de prédire en fonction du temps la composition détaillée de fluides produits par un gisement d'hydrocarbures en cours de production sont les suivantes: a) On représente tout d'abord le réservoir sous la forme d'un réseau de mailles (j) 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 sont 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. b) On établit, d'une façon connue en soi, le comportement thermodynamique de la (des) phase(s) non aqueuse(s) nécessaire à une représentation thermodynamique de type B.O.. Il s'agit de décrire la dépendance de paramètres de phase avec la pression, la température si celle varie en cours de simulation et, éventuellement, pour les modélisations de type B. O. non strict, un indicateur d'évolution de la composition (par exemple le contenu en condensat dissous dans la phase gaz) au cours d'un chemin (évolution) thermodynamique représentatif de celui que va suivre le fluide dans le gisement. Le chemin thermodynamique le plus courant est celui de l'opération dite différentielle au cours de laquelle la pression du mélange est progressivement abaissée à température constante: sous le point de saturation, l'une des phases est évacuée à chaque palier de pression (pour la vaporisation différentielle d'une phase liquide, la phase gaz est évacuée); durant cette opération, les compositions et propriétés des phases (viscosités, masses volumiques, facteur de volume, rapport de dissolution gaz-huile de la phase liquide, etc.) sont déterminées. Il est à noter que l'on peut également traiter les cas qui requièrent l'usage de plusieurs représentations thermodynamiques " Black Oil " par exemple, si la composition en début de simulation varie avec la profondeur ou latéralement à l'intérieur du réservoir, ou si l'on peut distinguer plusieurs chemins thermodynamiques locaux durant la modélisation. 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. L'origine des données sur le comportement du fluide de gisement peut être de nature expérimentale (expériences de laboratoire), elles peuvent être issues de mesures sur champ, ou prédites à partir de données frustres et d'hypothèses complémentaires en simulant numériquement le comportement du fluide de gisement. Parmi les propriétés de phase ainsi habituellement décrites en entrées des modélisations de type B.O. on trouve la viscosité des phases, le rapport de dissolution du gaz dans l'huile, le facteur de volume de l'huile, le facteur de compressibilité du gaz. c) On définit, avant la simulation B.O. et pour chaque zone ou domaine thermodynamique, le fluide par une représentation détaillée à Nrd composants et/ou pseudo-composants, d) Par zone thermodynamique, on détermine une équation d'état permettant de reproduire préalablement, avant la simulation B.O., le comportement thermodynamique du fluide au cours du chemin thermodynamique représentatif de celui que va suivre le fluide dans le gisement. e) On reproduit avec l'équation d'état les étapes (souvent des paliers de pression) du chemin thermodynamique de chaque zone, en particulier sous la pression de saturation o deux phases sont présentes. Lors de cette simulation thermodynamique préalable, on utilise, à la pression de saturation du fluide et en dessous, à chaque étape du chemin thermodynamique, les constantes d'équilibre obtenues à partir des résultats de la simulation, pour calculer le coefficient ADo et les n coefficients ADp de l'équation (1) correspondant à l'équation d'état choisie pour décrire le comportement du fluide hydrocarbure par exemple en minimisant la fonction: Nrd [ n (2) O(ADo, ADI,..,ADp,..,ADn) - = DO + ZADp 1-pi- Ln(ki) i=l p=l L'équation (1) ayant été obtenue à partir d'une équation sur la fugacité du constituant i dans une phase o apparaissent les coefficients Do et les n coefficients Dp propres à la phase: È(3) L{ =nDo+Z Dpfpn (3) Ln Xi= Do + Dp FI"pi, Ix+p) p=l on peut alternativement déterminer ces coefficients qui apparaissent comme des paramètres de phases, l'inconvénient étant d'augmenter le volume des données d'entrées. On dispose alors à chaque étape du chemin thermodynamique, prévu comme représentatif de celui que va suivre le fluide dans le gisement, de la dépendance en fonction de la pression, de la température et, le cas échéant, d'autres indicateurs nécessaires, des propriétés des phases nécessaires à la simulation B.O., et, sous la pression de saturation, du coefficient ADo et des n coefficients ADp (ou des n+l paramètres Do et Dp, p variant de 1 à n, par phase) qui seront utilisés pour calculer les constantes d'équilibre lors de l'étape de "delumping" décrite ciaprès. Les variations des différences ADp (p= 0 à 2) pour une opération différentielle sur une huile de gisement simulée à l'aide de l'équation de Peng- Robinson à deux paramètres, sont illustrées aux Fig. 1 à 3. On peut constater que le comportement de ces paramètres est suffisamment continu pour être aisément introduit sous forme de table, ou alternativement par une corrélation simple en fonction de la pression. Par exemple, dans le cas présent pour ADo: ADo = 1.406231 ( 1 - LN(P/Psat) [1 + 0.038 x LN(P/Psat)]) Au stockage pour utilisation ultérieure de l'évolution des valeurs de ces coefficients en fonction des abscisses de pression, température, et indicateur compositionnel, il peut être pratique d'ajouter l'évolution des valeurs de masse volumique de l'huile dans les conditions thermodynamiques de chaque étape du chemin thermodynamique suivi, tant au-dessus de la pression de saturation que au- dessous. On peut choisir également ou alternativement de stocker d'autres paramètres des phases huile et/ou gaz, mais ce n'est pas indispensable. f) On traduit, d'une façon connue en soi, le comportement thermodynamique de la (des) phase(s) hydrocarbure(s) en une ou plusieurs représentation(s) thermodynamique(s) de type B.O. dans les données d'entrée de la modélisation. Ces dépendances sont décrites au-dessus de la pression de saturation du mélange d'hydrocarbures et au-dessous de la pression de saturation, soit point par point dans des tables avec une méthode d'interpolation et, le cas échéant, une méthode d'extrapolation, soit par des corrélations analytiques, soit par un mélange de tables et de corrélations. g) On réalise, de façon connue en soi, une simulation de type B.O., ladite simulation permettant de calculer au moins dans chaque maille (j) et à des pas de m temps consécutifs (m, m+1, etc.), la pression (pj), la température (T7) (si celle-ci varie), les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire de mailles (j,h), les débits m Il volumiques des phases liquide (uojh) et vapeur (Ugih). h) On détermine, à chaque pas de temps (m) et dans chaque maille (j), m les constantes d'équilibre (kij) du composant i à partir des dépendances des coefficients (ADo, AD...ADp...ADn) de l'équation (1), ou des coefficients (Do, Dl,.,Dp,..Dn) de l'équation (3) écrite par phase, établies en e) avec la pression (pj), la température (Tj) et avec, le cas échéant en cas de modélisation B.O. non stricte un indicateur de composition comme le contenu en condensat dissous dans la phase gaz. i) On détermine, à chaque pas de temps (m), la fraction vaporisée (0j) dans chaque maille (j), soit à partir des résultats de la simulation de type B. O. soit, pour une meilleure précision en cas de coexistence des phases huile et gaz, 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.; "Procedure for Use of Electronic Digital Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium", J. Pet. Technol., 1952, 14, 19, à partir des fractions molaires de chaque composant i dans la composition détaillée m globale (zij) du fluide hydrocarbure dans la maille (j) au pas de temps (m) m m Nrd zii (kii-1) (4) m m- 0 i= + (kij -1) oj j) On estime la composition détaillée de chaque phase hydrocarbure, à chaque pas de temps (m) et dans chaque maille (j), en utilisant les relations suivantes m m (5) Xm Z1ii (kii - 1) Pur laphasehuile = 1 + (kj -1) ou m m m m kii Zii (kii - 1) (6) Yij - + (kj -1) O pour la phase gaz k) On estime la densité molaire de la phase huile,oj dans la maille (j) au pas de temps (m) par exemple à partir de la masse volumique de l'huile poj, dans la maille (j) au pas de temps (m), par la relation: m (7) oj = Mo, MMo étant la masse molaire de la phase liquide hydrocarbure (dite 'huile'). m Dans l'équation (7), la masse volumique poj peut donc être issue des résultats au pas de temps (m) de la simulation B.O. ou obtenues à partir des données d'entrée spécifiques à l'étape d'éclatement stockées lors de la simulation thermodynamique préalable. MMo peut être calculée à partir des résultats de la simulation B.O. si la masse molaire du constituant hydrocarbure lourd est définie par défaut ou dans les données d'entrées utilisateur, ou à partir de: Nid (8) MMo7 = x MM x,étant la masse molaire du composant ou i=l m pseudo-composant i, et les quantités xij ayant été calculées dans l'étape précédente. 1) On estime la densité molaire de la phase gaz ogj dans chaque maille (j) à chaque pas de temps (m), soit à partir de données d'entrée spécifiques à l'étape d'éclatement stockées lors de la simulation thermodynamique préalable, soit à partir m m des résultats de la simulation B.O. en pression Pj, température Tj, par exemple avec, lorsque le facteur de compressibilité du gaz fait partie des données d'entrée de la représentation thermodynamique de type B.O. la relation suivante: m (9), gj = Zgj R T, Zgj étant le facteur de compressibilité du gaz, R la constante des gaz parfaits ou, par exemple, si c'est le facteur de volume du gaz qui fait partie des données d'entrée de la représentation thermodynamique de type B.O. avec la relation suivante: m 1 (10) g j = m, Vmolst étant le volume molaire du gaz en conditions Bgj Vmolst standard et BgJ étant le facteur de volume du gaz. m) On évalue dans chaque maille (j) à chaque pas de temps (m) les débits molaires des phases liquide (Uoh) et vapeur (u.h) à partir des débits volumiques des phases liquide (uojh) et vapeur (ugh) obtenus lors de la simulation B.O. et les densités molaires obtenues dans les étapes k) et 1). n) On estime au pas de temps m+l la fraction molaire de chaque composant ou pseudo composant i dans la composition détaillée globale (z I) du composant ou pseudo composant i dans la composition détaillée globale (z ij) du fluide hydrocarbure dans la maille (j) connaissant la fraction molaire (zij) au pas de temps m à partir des équations ci-après: Zij Nj At (yij' Sg + xij, Soj)- At A (Yij' ugjh + xj uojh) h+1 J(j) (11)zi = m+l mI m m (12)NJ = NJ - At (S + Soj)- At A (Ugh + uOjh) h J(j) m l'équation (12) traduisant le bilan molaire total hydrocarbure sur la maille (j), NJ, étant la quantité molaire hydrocarbure contenue dans la maille (j) au pas de temps (m) en tenant compte des échanges de matière avec toutes les mailles (h) voisines de (j) qui forment l'ensemble J(j). m Dans l'équation (11), l'écriture des termes (yij) et (xij), - dans lesquels 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 -, suppose implicitement l'utilisation d'un schéma simple amont pour les flux compositionnels. Une écriture plus générale de ces termes est (yijh) et (Xijh) o xijh et Yih décrivent les compositions des phases liquide et gaz, obtenues d'une façon connue en soi, s'écoulant entre les mailles (j) et (h). Les termes source (débits molaires) (Sm) et (Somj) sont obtenus à partir des termes source (débits volumiques) (Si) et (Soj) dans la simulation B.O. en reproduisant pour les termes source les étapes j) à m). Les étapes h) à n) décrivent l'opération d'éclatement ou "delumping". Les détails fournis sont basés sur des équations de conservation en moles; on aurait pu tout aussi bien se baser sur des équations de conservation en masse. Dès lors que l'on sait décrire la composition détaillée dans chaque maille du modèle de simulation de type B.O. à l'instant t, qui peut être en particulier l'instant initial, on sait grâce à la présente méthode, 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 de type " Black Oil >". REVENDICATIONS 1) Méthode de modélisation permettant de prédire en fonction du temps la composition détaillée des fluides produits par un gisement souterrain en cours de production, dans au moins une zone thermodynamique définie dans la modélisation, caractérisée en ce que: - on discrétise le gisement par un maillage, chaque maille contenant une ou plusieurs phases dont au moins une phase non aqueuse; - on détermine la variation de paramètres thermodynamiques des phases non aqueuses nécessaires à une modélisation de type " Black Oil " (B.O.) au cours d'étapes d'un chemin thermodynamique suivi par les fluides dans le gisement; - on définit les fluides par une représentation détaillée à Nrd composants et/ou pseudo- composants; - on détermine au moins une fonction d'état à n paramètres permettant de simuler au moins dans la phase de préparation des données d'entrée, le comportement thermodynamique des fluides au cours des étapes du chemin thermodynamique suivi de façon à générer, dans chaque dite zone thermodynamique o l'on souhaite effectuer une opération d'éclatement, des données d'entrée supplémentaires; - on traduit le comportement thermodynamique de chaque phase non aqueuse en des données d'entrée convenant pour une représentation thermodynamique de type B.O., données d'entrée que l'on complète par les dites données d'entrée supplémentaires convenant à une opération d'éclatement; - on réalise une modélisation de type B.O., permettant de déterminer dans chaque maille et à des pas de temps consécutifs, des caractéristiques thermodynamiques de chaque phase non aqueuse et des données représentatives de déplacement des phases dans le gisement; et - on réalise une opération d'éclatement pour obtenir la composition détaillée des fluides dans chacune des mailles. 2) Méthode selon la revendication 1, caractérisée en ce que l'on caractérise le chemin thermodynamique en déterminant les variations de la composition des fluides et des propriétés des différentes phases. 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on détermine dans chaque maille au cours de la simulation de type B.O., des paramètres thermodynamiques tels que la pression (pj), la température (Tj) , si celle-ci varie, les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire de mailles(j,h), les débits m m volumiques des phases liquide (uojh) et vapeur (Ugih). 4) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'opération d'éclatement comporte une détermination des constantes d'équilibre à partir de données d'entrées spécifiques à l'opération d'éclatement, et une conversion de résultats exprimés en volume par des résultats utilisables dans des équations de conservation des quantités molaires ou massiques. ) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on utilise chaque fonction d'état pour générer n+1 relations fonctionnelles supplémentaires sous forme de tables de données ou de corrélations, que l'on inclut dans les données d'entrée.