FR 3059443 A1 20180601 FR 1661634 A 20161129 1 La présente invention concerne le domaine de l'exploration et de l'exploitation des gisements pétroliers ou de sites de stockage géologique de gaz. L'exploration et l'exploitation d'un gisement, notamment pétrolier, nécessitent d'acquérir une connaissance aussi précise que possible de la géologie souterraine, et ce, afin de permettre une estimation des réserves, une estimation de la production, ou encore la définition de la gestion de l'exploitation du gisement dans le temps. En effet, la détermination de l'emplacement d'un puits de production ou d'un puits d'injection au sein d'un gisement d'hydrocarbures, la constitution de la boue de forage, les caractéristiques de complétion, le choix d'un procédé de récupération des hydrocarbures (tel que l'injection d'eau par exemple) et des paramètres nécessaires à la mise en oeuvre de ce procédé (tels que la pression d'injection, le débit de production, etc.) nécessitent de bien connaître le gisement. Pour acquérir cette connaissance, l'industrie pétrolière allie les mesures sur champ (réalisées in situ, lors de campagnes sismiques, de mesures dans des puits, de carottages etc.) aux mesures réalisées en laboratoire ainsi qu'aux simulations numériques (réalisées au moyen de logiciels, tels qu'un simulateur d'écoulement). La formalisation de cette connaissance passe ensuite par l'établissement d'une maquette, aussi appelée modèle géologique, visant à rendre compte, de façon approchée, à la fois de la structure et du comportement d'un gisement. Ce type de maquette est généralement représenté sur un ordinateur, sous la forme d'une représentation maillée (qui peut être une grille régulière, ou bien un maillage structuré ou non structuré de façon plus générale), chacune des mailles de la représentation maillée comportant au moins une et préférentiellement plusieurs valeurs de propriétés pétrophysiques relatives au gisement lui-même (telles que porosité, perméabilité, faciès etc.) et/ou de propriétés relatives aux fluides présents dans le gisement (saturations en eau/huile/gaz, perméabilités relatives, etc.). Lorsqu'une représentation maillée respectant l'ensemble des données mesurées sur le terrain est finalement obtenue, cette représentation est utilisée pour prédire les déplacements de fluide dans la formation étudiée et planifier son développement futur, via la détermination de schémas d'exploitation. Pour ce faire, le spécialiste est amené 1) à sélectionner les zones d'intérêt (présentant les niveaux d'hydrocarbures les plus élevés par exemple) dans lesquelles forer de nouveaux puits, 2) à placer ces puits précisément dans les trois dimensions de l'espace, 3) à définir au moins un procédé de récupération (injection d'eau, de gaz, de vapeur etc.) et 4) à vérifier, par exemple via une simulation 3059443 2 d'écoulement, si la production d'hydrocarbures, avec cette implantation des puits et ce procédé de récupération, est optimale. Une simulation d'écoulement, aussi appelée simulation de réservoir, est une technique permettant de simuler les écoulements de fluides au moyen d'un logiciel appelé simulateur d'écoulement et d'une représentation 5 maillée de la formation étudiée. L'élaboration d'une représentation maillée réaliste en ce qui concerne les écoulements comprend généralement une étape de calage d'historique (« history matching » en anglais). Le calage d'historique consiste à modifier certains des paramètres 10 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. 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, 15 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.). Au cours du calage d'historique, on définit une fonction objectif qui mesure les écarts entre les données d'historique observées et les données d'historique simulées. Puis un processus de minimisation de la fonction objectif est mis en place. Cette 20 minimisation est généralement pilotée par un algorithme d'optimisation, basé sur un calcul de gradients de la fonction objectif. D'autres techniques permettent d'échantillonner l'espace des modèles respectant les données de production plutôt que de chercher à atteindre le modèle le plus probable parmi ceux-ci (par exemple, les algorithmes génétiques ou le recuit-simulé). 25 État de la technique Les documents suivants seront cités au cours de la description : 30 Barbier, M., Hamon, Y., Doligez, B., Callot, J.-P., Floquet, M., Daniel, J.-M., 2012, Stochastic joint simulation of facies and diagenesis,: A case study on early diagenesis of the Madison formation (Wyoming, USA), OGST, 67(1), 123-145 Caers, J., Hoffman, T., The probability perturbation method: A new look at Bayesian 35 inverse modeling, Math. Geol., 38(1), 81-100, 2006, D01:10.1007/s11004-005-9005-9 3059443 3 Doligez, B., Hamon, Y., Barbier, M., Nader, F., Lerat, 0., 2011, Advanced workflows for joint modelling of sedimentary facies and diagenetic overprint. Imact on reservoir quality, SPE 146621, SPE ATCE, Denver, CO, USA Gosselin, 0., Cominelli, A. van den Berg, S. et Chowdhury, S.D.: "A Gradient-Based 5 Approach for History Matching of Both Production and 4D Seismic Data", Proceeding 7th European Conference on the Math. of Oil Recovery, Baveno, Italy, 2000) Guardiano, F., Srivastava, M., 1993, Multivariate geostatistics: beyond bivariate moments, In: Soares, A. (Ed.), Geostatistics-Troia. Kluwer Academic Publications, Dordrecht, 133144, http://dx.doi.orq/10.1007/978-94-011-1739-5 12. 10 Holland, J.H., 1975, Adaptation In Natural And Artificial Systems, University Of Michigan Press, Ann Arbor Hu, L.-Y., 2000, Gradua) deformation and iterative calibration of Gaussian-related stochastic models, Mathematical Geology, 32(1), 87-108 Le Ravalec-Dupin, M., Da Veiga, S., Cosimulation as a perturbation method for calibrating 15 porosity and permeability fields to dynamic data, Computers and Geosciences, 37, 1400-1412, 2011, D01:10.1016/j.cageo.2010.10.013 Le Ravalec, M., Doligez, B., Lerat, 0., 2014, Integrated reservoir characterization and modeling,http://books.ifpeneraiesnouvelles.fr/ebooks/integrated reservoir characterizatio n nd modelire, doi: 10.2516/ifpen/2014001 20 Le Ravalec, M., Hu, L.-Y., Noetinger, B., 1999, Stochastic Reservoir Modeling Constrained to Dynamic Data: Local Calibration and Inference of the Structural Parameters, SPE 56556, SPE ATCE, Houston, TX, USA Pontiggia, M., Ortenzi, A., Ruvo, L., 2010, New integrated approach for diagenesis characterization and simulation, SPE 127236, SPE North Africa Technical Conference 25 and Exhibition, Cairo, Egypt Press, W.H., Teufolsky, S.A., Vetterling, W.T., Flannery, B.P., 1988, Numerical recipes in C, Cambridge University Press, New York, NY, USA Renard, D., Beucher, H., Doligez, B., 2008, Heterotopic bi-categorical variables in pluriGaussian truncated simulation, VIII Int. Geostat. Congress, GEOSTATS 2008, 1, pp. 28930 298 Roggero, F., Enchery, G., Ponsot-Jacquin, C., 2011, Method for gradually modifying lithologic facies proportions of a geological model, Patent number 7970593 Tarantola, A., 1987, Inverse Problem Theory - Methods For Data Fitting And Model Parameter Estimation. Elsevier Science Publishers, Amsterdam, The Netherlands, 613 pp 3059443 4 Tillier, E., Enchery, G., Le Ravalec, M., Gervais, V., 2010, Local and continuous updating of facies proportions for history-matching, Journal Petroleum Science and Engineering, 73(3-4), 194-203 5 On connaît des techniques calage d'historique pour lesquelles les paramètres d'ajustement sont relatifs aux faciès lithologiques rencontrés dans la formation souterraine étudiée. On peut en particulier citer la technique décrite dans le document EP 1978379 Al 10 (US 7970593 B2) qui consiste à modifier les proportions de faciès dans chacune des mailles de la représentation maillée représentative de la formation étudiée de façon à minimiser un écart entre données de production simulées et données de production mesurées, des propriétés d'écoulement (porosité, perméabilité par exemple) étant attribuées à chacun des faciès lithologiques considérés. Plus précisément, cette 15 technique décrit un ajustement itératif des proportions de faciès tout en respectant certaines règles de proportionnalité des faciès définies par l'ingénieur. On connaît également la technique décrite dans le document (Le Ravalec et al., 1999) qui consiste à ajuster la variabilité spatiale (paramétrée par la portée d'un variogramme) des faciès au sein de la représentation maillée de façon à minimiser un 20 écart entre données de production simulées et données de production mesurées. Par contre, il apparaît qu'à ce jour, les effets des processus diagénétiques ne sont pas pris en compte dans la mise à jour du modèle de réservoir par calage d'historique. Or, la diagénèse figure parmi les processus physico-chimiques majeurs de l'histoire d'un 25 bassin sédimentaire. La diagenèse englobe tous les processus de transformation chimique, biochimique et physique des sédiments s'étant déposés dans un bassin, en roches sédimentaires compactes. En effet, les sédiments qui se déposent dans un bassin sédimentaire sont meubles et riches en eau. Les processus affectent ces sédiments dès leurs dépôts et se poursuit tout au long de leur enfouissement, c'est-à-dire, à mesure que 30 d'autres sédiments viennent recouvrir les recouvrir progressivement sous plusieurs dizaines, centaines ou même milliers de mètres de matériel. À de telle profondeur, les processus diagénétiques sont influencés par les conditions de pression et de température. Ces transformations se déroulent selon différentes étapes, appelées aussi phases 35 diagénétiques, qui varient selon la nature des sédiments, de la nature des fluides circulant entre ces sédiments et des conditions d'enfouissement. 3059443 5 Ces processus comprennent entre autres la compaction (tassement mécanique lié au poids des terrains déposés par-dessus le sédiment ; ce processus tend à réduire la porosité dans la roche et à augmenter les points de contact entre les grains), la dissolution des grains à leurs points de contact, ou encore la cimentation (précipitation de 5 cristaux d'un minéral donné dans les espaces inter-granulaires liant les particules sédimentaires entre elles). L'analyse de lames minces extraites d'une carotte prélevée au sein d'une formation souterraine permet d'étudier l'histoire diagénétique d'une roche. Il est notamment question de définir la succession dans le temps de l'ensemble des transformations 10 minéralogiques ayant affecté les sédiments depuis leur dépôt. Il est possible de qualifier et quantifier la proportion de chacune de ces transformations minéralogiques dans un échantillon donné soit visuellement, en définissant par exemple un type et/ou un degré de diagenèse (correspondant à une proportion faible, moyenne ou importante d'une phase diagénétique donnée), soit en laboratoire par des techniques de laboratoire (DRX, 15 géochimie élémentaire, etc.), ou encore par analyse d'image en utilisant un logiciel de traitement d'image. Mais cette information est parcellaire, en quelques points le long du puits duquel a été extraite la carotte. L'extrapolation de ces informations parcellaires, entre les puits, est donc incertaine. Or les processus physico-chimiques subis par un dépôt sédimentaire au cours de la diagénèse affectent significativement les propriétés 20 d'écoulement de la roche résultant de ce dépôt sédimentaire. Par conséquent, il peut se produire que les représentations maillées utilisées dans les étapes classiques d'exploration et d'exploitation pétrolière, notamment lors d'un calage d'historique, soient basées sur des hypothèses erronées concernant l'état diagénétique des différents faciès présents dans une formation. 25 La présente invention vise à la prise en compte de l'état diagénétique d'une roche comme l'un des paramètres d'ajustement lors d'un calage d'historique. 30 3059443 6 Le procédé selon l'invention De façon générale, l'invention concerne un procédé pour exploiter une formation souterraine comportant des hydrocarbures, ladite formation comportant au moins un 5 faciès ayant subi au moins une phase diagénétique au cours des temps géologiques, à partir d'au moins un échantillon de roche issu de ladite formation, de données de production mesurées relatives à ladite formation et d'une représentation maillée représentative de ladite formation, ladite représentation maillée étant formée de mailles, chacune desdites mailles de ladite représentation maillée comportant au moins un 10 indicateur relatif à la répartition desdits faciès présents dans ladite maille. Le procédé selon l'invention comprend au moins les étapes suivantes : A) à partir au moins de mesures réalisées sur ledit échantillon, pour au moins un desdits faciès et pour chacune desdites mailles de ladite représentation maillée, on détermine une valeur d'au moins un paramètre de l'état diagénétique dudit 15 faciès ; B) à partir desdites données de production mesurées, pour au moins une maille de ladite représentation maillée, on met à jour au moins lesdites valeurs desdits paramètres de l'état diagénétique d'au moins un desdits faciès selon au moins les étapes suivantes : 20 i. pour chacune desdites mailles de ladite représentation maillée, à partir desdites valeurs desdits paramètres déterminés de l'état diagénétique desdits faciès, on détermine au moins une propriété pétrophysique relative à l'écoulement dans ladite formation ; ii. au moyen d'un simulateur d'écoulement et de ladite représentation maillée 25 remplie en propriétés pétrophysiques relatives à l'écoulement, on simule des données de production relatives à ladite formation ; iii. on détermine une fonction objectif mesurant au moins une différence entre lesdites données de production mesurées et lesdites données de production simulées ; 30 iv. pour au moins une desdites mailles de ladite représentation maillée, on modifie, à partir de ladite fonction objectif, au moins ladite valeur dudit paramètre de l'état diagénétique pour un desdits faciès ; et on réitère les étapes i) à iv) jusqu'à ce qu'un critère d'arrêt prédéterminé soit atteint ; 3059443 7 C- à partir d'au moins ladite représentation maillée mise à jour, on détermine un schéma d'exploitation desdits hydrocarbures de ladite formation, et on exploite lesdits hydrocarbures de ladite formation. 5 Selon un mode de mise en oeuvre de l'invention, ledit paramètre représentatif de l'état diagénétique d'un desdits faciès peut correspondre à la variabilité spatiale d'au moins un degré diagénétique d'au moins une desdites phases diagénétiques subies par ledit faciès. 10 Selon variante de l'invention, ladite variabilité spatiale peut être paramétrée par la portée d'un variogramme portant sur la distribution dudit degré diagénétique de ladite phase diagénétique. Selon un mode de mise en oeuvre de l'invention, ledit paramètre représentatif de 15 l'état diagénétique d'un desdits faciès peut correspondre à la proportion d'au moins un degré diagénétique d'au moins une desdites phases diagénétiques subies par ledit faciès. Selon une variante de l'invention, ledit paramètre peut correspondre à la proportion d'un desdits degrés diagénétiques, les proportions des autres degrés diagénétiques 20 restant constantes relativement. Selon un mode de mise en oeuvre de l'invention, à l'étape B, on peut mettre en outre à jour un paramètre représentatif de la distribution desdits faciès dans ladite formation et/ou une desdites propriétés pétrophysiques relatives à l'écoulement dans 25 ladite formation. Selon une réalisation de l'invention, on peut réaliser l'étape B) ii) au moyen d'une simulation bi-plurigaussienne. 30 Selon une variante de mise en oeuvre de l'invention, ledit critère d'arrêt peut correspondre à un minimum de ladite fonction objectif, à un nombre maximal d'itérations de l'étape B), ou bien à une différence maximale entre lesdites données de production mesurées et lesdites données de production simulées. 3059443 8 Avantageusement, ladite détermination dudit schéma d'exploitation peut comprendre la détermination de la position et de la géométrie de puits injecteurs et/ou producteurs, et/ou la détermination du type de récupération assistée. 5 En outre, l'invention concerne un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une quelconque des variantes de l'invention, lorsque ledit programme est exécuté sur un ordinateur. 10 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. 15 Présentation succincte des figures la Figure 1 a présente la distribution des faciès pour une section à z constant au sein 20 d'une représentation maillée d'une formation souterraine. la Figure 1 b présente la distribution des degrés diagénétiques par faciès pour la même section à z constant que celle présentée en Figure 1 a. La Figure 2 présente la distribution des degrés diagénétiques par faciès pour la même section à z constant que celle présentée en Figure 1 a, obtenue après une mise à jour 25 d'un paramètre représentatif de la distribution des degrés diagénétiques selon une première mise en oeuvre de l'invention. La Figure 3 présente la distribution des phases diagénétiques par faciès pour la même section à z constant que celle présentée en Figure 1 a, obtenue après une mise à jour d'un paramètre représentatif de la distribution des degrés diagénétiques selon une 30 seconde mise en oeuvre de l'invention. La Figure 4 présente la distribution des phases diagénétiques par faciès pour la même section à z constant que celle présentée en Figure 1 a, obtenue après une mise à jour d'un paramètre représentatif de la distribution des degrés diagénétiques selon une troisième mise en oeuvre de l'invention. 3059443 9 La Figure 5 présente l'évolution du « water out » mesuré au cours du temps, simulé à la première itération d'une mise en oeuvre du procédé selon l'invention, et simulé à l'itération finale de cette mise en oeuvre du procédé selon l'invention. La Figure 6 présente l'évolution, au cours des itérations du procédé selon l'invention, 5 de la fonction objectif et du paramètre représentatif de l'état diagénétique à ajuster dans le cas de la mise en oeuvre du procédé selon l'invention correspondant à celle de la Figure 5. 10 Description détaillée du procédé De façon générale, l'un des objets de l'invention concerne un procédé d'exploitation d'une formation souterraine comportant des hydrocarbures. Ce procédé comprend une étape de calage d'historique dont au moins l'un des paramètres d'ajustement est l'état 15 diagénétique d'au moins un faciès lithologique rencontré au sein de la formation souterraine. Les définitions suivantes sont utilisées au cours de la description de l'invention : faciès : on appelle faciès, ou encore faciès lithologique une classe de roches 20 présentant les mêmes caractéristiques. 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...). Ce type d'information est le plus souvent représenté par une propriété catégorielle, c'est-à-dire par un ensemble fini de 25 valeurs symboliques, représentées par des catégories (par exemple les catégories CO, C1, C2, etc.). Il s'agit d'une propriété particulièrement sensible pour l'exploitation d'un réservoir. En particulier, les proportions de faciès lithologiques peuvent avoir un impact majeur sur le comportement dynamique d'un champ pétrolier. Elles peuvent être identiques sur l'ensemble du domaine (cas stationnaire) ou varier selon la 30 position (cas non stationnaire). paramètre représentatif de l'état diagénétique d'un faciès : il s'agit d'un paramètre qui rend compte des processus physico-chimiques subis par un faciès. Notamment, il peut s'agir d'un paramètre relatifs aux différentes phases de diagénèse subies par un faciès donné au sein d'une formation souterraine, tel qu'un un paramètre décrivant la 35 variabilité spatiale des phases diagénétiques subies par un faciès donné, ou la 3059443 10 variabilité des proportions des phases diagénétiques observées au sein d'un même faciès. Il peut aussi s'agir d'un paramètre décrivant la variabilité spatiale des degrés diagénétiques de chacune des phases diagénétiques subies par un faciès donné, ou encore des proportions des degrés diagénétiques de chacune des phases 5 diagénétiques observées au sein d'un même faciès. Le procédé selon l'invention requiert de disposer : d'au moins un échantillon de roche de la formation souterraine étudiée : cet 10 échantillon peut avoir été prélevé in situ, par carottage par exemple ; de données de production mesurées relatives à la formation souterraine étudiée : il s'agit de données mesurées concrètement sur un gisement, relatives à la production de ce gisement. Il peut s'agir du volume cumulé d'huile produit, du débit d'huile, du « water cut » (rapport du volume d'eau produit sur volume total de liquide produit), de 15 la pression, du rapport du volume de gaz sur le volume d'huile, etc. Il peut aussi s'agir d'un essai de puits. Ces données de production sont mesurées au niveau de puits forés dans la formation souterraine étudiée. une représentation maillée de ladite formation : il s'agit d'une maquette de la formation souterraine étudiée, représentée généralement sur ordinateur sous la forme d'un 20 maillage ou grille, chacune des mailles de cette grille comportant une ou plusieurs valeurs de propriétés relatives à la formation étudiée. Selon l'invention, la représentation maillée comporte en chaque maille au moins une information sur la répartition des faciès présents dans la maille considérée. La répartition des faciès dans une maille peut être le pourcentage de chacun des faciès présents dans la maille 25 considérée. Un faciès peut être représenté sous la forme d'une propriété catégorielle, c'est-à-dire un ensemble fini de valeurs symboliques, représentées par des catégories (par exemple les catégories CO, C1, C2, etc.). Selon une mise en oeuvre, on utilise des données aux puits (données de diagraphie) pour déterminer une ou plusieurs courbes de proportions verticales de faciès. Ces courbes donnent la probabilité 30 d'occurrence de chaque faciès en fonction de la profondeur. Ces courbes sont ensuite utilisées pour construire, à partir de techniques de krigeage notamment, une matrice de proportion de faciès, qui détermine pour chaque maille de la grille représentative du réservoir les probabilités d'occurrence des différents faciès. Selon une mise en oeuvre de l'invention, la représentation maillée correspond au modèle géologique, 35 construit dans le but de décrire aussi précisément que possible la structure et les propriétés (dans le cas présent, au moins une répartition des différents faciès dans 3059443 11 chaque maille) de la formation étudiée. Un modèle géologique se doit de vérifier autant que possible les propriétés relatives à la formation étudiée collectées sur le terrain : les données de diagraphie mesurées le long des puits, les mesures réalisées sur des échantillons de roche prélevés dans les puits, les données déduites de 5 campagnes d'acquisition sismique, ou encore des données de production comme les débits d'huile, d'eau, les variations de pression etc. Selon une mise en oeuvre de l'invention, la représentation maillée comporte également des propriétés nécessaires à une simulation d'écoulement, telles que des valeurs de porosité, de perméabilité, de saturation en chacune des mailles de la représentation maillée, valeurs qui seront 10 mises à jour par le procédé selon l'invention. Le spécialiste en exploration et exploitation pétrolière a pleine connaissance de techniques pour construire une telle représentation maillée d'une formation. Le modèle géologique a généralement des mailles fines et par voie de conséquence un nombre de mailles important. Le calage d'historique nécessitant de réaliser plusieurs simulations d'écoulement, il peut s'avérer 15 nécessaire de déterminer une représentation maillée aux mailles plus grossières, afin de réduire le temps de calcul des simulations d'écoulement. Cette étape, bien connue dans le domaine de l'exploration et exploitation pétrolière consiste à réaliser un « upscaling » (« mise à l'échelle » en français), qui permet de passer d'une représentation à mailles fines à une représentation à mailles grossières, aussi appelée 20 modèle de réservoir. De façon générale, cette technique consiste à calculer un modèle grossier qui soit équivalent au modèle fin du point de vue des propriétés dynamiques (c'est-à-dire des propriétés liées à l'écoulement des fluides, telles que la perméabilité). un simulateur d'écoulement : un simulateur d'écoulement est un logiciel exécuté sur 25 ordinateur permettant de réaliser des simulations d'écoulement de fluides au sein d'un gisement. Ces simulations sont réalisées en résolvant un système d'équations aux dérivées partielles représentatives des écoulements, en utilisant par exemple des méthodes de volumes finis appliquées à une représentation maillée. De façon générale, lors d'une expérience de simulation de réservoir, l'ingénieur réservoir définit 30 une durée sur laquelle il souhaite simuler les écoulements dans la formation d'intérêt. Puis ces équations sont résolues par un schéma numérique, pas de temps après pas de temps, la taille du pas de temps (pouvant être une fonction ou une constante) étant intrinsèquement liée au schéma numérique choisi, et ce afin de garantir la stabilité de ce schéma numérique. Le résultat d'une simulation de réservoir est en particulier une 35 succession d'instantanés (comme des images prises à différents instants) représentatifs de l'état d'écoulement des fluides dans la formation étudiée. Par exemple, le logiciel PumaFlow (IFP Énergies nouvelles, France) est un simulateur d'écoulement. 3059443 12 Le procédé selon l'invention comporte au moins les étapes suivantes : 1. Détermination d'une représentation maillée renseignée sur l'état 5 diagénétique d'un faciès 1.1. Mesures expérimentales 1.2. Caractérisation de l'état diagénétique 1.3. Remplissage de la représentation maillée 10 2. Calage d'historique de production 2.1. Détermination de propriétés pétrophysiques 2.2. Simulation de données de production 2.3. Processus itératif d'optimisation 15 3. Exploitation des hydrocarbures de la formation Les différentes étapes du procédé selon l'invention sont détaillées ci-après. 20 1. Détermination d'une représentation maillée renseignée sur l'état diagénétique d'un faciès Au cours de cette étape, il s'agit, à partir au moins de mesures réalisées sur au moins un échantillon de roche issu de ladite formation, de déterminer une valeur d'au moins un paramètre représentatif de l'état diagénétique d'au moins un faciès donné présent dans 25 une maille donnée, et ce pour chacune des mailles de la représentation maillée représentative de la formation souterraine étudiée. Avantageusement, cette étape est réalisée au moyen de plusieurs échantillons de roche issus de la formation souterraine, prélevés en différentes localisations spatiales de la formation, de manière à caractériser de la manière la plus précise possible l'état 30 diagénétique des faciès rencontrés dans la formation. 3059443 13 1.1. Mesures expérimentales Au cours de cette sous-étape, il s'agit de caractériser l'état diagénétique d'au moins un faciès présent dans au moins un échantillon de roche issu de la formation souterraine étudiée. 5 Selon une mise en oeuvre de l'invention, cette caractérisation peut s'effectuer par une analyse de lames minces extraites d'une carotte. Cette analyse, effectuée sous un microscope conventionnel et/ou sous microscope électronique, donne accès à de nombreuses informations pétrographiques concernant l'échantillon de roche étudié (texture, type de grains, type de phases diagénétiques, porosité, etc.). On peut de cette 10 manière identifier par exemple les différentes phases diagénétiques subies par un faciès, et éventuellement de qualifier le degré diagénétique de chacune de ces phases diagénétiques (on peut par exemple qualifier si le faciès considéré a subi un degré faible, moyen ou fort de la phase diagénétique considérée). Selon une variante de réalisation de l'invention, cette caractérisation de l'état diagénétique d'un échantillon de roche à partir de 15 lame minces peut être réalisée au moyen d'une analyse d'image (par exemple au moyen d'un logiciel de traitement d'image par comptage de points). Avantageusement, on peut répéter cette sous-étape pour l'ensemble des échantillons de roche disponibles pour la formation souterraine étudiée. La répétition de ces étapes à l'ensemble des échantillons donne accès à l'évolution verticale et latérale qualitative et/ou 20 quantitative des informations désirées (par exemple le type et les proportions des grains, des phases diagénétiques, du réseau poreux, etc.). Avantageusement, on peut répéter cette sous-étape pour l'ensemble des faciès présents dans le ou les échantillons de roche considérés. 25 1.2. Caractérisation de l'état diagénétique A partir de l'analyse expérimentale décrite précédemment, on peut en déduire des informations sur l'état diagénétique du faciès étudié, comme par exemple sur la nature et sur la proportion (qualitative et/ou quantitative) des différentes phases diagénétiques subies par le faciès. 30 Selon l'invention, on caractérise cet état diagénétique par au moins un paramètre. Ce paramètre peut concerner : La variabilité spatiale des degrés diagénétiques d'au moins une phase diagénétique subie par un faciès donné : cette variabilité spatiale peut être évaluée au 35 moyen de techniques géostatistiques, bien connues du spécialiste, à partir des 3059443 14 informations sur l'état diagénétique obtenues au cours de la sous-étape précédente, par exemple à partir d'une l'analyse de lames minces. On peut par exemple construire un variogramme de la distribution spatiale des degrés diagénétiques pour chacune des phases diagénétiques d'un faciès donné, et paramétrer l'état diagénétique de cette phase 5 par la portée de ce variogramme (cf Le Ravalec et al., 1999). Cette variabilité spatiale peut aussi être représentée par des nombres aléatoires permettant de générer une réalisation d'une fonction aléatoire, ces nombres contrôlant la distribution spatiale des hétérogénéités. Selon une mise en oeuvre de l'invention, on peut utiliser pour ce faire la méthode de déformation graduelle (Hu, 2000), la méthode de perturbation des 10 probabilités (Caers et Hoffman, 2006) ou la méthode de déformation par co-simulation (Le Ravalec-Dupin et Da Veiga, 2011) pour modifier les nombres aléatoires et obtenir une nouvelle réalisation de la fonction aléatoire considérée, c'est à dire une réalisation montrant une distribution différente des degrés des phases diagénétiques dans l'espace. Les proportions des degrés diagénétiques d'au moins une phase 15 diagénétique subie par un faciès donné : selon une mise en oeuvre de l'invention, une paramétrisation basée sur les proportions des degrés diagénétiques d'une phase diagénétique observée au sein d'un même faciès est basée sur l'intensité de cette phase diagénétique. Un exemple d'une telle paramétrisation est donné dans le Tableau 1. Dans cet exemple, trois faciès sont considérés (faciès A, B et C) et sont affectés par une phase 20 diagénétique quelconque dont la proportion dans chacun des faciès est qualifiée de « faible » (notée L), « modérée » (notée M) et « forte » (notée H). Dans le cas de cet exemple, neuf paramètres sont à mettre à jour par le calage historique. Toutefois, dans la réalité, les faciès sédimentaires et les phases diagénétiques subies par un faciès donné sont plus nombreux. De plus, les proportions varient de façon significative dans l'espace. 25 On aurait alors à traiter un nombre considérable de paramètres, ce qui rendrait la résolution du problème difficile. Selon une mise en oeuvre de l'invention destinée à réduire le nombre de paramètres à déterminer, on peut piloter les variations des proportions de degré diagénétique d'au moins une phase diagénétique à partir d'un unique coefficient « a », 30 défini pour chacune des phases diagénétiques considérées, et variant entre 0 et 1. La valeur de ce coefficient, dit degré diagénétique de contrôle, correspond à la proportion d'un des degrés diagénétiques d'une phase diagénétique donnée, les proportions des autres degrés diagénétiques de la phase considérée étant constantes relativement entre elles. Par exemple, dans le cas de trois proportions P1, P2 et P3 (correspondant 35 respectivement à une proportion « faible », « moyenne », « importante ») d'une phase diagénétique donnée d'un faciès donné, on peut définir une telle paramétrisation de la manière suivante : 3059443 15 pP1 new - a pP2new = (1+( pP1init-a)/(pP2init - - - 11nP2 init pP3new = (1+( pPlinit -a)/(pP2 - )11/P':i init - -init où p est la notation utilisée pour indiquer une proportion, finit indique les valeurs initiales 5 des proportions de degré diagénétique et new celles après ajustement. Ainsi, dans cet exemple de paramétrisation, le degré diagénétique de contrôle est le degré P1, la proportion du degré de contrôle est donnée par le paramètre a, et les proportions des autres degrés diagénétiques en découlent directement. Ainsi, selon cette paramétrisation, le coefficient a devient le seul paramètre à ajuster : en le modifiant, on modifie l'ensemble 10 des proportions des degrés diagénétiques pour la phase diagénétique considérée du faciès considéré. Ainsi, si on augmente la proportion de degré diagénétique P1 de 10`)/0, on réduit d'autant les proportions des degrés diagénétiques P2 et P3. Toutefois, on préserve le rapport entre les proportions des degrés diagénétiques P2 et P3. Selon une variante de l'invention, on peut mettre en oeuvre cette paramétrisation pour chacune des 15 N phases diagénétiques d'un faciès donné. Selon une autre mise en oeuvre de l'invention, on utilise non pas un unique degré de contrôle par phase diagénétique, mais un groupe de degrés diagénétiques (par exemple on forme le groupe de contrôle L+M dans l'exemple précédemment décrit). Selon une autre alternative, on fixe comme constante la proportion d'un degré diagénétique donné ou d'un groupe de degrés diagénétiques, et on 20 modifie uniquement les proportions des autres degrés diagénétiques. Par la suite, cette paramétrisation sera appelée « paramétrisation par degré(s) de contrôle ». Tableau 1 L M H Faciès A 33,3% 16,7% 50% Faciès B 50% 50% 0 Faciès C 100% 0 0 25 A l'issue de cette sous-étape, on dispose d'au moins une valeur d'au moins un paramètre reflétant l'état diagénétique d'au moins un faciès, et ce pour au moins un échantillon de roche prélevé au sein de la formation. 1.3. Remplissage de la représentation maillée 30 Au cours de cette sous-étape, il s'agit de compléter la représentation maillée (qui contient à ce stade au moins les informations relatives aux proportions de faciès) avec les 3059443 16 valeurs d'au moins un paramètre reflétant l'état diagénétique d'au moins un des faciès, valeurs déterminées au moyen de la mise en oeuvre des deux sous-étapes précédentes. Du fait du nombre restreint d'échantillons de roche, les informations concernant le paramètre représentatif de l'état diagénétique d'un faciès peuvent être limitées dans 5 l'espace. Selon un mode de réalisation de l'invention, le remplissage des mailles de la représentation maillée peut nécessiter de recourir à des simulations de type géostatistique. Avantageusement, le remplissage des mailles tient compte de la continuité et de la connectivité des structures géologiques étudiées. Selon une mise en oeuvre de l'invention, on utilise la technique décrite dans le 10 document (Renard et al., 2008) pour générer une représentation maillée rendant compte de la distribution des faciès et des phases diagénétiques par faciès. Cette technique consiste en une simulation bi-plurigaussienne, qui considère deux modèles plurigaussiens, un pour les faciès et un pour la diagenèse. Doligez et al. (2011) ou Barbier et al. (2012) présentent des cas d'application et appliquent des techniques de simulation 15 plurigaussienne, bi-plurigaussienne, ou encore emboîtées. Ainsi, à l'issue de cette étape, on dispose d'une représentation maillée pour laquelle, dans chaque maille, est renseignée au moins une valeur d'au moins un paramètre représentatif de l'état diagénétique d'au moins un faciès présent dans la 20 formation souterraine étudiée. 2. Calage d'historique de production Au cours de cette étape, il s'agit, à partir de données de production mesurées, de 25 mettre à jour les valeurs présentes dans les mailles de la représentation maillée relatives à au moins un paramètre représentatif de l'état diagénétique d'au moins un faciès. 2.1. Détermination de propriétés pétrophysiques Au cours de cette sous-étape, il s'agit dans un premier temps de convertir les 30 valeurs relatives aux paramètres représentatifs de l'état diagénétique en propriétés pétrophysiques, telles que la porosité et la perméabilité. De cette manière, la représentation maillée comporte des propriétés relatives aux écoulements dans la formation. 3059443 17 Le spécialiste a parfaite connaissance de techniques pour une telle conversion. On pourra se référer par exemple aux documents (Doligez et al., 2011 ; Barbier et al., 2012). 2.2. Simulation de données de production 5 Au cours de cette sous-étape, il s'agit, au moyen d'un simulateur d'écoulement et de la représentation maillée remplie en propriétés pétrophysiques à l'issue de la sous-étape précédente, de simuler des données de production relatives à ladite formation. Cette étape est résolue au moyen d'un simulateur d'écoulement, tel que le logiciel PumaFlow ® (IFP Énergies nouvelles, France). 10 Très préférentiellement et afin d'avoir des données de production mesurées et simulées pouvant être comparées, la simulation de données de production est effectuée aux emplacements dans la représentation maillée correspondant aux emplacements réels des puits traversant la formation étudiée et pour lesquels on dispose de données de production mesurées. 15 Très préférentiellement et afin d'avoir des données de production mesurées et simulées pouvant être comparées, on simule à l'aide du simulateur d'écoulement le même type de données de production que le type de données de production mesurées. Par type de donnée de production on entend par exemple le volume cumulé d'huile produit, le débit d'huile, le « water cut » (rapport du volume d'eau produit sur volume total de liquide 20 produit), la pression, le rapport du volume de gaz sur le volume d'huile, etc. Il peut aussi d'agir d'un essai de puits. 2.3. Processus itératif d'optimisation 25 Au cours de cette sous-étape, il s'agit de déterminer une fonction objectif mesurant au moins la différence entre les données de production mesurées et les données de production simulées au cours de la sous-étape précédente, puis de déterminer de manière itérative, les valeurs des paramètres représentatifs de l'état diagénétique permettant de satisfaire un certain critère d'arrêt prédéfini par le 30 spécialiste. Selon un mode de mise en oeuvre de l'invention, on estime une fonction objectif en mesurant la différence entre les données de production mesurées et les données de production simulées au sens de la norme L2. Selon une réalisation de l'invention, la fonction objectif s'écrit selon une formule de la forme : 35 3059443 18 où F(2)- Nwell N time (a Donbs _ ( ) )2 1 1( 1 0- ( ) me a) 19'1 N () taire i=1 n=1 Nwell est le nombre de puits P°bs est l'ensemble des données de production mesurées 5 P5'm est l'ensemble des données de production simulées, op sont les poids associés aux données de production, o-p(i) est l'intervalle de confiance sur les données de production au puits i, Mme représente le nombre de mesures effectuées dans le temps pour le puits i, 10 2 = (2v, 112, IlAi) sont les M paramètres représentatifs de l'état diagénétique d'un faciès. Avantageusement, la fonction objectif peut comporter d'autres termes, notamment un terme visant à minimiser un écart entre des données sismiques simulées et des 15 données sismiques mesurées et/ou un terme dit de régularisation permettant de réduire l'espace des solutions possibles. Selon une mise en oeuvre de l'invention, on met à jour, en plus des paramètres représentatifs de l'état diagénétique d'au moins un faciès, un paramètre représentatif de la distribution de faciès au sein de la formation souterraine étudiée. Un tel paramètre peut 20 concerner la distribution spatiale des faciès (cf par exemple le document Le Ravalec et al., 1999), ou encore les proportions de faciès (cf par exemple le document EP 1978379 Al (US 7970593 B2)). Selon une autre mise en oeuvre de l'invention, on met à jour, en plus des paramètres représentatifs de l'état diagénétique d'au moins un faciès, les propriétés 25 pétrophysiques relatives à l'écoulement au sein de la formation souterraine étudiée. Un tel paramètre peut concerner la distribution spatiale des valeurs de porosités ou perméabilités dans les faciès. On pourra se référer au document (Le Ravalec et al., 2014) pour un tel calage. Différentes techniques d'optimisation d'une fonction objectif existent. On connaît par 30 exemple la méthode basée sur les gradients, telle que décrite par exemple dans (Gosselin et al., 2000). De façon générale, la fonction objectif relative au problème étudié étant non linéaire, ces techniques d'optimisation résolvent le problème posé de manière itérative. A partir de l'estimation de la fonction objectif à une itération donnée, c'est-à-dire pour un jeu de paramètres à une itération donnée, on estime un jeu de paramètres 35 permettant d'abaisser la valeur de la fonction objectif à l'itération suivante, et ce, en tenant 3059443 19 compte de caractéristiques de la fonction objectif (comme par exemple, le gradient de la fonction objectif dans le cas de la méthode des gradients précitée). A chaque itération, une estimation de la fonction objectif est effectuée. Notamment, à chaque itération, des données de production sont simulées sur la représentation maillée de l'itération en cours. 5 Ce processus itératif est arrêté lorsqu'un critère d'arrêt prédéfini est atteint. Ce critère peut être basé sur la valeur de la fonction objectif (par exemple lorsque l'estimation de la fonction coût a atteint un minimum, ou bien est en deçà d'un seuil prédéfini par le spécialiste), sur la mesure de la différence entre données de production simulées et données de production mesurées (par exemple lorsque la différence quadratique 10 moyenne est en deçà d'un seuil prédéfini par le spécialiste), ou encore lorsqu'un nombre maximum d'itérations du processus itératif, défini par le spécialiste, a été atteint. Le spécialiste a une parfaite connaissance de la manière de prédéfinir de tels critères d'arrêt de processus d'optimisation itératif. 15 La mise à jour selon l'invention des valeurs des paramètres représentatifs de l'état diagénétique d'un faciès donné par calage de données d'historique peut se faire pour une seule maille de la représentation maillée, sur l'ensemble de la représentation maillée (soit pour toutes les mailles de la représentation maillée) ou bien sur certaines zones identifiées (un ou plusieurs ensembles de mailles voisines) au préalable par l'ingénieur 20 réservoir. Selon une mise en oeuvre de l'invention, la mise à jour selon l'invention des valeurs des paramètres représentatifs de l'état diagénétique d'un faciès donné par calage de données d'historique de production est opérée pour l'ensemble des faciès lithologiques d'une représentation souterraine. Avantageusement, la mise à jour selon l'invention des 25 valeurs des paramètres représentatifs de l'état diagénétique d'un faciès donné par calage de données d'historique est opéré pour les faciès de la formation étudiée ayant subi des phases diagénétiques impactant le plus les propriétés pétrophysiques relatives à l'écoulement. 30 3. Exploitation des hydrocarbures de la formation Au cours de cette étape, il s'agit dans un premier temps de définir au moins un schéma d'exploitation des hydrocarbures contenus dans la formation, en tenant compte de critères technico-économiques prédéfinis par le spécialiste. Il peut s'agir d'un scénario 35 offrant un fort taux de récupération du fluide, sur une longue durée d'exploitation, et nécessitant un nombre de puits limité. Puis, selon l'invention, on exploite le fluide de la formation étudiée en fonction de ce schéma d'exploitation optimal. 3059443 20 Selon l'invention, la détermination du schéma d'exploitation est réalisée à l'aide d'une simulation d'écoulement exploitant au moins la représentation maillée de la formation étudiée mise à jour à l'issue de l'étape de calage d'historique selon l'invention. En effet, à l'issue de l'étape 2 décrite précédemment, on dispose d'une représentation 5 maillée satisfaisant les données de production mesurées. Un exemple de simulateur d'écoulement est le logiciel PumaFlow (IFP Energies nouvelles, France). Selon l'invention, à tout instant t de la simulation, le simulateur d'écoulement résout l'ensemble des équations d'écoulement propres à chaque maille et délivre des valeurs solutions des inconnues (saturations, pressions, concentrations, température,...) prédites 10 à cet instant t. De cette résolution, découle la connaissance des quantités d'huile produites et de l'état du gisement (distribution des pressions, saturations, etc.) à l'instant considéré. A partir d'un scénario d'exploitation donné, de la représentation maillée du gisement, on est alors capable de simuler la production d'hydrocarbures escomptée, au moyen du simulateur d'écoulement. 15 La définition d'un schéma d'exploitation du fluide de la formation étudiée peut consister à choisir un type de procédé de récupération (par exemple un procédé de récupération par injection d'eau) puis à déterminer, pour ce type de procédé, le nombre, la géométrie et l'implantation (position et espacement) des puits injecteurs et producteurs. 20 Afin de définir un schéma d'exploitation optimal, différents essais de différents scenarii de production peuvent être réalisés à l'aide d'un simulateur d'écoulement. Le schéma d'exploitation offrant le meilleur taux de récupération de fluide pour un coût le plus faible sera par exemple préféré. 25 En sélectionnant divers scénarii, caractérisés par exemple par diverses implantations respectives des puits injecteurs et producteurs, et en simulant la production de fluide pour chacun d'eux, on peut sélectionner le scénario permettant d'optimiser la production de la formation considérée suivant des critères technico-économiques prédéfinis par le spécialiste. 30 Les spécialistes exploitent alors le fluide de la formation considérée selon le scénario permettant d'optimiser la production du gisement, par exemple en forant les puits injecteurs et producteurs définis par ledit schéma d'exploitation optimal et/ou en mettant en oeuvre le procédé de récupération défini par le schéma d'exploitation optimal. 35 3059443 21 L'invention concerne en outre un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des caractéristiques précédentes, lorsque 5 ledit programme est exécuté sur un ordinateur. Exemple de réalisation 10 Les caractéristiques et avantages du procédé selon l'invention apparaîtront plus clairement à la lecture de l'exemple d'application ci-après. Le gisement étudié pour cette illustration correspond à une formation souterraine 15 caractérisée par 3 faciès (notés A, B et C), chacun de ces faciès étant affecté par une (unique) phase diagénétique dont le degré de transformation ou degré diagénétique varie de faible (noté « L » pour low), moyen (noté « M » pour medium) à important (noté « H pour high). Les proportions de ces degrés diagénétiques ont été déterminées par l'analyse des lames minces réalisées sur un échantillon de roche provenant de la 20 formation étudiée (cf sous-étape étape 1.1 décrite ci-dessus) et sont reportées dans le Tableau 1 présenté précédemment. La représentation maillée représentative de cette formation est composée de 200x200x100 mailles, chaque maille comprenant au moins un indicateur relatif aux faciès rencontrés (faciès notés A, B, ou C) dans la maille considérée et les proportions relatives aux degrés diagénétiques de la phase diagénétique 25 considérée (telles que présentées dans le Tableau 1), et ce pour chacun des faciès rencontrés dans la maille. La distribution des faciès et la distribution des proportions de degré diagénétique pour chaque faciès pour une section à profondeur constante donnée sont présentées respectivement en Figure 1 a et en Figure 1 b (les contours superposés sur la carte de distribution des degrés diagénétiques rappellent les limites des faciès A, 30 B et C présentés en Figure 1 a). Selon l'invention, on définit une paramétrisation représentative de l'état diagénétique de chacun des faciès rencontrés dans la formation (cf sous-étape 1.2 décrite précédemment). A titre d'illustration, on présente ci-après trois modes de mise en oeuvre 35 de l'invention, suivant trois paramétrisations différentes : 3059443 22 a) Modèle de variabilité spatiale et portée du variogramme Pour cette première mise en oeuvre, on établit un variogramme décrivant la variabilité spatiale des proportions de degré diagénétique par faciès, et on en détermine la portée. A titre illustratif, on modifie la portée de ce variogramme en la divisant artificiellement par 2. 5 La distribution des proportions de degré diagénétique par faciès en résultant est présentée en Figure 2 et est à comparer avec la distribution en Figure 1 b. Pour le faciès C, aucune modification n'est visible puisque ce faciès ne contient que le degré diagénétique faible (valeur à 100% du degré L dans le Tableau 1). Pour les 2 autres faciès, on peut observer que les hétérogénéités des degrés diagénétiques sont bien plus 10 petites. b) Variation de la distribution spatiale des hétérogénéités par déformation graduelle Pour cette deuxième mise en oeuvre, on établit un variogramme décrivant la variabilité spatiale des degrés diagénétiques subis par un faciès donné, et on génère des nombres 15 aléatoires pour produire une réalisation de la fonction aléatoire définie par le variogramme. A titre illustratif, on utilise la méthode de déformation graduelle pour modifier les nombres aléatoires. Le modèle des degrés diagénétiques ainsi obtenu est présenté en Figure 3 et est à comparer à la Figure 1 b. Le paramètre de déformation graduel appliqué est de 0.1, ce qui implique une faible perturbation des nombres 20 aléatoires et donc une faible perturbation de la position des hétérogénéités. Cette perturbation suffit néanmoins pour faire apparaître une modification de la distribution des hétérogénéités des degrés diagénétiques dans les faciès A et B. Pour le faciès C, aucune modification n'est visible puisque ce faciès ne contient que le degré diagénétique faible (valeur à 100% du degré L dans le Tableau 1). 25 c) Variation des proportions des degrés diagénétiques Pour cette troisième mise en oeuvre, on utilise la paramétrisation particulière présentée en section 1.2 concernant la représentation des proportions des degrés diagénétiques observés au sein d'un même faciès, permettant un contrôle par un paramètre a attaché à 30 au moins un degré diagénétique de contrôle. Pour cet exemple, on définit le degré diagénétique de contrôle comme étant le degré L dans le faciès A, et on décrit sa proportion par le paramètre a. A titre illustratif, on augmente la proportion du degré diagénétique faible (L) dans le faciès A, en passant d'une proportion de 1/3 à 1/2 pour ce degré diagénétique. La Figure 4 présente la distribution des degrés diagénétiques en 35 résultant. On observe bien une augmentation de la proportion du degré L, les proportions des différents degrés dans les autres faciès restant inchangées relativement. 3059443 23 On applique ensuite le procédé de calage d'historique selon l'invention, à partir de données de production mesurées sur la formation souterraine telle que décrite précédemment Pour cette mise en oeuvre du procédé selon l'invention, on met à jour uniquement 5 les valeurs des paramètres représentatifs de l'état diagénétique du faciès A, les valeurs des paramètres représentatifs de l'état diagénétique des autres faciès ayant été obtenues par ailleurs et étant fixées pour la mise en oeuvre de l'invention. Au cours de cette mise en oeuvre du procédé selon l'invention, on paramétrise l'état diagénétique du faciès A par les proportions des degrés diagénétiques subis par ce 10 faciès, en utilisant la paramétrisation dite par degré diagénétique de contrôle (telle que présentée à la sous-étape 1.2). Pour cette mise en oeuvre, le paramètre a correspond à la proportion du degré diagénétique L dans le faciès A. La formation souterraine étudiée est exploitée à partir de 2 puits : un puits injecteur 15 d'eau dans le coin inférieur gauche et un puits producteur dans le coin supérieur droit. Les débits d'injection et de production sont de 200 m3/jour. Les valeurs de water cut (ratio de production d'eau par à rapport au volume total de liquide) mesurées MES au puits producteur à différentes dates sont présentées par les losanges de la Figure 5. Le temps de percée, c'est-à-dire le moment où l'eau perce au puits producteur, est de 2250 jours. 20 Par ailleurs, les mesures effectuées sur la formation (par carottage notamment) permettent d'établir des propriétés pétrophysiques associées aux différents couples (faciès, degré diagénétique) telles que reportées dans le Tableau 2. Tableau 2 Faciès Degré diagénétique Porosité (%) Perméabilité (mD) A L 0,30 500 A M 0,15 100 A H 0,05 1 B L 0,20 200 B M 0,10 50 C L 0,05 2 25 Afin de débuter le processus itératif de mise à jour selon l'invention, on définit une valeur d'un paramètre a contrôlant le degré diagénétique L égal à 0.38. Autrement dit, on démarre le processus itératif de calage d'historique selon l'invention en considérant une 30 distribution des degrés diagénétiques caractérisée par une proportion du degré L (iow) dans le faciès A de 38%. 3059443 24 Au moyen d'un simulateur d'écoulement, on simule l'évolution au cours du temps du water cut pour cette valeur du paramètre a. Le water cut simulé WCS sur la représentation maillée comportant la valeur de 0,38 pour le paramètre a est présenté en Figure 5 par la courbe comportant des carrés, alors que le water cut mesuré MES est 5 présenté par les points symbolisés par des losanges sur cette même figure. On peut observer que ces courbes étant relativement éloignées l'une de l'autre, une mise à jour significative du paramètre a est à réaliser par le procédé selon l'invention. Par ailleurs, on peut observer que le temps de percée simulé, de 1400 jours, est nettement plus faible que le temps de percée mesuré, de 2250 jours, sur le gisement étudié. 10 Pour cette mise en oeuvre de l'invention, on considère une fonction objectif composée uniquement d'une mesure de la différence entre le water cut mesuré et le water cut simulé, ainsi qu'un critère d'arrêt basé sur la minimisation de cette fonction objectif. On modifie alors de manière itérative le paramètre a représentatif de l'état 15 diagénétique du faciès A, au moyen de cette fonction objectif évaluée à chaque itération, un water cut simulé étant évalué à chacune des itérations du processus itératif, sur la représentation maillée de l'itération courante. Pour cet exemple, la fonction objectif FO atteint un minimum au bout de 10 itérations (cf courbe avec des carrés en Figure 6). Le paramètre a dit de degré de contrôle 20 obtenu à l'itération vaut alors 0,8038 (cf courbe avec des losanges en Figure 6). L'évolution du water cut dans le temps, simulé pour la représentation maillée déterminée à l'itération 10 est présenté en Figure 5 par la courbe pleine notée SIM. On peut observer que cette courbe se superpose parfaitement aux valeurs discrètes du water cut mesurées MES (losanges en Figure 5). 25 Ainsi, le procédé selon l'invention permet de prendre en compte l'état diagénétique d'un faciès comme paramètre d'ajustement de données de production mesurées sur une 30 formation souterraine comportant des hydrocarbures. Un tel procédé permet de prendre en compte de façon plus exhaustive la complexité des phénomènes géologiques subis par une formation souterraine. En effet, la présente invention permet de prendre en compte, entre autres, les processus diagénétiques subis par des sédiments dès leur dépôt et tout au long de leur enfouissement, ces processus affectant significativement les 35 propriétés d'écoulement de la roche résultant de ce dépôt sédimentaire. La présente invention permet donc d'obtenir une représentation maillée d'une formation souterraine remplie en propriétés d'écoulement plus fiables et/ou réalistes, et permet de fait d'améliorer l'exploitation des hydrocarbures de la formation. REVENDICATIONS1. Procédé pour exploiter une formation souterraine comportant des hydrocarbures, ladite formation comportant au moins un faciès ayant subi au moins une phase diagénétique au cours des temps géologiques, à partir d'au moins un échantillon de roche issu de ladite formation, de données de production mesurées relatives à ladite formation et d'une représentation maillée représentative de ladite formation, ladite représentation maillée étant formée de mailles, chacune desdites mailles de ladite représentation maillée comportant au moins un indicateur relatif à la répartition desdits faciès présents dans ladite maille, caractérisé en ce qu'on réalise au moins les étapes suivantes : A. à partir au moins de mesures réalisées sur ledit échantillon, pour au moins un desdits faciès et pour chacune desdites mailles de ladite représentation maillée, on détermine une valeur d'au moins un paramètre de l'état diagénétique dudit faciès ; B. à partir desdites données de production mesurées, pour au moins une maille de ladite représentation maillée, on met à jour au moins lesdites valeurs desdits paramètres de l'état diagénétique d'au moins un desdits faciès selon les étapes suivantes : i. pour chacune desdites mailles de ladite représentation maillée, à partir desdites valeurs desdits paramètres déterminés de l'état diagénétique desdits faciès, on détermine au moins une propriété pétrophysique relative à l'écoulement dans ladite formation ; ii. au moyen d'un simulateur d'écoulement et de ladite représentation maillée remplie en propriétés pétrophysiques relatives à l'écoulement, on simule des données de production relatives à ladite formation ; iii. on détermine une fonction objectif mesurant au moins une différence entre lesdites données de production mesurées et lesdites données de production simulées ; iv. pour au moins une desdites mailles de ladite représentation maillée, on modifie, à partir de ladite fonction objectif, au moins ladite valeur dudit paramètre de l'état diagénétique pour un desdits faciès ; et on réitère les étapes i) à iv) jusqu'à ce qu'un critère d'arrêt prédéterminé soit atteint ; 3059443 26 C- à partir d'au moins ladite représentation maillée mise à jour, on détermine un schéma d'exploitation desdits hydrocarbures de ladite formation, et on exploite lesdits hydrocarbures de ladite formation. 5 2. Procédé selon la revendication 1, dans lequel ledit paramètre représentatif de l'état diagénétique d'un desdits faciès correspond à la variabilité spatiale d'au moins un degré diagénétique d'au moins une desdites phases diagénétiques subies par ledit faciès. 10 3. Procédé selon la revendication 2, dans lequel ladite variabilité spatiale est paramétrée par la portée d'un variogramme portant sur la distribution dudit degré diagénétique de ladite phase diagénétique. 4. Procédé selon la revendication 1, dans lequel ledit paramètre représentatif de 15 l'état diagénétique d'un desdits faciès correspond à la proportion d'au moins un degré diagénétique d'au moins une desdites phases diagénétiques subies par ledit faciès. 5. Procédé selon la revendication 4, dans lequel ledit paramètre correspond à la 20 proportion d'un desdits degrés diagénétiques, les proportions des autres degrés diagénétiques restant constantes relativement. 6. Procédé selon l'une des revendications précédentes, dans lequel, à l'étape B), on met en outre à jour un paramètre représentatif de la distribution desdits faciès 25 dans ladite formation et/ou une desdites propriétés pétrophysiques relatives à l'écoulement dans ladite formation. 7. Procédé selon l'une des revendications précédentes, dans lequel on réalise l'étape B) ii) au moyen d'une simulation bi-plurigaussienne. 30 8. Procédé selon l'une des revendications précédentes, dans lequel ledit critère d'arrêt correspond à un minimum de ladite fonction objectif, à un nombre maximal d'itérations de l'étape B), ou bien à une différence maximale entre lesdites données de production mesurées et lesdites données de production simulées. 3059443 27 9. Procédé selon l'une des revendications précédentes, dans lequel ladite détermination dudit schéma d'exploitation comprend la détermination de la position et de la géométrie de puits injecteurs et/ou producteurs, et/ou la détermination du type de récupération assistée. 5 10. Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des revendications 10 précédentes, lorsque ledit programme est exécuté sur un ordinateur.