FR 2949586 A1 20110304 FR 0904083 A 20090826 10 La présente invention concerne le domaine technique de l'industrie pétrolière, et plus particulièrement la caractérisation des réservoirs pétroliers par construction d'une représentation du réservoir, appelée modèle de réservoir. L'optimisation et l'exploitation des gisements pétroliers reposent sur une 15 description aussi précise que possible de la structure, des propriétés pétrophysiques, des propriétés fluides, etc., du gisement étudié. Pour ce faire, les spécialistes utilisent un outil qui permet de rendre compte de ces aspects de façon approchée : le modèle de réservoir. Un tel modèle constitue une maquette du sous-sol, représentative à la fois de sa structure et de son comportement. Généralement, ce type de maquette est 20 représenté sur un ordinateur, et l'on parle alors de modèle numérique. Ces modèles bien connus et largement utilisés dans l'industrie pétrolière, permettent de déterminer de nombreux paramètres techniques relatifs à la recherche, l'étude ou l'exploitation d'un réservoir, d'hydrocarbures par exemple. En effet, un modèle de réservoir est représentatif de la structure du réservoir ainsi que de son 25 comportement. On peut ainsi par exemple déterminer quelles sont les zones qui ont le plus de chances de contenir des hydrocarbures, les zones dans lesquelles il peut être intéressant/nécessaire de forer un puits d'injection pour améliorer la récupération des hydrocarbures, le type d'outils à utiliser, les propriétés des fluides utilisés et récupérés.... Ces interprétations de modèles de réservoir en termes de � paramètres 30 techniques d'exploitation � sont bien connues des spécialistes, même si de nouvelles méthodes sont régulièrement développées. II est donc crucial, dans le domaine pétrolier, de construire un modèle de réservoir aussi précis que possible. Pour ce faire, l'intégration de toutes les données disponibles est indispensable. Un modèle de réservoir a donc pour vocation de rendre compte, au mieux, de toutes les informations d'un gisement. Un modèle de réservoir est représentatif lorsqu'une simulation de réservoir fournit des estimations de données d'historique très proches des données observées. On appelle données d'historique, les données de production issues de mesures aux puits en réponse à la production du réservoir (production d'huile, production d'eau d'un ou plusieurs puits, ratio gaz/huile (GOR), proportion d'eau de production ("water cut")), et/ou les données de sismiques répétitives (impédances sismiques 4D dans une ou plusieurs régions, etc.). Une simulation de réservoir est une technique permettant de simuler les écoulements de fluides au sein d'un réservoir au moyen d'un logiciel appelé simulateur d'écoulement. Le calage d'historique ("history matching") consiste à modifier les paramètres d'un modèle de réservoir, tels que les perméabilités, les porosités ou les skins de puits (représentant les endommagements autour du puits), les connections de failles..., pour minimiser les écarts entre les données d'historique simulées et mesurées. Les paramètres peuvent être liés à des régions géographiques comme les perméabilités ou porosités autour d'un puits ou plusieurs puits. Les techniques de calage d'historique assisté sont largement utilisées pour caractériser un réservoir, en intégrant les données de puits et les données sismiques. On connaît par exemple les techniques décrites dans les documents suivants : Roggero, F. et Hu, L.Y.: "Graduai Deformation of Continuous Geostatistical Models for History Matching", paper SPE 49004, Proc. SPE Annual Technical Conference 25 and Exhibition, New Orleans, USA, 1998; Gosselin, O., Cominelli, A. van den Berg, S. et Chowdhury, S.D.: "A Gradient-Based Approach for History Matching of Both Production and 4D Seismic Data", Proceeding 7th European Conference on the Math. of Oil Recovery, Baveno, Italy, 2000; 30 - Cheng, H., Wen, X., Milliken, W.J. et Datta-Gupta, A.: "Field Experiences with Assisted and Automated History Matching", papier SPE 89857, SPE ATC&E, Houston, TX, USA, 2004; Roggero, F., Ding, D.Y., Berthet, P., Lerat, O., Cap, J. et Schreiber, P.E.: "Matching of Production History and 4D Seismic Data - Application to the Girassol Field, Offshore Angola", papier SPE 109929, SPE ATC&E, Anaheim, California, USA, 2007. Au cours du calage d'historique, on définit d'abord une fonction objectif qui mesure les écarts entre les données d'historique observées et les données d'historique simulées. Plus le nombre de paramètres du modèle de réservoir est grand, plus le calage est difficile, car il faut faire davantage de simulations pour évaluer la fonction objectif afin de trouver un meilleur résultat. Différentes techniques d'optimisation sont développées pour minimiser la fonction objectif. Les approches basées sur la méthode des gradients (Gosselin, O., Cominelli, A. van den Berg, S. et Chowdhury, S.D.: "A Gradient-Based Approach for History Matching of Both Production and 4D Seismic Data", Proceeding 7th European Conference on the Math. of Oil Recovery, Baveno, Italy, 2000) sont largement utilisées dans le calage historique. D'autres approches comme les optimisations stochastiques (Mohamed, L., Christie, M. et Demyanov, V.: "Comparison of Stochastic Sampling Algorithms for Uncertainty Quantification", paper SPE119139 presented at SPE Reservoir Simulation Symposium, Houston, Feb. 2-4, 2009), ou les algorithmes d'apprentissage adaptatif (Feraille, M. et Roggero, F.: "Uncertainty Quantification for Mature Field Combining the Bayesian Inversion Formalism and Experimental Design Approach", 9th European Conf. on Math. of Oil Recovery, Cannes, France, 30 August - 2 Sept. 2004) sont parfois aussi utilisées. Dans toutes ces approches, on applique directement les méthodes d'optimisation à la fonction objectif qui contient l'ensemble des données du réservoir. Ainsi, les techniques de calage d'historique assisté sont développées pour aider les ingénieurs de réservoir à améliorer l'efficacité du calage. Mais ce calage d'historique est généralement un processus long et fastidieux, qui nécessite aussi beaucoup d'efforts et d'expertise de la part des ingénieurs de réservoir. Ces méthodes nécessitent de nombreuses simulations de réservoir pour évaluer la fonction objectif, d'autant plus que le nombre de paramètres est grand. Ces techniques sont donc très coûteuses en temps CPU, et ne sont donc pas toujours adaptées directement aux besoins des spécialistes en charge de l'exploitation du réservoir pétrolier. II est donc important de réduire le nombre de simulations dans le processus d'optimisation. Pour réduire le nombre de simulations, Maschio, C. et Schiozer, D.J.: "A new methodology for assisted history matching using independent objective functions" Petroleum Science and Technology, v 26, n 9, p 1047-1062, June 2008, proposent d'utiliser des fonctions objectifs indépendantes. Dans cette approche, la fonction objectif globale est décomposée en plusieurs fonctions objectif totalement indépendantes, et il est supposé qu'un paramètre qui influence une fonction objectif indépendante ne doit pas avoir d'impact sur d'autres fonctions objectif indépendantes. En pratique, il est très difficile de trouver des cas vérifiant cette hypothèse. De plus, cette méthode consiste à minimiser les fonctions objectif indépendantes, sans considérer la minimisation de la fonction objectif globale. Ainsi, l'objet de l'invention concerne un procédé pour exploiter un réservoir pétrolier à partir d'un modèle de réservoir, dans lequel on réalise un calage d'historique où le nombre de simulations de réservoir est réduit. Le procédé utilise une approche d'optimisation de fonctions objectif locales définies par régions, tout en considérant la cohérence globale du problème en minimisant la fonction objectif globale sur l'ensemble des données du réservoir. La réduction du nombre de simulations de réservoir est obtenue d'une part par une réduction du nombre de paramètres pour les fonctions objectif locales, via un choix de région pour la décomposition, et d'autre part, par une planification des perturbations de paramètres à réaliser pour minimiser la fonction objectif globale. Le procédé selon l'invention L'invention concerne un procédé (ou méthode) pour exploiter un réservoir pétrolier à partir d'un modèle de réservoir constitué d'un ensemble de mailles, auxquelles sont affectées des valeurs d'au moins M paramètres issues de mesures au sein du réservoir pétrolier. Le procédé comporte les étape suivantes : - on définit, en fonction desdits M paramètres, une fonction objectif globale mesurant des écarts entre des valeurs de données d'historique mesurées et des 30 valeurs de données d'historique simulées au moyen dudit modèle de réservoir ; - on décompose ladite fonction objectif globale en une somme de k fonctions objectifs locales, chaque fonction objectif locale mesurant lesdits écarts sur une région géographique dudit réservoir, on définit chaque fonction objectif locale à partir de mk paramètres ayant un impact sur lesdites données d'historique correspondant à chaque région, lesdits mk paramètres étant choisis parmi lesdits M paramètres, et chaque région étant déterminée de façon à minimiser ledit nombre mk de paramètres ; - on modifie des valeurs de paramètres en minimisant ladite fonction objectif globale au moyen d'une méthode de gradient, dans laquelle on estime des dérivées des fonctions objectif locales par rapport aux paramètres mk au moyen d'une technique de perturbation de paramètres basée sur une matrice de perturbations ayant un rang égal au nombre de paramètres mk de ladite fonction objectif locale ; et on utilise ledit modèle ainsi modifié pour exploiter ledit réservoir. Selon l'invention, les mk paramètres peuvent être sélectionnés au moyen d'une analyse de sensibilité pour chaque fonction objectif locale, et les perturbations peuvent être effectuées avec différentes valeurs pour un même paramètre. Les données d'historique sont de préférence des données de production issues de mesures aux puits en réponse à une production du réservoir, et/ou des données de sismiques répétitives. Les données de production peuvent être choisies parmi les données suivantes : production d'huile, production de gaz, production d'eau de puits, ratio gaz/huile (GOR), proportion d'eau de production, pression, température, composition des hydrocarbures. Les M paramètres peuvent être choisis parmi les paramètres suivants : propriétés pétrophysiques telles que perméabilité, porosité, perméabilités relatives et pressions capillaires, propriétés pétroélastiques, distribution spatiale de faciès, endommagement autour de puits, indice de connectivité de puits, connexion ou rejets de failles, paramètres structuraux. Le procédé (ou méthode) de calage d'historique selon l'invention, permet donc de réduire le temps CPU nécessaire au calage des données d'historique, et donc, de diminuer significativement le temps passé dans les études de réservoir, permettant ainsi de répondre aux attentes des ingénieurs réservoirs en charge de l'exploitation d'un réservoir pétrolier. Cette approche est générale, car les interférences entre les fonctions objectifs locales et les paramètres peuvent être prises en compte, puisqu'un paramètre peut avoir un impact sur plusieurs fonctions objectifs locales. D'autres caractéristiques et avantages du procédé 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 des figures La figure 1 montre un réservoir hétérogène avec 4 puits producteurs P,, P2, P3 et P4. Les paramètres K,, K2, K3 et K4 représentent les perméabilités (ou multiplicateurs de 5 perméabilités) autour de ces 4 puits. - La figure 2 montre le réservoir hétérogène avec 25 puits dont 12 injecteurs (I1 à 112) et 13 producteurs (P1 à P13) dans l'exemple. Ce réservoir est maillé avec 50 mailles dans la direction x et 50 mailles dans la direction y. - La figure 3 montre le schéma de perturbations. Les numéros dans les cercles 10 correspondent aux numéros de perturbations. - La figure 4 montre la variation de la fonction objectif dans la procédure d'optimisation avec l'approche conventionnelle. L'axe X correspond au nombre de simulations pour évaluer la fonction objectif, et l'axe Y représente la valeur de la fonction objectif. 15 La figure 5 montre la variation de la fonction objectif dans la procédure d'optimisation avec le procédé selon l'invention. L'axe X correspond au nombre de simulations pour évaluer la fonction objectif, et l'axe Y représente la valeur de la fonction objectif. La figure 6 montre les erreurs sur les paramètres de différents modèles. INI 20 représente le modèle initial. CONVENTION_287 représente le modèle optimal (287ème modèle) avec l'optimisation conventionnelle. INVENTION_55 est le 55ème modèle dans l'optimisation avec le procédé selon l'invention. INVENTION_187 est le modèle optimal (187ème modèle) avec le procédé selon l'invention. L'axe X est le numéro de paramètre, et l'axe Y est l'erreur sur le paramètre (valeur de skin). 25 - La figure 7a montre les comparaisons de production d'huile au puits P2 de différents modèles. REF est la solution de référence, représentant les données mesurées. INI représente le modèle initial. CONVENTION_287 représente le modèle optimal (287ème modèle) avec l'optimisation conventionnelle. INVENTION_55 est le 55ème modèle dans l'optimisation avec le procédé selon l'invention. INVENTION_187 est 30 le modèle optimal (187ème modèle) avec le procédé selon l'invention. L'axe X est le temps en jour. L'axe Y est la production d'huile en m3/jour. La figure 7b montre les comparaisons de water-cut au puits P2 de différents modèles. L'axe X est le temps en jour. L'axe Y est water-cut sans unité. La figure 8a montre les comparaisons de production d'huile au puits P6 de différents modèles. L'axe X est le temps en jour. L'axe Y est la production d'huile en m3/jour. La figure 8b montre les comparaisons de water-cut au puits P6. L'axe X est le temps en jour. L'axe Y est une donnée de water-cut, sans unité. - La figure 9a montre les comparaisons de production d'huile au puits P12. L'axe X est le temps en jour. L'axe Y est la production d'huile en m3/jour. La figure 9b montre les comparaisons de water-cut au puits P12. L'axe X est le temps en jour. L'axe Y est le water-cut sans unité. - La figure 10 montre les comparaisons d'injection d'eau puits 16 de différents 10 modèles. L'axe X est le temps en jour. L'axe Y est la production d'eau en m3/jour. - La figure 11 montre les comparaisons d'injection d'eau puits 112 de différents modèles. L'axe X est le temps en jour. L'axe Y est la production d'eau en m3/jour. Description détaillée du procédé 15 Le procédé (ou méthode) selon l'invention permet d'exploiter un réservoir pétrolier à partir d'un calage d'historique optimisé d'un modèle de réservoir. La calage d'historique est optimisé en ce qu'il permet de réduire significativement le nombre de simulations de réservoir, qui sont très coûteuses en temps CPU, ainsi que le temps passé dans une étude de réservoir. 20 Dans la présente invention, on propose de décomposer le calage par régions. Cette décomposition représente un partitionnement des données à caler. Les régions sont définies de telle sorte qu'une partie seulement des paramètres ait un impact significatif sur les données d'historique correspondant à chaque région. En procédant de la sorte, le calage peut être réalisé de façon relativement indépendante dans 25 chaque région, tout en obtenant une amélioration globale. Le nombre de paramètres influents étant réduit par région, il est ainsi possible de limiter fortement le nombre de simulations nécessaires au calage. Pour assurer la cohérence du calage sur l'ensemble du réservoir, on associe une technique de perturbation de paramètres qui tient compte des interférences entre les régions. Le procédé comporte principalement 30 quatre étapes : 1. on génère un modèle de réservoir initial et on choisit les paramètres à modifier pour améliorer ce modèle ; 2. on définit une fonction objectif globale pour le calage d'historique et on décompose cette fonction objectif en composantes locales par région, basées sur le partitionnement des données ; 3. on planifie de façon optimale les simulations à réaliser, en perturbant les valeurs de paramètres, pour calculer les gradients de la fonction objectif ; 4. on construit un nouveau modèle de réservoir en modifiant les paramètres de façon à minimiser la fonction objectif par une procédure d'optimisation basée sur la méthode des gradients. 1- Génération d'un modèle initial de réservoir On commence par générer un modèle de réservoir initial, qui prend en compte différents aspects tels que la structure géologique, les propriétés pétrophysiques, les propriétés fluides, les puits, etc. Un modèle de réservoir est constitué d'un ensemble de mailles, auxquelles sont affectées des valeurs d'au moins M paramètres issues de mesures au sein du réservoir pétrolier. Les paramètres, qui sont utilisés pour décrire la structure, la pétrophysique, etc., sont estimés en utilisant des données statiques. Certains paramètres ont de grandes incertitudes et peuvent être éloignés de la réalité. Néanmoins, nous pouvons utiliser les données dynamiques telles que les données de production ou de sismique 4D pour réduire les incertitudes et améliorer la caractérisation de ces paramètres via la procédure de calage d'historique. Dans cette étape, nous choisissons tous les paramètres potentiels à modifier pour améliorer le modèle de réservoir. La génération d'un modèle de réservoir initial est bien connue des spécialistes. 2- Définition et décomposition d'une fonction obiectif globale La fonction objectif globale, qui mesure des écarts entre des valeurs de données d'historique mesurées et des valeurs de données d'historique simulées au moyen du modèle de réservoir, est généralement donnée par la formule suivante : Nwell Ntime (r) obs srm 1 P n ùPin (~) F(2) _ 2 %,i Ntime(i) 6p(i) i=1 n=1 Nzone 1 Nrime (J) s>obs ùsj n simn (4)2 ~( +ù ws 2 J Ntime(j) 6 (J) 1 s j= où NWeä est le nombre de puits Nzone est le nombre de zone sismiques 4D, P°bs est l'ensemble des données de production observées, Ps'm est l'ensemble des données de production simulées avec le modèle, es représente les données de sismique 4D observées, Ss'm représente les données de sismique 4D simulées avec le modèle, wp sont les poids associés aux données de production, ws sont les poids associés aux données sismiques 4D, 6p(i) est l'intervalle de confiance sur les données de production au puits i, o (j) est l'intervalle de confiance sur les données sismiques à la zone j, Ntme représente le nombre de point en temps pour effectuer les mesures, 2 = (21, 22, ..., 2M) sont les M paramètres à optimiser. Les données de production peuvent être choisies parmi les données suivantes : production d'huile, production de gaz, production d'eau de puits, ratio gaz/huile (GOR), proportion d'eau de production, pression, température, composition des hydrocarbures. Cette liste n'est bien sur non exhaustive. 20 Les paramètres peuvent être choisis parmi les paramètres suivants : propriétés pétrophysiques telles que perméabilité, porosité, perméabilités relatives et pressions capillaires, propriétés pétroélastiques, distribution spatiale de faciès, endommagement autour de puits, indice de connectivité de puits, connexion ou rejets de failles, paramètres structuraux. Cette liste n'est bien sur non exhaustive. )2 n=1 10 15 25 Cette fonction objectif globale peut être directement décomposée sur le partitionnement des données par région sous la forme suivante : Nreg F(21,...,/,M)= EFk(~11,...,2M) k=1 où Fk(1 ..., /,M) est une composante (fonction objectif locale) qui estime les écarts entre les données observées et les données simulées sur les puits et les zones sismiques 4D dans la région k (k = 1, ..., N1eg), et Nfeg est le nombre de régions. A noter que la fonction objectif locale dépend théoriquement toujours de M paramètres, mais les sensibilités de ces paramètres sur chaque fonction objectif locale peuvent être très différentes. Nous pouvons ainsi réduire, par une analyse de sensibilité par exemple, le nombre de paramètres principaux (significatifs) pour les fonctions objectif locales. Ainsi, en partitionnant convenablement les composantes de la fonction objectif, c'est-à-dire en définissant convenablement les régions, et les paramètres par régions, nous pouvons réduire au minimum le nombre de paramètres principaux pour les fonctions objectif locales. Notons que le principe de partitionnement de la fonction objectif, illustré ci-dessus avec une fonction objectif de type "moindres carrés", peut être généralisé à d'autres types de formulation. Cette méthode peut être illustrée sur un exemple de réservoir hétérogène comportant 4 puits producteurs (figure 1). Nous supposons sur ce cas avoir une grande incertitude sur les perméabilités. Les productions de ces 4 puits producteurs P,, P2, P3 et P4 dépendent, bien entendu, des valeurs des perméabilités du champ. Nous voulons modifier les valeurs des perméabilités (ou des multiplicateurs de perméabilités) K,, K2, K3 et K4 autour de ces 4 puits par minimisation de la fonction objectif pour améliorer le modèle. Sur cet exemple, nous pouvons décomposer la fonction objectif par puits de la façon suivante : 4 F(KI,K2,K3,K4)= EFk(KI,K2,K3,K4) k=1 où Fk(K,,K2,K3,K4) correspond aux écarts entre données de production observées et simulées pour le puits Pk (K = 1, 2, 3, 4). La production du puits P, dépend principalement de la valeur du paramètre K,, et les impacts des autres paramètres (K2, K3 et K4) sur la production du P, sont très faibles. La dépendance des paramètres sur la fonction objectif locale FI peut être approchée par : F1(K1,K2,K3,K4) F1(K1) De même, nous avons, F2(K1,K2,K3,K4) z'F2(K2) F3(K1,K2,K3,K4) F3(K3) F4(K1,K2,K3,K4) F4(K4) Dans les cas plus complexes, les dépendances de paramètres peuvent être 10 déterminées par une analyse de sensibilité. Cette technique est bien connue des spécialistes, un exemple est décrit dans le document suivant : A. Saltelli, S. Tarantola, F. Campolongo. Sensitivity Analysis in Practice : A Guide to Assessing Scientific Models. Wiley, 2004 Ainsi, la décomposition de la fonction objectif globale selon l'invention conduit à 15 ce que chaque composante (fonction objectif locale) dépende du plus petit nombre possible de paramètres : Fk(21,22,...,/1M)Fk(2k1,2k2,...,2kmk ) avec mk le plus petit possible (k = 1, ..., Nreg). Soit Mm = max(m,, ..., m4), la valeur maximum parmi m,, ..., mNreg, l'objectif de cette étape est de chercher une 20 décomposition de la fonction objectif telle que Mn, soit le plus petit possible avec Mm D.F = aF de la fonction objectif F(X1, ..., ..., %,M) par rapport aux paramètres X; (i = a/îi 1, ..., M) sont calculées au moyen d'une technique de perturbation de paramètres, au 5 cours de laquelle on évalue M+1 valeurs de la fonction objectif. Par exemple, la dérivée DiF = aF peut être calculée par : aa,i DiF(21,...,2i, 1M) F(21,...,2i +zl/li,...,/,M)ùF(21,...,2i,...,/lM) A2i (i=1,...,M) Pour connaître toutes les M dérivées, il faut évaluer la fonction objectif F sur le 10 point initial (X, ..., XM) et pour M points perturbés (~.1,..., 2i + ~M) pour i = 1, ..., M, c'est-à-dire, faire M+1 simulations de réservoir. Or la simulation de réservoir est très coûteuse en temps CPU. De même, pour une fonction objectif locale Fk (k = 1, ..., Nfeg), nous avons besoin de mk perturbations (ou mk+1 simulations) pour calculer ses dérivées. Si Mn, 15 à-dire, mk objectif locale de façon à pouvoir calculer ses dérivées par rapport à tous ses paramètres. Pour réaliser cette planification, on construit une matrice, appelée matrice de perturbations. Cette matrice contient autant de colonnes que de paramètres utilisés pour définir la fonction objectif globale, et autant de lignes que de perturbations à réaliser. Chaque case de la matrice comporte la variation du paramètre à réaliser. Cette variation est également appelée � pas de perturbation � et est notée a;,. Ainsi, dans le tableau 1, qui correspond à une matrice de perturbations, a;, correspond à la variation du paramètre X pour la perturbation i, L étant le nombre de perturbations. Tableau la - Planification des perturbations perturbation X1 X2 X3 X4 XM 1 a11 a12 a13 a14 a1M 2 a21 a22 a23 a24 a2M 3 a31 a32 a33 a34 a3M L au aL2 aL3 aL4 aLM Une planification classique correspond à L = M, = 0 (si i ~ j) et a;; = (i=1, L;j=1,...,M): Tableau 1 b - Planification des perturbations classique perturbation X2 X3 X4 • • • XM 1 A.i,l 0 0 0 0 2 0 d22 0 0 0 3 0 0 d~,3 0 0 L(=M) 0 0 0 0 A/lm Il faut donc, dans un cas classique, réaliser L (L=M) perturbations, pour appliquer 15 une variation indépendamment à chaque paramètre, et pouvoir ainsi calculer les M dérivées de la fonction objectif locale. Planification des perturbations Selon l'invention, on applique une technique de planification des perturbations dans laquelle on applique une variation à plusieurs paramètres en même temps, de façon à minimiser le nombre de perturbations (nombre de lignes de la matrice de perturbation) tout en permettant de calculer les M dérivées de la fonction objectif locale. Pour ce faire, on étudie le nombre minimum de perturbations et les pas des perturbations a;; nécessaires au calcul des dérivées d'une fonction objectif locale. Nous prenons un exemple dont la fonction objectif locale F, a 3 paramètres principaux. Sans perte de généralité, nous supposons que les 3 paramètres principaux de la fonction objectif locale F, sont 2 ,1, 22 et 2,3. Soit (29,202,203) le point initial, en utilisant le développement de Taylor, la fonction objectif locale F, correspondant à la perturbation i peut s'exprimer par : 3 FI(29 +ai1,+ai2,2 +ai3) F1(201,22,/103)+ .FI(201 '222'/103 )a j=1 (i=1,...,L) Dans le système ci-dessus, D;FI (j = 1, 2, 3) sont inconnues, et tous les autres termes sont connus en utilisant les simulations de réservoir sur le point initial .102 , 2o) et les points perturbés (29+ail,22+ai2,2 +ai3) (j = 1, ..., L). Pour calculer les dérivées D;Fi (j = 1, 2, 3), nous considérons le système d'équations suivant : 3 EaijD1F1 =F1(21 +ail, /12 +ai2,a,0 +ai3)ùFI(~i,9,I.o,/L3) j=1 (i=1,...,L) Il existe une infinité de solutions si le nombre de perturbations est inférieur au nombre de paramètres qui est fixé à 3 dans cet exemple. Il faut donc au minimum 3 25 perturbations, et les coefficients a;; des perturbations doivent vérifier que le rang de la 14 20 15 matrice C= a11 a12 a13 a21 a22 a23 est 3. Cette matrice est appelée la matrice de aL1 aL2 aL3i perturbation de la fonction objectif locale F,. Nous pouvons utiliser différents pas de perturbations a;; pour le même paramètre 4; dans différentes perturbations i (i = 1, ..., L), afin que le rang de C soit égale au 5 nombre de paramètres mk. Il est important de minimiser le nombre de perturbations L. Pour chaque fonction objectif locale, il est possible d'utiliser seulement mk perturbations pour calculer ses dérivées. L'objectif de cette étape est d'arranger les coefficients des perturbations pour minimiser le nombre de perturbations L. En particulier, nous souhaitons avoir L = Mn,. Calcul des dérivées Pour calculer les dérivées D;F, (j = 1, 2, 3) dans l'exemple précédent, on résout le problème de minimisation suivant : L 3 trouver (x,, x2, x3) tel que la fonction G(xl,x2,x3) _ E(Eaijx j ùbi )2 avec i=1 j=1 15 bi =F1(2 +ail,.1,02 +ai2,23 +ai3)ùF1(21 ,202 ,.,3) (i = 1, 2, 3) atteigne le minimum. Les x,, x2, x3 correspondent aux dérivées D;F, (j = 1, 2, 3) que l'on souhaite calculer. D'une manière générale, pour déterminer les dérivées de la fonction objectif locale Fk qui dépend de mk paramètres principaux, il faut au moins mk perturbations et 20 le rang de la matrice de perturbation de Fk est mk. Soit C la matrice de perturbation de la fonction objectif Fk, C est une matrice de L lignes et mk colonnes avec L mk . Le rang de C doit être mk, et les dérivées de la fonction objectif locale Fk peuvent être calculées en résolvant le problème de minimisation suivant : L mk Trouver (xl,...,xmk) tels que la fonction G(xl,...,xmk )= E(Icijx j ùbi )2 i=1 j=1 25 soit minimale, où c;; (i = 1, ..., L ; j = 1, ..., mk) sont les coefficients de la matrice de 10 perturbation C, bi =Fk(P)ûFk(Po) (i = 1, ..., L) avec Po le point initial et Pi le point correspondant à la perturbation i. Ainsi, on estime les dérivées des fonctions objectif locales par rapport aux paramètres mk, au moyen d'une technique de perturbation de paramètres au cours de laquelle on évalue R valeurs de la fonction objectif locale. La valeur de R est diminuée en construisant pour chaque fonction objectif une matrice de perturbations ayant un rang égal au nombre de paramètres mkde la fonction objectif locale. Dans l'exemple de réservoir hétérogène avec 4 puits ci-dessus, nous pouvons 10 utiliser une perturbation (2 simulations) au lieu de 4 perturbations (5 simulations) pour calculer toutes les dérivées. 4- Construction du nouveau modèle par optimisation Le modèle de réservoir initial est ainsi mis à jour en modifiant les paramètres de 15 ce modèle. Les nouveaux paramètres étant ceux qui minimisent la fonction objectif globale. Plusieurs méthodes d'optimisation basées sur le calcul de gradients sont disponibles (voir D. Sinoquet et F. Delbos :"Adapted Nonlinear Optimization Method for Production Data and 4D Seismic Inversion", proceeding of European Conference on Mathematics of Oil Recovery, Bergen, Norvège, 8-11 Sept. 2008) pour le calage 20 d'historique. Nous pouvons mettre en oeuvre une méthode d'optimisation en utilisant la décomposition de la fonction objectif et la planification de perturbations décrites dans les étapes précédentes pour minimiser la fonction objectif (globale) en modifiant les paramètres. Nous pouvons ainsi obtenir un nouveau modèle de réservoir avec les 25 paramètres optimaux qui minimisent la fonction objectif. Le modèle ainsi calé aux données d'historique permet une optimisation de l'exploitation du réservoir pétrolier. Ce modèle permet de déterminer de façon plus précise et plus rapide de nombreux paramètres techniques relatifs à la recherche, l'étude ou l'exploitation d'un réservoir, d'hydrocarbures par exemple. On peut ainsi par 30 exemple déterminer quelles zones ont le plus de chances de contenir des hydrocarbures, les zones dans lesquelles il peut être intéressant/nécessaire de forer un puits d'injection pour améliorer la récupération des hydrocarbures, le type d'outils à utiliser, les propriétés des fluides utilisés et récupérés.... Exemple d'application Soit un réservoir hétérogène de taille 2500 m dans la direction x, 2500 m dans la direction y et 10 m dans la direction z. Ce réservoir est maillé avec 50 mailles dans la direction x et 50 mailles dans la direction y, et ses perméabilités (K) sont présentées sur la figure 2. Vingt cinq puits sont forés dans ce réservoir dont douze puits injecteurs et treize puits producteurs. La pression initiale du réservoir est 250 bars. Les puits fonctionnent à pressions imposées avec 320 bars aux puits injecteurs et 180 bars aux producteurs. Les skins (les endommagements ou améliorations de perméabilités autour du puits) sont nuls pour tous les puits. Un historique de 4000 jours est disponible. Les données de production mesurées sont : les débits d'huile aux 13 puits producteurs, le débit d'eau aux 12 puits injecteurs et les water-cuts aux 13 puits producteurs. La fonction objectif globale, qui mesure les écarts entre les données d'historique simulées et les données d'historique mesurées, est définie de la façon suivante : F(2) = (Nprod Ntime (Qoo b n sim (~))2 + 2 i=1 Ntime n=1 6° t Nprod Ntime obs _ sim j (FWi n FWi n 0) i=1 n=1 Ninj Ntime obs sim I (Qw,i,n ùQw,i,n(2))2) i=1 Ntime n=1 6w(t) où Nprod = 13 est le nombre de puits producteurs, Nin; = 12 est le nombre d'injecteurs, Ntime = 40 est le nombre de points mesurés avec une mesure tous les 100 jours, Q0 est le débit d'huile aux producteurs, QW est le débit d'eau aux injecteurs, FW est le water-cut, 25 60, 6f,,,, et a w, qui sont les intervalles de confiance pour le débit d'huile, le watercut et le débit d'eau, ont été calculés en fonction de la moyenne et la variance des données à chaque puits sur cet exemple, X = (X1, •••, X 13, 2i1, •••, X 12) correspond aux valeurs de skins aux puits, les premières valeurs 41, ..., X13 étant attachées aux treize puits producteurs et 2^,11, ..., 2412 aux douze puits injecteurs. Cette fonction objectif globale est décomposée en vingt cinq fonctions objectifs locales, liées aux producteurs et aux injecteurs : 25 F(ap1,...,4'13,211,...,2112) = Fk(2PI,...,2P13,/,I1,....2112) k=1 avec Fk(2P1 1'13,/I71,...,2112)= Ntime obs sim F.Wobs _ FWsim 1 ((Qo,k,n ùQo,k,n )2 +( k,n k,n )2) 2 Ntime 71=1 ao(i) 6fw(i) (k = 1, ..., 13) correspondant aux 13 producteurs, et Ntime obis; sim Fk(SPI 1 1 Qw,kù13,n ùQw,kù13,n 2 ~P13 ~I1, 2112 ( (k = 14, 2 Ntime n=1 6w(i) ..., 25) correspondant aux 12 puits injecteurs. Dans la construction du modèle de réservoir initial, nous ne connaissons pas les véritables valeurs des skins. Nous adoptons d'abord des valeurs initiales de skins estimées à 10 pour tous les puits. Nous supposons que les valeurs possibles de skins sont comprises dans un intervalle entre 0 et 20. La localisation des puits dans ce réservoir (figure 2) nous montre que le skin d'un puits injecteur a un fort impact uniquement sur le puits lui-même et sur les producteurs environnants (maximum 4 producteurs), et son impact sur les autres puits doit être relativement faible. De même, le skin d'un puits producteur a un fort impact seulement sur le puits lui-même et sur les injecteurs environnants, et peu d'impact sur les autres puits. Donc, pour chaque fonction objectif locale, nous avons au maximum 5 paramètres principaux (le skin sur le puits considéré dans la fonction objectif locale et les skins sur les puits voisins). Les perturbations à réaliser sur les paramètres, qui doivent garantir le rang de la matrice de perturbation de chaque fonction objectif locale égale au nombre de paramètres principaux de cette fonction objectif locale, peuvent être planifiées de la façon suivante (figure 3) : Perturbation 1 : on perturbe les puits dans les colonnes 1 et 4; Perturbation 2 : on perturbe les puits dans les colonnes 2 et 5; Perturbation 3 : on perturbe les puits dans les lignes 1 et 4; Perturbation 4 : on perturbe les puits dans les lignes 2 et 5; Perturbation 5 : on perturbe les puits dans la ligne 3 et la colonne 3; Ces perturbations sont résumées dans le Tableau 2. Les cellules avec une croix correspondent aux perturbations avec des coefficients non nuls. Les cellules non remplies signifient qu'il n'y a pas de perturbation pour les paramètres correspondants, ce qui se traduit par des coefficients nuls pour ces cellules. Tableau 2 - Plan de perturbations pertur- ~`P1 %P2 !~P3 %P4 %P5 XP6 XP7 !P8 %P9 !~P10 XP11 %P12 !~P13 %11 XI2 XI3 XI4 XI5 %I6 XI7 %18 X19 X110 X111 %112 bation 1 X X X X X X X X X X 2 X X X X X X X X X X 3 X X X X X X X X X X 4 X X X X X X X X X X 5 X X X X X X X X X La fonction objectif globale est minimisée en utilisant une procédure d'optimisation basée sur la méthode des gradients pour optimiser les 25 paramètres. La valeur de fonction objectif pour le modèle initial est 1031. La variation de la fonction objectif avec l'approche conventionnelle, qui calcule les dérivées avec 25 perturbations, est illustrée dans la figure 4. L'axe X correspond au nombre de 25 simulations. Au total, 287 simulations sont effectuées, et la fonction objectif optimale est 422 pour la 287ième simulation (on note ce modèle par CONVENTION_287). La fonction objectif en fonction du nombre de simulations avec la méthode proposée dans cette invention est illustrée dans la figure 5. Au total, 197 simulations sont effectuées, et la valeur minimum de 54 est obtenue pour la 187ième simulation (on note ce modèle par INVENTION_187). Nous notons aussi que la fonction objectif est réduite à 418 à la 55ième simulations (on note ce modèle par INVENTION_55), ce qui signifie que nous pouvons obtenir un modèle de qualité identique (mêmes valeurs de fonctions objectifs) par rapport au modèle optimal obtenu avec l'approche conventionnelle (CONVENTION_ 287) avec 5 fois moins de simulations. Ce résultat est cohérent avec l'analyse théorique, car nous utilisons 5 perturbations au lieu de 25 pour l'approche conventionnelle et nous réduisons d'un facteur 5 le nombre de perturbations (simulations) pour le calcul des dérivées. Avec la méthode proposée, le gain en temps CPU et en temps ingénieur consacré à l'étude d'un réservoir est conséquent. Les erreurs absolues sur les paramètres (skins) sont présentées sur la figure 6 pour le modèle initial et pour quelques modèles améliorés (CONVENTION_287, INVENTION_55 et INVENTION_187). Les modèles CONVENTION_287 et INVENTION_55 donnent des résultats assez proches et réduisent les erreurs par rapport au modèle initial. Cependant, les réductions d'erreurs sont limitées. Avec le modèle optimal du procédé selon l'invention INVENTION_187, les résultats sont nettement meilleurs. Les figures 7a à 11 montrent les résultats des simulations sur quelques puits (P2, P6, P12 et 16, 112). Nous pouvons considérer que les résultats sont équivalents pour les modèles CONVENTION_287 et INVENTION_55, qui réduisent légèrement les erreurs. Néanmoins, le modèle optimal obtenu par le procédé selon l'invention INVENTION_187 donne toujours des résultats nettement meilleurs. Avec ce procédé, nous pouvons obtenir un meilleur modèle avec moins de simulations. REVENDICATIONS1. Procédé pour exploiter un réservoir pétrolier à partir d'un modèle de réservoir constitué d'un ensemble de mailles, auxquelles sont affectées des valeurs d'au moins M paramètres issues de mesures au sein dudit réservoir pétrolier, caractérisée en ce que l'on réalise les étape suivantes : - on définit, en fonction desdits M paramètres, une fonction objectif globale mesurant des écarts entre des valeurs de données d'historique mesurées et des valeurs de données d'historique simulées au moyen dudit modèle de réservoir ; - on décompose ladite fonction objectif globale en une somme de k fonctions objectifs locales, chaque fonction objectif locale mesurant lesdits écarts sur une région géographique dudit réservoir, on définit chaque fonction objectif locale à partir de mk paramètres ayant un impact sur lesdites données d'historique correspondant à chaque région, lesdits mk paramètres étant choisis parmi lesdits M paramètres, et chaque région étant déterminée de façon à minimiser ledit nombre mk de paramètres ; - on modifie des valeurs de paramètres en minimisant ladite fonction objectif globale au moyen d'une méthode de gradient, dans laquelle on estime des dérivées des fonctions objectif locales par rapport aux paramètres mk au moyen d'une technique de perturbation de paramètres basée sur une matrice de perturbations ayant un rang égal au nombre de paramètres mk de ladite fonction objectif locale ; et - on utilise ledit modèle ainsi modifié pour exploiter ledit réservoir. 2. Procédé selon la revendication 1, dans lequel lesdits mk paramètres sont sélectionnés au moyen d'une analyse de sensibilité pour chaque fonction objectif locale. 3. Procédé selon l'une des revendications précédentes, dans lequel les perturbations sont effectuées avec différentes valeurs pour un même paramètre. 4. Procédé selon l'une des revendications précédentes, dans lequel les données d'historique sont des données de production issues de mesures aux puits en réponse à une production du réservoir, et/ou des données de sismiques répétitives. 5. Procédé selon la revendication 4, dans lequel les données de production sont choisies parmi les données suivantes : production d'huile, production de gaz, production d'eau de puits, ratio gaz/huile (GOR), proportion d'eau de production, pression, température, composition des hydrocarbures. 6. Procédé selon l'une des revendications précédentes, dans lequel les M paramètres sont choisis parmi les paramètres suivants : propriétés pétrophysiques telles que perméabilité, porosité, perméabilités relatives et pressions capillaires, propriétés pétroélastiques, distribution spatiale de faciès, endommagement autour de puits, indice de connectivité de puits, connexion ou rejets de failles, paramètres structuraux.