FR 3004270 A1 20141010 FR 1353066 A 20130405 Domaine technique La présente divulgation s'inscrit dans les domaines de l'exploration et de l'exploitation des réservoirs de fluides compressibles, notamment des hydrocarbures. Le réservoir peut être un gisement naturel en milieu poreux, un gisement de gaz de schiste (« shale gas ») ou un stockage artificiel comme par 10 exemple des stockages en aquifères ou en gisements déplétés. L'exploitation d'un tel réservoir nécessite d'être capable de prévoir son fonctionnement et ses performances futures, en termes de volumes ou de débits. Or les connaissances de la structure du réservoir et de son comportement vis-à-vis du fluide sont limitées. Elles proviennent, dans un premier temps, de forages 15 réalisés en nombre limité et de données d'imagerie sismique, qui n'offrent qu'une connaissance partielle du réservoir. Dans un deuxième temps, elles proviennent de mesures de production réalisées par des puits d'exploitation ou de surveillance du réservoir, telles le débit de fluide, la pression du réservoir, la production d'eau ou la position des interfaces des fluides. 20 Pour appréhender et modéliser le réservoir, les géologues peuvent utiliser des cartes de propriétés, représentant par exemple dans un plan horizontal une propriété telle que la profondeur du réservoir, son épaisseur, sa porosité ou sa perméabilité. Chaque carte présente alors une unique propriété, représentée par une intensité de gris ou un code couleur. Pour rendre compte des différents 25 aspects du réservoir, un modèle composé de plusieurs cartes représentant chacune une propriété locale différente, ainsi que d'autres paramètres représentant des propriétés globales sont utilisés. Dans le cas d'un réservoir à symétrie de révolution autour d'un axe vertical - on parle de géométrie radiale - il est possible d'utiliser des graphiques présentant, en 30 fonction de la coordonnée radiale, plusieurs propriétés, comme par exemple l'épaisseur, la porosité et la perméabilité. Ces informations peuvent être représentées par des courbes, ou des codes couleur ou des intensités de gris. Dans les deux cas, que l'on utilise un graphique dans le cadre d'une modélisation radiale, ou une carte à deux dimensions, la ou les représentations graphiques utilisées, constituent, avec éventuellement d'autres paramètres numériques non figurés sur la représentation graphique, un modèle du réservoir. Le modèle de réservoir doit faire l'objet d'un processus de calage au vu des mesures de production. Cela implique de modifier les paramètres numériques ainsi que les cartes ou graphiques, en partant de réalisations envisageables, mais naturellement, plus le modèle est sophistiqué, plus le nombre d'éléments à optimiser est élevé, et plus délicate est l'opération de calage, à la fois à cause du temps et des moyens de calcul nécessaires et des risques d'absence de convergence ou de convergence vers une solution erronée. On connait par exemple une méthode, décrite dans le document W02004/086280, pour former un modèle par une sélection de réalisations et des déformations graduelles. Cette méthode ne permet pas, aux différentes étapes de sélection, un contrôle du géologue, et est donc susceptible, dans certaines situations, de mener à des solutions non réalistes. Définition de l'invention et avantages associés Pour surmonter ces difficultés et obtenir des avantages qui apparaîtront dans la présente divulgation, il est proposé une méthode et un outil de détermination d'un modèle calé pour un réservoir souterrain de fluide, basés sur, - la génération d'une bibliothèque de représentations graphiques de caractéristiques du réservoir pour des réalisations différentes du réservoir, - la classification topologique auto-organisatrice des modèles graphiques de la bibliothèque dans un espace de classes, - et l'optimisation d'une fonction objectif au moins dans ledit espace de classes, pour déterminer une réalisation du réservoir calée sur des mesures d'exploitation. Grâce à ces caractéristiques, on permet un contrôle par le praticien de la validité du modèle calé au cours du processus, ainsi qu'une amélioration de la convergence l'optimisation étant menée sur l'espace des classes, qui est de plus faible dimension que celui des modèles ou des représentations graphiques, tout en offrant la possibilité au praticien d'intégrer dans le processus des représentations graphiques d'origine variées. Eventuellement la méthode comprend de plus une étape d'extraction d'un attribut pour des représentations graphiques de la bibliothèque, la classification étant faite sur la base d'au moins l'attribut. Les principes divulgués vont maintenant être décrits en relation avec des figures dont la liste est donnée ci-après. Liste des figures La figure 1 représente un stockage de gaz naturel. 10 La figure 2 présente une carte telle qu'utilisée dans un exemple de réalisation. La figure 3 présente les principales étapes d'un exemple. La figure 4 présente un aspect d'un autre exemple. La figure 5 présente un graphique tel qu'utilisé dans un autre exemple. La figure 6 présente un aspect d'un exemple utilisant le graphique de la figure 5. 15 La figure 7 présente un autre exemple. La figure 8 présente encore un autre exemple. Description détaillée de 'Invention En figure 1, on a représenté un stockage constituant un réservoir auquel on 20 s'intéresse dans la présente divulgation. Il comprend un réservoir 1, et on a figuré le puits de gisement 2, qui permet, via un manifold, de mesurer une pression de gisement, ainsi que le puits d'exploitation 3, par lequel circule un fluide, entrant ou sortant. La structure interne et les propriétés internes du réservoir ne sont pas connues, puisque seules quelques informations, comme évoqué plus haut, sont 25 disponibles pour le géologue. En figure 2, on a représenté une représentation graphique, ici une carte 10 constituant, avec quelques paramètres P1, P2, P3 et P4, une réalisation possible du réservoir 1. La carte 10 est un diagramme représentant un plan x,y de coupe horizontal, et figure une propriété du réservoir, par un code couleur appliqué en 30 chaque point. La propriété, dans l'exemple choisi est la porosité, mais la perméabilité, la hauteur, le faciès géologique ou la profondeur pourraient, par exemple être choisis. La carte est ici générée de manière automatique, par ordinateur, sans intervention d'un spécialiste géologue, pour avoir un niveau minimal de compatibilité avec les connaissances géologiques CG disponibles au sujet du réservoir 1. Elle est obtenue par simulation géostatistique, par modélisation stochastique non géostatistique, ou par interprétation déterministe. En figure 3, on a figuré un schéma représentant les étapes d'un exemple de mise en oeuvre de l'invention. Sur la base de connaissances géologiques CG d'un réservoir, une bibliothèque 100 formée d'un grand nombre de cartes du réservoir est générée, par ordinateur au cours d'une étape S100. Au sein de la bibliothèque, les cartes sont obtenues, toutes, par simulations géostatistiques, toutes par modélisation stochastique non géostatistique ou toutes par interprétation déterministe ou encore, certaines par une première de ces méthodes et d'autres par une deuxième de ces méthodes, et éventuellement encore certaines par la troisième de ces familles de méthodes. Chaque carte de la bibliothèque 100 présente une même propriété du réservoir en chaque point. La propriété, dans l'exemple choisi, est la porosité. Chaque carte diffère des autres cartes de la bibliothèque, car elles ont été générées en utilisant des simulations géostatistiques distinctes. Elles respectent néanmoins toutes, a minima, les connaissances géologiques CG disponibles pour le réservoir. Chaque carte peut être rapportée à un ou plusieurs attributs calculés à partir de la carte ou extraits de celle-ci. Ainsi, pour des cartes de profondeur, la surface, la forme ou le volume de la surface fermée, extraits de la carte, peuvent être utilisés comme attribut pour la réalisation. Pour d'autres cartes ou pour les mêmes cartes de profondeur, un vecteur constitué des valeurs de la propriété représentée aux noeuds d'une grille peut être utilisé comme attribut de carte. Au cours d'une étape 5200, des classes de cartes sont définies et les réalisations de la bibliothèque 100 sont classées dans ces classes. Les classes sont repérées par deux indices qui sont des nombres entiers i et j, et on peut donc représenter l'ensemble des classes, organisé, sous la forme d'un tableau à deux entrées, de dimensions n et m. Toutes les cartes sont reliées à une classe, ce qui, sur la figure est représenté par leur placement dans des cases ordonnées en deux dimensions figurant les classes. L'étape de classification S200 est effectuée avec une technique de classification topologique auto-organisatrice telle que les cartes auto-organisatrices de Kohonen (SOM pour « Self organizing map » en anglais) ou GTM (« Generative Topographic Maping » en anglais). Grâce à cette méthode, une distance peut être calculée entre chaque classe, sur la base, dans le cas présenté ici, des indices i et j. Les classes sont organisées sur au moins une dimension, et ici, elles le sont sur deux dimensions. Elles peuvent l'être sur plus de deux dimensions. Dans chaque classe, une carte représentative est déterminée. Il peut s'agir de l'une des cartes classées dans la classe, par exemple celle qui est la plus proche d'un barycentre des cartes de la classe, ou d'une carte déterminée sur la base des cartes de la classe, comme par exemple une carte moyenne, calculée sur l'ensemble des cartes de la classe, ou sur une sélection de celles-ci, les cartes les plus éloignées du barycentre étant écartées pour le calcul de la carte moyenne. A ce stade, une visualisation à l'écran des cartes représentatives pour les n fois m classes permet à un géologue d'écarter des classes qui ne seraient pas réalistes au vu de l'expérience du praticien, portant sur les propriétés courantes des réservoirs, ou les propriétés particulières du réservoir étudié. Les classes écartées sont remplacées à l'aide d'une nouvelle application de la méthode de classification S200, sur une bibliothèque comprenant toutes les cartes initiales, à l'exception des cartes de la classe supprimée. Les classes sont donc refondées, les cartes non éliminées redistribuées dans des nouvelles classes et une nouvelle visualisation à l'écran permet au géologue de juger de la pertinence des cartes représentatives des classes qui lui sont proposées, jusqu'à ce que chacune d'entre elles soit raisonnablement réaliste et puisse être conservée pour l'étape suivante. Au cours d'une étape S300, une fonction objectif, calculée pour chaque classe à l'aide de la carte représentative correspondante est optimisée dans l'espace des classes, c'est-à-dire l'espace des couples i,j. La méthode d'optimisation peut utiliser comme variable les coordonnées entières et discrètes de classes dans l'espace i,j. Mais dans une variante, il est prévu de déterminer des cartes intermédiaires entre plusieurs classes voisines de l'espace i,j, de manière continue, les cartes intermédiaires étant obtenues par interpolation des cartes voisines avec différentes pondérations, choisies pour optimiser au final la fonction objectif. La fonction objectif est construite de manière à mesurer la qualité du calage, c'est-à-dire la capacité de la carte choisie, qui est appelée à constituer le modèle ou un 5 élément du modèle de réservoir, à rendre compte des mesures d'exploitation ME, comme la pression dans les puits en fonction du temps ou le débit d'eau aux puits en fonction du temps. La fonction objectif est optimisée dans un processus mis en oeuvre par un solveur, qui prend notamment en compte la qualité de reproduction des bilans de matière constatés au cours de l'exploitation du réservoir. Différentes 10 méthodes d'optimisation numérique peuvent être mises en oeuvre. Dans une variante, la fonction objectif est calculée, non plus seulement sur les variables i,j mais aussi sur les valeurs qui peuvent être prises par un ou plusieurs paramètres P1, P2, P3 ou P4 tels qu'évoqués en relation avec la figure 2. Cette variante est représentée en figure 4, où la fonction objectif est calculée sur la 15 base de 15 paramètres qui sont des nombres réels indépendants des cartes c'est à dire non dépendants directement des cartes, et des coordonnées i,j des classes de cartes organisées comme présenté plus haut. La fonction objectif est la minimisation de la somme de deux fonctions, à savoir Dpg et Dzrio. Ainsi, la fonction objectif est minimisée sur 17 paramètres. 20 En figure 5, on a figuré la représentation graphique d'un modèle radial, comme évoqué plus haut. Il s'agit d'un graphique, où l'axe des abscisses présente la coordonnée radiale autour d'un axe vertical, le réservoir étant considéré posséder une géométrie de révolution autour de cet axe. Diverses informations sont présentées sur ce graphique, par des courbes et des codes couleur. 25 On peut utiliser des cartes en coordonnées polaires (rayon azimuth). En figure 6, on a représenté une classification d'une bibliothèque de graphiques tels que présentés en figure 5. On ne présente ici, par simplification, pour chaque classe, qu'un seul graphique représentatif, choisi ou calculé comme présenté plus haut pour les cartes représentatives de classes. Mais comme précédemment, 30 chaque classe comprend de nombreux graphiques. Les classes sont organisées suivant deux axes, et on peut calculer une distance entre deux classes, à l'aide des coordonnées de la classe sur les deux axes. En figure 7, on a représenté un aspect d'une variante de mise en oeuvre impliquant plusieurs cartes de propriétés comme celle représentée en figure 2. Dans cet exemple, on génère indépendamment plusieurs bibliothèques de cartes, les cartes de chaque bibliothèque présentant une propriété donnée. On effectue, indépendamment, des classifications de chaque bibliothèque ordonnées sur deux dimensions puis le géologue peut visualiser les cartes représentatives des classes, présentées sur deux dimensions. Ainsi, on a présenté en figure 7, les cartes représentatives des classes ordonnées pour une bibliothèque de cartes de la propriété HT (hauteur totale), les cartes représentatives des classes ordonnées pour une bibliothèque de cartes de la propriété HU (hauteur utile), et les cartes représentatives des classes ordonnées pour une bibliothèque de cartes de la propriété PHIU (porosité utile). Qui plus est, deux zones ou niveaux de réservoir (ou profondeurs de réservoir mesurées à partir de la surface du sol) sont étudiés dans le cas présenté en figure 6. Ces deux niveaux de réservoir sont notés R1B et R1C. Deux bibliothèques sont donc préparées pour chaque propriété, les cartes de la première bibliothèque figurant la propriété à la profondeur du niveau R1B, et les cartes de la deuxième bibliothèque figurant la propriété à la profondeur R1C. Le nombre de profondeurs peut naturellement être plus grand que deux, et il y a dans ce cas autant de bibliothèques par propriété que de profondeurs étudiées. Chacune des bibliothèques fait l'objet, indépendamment des autres, d'une classification organisée sur deux axes i,j, et les cartes représentatives des classes peuvent être visualisées par un géologue pour éliminer les classes non réalistes. Une fonction objectif prenant en variables les coordonnées i,j de chacune des classifications de bibliothèques de cartes, soit en figure 6, six couples i,j, est ensuite optimisée pour déterminer des cartes, extraites ou calculées à partie des cartes représentatives des classes pour chaque classification de bibliothèque. La réalisation de réservoir retenue pour modéliser celui-ci est ensuite une combinaison des cartes retenues. La combinaison des cartes de deux niveaux réservoir différent est une juxtaposition, alors que la combinaison des cartes de deux propriétés différentes pour un même niveau de réservoir est une superposition des cartes. Comme évoqué précédemment, des paramètres numériques P1, P2, P3 et P4 non représentés sur des cartes peuvent aussi être utilisés comme variables de la fonction objectif, et celle-ci peut être minimisée en faisant varier non seulement les couples i,j de coordonnées de classes dans les classifications de cartes, mais aussi les paramètres P1, P2... P4. Ainsi, on a représenté en figure 8, le processus d'optimisation utilisant six bibliothèques de cartes classifiées et 15 paramètres hors des cartes, ce qui définit un ensemble de 21 paramètres à faire varier pour optimiser la fonction objectif. L'optimisation de la fonction objectif est faite, dans une variante, en imposant une corrélation minimale entre au moins deux propriétés différentes du réservoir, par exemple les propriétés HU et PHIU de la figure 7. L'invention ne se limite pas aux modes de réalisation décrits mais s'étend à toutes les variantes dans le cadre de la portée des revendications. REVENDICATIONS1. Méthode de détermination d'un modèle calé pour un réservoir souterrain de fluide comprenant, 10 - une étape de génération (S100) d'une bibliothèque (100) de représentations graphiques décrivant le réservoir pour des réalisations différentes de celui-ci, - une étape de classification auto-organisatrice (S200) des représentations graphiques de la bibliothèque dans un espace (200) de classes organisées 15 sur au moins une dimension, - et une étape d'optimisation (S300) d'une fonction objectif au moins dans ledit espace de classes, pour déterminer une réalisation (300) du réservoir calée sur des mesures d'exploitation (ME). 2. Méthode de détermination selon la revendication 1, comprenant de plus 20 une étape d'extraction d'un attribut pour des représentations graphiques de la bibliothèque, la classification étant faite sur la base d'au moins l'attribut. 3. Méthode de détermination selon la revendication 1 ou la revendication 2, dans laquelle l'optimisation est faite en calculant la fonction objectif pour une classe en utilisant une représentation graphique représentative de la 25 classe. 4. Méthode de détermination selon l'une des revendications 1 à 3, dans laquelle, entre l'étape de classification et l'étape d'optimisation, la possibilité d'éliminer une ou plusieurs classes apparaissant comme non réalistes par observation visuelle d'au moins une représentation graphique 30 représentative de chaque classe est offerte à un opérateur. 5. Méthode de détermination selon l'une des revendications 1 à 4, dans laquelle la classification est effectuée sur deux dimensions ou sur plus de deux dimensions. 6. Méthode de détermination selon l'une des revendications 1 à 5, dans laquelle la fonction objectif est optimisée sur un ou plusieurs paramètres qui sont des nombres réels non dépendants directement des cartes, et sur de coordonnées des classes de cartes sur la au moins une dimension. 7. Méthode de détermination selon l'une des revendications 1 à 6, dans laquelle plusieurs bibliothèques de représentations graphiques sont générées, les représentations graphiques de chaque bibliothèque faisant l'objet, indépendamment, de classifications auto-organisatrices dans des espaces de classes respectifs, et la fonction objectif étant optimisée dans des espaces de classes respectifs, une combinaison de représentations graphiques représentatives des classes optimisant la fonction étant utilisée pour la détermination de la réalisation calée. 8. Méthode de détermination selon la revendication 7, dans laquelle des bibliothèques de représentations graphiques portant sur des zones distinctes du réservoir sont générées, la fonction objectif étant optimisée dans les espaces de classes comprenant lesdites bibliothèques. 9. Méthode de détermination selon la revendication 7 ou la revendication 8, dans laquelle des bibliothèques de représentations graphiques portant sur des propriétés distinctes du réservoir sont générées, la fonction objectif étant optimisée dans les espaces de classes comprenant lesdites bibliothèques. 10. Méthode de détermination selon la revendication 9, dans laquelle l'optimisation de la fonction objectif est faite en imposant une corrélation minimale entre au moins deux propriétés différentes du réservoir décrites dans deux bibliothèques différentes. 11.Méthode de détermination selon l'une des revendications 1 à 10, dans laquelle les représentations graphiques comprennent des cartes représentant une propriété dans un plan du réservoir 12.Méthode de détermination selon l'une des revendications 1 à 11, dans laquelle les représentations graphiques comprennent des représentations d'au moins une propriété le long d'un rayon du réservoir, celui-ci étant considéré avoir une symétrie de révolution autour d'un axe vertical. 13.Méthode de détermination selon l'une des revendications 1 à 12, dans laquelle, l'optimisation est faite sur des coordonnées discrètes des classes. 14. Méthode de détermination selon l'une des revendications 1 à 13, dans laquelle l'optimisation est faite sur des coordonnées continues dans l'espace des classes, en interpolant des représentations graphiques représentatives des classes. 15. Méthode de détermination selon la revendication 1 à 14, avec une bibliothèque comprenant au moins des cartes obtenues par simulations géostatistiques, des cartes obtenues par modélisation stochastique non géostatistique, ou au moins des cartes obtenues par interprétation déterministe, ou au moins des cartes obtenues par une première de ces méthodes et des cartes obtenues par une deuxième de ces méthodes, ou au moins des cartes obtenues par respectivement chacune des trois méthodes. 16.Méthode de détermination selon l'une des revendications 1 à 15, la propriété comprenant la porosité, la perméabilité, la hauteur, la profondeur ou le faciès géologique.