Brevet FR3045098A1 — publié le 20170616

La présente invention concerne le domaine de l'exploration pétrolière. L'exploration pétrolière consiste à rechercher des gisements d'hydrocarbures dans un bassin sédimentaire. La démarche générale comporte des aller-retour entre : des expertises et des études géologiques d'évaluation du potentiel pétrolier du bassin sédimentaire, réalisées à partir des données disponibles (affleurements, campagnes sismiques, forages). Cette expertise vise : o à mieux comprendre l'architecture et l'histoire géologique du sous-sol, notamment à étudier si des processus de maturation et de migration d'hydrocarbures ont pu se mettre en place ; o à identifier les zones du sous-sol dans lesquelles ces hydrocarbures ont pu s'accumuler ; o à établir quelles zones présentent le meilleur potentiel économique, évalué à partir du volume et de la nature des hydrocarbures probablement piégés (viscosité, taux de mélange avec de l'eau, composition chimique ...), ainsi que de leur coût d'exploitation (contrôlé par exemple par la profondeur et la pression de fluide). des forages d'exploration dans les différentes zones 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. Dans certains bassins sédimentaires, ayant subi une histoire géologique compliquée faisant interagir de nombreux processus physiques, ou lorsque le volume de données est très important, une simple intervention humaine d'expertise ne suffit pas à prédire l'emplacement et l'intérêt économique des gisements. On met alors en oeuvre une démarche faisant appel à des outils informatiques de synthèse des données disponibles et de 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". Elle fournit une cartographie prédictive du sous-sol indiquant l'emplacement probable des gisements, ainsi que la teneur, la nature et la pression des hydrocarbures qui y sont piégés. La compréhension des principes de la genèse des hydrocarbures et de leurs liens avec l'histoire géologique du sous-sol a permis la mise au point de méthodes de prédiction du 3045098 2 potentiel pétrolier et de l'emplacement des gisements dans les bassins sédimentaires. Ces méthodes s'appuyaient sur des observations géologiques effectuées sur le terrain, et des géologues experts intégraient leurs connaissances dans un raisonnement visant à fournir des hypothèses de scénarios de formation, migration et piégeage des hydrocarbures pour 5 chaque bassin étudié. Les forages d'exploration étaient ensuite réalisés suivant leurs recommandations, avec un taux de succès de l'ordre d'une découverte de gisement pour 10 puits forés. Les informations fournies par cette phase d'expertise étaient essentiellement qualitatives. Dans les années 80, le développement de l'informatique a fourni de nouveaux outils 10 pour assister les géologues dans cette démarche, les aider à fournir des informations quantitatives et plus fiables, et leur permettre ainsi d'augmenter ce taux de succès. De nombreux codes informatiques ont été développés pour mieux appréhender et quantifier chacun des phénomènes géologiques contrôlant la formation, la migration et le piégeage des hydrocarbures. Parmi ces outils, on peut citer les logiciels dits "de modélisation (ou de 15 simulation) de bassin". La modélisation de bassin Les logiciels de modélisation de bassin permettent de simuler en une, deux ou trois dimensions, l'ensemble des processus sédimentaires, tectoniques, thermiques, 20 hydrodynamiques et de chimie organique et inorganique qui interviennent lors de la formation d'un bassin pétrolier. Ces logiciels constituent un outil fondamental pour l'exploration pétrolière, car ils fournissent en sortie l'ensemble des données nécessaires à la prédiction de la localisation des gisements et de leur intérêt économique. La démarche suivie par la plupart des modèles de bassin comporte trois phases : 25 1. une phase de construction d'un maillage du sous-sol suivant une hypothèse sur son architecture interne et sur les propriétés qui caractérisent chaque maille : par exemple leur porosité, leur nature sédimentaire (argile, sable ...) ou encore leur teneur en matière organique au moment de leur sédimentation. La construction de ce modèle se base sur des données acquises par campagnes sismiques ou 30 mesures lors du forage. Ce maillage est structuré en couches : un groupe de mailles est affecté à chaque couche géologique du bassin modélisé. 2. une phase de reconstruction du maillage représentant des états antérieurs de l'architecture du bassin. Cette étape est réalisée à l'aide, par exemple, d'une méthode de "backstripping" (Steckler, M.S., and A.B. Watts, Subsidence of the Atlantic-type continental margin off New York, Earth Planet. Sci. Lett., 41, 1-13, 3045098 3 1978.) ou une méthode de restauration structurale (décrit par exemple dans la demande de brevet FR 2930350 - US 8150669). 3. une étape de simulation numérique d'une sélection de phénomènes physiques se déroulant au cours de l'évolution du bassin et contribuant à la formation des 5 pièges pétroliers. Cette étape s'appuie sur une représentation discrétisée du temps en « évènements », chaque évènement étant simulé par une succession de pas de temps. Le début et la fin d'un évènement correspondent à deux états successifs de l'évolution de l'architecture du bassin fournis à l'étape 2 précédente. 10 Genèse des hydrocarbures Parmi les phénomènes physiques pris en compte dans l'étape 3 ci-dessus, on simule notamment la transformation en hydrocarbures de la matière organique initialement enfouie avec les sédiments, et le cheminement de ces hydrocarbures depuis les roches, dans lesquelles ils se sont formés jusqu'à celles où ils sont piégés. 15 Les procédés de modélisation de bassin proposés actuellement permettent de simuler trois processus possibles de formation d'hydrocarbures : par action directe de certains micro-organismes présents en profondeur dans les sédiments qui dégradent la matière organique initiale en produisant notamment du méthane ; il s'agit de la production biogénique primaire (décrit par exemple dans la 20 demande de brevet WO 2014/040622 Al). par transformation chimique et/ou thermodynamique de la matière organique sous l'effet des fortes températures qui règnent dans le sous-sol profond (décrit par exemple dans la demande de brevet WO 2014/040622 Al ). par action de certains micro-organismes du sous-sol profond qui dégradent 25 certains hydrocarbures en produisant notamment du méthane : il s'agit de la production biogénique à partir d'hydrocarbures fluides (décrit par exemple dans la demande de brevet FR 2888251 - US 8753867) Migration des hydrocarbures 30 Ces mêmes logiciels permettent également de simuler le transport des hydrocarbures produits, notamment le méthane biogénique primaire et secondaire, en faisant l'hypothèse que ceux-ci se déplacent sous la forme d'une ou plusieurs phases, séparées de la phase aqueuse, que l'on nommera « phases hydrocarbures » ou « gaz libre ». Plusieurs méthodes 3045098 4 sont disponibles pour modéliser la migration des phases hydrocarbures : la méthode par lancer de rayon (Sylta, Modeling techniques for hydrocarbon migration, In EAGE 64th Conference and Exhibition, Florence, 2002), la méthode par invasion percolation (Carruthers, Transport of secondary oil migration using gradient-driven invasion percolation 5 techniques. PhD thesis, Heriot-Watt University, Edinburgh, Scotland, UK, 1998), et la résolution des équations de Darcy généralisées (Schneider, Modelling multi-phase flow of petroleum at the sedimentary basin scale. Journal of Geochemical exploration 78-79 (2003) 693-696). 10 Les hydrocarbures piégés dans les réservoirs pétroliers sont majoritairement issus de la transformation en hydrocarbures de la matière organique solide et ce principalement sous l'effet de la chaleur. Cette dernière est souvent caractérisée par sa teneur en carbone organique (TOC) et son indice d'hydrogène (IH) qui permettent de connaître le potentiel pétrolier de la roche mère (i.e. sa capacité à produire des hydrocarbures, sa richesse). Les 15 valeurs de IH et de TOC sont classiquement évaluées à l'état actuel (i.e. après possible altération thermique subie dans les conditions géologiques), grâce à des outils tels que le Rock-Eval ® (IFP Energies nouvelles, France). Les modèles de bassin ont parmi leurs objectifs, l'estimation des quantités d'hydrocarbures générés par les roches-mères et piégés dans les réservoirs pétroliers. Ces 20 quantités dépendent ainsi directement de la distribution spatiale de la matière organique et de sa richesse. Il est toutefois très difficile d'avoir une représentation juste de ces paramètres, car les roches-mères (qui contiennent la matière organique) ne sont connues que lorsqu'elles affleurent (elles peuvent alors être relativement bien caractérisées, mais leur caractérisation revêt une importance moindre, car il est alors peu probable qu'elles aient 25 généré des hydrocarbures, ou que ces derniers soient encore exploitables), ou lorsque des puits de forage pétrolier ont pu l'atteindre. Pour des réservoirs dits « conventionnels », les forages pétroliers ne vont malheureusement que très peu souvent jusqu'à la roche-mère, ce qui rend l'estimation de la distribution de la richesse de la matière organique difficile. Même dans le cas des réservoirs dits « non-conventionnels », pour lequel le réservoir est 30 généralement la roche-mère elle-même, il reste difficile de l'estimer correctement, la densité des puits n'étant pas toujours suffisante pour avoir une image globale de la distribution de la matière organique. En modélisation de bassin, les données concernant la richesse de la matière organique sont indispensables à la simulation de la génération et de la migration des 35 hydrocarbures, les utilisateurs sont donc obligés de construire des cartes de TOC et de IH à 3045098 5 partir des quelques données en leur possession. Les cartes sont alors souvent uniformes (même valeur sur l'ensemble de la carte) ou obtenues par interpolation / extrapolation mathématique (méthodes de krigeage, de « smoothing »...) à partir de quelques points. Les processus responsables du dépôt de la matière organique ne sont ainsi pas pris en compte 5 dans l'estimation des cartes de richesse. Une autre solution tenant compte des processus liés au dépôt de la matière organique consiste à utiliser des outils de modélisation des dépôts sédimentaires tels que DionisosFlow 0 (IFP Energies nouvelles). Ces outils demandent toutefois de gros investissements en temps car ils requièrent de nombreuses données additionnelles, des simulations 10 supplémentaires et des compétences que n'ont pas nécessairement les utilisateurs de modélisation de bassin. Pour pallier ces inconvénients, la présente invention concerne un procédé d'exploitation d'un bassin sédimentaire. L'exploitation du bassin sédimentaire est permise 15 par le choix de zones d'intérêt au sein du bassin, les zones étant déterminées au moyen de cartes de teneur en carbone organique TOC et d'indice d'hydrogène IH d'au moins une couche sédimentaire. Selon l'invention, les cartes de teneur en carbone organique TOC et d'indice d'hydrogène IH sont obtenues en prenant en compte des valeurs actuelles de TOC et IH au niveau du puits et en prenant en compte les cartes de bathymétrie et de vitesse de 20 sédimentation. Ainsi, le procédé selon l'invention permet de déterminer des cartes de richesse TOC et IH avant que ces propriétés ne soient affectés par les processus thermiques sur la base des informations dont dispose le modélisateur de bassin (généralement TOC et IH en quelques puits, généralement affectés par les processus thermiques) et des processus physiques et chimiques (cartes de sédimentation et de bathymétrie) sans avoir recours à des 25 simulations de dépôt sédimentaire, qui nécessiteraient plus de temps, de données et de compétences. Le procédé selon l'invention L'invention concerne un procédé d'exploitation d'un bassin sédimentaire, ledit bassin 30 sédimentaire comprenant au moins une couche sédimentaire, et étant traversé par au moins un puits, dans lequel on connaît la valeur de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au temps actuel au niveau d'au moins un puits, une carte de bathymétrie et une carte de vitesse de sédimentation de ladite couche sédimentaire. Pour ce procédé, on réalise les étapes suivantes : 3045098 6 a) on construit pour au moins une couche sédimentaire des cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire à partir desdites valeurs de teneur en carbone organique TOCA(P) et de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits, et au moyen 5 desdites cartes de bathymétrie et de vitesse de sédimentation de ladite couche sédimentaire ; b) on détermine au moins une zone dudit bassin sédimentaire comportant des hydrocarbures au moyen desdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire ; et 10 c) on exploite ladite zone dudit bassin sédimentaire. Selon un mode de réalisation de l'invention, on construit pour au moins une couche sédimentaire des cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire en mettant en oeuvre les étapes suivantes : 15 i) on construit des modèles de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 en tout point de ladite couche sédimentaire au temps du dépôt de ladite couche sédimentaire, lesdits modèles reliant ladite teneur en carbone organique TOCS et ledit indice d'hydrogène IH0 en tout point de ladite couche sédimentaire au temps du dépôt de ladite couche sédimentaire, avec lesdites 20 cartes de bathymétrie et de vitesse de sédimentation de ladite couche sédimentaire et avec des paramètres chimiques et biologiques de ladite couche sédimentaire ; ii) on cale lesdits modèles en déterminant lesdits paramètres chimiques et biologiques au moyen desdites valeurs actuelles de teneur en carbone organique 25 TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au niveau dudit puits ; et iii) on construit lesdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire en appliquant lesdits modèles calés en toute maille de ladite couche sédimentaire au moyen des paramètres chimiques et biologiques déterminés et desdites cartes de vitesse 30 de sédimentation et de bathymétrie. Avantageusement, on cale lesdits modèles en mettant en oeuvre les étapes suivantes : (1) on choisit des paramètres biologiques et chimiques de ladite couche sédimentaire ; 35 (2) on détermine des valeurs de teneur en carbone organique T000s(P) et d'indice d'hydrogène IH0s(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire au moyen desdits modèles, desdits paramètres 3045098 7 biologiques et chimiques choisis et desdites cartes de bathymétrie et de vitesse de sédimentation ; (3) on calcule une valeur de teneur en carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire au moyen d'un 5 bilan carbone et desdites valeurs de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits ; (4) on compare lesdites valeurs déterminées TOCos(P) et calculées T000Bc(P) de teneur en carbone organique au niveau dudit puits au temps de dépôt de ladite couche sédimentaire, et on réitère les étapes (2) à (4) en modifiant 10 lesdits paramètres biologiques et chimiques pour minimiser ladite comparaison. Selon une variante de réalisation, on calcule la valeur de teneur en carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire au moyen TOCA(P) d'une formule du type : TOCOBC (P) =1 IHos(P)-1HA(B)>(0.84 1000 15 Alternativement, on calcule la valeur de teneur en carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire au moyen d'une formule TOCA(P)-0.84XS1 du type : TOCOBC (P) =1 /Hos(P)-/HA(B) avec 51 la quantité d'huiles libres. >(0.84 1000 De préférence, ladite valeur de teneur en carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire est fonction d'un paramètre 20 relatif à la formation de gaz biogénique au sein de ladite couche sédimentaire. Selon une mise en oeuvre de l'invention, lesdits paramètres biologiques et chimiques sont des paramètres représentatifs du dépôt de matière organique et/ou des réactions de réduction et/ou d'oxydation. Conformément à un mode de réalisation de l'invention, lesdits modèles sont formés par 25 interpolation desdites valeurs de teneur en carbone organique TOC° (P) et d'indice d'hydrogène IHA(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire, l'interpolation étant contrainte par lesdites cartes de vitesse de sédimentation et de bathymétrie. De manière avantageuse, lesdites valeurs de teneur en carbone organique TOCA(P) et 30 de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits sont déterminées par analyse d'une pyrolyse d'une roche de ladite couche sédimentaire. Conformément à un mode de réalisation de l'invention, on détermine au moins une zone dudit bassin sédimentaire comportant des hydrocarbures au moyen d'un modèle dudit 3045098 8 bassin sédimentaire construit au moyen desdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire. En outre, l'invention concerne un produit programme d'ordinateur téléchargeable 5 depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des caractéristiques précédentes, lorsque ledit programme est exécuté sur un ordinateur. 10 Présentation succincte des figures 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. La figure 1 illustre les étapes du procédé d'exploitation selon un mode de réalisation de 15 l'invention. La figure 2 illustre les étapes du calage du modèle selon un mode de réalisation de l'invention. La figure 3 illustre une carte de bathymétrie pour un exemple. La figure 4 illustre une carte de vitesse de sédimentation pour l'exemple de la figure 3. 20 La figure 5 illustre une carte de teneur en carbone organique TOCS au niveau de la couche sédimentaire pour l'exemple de la figure 3, la carte étant obtenue par le procédé selon l'invention. La figure 6 illustre une carte d'indice d'hydrogène IH0 au niveau de la couche sédimentaire pour l'exemple de la figure 3, la carte étant obtenue par le procédé selon 25 l'invention. Description détaillée de l'invention La présente invention concerne un procédé d'exploitation d'un bassin sédimentaire. Le bassin sédimentaire comprend au moins une couche sédimentaire, et est traversé par au 30 moins un puits, par exemple un puits d'exploration, un puits injecteur ou un puits producteur. Notations : Au cours de la description, les notations suivantes seront utilisées : TOC : teneur en carbone, avec : 35 - TOCA : valeur actuelle de la teneur en carbone, 3045098 9 TOC° : valeur de teneur en carbone lors du dépôt de la couche sédimentaire, TOCos : valeur de teneur en carbone simulée obtenue par le modèle selon l'invention, la valeur étant celle lors du dépôt de la couche sédimentaire, TOCosc : valeur de teneur en carbone calculée obtenue par un bilan carbone, 5 la valeur étant celle lors du dépôt de la couche sédimentaire, IH : indice d'hydrogène, avec : IHA : valeur actuelle de l'indice d'hydrogène, IH0 : valeur de l'indice d'hydrogène lors du dépôt de la couche sédimentaire, IHos : valeur de l'indice d'hydrogène simulée obtenue par le modèle selon 10 l'invention, la valeur étant celle lors du dépôt de la couche sédimentaire. Ces notations, lorsqu'elles sont suivies de la notation (P), désignent les valeurs au niveau d'un puits. Lorsqu'elles ne sont pas suivies de la notation (P), désignent les valeurs en tout point de la couche sédimentaire. On appelle valeur actuelle ou valeur au temps actuelle, la valeur obtenue par mesure à 15 l'instant présent (au pas de temps actuel), auquel on met en oeuvre le procédé selon l'invention, la roche peut avoir subie des altérations thermiques dans les conditions géologiques. A contrario, la valeur lors du dépôt ou au temps du dépôt de la couche sédimentaire, correspond à la valeur au temps géologique (c'est-à-dire au pas de temps géologique), pour lequel la couche sédimentaire s'est formée. 20 Le procédé d'exploitation selon l'invention utilise les informations suivantes : TOCA(P) : valeur actuelle de la teneur en carbone au niveau d'au moins un puits, cette valeur peut être obtenue grâce à des outils d'analyse de roche prélevée dans le bassin sédimentaire, tel que le Rock-Eval 0 (IFP Energies 25 nouvelles, France) ou tout outil analogue, IHA(P) : valeur actuelle de l'indice d'hydrogène au niveau d'au moins un puits, cette valeur peut être obtenue grâce à des outils d'analyse de roche prélevée dans le bassin sédimentaire, tel que le Rock-Eval 0 (IFP Energies nouvelles, France) ou tout outil analogue, 30 une carte de bathymétrie de la couche sédimentaire, qui est une carte qui indique la profondeur des mailles de la carte par rapport au niveau de la mer, cette carte est généralement utilisée lors de la construction de modèle de bassin, la figure 3, qui sera détaillée dans la suite de la description, est un exemple de carte de bathymétrie, et 35 une carte de vitesse de sédimentation, qui est une carte qui est représentative de la durée du dépôt et de l'épaisseur de l'évènement, cette carte est généralement utilisée lors de la construction de modèle de bassin, la figure 4, 3045098 10 qui sera détaillée dans la suite de la description, est un exemple de carte de vitesse de sédimentation. On appelle carte, une représentation maillée de la couche sédimentaire, elle 5 représente une maquette de la couche sédimentaire. A chaque maille est associée une propriété de la maille. Par exemple, la propriété peut être la bathymétrie, la vitesse de sédimentation, la teneur en carbone organique, l'indice d'hydrogène, etc. Classiquement, la carte peut être une représentation numérique, réalisée sur ordinateur. On rappelle que le Rock-Eval 0 (IFP Energies nouvelles, France) est un système 10 mettant en oeuvre une pyrolyse d'une roche issue du sous-sol. Elle fournit, en effet, et d'une façon rapide, différentes informations sur le contenu organique des roches, telles que le potentiel pétrolier des séries rencontrées, la nature des kérogènes, leur état de maturation. Le système est la méthode sont notamment décrits dans la demande de brevet FR 2722296 (US 5843787). Par exemple, le Rock-Eval 0 (IFP Energies nouvelles, France) peut délivrer 15 des paramètres classiques tels que : - la quantité de gaz libres: 50 ; - la quantité d'huiles libres: 51' ; - le potentiel pétrolier, ou hydrocarbures issus de pyrolyse: S2 ; - la température du sommet du pic S2: TMax ; 20 - le carbone organique total, TOC ; et - l'index d'hydrogène IH Le procédé selon l'invention comprend les étapes suivantes : 1) construction des cartes de TOC et d'IH pour au moins une couche sédimentaire, 25 2) détermination d'une zone d'intérêt, et 3) exploitation du bassin sédimentaire. La figure 1 illustre de manière schématique et non limitative les étapes du procédé selon un mode de réalisation de l'invention. Pour le procédé illustré, les informations 30 suivantes sont connues : la valeur actuelle de la teneur en carbone au niveau d'un puits TOCA(P), valeur actuelle de l'indice d'hydrogène au niveau d'un puits IHA(P), la carte de bathymétrie CB et la carte de vitesse de sédimentation CS. Une première étape consiste en une étape de construction CON des cartes de teneur en carbone organique TOC et d'indice hydrogène IH pour au moins une couche sédimentaire. Selon un mode de réalisation de 35 l'invention, détaillé dans la suite de la description, l'étape de construction CON peut comprendre une étape de modélisation MOD, une étape de calage des modèles CAL, et une étape d'application des modèles pour former les cartes CAR. Après l'étape de construction 3045098 11 CON, le procédé selon l'invention peut comprendre une étape détermination d'au moins une zone d'intérêt ZI, en fonction des cartes de TOC et d'IH construites, par exemple au moyen d'une modélisation de réservoir. Enfin, le procédé selon l'invention comprend une étape d'exploitation du bassin sédimentaire EXP en fonction de la zone d'intérêt déterminée ZI. 5 Etape 1) Construction des cartes de TOC d'IH Lors de cette étape, on construit les cartes de teneur en carbone organique TOC, et d'indice d'hydrogène IH pour au moins une couche sédimentaire du bassin sédimentaire. Les valeurs de TOC et IH des cartes étant celles déterminées au temps du dépôt de la 10 couche sédimentaire. De préférence, on construit les cartes de TOC et d'IH pour chaque couche sédimentaire du bassin sédimentaire pouvant comprendre des hydrocarbures. Cette étape 1) est alors répétée pour chaque couche du bassin sédimentaire à étudier. Cette étape peut être réalisée numériquement. 15 Selon une mise en oeuvre de cette étape, on peut construire les cartes de TOC et d'IH en mettant en oeuvre les étapes suivantes : i) on construit des modèles de teneur en carbone organique TOC° et d'indice d'hydrogène IH0 en tout point de la couche sédimentaire au temps du dépôt de la couche sédimentaire, les modèles reliant la teneur en carbone organique TOC° et 20 l'indice d'hydrogène IH0 en tout point de ladite couche sédimentaire au temps du dépôt de ladite couche sédimentaire, avec les cartes de bathymétrie et de vitesse de sédimentation de la couche sédimentaire, et avec des paramètres chimiques et biologiques de la couche sédimentaire ; ii) on cale les modèles en déterminant les paramètres chimiques et biologiques au 25 moyen des valeurs actuelles de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au niveau du puits ; et iii) on construit les cartes de teneur en carbone organique TOC° et d'indice d'hydrogène IH0 au temps du dépôt de la couche sédimentaire en appliquant les modèles calés en tout point du bassin sédimentaire au moyen des paramètres 30 chimiques et biologiques déterminés et des cartes de vitesse de sédimentation et de bathymétrie. Les paramètres biologiques et chimiques des modèles de TOC et d'IH peuvent être des paramètres représentatifs du dépôt de matière organique et/ou des réactions de réduction et/ou d'oxydation de la matière organique. 35 Les modèles prennent donc en compte des données relatives à la bathymétrie et à la vitesse de sédimentation, qui sont variables au sein de la couche sédimentaire. Par conséquent, les modèles peuvent donc être représentatifs de la couche sédimentaire, en 3045098 12 prenant en compte ces phénomènes géologiques. De plus, les modèles prennent en compte les conditions biologiques et chimiques, ce qui permet de rendre ces modèles cohérents avec les phénomènes biologiques et chimiques au sein de la couche sédimentaire. Ainsi, le modèle de teneur en carbone organique TOCS en tout point de la couche 5 sédimentaire au temps du dépôt de la couche sédimentaire peut s'écrire sous la forme d'une fonction f : T000 = f(CB,CS,PB,PC) avec CB la carte de bathymétrie (connue), CS la carte de vitesse de sédimentation (connue), PB les paramètres biologiques (à caler) et PC les paramètres chimiques (à caler). De la même manière, le modèle d'indice d'hydrogène IH0 en tout point de la couche 10 sédimentaire au temps du dépôt de la couche sédimentaire peut s'écrire sous la forme d'une fonction g : 1H0 = g (CB,CS,PB,PC) avec CB la carte de bathymétrie (connue), CS la carte de vitesse de sédimentation (connue), PB les paramètres biologiques (à caler) et PC les paramètres chimiques (à caler). 15 Le but de calage est de déterminer les paramètres biologiques et chimiques qui permettent de rendre cohérent le modèle avec les valeurs connues au niveau du puits. Le calage des modèles est réalisé au moyen des valeurs actuelles de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au niveau du puits. Ainsi, les modèles sont calés avec des données mesurées et permettent une cohérence entre les 20 modèles et les données réelles. Selon un mode de réalisation, pour lequel des valeurs actuelles de teneur en carbone organique TOCA(P) et de l'indice d'hydrogène IHA(P) sont connues au niveau de plusieurs puits, le calage peut être mis en oeuvre pour déterminer les paramètres biologiques et chimiques qui permettent de rendre cohérents le modèle avec les valeurs actuelles pour tous 25 les puits. La figure 2 illustre de manière schématique et non limitative, les étapes du calage du modèle. Dans un premier temps, on choisit des paramètres biologiques et chimiques PAR. Ces paramètres PAR sont utilisés pour simuler SIM, c'est-à-dire pour déterminer à l'aide des modèles, des valeurs de TOC et d'IH au niveau du puits lors du dépôt de la couche 30 sédimentaire. Ensuite, on calcule au moyen d'un bilan carbone BC, une valeur de TOC au niveau du puits lors du dépôt de la couche sédimentaire. Enfin, la valeur de TOC simulée et la valeur de TOC calculée sont comparées COMP. Si la différence entre la valeur simulée et la valeur calculée n'est pas minimale, alors les étapes précédentes sont répétées (flèche en pointillés) en modifiant les paramètres biologiques et chimiques. Dans le cas inverse, on 35 conserve les paramètres biologiques et chimiques qui minimisent la différence entre la valeur 3045098 13 simulée et la valeur calculée. Dans le cas de plusieurs puits, la sélection des paramètres biologiques et chimiques optimaux sont déduits en minimisant la somme des écarts, selon une méthode classique (par exemple les moindres carrés). 5 Selon une mise en oeuvre préférée de l'invention, on peut réaliser le calage en mettant en oeuvre les étapes suivantes : (1) on choisit des paramètres biologiques et chimiques de la couche sédimentaire ; (2) on détermine des valeurs de teneur en carbone organique T000s(P) et 10 d'indice d'hydrogène IHos(P) au niveau du puits au temps du dépôt de la couche sédimentaire au moyen des modèles, des paramètres biologiques et chimiques choisis à l'étape précédente ou modifiés pour la nouvelle itération, et des cartes de bathymétrie et de vitesse de sédimentation ; (3) on calcule une valeur de teneur en carbone organique T000Bc(P) au niveau 15 dudit puits au temps du dépôt de ladite couche sédimentaire au moyen d'un bilan carbone et des valeurs de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au temps actuel au niveau du puits ; (4) on compare les valeurs déterminées TOCos(P) et calculées TOCoBc(P) de teneur en carbone organique au niveau dudit puits au temps de dépôt de la 20 couche sédimentaire, et on réitère les étapes (2) à (4) en modifiant lesdits paramètres biologiques et chimiques pour minimiser la comparaison. Par exemple, si la comparaison est supérieure à un seuil, les paramètres biologiques et chimiques sont modifiés, et si la comparaison est inférieure à ce seuil, alors les paramètres biologiques et chimiques utilisés dans les 25 modèles sont ceux qui permettent cette comparaison. Pour la première itération, le choix des paramètres biologiques et chimiques peut être quelconque, ou dépendre de mesures. Pour les itérations suivantes, les paramètres biologiques et chimiques sont déterminés de manière à minimiser la différence entre la valeur simulée et la valeur calculée. 30 Pour l'étape de détermination des valeurs T000s(P) et IHos(P), on applique les modèles, au niveau du puits, avec les paramètres choisis et les cartes de vitesse de sédimentation et de bathymétrie connues. Pour l'étape de calcul de la valeur de T000Bc(P), on applique un bilan carbone au niveau du puits. Selon une première alternative, on peut calculer la valeur de teneur en 35 carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche 3045098 14 TOCA(P) sédimentaire au moyen d'une formule du type : T 0 COBC (P) - -IHos(p)-IFIA(B) 1 xC 1000 avec C une constante, correspondant au ratio en pourcent de carbone organique dans les hydrocarbures. Classiquement, la constante C peut valoir 0,84. Selon une deuxième alternative, on peut calculer la valeur de teneur en carbone organique T000Bc(P) au niveau 5 dudit puits au temps du dépôt de ladite couche sédimentaire au moyen d'une formule du TOCA (P)- 0.84 xS1 type : T000BC (P) =1 IHos(P)-IHA(B) avec 51 la quantité d'huiles libres, cette XC 1000 quantité peut être mesurée à partir de l'outil Rock-Eval 0 (IFP Energies nouvelles, France), et C une constante, correspondant au ratio en pourcent de carbone organique dans les hydrocarbures. Classiquement, la constante C peut valoir 0,84. Pour ces deux alternatives, 10 le bilan carbone dépend de la valeur déterminée de l'indice d'hydrogène II-Ios(P) et des données au niveau du puits. En outre, le bilan carbone utilisé pour le calcul de la valeur de T00000(P) peut être fonction de la formation du gaz biogénique, de manière à être plus complet. La comparaison de la valeur déterminée par les modèles et la valeur calculée peut 15 consister en une comparaison mathématique des valeurs, le but est de minimiser l'écart entre la valeur déterminée par les modèles et la valeur déterminée par bilan carbone. Si ces valeurs ne sont pas cohérentes alors les étapes sont réitérées en optimisant les paramètres biologiques et chimiques. Pour le mode de réalisation, pour lequel plusieurs valeurs actuelles de teneur en carbone organique TOCA(P) et de l'indice d'hydrogène IHA(P) sont connues au 20 niveau de plusieurs puits, la comparaison peut être fonction des comparaisons au niveau de chaque puits, par exemple au moyen d'une méthode des moindres carrés. A la fin du calage, les modèles de valeur de teneur en carbone organique et d'indice d'hydrogène sont entièrement paramétrés. On peut donc construire les cartes de TOC et 25 d'IH pour la couche sédimentaire en appliquant le modèle calé (avec les paramètres biologiques et chimiques calés et les cartes de bathymétrie et de vitesse de sédimentation connues) à toutes les mailles représentatives du bassin sédimentaire. L'avantage du procédé selon l'invention est qu'il permet ensuite de déterminer les cartes de TOC° et de IH0 à l'échelle du bassin sur la base des processus chimiques et 30 biologiques. Le procédé selon l'invention ne fait aucune hypothèse a priori sur la valeur de l'indice d'hydrogène IH0, qui est déterminée dans la procédure d'optimisation des paramètres chimiques et biologiques. Dans les approches classiques de l'art antérieur, une hypothèse sur le IH0 est nécessaire au calcul de la valeur de teneur en carbone organique TOC° à partir 3045098 15 des équations tirées des bilan carbones. Par manque d'information, un IH0 homogène spatialement dans le bassin est utilisé. Or, les conditions de vitesse de sédimentation et de bathymétrie variables dans le bassin sont connues et on peut voir qu'elles permettent d'aboutir à des variations spatiales de IH0 (comme illustré dans la fig. 4) si elles sont prises 5 en compte. Etape 2) Détermination d'une zone d'intérêt Au moyen des cartes de TOC et d'IH de la couche sédimentaire, on détermine au moins une zone d'intérêt du bassin sédimentaire, qui peut comporter des hydrocarbures. En 10 effet, la teneur en carbone organique et en hydrogène sont des indices de zones dans lesquelles, des hydrocarbures ont pu se former et/ou s'accumuler. De préférence, on met en oeuvre cette étape au moyen d'une modélisation de bassin, qui intègre les cartes de TOC et d'IH construites. Elle fournit une cartographie prédictive du 15 sous-sol indiquant l'emplacement probable des gisements, ainsi que la teneur, la nature et la pression des hydrocarbures qui y sont piégés, en prenant en compte la genèse, la migration et éventuellement la dégradation des hydrocarbures. Le modèle de bassin s'appuie sur une discrétisation spatiale et temporelle du bassin, il représente une maquette du bassin sédimentaire. La discrétisation spatiale est organisée en couches de mailles représentant les 20 couches géologiques. A chaque maille est associée au moins une propriété de la maille, par exemple la densité de roche, la porosité, la teneur en carbone organique TOC, l'indice d'hydrogène IH, etc. La discrétisation temporelle repose sur le découpage du temps en intervalles de temps élémentaires : les « pas de temps ». Une simulation complète est constituée d'un ensemble de séquences temporelles que l'on nomme « évènements », qui 25 sont elles-mêmes subdivisées en pas de temps. Un évènement dure donc au moins un pas de temps, et correspond à la mise en place d'une couche géologique. Classiquement, le modèle de bassin peut être une représentation numérique, réalisée sur ordinateur. Les logiciels de modélisation de bassin permettent de simuler en une, deux ou trois dimensions, l'ensemble des processus sédimentaires, tectoniques, thermiques, 30 hydrodynamiques et de chimie organique et inorganique qui interviennent lors de la formation d'un bassin pétrolier. Ces logiciels constituent un outil fondamental pour l'exploration pétrolière, car ils fournissent en sortie l'ensemble des données utiles à la prédiction de la localisation des gisements et de leur intérêt économique. La démarche suivie par la plupart des modèles de bassin comporte trois phases : 3045098 16 1. une phase de construction d'un maillage du sous-sol suivant une hypothèse sur son architecture interne et sur les propriétés qui caractérisent chaque maille : par exemple leur porosité, leur nature sédimentaire (argile, sable ...) ou encore leur teneur en matière organique au moment de leur sédimentation. La construction 5 de ce modèle se base sur des données acquises par campagnes sismiques ou mesures en forages. Ce maillage est structuré en couches : un groupe de mailles est affecté à chaque couche géologique du bassin modélisé. 2. une phase de reconstruction du maillage représentant des états antérieurs de l'architecture du bassin. Cette étape est réalisée à l'aide, par exemple, d'une 10 méthode de "backstripping" (Steckler, M.S., and A.B. Watts, Subsidence of the Atlantic-type continental margin off New York, Earth Planet. Sci. Lett., 41, 1-13, 1978.) ou une méthode de restauration structurale (décrit par exemple dans la demande de brevet FR 2930350 - US 8150669). Lors de cette phase, on intègre les cartes de TOC et d'IH construites pour la(es) couche(s) sédimentaire(s). 15 3. une étape de simulation numérique d'une sélection de phénomènes physiques se déroulant au cours de l'évolution du bassin et contribuant à la formation des pièges pétroliers. Cette étape s'appuie sur une représentation discrétisée du temps en « évènements », chaque évènement étant simulé par une succession de pas de temps. Le début et la fin d'un évènement correspondent à deux états 20 successifs de l'évolution de l'architecture du bassin fournis à l'étape 2 précédente. A partir la modélisation de bassin, le spécialiste peut déterminer au moins une zone du bassin, correspondant à des mailles de ladite représentation maillée au temps actuel 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 les 25 zones du bassin étudié présentant le meilleur potentiel pétrolier. Etape 3) Exploitation du bassin sédimentaire Lors de cette étape, on exploite la zone d'intérêt déterminée du bassin sédimentaire. L'exploitation pétrolière du bassin peut alors prendre plusieurs formes, notamment : 30 - la réalisation de forages d'exploration dans les différentes zones sélectionnées 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, 3045098 17 - 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. 5 En outre, l'invention concerne un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé tel que décrit précédemment, lorsque ledit programme est 10 exécuté sur un ordinateur. Exemple illustratif Le procédé selon l'invention est appliqué pour déterminer la richesse de la roche mère dans une couche sédimentaire du bassin Ouest Canadien. 15 Ce bassin fait actuellement l'objet d'une exploration pétrolière intense notamment dans les réservoirs « non-conventionnels » des formations Montney et Doig. La rentabilité des prospects non-conventionnels étant en partie liée à la richesse de la roche-mère (qui constitue également le réservoir des hydrocarbures), il est essentiel pour les producteurs pétroliers de correctement estimer la distribution du potentiel pétrolier dans le bassin avant de réaliser les 20 modélisations de bassin. Le procédé de détermination du TOC selon l'invention a donc été appliqué à ce bassin. La zone étudiée représente une surface de 500 000 km2 environ (500 km x 1000 km). Elle correspond aux couches sédimentaires de la base de la formation Montney, déposée au Trias sur une période de 2 millions d'années. Le bassin est discrétisé à l'aide de 5000 mailles. 25 Pour la couche considérée, on dispose de mesures de TOC et de IH en 7 puits dans le bassin (à comparer aux 50 x 100 soit 5000 mailles du modèle) et des cartes de bathymétrie (figure 3) et de vitesse de sédimentation (figure 4). Sur la figure 3, les zones Z1 ont une bathymétrie comprise entre -75 et -50 m, la zone Z2 une bathymétrie entre - 50 et -25 m, la zone Z3 une bathymétrie entre - 25 et 0 m, la zone 30 Z4 une bathymétrie entre 0 et 75 m, la zone Z5 entre 75 et 100 m, la zone Z6 entre 100 et 150 m, et la zone Z7 au-delà de 150 m. Sur la figure 4, les zones Al ont une vitesse de sédimentation comprise entre 0 et 8 m/Ma, la zone A2 entre 8 et 16 m/Ma, la zone A3 entre 16 et 32 m/Ma, la zone A4 entre 32 et 56 m/Ma, la zone A5 entre 56 et 64 m/Ma, la zone A6 entre 64 et 72 m/Ma. 3045098 18 Les résultats de mesures sur les échantillons de roche par Rock-Eval 6 ® (IFP Energies nouvelles, France) sont reportés dans le tableau 1 ci-dessous avec la position des puits utilisés (les positions I et J représentent les numéros des mailles dans le domaine). 5 Tableau 1 - TOC et IH mesurés Position en I Position en J TOCA(P) IHA(P) 10 58 1.36 0 10 57 1.39 52 12 57 1.49 143 15 48 1.52 21 17 44 1.66 59 21 39 1.02 105 29 39 1.00 189 23 43 0.64 182 A partir de ces informations, au moyen du procédé selon l'invention, les valeurs des paramètres biologiques et chimiques sont déterminées par calage, puis appliquées à l'ensemble du bassin pour obtenir des cartes à l'échelle du bassin des valeurs de TOC° (figure 10 5) et de IH0 (figure 6). Sur les figures 5 et 6, les zones blanches correspondent à une zone continentale, non lacustre ou marine. Par conséquent, en ces zones il n'est pas utile de déterminer les teneurs en carbone organique et l'indice d'hydrogène. Sur la figure 5, la zone B1 a un TOC compris entre 0 et 1,2 %, la zone B2 un TOC entre 1,2 et 2 %, la zone B3 un TOC entre 2 et 2.4 %, la zone B4 un TOC entre 2.4 et 2.8 %, la 15 zone B5 un TOC au-delà de 2.8 `Vo. Sur la figure 6, la zone Cl a un IH supérieur à 520, la zone C2 un IH entre 480 et 520, la zone C3 un IH entre 460 et 480, la zone C4 un IH entre 440 et 460, la zone C5 un IH inférieur à 440. 20 Pour cet exemple comparatif, on compare les valeurs de TOCA(P) déterminés au moyen du procédé selon l'invention, avec des valeurs de TOCO(P) obtenues, sur la base d'une 3045098 19 méthode de détermination du TOCS, à partir des valeurs de IH, IH0 et TOC actuelle, couramment admise (ex : Jarvie, 2012). Pour la comparaison, on utilise la formule suivante de l'art antérieur : TOC° (P) = TOCA(P) /Hos(P) - /HA (B) 1 x 0.84 1000 5 Tableau 2 - TOC simulés et calculés Position en I Position en J T000(P) T000(P) selon l'art antérieur selon l'invention 10 58 2.199239913 2.202300298 10 57 2.1640443 2.079777204 12 57 2.093971256 2.024799477 15 48 2.659069234 2.522195442 17 44 2.540686772 2.758488893 21 39 1.232681147 1.626519423 29 39 1.17093617 0.923460331 23 43 2.199239913 2.202300298 On observe donc que les valeurs de teneur organiques au temps du dépôt sont très similaires au niveau des puits. La construction de la carte de TOC selon l'invention est donc bien représentative du bassin sédimentaire au niveau des puits. 10 REVENDICATIONS1) Procédé d'exploitation d'un bassin sédimentaire, ledit bassin sédimentaire comprenant au moins une couche sédimentaire, et étant traversé par au moins un puits, dans lequel on connaît la valeur de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au temps actuel au niveau d'au moins un puits, une carte de bathymétrie et une carte de vitesse de sédimentation de ladite couche sédimentaire, caractérisé en ce qu'on réalise les étapes suivantes : a) on construit pour au moins une couche sédimentaire des cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire à partir desdites valeurs de teneur en carbone organique TOCA(P) et de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits, et au moyen desdites cartes de bathymétrie et de vitesse de sédimentation de ladite couche sédimentaire ; b) on détermine au moins une zone dudit bassin sédimentaire comportant des hydrocarbures au moyen desdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire ; et c) on exploite ladite zone dudit bassin sédimentaire. 2) Procédé selon la revendication 1, dans lequel on construit pour au moins une couche sédimentaire des cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire en mettant en oeuvre les étapes suivantes : i) on construit des modèles de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 en tout point de ladite couche sédimentaire au temps du dépôt de ladite couche sédimentaire, lesdits modèles reliant ladite teneur en carbone organique TOCS et ledit indice d'hydrogène IH0 en tout point de ladite couche sédimentaire au temps du dépôt de ladite couche sédimentaire, avec lesdites cartes de bathymétrie et de vitesse de sédimentation de ladite couche sédimentaire et avec des paramètres chimiques et biologiques de ladite couche sédimentaire ; ii) on cale lesdits modèles en déterminant lesdits paramètres chimiques et biologiques au moyen desdites valeurs actuelles de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au niveau dudit puits ; et iii) on construit lesdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire en appliquant lesdits modèles calés en toute maille de ladite couche sédimentaire au moyen 3045098 21 des paramètres chimiques et biologiques déterminés et desdites cartes de vitesse de sédimentation et de bathymétrie. 3) Procédé selon la revendication 2, dans lequel on cale lesdits modèles en mettant en 5 oeuvre les étapes suivantes : (1) on choisit des paramètres biologiques et chimiques de ladite couche sédimentaire ; (2) on détermine des valeurs de teneur en carbone organique T000s(P) et d'indice d'hydrogène IHos(P) au niveau dudit puits au temps du dépôt de ladite 10 couche sédimentaire au moyen desdits modèles, desdits paramètres biologiques et chimiques choisis et desdites cartes de bathymétrie et de vitesse de sédimentation ; (3) on calcule une valeur de teneur en carbone organique TOCosc(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire au moyen d'un 15 bilan carbone et desdites valeurs de teneur en carbone organique TOCA(P) et la valeur de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits ; (4) on compare lesdites valeurs déterminées TOCos(P) et calculées TOCoBc(P) de teneur en carbone organique au niveau dudit puits au temps de dépôt de ladite couche sédimentaire, et on réitère les étapes (2) à (4) en modifiant 20 lesdits paramètres biologiques et chimiques pour minimiser ladite comparaison. 4) Procédé selon la revendication 3, dans lequel on calcule la valeur de teneur en carbone organique TOCosc(P) au niveau dudit puits au temps du dépôt de ladite couche 25 sédimentaire au moyen d'une formule du type : TOCA(P) C A(P) TOC Pl = 1 IHos(P)-IHA(B) X0.84 1000 5) Procédé selon la revendication 3, dans lequel on calcule la valeur de teneur en carbone organique TOCosc(P) au niveau dudit puits au temps du dépôt de ladite couche 30 sédimentaire au moyen d'une formule du type : TOCA(P)-0.84xS1 T000Bc(P) = 1 /Hos(P)-/HA(B)x0.84 avec S1 la quantité d'huiles libres. l000 6) Procédé selon l'une des revendications 3 à 5, dans lequel ladite valeur de teneur en carbone organique T000Bc(P) au niveau dudit puits au temps du dépôt de ladite couche 3045098 22 sédimentaire est fonction d'un paramètre relatif à la formation de gaz biogénique au sein de ladite couche sédimentaire. 7) Procédé selon l'une des revendications 2 à 6, dans lequel lesdits paramètres biologiques 5 et chimiques sont des paramètres représentatifs du dépôt de matière organique et/ou des réactions de réduction et/ou d'oxydation. 8) Procédé selon l'une des revendications 2 à 7, dans lequel lesdits modèles sont formés par interpolation desdites valeurs de teneur en carbone organique TOCS (P) et d'indice 10 d'hydrogène IHA(P) au niveau dudit puits au temps du dépôt de ladite couche sédimentaire, l'interpolation étant contrainte par lesdites cartes de vitesse de sédimentation et de bathymétrie. 9) Procédé selon l'une des revendications précédentes, dans lequel lesdites valeurs de 15 teneur en carbone organique TOCA(P) et de l'indice d'hydrogène IHA(P) au temps actuel au niveau dudit puits sont déterminées par analyse d'une pyrolyse d'une roche de ladite couche sédimentaire. 10) Procédé selon l'une des revendications précédentes, dans lequel on détermine au moins 20 une zone dudit bassin sédimentaire comportant des hydrocarbures au moyen d'un modèle dudit bassin sédimentaire construit au moyen desdites cartes de teneur en carbone organique TOCS et d'indice d'hydrogène IH0 au temps du dépôt de ladite couche sédimentaire. 25 11) 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 précédentes, lorsque ledit programme est exécuté sur un ordinateur. 30