FR 2857764 A1 20050121 FR 0308841 A 20030718 La présente invention concerne une méthode de déformation graduelle d'un modèle booléen de simulation d'un milieu hétérogène, contraint par des données dynamiques. La méthode selon l'invention s'applique notamment à la construction d'un modèle servant à représenter la distribution spatiale dans une zone souterraine de certaines propriétés physiques du milieu comme la perméabilité, la porosité, etc. Etat de la technique Les documents suivants que l'on va citer dans le cours de la description ci-après, illustrent l'état de la technique, : - Chilès, J. P., and Delfiner, P., 1999, Geostatistics : spatial uncertainty, Wiley, New York, 695p. - Hu, L. -Y., 2000a, Graduai déformation and itérative calibration of Gaussian-related stochastic models, Math. Geol., 32(1). Hu, L.-Y., 2000b, Graduai déformation of non-Gaussian stochastic models, Geostats 2000 Cape Town, WJ Kleingeld and DG Krige (eds.), 1, 94-103. - Hu, L.-Y., 2003, History matching of object-based stochastic réservoir models. SPE 81503. - Journel, A., and Huijbregts, C.J., 1978, Mining Geostatistics, Academic Press, London, 600p. - Le Ravalée, M., Noetinger, B., and Hu, L. -Y., 2000, The FFT moving average (FFT- MA) generator : an efficient numerical method for generating and conditioning Gaussian simulations, Math. Geol., 32 (6), 701-723. - Le Ravalec, M., Hu, L. -Y. and Noetinger, B., 2001, Stochastic Réservoir modeling constrained to dynamic data : local calibration and inference of the structural parameters, SPE Journal, 25-31. Différentes méthodes basées sur un schéma de déformation graduelle, ont fait l'objet de brevets ou demandes de brevet suivants du demandeur : 2. 780.798, FR 2.795.841, FR 2. 821.946, FR 02/13. 632 ou FR 03/02.199. Les réalisations de modèles stochastiques gaussiens, ou apparentés, sont souvent utilisées pour représenter la distribution spatiale de certaines propriétés physiques, comme la perméabilité ou la porosité, dans les réservoirs souterrains. On recourt alors couramment aux méthodes inverses pour contraindre ces réalisations à des données dont elles dépendent de façon non linéaire. C'est le cas notamment en hydrologie ou en ingénierie pétrolière. Ces méthodes reposent sur la minimisation d'une fonction objectif, encore appelée fonction coût, qui mesure l'écart entre les données mesurées sur le terrain et les réponses correspondantes simulées numériquement pour les réalisations représentant le milieu que l'on cherche à caractériser. Le but affiché est d'identifier les réalisations associées aux valeurs de fonctions objectifs les plus faibles, c'est-à-dire les réalisations les plus cohérentes au regard des données. La méthode de déformation graduelle a été introduite dans ce cadre (cf. Hu, 2000a; Le Ravalec et al., 2000). Cette technique de paramétrage géostatistique permet de modifier graduellement les réalisations à partir d'un nombre restreint de paramètres. Elle est particulièrement adaptée aux problèmes de minimisation, car appliquée aux réalisations, elle induit une variation continue et régulière de la fonction objectif. De fait, la minimisation peut être conduite à partir des techniques les plus avancées, à savoir les techniques de gradients. La méthode de déformation graduelle s'est avérée efficace pour contraindre des modèles de réservoirs pétroliers à des données de production (cf. Le Ravalec et al., 2001). Par la suite, la méthode de déformation graduelle, initialement mise en place pour des modèles gaussiens, a été étendue à des modèles non gaussiens (cf. Hu, 2000b), et plus particulièrement aux modèles objets ou booléens. Ces modèles sont utilisés pour décrire des milieux comprenant par exemple des chenaux ou des fractures. Les chenaux ou fractures sont alors assimilés à des objets. Un algorithme a été proposé pour simuler la migration graduelle des objets dans l'espace, c'est à dire le déplacement graduel des objets dans l'espace. Pour une réalisation d'un modèle objet, ce type de perturbation se traduit par une variation lissée de la fonction objectif, comme pour les modèles gaussiens. Cet algorithme a ensuite été généralisé au modèle booléen non-stationnaire et conditionnel (cf. Hu, 2003). Par ailleurs, toujours en s'appuyant sur la déformation graduelle des lois gaussiennes, (cf. Hu (2003) propose des solutions permettant de modifier progressivement le nombre d'objets peuplant un modèle, ce nombre étant représentatif d'une loi de Poisson. Néanmoins, une des limites de ces développements est que les objets apparaissent ou disparaissent soudainement, ce qui peut engendrer des discontinuités sévères de la fonction objectif et rendre difficilement praticables les techniques de gradients utilisées classiquement pour mener à bien le processus de minimisation. Déformation graduelle : rappelsFonction aléatoire multigaussienne Jusqu' à présent, les principes de déformation graduelle qui ont été proposés s'appliquent aux fonctions aléatoires multigaussiennes. Soient, par exemple, deux fonctions aléatoires indépendantes, Yl(x) et Y2(x), multigaussiennes et stationnaires d'ordre 2. x est le vecteur position. On suppose que ces deux fonctions ont même moyennes et variances, à savoir 0 et 1, et même fonction de covariance. On construit alors une nouvelle fonction aléatoire Y(t) en combinant Y1 et Y2 suivant l'expression: Y(t) = Y1 cos(t)+ Y2 sin(t). On peut montrer que quel que soit t, Y a même moyenne, variance et modèle de covariance que Y1 et Y2. En outre, Y(t) est aussi une fonction aléatoire multigaussienne, car elle est la somme de deux fonctions aléatoires multigaussiennes. Suivant ce principe de combinaison, on peut à partir de deux réalisations indépendantes y1 et y2 de Y1 et Y2, construire une chaîne de réalisations y(t) qui dépendent uniquement du paramètre de déformation t. L'idée de base des processus de minimisation par déformation graduelle est d'explorer cette chaîne de réalisations et de déterminer le paramètre de déformation qui fournit la réalisation la plus compatible avec les données mesurées sur le terrain, c'est-à-dire les pressions, les débits de production, les temps de percée, etc. Comme l'exploration d'une unique chaîne de réalisation ne permet pas en général d'identifier une réalisation fournissant une fonction objectif suffisamment faible, on itère le processus de recherche. On combine alors la réalisation optimale déterminée pour la 1ère chaîne avec une nouvelle réalisation indépendante de Y2 et on en déduit une nouvelle chaîne de réalisations dont l'exploration peut fournir une réalisation qui diminue davantage la fonction objectif, etc. Processus ponctuel de Poisson L'élément clé des modèles booléens est un processus ponctuel de Poisson qui caractérise l'agencement spatial des objets. Considérons un modèle booléen de base pour lequel les objets ont même forme et sont répartis aléatoirement et uniformément dans l'espace. Les positions de ces objets sont distribuées suivant un processus ponctuel de Poisson de densité constante. En d'autres termes, la position d'un objet dans l'espace à n dimensions [0,1]" est définie par le vecteur x dont les n composantes sont des nombres uniformes tirés indépendamment selon la loi de distribution uniforme entre 0 et 1. La technique de migration (cf. Hu, 2000b) des objets consiste à déformer graduellement la position des objets. La position d'un objet est déterminée par des nombres uniformes. Avant toute chose, on transforme ces nombres uniformes en nombres gaussiens : Y=G-1 (x). G est la fonction de répartition normale standard. Soient x1 la position initiale d'un objet donné et X2 une autre position possible, indépendante de x1. On définit une trajectoire pour l'objet en combinant les transformées gaussiennes de ces deux positions suivant la méthode de déformation graduelle : x) = GlG-I (XI )cos(t)+ G-I (X2 )sin(t)1 On peut montrer que pour toute valeur du paramètre de déformation t, x est un point uniforme de [0,1]". Quand les deux positions XI et X2 sont fixées, la trajectoire est complètement déterminée. Un exemple bidimensionnel est visualisé sur la Fig.1. La technique de migration des objets constitue un premier pas vers la déformation graduelle des simulations booléennes. Une de ses limites est que le nombre d'objets est supposé constant au cours de la déformation. Des solutions sont proposées pour modifier progressivement le nombre d'objets peuplant un modèle (cf. brevet N 01/03 194 ou Hu, 2003 précités). Néanmoins, une des limites de ces développements est que les objets apparaissent ou disparaissent soudainement, ce qui peut engendrer des discontinuités sévères de la fonction objectif. La méthode de déformation graduelle selon l'invention permet, on va le voir, de réduire cette discontinuité et ainsi de faciliter la mise en #uvre des technique d'optimisation basées sur le gradient. LA METHODE SELON L'INVENTION La méthode selon l'invention permet de déformer graduellement un modèle booléen permettant de simuler au mieux la configuration spatiale dans une zone souterraine hétérogène, d'objets géologiques définis par des grandeurs physiques, qui soit contraint par des données dynamiques mesurées, représentatives de déplacements de fluides dans le milieu, et en imposant une continuité dans l'évolution des dits objets. Elle comporte la réalisation d'un processus itératif d'optimisation à partir de réalisations incluant chacune des objets dont le nombre est une variable aléatoire de Poisson de moyenne définie, et une minimisation d'une fonction objectif mesurant l'écart entre des données dynamiques réelles et les données dynamiques simulées au moyen d'un simulateur d'écoulement à partir de la réalisation combinée, par ajustement de coefficients de combinaison, le processus d'ajustement itératif étant poursuivi jusqu'à l'obtention d'une réalisation optimale du modèle stochastique. La méthode se distingue en ce que, à chaque itération : - on construit une réalisation combinée obtenue par combinaison d'une part d'une réalisation initiale formée de Nl(t) objets correspondant à une première valeur moyenne et d'au moins une deuxième réalisation indépendante du même modèle formée de N2 (t) correspondant à une deuxième valeur moyenne, cette combinaison étant telle que le nombre N (t) de la combinaison ait une valeur moyenne constamment égale à la somme de la première et de la deuxième valeur moyenne. Suivant un mode de mise en #uvre, à chaque itération et pour une même valeur moyenne de la combinaison, on fait varier de façon concomitante la première et la deuxième valeur moyenne de manière à faire varier graduellement le nombre et la taille d'objets issus de chacune des réalisations combinées. En d'autres termes, la méthode proposée permet non seulement de modifier graduellement le nombre d'objets peuplant un modèle booléen mais aussi et surtout d'envisager une nouvelle technique de construction des modèles booléens avec, à présent, apparition et disparition progressives des objets. Cette nouvelle possibilité permet d'atténuer les discontinuités de la fonction objectif dans une procédure d'inversion et par conséquent de faciliter la mise en place des algorithmes d'optimisation basés sur le gradient. Le schéma proposé de déformation graduelle appliqué aux variables aléatoires de Poisson rend possible la déformation graduelle du nombre d'objets peuplant des simulations booléennes. De plus, nous proposons d'associer la taille des objets au processus de Poisson utilisé pour générer le nombre d'objets d'une simulation booléenne, ce qui nous conduit à un algorithme permettant de faire apparaître ou disparaître un objet de façon progressive lors de la déformation d'une réalisation booléenne. Ce faisant, on réduit les discontinuités éventuelles de la fonction objectif. Ces algorithmes sont aisés à mettre en #uvre et peuvent être combinés avec un algorithme d'optimisation pour caler des modèles de réservoir de type booléen à des données de production, pression, temps de percée, etc. Contrairement au modèle de réservoir de type gaussien, la discontinuité d'une fonction objectif est intrinsèque au modèle de réservoir de type objet. Dans certains cas, le déplacement d'un objet d'une seule maille peut provoquer un changement brutal de la connexité du modèle et par conséquent un changement radical du comportement hydrodynamique du modèle. En faisant apparaître ou disparaître progressivement un objet dans une simulation booléenne, nous n'éliminons pas, mais atténuons cette discontinuité. Par conséquent, nous facilitons la mise en #uvre des algorithmes d'optimisation basés sur les gradients. PRESENTATION SUCCINCTE DES FIGURES Les caractéristiques et avantages de la méthode et du dispositif selon l'invention, apparaîtront plus clairement à 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ù : - la figure 1 montre un exemple de trajectoire définie à partir de deux points dans [0,1]2 ; la figure 2 montre une simulation d'une réalisation n pour une variable de Poisson de paramètre # ; la figure 3 montre un exemple de déformation graduelle d'une variable de Poisson de moyenne À en combinant deux variables de Poisson indépendantes, N1 et N2 ; la figure 4 montre un exemple de déformation graduelle d'une réalisation de Poisson pour une variable de paramètre 10 ; la figure 5 montre un exemple de déformation graduelle du nombre et de la position des objets avec apparition et disparition progressive. La petite ellipse grise disparaît et une ellipse noire apparaît ; la figure 6 montre un exemple de modèle de référence au centre. Pression au puits injecteur à gauche et flux fractionnaire au puits producteur à droite ; la figure 7 montre un exemple d'invasion du réservoir par l'eau injectée ; la figure 8 montre l'évolution de la fonction objectif lorsque le paramètre de déformation décrit une chaîne ; la figure 9 montre un exemple de déformation graduelle du nombre et de la position des objets avec apparition et disparition soudaine ; cette chaîne correspond aux quatre points encerclés de la figure 8 ; la figure 10 montre un autre exemple de déformation graduelle du nombre et de la position des objets avec apparition et disparition progressive ; il s'agit de la même chaîne que pour la figure 9. DESCRIPTION DETAILLEE DE LA METHODE Déformation graduelle du nombre d'objets La technique de migration des objets constitue un premier pas vers la déformation graduelle des simulations booléennes. Une de ses limites est que le nombre d'objets est upposé constant au cours de la déformation. Des solutions sont proposées pour modifier )rogressivement le nombre d'objets peuplant un modèle (cf. Brevet N 01/03 194 ou Hu, @003 précités). Néanmoins, une des limites de ces développements est que les objets apparaissent ou disparaissent soudainement, ce qui peut engendrer des discontinuités sévères de la fonction objectif. La méthode de déformation graduelle selon l'invention permet, on va le voir, de réduire cette discontinuité et ainsi de faciliter la mise en #uvre des technique d'optimisation basées sur le gradient. Loi de Poisson Notre objectif est de déformer graduellement le nombre des objets pour une simulation booléenne. Ce nombre est aléatoire et suit une loi de Poisson, c'est-à-dire que la probabilité pour que ce nombre (non négatif) soit n est : où # est la moyenne et la variance de la variable de Poisson N . Une technique possible pour simuler une variable de Poisson est de générer un processus de Poisson de moyenne 1 sur une période égale à #, puisque le nombre d'événements survenant sur cette période suit une loi de Poisson de paramètre #. Pour ce faire, on génère successivement des intervalles OEI, E1E2, ...EnEn+1 (figure 2) indépendants suivant une loi exponentielle de moyenne égale à 1 (loi [gamma]1), et on ajoute leurs valeurs jusqu'à dépasser #. La fonction de répartition d'une variable X suivant la loi exponentielle de moyenne 1 s'écrit : F (x)=1-exp(-x). Si r est un nombre aléatoire, uniformément distribué entre 0 et 1, 1-r est aussi uniformément distribué sur [0 ;1]. On pose alors 1- r = 1 - exp(- x) . De fait, pour simuler les segments successifs, il suffit de répéter x=-Log(r) pour différents r. La réalisation n de la variable de Poisson N est alors le plus grand entier n tel que : ou, ce qui est équivalent, mais plus économique du point de vue calcul : Par contre, cette seconde formulation est sans doute plus numériquement instable quand # est grand. Déformation graduelle de variables de Poisson Les principes de déformation graduelle présentés plus haut s'appliquent à des nombres gaussiens qui varient continûment dans 9\. Les variables de Poisson étant des entiers, ces principes ne s'appliquent pas en l'état. On rappelle que la somme de deux variables de Poisson, indépendantes, de paramètres l'une # et l'autre , est encore une variable de Poisson, mais de paramètre # + . Ce résultat fondamental est à la base de l'algorithme de déformation graduelle que nous proposons. Soient NI et N2 deux variables de Poisson, indépendantes et de paramètre # (cf. figure 3). On propose de déformer graduellement les paramètres des deux variables additionnées, mais en respectant à chaque instant le fait que la somme des paramètres vaut #. Pour éviter toute confusion entre le paramètre des variables de Poisson et le paramètre de déformation, on parlera respectivement de moyenne et paramètre de déformation. On construit ainsi une nouvelle variable de Poisson N(t) de moyenne . à partir de : Ntlll = Nllaltl+ N21a2t1'1 2 avec La, = 1. a, (t)# est la moyenne de la variable Ni au moment de la combinaison(cf. @=1 figure 3). On choisira, par exemple, al = cos2 (t) et a2 = sin2 (t) . Ce paramétrage s'étend facilement à n dimensions : n n N(t XÀ} = L. N, la, (t)À} avec L. a, = 1 , et présente l'avantage d'être périodique. 1=1 1=1 Pour les coefficients ai, on peut choisir par exemple des fonctions trigonométriques. En variant le paramètre de déformation t, on modifie la moyenne des variables additionnées, ce qui affecte les réalisations de ces variables. Leur somme change donc aussi. En suivant ce principe de déformation, on peut à partir de deux réalisations n1 et n2 de Ni et N2, construire une chaîne de réalisations n(t) dépendant uniquement du paramètre de déformation t. Quand t égale 0, la réalisation de N vaut n1 ; quand t égale ~#/2, la réalisation de N vaut n2 . On peut montrer que N(t) est périodique de période 1[. La figure 4 illustre ce principe de déformation pour une variable de Poisson de moyenne 10. Quand le paramètre de déformation divisé par 1[ vaut 0, la réalisation est la même que la réalisation considérée initialement (nI=12). Quand il vaut 0.5, la réalisation est la même que la deuxième réalisation considérée pour la combinaison (n2=9). Une autre approche pourrait être envisagée pour déformer graduellement une variable de Poisson de paramètre #. Elle consisterait à déplacer un segment de longueur # sur une droite construite à partir de l'addition d'une infinité de segments de longueurs obtenues par une loi [gamma]1. Le nombre de segments complets rentrant dans le segment de longueur A serait une réalisation de la variable de Poisson et varierait avec la position de ce segment. Construction successive de chaînes de déformation Dans le cas des réalisations gaussiennes, la méthode de déformation graduelle s'intègre naturellement dans les processus de minimisation. L'ensemble se traduit par l'exploration successive de chaînes de réalisations gaussiennes fabriquées à partir d'une réalisation initiale qui correspond à l'optimum déterminé pour la chaîne précédente et d'une seconde réalisation gaussienne générée aléatoirement pour chaque chaîne. Pour les nombres de Poisson, on suit une démarche analogue. La déformation graduelle selon les principes précisés dans la section précédente de deux réalisations d'une loi de Poisson, une, dite initiale, et une seconde, dite complémentaire, fournit une première chaîne. L'exploration de cette chaîne conduit à l'identification d'une réalisation optimale qui minimise la fonction objectif. Cette réalisation est alors utilisée comme réalisation initiale pour la création d'une nouvelle chaîne. Cette nouvelle chaîne requiert en outre une nouvelle réalisation complémentaire, générée indépendamment. A nouveau, l'exploration de cette nouvelle chaîne peut permettre de réduire davantage la fonction objectif. Ce processus de recherche est itéré jusqu'à ce que la fonction objectif soit jugée suffisamment faible. Apparition et disparition progressive d'un objet Le principe de déformation graduelle suggéré plus haut permet de perturber au :ours d'une optimisation le nombre d'objets peuplant une simulation booléenne. clairement, l'apparition ou la disparition brutale d'objets est susceptible d'induire une variation brusque de la fonction objectif. Or, les optimisations s'appuyant sur des calculs le gradients nécessitent une évolution continue de la fonction objectif. Nous proposons lonc un nouveau type de modèles booléens assurant l'apparition et la disparition régressives des objets. Principe Revenons à la génération d'un nombre n de Poisson à partir d'un processus de 'oisson. La méthode rappelée plus haut (cf. figure 2) consiste à mettre bout à bout des segments jusqu'à l'obtention d'une longueur excédant le paramètre # de la variable de ?oisson. n est l'entier tel que OEn .. Appelons Li le point tel que la longueur du segment OL1 égale la moyenne de la variable de Poisson (cf. figure 3). Au cours du processus de déformation, ce point se iéplace. A t=0, pour la variable NI, OLi=, . Le segment OL1 comprend alors n1 segments complets (OEI, EIE2, ...Em-lEni) et un (+1)' segment tronqué (EniLi). On en déduit lue la simulation booléenne est peuplée de n1+1 objets de taille, surface ou volume définis i partir des longueurs anamorphosées des segments OEi, E1E2, ... E,,Ll. Par souci de simplification, on parlera dans la suite de taille d'objets. Si on augmente le paramètre de iéformation t, le segment OL1 se réduit ; il en va de même pour En1L1. En conséquence, la :aille du (n1+1)ème objet diminue: l'objet disparaît progressivement. Si la moyenne de la variable Ni continue de décroître, le n1ème objet va aussi commencer à se réduire, etc. Dans e même temps, la moyenne de la variable N2 augmente. La réalisation pour cette variable vaut d'abord 0, puis un premier objet apparaît progressivement. Sa taille est contrôlée par a longueur du segment PL2, où le point L2 est pour N2 l'équivalent de L1 pour Ni. Quand PL2=PI1, le premier objet est complet. Si la moyenne de N2 continue d'augmenter, un deuxième objet apparaît. Sa taille dépend à présent de la longueur du segment I1L2. La figure 5 illustre la méthode proposée. On distingue deux familles d'objets, des ellipses noires et des ellipses grises. La différence de couleur tient au fait que les objets considérés sont obtenus à partir de deux variables de Poisson, Ni et N2. La position des ellipses varie graduellement ; ellipses apparaissent et disparaissent progressivement. Sur l'image 1, on observe 4 ellipses grises dont 3 grandes (deux se superposent plus ou moins) et une petite : dernière est en train d'apparaître. On remarque aussi une ellipse noire et un point noir, qui traduit l'apparition d'une nouvelle ellipse noire. Sur les images 2 à 5, la petite ellipse grise devient de plus en plus petite jusqu'à n'être plus qu'un point. Dans le même temps, la petite ellipse noire croît. Sur l'image 6, la petite ellipse grise a totalement disparu. Une autre ellipse grise commence à présent à se réduire. La petite ellipse noire devient un peu plus grande. Comme plus haut, on pourrait simuler le même type de déformation graduelle en déplaçant un segment de longueur # sur une droite construite à partir de l'addition d'une infinité de segments de longueurs obtenues par une loi [gamma]1. Ce segment de longueur , comprendrait des segments complets au centre et, aux extrémités, un segment qui tendrait à grandir et un autre qui afficherait le comportement inverse. Impact sur lafonction objectif A titre d'exemple, nous considérons un modèle de réservoir synthétique figure 6. Il s'agit d'une roche réservoir de perméabilité 500 mD et de lentilles peu perméables de perméabilité 50 mD. Les lentilles sont représentées par des ellipses. Par simplicité, on suppose la porosité constante et égale à 30% partout. Le modèle est discrétisé sur une grille de 200x200 mailles, de longueur 1 m suivant l'axe X et 0. 8 m suivant l'axe Y. Pour ce réservoir initialement saturé en huile, on effectue le test suivant à l'aide du simulateur d'écoulement 3DSL. On injecte de l'eau à 100 m3/jour pendant 100 jours au puits I et on produit de l'huile à pression constante au puits P. La pression au puits injecteur et le flux fractionnaire en eau au puits producteur sont reportés sur la figure 6. Ces données sont appelées données de référence car elles se rapportent au modèle de référence. L'invasion en eau est illustrée par la figure 7 pour différents temps. Les lentilles forment initialement une barrière à l'écoulement, puis peu à peu se saturent en eau. Supposons à présent que l'on ignore tout de la position et du nombre des lentilles. Ne connaissant que la réponse en pression au puits injecteur et la réponse en eau au puits producteur, on cherche à identifier un modèle de réservoir aussi cohérent que possible avec ces données. On suppose que des informations issues de la géologie, par exemple, permettent d'approcher le nombre des lentilles à partir d'une loi de Poisson de paramètre 10. En appliquant les processus de déformation graduelle suggérés plus avant, on construit une chaîne de simulations booléennes caractérisées par la variation graduelle du nombre d'ellipses et leur position. On aurait pu aussi considérer la déformation graduelle de la taille des ellipses. Pour toutes ces simulations, on simule numériquement le test d'injection et on en déduit une fonction objectif mesurant l'écart entre les données de référence et les données simulées (figure 8). Le paramètre de déformation t contrôle à la fois la déformation du nombre d'ellipses et leur position. Les relations de déformation graduelle sont telles qu'on explore une période complète avec t variant entre -#/2 et +#/2. On examine deux cas. Les ellipses apparaissent et disparaissent soudainement pour le premier et progressivement pour le deuxième. Pour le premier cas, la fonction objectif évolue parfois par sauts. On vérifie que ces derniers sont absents, ou au moins très atténués, pour le deuxième. La variation du nombre d'objets en fonction du paramètre de déformation tout au long de la chaîne est décrite sur la figure 4 : nombre est celui que l'on a pour des apparitions et disparitions soudaines d'ellipses. Sur la figure 8, pour la courbe en trait épais, le nombre d'objets varie avec apparition et disparition progressive des objets. Pour les losanges, le nombre d'objets varie avec apparition et disparition soudaine des objets - il est spécifié par la figure 4. Concentrons nous sur les 4 réalisations dont les valeurs de fonction objectif sont encerclées sur la figure 8. On se trouve justement dans un cas où la fonction objectif affiche un comportement discontinu lorsque les ellipses apparaissent et disparaissent soudainement. Au contraire, l'évolution est continue pour des apparitions et disparitions progressives. Les réalisations en question sont représentées sur la figure 9 et la figure 10. Sur la figure 9 2), on observe que la disparition brusque d'une ellipse agrandit significativement le couloir d'écoulement du puits injecteur vers le puits producteur. Cette réalisation est d'ailleurs celle qui, de toute la chaîne, fournit la fonction objectif la plus faible. Ce n'est pas le cas sur la figure 10 2), même si l'ellipse concernée a diminué en taille. Cette dernière continue de se réduire sur la figure 10 3). Sur la figure 10 4), elle devient si petite que son effet est négligeable et que les scenarii 9-4 et 10-4 sont finalement équivalents. Construction successive de chaînes de déformation Comme explicité à la fin de la section 2 ci-dessus, on peut, au cours d'un processus de minimisation, explorer successivement plusieurs chaînes, sachant que la réalisation de départ pour une chaîne est la réalisation optimale identifiée pour la chaîne précédente. La seule différence par rapport à ce qui a été dit dans la section 3, est qu'on intègre aussi les segments tronqués dans la construction du segment de longueur #. REVENDICATIONS 1) Méthode pour déformer graduellement un modèle booléen permettant de simuler au mieux la configuration spatiale dans une zone souterraine hétérogène, d'objets géologiques définis par des grandeurs physiques, qui soit contraint par des données dynamiques mesurées, représentatives de déplacements de fluides dans le milieu, et en imposant une continuité dans l'évolution des dits objets, dans laquelle on réalise un processus itératif d'optimisation à partir de réalisations incluant chacune des objets dont le nombre est tiré d'une variable aléatoire de Poisson de moyenne définie, et on minimise une fonction objectif mesurant l'écart entre des données dynamiques réelles et les données dynamiques simulées au moyen d'un simulateur d'écoulement à partir de la réalisation combinée, par ajustement de coefficients de combinaison, le processus d'ajustement itératif étant poursuivi jusqu'à l'obtention d'une réalisation optimale du modèle stochastique, caractérisée en ce que, à chaque itération : - on construit une réalisation combinée obtenue par combinaison d'une part d'une réalisation initiale formée de Nl(t) objets correspondant à une première valeur moyenne et d'au moins une deuxième réalisation indépendante du même modèle formée de N2 (t) correspondant à une deuxième valeur moyenne, cette combinaison étant telle que le nombre N (t) de la combinaison ait une valeur moyenne égale à la somme de la première et de la deuxième valeur moyenne. 2) Méthode selon la revendication 1, dans laquelle, à chaque itération, pour une même valeur moyenne de la combinaison, on fait varier de façon concomitante la première et la deuxième valeur moyenne de manière à faire varier graduellement le nombre d'objets issus de chacune des réalisations combinées. 3) Méthode selon la revendication 1 ou 2, dans laquelle on associe la taille d'objets à la procédure de génération du nombre d'objets de manière à faire apparaître ou disparaître un objet progressivement. 4) Méthode selon la revendication 1 ou 2, dans laquelle le nombre N (t) de la combinaison de moyenne .iL est relié aux nombres d'objets respectifs Ni (i=l...n) des n n réalisations combinées par la relation : N(t XA} = L N, {a, (t)A} avec La, =1, relation ,=1 ,=1 où at (t)# est la moyenne de la variable Ni au moment de la combinaison. 5) Méthode selon la revendication 4, dans laquelle pour une combinaison impliquant seulement deux réalisations, on choisit pour les coefficients ai, a2 des fonctions trigonométriques telles que par exemple les valeurs al = cos2(t) et a2 = sin2(t)