FR 3045842 A1 20170623 FR 1562666 A 20151217 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. Plus particulièrement, la présente invention concerne la modélisation de bassin, appliquée notamment à l'évaluation du potentiel pétrolier d'un bassin sédimentaire ayant subi une histoire géologique complexe. 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. Un bassin sédimentaire résulte du dépôt, au cours des temps géologiques, de sédiments au sein d'une dépression de la croûte terrestre. Ces sédiments, meubles et riches en eau vont être soumis, lors de leur enfouissement progressif dans le bassin, à des conditions de pression et de température qui vont les transformer en roches sédimentaires compactes, dites couches géologiques. L'architecture au temps actuel d'un bassin sédimentaire résulte notamment d'une déformation mécanique du sous-sol au cours des temps géologiques. Cette déformation comporte a minima une compaction des couches géologiques, due à l'enfouissement 3045842 2 progressif de ces couches au sein du bassin, sous l'effet de l'arrivée de nouveaux sédiments. Mais un bassin sédimentaire est également le plus souvent soumis à des mouvements tectoniques de grande ampleur, générant par exemple des plissements des couches géologiques, ou encore des failles venant créer une rupture dans les couches 5 géologiques. La nature des hydrocarbures présents dans un bassin sédimentaire quant à elle résulte notamment du type de matière organique présente dans les sédiments s'étant déposés, mais aussi des conditions de pression et de température auxquelles est soumis le bassin au cours des temps géologiques. 10 La Figure 1 montre une représentation schématique d'un bassin sédimentaire, comportant plusieurs couches géologiques (a,c), délimitées par des interfaces sédimentaires (b) traversées par une faille (e), et une accumulation d'hydrocarbures (d) dans une des couches géologiques du bassin considéré (c). 15 La formation d'un bassin sédimentaire implique donc un grand nombre de processus physiques et chimiques complexes, pouvant qui plus est interagir. Face à cette complexité, la prédiction du potentiel pétrolier d'un bassin sédimentaire requiert de disposer d'outils informatiques permettant de simuler, de la façon la plus réaliste possible, les phénomènes physiques et chimiques impliqués lors de la formation du bassin étudié. 20 Ce type de reconstitution de l'histoire de la formation d'un bassin sédimentaire, appelée aussi modélisation de bassin, est le plus souvent réalisée à l'aide d'une famille d'outils informatiques permettant 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. 25 La modélisation de bassin comporte classiquement trois étapes : une étape de géo-modélisation, qui consiste à construire une représentation maillée du bassin étudié. Cette représentation maillée est le plus souvent structurée en couches, c'est-à-dire qu'un groupe de mailles est affecté à chaque couche géologique du bassin modélisé. Puis, chaque maille de cette représentation maillée est remplie par une ou 30 plusieurs propriétés pétrophysiques, telles que la porosité, le faciès (argile, sable, etc) ou encore la teneur en matière organique. La construction de ce modèle se base sur des 3045842 3 données acquises au cours de campagnes sismiques, de mesures dans des puits, de carottages, etc. une étape de reconstruction structurale de l'architecture du bassin : il s'agit de reconstruire les architectures passées du bassin. Pour ce faire, on déforme la 5 représentation maillée construite à l'étape précédente afin de représenter l'évolution anti- chronologique de l'architecture du sous-sol au cours des temps géologiques, et ce pour différentes périodes temporelles appelées aussi pas de temps. une étape de simulation de bassin : il s'agit d'une simulation numérique d'une sélection de phénomènes physiques et chimiques se déroulant au cours de l'évolution du bassin et 10 contribuant à la formation des pièges pétroliers. Cette simulation de bassin s'effectue pas de temps par pas de temps, et s'appuie pour chaque pas de temps sur les représentations maillées construites pour chaque pas de temps à l'étape précédente. En particulier, une simulation de bassin fournit une cartographie prédictive du sous-sol, indiquant l'emplacement probable des gisements, ainsi que la teneur, la nature et la 15 pression des hydrocarbures qui y sont piégés. En fournissant des informations quantitatives et fiables, cette démarche intégrée de modélisation de bassin permet d'augmenter le taux de succès lors d'un forage d'un puits d'exploration. 20 De façon générale, un bassin sédimentaire peut être soumis au cours de son histoire à des contraintes mécaniques caractérisées par des composantes dans les trois dimensions de l'espace, ces contraintes pouvant être locales ou régionales, et variables au cours du temps. Ces contraintes mécaniques sont d'une part induites par les dépôts de sédiments eux-mêmes. Dans ce cas, les contraintes mécaniques comportent une composante verticale, 25 liée au poids des sédiments sur les couches déjà déposées, mais comportent le plus souvent également une composante horizontale, les dépôts sédimentaires n'étant généralement pas invariants latéralement. D'autre part, un bassin sédimentaire subit, tout au long de sa formation, des contraintes mécaniques induites par des mouvements tectoniques liés à la géodynamique terrestre, tels que des mouvements en extension (provoquant 30 l'ouverture du bassin avec par exemple la formation de rift) ou encore des mouvements en compression (provoquant des plissements, des chevauchements, des fractures au sein du bassin etc). Ces mouvements tectoniques induisent le plus souvent des variations de contraintes mécaniques dans les trois dimensions de l'espace. A noter qu'une couche déjà 3045842 4 déposée va subir des variations de contraintes induites par les mouvements tectoniques auxquels est soumis un bassin sédimentaire tout au long de sa formation. Très classiquement, comme décrit par exemple le document (Schneider F., 2003), les logiciels de simulation de bassin font une hypothèse de variations uniquement verticales des 5 contraintes mécaniques affectant un bassin sédimentaire. Plus précisément, les logiciels de simulation de bassin ne prennent en compte que la composante verticale des variations de contraintes mécaniques induites par le poids des dépôts sédimentaires successifs au cours du temps. On parle alors de simulation 1D des effets mécaniques. La non prise en compte des variations tridimensionnelles des contraintes mécaniques 10 peut avoir de lourdes conséquences dans l'évaluation du potentiel pétrolier d'un bassin sédimentaire. En effet, les contraintes mécaniques horizontales peuvent par exemple engendrer des fractures ou des plissements dans les couches géologiques du bassin, ce qui peut fortement modifier les caractéristiques des gisements pétroliers et des roches couvertures de ces gisements, et par voie de conséquence les chemins préférentiels 15 d'écoulement, les niveaux de pression dans le bassin, la localisation des pièges à hydrocarbures, etc. Une prise en compte en trois dimensions des variations de contraintes mécaniques subies par un bassin sédimentaire au cours de sa genèse apparaît donc importante pour une 20 prédiction réaliste du potentiel pétrolier dudit bassin, et ce d'autant plus lorsque ce bassin a subi une histoire géologique complexe. État de la technique 25 Les documents suivants seront cités au cours de la description : Nayroles, B., G. Touzot and P. Villon. 1991. The diffuse approximation. C. R. Acad. Sci., Paris, série II. 313: 293-296. Scheichl, R., Masson, R., Wendebourg, J., Decoupling and Block Preconditioning for 30 Sedimentary Basin Simulations, Computational Geosciences 7(4), pp. 295-318, 2003. Schneider F., Modelling multi-phase flow of petroleum at the sedimentary basin scale. Journal of Geochemical exploration 78-79 (2003) 693-696. 3045842 5 Schneider, F., S. Wolf, I. Faille, D. Pot, A 3D Basin Model for Hydrocarbon Potential Evaluation : Application to Congo Offshore , Oil & Gas Science and Technology - Rev. IFP, Vol. 55 (2000), No. 1, pp. 3-13. Steckler, M.S., and A.B. Watts, Subsidence of the Atlantic-type continental margin off 5 New York, Earth Planet. Sci. Lett., 41, 1-13, 1978. Zoback ,M. D., Reservoir Geomechanics, 2010. Zienkiewicz, O.C., RL Taylor et JZ Zhu, The Finite Element Method: Ifs Basis and Fundamentals (seventh ed.), 2013. 10 On connaît le brevet US 8271243 B2 qui divulgue une coopération entre une simulation de bassin et une simulation géomécanique, en vue d'une prise en compte des effets géomécaniques tridimensionnels au cours de la modélisation de bassin. La modélisation géomécanique, appliquée au domaine de l'évaluation du potentiel 15 pétrolier d'un bassin sédimentaire, utilise classiquement des historiques de pression, de température et de saturation, des quantités de sédiments déposés ou érodés, d'éventuelles contraintes tectoniques, ainsi que des lois de comportement associées aux différentes lithologies du domaine modélisé, et ce afin de décrire le comportement géomécanique d'un bassin au cours des temps géologiques. 20 La mise en oeuvre d'une modélisation géomécanique comporte généralement deux phases : une phase de géo-modélisation, qui consiste en la construction d'un maillage du bassin sédimentaire à étudier. Comme pour la simulation de bassin, cette représentation maillée est le plus souvent structurée en couches, c'est-à-dire qu'un groupe de mailles est 25 affecté à chaque couche géologique du bassin modélisé. Puis, chaque maille de cette représentation maillée est associée par exemple à une loi de comportement géomécanique, une porosité, une pression, une contrainte ou encore une densité. La construction de ce modèle se base sur des données acquises au cours de campagnes sismiques, de mesures dans des puits, de carottages, etc. Pour des raisons de résolution 30 numérique des équations impliquées en simulation de bassin et en simulation géomécanique, une représentation maillée adaptée à une simulation de bassin peut être différente d'une représentation maillée adaptée à une simulation géomécanique. 3045842 6 une étape de simulation numérique géomécanique, qui permet de calculer l'évolution au cours du temps de la distribution des contraintes dans chacune des mailles de la représentation maillée, ainsi que les déformations résultantes. Ce type de technique est utilisé dans le domaine de l'industrie pétrolière, mais il est aussi utilisé par exemple dans 5 le domaine de la géotechnique. Les simulateurs géomécaniques utilisent pour ce faire une représentation maillée telle que décrite précédemment, et vont résoudre de manière discrète l'équation de la conservation de la quantité de mouvement par la méthode des éléments finis. Cette technique permet notamment de prédire les variations de contraintes dans les trois dimensions de l'espace. 10 Le brevet US 8271243 B2 décrit un procédé basé sur une coopération entre une simulation de bassin et une simulation géomécanique. Dans l'approche décrite dans le brevet US 8271243 B2, on ne détermine les variations de contraintes géomécaniques tridimensionnelles auxquelles est soumise une couche donnée qu'une seule fois, uniquement au cours de son dépôt. Autrement dit, dans l'approche décrite dans le brevet US 15 8271243 B2, on néglige les effets géomécaniques tridimensionnels résultant du dépôt de couches ultérieures ou de mouvements tectoniques affectant le bassin après le dépôt de cette couche. Une telle mise en oeuvre de la prise en compte de l'évolution des contraintes géomécaniques tridimensionnelles au fur et à mesure de la mise en place d'un bassin 20 sédimentaire est donc insatisfaisante, dès que ce bassin est caractérisé par des dépôts variables latéralement (induisant des variations de contraintes horizontales, et non pas uniquement verticales) et est soumis à des contraintes tectoniques (comprenant généralement une composante horizontale non nulle). Il peut notamment en résulter une erreur dans l'estimation finale des pressions de pores et des saturations en fluide dans 25 chacune des mailles de la représentation maillée représentative du temps actuel et obtenue après la mise en oeuvre du procédé selon le brevet US 8271243 B2. Or une connaissance précise de ces informations est cruciale pour l'évaluation pétrolière d'un bassin sédimentaire (par exemple la décision d'exploiter un bassin sédimentaire donné se base notamment sur la connaissance de la saturation en hydrocarbures dans ce bassin) et son exploitation (par 30 exemple la décision d'utiliser tel ou tel schéma d'exploitation pétrolière se base sur la connaissance des pressions de pores au sein de ce bassin). 3045842 7 L'objet de la présente invention concerne la prise en compte dans une simulation de bassin des effets géomécaniques tridimensionnels auxquels est soumis une couche lors de son dépôt, mais aussi des effets géomécaniques tridimensionnels engendrés par le dépôt ou l'érosion de cette couche sur les couches sous-jacentes, déjà déposées. De plus, l'invention 5 peut permettre la prise en compte de contraintes tectoniques. Ainsi, la présente invention vise à une meilleure prédiction des quantités prédites par une simulation de bassin, notamment les pressions de pores et les saturations en fluide. 10 Le procédé selon l'invention Ainsi, la présente invention concerne un procédé d'exploitation d'un bassin sédimentaire comportant des hydrocarbures, comprenant une reconstitution de la formation dudit bassin, ladite formation comprenant au moins le dépôt de deux couches sédimentaires, ladite reconstitution étant réalisée au moyen d'un simulateur de bassin et d'un simulateur 15 géomécanique coopérant entre eux, au moyen de mesures de propriétés relatives audit bassin et de représentations maillées représentatives dudit bassin au moment du dépôt de chacune desdites couches sédimentaires. A partir desdites mesures et d'une reconstitution de ladite formation de la première desdites couches, ladite reconstitution pour au moins une desdites couches sus-jacentes à ladite première couche est mise en oeuvre selon au moins 20 les étapes suivantes : A. au moyen dudit simulateur de bassin, on réalise une simulation de bassin conjointe de ladite couche et d'au moins une couche sous-jacente à ladite couche, et on détermine un premier jeu de paramètres ; B. à partir d'au moins une partie dudit premier jeu desdits paramètres et dudit 25 simulateur géomécanique, on réalise une simulation géomécanique conjointe de ladite couche et d'au moins une couche sous-jacente à ladite couche, et on détermine un deuxième jeu de paramètres ; C. on mesure un écart entre au moins une partie desdits paramètres dudit premier jeu et au moins une partie desdits paramètres dudit deuxième jeu, et 30 on réitère les étapes A) à C) en appliquant une correction à ladite simulation de bassin si ledit écart est supérieur à un seuil prédéfini ; 3045842 8 Puis, au moyen de ladite reconstitution réalisée pour lesdites couches, on sélectionne au moins une zone dudit bassin comportant lesdits hydrocarbures, et on exploite ledit bassin en fonction de ladite zone sélectionnée. 5 Selon un mode de mise en oeuvre de l'invention, ladite reconstitution relative à ladite première couche peut être réalisée selon au moins les étapes suivantes : i. au moyen dudit simulateur de bassin, on réalise une simulation de bassin relative à ladite première couche, et on détermine un premier jeu de paramètres ; 10 ii. à partir d'au moins une partie dudit premier jeu desdits paramètres et dudit simulateur géomécanique, on réalise une simulation relative à ladite première couche, et on détermine un deuxième jeu de paramètres ; iii. on mesure un écart entre au moins une partie desdits paramètres dudit premier jeu et au moins une partie desdits paramètres dudit deuxième jeu, et 15 on réitère les étapes i) à iii) en appliquant une correction à ladite simulation de bassin si ledit écart est supérieur à un seuil prédéfini. Selon un mode de mise en oeuvre de l'invention, ledit premier jeu de paramètres peut comprendre au moins la porosité et la pression en chacune des mailles de ladite 20 représentation maillée. Selon un mode de mise en oeuvre de l'invention, ledit deuxième jeu de paramètres peut comprendre au moins la porosité en chacune des mailles de ladite représentation maillée. 25 Préférentiellement, ledit écart peut être fondé sur une mesure de la différence entre ladite porosité issue dudit premier jeu de et ladite porosité issue dudit deuxième jeu. Avantageusement, ledit écart peut être un écart absolu MESabs défini selon une 30 formule du type : 3045842 9 MESabs = MaXnENO(Pe - (PP2.1) où (/),.: (respectivement qqn) correspond à ladite porosité déterminée par ledit simulateur de bassin (respectivement par ledit simulateur géomécanique) en une maille n de ladite représentation maillée comprenant N mailles. 5 Selon un mode de mise en oeuvre de l'invention, ledit écart peut être un écart relatif MESrel défini selon une formule du type : mEsrel 2 xmaxnENG (PZ -(Pri,l) = (Pn9-E(Pnb où (/),-, (respectivement (png) correspond à ladite porosité déterminée par ledit simulateur de 10 bassin (respectivement par ledit simulateur géomécanique) en une maille n de ladite représentation maillée comprenant N mailles. Selon un mode de mise en oeuvre de l'invention, ledit procédé peut comprendre une étape de simulation de l'érosion d'au moins une desdites couches et/ou de simulation d'une 15 lacune géologique. 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 20 pour la mise en oeuvre du procédé selon la description ci-dessus, lorsque ledit programme est exécuté sur un ordinateur. 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 25 figures annexées et décrites ci-après. 3045842 10 Présentation succincte des figures La Figure 1 montre une représentation schématique du sous-sol d'un bassin pétrolier. La Figure 2 montre un exemple de bassin sédimentaire (à gauche) et un exemple de 5 représentation maillée (à droite) de ce bassin. La Figure 3 représente la reconstruction structurale, représenté par 3 états de déformations à trois périodes de temps différentes, d'un bassin sédimentaire. La Figure 4 présente une courbe représentant l'évolution de la contrainte effective en fonction de la porosité. 10 La Figure 5 compare la prédiction de l'évolution de la pression de pores au sein d'un bassin sédimentaire, obtenue par le procédé selon l'invention et obtenue par un procédé selon l'art antérieur. 15 Description détaillée du procédé L'invention concerne un procédé d'exploitation d'un bassin sédimentaire comprenant des hydrocarbures, notamment l'identification d'au moins une zone dudit bassin dans laquelle des hydrocarbures ont pu s'accumuler, en vue de l'extraction de ces hydrocarbures. 20 L'un des objets de l'invention concerne une reconstitution réaliste de la formation d'un bassin sédimentaire résultant du dépôt d'au moins deux couches sédimentaires, ladite reconstitution d'au moins la deuxième couche étant réalisée via un simulateur de bassin et un simulateur géomécanique coopérant entre eux. Par la suite, on considère que le dépôt 25 d'une couche sédimentaire correspond à une période de temps (ou pas de temps) d'une simulation. A noter qu'une période de temps d'une simulation peut toutefois également correspondre à l'érosion d'une couche sédimentaire, ou encore à une lacune géologique. Selon l'invention, à partir de la reconstitution de la formation d'une première couche (par exemple la couche s'étant déposée en premier dans le bassin), on réalise une reconstitution 30 pour une couche sus-jacente à ladite première couche, en appliquant au moins une simulation de bassin et une simulation géomécanique pour ladite couche sus-jacente de 3045842 11 façon conjointe à au moins une couche sous-jacente (c'est-à-dire une couche déjà déposée, comme par exemple la première couche) de ce bassin. La présente invention requiert de disposer : 5 de mesures de propriétés relatives au bassin : il s'agit de mesures réalisées in situ (par exemple par carottage, via des diagraphies réalisées dans des puits, par des campagnes d'acquisition sismique, etc), en différents points du bassin étudié, nécessaires à une simulation de bassin et à une simulation géomécanique, telles que la porosité, la perméabilité, ou encore la lithologie au temps actuel. Selon un mode de mise en oeuvre de 10 l'invention visant à prendre en compte les contraintes tridimensionnelles induites par des mouvements tectoniques, il convient de disposer d'estimations de la valeur des contraintes in situ et de leur orientation dans l'espace. Ces estimations peuvent être obtenues par exemple à partir de « leak-off tests » (tests d'intégrité de puits), de mini-frac (par exemple Zoback, 2010) et d'analyse de l'ovalisation des puits. Une autre source d'information peut 15 être l'analyse géologique de la géométrie actuelle du bassin, qui permet d'estimer les raccourcissements et extension tectonique éventuels, pour aboutir idéalement à un scénario cinématique pertinent. d'un simulateur de bassin, c'est-à-dire d'un logiciel permettant d'effectuer une simulation de bassin de façon numérique, à l'aide d'un ordinateur. Plus précisément, un 20 simulateur de bassin permet de simuler numériquement l'évolution (ce qui inclut leur genèse et leur migration) des fluides (hydrocarbures, mais également l'eau de formation) au sein du bassin étudié et de leurs propriétés (évolution des pressions de fluide, des saturations et des températures), ainsi que l'évolution de propriétés pétrophysiques des roches constituant les couches sédimentaires du bassin étudié (notamment la porosité et la perméabilité). Selon 25 l'invention, une simulation de bassin est réalisée à l'aide du simulateur de bassin pour une succession de périodes temporelles (appelées aussi pas de temps), chaque période correspondant à un événement géologique, tel que le dépôt ou l'érosion d'une couche sédimentaire, ou encore une lacune géologique. Selon l'invention, le simulateur de bassin requiert de disposer de représentations maillées du bassin pour chaque pas de temps de la 30 simulation (la construction de telles représentations maillées sera décrite au cours des étapes 1 et 2 décrites ci-après). Ainsi, la simulation de bassin consiste en la résolution d'un système d'équations différentielles décrivant l'évolution au cours du temps des grandeurs physiques étudiées. Pour ce faire, on peut par exemple utiliser une discrétisation par la 3045842 12 méthode des volumes finis, comme décrit par exemple dans (Scheichl et al., 2003). Conformément au principe des méthodes de volumes finis centrés sur les mailles, les inconnues sont discrétisées par une valeur constante par maille et les équations de conservation (masse ou chaleur) sont intégrées en espace sur chaque maille et en temps 5 entre deux pas de temps successifs. Les équations discrètes expriment alors que la quantité conservée dans une maille à un pas de temps donné est égale à la quantité contenue dans la maille au pas de temps précédent, augmentée des flux de quantités entrés dans la maille et réduite des flux de quantités sorties de la maille par ses faces, plus les apports extérieurs. A chaque pas de temps et dans chaque maille de la représentation maillée du bassin au pas 10 de temps considéré, le simulateur de bassin requis pour la mise en oeuvre de l'invention permet au moins de calculer les quantités physiques suivantes : la pression de pores et la porosité. Préférentiellement, on détermine également par simulation de bassin les températures et les saturations. Un exemple d'un tel simulateur de bassin est le logiciel TemisFlowTM (IFP Énergies nouvelles, France). 15 - d'un simulateur géomécanique, c'est-à-dire d'un logiciel permettant d'effectuer une simulation géomécanique de façon numérique, à l'aide d'un ordinateur. Plus précisément, un simulateur géomécanique permet de simuler l'évolution des contraintes et déformations au sein d'un bassin sédimentaire au cours de sa formation. La mise en oeuvre de la simulation géomécanique nécessite de disposer d'une représentation maillée du bassin étudié pour 20 chaque période (ou pas de temps) de l'histoire du bassin pour laquelle on souhaite faire une estimation des effets géomécaniques. Selon un mode de réalisation de la présente invention, la simulation géomécanique repose sur une discrétisation et une résolution de l'équation de la conservation de la quantité de mouvement par la méthode des éléments finis (cf par exemple Zienkiewicz et al., 2013). Dans ce cas, et si la représentation maillée pour la 25 simulation de bassin est adaptée à une résolution par la méthode des volumes finis, on peut déterminer une représentation maillée adaptée à la simulation géomécanique à partir d'une représentation maillée adaptée à la simulation de bassin par un remaillage. L'étape de remaillage peut consister en le découpage des hexaèdres dégénérés, éventuellement présents dans les maillages adaptés à une résolution par la méthode des volumes finis, en 30 des tétraèdres adaptés à une résolution par la méthode des éléments finis. Dans le cas de représentations maillées différentes entre simulation de bassin et simulation géomécanique, on peut aussi assurer la correspondance entre les deux représentations maillées via un algorithme dit de « mapping » (« cartographie » en français). Le simulateur géomécanique selon l'invention requiert au moins une valeur de pression de pores pour chaque pas de 3045842 13 temps et en chaque maille de la représentation maillée relative au pas de temps considéré, et permet de calculer au moins le tenseur des contraintes et le tenseur des déformations en chacune des mailles de la représentation maillée, à partir d'un état initial défini en termes de contraintes, pressions et porosités. En outre, chaque maille de la représentation maillée 5 associée à un pas de temps donné est associée au moins à une loi de comportement géomécanique. La loi de comportement géomécanique peut être exprimée par un module d'Young, un coefficient de Poisson, une limite d'élasticité et/ou encore une loi d'écrouissage, ces paramètres peuvent être estimés à partir du faciès (argile, sable, etc) présent dans la maille considérée. Selon un mode de mise en oeuvre de l'invention, le simulateur 10 géomécanique peut en outre prendre en compte des conditions aux limites, permettant notamment de prendre en compte des mouvements tectoniques régionaux. Un exemple d'un tel simulateur géomécanique est le logiciel ABAQUSTM (Dassault Systèmes, France). 15 La présente invention comporte alors au moins les étapes suivantes 1. Construction de représentations maillées du bassin 2. Reconstitution de la formation du bassin 2.1. Application d'une simulation de bassin 20 2.2. Application d'une simulation géomécanique 2.3. Contrôle de cohérence 3. Exploitation du bassin sédimentaire 25 Les principales étapes de la présente invention sont détaillées ci-après. 1. Construction de représentations maillées du bassin La mise en oeuvre de l'invention repose sur une simulation de bassin coopérant avec une simulation géomécanique. Par conséquent le procédé selon l'invention requiert de 3045842 14 disposer d'une représentation maillée pour chaque pas de temps de simulation considéré pour la mise en oeuvre de l'invention, un pas de temps correspondant au dépôt d'une couche sédimentaire, à une lacune géologique, ou encore à l'érosion d'une couche sédimentaire. Selon un mode de mise en oeuvre de l'invention, cette étape est réalisée tout d'abord 5 en établissant une représentation maillée du bassin étudié au temps actuel (étape 1.1), puis en reconstruisant les architectures passées du bassin (étape 1.2), en partant du temps actuel jusqu'à un temps géologique t, antérieur à l'actuel. Ce mode particulier de mise en oeuvre de l'invention est décrit plus en détail ci-après. 10 1.1 Construction d'une représentation maillée au temps actuel Une représentation maillée est une maquette du bassin sédimentaire généralement représentée sur ordinateur, sous la forme d'un maillage ou grille, chaque maille étant caractérisée par une ou plusieurs propriétés pétrophysiques relatives au bassin, telles que porosité, faciès, perméabilité etc. La construction de ce modèle se base sur des données 15 acquises au cours de campagnes sismiques, de mesures dans des puits, de carottages, etc. Plus précisément, la construction d'une représentation maillée d'un bassin consiste à discrétiser en trois dimensions l'architecture du bassin, à attribuer des propriétés à chacune des mailles de cette représentation maillée, et à ajouter des conditions aux limites de cette représentation pour rendre compte de l'interaction de la zone modélisée avec son 20 environnement. Pour ce faire, on exploite notamment les mesures de propriétés réalisées en différents points du bassin décrites ci-dessus, que l'on extrapole et/ou interpole, en les différentes mailles de la représentation maillée, suivant des hypothèses plus ou moins restrictives. Le plus souvent la discrétisation spatiale d'un bassin sédimentaire est organisée en 25 couches de mailles, représentant chacune les différentes couches géologiques du bassin étudié. La Figure 2 illustre à gauche un exemple de bassin sédimentaire, et à droite un exemple de représentation maillée de ce bassin. La représentation maillée pouvant être utilisée pour la mise en oeuvre de l'invention comprend notamment en chaque maille une information sur la lithologie, une valeur de 30 porosité, une valeur de perméabilité, ainsi que des propriétés relatives aux fluides (notamment la saturation). A noter que le spécialiste est à même de déduire, à partir de ces 3045842 15 quantités, des informations relatives à la compaction dans chacune des mailles de la représentation maillée. 5 1.2. Reconstruction structurale Au cours de cette sous-étape, il s'agit de reconstruire les architectures passées du bassin, en partant du temps actuel jusqu'à un temps géologique t, antérieur à l'actuel. Pour ce faire, on déforme la représentation maillée construite à l'étape précédente afin de représenter l'évolution anti-chronologique de l'architecture du sous-sol au cours des temps 10 géologiques, et ce pour chaque pas de temps de la simulation. On obtient ainsi une représentation maillée pour chaque pas de temps de la simulation, depuis le temps actuel jusqu'au temps géologique t. Selon un mode de réalisation de la présente invention, la reconstruction structurale peut être particulièrement simple si elle se base sur l'hypothèse que sa déformation résulte 15 uniquement d'une combinaison de mouvements verticaux par compaction du sédiment ou par surrection ou affaissement de son socle. Cette technique, connue sous le terme de « backstripping » (ou «décompaction progressive du bassin » en français) est décrite par exemple dans (Steckler et Watts, 1978). Selon un autre mode de réalisation de la présente invention, dans le cas de bassins 20 dont l'histoire tectonique est complexe, notamment dans le cas de bassins présentant des failles, il convient d'utiliser des techniques aux hypothèses moins restrictives, telles que la restauration structurale. Une telle restauration structurale est décrite par exemple dans le document FR 2 930 350 A (US 2009/0265152 A). La restauration structurale consiste à calculer les déformations successives que le bassin a subies, en intégrant les déformations 25 dues à la compaction et celles qui résultent des forces tectoniques. Dans l'exemple de la Figure 3, trois états sont utilisés pour représenter la déformation du sous-sol au cours des temps géologiques. La représentation maillée de gauche représente l'état actuel, où l'on peut observer une interface de glissement (ici une faille). La représentation maillée de droite représente le même bassin sédimentaire à un temps géologique t, antérieur à l'actuel. A ce 30 temps t, les couches sédimentaires n'étaient pas encore fracturées. La représentation maillée centrale est un état intermédiaire, c'est-à-dire qu'elle représente le bassin sédimentaire à un temps t', compris entre le temps t et l'actuel. On constate que le glissement a commencé à modifier l'architecture du bassin. 3045842 16 2. Reconstitution de la formation du bassin Au cours de cette étape, il s'agit de reconstituer la formation du bassin étudié, à partir des représentations maillées établies aux étapes précédentes, au moyen d'un simulateur de 5 bassin et d'un simulateur géomécanique coopérant entre eux. Pour ce faire, il est nécessaire de disposer préalablement d'une reconstitution de la formation de la première des couches du bassin sédimentaire étudié, c'est-à-dire celle s'étant déposée la première dans le bassin, ou encore la couche la plus ancienne. Cette reconstitution de la première des couches peut avoir été obtenue par tout moyen. Selon un 10 mode de mise en oeuvre de la présente invention, cette reconstitution est obtenue par une simulation de bassin. Puis, pour au moins une couche s'étant déposée ultérieurement à la première couche (c'est-à-dire une couche située au-dessus de la première des couches), on applique au moins une simulation de bassin suivie d'une simulation géomécanique, lesdites simulations 15 prenant en compte conjointement la couche considérée et au moins une couche sous- jacente à la couche considérée, la coopération étant assurée par le transfert d'une partie des paramètres en sortie de la simulation de bassin à la simulation géomécanique. Préférentiellement, pour une couche donnée, la simulation de bassin est appliquée en coopération à la simulation géomécanique pour la couche considérée, conjointement à 20 toutes les couches précédemment déposées (c'est-à-dire à toutes les couches sous- jacentes à la couche considérée). Ainsi selon l'invention, on met en oeuvre les sous-étapes suivantes pour chaque couche considérée, les sous-étapes 2.1 à 2.3 étant répétées pour au moins chaque couche 25 considérée : 2.1. Application d'une simulation de bassin Au cours de cette sous-étape, il s'agit d'appliquer une simulation de bassin, au moyen d'un simulateur de bassin, à la couche considérée, conjointement à au moins une couche 30 sous-jacente, c'est-à-dire une couche déjà déposée. La simulation de bassin conjointe de 3045842 17 plusieurs couches sédimentaires permet de prendre en compte les contraintes imposées par le dépôt des couches sus-jacentes sur les couches sous-jacentes. Ainsi, la simulation de bassin appliquée selon l'invention permet de calculer des paramètres pour la couche considérée, mais permet aussi de mettre à jour des paramètres 5 calculés pour la couche précédemment déposée et considérée conjointement avec la couche en cours. Selon l'invention, la simulation de bassin pour la couche considérée est mise en oeuvre en utilisant la représentation maillée déterminée au cours de l'étape 2 précédente pour la période de temps correspondant au dépôt de la couche courante. Selon un mode de mise en 10 oeuvre de l'invention, ladite représentation maillée requise contient en chaque maille des informations usuelles en simulation de bassin, telles que des propriétés pétrophysiques relatives à la lithologie, des informations caractérisant le type de compaction, des propriétés relatives aux fluides présents dans la formation à la période de temps considérée, ainsi que les limites et géométries des couches déjà déposées. 15 Selon un mode de mise en oeuvre de l'invention, le simulateur de bassin utilisé pour la mise en oeuvre de la présente invention permet de discrétiser et résoudre les équations décrite dans le document (Schneider et al., 2000). Le logiciel TemisFlowTM (IFP Énergies nouvelles, France) est un exemple d'un tel simulateur de bassin. A l'issue de l'application de la simulation de bassin pour la période de temps 20 considérée, on obtient un premier jeu de paramètres, tels que des pressions de pores, des saturations, des températures, des porosités dans chacune des mailles de la représentation maillée représentative de la période considérée. Selon l'invention, ce premier jeu de paramètres comprend au moins les porosités en chaque maille de la représentation maillée pour la période considérée. Selon un mode de mise en oeuvre de l'invention, ce premier jeu 25 de paramètres comprend au moins les pressions de pores et les porosités en chaque maille de la représentation maillée pour la période considérée. 2.2 Application d'une simulation géomécanique Au cours de cette sous-étape, il s'agit de réaliser une simulation géomécanique, pour 30 la même période de temps que celle considérée pour la sous-étape précédente. Selon l'invention, la simulation géomécanique est appliquée conjointement à la couche s'étant déposée lors de la période en cours et à au moins à une autre couche s'étant déposée à une 3045842 18 période de temps antérieure. La simulation géomécanique conjointe de plusieurs couches sédimentaires permet de prendre en compte les contraintes mécaniques imposées par le dépôt des couches sus-jacentes sur les couches sous-jacentes. Préférentiellement, on appliquera la simulation géomécanique à l'ensemble des 5 couches déposées jusqu'à la période de temps considérée. Autrement dit, la simulation géomécanique à une période de temps donnée est mise en oeuvre en incluant toute la colonne de couches sédimentaires déjà déposées jusqu'à cette période. De cette façon, le procédé selon l'invention permet de non seulement simuler les effets géomécaniques tridimensionnels qui se sont produits dans la couche déposée pendant la période de temps 10 considérée, mais aussi les effets produits sur toutes les couches déjà déposées jusqu'à la période de temps considérée. Selon l'invention, la simulation géomécanique est réalisée à partir d'une représentation maillée représentative du bassin pour la période de temps considérée. Selon un mode de mise en oeuvre de l'invention, il peut s'agir de la représentation maillée utilisée pour la 15 simulation de bassin lors de la mise en oeuvre de la sous-étape précédente, éventuellement adaptée à la méthode de résolution des équations requises par la simulation géomécanique comme décrit précédemment. Selon l'invention, une partie du premier jeu de paramètres, calculé à la sous-étape précédente pour la même période de temps, est transmis au simulateur géomécanique. 20 Selon un mode de mise en oeuvre de l'invention, une partie du premier jeu de paramètres calculé à la sous-étape précédente correspond à la pression au début et à la fin de la période de temps considérée (on parle d'historique de pressions) par la simulation de bassin. Cet historique de pressions peut alors être utilisé comme une contrainte pour la simulation géomécanique relative à la période de temps considérée. Selon un autre mode de 25 mise en oeuvre, on transmet en outre à la simulation géomécanique les températures et/ou les saturations déterminées par simulation de bassin à la sous-étape précédente. Selon un mode de mise en oeuvre de l'invention, dans le cas de représentations maillées différentes entre simulation de bassin et simulation géomécanique, on peut assurer la correspondance entre les deux représentations maillées via un algorithme dit de « mapping » 30 (« cartographie » en français). Préférentiellement, la simulation géomécanique prend aussi en compte des contraintes mécaniques liées aux mouvements tectoniques régionaux subis par le bassin sédimentaire étudié. Selon un mode de mise en oeuvre de l'invention, ces contraintes mécaniques sont 3045842 19 prises en compte comme conditions aux limites lors des simulations géomécaniques, imposant par exemple à la simulation géomécanique de vérifier, sur les bords du domaine de simulation, des déplacements horizontaux (traduisant une compression ou une extension). La prise en compte de telles conditions aux limites peut être réalisée selon la méthode des 5 éléments finis via l'utilisation de conditions aux limites en déplacement ou contraintes. Ainsi, selon ce mode de réalisation, la représentation maillée pour la période de temps considérée se déforme au cours de l'étape de simulation géomécanique, en respectant les conditions aux limites établies. A l'issue de cette sous-étape, on obtient un deuxième jeu de paramètres relatif à la 10 période de temps considérée. Selon un mode de mise en oeuvre de l'invention, ce deuxième jeu de paramètres comprend au moins les porosités en chaque maille de la représentation maillée pour la période considérée. A noter que la représentation maillée à la fin d'un pas de temps d'une simulation géomécanique n'est pas nécessairement identique à la représentation maillée issue de l'étape de reconstruction structurale et déterminée au cours 15 de l'étape 2 pour le pas de temps suivant. 2.3 Contrôle de cohérence Au cours de cette étape, il s'agit d'effectuer un contrôle de cohérence entre au moins une partie des paramètres du premier jeu de paramètres, déterminé à l'étape 2.1 par 20 simulation de bassin, et au moins une partie des paramètres du deuxième jeu de paramètres, déterminé à l'étape 2.2 par simulation géomécanique. Pour ce faire, selon l'invention, on mesure un écart entre au moins une partie des paramètres du premier jeu (issu de la simulation de bassin) et au moins une partie des paramètres du deuxième jeu (issu de la simulation géomécanique) et on vérifie si cet écart s'inscrit ou non au-delà d'un 25 certain seuil prédéfini par le spécialiste. Selon un mode de mise en oeuvre de l'invention, ledit écart est évalué en comparant les porosités estimées par simulation de bassin et les porosités estimées par simulation géomécanique. L'écart peut être estimé maille à maille si les modèles bassin et géomécaniques ont des maillages identiques, sinon, on a recours à une technique de 30 transfert de champs, bien connue du spécialiste. On trouvera par exemple une référence concernant une telle technique dans le document (Nayroles et al., 1991) Selon un mode de mise en oeuvre de l'invention, on calcule, pour le pas de temps t considéré, un écart MES, selon une formule du type 3045842 20 MESrs = maxnENO(Pn,t etl) (1) où (pt (respectivement cent) est la porosité déterminée par simulation de bassin (respectivement pas simulation géomécanique) en une maille n de la représentation maillée relative au pas de temps t., la représentation maillée comprenant N mailles. A noter que si le 5 simulateur géomécanique ne détermine pas directement la porosité, celle-ci peut être obtenue par la formule : Mme (177. - Vs )/177. avec avec T correspondant à au volume total de la maille et Vs au volume solide dans la maille. Selon ce mode de mise en oeuvre, on considère que les prédictions de porosité sont cohérentes entre simulation géomécanique et simulation de bassin pour le pas de temps t considéré lorsque l'écart entre ces deux 10 prédictions est inférieure à une valeur de seuil prédéfinie par le spécialiste. Avantageusement, la valeur de seuil prédéfinie peut être comprise en 0,0001 et 0,02. Préférentiellement, on considère que les prédictions de porosité sont cohérentes entre simulation géomécanique et simulation de bassin pour le pas de temps t considéré lorsque MESt REVENDICATIONS1- Procédé d'exploitation d'un bassin sédimentaire comportant des hydrocarbures, comprenant une reconstitution de la formation dudit bassin, ladite formation comprenant au moins le dépôt de deux couches sédimentaires, ladite reconstitution étant réalisée au moyen d'un simulateur de bassin et d'un simulateur géomécanique coopérant entre eux, au moyen de mesures de propriétés relatives audit bassin et de représentations maillées représentatives dudit bassin au moment du dépôt de chacune desdites couches sédimentaires, caractérisé en ce que, à partir desdites mesures et d'une reconstitution de ladite formation de la première desdites couches, ladite reconstitution pour au moins une desdites couches sus-jacentes à ladite première couche est mise en oeuvre selon au moins les étapes suivantes : A. au moyen dudit simulateur de bassin, on réalise une simulation de bassin conjointe de ladite couche et d'au moins une couche sous-jacente à ladite couche, et on détermine un premier jeu de paramètres ; B. à partir d'au moins une partie dudit premier jeu desdits paramètres et dudit simulateur géomécanique, on réalise une simulation géomécanique conjointe de ladite couche et d'au moins une couche sous-jacente à ladite couche, et on détermine un deuxième jeu de paramètres ; C. on mesure un écart entre au moins une partie desdits paramètres dudit premier jeu et au moins une partie desdits paramètres dudit deuxième jeu, et on réitère les étapes A) à C) en appliquant une correction à ladite simulation de bassin si ledit écart est supérieur à un seuil prédéfini ; et en ce que, au moyen de ladite reconstitution réalisée pour lesdites couches, on sélectionne au moins une zone dudit bassin comportant lesdits hydrocarbures, et on exploite ledit bassin en fonction de ladite zone sélectionnée. 2- Procédé selon la revendication 1, dans lequel ladite reconstitution relative à ladite première couche est réalisée selon au moins les étapes suivantes : 3045842 28 i. au moyen dudit simulateur de bassin, on réalise une simulation de bassin relative à ladite première couche, et on détermine un premier jeu de paramètres ; ii. à partir d'au moins une partie dudit premier jeu desdits paramètres et dudit 5 simulateur géomécanique, on réalise une simulation relative à ladite première couche, et on détermine un deuxième jeu de paramètres ; iii. on mesure un écart entre au moins une partie desdits paramètres dudit premier jeu et au moins une partie desdits paramètres dudit deuxième jeu, et on réitère les étapes i) à iii) en appliquant une correction à ladite simulation de 10 bassin si ledit écart est supérieur à un seuil prédéfini. 3- Procédé selon l'une des revendications précédentes, dans lequel ledit premier jeu de paramètres comprend au moins la porosité et la pression en chacune des mailles de ladite représentation maillée. 15 4- Procédé selon la revendication 3, dans lequel ledit deuxième jeu de paramètres comprend au moins la porosité en chacune des mailles de ladite représentation maillée. 20 5- Procédé selon l'une des revendications 3 à 4, dans lequel ledit écart est fondé sur une mesure de la différence entre ladite porosité issue dudit premier jeu de et ladite porosité issue dudit deuxième jeu. 6- Procédé selon la revendication 5, dans lequel ledit écart est un écart absolu 25 MESabs défini selon une formule du type : MESabs = MaXnENO(Pe (PP2.1) où (/),,P (respectivement qqn) correspond à ladite porosité déterminée par ledit simulateur de bassin (respectivement par ledit simulateur géomécanique) en une maille n de ladite représentation maillée comprenant N mailles. 30 3045842 29 7- Procédé selon la revendication 5, dans lequel ledit écart est un écart relatif MESrel défini selon une formule du type : mEsrel 2 xmaxnENG (PZ -(Pri,l) = (Pn9-E(Pnb où (/),,P (respectivement (pe) correspond à ladite porosité déterminée par ledit 5 simulateur de bassin (respectivement par ledit simulateur géomécanique) en une maille n de ladite représentation maillée comprenant N mailles. 8- Procédé selon l'une des revendications précédentes, dans lequel ledit procédé comprend une étape de simulation de l'érosion d'au moins une desdites couches 10 et/ou de simulation d'une lacune géologique. 9- 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 15 oeuvre du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.