FR 2775094 A1 19990820 FR 9801970 A 19980218 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, la composition détaillée d'un fluide produit par un réservoir et plus particulièrement la 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étenniner 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 des pseudo-composants regroupant l'azote et le méthane, des pseudo-composants 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 à 8 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. 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 peuvent être calculées en utilisant une des équations d'état bien connues des spécialistes et qui ne seront donc pas rappelées. 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 20 à 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 la 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. 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 (el). les constantes d'équilibre liquide-vapeur (k") ) de chaque composant (i), les débits d'injection ou de production (sijm) et pour chaque paire de mailles (j, h) les débits des phases liquide (um,h) et vapeur (U,3,h) 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 détaillé de la composition Q du fluide, à partir des valeurs de la fraction vaporisée (#jm) et des constantes d'équilibre (ksi) ) 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 même 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 à 8, défais 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. Le fluide regroupé à N composants et pseudo-composants est décrit par une équation d'état qui pourrait être par exemple celle de PENG ROBINS ON adaptée à ladite composition regroupée. Au début de l'exploitation du gisement, on connaît, par des mesures effectuées au préalable, pour chaque maille j du réseau, les variables principales à l'instant t=O ou pas de temps m=O. Ces variables principales sont: la quantité de fluide: fJ la fraction molaire du composant i: z je avec i zizi 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 0, la fraction molaire de chaque composant i dans la phase huile (x" avec i = 1...N) et la fraction molaire dudit composant i dans la phase gaz (yt avec i = Ces diverses fractions sont déterminées avec le système d'équations suivant: dans lesquelles bi est la fugacité du composant i. Selon une première caractéristique de l'invention, on mémorise la fraction vaporisée 0T? et les constantes d'équilibre avec i variant entre 1 et N, afin de les utiliser dans une autre étape dénommée "delumping" qui sera explicitée ultérieurement. Dans une deuxième étape et pour chaque maille, on évalue certaines des propriétés de chaque phase liquide et gaz, à savoir masse volumique (pOm, et :) z 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 ml, sont: Pour résoudre ces équations, les coefficients à calculer sont: - la compressibilité totale CTS (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), a,m est le volume molaire partiel du composant i dans la maille j au pas de temps m, 7i,mJ et h,m sont les mobilités des phases huile (o) et gaz (g), t,, est la transmissivité entre les centres des mailles j et h, I é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. p=B 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. 1 étant la maille amont avec l=j ou h. Selon une autre 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+ 1 à 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 à l'étape 1. Pour ce faire, on applique les lois de conservation défPmie par les équations suivantes Par un processus itératif, on répète les étapes ci-dessus pour calculer les variables principales aux temps t+2At, t+3At, etc... ou encore aux pas de temps m+l, m+2, m+3, etc..., l'incrément At pouvant être constant ou variable. 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 20 à 30. On connaît la composition du fluide détaillé à Q composants et/ou pseudo-composants au pas de temps m=O. 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. Comme l'écoulement peut s'effectuer dans n'importe quelle direction entre les mailles voisines, chaque maille ne peut pas être traitée indépendamment. C'est pourquoi, à chaque pas de temps, toutes les mailles sont traitées ensemble. Pour réaliser le delumping, on fait 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 du fluide, #jm=#jm - le nombre de moles de liquide et de vapeur s'écoulant entre les mailles et dans ou hors des puits est indépendant du nombre de constituants du fluide, UoJh = Uouh et Ugh = Sojm=sojm et Sgjm=sgjm 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é. Le calcul des constantes d'équilibre peut être effectué à partir de l'expression 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 ROBIN SON, l'équation généralisée de LEIBOVICI peut être écrite pour le fluide détaillé de la manière suivante: 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 calculés à partir des constantes d'équilibre pour le fluide regroupé (kt) en minimisant la fonction: De ce fait, lorsqu'on connaît la composition détaillée ZIJ 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 à l 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 zlmJ'' au pas de temps consécutif m+l à partir des équations ci-après dans lesquelles j'=j pour un débit positif de la maille j vers la maille h ou dans le puits et j'=h pour un débit négatif, j' correspondant au fluide injecté pour les puits d'injection où S sera négatif. 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. REVENDICATIONS 1. Méthode de simulation pour prédire, en fonction du temps, la 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 au moins pour chaque maille (j) et à des pas de temps consécutifs (m, m+l, ...) la fraction vaporisé (û, ), les constantes d'équilibre liquide-vapeur (k;) de chaque composant (i), les débits d'injection ou de production (sp) et pour chaque paire de mailles (j, h) les débits des phases liquide (Uojh) et vapeur (UéJh) ) 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 détaillé de la composition Q du fluide, à partir des valeurs de la fraction vaporisée (0T?) et des constantes d'équilibre (kit) ) 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). 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.