FR 2914434 A1 20081003 FR 0702405 A 20070330 La présente invention concerne un domaine technique de l'industrie pétrolière : la modélisation géologique de réservoir pétrolier. Plus particulièrement, l'invention concerne une méthode pour modifier graduellement un modèle géologique représentatif d'un réservoir souterrain de façon à respecter des 5 proportions moyennes fixées de ces faciès lithologiques. L'invention peut ainsi être directement appliquée au calage d'historique ( history matching ) d'un modèle géologique. L'objectif d'un tel calage est d'ajuster certains paramètres d'un modèle géologique de réservoir de façon à reproduire des données dynamiques mesurées sur le terrain. Il peut s'agir de données de production ou de données 10 issues de sismiques répétées (4D). L'optimisation et l'exploitation des gisements pétroliers reposent sur une description aussi précise que possible de la structure et des propriétés pétrophysiques du gisement étudié. Pour ce faire, les spécialistes utilisent un outil qui permet de rendre compte de ces deux 15 aspects de façon approchée : le modèle géologique. Un modèle géologique a donc pour vocation de rendre compte, au mieux, de la structure et des propriétés pétrophysiques d'un gisement. Le modèle géologique est constitué de différents éléments : un maillage qui constitue l'ossature du gisement et qui doit être représentatif de la structure, et des cartes tri ou bidimensionnelles de propriétés pétrophysiques associées à ce maillage et qui doivent être 20 représentatives du comportement statique ou dynamique du gisement. Dans le cadre de la modélisation géologique de réservoir pétrolier, la modélisation géostatistique est de plus en plus utilisée. Afin de prédire le comportement dynamique d'un gisement en production, des simulations numériques d'écoulement sont réalisées à partir du modèle géologique. Pour réduire les incertitudes et obtenir un modèle représentatif de la 25 réalité, le modèle géologique doit être conditionné par l'ensemble des données disponibles. Pour appliquer ce type de modélisation en intégrant différentes sources de données, on paramètre le modèle, c'est-à-dire que l'on définit des paramètres, que l'on appelle paramètres de calage du modèle . A titre d'exemple, les paramètres de calage peuvent être les propriétés pétrophysiques moyennes des faciès (perméabilité, porosité, etc.), des propriétés 30 géostatistiques telles que certains paramètres du variogramme (portée, direction d'anisotropie, etc.) ou les coefficients des lois liant certaines propriétés (rapport d'anisotropie verticale, loi porosité/perméabilité, etc.). La prise en compte des données statiques à l'aide de ce type de modélisation est bien connue. Elle permet de construire des modèles géologiques détaillés, appelés alors modèles géostatistiques, représentant au mieux l'ensemble des données statiques disponibles. Grâce à cette technique, les données de puits, les données de sismique et les observations géologiques (analogues, etc.) peuvent être intégrées dans un modèle géologique cohérent en tenant compte des incertitudes. Les données de puits permettent de contraindre localement les modélisations géostatistiques par des observations de faciès ou par des propriétés pétrophysiques mesurées (porosité, perméabilité, etc.). Les données de sismique sont utilisées pour caractériser la structure du réservoir. Dans certains cas, une analyse des attributs sismiques peut être réalisée pour définir les variations spatiales des proportions de faciès, de façon à introduire une contrainte non stationnaire dans le processus de modélisation géostatistique. La prise en compte des données dynamiques requiert une méthodologie adaptée pour modifier les modèles issus de modélisations géostatistiques appliquées aux données statiques. Cette étape consistant à ajuster un modèle existant pour qu'il soit cohérent avec des données dynamiques, est appelée calage d'historique . Ce processus de calage doit permettre d'ajuster certains paramètres du modèle géostatistique en conservant la cohérence avec les données statiques et les propriétés géostatistiques. État de la technique De nombreuses techniques de calage ont été développées au cours des dernières années pour paramétrer puis ajuster les modèles géostatistiques aux données dynamiques. Il s'agit de modifier itérativement le modèle géostatistique initial jusqu'à ce que le comportement dynamique simulé reproduise les données dynamiques collectées. Ce problème se résout en minimisant une fonction objectif qui quantifie la différence entre les données dynamiques et les réponses correspondantes simulées. On connaît des techniques de calage des données de production visant à préserver la structure spatiale du modèle géologique tout au long du processus de minimisation. En d'autres termes, il est possible de modifier les propriétés pétrophysiques attribuées aux mailles d'un modèle (un modèle est constitué d'un ensemble de mailles ayant chacune des coordonnées spatiales bien définies et relatives au réservoir modélisé) sans corrompre les propriétés statistiques globales de ces propriétés. On peut citer par exemple la méthode des points pilotes qui permet de modifier le modèle géologique localement autour des points pilotes : RamaRao, B.S., LaVenue, A.M., Marsilly, G. de, et Marietta, M.G., 1995, "Pilot point methodology for automated calibration of an ensemble of conditionally simulated transmissivity fields. 1. Theory and computational experiments", WRR, 31(3), 475-493. On connaît également la méthode des déformations graduelles, décrite dans les documents suivants : FR-2,780,798 (US 6.618.695), FR-2,795,841 (US 6.813.565), FR-2,821,946(US 7.151.845), FR 2.846.767 (US 2004/0138862) et FR 2.851.670 (US 2006/0149520). Cette méthode permet de perturber une réalisation géostatistique d'une fonction aléatoire (c'est à dire une image possible du modèle de réservoir) à partir d'un nombre réduit de paramètres tout en préservant la structure spatiale de cette même réalisation. L'avantage de cette méthode est qu'elle permet de modifier la réalisation globalement ou localement. Lorsque la déformation graduelle est globale, on considère un jeu de coefficients de déformation qui tous s'appliquent de la même manière à l'ensemble des mailles du modèle. Lorsque la déformation graduelle est locale, des jeux de coefficients de déformation sont affectés à des zones du modèle et toutes les mailles d'une même zone sont: identiquement affectées par le jeu de coefficient correspondant. Dans le cas d'une simulation en faciès, les déformations graduelles conduisent à modifier la distribution spatiale des hétérogénéités (décrites par la position spatiale des faciès), tout en respectant les propriétés géostatistiques moyennes (variogramme et proportions de faciès). Les déformations graduelles peuvent s'appliquer sur l'ensemble du modèle ou de façon régionale. Cependant, le changement de réalisation géostatistique ou des propriétés pétrophysiques ne suffit pas toujours à expliquer à lui seul les données dynamiques. Les proportions moyennes de faciès peuvent en particulier avoir un impact majeur sur le comportement dynamique du champ. Ces données de proportions peuvent être globales sur l'ensemble du réservoir dans le cas d'un modèle stationnaire (proportion moyenne de sables par exemple), ou déterminer les tendances ou variations spatiales dans le cas de modèles non stationnaires. Les proportions initiales sont en général déterminées à partir des données de puits et/ou de données de sismique. Elles peuvent être entachées d'incertitudes, aussi bien dans leurs valeurs globales que dans les tendances spatiales. Un objet de l'invention concerne une méthode permettant de modifier graduellement un modèle géologique en transformant localement ou globalement les proportions de faciès lithologiques, afin de contraindre ce modèle géologique par des données dynamiques dans un processus de calage d'historique par exemple. La méthode selon l'invention Ainsi, l'invention concerne une méthode pour modifier graduellement un modèle géologique représentatif d'un réservoir souterrain et constitué d'un ensemble de mailles renseignées en faciès lithologiques, de façon à respecter des proportions moyennes fixées de 10 ces faciès lithologiques, ces proportions moyennes étant calculées en moyennant des proportions de faciès sur au moins une zone géographique du réservoir. La méthode comporte les étapes suivantes : - on définit la zone géographique du réservoir en sélectionnant un ensemble de mailles du modèle géologique ; 15 on regroupe p faciès lithologiques dans un premier groupe de faciès ; au sein de ce premier groupe, on regroupe q faciès lithologiques dans un second groupe de faciès ; - on définit un paramètre de transformation par le rapport Rg de la proportion moyenne Xa,g des faciès du second groupe sur la proportion moyenne Xs g des 20 faciès du premier groupe ; - on transforme les proportions de faciès lithologiques du second groupe et de son complémentaire dans le premier groupe, en chaque maille de la zone géographique, de façon à imposer une valeur fixée au rapport Rg ; et - on simule un nouveau modèle géologique contraint par les proportions de faciès 25 lithologiques modifiées. Les proportions de faciès lithologiques du second groupe peuvent être transformées en appliquant une transformation du type : Vk e Fs, pk(m) = a(m,k).pk'O(m) avec: k un faciès appartenant à l'ensemble des p faciès lithologiques du premier groupe de faciès, pk (m) , la proportion du faciès k dans une maille m après transformation k,o p (m) , la proportion du faciès k dans une maille m avant transformation - a(m, k) , un coefficient fonction du paramètre de transformation Rg et des proportions de faciès avant transformation. Pour la simulation du nouveau modèle géologique on peut avantageusement utiliser une méthode de simulation par gaussiennes seuillées. Selon l'invention on peut contraindre les valeurs moyennes de proportions à l'aide d'une technique de krigeage, ce qui procure une continuité entre la zone subissant des transformations et l'extérieur de la zone où aucune transformation n'est appliquée. La méthode permet d'appliquer successivement plusieurs transformations. On peut par exemple modifier le contraste entre les faciès lithologiques en appliquant au modèle 15 géologique une transformation en puissance définie comme suit : dm E Zone, b'k e Facies, pk,final (in) = (pk (m))'' ~(p' (m))fi iEFacies Facies : ensemble des faciès lithologiques du modèle géologique. Zone : une zone du domaine définie par l'utilisateur et correspondant à un ensemble de mailles. pk (m) : proportion locale du faciès k dans la maille m pk,fnar (m) proportion du faciès k dans la maille m après application de ladite transformation en puissance. paramètre de contraste, nombre réel compris entre 0 et 1. 25 Enfin, selon l'invention, la méthode peut être appliquer au calage d'un modèle géologique d'un réservoir souterrain de façon à reproduire des données dynamiques mesurées, où : 20 ce modèle étant constitué d'un ensemble de mailles renseignées en faciès lithologiques. Selon une telle application, on mesure une différence entre les valeurs de données dynamiques mesurées et des valeurs de données dynamiques simulées à l'aide d"un simulateur d'écoulement appliqué au modèle géologique. Puis on modifie graduellement le modèle géologique de façon à minimiser cette erreur. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation succincte des figures La figure 1 est un schéma illustrant les étapes de la méthode de calage d'historique selon l'invention. La figure 2 illustre le principe de la méthode des gaussiennes seuillées. 15 La figure 3 illustre un modèle géologique initial (à gauche) et une définition de régions chenaux et levées (à droite). La figure 4 représente la distribution spatiale des faciès sur le modèle réservoir avant et après calage des proportions de faciès pour les couches 1 et 2. La figure 5 représente la distribution spatiale des faciès sur le modèle réservoir avant et 20 après calage des proportions de faciès pour les couches 3 et 4. Description détaillée de la méthode Un objet de l'invention concerne une méthode permettant de réaliser un calage d'un modèle géologique par rapport à des données dynamiques en modifiant localement ou 25 globalement les proportions de faciès lithologiques. On appelle faciès lithologique une propriété d'une roche. Par exemple, le faciès lithologique peut se référer à la nature géologique de la roche (argile ; grès ; calcaire...), à son type de porosité (roche non consolidée et très poreuse; roche peu poreuse...), à la nature du fluide piégé dans sa porosité (saumure ; huile ; gaz...).10 On appelle modèle géologique une carte (représentation en deux dimensions) ou un cube (représentation en trois dimensions) constitué de mailles associées à des propriétés liées à un réservoir souterrain étudié. Un tel modèle rend compte également des caractéristiques structurales de ce réservoir. La méthode est décrite à travers un exemple particulier de réalisation. concernant le calage d'historique des données de production, mesurées au niveau de puits producteurs d'hydrocarbures. De nos jours, les boucles de calage d'historique se résument comme suit. Tout d'abord, 10 un modèle géologique représentatif d'un réservoir pétrolier est élaboré à une échelle fine. Ce modèle respecte les données statiques, c'est-à-dire les données qui ne dépendent pas des mouvements de fluide au sein du réservoir modélisé. Puis, ce modèle géologique est transféré vers une échelle grossière, appelée échelle réservoir . Ces notions d'échelles sont parfaitement connues des spécialistes. On obtient alors un modèle, dit modèle de réservoir sur 15 lequel on effectue des simulations numériques d'écoulement à l'aide d'un simulateur d'écoulement, type de logiciel bien connu des spécialistes. Le passage à l'échelle du modèle de réservoir permet de réduire le nombre de mailles par rapport au modèle géologique fin. De fait, le temps CPU requis pour une simulation d'écoulement devient acceptable en pratique. La simulation d'écoulement fournit des données de production simulées, qui sont alors 20 comparées aux données de production effectivement mesurées. Si l'écart entre les données dynamiques simulées et les données dynamiques mesurées est suffisamment faible, le processus de calage s'arrête. Dans le cas contraire, on revient au modèle géologique fin que l'on modifie de façon à améliorer le calage des données de production. Le processus itératif est alors répété jusqu'à ce que le calage soit jugé satisfaisant. 25 Selon la méthode, on intègre à cette boucle de calage d'historique classique une paramétrisation adaptée permettant d'ajuster les proportions de faciès lithologiques du modèle géologique d'un réservoir souterrain de façon à reproduire des données dynamiques mesurées sur le terrain. Ces données dynamiques peuvent être des données de production ou des données issues de sismiques répétées (4D). 30 Selon l'invention, on dispose d'un modèle géologique initial décrivant le réservoir étudié en faciès. De nombreuses techniques permettant de construire un modèle géologique contraint sont connues des spécialistes. On peut par exemple citer une technique de5 modélisation non stationnaire appelée simulation par gaussiennes seuillées. Cette technique est décrite dans le document suivant : Le Loc'h G. and Galli A., 1997 : Truncated plurigaussian method : theoretical and practical points of view. In : Geostatistics Wollongong'96, E.Y. Baafi and N.A Schofield eds, Kluwer, p.211-222 Selon l'invention, chacune des mailles du modèle géologique est renseignée au moins par les valeurs de proportions de faciès lithologiques, par exemple : 15% de sable, 30% de grès, 5% d'argile, etc. Les proportions de faciès peuvent être identiques sur toutes les mailles d'une couche du modèle (on parle alors d'un modèle stationnaire) ou être variables d'une maille à une autre (cas d'un modèle non stationnaire). Ces proportions initiales de faciès géologiques peuvent être déterminées à partir des données de puits et/ou d'attributs ou de données sismiques. Une telle étape est bien connue des spécialistes et l'on peut citer par exemple les documents suivants : Doligez B., Fournier F., Jolivet G., Gancarski S., Beucher H., 2002: Seismic facies map integration in geostatistical geological model: a field case. EAGE, conference & technical exhibition of European Association of Geoscientists & Engineers, 64th, Florence, 27-30 May 2002, Extended abstracts, Vol. 2, P215-219. Barens L., Biver P., 2004: Reservoir Facies Prediction from Geostatistical Inverted Seismic Data, Abu Dhabi International Conference and Exhibition, 10-13 October, SPE 88690-MS. - la demande de brevet française FR 06/10.792. A partir des proportions de faciès et des données de puits, une technique de simulation géostatistique, telle que la méthode des gaussiennes seuillées, est utilisée pour générer une ou plusieurs réalisations de distributions de faciès lithologiques sur le maillage géologique. Pour chaque réalisation, la probabilité de simulation d'un faciès lithologique sur une maille est conditionnée par les données de proportions de faciès. De façon avantageuse ce modèle géologique est également renseigné en propriétés pétrophysiques (porosité, perméabilité, ...), qui sont assignées aux faciès lithologiques. Le réservoir pétrolier est équipé de puits pour lesquels on connaît des données de production, à partir desquelles on souhaite caler le modèle géologique. Partant d'un modèle géologique initial (MGt), la méthode comporte principalement les étapes suivantes (figure 1) : a) on définit un paramètre de transformation (Rg) des proportions moyennes de faciès 5 lithologiques ; b) on détermine la valeur que doit avoir ce paramètre pour que le modèle cale au mieux aux données de production. Dans le cadre du calage d'historique ceci est réalisé à l'aide d'un simulateur d'écoulement (SIM), type de logiciel bien connu des spécialistes. A partir d'un modèle géologique représentatif du gisement réel, un tel simulateur permet de calculer 10 une prédiction des données de production (DDsim). Si l'on modifie le modèle, et notamment les proportions de faciès, la réponse, c'est-à-dire la simulation des données de production, est différente. On peut ainsi, en comparant les données de production simulées et celles mesurées (DDmes), déterminer une valeur du paramètre de transformation, tel qu'une fois modifié, le modèle géologique offre une réponse plus 15 proche de la réalité ; c) on modifie (Modif) les proportions de faciès lithologiques du modèle géologique, de façon à ce que le paramètre ait la valeur de calage déterminée ; d) on simule (GS) un nouveau modèle en faciès en fonction de ces nouvelles proportions. Les étapes b) à d) peuvent être répétées à chaque boucle de calage telles que décrite 20 précédemment. Modification des proportions de faciès lithologiques d'un modèle géologique La méthode décrite ci-après correspond aux étapes a) et c) du calage. Elle constitue un autre objet de l'invention qui concerne une méthode permettant de modifier les proportions de 25 faciès lithologiques d'un modèle géologique de façon à respecter des proportions moyennes fixées de ces faciès lithologiques. On appelle proportions moyennes de faciès lithologiques , l'estimation globale des proportions sur un volume de discrétisation formé par un ensemble donné de mailles. On appelle proportions locales de faciès lithologiques les données de proportions dans un 30 volume réduit à une maille du modèle. La méthode propose de sélectionner un groupe de faciès, appelé sélection, et de contraindre dans cette sélection la proportion moyenne d'un sous-groupe de faciès, appelé association. Autrement dit, on sélectionne un groupe de faciès (la sélection), puis un sous-groupe (l'association) et on définit un paramètre par le rapport entre les proportions moyennes des faciès de l'association et celles de la sélection. Ce paramètre permet de modifier les proportions des faciès de la sélection, et ceci sur l'ensemble du modèle ou sur une zone donnée formée par un ensemble de mailles. Lorsque le modèle est stationnaire, la modification des proportions est appliquée de façon uniforme sur l'ensemble des mailles de la zone. Lorsque le modèle est non-stationnaire, les proportions de chaque maille sont modifiées en tenant compte des proportions locales initiales, de façon à respecter le paramètre de proportion moyenne. La méthode propose également de définir un autre paramètre d'ajustement des proportions, appelé paramètre de contraste. Le principe de ce paramètre de contraste est d'agir sur la proportion du faciès lithologique le plus probable par rapport à la proportion des autres faciès. Ce paramètre de contraste peut également être appliqué à l'ensemble du modèle ou sur une zone correspondant à un ensemble de mailles. Une loi de modification en puissance est proposée ici, mais le principe de paramètre de contraste peut être étendu à d'autres lois. Pour exposer de façon plus détaillée la méthode de modification des proportions de 20 faciès, on introduit tout d'abord les notations et définitions suivantes : Facies : ensemble des faciès lithologiques du modèle géologique. Zone une zone du domaine définie par l'utilisateur (que nous supposerons unique pour simplifier) et correspondant à un ensemble de mailles. pk (m) : proportion locale du faciès k dans la maille m. 25 pgk proportion moyenne du faciès k. Dans le cas stationnaire, les proportions moyennes et locales coïncident, et dans le cas non stationnaire : pgk = 1 Epk(m) tard (Zone) m{zone card(Zone) est le cardinal de Zone, c'est-à-dire le nombre de mailles contenues dans la zone Zone. FS ensemble des faciès lithologiques contenus dans la sélection. Fa ensemble des faciès lithologiques de l'association. Xs g : proportion moyenne de la sélection, Xs g = l pg k keF X, ,,g : proportion moyenne de l'association, Xa g = E pg k kEF X rapport de proportion moyen a'g S,g X, (m) : proportion locale de la sélection dans la maille m, Xs (m) = I pk (r) kEF Xa (m) : proportion locale de l'association dans la maille m, Xa (m) = pk (m) kEF R(m) : rapport de proportion local, X a (m) XS(m) Q paramètre de contraste correspondant à une loi en puissance (exposant des 10 proportions de faciès) Pour toutes ces notations, un exposant indique une valeur initiale, c"est-à-dire avant modification par la méthode. La méthode proposée permet de contraindre la proportion moyenne des faciès associés dans la sélection en agissant sur le paramètre Rg . Toutes les proportions de faciès de la 15 sélection sont affectées par cette contrainte, tout en gardant les proportions des faciès restants (i.e. hors sélection) inchangées. Cette méthode comporte principalement les étapes suivantes : on définit au moins une zone géographique du réservoir dans laquelle on souhaite modifier des proportions de faciès pk (m) , en sélectionnant un ensemble de mailles du modèle géologique (la zone peut être étendue à l'ensemble du réservoir) ; 20 on regroupe p faciès lithologiques dans un premier groupe de faciès, appelé sélection ; Rg - au sein de ce premier groupe, on regroupe q faciès lithologiques dans un second groupe de faciès, appelé association ; - on définit un paramètre de transformation par le rapport Rg de la proportion moyenne Xa g des faciès du second groupe sur la proportion moyenne XS g des faciès du premier groupe : R _ Xa,g g Xs,g On calcule Rg pour le modèle initial. Puis on fixe une valeur à respecter : Rg = Rg R. on modifie les proportions de faciès lithologiques de l'association (Fa), et de son complémentaire dans la sélection (FS \ Fa), en chaque maille de la zone géographique et au prorata des proportions initiales pk'0 (m) , de façon à satisfaire le rapport Rg . La valeur Rg peut être définie en fonction de données additionnelles (données de puits, interprétation de tests de puits, contrainte géologique, contrainte sismique, ,etc.) ou dans le processus itératif de calage d'historique. Cas stationnaire Dans un cas stationnaire, les proportions de faciès pk (m) sont identiques en tout point de la grille. Par conséquent : pgk = pk(m) . R En notant A= g R g1ùR et B = g , on modifie les proportions 1ùR g de faciès de l'association, et de son complémentaire dans la sélection, de la façon suivante : Vk E Fa , k (Eq. 1) p (m) = pk = A x pg•o VkEFS\Fa, pk (m) = pk = B x pg o (Eq. 2) Ainsi : Xa (m) = Epk (m) _ >2 A x pg keF keF Trois cas particuliers peuvent se produire : • Si Xs,g = 0, la transformation ne peut être appliquée et les proportions de faciès doivent rester inchangées. • Si R g =1, ceci signifie qu'initialement toutes les proportions des faciès de l'ensemble Xs \ Xa sont nulles. Dans ce cas, la formule (Eq. 2) ne peut être appliquée et on supposera que les faciès du groupe XS \ Xa apparaissent en proportions identiques. Ainsi, dans ce cas : Vk E F a , pk (m) = pk = A x pk O (Eq. 1) 13 b'k E Fs \ Fa , k k (1ùRg)X g p (m) pg ù tard (FS \F1) (Eq. 2') • Si R0 g = 0 , on procède de la même façon : pk (m) = pg = RgX s,g tard (Fa ) VkEFs\Fa, pk(m)g=Bxpk'0 (Eq. 2) Cas non stationnaire Dans un cas non stationnaire, la modification est analogue au cas précédent et préserve, dans une certaine mesure, la répartition locale des faciès observée initialement. En effet, la modification globale des proportions de faciès est distribuée localement en fonction des proportions initiales de chaque maille. Ainsi, en utilisant les mêmes notations : b'k e Fa , k,0 m (Eq. 3) Aroc (m) pk (ni) ù pk,o (m) + x Xo VkEFs\Fa, k,0(m) (Eq. 4) avec : P Broc pk (m) = pk,O (m) + x o XS (m) ù Xa (m) dkEFa, (Eq. 1') 10 (A -1)X, (m), si Rg Ar o ù (B -1)(X (m) ù X (m)), sinon. Broa = ùAroc Comme pour le cas stationnaire, la transformation ne peut être appliquée lorsque 5 XS g = 0 . Remarquons que les grandeurs Aron et Bk), sont toujours bien définies, y compris pour R g = 0 et R g = 1 . Localement, des cas particuliers peuvent apparaître. En effet, pour Xs (m) = 0 , les proportions de faciès doivent rester inchangées pour la maille in, pour X 0a (m) = 0 , alors Vk E Fa, pk (m) = A` ` card (Fa ) pour x s m x a m, alors ~/k E F\ F k m B10 ( ) = ( ) a' p ( ) = card (Fs \ Fa) En somme, quelque soit le cas, on transforme les proportions de faciès lithologiques de l'association, en appliquant une transformation du type : Vk E FS , pk (m) = a(m, k)•pk'0 (m) avec: 15 pk (m) , la proportion du faciès k dans une maille m après transformation, - pk'O (m) , la proportion du faciès k dans une maille m avant transformation, - a(m, k) , un coefficient fonction du paramètre de transformation Rg et des proportions de faciès initiales. 20 Ainsi, à l'aide d'un seul paramètre global Rg , on parvient à modifier les proportions locales des faciès de la sélection sur l'ensemble de la zone à modifier. Quelque soit le cas, la méthode permet de tenir compte des données observées aux puits (données conditionnantes) : les mailles traversées par des puits sont exclues de la zone de transformation et les faciès observés leur sont directement attribués. Les transformations définies ci-dessus peuvent être appliquées sur des régions préalablement définies par unensemble de mailles (transformation régionale) ou sur l'ensemble du modèle (transformation globale, la région est le modèle complet). Dans le cas de transformations régionales, on définit une zone particulière du réservoir dans laquelle on souhaite modifier les proportions de faciès. Ces proportions moyennes de faciès sont calculées de façon indépendante sur chaque région. La méthode de transformation n'affecte que les proportions de faciès de la zone considérée. Selon un autre mode de réalisation, on utilise un second paramètre de transformation, appelé paramètre de contraste, pour calculer les proportions de faciès finales sur l'ensemble des mailles du modèle géologique. Ce paramètre de contraste permet d'accentuer ou de réduire la probabilité d'apparition des faciès lithologiques dont la probabilité est dominante. Dans un processus de simulation géostatistique, un paramètre de contraste élevé se traduira par des réalisations géostatistiques se rapprochant du faciès le plus probable. A l'inverse, un paramètre de contraste faible tendra vers une apparition équiprobable de faciès sur un ensemble de réalisations indépendantes, se traduisant par des proportions également réparties entre faciès. Pour faire varier le contraste entre les faciès, une transformation en puissance peut être appliquée, avec la formulation suivante : Vm E Zone,Vk e Facies, pk,final (m) _ (pk (m))Q (p' (m iEFacies Le paramètre f est le paramètre de contraste, il s'agit d'un nombre réel compris entre 0 et 1. Il permet d'accentuer la proportion des faciès les plus représentés par des valeurs d'exposant supérieures à un. Lorsque le paramètre est inférieur à un, l'influence du faciès le plus représenté est diminuée et on tend vers des proportions également réparties de faciès pour des valeurs faibles. Lorsque le paramètre vaut un, les proportions initiales sont inchangées. Les nouvelles valeurs obtenues sont normalisées pour que la somme des 30 proportions transfotuiées soit unitaire. Ce paramètre de contraste peut être appliqué, de façon globale ou par régions, une fois toutes les autres transformations effectuées ou bien avant les autres transformations. Dans le cas où les proportions de faciès initiales sont déterminées à partir d'attributs ou de données sismiques, l'augmentation du paramètre de contraste permet d'accentuer l'impact de la contrainte sismique par rapport aux réalisations aléatoires. A l'extrême, un paramètre de contraste très élevé se traduira par des réalisations déterministes, correspondant à la distribution du faciès le plus probable. Modification des faciès lithologiques d'un modèle géologique La méthode décrite ci-après correspond à l'étape d) du calage. Elle constitue un autre objet de l'invention qui concerne une méthode permettant de modifier graduellement un modèle géostatistique en modifiant localement ou globalement les proportions de faciès lithologiques. La méthode comporte principalement les étapes suivantes : - on définit au moins une zone géographique du réservoir dans laquelle on souhaite modifier des proportions de faciès p' (m), en sélectionnant un ensemble de mailles du modèle géologique ; - on modifie des proportions de faciès lithologiques en appliquant la méthode précédente de modification des proportions de faciès lithologiques, et notamment en modifiant les proportions de faciès lithologiques d'une association (Fa ), et de son complémentaire dans sa sélection (FS \ Fa ), en chaque maille de la zone géographique, de façon à imposer une valeur fixée au paramètre de transformation Rg. Un paramètre de contraste peut également être appliqué ; on simule un nouveau modèle géologique contraint par lesdites proportions de faciès lithologiques modifiées. De préférence on utilise une méthode de simulation qui permet de conférer à cette transformation un caractère graduel et continu. On peut par exemple utiliser une simulation basée sur les gaussiennes seuillées déjà citée. En effet cette méthode consiste à générer une réalisation gaussienne d'une fonction 30 aléatoire. Cette réalisation gaussienne d'une variable continue est transformée en catégories,25 correspondant aux faciès lithologiques, en la tronquant par des valeurs de seuils. Les seuils définissent des intervalles d'appartenance à un faciès lithologique. La figure 2 illustre le principe de cette méthode des gaussiennes seuillées. Elle montre une réalisation gaussienne lD d'une fonction aléatoire (GRF) et deux seuils (SI et S2) permettant de définir trois faciès (1, 2 et 3). Ainsi à chaque valeur simulée sur une grille, on peut faire correspondre un faciès. Les seuils sont définis en fonction des proportions de faciès lithologiques. Un intervalle entre 2 seuils est faible si le faciès correspondant à une faible proportion. En d'autres termes, ces proportions permettent de déduire la position des seuils. Selon la méthode on modifie les proportions de faciès. On modifie donc la position des seuils. En réalisant une simulation par la méthode des gaussiennes seuillées on obtient un nouveau modèle en faciès. Comme la réalisation gaussienne est une variable continue, si l'on déplace graduellement un seuil, la modification de la distribution de faciès est nécessairement graduel. Selon un mode particulier de réalisation, plusieurs transformations peuvent être 15 appliquées successivement de façon à introduire différents paramètres de calage. Ces transformations peuvent être différentes. Cette méthode est peu coûteuse en temps de calcul, mais elle ne tient pas compte des données de contraintes de continuités pour modifier la répartition des proportions de faciès : les transformations des proportions sont effectuées indépendamment les unes des autres sur 20 chaque région. Notons cependant que la continuité des réalisations en faciès est maintenue du fait qu'un variogramme est appliqué dans le processus de simulation géostatistique. Lorsqu'une continuité doit être maintenue sur la répartition des proportions de faciès, une approche par krigeage peut être avantageusement combinée, de façon à lisser les transformations et assurer une continuité entre régions. Dans ce cas, la technique de krigeage 25 permet de contraindre les valeurs moyennes de proportions par région, tout en gardant une évolution progressive de la transformation déterminée par un variogramme. :Pour ce faire, on procède comme suit. En fonction de la transformation choisie (par l'utilisateur), on calcule les valeurs R (m) en chaque maille m du modèle y compris en dehors de la zone de 30 transformation. On considère cette fonction discrète R comme une première distribution que l'on conditionne ensuite aux puits en fonction des faciès observés mais aussi sur la zone de transformation en imposant une valeur moyenne Rg rentrée par l'utilisateur ou déduite d'une optimisation. On obtient alors une seconde distribution Rond . Aux puits les valeurs de Rcond sont égales à 1 ou 0 si le faciès observé appartient à l'ensemble Fa ou FS \ Fa . Il peut aussi arriver que ce faciès ne soit pas présent dans la sélection, dans ce cas, Rond n'est pas définie, on ne cherchera donc pas à conditionner Rcond en ce point. Le conditionnement peut être réalisé à l'aide d'un algorithme de krigeage simple où des données de variogramme doivent être spécifiées en plus des données relatives aux transformations. Le conditionnement par krigeage est une technique parfaitement connue des spécialistes. Le krigeage simple et le krigeage par bloc sont, en particulier détaillés dans : P. Goovaerts : Geostatistics for Natural Resources Evaluation, Oxford University Press, 10 1997, pp. 125-185. A partir des nouvelles valeurs Rcond (m) , on recalcule les proportions des faciès de la sélection en chaque maille du modèle de la façon suivante : pk (m) = A(m) x pk•0 (m) , avec A(m) = Rc nd (m) R (m) ' • b'kEFS\Fa, pk( )ù pk,o( )~ -1ùRcood(m). mù B(m) x m avec B(m) - 1-R (m) 15 Des cas particuliers apparaissent aux mailles m où X (m) = 0, R (m) = 1 ou encore R (m) = 0 . Ces cas particuliers peuvent être résolus d'une façon analogue au cas stationnaire. Exemple d'application L'exemple d'application présenté dans cette section est basé sur un modèle synthétique 20 de réservoir, représentant une simulation de test de puits dans un chenal. Un modèle initial est construit à l'aide d'une simulation géostatistique par la méthode des gaussiennes seuillées. Les proportions de faciès initiales, qui peuvent être issues de données de sismique 3D, donnent les tendances principales de répartition spatiale des faciès. Le modèle réservoir est supposé décrit avec quatre faciès, qui sont par ordre de qualités réservoir croissantes : "argiles" (Fl), "argiles 25 laminées" (F2), "sables laminés" (F3) et "sables" (F4). Les faciès "argiles" et "argiles laminés" (Fl+F2) sont présents dans les régions levées, et les faciès "sables laminés" et "sables" (F3+F4) dans les chenaux. Les proportions moyennes d'argiles laminées dans les • dkEFa, levées et de sables dans les chenaux sont entachées d'une marge d'incertitude que l'on cherche à réduire en contraignant les proportions moyennes par les données de tests de puits. Le modèle géologique initial, donné par la distribution spatiale des faciès, est présenté sur la Figure 3. Le modèle géologique initial est présenté à gauche. Selon une vue de dessus (axes X et Y) pour un niveau C3 en haut de la figure 3, et selon une section verticale en fonction de la profondeur Z pour Y=26, en bas de la figure 3. On distingue la géométrie du chenal, dans lequel est implanté un puits de production (PW). Le puits est perforé sur l'ensemble de la hauteur du réservoir. La définition des régions chenaux et levées est présentée à droite de la figure 3. Les régions sont utilisées pour calculer et: paramétrer les proportions moyennes de faciès en vue du calage des données de test de puits. Les quatre faciès F1 à F4 sont représentés. Les propriétés pétrophysiques (porosité et perméabilité) sont constantes par faciès. Les valeurs suivantes (Tableau 1) ont été adoptées pour les simulations d'écoulement : Faciès Porosité Perméabilité Rapport d'anisotropie (adim) horizontale verticale Kv/Kh (mD) (adirn) F1 0.02 10 0.1 F2 0.08 50 0.1 F3 0.15 300 0.1 F4 0.25 800 0.1 Tableau 1 - Valeurs des propriétés pétrophysiques par faciès. Les régions chenaux et levées ont été délimitées à partir des faciès (Figure 3) : les régions levees correspondent à la présence de faciès FI et F2, les régions chenaux à la présence de faciès F3 et F4. Ces régions sont utilisées pour paramétrer les proportions moyennes de faciès dans les chenaux et levées, en vue du calage de ces paramètres par les données de test de puits: • un premier paramètre de transformation correspond à la proportion moyenne d'argiles laminées dans les levées: %F2 • un second paramètre de transformation correspond à la proportion moyenne de sables dans les chenaux: %F3 %(F3 + F4) %(F1 + F2) Un test de puits synthétique a été simulé à partir de proportions de faciès de référence dans un cadre non-stationnaire. Le modèle de référence a été généré avec les valeurs suivantes de paramètres de transformation : 0.46 pour le premier paramètre et 0.51 pour le second paramètre. Le test de puits (draw-down) est simulé sur une période de 100 heures, avec un débit constant de 650 m3/jour. Les données synthétiques de tests de puits comprennent la différence de pression et sa dérivée. Pour le calage du test de puits, les valeurs initiales ont été fixées à 0.2 pour les deux paramètres de transformation, qui peuvent varier entre 0 et 1 au cours du calage. Le Tableau 2 ci-dessous synthétise les valeurs des paramètres de proportions pour le modèle de référence et pour le modèle initial. Les proportions de référence (colonne de droite) sont les valeurs utilisées pour générer les données synthétiques de test de puits. Les valeurs initiales constituent le point de départ de l'optimisation réalisée pour le calage du test de puits. Région Faciès de la sélection Faciès de Proportions initiales Valeurs des proportions l'association des faciès de de référence l'association au sein de la sélection Levées F1 + F2 F2 0.2 0.46 Chenaux F3 + F4 F4 0.2 0.51 Tableau 2 - Définition et valeurs des paramètres de proportions moyennes de faciès. Un algorithme d'optimisation, basé sur une méthode de gradients, a été mis en oeuvre pour caler le test de puits en ajustant les deux paramètres de proportions., en partant des valeurs initiales. L'algorithme réalise une minimisation itérative d'une fonction objectif du type moindres carrés, qui quantifie l'écart entre la simulation numérique du test de puits et les données. L'optimisation converge en 7 simulations et permet de retrouver des valeurs moyennes de proportions de faciès très proches des valeurs de référence: 0.453 pour le paramètre 1 (pour une valeur de référence de 0.46) et 0.509 pour le paramètre 2 (pour une valeur de référence de 0.51). Les valeurs optimales de paramètres permettent une très bonne restitution des données de test de puits. La Figure 4 et la Figure 5 montrent l'impact du calage des proportions de faciès sur le modèle géologique. La figure 4 montre la distribution spatiale des faciès sur le modèle réservoir avant et après calage des proportions de faciès. Les images de gauche représentent les couches 1 et 2 (Cl, C2) obtenues avec les proportions initiales, les images de droite les mêmes couches avec les proportions finales. La figure 5 montre la distribution spatiale des faciès sur le modèle réservoir avant et après calage des proportions de faciès. Les images de gauche représentent les couches 3 et 4 (C3, C4) obtenues avec les proportions initiales, les images de droite les mêmes couches avec les proportions finales. Le calage du test de puits a été réalisé en augmentant la proportion moyenne du faciès F2 (argiles laminées) dans les levées et la proportion moyenne de faciès F4 (sables) dans les chenaux. La répartition spatiale des faciès montre une étendue plus importante de ces faciès sur le modèle optimal par rapport au modèle initial. On note que l'optimisation est réalisée sur une réalisation géostatistique fixée (germe figé au cours de l'optimisation). Il est également possible de caler simultanément les proportions de faciès et la réalisation en utilisant la méthode des déformations graduelles. Dans ce cas, le calage des proportions de faciès permet de modifier les tendances principales, tandis que la méthode des déformations graduelle permet de modifier la répartition locale des hétérogénéités à proportions fixées. Cette méthodologie est donc complémentaire par rapport à la méthode des déformations graduelles : le calage des proportions de faciès introduit un contrôle sur les tendances spatiales moyennes tandis que les déformations graduelles agissent sur la distributions locale des hétérogénéités à proportions moyennes fixées.25 REVENDICATIONS 1) Méthode pour modifier graduellement un modèle géologique représentatif d'un réservoir souterrain et constitué d'un ensemble de mailles renseignées en faciès lithologiques, de façon à respecter des proportions moyennes fixées de ces faciès lithologiques, lesdites proportions moyennes étant calculées en moyennant des proportions de faciès sur au moins une zone géographique du réservoir, caractérisée en ce qu'elle comporte les étapes suivantes : - on définit ladite zone géographique du réservoir en sélectionnant un ensemble de mailles dudit modèle géologique ; on regroupe p faciès lithologiques dans un premier groupe de faciès ; au sein de ce premier groupe, on regroupe q faciès lithologiques dans un second groupe de faciès ; on définit un paramètre de transformation par le rapport Rg de la proportion moyenne Xa,g des faciès du second groupe sur la proportion moyenne Xs g des 15 faciès du premier groupe ; on transforme lesdites proportions de faciès lithologiques du second groupe et de son complémentaire dans le premier groupe, en chaque maille de la zone géographique, de façon à imposer une valeur fixée au rapport Rg ; et on simule un nouveau modèle géologique contraint par lesdites proportions de faciès 20 lithologiques modifiées. 2) Méthode selon la revendication 1, dans laquelle on transforme lesdites proportions de faciès lithologiques du second groupe, en appliquant une transformation du type : `dk E FS , pk (m) = a(m, k)•pk'O (m) avec: 25 - k un faciès appartenant à l'ensemble des p faciès lithologiques dudit premier groupe de faciès, pk(m) , la proportion du faciès k dans une maille m après transformation, k,o p (m) , la proportion du faciès k dans une maille m avant transformation, a(m, k) , un coefficient fonction du paramètre de transformation Rg et des proportions de faciès avant transformation. 3) Méthode selon l'une des revendications 1 et 2, dans laquelle on simule un nouveau modèle géologique par la méthode des gaussiennes seuillées. 4) Méthode selon l'une des revendications 1 à 3, dans laquelle on contraint les valeurs moyennes de proportions à l'aide d'une technique de krigeage. 5) Méthode selon l'une des revendications précédentes, dans laquelle on applique successivement plusieurs transformations. 6) Méthode selon la revendication 5, dans laquelle on modifie un contraste entre lesdits faciès lithologiques en appliquant au modèle géologique une transformation en puissance définie comme suit : Vm E Zone,Vk e Facies, pk,final (m) _ (pk (m))R E(p` (m))a iEFacies ou : Facies : ensemble des faciès lithologiques du modèle géologique. Zone : une zone du domaine définie par l'utilisateur et correspondant à un ensemble de mailles. pk (m) : proportion locale du faciès k dans la maille m. proportion du faciès k dans la maille m après application de ladite transformation en puissance. paramètre de contraste, nombre réel compris entre 0 et 1. 7) Application au calage d'un modèle géologique d'un réservoir souterrain de façon à reproduire des données dynamiques mesurées, dans laquelle on mesure une différence entre lesdites valeurs de données dynamiques mesurées et des valeurs de données dynamiques simulées à l'aide d'un simulateur d'écoulement appliqué audit modèle géologique, et l'on 25 modifie graduellement ledit modèle géologique selon l'une des revendications précédentes, de façon à minimiser ladite erreur. pk,final (m) :20