FR 3037353 A1 20161216 FR 1555389 A 20150612 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 pétrolière consiste à rechercher des gisements d'hydrocarbures au sein d'un bassin sédimentaire. La compréhension des principes de la genèse des hydrocarbures et leurs liens avec l'histoire géologique du sous-sol, a permis la mise au point de méthodes d'évaluation du potentiel pétrolier d'un bassin sédimentaire. La démarche générale de l'évaluation du potentiel pétrolier d'un bassin sédimentaire comporte des allers-retours entre une prédiction du potentiel pétrolier du bassin sédimentaire, réalisée à partir d'informations disponibles concernant le bassin étudié (affleurements, campagnes sismiques, forages par exemple), et des forages d'exploration dans les différentes zones présentant le meilleur potentiel, afin de confirmer ou infirmer le potentiel prédit préalablement, et d'acquérir de nouvelles données pour alimenter de nouvelles études plus précises. L'exploitation pétrolière d'un gisement consiste, à partir des informations récoltées lors de la phase d'exploration pétrolière, à sélectionner les zones du gisement présentant le meilleur potentiel pétrolier, à définir des schémas d'exploitation optimaux pour ces zones (par exemple à l'aide d'une simulation de réservoir, afin de définir les nombre et positions des puits d'exploitation permettant une récupération d'hydrocarbures optimale), à forer des puits d'exploitation et, de façon générale, à mettre en place les infrastructures de production nécessaires au développement du gisement. Dans certains bassins sédimentaires, ayant subi une histoire géologique complexe, faisant interagir de nombreux processus physiques, ou lorsque le volume de données est très important, l'évaluation du potentiel pétrolier d'un bassin sédimentaire requiert, généralement, de disposer d'outils informatiques (logiciels exécutés sur ordinateur) permettant la synthèse des données disponibles, et d'outils informatiques permettant la simulation de l'histoire géologique et des multiples processus physiques qui la contrôlent. Il s'agit d'une démarche dite de "modélisation de bassin". La famille des logiciels dits de modélisation de bassin permettent de simuler en une, deux ou trois dimensions, les processus sédimentaires, tectoniques, thermiques, hydrodynamiques et de chimie organique et inorganique qui interviennent lors de la formation d'un bassin pétrolier. Concernant plus particulièrement les processus sédimentaires, les spécialistes utilisent des outils mettant en oeuvre un ensemble d'équations simulant l'évolution 3037353 2 sédimentaire d'un bassin au cours des temps géologiques, c'est-à-dire depuis le dépôt des sédiments jusqu'à un temps actuel. La simulation de l'histoire sédimentaire d'un bassin requiert la prise en compte de différents paramètres : (1) l'évaluation de l'espace disponible pour la sédimentation, lié à des mouvements tectoniques et/ou 5 eustatiques, (2) l'apport en sédiments dans le bassin, soit par les frontières, soit par l'intermédiaire d'une production ou d'une précipitation in situ, (3) le transport de ces sédiments dans l'espace disponible créé, (4) l'évolution de ces sédiments au cours de l'enfouissement, ce qu'on appelle la diagenèse. Ce type de simulation, dite simulation stratigraphique, permet notamment au spécialiste de tester différentes hypothèses 10 sur les processus sédimentaires ayant affecté un bassin, et de mettre à jour ces hypothèses par comparaison du résultat obtenu par simulation avec l'état actuel observé des dépôts sédimentaires d'un bassin. Le logiciel DionisosFlowe (IFP Energies nouvelles, France) est un exemple d'un logiciel, dit simulateur stratigraphique, mettant en oeuvre une simulation stratigraphique. 15 La diagénèse figure donc parmi les processus sédimentaires majeurs de l'histoire d'un bassin sédimentaire. La diagenèse consiste en la 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. Ces sédiments vont être soumis, 20 lors de leur enfouissement progressif dans le bassin, à des conditions de pression et de température qui vont les transformer. Cette transformation survient généralement à faible profondeur et se déroule selon différentes étapes qui varient selon la nature des sédiments et des conditions d'enfouissement. En comparaison des roches sédimentaires clastiques de type grès ou argiles, la 25 diagénèse des roches carbonatées est généralement complexe, et peut notamment consister en de nombreux processus chimiques et/ou biologiques qui s'enchainent les uns après les autres. La diagenèse augmentant avec le temps et la profondeur, elle est marquée par (1) la compaction des sédiments avec perte d'eau (tassement mécanique lié au poids des terrains déposés par-dessus le sédiment ; ce processus 30 tend à réduire la porosité dans la roche et à augmenter les points de contact entre les grains), (2) une augmentation de la température par enfouissement, ce qui favorise les réactions chimiques et (3) une multiplication de réactions variées et complexes telles que : la transformation (ou épigénisation) de certains minéraux en d'autres minéraux (par exemple dolomitisation), la dissolution des grains à leurs points de 35 contact, et la précipitation (cimentation) dans les espaces inter-granulaires. Par ailleurs, chaque roche carbonatée de chaque bassin subit des étapes diagénétiques propres, l'intensité de chaque étape pouvant même varier d'un point à un autre du 3037353 3 bassin considéré. On parle alors de chemin diagénétique, qui peut être plus ou moins complexe. État de la technique 5 Les documents suivants seront cités au cours de la description : Adelinet, M., Fortin, J., & Guéguen, Y., 2011a. Dispersion of elastic moduli in a porous-cracked rock: Theoretical predictions for squirt-flow. Tectonophysics, 503(1), 10 173-181. Adelinet, M., Dorbath, C., Le Ravalec, M., Fortin, J., & Guéguen, Y., 2011b. Deriving microstructure and fluid state within the lcelandic crust from the inversion of tomography data. Geophysical Research Letters, 38(3). Granjeon, D. & Joseph, P., 1999. Concepts and applications of a 3-D multiple lithology, diffusive model in stratigraphic modeling. Numerical Experiments in Stratigraphy Recent Advances in Stratigraphic and Sedimentologic Computer Simulations SEPM Special Publications No 62. Xu, S., & Payne, M. A. (2009). Modeling elastic properties in carbonate rocks. The Leading Edge, 28(1), 66-74. Les processus intervenant au cours de la diagénèse d'une roche carbonatée 25 ont pour effet de modifier les paramètres microstructuraux d'une roche (nature et géométrie des grains formant la matrice de la roche ; nature et géométrie des pores de la roche). De fait, les propriétés mécaniques de la roche carbonatée se trouvent affectées par les processus chimiques et/ou biologiques intervenant au cours de la diagénèse. Ainsi, les transformations diagénétiques subies par une roche au cours du 30 temps ont pour conséquence la variation au cours des temps géologiques des propriétés mécaniques des roches (modules élastiques), et a fortiori de leurs propriétés pétrophysiques (porosité, perméabilité par exemple). La prise en compte de la diagénèse dans les outils informatiques de simulation 35 stratigraphique est à ce jour restreinte, puisque seul l'impact de la compaction sédimentaire sur les paramètres mécaniques est simulé numériquement. Ainsi, dans 15 20 3037353 4 le document (Granjeon & Joseph, 1999), on décrit des lois de compaction classique qui relient la porosité du sédiment à l'enfouissement, permettant ainsi de quantifier le volume des couches sédimentaires. Si une telle restriction peut être satisfaisante (c'est-à-dire produisant un résultat de simulation suffisamment proche de la réalité) 5 pour des roches sédimentaires de type clastique, elle ne peut conduire à une simulation satisfaisante de la diagénèse des roches carbonatées. Or les roches carbonatées représentent plus de 50 (3/0 des roches réservoirs actuellement exploitées dans le monde. Il apparaît donc important de pouvoir prendre correctement en compte le phénomène de diagénèse, dans sa complexité, dans le 10 cas de bassins sédimentaires comportant des roches carbonatées. Notamment, il apparaît important de prendre en compte, dans une simulation stratigraphique, l'aspect évolutif dans le temps, induit par la diagénèse, des paramètres mécaniques d'une roche carbonatée. 15 On connaît le document (Xu and Paine, 2009) qui présente un procédé pour déterminer des propriétés mécaniques d'une roche carbonatée à partir de mesures expérimentales. Plus précisément, les propriétés mécaniques sont déterminées à partir d'un modèle de microstructure, en testant différentes valeurs des paramètres de celui-ci (aplatissement et augmentation de la porosité). Ces tests ne considèrent 20 toutefois pas une évolution dans le temps des paramètres du modèle microstructural, et donc des paramètres mécaniques des roches carbonatées. On connaît également le document (Adelinet et al, 2011a) qui concerne un procédé pour déterminer des propriétés structurales d'une roche basaltique à partir de mesures réalisées sur le terrain et d'une représentation par milieu effectif. La 25 modélisation en milieu effectif permet, à partir d'une description fine de la microstructure d'une roche à l'échelle d'un Volume Elémentaire Représentatif (VER), de calculer des propriétés mécaniques homogénéisées à l'échelle de ce volume. Dans ce document, des données de tomographie sismique sont utilisées pour inverser deux paramètres microstructuraux du modèle effectif : la densité de fissure 30 et le module d'incompressibilité du fluide remplissant les inclusions de porosité. Ce document ne considère pas une évolution dans le temps des paramètres d'un modèle microstructural, et donc des paramètres mécaniques de la roche considérée. 35 3037353 5 La présente invention a pour objet un procédé permettant de déterminer une évolution des paramètres mécaniques d'une roche carbonatée au cours des différentes étapes de la diagénèse subies par cette roche au sein d'un bassin sédimentaire. Ces paramètres sont ensuite pris en compte dans une simulation 5 stratigraphique, afin de contribuer à une meilleure appréhension d'un bassin sédimentaire comportant des roches carbonatées, et ainsi à une évaluation pétrolière plus fiable de ce type de bassin. 10 Le procédé selon l'invention Ainsi, la présente invention concerne un procédé d'exploitation pétrolière d'un bassin sédimentaire, ledit bassin comportant au moins une couche de sédiments carbonatés. Au moyen d'un simulateur stratigraphique permettant de reconstituer 15 l'histoire sédimentaire dudit bassin depuis un temps géologique t jusqu'à un temps actuel, au moyen d'au moins un échantillon de roche de ladite couche et d'une échelle d'un Volume Elémentaire Représentatif, ladite échelle étant déterminée en fonction dudit échantillon, le procédé comporte les étapes suivantes pour ladite couche : 20 A. à partir de mesures réalisées sur ledit échantillon, on détermine des paramètres d'un modèle microstructural représentatif de l'état diagénétique de ladite couche audit temps actuel, lesdits paramètres dudit modèle microstructural étant définis à ladite échelle ; B. à partir de mesures réalisées sur ledit échantillon, on identifie au moins une 25 étape diagénétique subie par lesdits sédiments de ladite couche depuis ledit temps géologique t jusqu'audit temps actuel, et on détermine des variations minimales et maximales desdits paramètres dudit modèle microstructural pour chacune desdites étapes diagénétiques ; C. on détermine au moins un paramètre mécanique desdits sédiments de ladite 30 couche pour chacune desdites étapes diagénétiques, au moyen d'une modélisation en milieu effectif et desdites variations desdits paramètres dudit modèle microstructural déterminées pour chacune desdites étapes diagénétiques ; ainsi que les étapes suivantes : 3037353 6 D. on évalue le potentiel pétrolier dudit bassin au moins au moyen dudit simulateur et desdits paramètres mécaniques déterminés pour chacune desdites étapes diagénétiques, et on sélectionne au moins une zone dudit bassin présentant ledit potentiel pétrolier le plus élevé ; 5 E. on exploite ledit bassin en fonction de ladite zone sélectionnée. Préférentiellement, lesdites mesures peuvent consister en des mesures de caractérisation de ladite roche réalisées au microscope, par diffractométrie de Rayons X, ou par porosimétrie. 10 Selon un mode de réalisation de l'invention, lesdits paramètres microstructuraux peuvent comprendre la souplesse des interfaces entre grains de ladite roche. 15 Selon un mode de réalisation de l'invention, ladite souplesse peut être supposée invariante au cours desdites étapes diagénétiques. Selon un mode de réalisation de l'invention, on peut déterminer au moins un desdits paramètres dudit modèle microstructural par modélisation inverse. 20 Préférentiellement, lesdites variations minimales et maximales peuvent être déterminées à partir de mesures réalisées sur un nombre d'échantillons de ladite roche au moins égal au nombre desdites étapes diagénétiques. 25 Avantageusement, on peut déterminer lesdites variations minimales et maximales de la micro-porosité, de la macro-porosité, et de la composition minéralogique. Selon un mode de réalisation de l'invention, à partir des paramètres 30 mécaniques déterminés pour chacune desdites étapes diagénétiques, on peut déterminer la perméabilité de ladite couche pour chacune desdites étapes diagénétiques. 3037353 7 Selon un mode de réalisation de l'invention, à partir des paramètres mécaniques déterminés pour chacune desdites étapes diagénétiques, on peut déterminer un cube de données sismiques synthétiques pour chacune desdites étapes diagénétiques. 5 Préférentiellement, à l'étape D, on peut simuler en outre au moins un processus choisi parmi les processus tectoniques, thermiques, hydrodynamiques, et de chimie organique et inorganique ayant affecté ledit bassin. 10 Avantageusement, à l'étape E, on peut réaliser au moins un forage d'exploitation et/ou d'exploration dans lesdites zones sélectionnées pour la récupération des hydrocarbures présents au sein dudit bassin. En outre, l'invention concerne un produit programme d'ordinateur 15 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 la description ci-dessus, lorsque ledit programme est exécuté sur un ordinateur. 20 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. 25 Présentation succincte des figures - La Figure 1 présente un exemple illustratif d'un chemin diagénétique constitué de quatre étapes diagénétiques différentes. - la Figure 2 présente l'évolution au cours de la diagénèse des paramètres 30 microstructuraux associés à l'exemple présenté en Figure 1. - la Figure 3 présente l'évolution des modules élastiques et des vitesses élastiques en fonction des étapes diagénétiques déterminées pour l'exemple présenté en Figure 1. 3037353 8 - La Figure 4 présente l'évolution de la perméabilité au cours des étapes diagénétiques déterminées pour l'exemple présenté en Figure 1. 5 Description détaillée du procédé Les définitions suivantes sont utilisées au cours de la description de l'invention : Modélisation en milieu effectif : il s'agit d'une modélisation physique permettant d'estimer les propriétés effectives d'un milieu à partir des 10 propriétés locales de ses constituants. Volume Elémentaire Représentatif d'un milieu: il s'agit d'un volume de taille suffisamment grande pour pouvoir définir des propriétés homogènes et représentatives du milieu étudié. Inclusions carbonatées : il s'agit d'éléments solides composants une roche 15 carbonatée. Il peut s'agit de bioclastes (éléments fossiles d'origine animale ou végétale, le plus souvent en fragments) ou d'oolithes (sphères formées d'un nucleus et de différentes enveloppes). Etude de lames minces : à partir d'un échantillon de roche amincie jusqu'à le rendre transparent, on observe la microstructure de la roche en lumière 20 transmise au microscope optique. Modèle microstructural : il s'agit de la simplification de la microstructure, vue en microscopie par exemple, afin de la traduire en modèle de milieu effectif. La description d'un modèle microstructural comprend au moins la caractérisation de la matrice minérale (nature et géométrie des grains 25 formant la matrice) et la caractérisation des inclusions de porosité (nature et géométrie des pores). Par nature de la matrice minérale, on entend la nature lithologique des grains (par exemple quartz, argiles, calcaire). Par géométrie des grains formant la matrice, on entend la forme des grains de la matrice (ellipsoïdes plus ou moins aplaties). Par nature des pores, on 30 entend distinguer la porosité de fissure, la porosité issue d'une dissolution, etc. Et par géométrie des pores, on entend la forme des pores (pièce de monnaie pour les fissures par exemple, ou ellipsoïdes pour les pores équants). 3037353 9 Paramètres ou propriétés mécaniques : il s'agit des modules élastiques définis en mécanique des milieux continus. Le comportement élastique d'un matériau homogène isotrope et linéaire est caractérisé par deux modules élastiques indépendants (modules d'incompressibilité et de cisaillement) 5 qui sont des constantes intrinsèques du matériau. L'invention concerne un procédé d'exploitation pétrolière d'un bassin sédimentaire comportant au moins une couche de sédiments carbonatés. En particulier, l'invention concerne la modélisation, au sein d'une simulation 10 stratigraphique, du phénomène de diagénèse subi par des sédiments carbonatés. Une étape importante du procédé selon l'invention est la caractérisation mécanique du phénomène de diagénèse subi par les sédiments carbonatés du bassin étudié. Par procédé d'exploitation pétrolière d'un bassin sédimentaire, on entend un procédé permettant l'exploitation des hydrocarbures présents au sein dudit bassin 15 sédimentaire. La présente invention requiert : de disposer d'un simulateur stratigraphique selon l'art antérieur : un simulateur stratigraphique est un logiciel visant à reconstituer les 20 processus sédimentaires ayant affecté le bassin depuis un temps géologique t jusqu'au temps actuel. La simulation de l'histoire sédimentaire d'un bassin requiert l'élaboration de systèmes d'équations permettant d'estimer : (1) l'évaluation de l'espace disponible pour la sédimentation, lié à des mouvements tectoniques et/ou eustatiques, (2) 25 l'apport en sédiments dans le bassin, soit par les frontières, soit par l'intermédiaire d'une production ou d'une précipitation in situ, (3) le transport de ces sédiments dans l'espace disponible créé, (4) l'évolution de ces sédiments au cours de l'enfouissement, c'est-à-dire la diagenèse. Par simulateur stratigraphique selon l'art antérieur, on entend un 30 simulateur stratigraphique modélisant la diagénèse par le seul phénomène de compaction des sédiments ; de disposer d'au moins un échantillon de roche de chacune des couches de sédiments carbonatés du bassin sédimentaire étudié : cet échantillon peut avoir été prélevé in situ, par carottage par exemple ; 3037353 10 de définir une échelle d'un Volume Elémentaire Représentatif (VER) : l'échelle d'un VER est fonction de la taille des échantillons de roche mis à disposition. Il s'agit de s'affranchir de tout élément microstructural perturbant la représentativité du volume (grande fissure 5 parcourant tout l'échantillon, trous non associés à une porosité dans tout l'échantillon, etc...). 10 1. 15 2. 20 3. La présente invention comporte au moins les étapes suivantes : Caractérisation mécanique de la diagénèse 1.1. Détermination des paramètres d'un modèle microstructural de l'état diagénétique actuel 1.2. Identification des différentes étapes diagénétiques 1.3. Détermination des variations minimales et maximales des paramètres du modèle microstructural pour chacune des étapes diagénétiques 1.4. Détermination des paramètres mécaniques par modélisation en milieu effectif pour chacune des étapes diagénétiques Evaluation du potentiel pétrolier Exploitation du bassin sédimentaire Les principales étapes de la présente invention sont détaillées ci-après. Elles sont illustrées sur un exemple (non limitatif) de chemin diagénétique subi par une 25 couche composée de sédiments carbonatée donnée. 1. Caractérisation mécanique de la diagénèse L'objet de cette première étape consiste en la caractérisation mécanique de la 30 diagenèse ayant affecté les couches de sédiments carbonatés du bassin sédimentaire étudié. Cette étape peut être découpée, non limitativement, en quatre sous-étapes, applicables parallèlement ou séquentiellement à chacune des couches de sédiments carbonatés du bassin étudié. Les quatre sous-étapes en question sont détaillées pour une couche de sédiments carbonatés donnée. 3037353 11 1.1. Détermination des paramètres d'un modèle microstructural de l'état diagénétique actuel 5 Au cours de cette sous-étape, il s'agit de déterminer les paramètres d'un modèle microstructural représentatif de l'état diagénétique actuel de la couche de sédiments carbonatés considérée, à partir de mesures expérimentales réalisées sur au moins un échantillon de roche de la couche considérée. Selon l'invention, les paramètres du modèle microstructural sont définis à l'échelle d'un Volume 10 Elémentaire Représentatif (VER), de sorte à pouvoir exploiter ces paramètres par une approche en milieux effectifs dans la sous-étape 1.4 décrite ultérieurement. Un échantillon de roche, prélevé par exemple par un carottage, permet d'accéder aux paramètres du modèle microstructural de l'état diagénétique final de la roche étudiée. En effet, certaines étapes du processus de formation de la roche 15 prélevée au temps actuel ne sont plus visibles que sous la forme de traces ou d'éléments géométriques (phase minéralogique incluse ou entourant une autre par exemple). Ainsi, au cours de cette étape, on réalise des mesures sur un échantillon prélevé au temps actuel afin de déterminer les paramètres microstructuraux 20 représentatifs de la matrice minérale (nature et géométrie des grains formant la matrice) et les paramètres microstructuraux représentatifs des inclusions de porosité (nature et géométrie des pores). Selon l'invention, on a recours à des mesures directes sur échantillon par au moins une des techniques suivantes : 25 une étude microscopique, réalisée par exemple à l'aide d'un microscope optique ou électronique à balayage : une étude microscopique permet de caractériser la matrice de la roche étudiée ainsi que sa porosité. Ainsi, par une étude microscopique, on peut accéder à la géométrie et la disposition des phases solides (la matrice), et la géométrie et la disposition de la porosité de 30 la roche. On peut alors renseigner la matrice, les inclusions cristallines et les supports de porosité (pores sphériques ou fissures par exemple) à l'échelle d'un Volume Elémentaire Représentatif (VER) ; une diffractométrie des Rayons X (DRX), réalisée à l'aide d'un diffractomètre : une diffractométrie des Rayons X permet de quantifier les différentes phases 35 minéralogiques d'un échantillon donné, ce qui permet de renseigner une fraction volumique des inclusions solides à l'échelle d'un Volume Elémentaire Représentatif (VER) ; 3037353 12 une porosimétrie, réalisée à l'aide d'un spectromètre à résonance magnétique nucléaire (RMN), d'un porosimétre au Mercure ou d'un porosimétre à l'Hélium : ce type de mesures permet de quantifier le ratio entre la micro-porosité et la macro-porosité à l'échelle d'un Volume Elémentaire 5 Représentatif (VER). Selon un mode de mise en oeuvre de l'invention, afin de compléter le modèle microstructural, on peut avoir recours à une modélisation inverse réalisée à partir de mesures sur au moins un échantillon de la roche carbonatée étudiée. 10 En effet, les roches carbonatées sont bien souvent caractérisées par un agencement minéralogique hétérogène, ce qui entraîne une complexification dans la réponse mécanique de ces roches. Certains paramètres du modèle microstructural, notamment la souplesse existante entre les différentes inclusions carbonatées (bioclastes, oolithes par exemple), ne peuvent alors être directement approchés par 15 des mesures. Afin de quantifier ces paramètres, on peut alors avoir recours à une modélisation inverse. On appelle modélisation inverse une technique d'inversion itérative. Plus précisément, on construit une fonction objectif mesurant l'écart entre des données expérimentales et des données théoriques, calculées à partir de valeurs initiales pour 20 les paramètres à déterminer, puis on modifie les valeurs de ces paramètres, itération après itération, jusqu'à trouver un minimum de la fonction objectif. De nombreux algorithmes de minimisation de fonction objectif sont connus du spécialiste, tels que la méthode de Gauss-Newton, la méthode de Newton-Raphson ou encore le gradient conjugué. Selon un mode de réalisation préféré de la présente invention, on utilise la 25 méthode de Gauss-Newton. Selon un mode de mise en oeuvre de la présente invention, les données expérimentales de la fonction objectif sont des mesures ultrasoniques des vitesses des ondes sismiques de compression (ondes P) et des ondes sismiques de cisaillement (ondes S). Ces mesures peuvent avoir été obtenues en laboratoire, ou 30 bien avoir été obtenues par une campagne d'acquisition sismique suivie d'un traitement sismique, et d'une mise à l'échelle comme présentée dans la demande FR 2951555 (US 12/908130). Selon un mode de mise en oeuvre de la présente invention, le calcul des vitesses théoriques, à partir de valeurs des paramètres microstructuraux à l'itération 35 courante, peut être obtenu par une modélisation en milieu effectif telle que décrite dans (Adelinet et al, 2011b). 3037353 13 Les données théoriques sont alors comparées aux données expérimentales, et la modélisation inverse permet de minimiser l'écart entre ces deux jeux de données, en ajustant les valeurs des paramètres microstructuraux recherchés. 5 Ainsi, à l'issue de cette première sous-étape, on obtient un modèle microstructural, renseigné à l'échelle d'un Volume Elémentaire Représentatif (VER), représentatif de l'état diagénétique au temps actuel de la couche de sédiments carbonatés considérée. 10 1.2. Identification des différentes étapes diagénétiques Au cours de cette sous-étape, il s'agit d'identifier les différents stades diagénétiques par lesquels est passé la roche carbonatée du bassin sédimentaire étudié, à partir de mesures réalisées sur au moins un échantillon de la roche 15 carbonatée étudiée. Selon l'invention, on identifie au moins une étape diagénétique. Préférentiellement, on identifie plusieurs étapes diagénétiques. Si les échantillons de roche prélevés au temps actuel renseignent sur un état diagénétique final d'un point de vue des propriétés mécaniques, une étude fine de la microstructure via des lames minces renseigne sur l'histoire diagénétique subie par la 20 roche. En effet, certaines étapes précoces ont été seulement partiellement effacées par les étapes tardives et peuvent donc être encore identifiées. A partir de différentes lames minces, le géologue carbonatiste est capable d'identifier et d'ordonner les différents processus subis par la roche carbonatée au cours de la diagenèse, telles que la cimentation, la dolomitisation, l'aragonitisation ou encore la dissolution. La 25 Figure 1 présente un schéma à but illustratif du chemin diagénétique suivi par une roche carbonatée donnée. Ainsi, cette Figure présente une succession d'images, chaque image simulant une visualisation au microscope d'un échantillon de la roche considérée pour une étape diagénétique donnée. Le chemin diagénétique de la roche considérée est constitué de quatre étapes diagénétiques : une étape de cimentation 30 50 (les grains étant représentés en gris moyen), une étape de dissolution 51 (entrainant la formation de macro-porosité, représentée par des ellipses blanches), une étape de dolomitisation S2 (entrainant le remplacement des minéraux de calcite par de la dolomite), et une étape de dissolution S3 (entrainant la formation de micro-porosité, représentée par des ellipses blanches intra-grains). 35 3037353 14 1.3. Détermination des variations minimales et maximales des paramètres du modèle microstructural pour chacune des étapes diagénétiques Au cours de cette sous-étape, il s'agit de déterminer les bornes de variation des 5 paramètres du modèle microstructural déterminés à l'étape 1.1, et ce, pour chacune des étapes diagénétiques identifiées à l'étape 1.2. Pour cette étape, on part du principe qu'un échantillon de la couche carbonatée considérée, voire une partie d'un échantillon, n'a pas subi le même état d'avancement diagénétique que respectivement un autre échantillon de cette couche, 10 ou qu'une autre partie d'un même échantillon. Ainsi, les mesures de paramètres microstructuraux peuvent être différentes d'un échantillon à l'autre, voire d'une partie d'un échantillon à l'autre. Selon l'invention, les valeurs minimales et maximales des paramètres du modèle microstructural sont déterminées à partir de mesures réalisées sur au moins un échantillon de roche de la couche carbonatée. Préférentiellement, on 15 évalue les valeurs minimales et maximales des paramètres du modèle microstructural sur plusieurs échantillons, afin de bénéficier de la dispersion des mesures. Préférentiellement, les valeurs minimales et maximales des paramètres du modèle microstructural sont déterminées à partir d'un nombre d'échantillons de roche au moins égal au nombre d'étapes diagénétiques identifiées à l'étape 1.2. 20 Selon l'invention, les mesures utilisées pour la détermination des valeurs minimales et maximales des paramètres mécaniques sont réalisées par au moins une des techniques décrites dans l'étape 1.1 (c'est-à-dire une étude microscopique, une diffractométrie, une porosimétrie). En particulier, l'étude de lames minces permet une quantification des paramètres microstructuraux tels que le remplacement de 25 cristaux de calcite par des cristaux de dolomite, le remplissage incomplet d'une porosité par une phase minérale. Selon un mode de mise en oeuvre de l'invention, le spécialiste fixe au moins une borne de variations d'au moins un paramètre mécanique à une valeur prédéfinie. Par exemple, si un échantillon ou le nombre d'échantillons disponibles ne permettent 30 pas d'accéder à une borne de variations d'un des paramètres mécaniques, le spécialiste peut fixer cette borne à partir de base de données existantes dans le domaine, de ces connaissances générales, etc. Selon un mode préféré de mise en oeuvre de l'invention, on identifie au moins les valeurs minimales et maximales de la micro-porosité, de la macro-porosité et de 35 la composition minéralogique. En faisant varier les paramètres du modèle microstructural entre ces valeurs minimales et maximales au cours du temps (la durée de chaque étape diagénétique 3037353 15 pouvant être choisie arbitrairement), on dispose alors de l'évolution des paramètres du modèle microstructural au cours des différentes étapes diagénétiques du chemin diagénétique établi à l'étape 1.2 précédente. Ainsi, la Figure 2 montre l'évolution, au cours de la diagénèse, des paramètres microstructuraux associés à l'exemple 5 présenté en Figure 1. Plus précisément, la Figure 2A présente l'évolution de la (macro-)porosité d (grandeur sans unité comprise entre 0 et 1) au cours de l'étape diagénétique Si, la Figure 2B présente l'évolution du ratio R (grandeur sans unité comprise entre 0 et 1) de remplacement de la calcite par la dolomite au cours de l'étape diagénétique S2, et la Figure 20 présente l'évolution de la (micro-)porosité 10 au cours de l'étape diagénétique S3. Selon un mode de réalisation dans lequel le modèle microstructural comprend la souplesse des interfaces entre grains (voir l'étape 1.1), on suppose que cette souplesse est invariante au cours des étapes diagénétiques. 15 1.4. Détermination des paramètres mécaniques par modélisation en milieu effectif pour chacune des étapes diagénétiques Au cours de cette étape, à partir des paramètres microstructuraux définis à 20 l'étape 1.1 et de l'évolution de ces paramètres au cours des différentes étapes diagénétiques déterminée à l'étape 1.3, on détermine au moins un paramètre mécanique de la roche carbonatée étudiée par modélisation en milieu effectif, et ce pour chacune des étapes diagénétiques identifiées à l'étape 1.2. Préférentiellement, et dans le cas d'une roche isotrope, on détermine deux paramètres mécaniques : le 25 module d'incompressibilité et le module de cisaillement. La modélisation en milieu effectif permet, à partir d'une description fine de la microstructure d'une roche à l'échelle d'un Volume Elémentaire Représentatif (VER), de calculer les propriétés mécaniques homogénéisées. L'évolution des paramètres microstructuraux ayant été déterminée au cours des différentes étapes diagénétiques 30 identifiées, les propriétés mécaniques sont calculées directement par homogénéisation pour chacune des étapes du chemin diagénétique. Ce calcul repose sur la résolution du problème d'Eshelby, c'est-à-dire sur la résolution de la perturbation induite au premier ordre par la présence d'une inclusion ellipsoïdale dans une matrice. 35 Dans l'exemple présenté en Figure 1, plusieurs inclusions étant présentes dans la matrice, on parle alors de problème auxiliaire d'Eshelby. Dans ce cas, des méthodes de calcul des moyennes sont utilisées pour calculer les paramètres 3037353 16 mécaniques du milieu. Ces calculs se font dans un cas isotrope et deux modules élastiques EM indépendants sont calculés au cours du chemin diagénétique. Ainsi, la Figure 3A présente en trait plein l'évolution au cours de la diagénèse du module d'incompressibilité (« bulk modulus » en anglais), qui est une constante propre au 5 matériau étudié, reliant la contrainte au taux de déformation d'un matériau isotrope soumis à une compression isostatique. La Figure 3A présente également en trait pointillé l'évolution au cours de la diagénèse du module de cisaillement (« shear modulus » en anglais), qui est une constante intrinsèque au matériau étudié, qui intervient dans la caractérisation des déformations causées par des efforts de 10 cisaillement reliant la contrainte au taux de déformation d'un matériau isotrope soumis à une compression isostatique. Par ailleurs, à partir de l'évolution des modules élastiques au cours de la diagenèse et en résolvant l'équation de Christoffel, on peut en déduire l'évolution des vitesses V des ondes sismiques P (courbe en trait plein en Figure 3B) et S (courbe en trait pointillé en Figure 3B) au cours de la 15 diagenèse. Selon un mode de mise en oeuvre de la présente invention, on vérifie au cours de cette étape la cohérence entre le résultat de la modélisation obtenue par la présente invention et les mesures expérimentales représentatives de l'étape diagénétique terminale. Ainsi, la Figure 3B représente par des triangles les mesures des vitesses des ondes sismiques P et S réalisées en laboratoire (mesures 20 ultrasoniques par exemple) sur des échantillons de roche. Selon un autre mode de mise en oeuvre de la présente invention, à partir des propriétés mécaniques évaluées pour chaque étape du chemin diagénétique, on 25 détermine la perméabilité de la roche carbonatée pour chaque étape du chemin diagénétique. Pour ce faire, les modèles de milieu effectif, initialement renseignés en propriétés mécaniques, sont convertis en perméabilité. On a ainsi accès à l'évolution de la perméabilité k au cours des différentes étapes de la diagenèse, telle que présentée en Figure 4. Le spécialiste a parfaite connaissance de méthodes pour 30 transformer des paramètres mécaniques en perméabilité. Selon un autre mode de mise en oeuvre de la présente invention, à partir des propriétés mécaniques déterminées pour chaque étape du chemin diagénétique, on construit des cubes de données sismiques, dits synthétiques, pour chacune des étapes diagénétiques identifiées. Pour ce faire, on construit un cube d'impédance 35 sismique pour chaque étape diagénétique, à partir des propriétés mécaniques déterminées pour l'étape considérée. Puis on utilise une technique de simulation de données sismiques qui permet, à partir d'une ondelette sismique, de transformer ces 3037353 17 cubes d'impédances en données sismiques synthétiques. Le spécialiste a parfaite connaissance de méthodes pour transformer des propriétés mécaniques en impédances sismiques, et transformer des impédances sismiques en données sismiques synthétiques. On peut alors comparer le cube de données sismiques 5 synthétiques obtenu par la présente invention pour l'étape diagénétique finale avec un cube de données sismiques réelles. En fonction des conclusions de cette comparaison, le spécialiste peut déduire si certaines hypothèses faites sur les paramètres de la simulation stratigraphique sont pertinentes ou non, et, en conséquence, modifier ou non les paramètres en question. 10 Ainsi la présente invention permet de faire le lien entre la description géologique et sédimentologique des différentes étapes de diagenèse subies par une roche carbonatée et les propriétés mécaniques, et éventuellement pétrophysiques et/ou sismique, de la roche au cours de ces différentes étapes diagénétiques. 15 2. Evaluation du potentiel pétrolier A l'issue de l'étape précédente, on obtient une modélisation de l'évolution des 20 paramètres mécaniques d'une roche carbonatée au cours du temps. Selon l'invention, cette modélisation est prise en compte dans une simulation stratigraphique, permettant ainsi de contribuer à une meilleure compréhension de l'histoire sédimentaire du bassin étudié. On peut en outre avoir recours à d'autres outils de la famille de la modélisation 25 de bassin afin de simuler les processus tectoniques, thermiques, hydrodynamiques et de chimie organique et inorganique ayant affecté le bassin étudié. Un exemple d'un tel outil de modélisation de bassin est le logiciel TEMISFLOW (IFP Energies nouvelles, France). Ainsi, à l'issue de cette étape, le spécialiste peut disposer d'informations sur : 30 i. la mise en place des couches sédimentaires, ii. les effets de la diagénèse sur les sédiments ainsi déposés, iii. leur réchauffement au cours de leur enfouissement, iv. les modifications de pressions de fluides résultant de cet enfouissement, v. la formation des hydrocarbures formés par thermogénèse, 3037353 18 vi. le déplacement de ces hydrocarbures dans le bassin sous l'effet de la flottabilité, de la capillarité, des différences de gradients de pression, vii. la quantité d'hydrocarbures issus de la thermogénèse. A partir de telles informations, le spécialiste a alors connaissance des zones 5 dudit bassin comportant des hydrocarbures, ainsi que la teneur, la nature et la pression des hydrocarbures qui y sont piégés. Le spécialiste est alors en mesure de sélectionner la ou les zones du bassin étudié présentant le meilleur potentiel pétrolier. 10 3. Exploitation du bassin sédimentaire L'exploitation pétrolière du bassin peut alors prendre plusieurs formes, notamment : la réalisation de forages d'exploration dans les différentes zones sélectionnées 15 comme présentant le meilleur potentiel, afin de confirmer ou infirmer le potentiel estimé préalablement, et d'acquérir de nouvelles données pour alimenter de nouvelles études plus précises, la définition de schémas d'exploitation optimaux pour les zones sélectionnées, par exemple à l'aide d'une simulation de réservoir, afin de définir les nombre et 20 positions des puits d'exploitation permettant une récupération d'hydrocarbures optimale, la réalisation de forages d'exploitation (puits producteurs ou injecteurs) pour la récupération des hydrocarbures présents au sein du bassin sédimentaire dans les zones sélectionnées comme présentant le meilleur potentiel, 25 la mise en place des infrastructures de production nécessaires au développement du gisement. Produit programme d'ordinateur En outre, l'invention concerne un produit programme d'ordinateur 30 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é tel que décrit précédemment, lorsque ledit programme est exécuté sur un ordinateur. REVENDICATIONS1) Procédé d'exploitation pétrolière d'un bassin sédimentaire, ledit bassin comportant au moins une couche de sédiments carbonatés, au moyen d'un simulateur stratigraphique permettant de reconstituer l'histoire sédimentaire dudit bassin depuis un temps géologique t jusqu'à un temps actuel, au moyen d'au moins un échantillon de roche de ladite couche et d'une échelle d'un Volume Elémentaire Représentatif, ladite échelle étant déterminée en fonction dudit échantillon, caractérisé en ce qu'on réalise les étapes suivantes pour ladite couche : A. à partir de mesures réalisées sur ledit échantillon, on détermine des paramètres d'un modèle microstructural représentatif de l'état diagénétique de ladite couche audit temps actuel, lesdits paramètres dudit modèle microstructural étant définis à ladite échelle ; B. à partir de mesures réalisées sur ledit échantillon, on identifie au moins une étape diagénétique subie par lesdits sédiments de ladite couche depuis ledit temps géologique t jusqu'audit temps actuel, et on détermine des variations minimales et maximales desdits paramètres dudit modèle microstructural pour chacune desdites étapes diagénétiques ; C. on détermine au moins un paramètre mécanique desdits sédiments de ladite couche pour chacune desdites étapes diagénétiques, au moyen d'une modélisation en milieu effectif et desdites variations desdits paramètres dudit modèle microstructural déterminées pour chacune desdites étapes diagénétiques ; et en ce qu'on réalise les étapes suivantes : D. on évalue le potentiel pétrolier dudit bassin au moins au moyen dudit simulateur et desdits paramètres mécaniques déterminés pour chacune desdites étapes diagénétiques, et on sélectionne au moins une zone dudit bassin présentant ledit potentiel pétrolier le plus élevé ; E. on exploite ledit bassin en fonction de ladite zone sélectionnée. 3037353 20 2) Procédé selon la revendication précédente, dans lequel lesdites mesures consistent en des mesures de caractérisation de ladite roche réalisées au microscope, par diffractométrie de Rayons X, ou par porosimétrie. 5 3) Procédé selon l'une des revendications précédentes, dans lequel lesdits paramètres microstructuraux comprennent la souplesse des interfaces entre grains de ladite roche. 4) Procédé selon la revendication 3, dans lequel ladite souplesse est supposée 10 invariante au cours desdites étapes diagénétiques. 5) Procédé selon l'une des revendications précédentes, dans lequel on détermine au moins un desdits paramètres dudit modèle microstructural par modélisation inverse. 15 6) Procédé selon l'une des revendications précédentes, dans lequel lesdites variations minimales et maximales sont déterminées à partir de mesures réalisées sur un nombre d'échantillons de ladite roche au moins égal au nombre desdites étapes diagénétiques. 20 7) Procédé selon l'une des revendications précédentes, dans lequel on détermine lesdites variations minimales et maximales de la micro-porosité, de la macroporosité, et de la composition minéralogique. 25 8) Procédé selon l'une des revendications précédentes, dans lequel à partir des paramètres mécaniques déterminés pour chacune desdites étapes diagénétiques, on détermine la perméabilité de ladite couche pour chacune desdites étapes diag én étiques. 30 9) Procédé selon l'une des revendications précédentes, dans lequel à partir des paramètres mécaniques déterminés pour chacune desdites étapes diagénétiques, on détermine un cube de données sismiques synthétiques pour chacune desdites étapes diagénétiques. 3037353 21 10) Procédé selon l'une des revendications précédentes, dans lequel, à l'étape D, on simule en outre au moins un processus choisi parmi les processus tectoniques, thermiques, hydrodynamiques, et de chimie organique et inorganique ayant affecté ledit bassin. 5 11) Procédé selon l'une des revendications précédentes, dans lequel, à l'étape E, on réalise au moins un forage d'exploitation et/ou d'exploration dans lesdites zones sélectionnées pour la récupération des hydrocarbures présents au sein dudit bassin. 10 12) 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 15 précédentes, lorsque ledit programme est exécuté sur un ordinateur.