Brevet FR3067500A1 — publié le 20181214

La présente invention concerne le domaine de l'exploration et de l'exploitation des gisements pétroliers ou de sites de stockage géologique de gaz. L'exploration pétrolière consiste à rechercher des gisements d'hydrocarbures au sein d'un bassin sédimentaire. La compréhension des principes de la genèse des hydrocarbures et leurs liens avec l'histoire géologique du sous-sol, a permis la mise au point de méthodes d'évaluation du potentiel pétrolier d'un bassin sédimentaire. La démarche générale de l'évaluation du potentiel pétrolier d'un bassin sédimentaire comporte des allers-retours entre une prédiction du potentiel pétrolier du bassin sédimentaire, réalisée à partir d'informations mesurées relatives au bassin étudié (analyse d'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 informations permettant de préciser les prédictions de potentiel pétrolier du bassin étudié. 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 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. L'évaluation du potentiel pétrolier d'un bassin sédimentaire peut se faire soit de manière probabiliste soit de manière analytique. L'approche analytique est le plus souvent mise en oeuvre au moyen de logiciels exécutés sur ordinateur, permettant la synthèse des données disponibles et la simulation de l'histoire géologique du bassin étudié. En fonction de la complexité de l'histoire géologique du bassin étudié, la famille des logiciels mis en oeuvre pour l'approche analytique peut comprendre la simulation à une, deux ou trois dimensions, des processus sédimentaires, tectoniques, thermiques, hydrodynamiques, de chimie organique et inorganique qui interviennent lors de la formation d'un bassin pétrolier. On parle alors d'une démarche dite de « modélisation de bassin ». Concernant plus particulièrement la modélisation des processus sédimentaires, on met en oeuvre des logiciels dits de modélisation stratigraphique. La modélisation stratigraphique a pour objectif de simuler l'évolution d'un bassin sédimentaire au cours des temps géologiques afin de quantifier notamment la géométrie des couches sédimentaires, le type de sédiments déposés, la profondeur d'eau au moment du dépôt, etc. On connaît des modèles stratigraphiques 3D, résolus de manière numérique, via des logiciels exécutés sur ordinateur, tels que le simulateur 2 stratigraphique DIONISOSFLOWTM (IFP Energies nouvelles, France). De manière générale, les différents modèles stratigraphiques prennent en compte trois grands processus : - la création (ou la suppression), d'espace disponible à la sédimentation par les mouvements tectoniques, eustatiques et flexuraux, 5 - l'apport en sédiment dans le bassin, soit par les frontières, soit par l'intermédiaire d'une production ou d'une précipitation in situ, et - le transport de ces sédiments dans l'espace disponible créé, ces transports peuvent être le transport permanent des sédiments, le transport épisodique des sédiments et/ou le transport catastrophique de ces sédiments. 10 Une nouvelle stratigraphie est ensuite déterminée sur la base du principe de conservation de la masse. Etat de la technique 15 Les documents suivants seront cités dans la suite de la description : Bohacs, K and Suter, J. (1997). Sequence Stratigraphic Distribution of Coaly Rocks: Fundamental Controls and Paralic Examples. AAPG Bulletin, V. 81, No. 10 (October 1997), P. 20 1612-1639. Granjeon D. et Joseph P. (1999) Concepts and applications of a 3D multiple lithology, diffusive model in stratigraphic modeling. In: Harbaugh J. W., Watney W. L., Rankey E. C., Slingerland R., Goldstein R. H., Franseen E. K. (Eds.). Numerical Experiments in Stratigraphy: Recent Advances in Stratigraphic and Sedimentologic Computer Simulations. SEPM Special 25 Publications, 62, p. 197-210. Mann, U. and Zweigel, J., (2008). Modelling source rock distribution and quality variations: The OF-Mod approach. In: Analogue and Numerical Forward Modelling of Sedimentary Systems; from Understanding to Prediction (P.L. de Boer, G. Postma, C.J. van der Zwan, P.M. Burgess and P. Kukla, eds.). Special Publication 40 of the International Association of 30 Sedimentologists, pp. 239-274. L'évaluation du potentiel pétrolier d'un bassin sédimentaire nécessite de caractériser au mieux la distribution des sédiments au sein d'un bassin sédimentaire, et plus particulièrement la distribution et la qualité de la matière organique au sein de ces sédiments, contrôlées par les 35 processus de dépôt, de dégradation et/ou de préservation de cette matière organique. Une 3 bonne caractérisation des couches sédimentaires riches en matière organique passe ainsi notamment par une bonne compréhension de la stratigraphie du bassin, c'est-à-dire de la géométrie des couches sédimentaires, de la répartition des sédiments fins et grossiers au sein de ces couches, et de l'impact qu'aura la dynamique de ce remplissage sur la distribution finale 5 et la qualité de la matière organique. Majoritairement, les modèles stratigraphiques actuels permettent de retracer l'histoire des dépôts de sédiments inorganiques au sein d'un bassin. Toutefois, l'industrie pétrolière souhaite en plus connaître si, au cours des temps géologiques, le bassin a connu des conditions favorables au dépôt de matière organique et à sa préservation, et la localisation et la qualité de 10 ces dépôts le cas échéant. Ainsi, il apparaît important de prédire la distribution et la qualité de la matière organique au sein d'un bassin sédimentaire. De manière générale, on peut distinguer trois types de matière organique : - la matière organique d'origine marine, qui résulte de la production in situ de 15 matière organique, initialement produite en zone photique, puis préservée au cours des temps géologiques ; - la matière organique lacustre, qui résulte de la production in situ de matière organique, initialement produite en zone photique ou benthique, puis préservée au cours des temps géologiques ; 20 la matière organique d'origine terrestre, qui résulte de la production in situ de la matière organique à la surface émergée d'un bassin sédimentaire, puis préservée au cours des temps géologiques. On connaît le modèle Of-Mod TM (SINTEF, Norvège) qui permet de prédire la distribution 25 et la qualité de la matière organique déposée en environnement marin dans des séquences stratigraphiques. Ce modèle est notamment décrit dans le document (Mann and Zweigel, 2008). Ce modèle nécessite un modèle numérique 3D du bassin sédimentaire étudié, fourni par un modèle stratigraphique ou un géomodeleur, et quantifie la teneur en carbone organique dans les sédiments en simulant la production et la dégradation de la matière organique. Toutefois, si 30 ce modèle peut être adapté à la prédiction de la distribution et de la qualité de la matière organique marine, la gestion de la matière organique terrestre dans ce modèle est relativement éloignée de la réalité physique. En effet, l'approche décrite dans cette demande repose sur une hypothèse de corrélation entre les apports en sédiments inorganiques et la matière organique terrestre. De plus, ce modèle ne détermine pas de zones de production/préservation de la 35 matière organique terrestre en fonction notamment d'un bilan hydrologique et/ou des circulations de fluide au sein du bassin. 4 On connaît également la demande EP 2743738 (US 2014/0163883) qui concerne la détermination de la distribution et de la qualité de la matière organique dans un bassin sédimentaire à partir d'un couplage entre un modèle stratigraphique et un modèle de 5 production, dégradation et transport de la matière organique. Toutefois, le modèle de production, dégradation et transport décrit dans cette demande est relativement simpliste dans le cas du dépôt de matière organique terrestre car il ne comprend pas une modélisation réaliste de l'état hydrologique du bassin et/ou des circulations de fluide au sein du bassin. En effet, dans le cas de cette approche, le niveau des nappes phréatiques est supposé proportionnel à 10 un flux d'eau à la surface du bassin, ce qui est une hypothèse qui dans de nombreux cas est éloignée de la réalité. Ainsi, par exemple, cette hypothèse néglige : (1) l'impact très important que peut avoir la présence d'un lac sur la nappe phréatique et donc sur la zone vadose ; (2) l'existence des systèmes de rivières déconnectés des nappes. Ainsi, une telle hypothèse pourrait conduire à ne pas reproduire correctement les environnements palustres (zones de 15 production de matière organique terrestre autour des lacs) ou à surestimer la production de matière organique terrestre dans les hauts topographiques (où rivières et nappes peuvent apparaitre souvent déconnectées). 20 L'invention a pour objet une méthode permettant d'estimer de façon la plus déterministe possible la distribution de la matière organique d'origine terrestre dans un bassin sédimentaire. Pour ce faire, la présente invention repose sur un couplage entre un modèle stratigraphique, un modèle hydrologique, et un modèle d'accumulation in situ de la matière organique terrestre plus proche de la réalité physique. En particulier, le modèle d'accumulation in situ de la matière 25 organique terrestre tient compte des conditions hydriques au sein du bassin sédimentaire, celles-ci étant déterminées à partir d'un modèle hydrologique plus représentatif de la réalité physique. Par ailleurs, le couplage entre le modèle d'accumulation de la matière organique terrestre et du modèle stratigraphique, permet que la matière organique déposée à un pas de temps donné soit transportée, via la modélisation stratigraphique, à un pas de temps suivant. 30 Selon une variante de réalisation, la présente invention permet de modéliser, en sus de l'accumulation de la matière organique terrestre préservée à l'endroit de sa zone de production, l'accumulation de la partie non préservée in situ de la matière organique terrestre produite. 35 5 Le procédé selon l'invention Ainsi, la présente invention concerne un procédé pour exploiter des hydrocarbures au 5 sein d'un bassin sédimentaire, ledit bassin résultant notamment du dépôt de sédiments inorganiques et de matière organique terrestre, ledit procédé étant mis en oeuvre à partir de mesures de propriétés relatives auxdits dépôts et à l'hydrologie dudit bassin, au moyen d'un simulateur stratigraphique exécutant sur ordinateur une modélisation stratigraphique pour une succession de pas de temps. Selon l'invention, pour au moins un pas de temps : 10 A. au moyen dudit simulateur et de paramètres dudit simulateur pour ledit pas de temps, lesdits paramètres dudit simulateur étant construits à partir desdites mesures, on détermine une représentation maillée représentative au moins de la distribution desdits sédiments inorganiques au sein dudit bassin pour ledit pas de temps ; B. au moyen d'un modèle hydrologique et de paramètres dudit modèle 15 hydrologique, lesdits paramètres dudit modèle hydrologique étant construits à partir desdites mesures et de ladite représentation maillée, on délimite au moins une zone insaturée au sein dudit bassin ; C. au moyen d'un modèle d'accumulation de ladite matière organique et de paramètres dudit modèle d'accumulation, lesdits paramètres dudit modèle 20 d'accumulation comprenant l'épaisseur de ladite zone insaturée audit pas de temps, la stabilité de ladite épaisseur de ladite zone insaturée au cours dudit pas de temps, et le flux d'eau de surface audit pas de temps, on détermine au moins une zone d'accumulation de ladite matière organique terrestre au sein dudit bassin et une quantité de matière organique terrestre accumulée dans ladite zone d'accumulation. 25 Puis, à partir au moins de ladite quantité de matière organique terrestre accumulée et de ladite représentation maillée, on quantifie le potentiel pétrolier dudit bassin, on définit un schéma d'exploitation dudit bassin, et on exploite ledit bassin en fonction dudit schéma. Selon une réalisation de l'invention, chacune des mailles de ladite représentation maillée 30 déterminée au moyen dudit simulateur stratigraphique peut comprendre au moins une propriété relative au faciès lithologique desdits sédiments inorganiques de ladite maille. Selon une mise en oeuvre de l'invention, lesdits paramètres dudit modèle hydrologique peuvent comprendre la conductivité hydraulique en chaque maille de ladite représentation 35 maillée, ladite conductivité hydraulique en une desdites mailles pouvant être déterminée au 6 moins à partir de ladite propriété relative au faciès lithologique desdits sédiments organiques en ladite maille. Selon une variante de réalisation de l'invention, ledit modèle hydrologique peut délimiter 5 une zone saturée dudit bassin en tenant compte de gradients de pression dans ledit bassin, et on définit alors la base de ladite zone insaturée par le toit de ladite zone saturée. Selon une mise en oeuvre de l'invention, une maille de ladite représentation maillée peut être comprise dans une zone d'accumulation potentielle de ladite matière organique terrestre 10 si. - ladite épaisseur de ladite zone insaturée audit pas de temps est au plus égale à lm ; et - ledit flux d'eau de surface est au plus égal à 10 m3/s. Selon une mise en oeuvre de l'invention, une maille de ladite représentation maillée située 15 dans ladite zone d'accumulation potentielle peut faire partie de ladite zone d'accumulation de ladite matière organique terrestre si ladite épaisseur de ladite zone insaturée dans ladite zone d'accumulation potentielle est invariante au cours dudit pas de temps. Selon une variante de réalisation de l'invention, ladite zone insaturée peut être considérée 20 comme stable si une égalité de la forme suivante est vérifiée : Su = SR + SL où SL est un taux de variations eustatiques, SR est un taux de sédimentation relatif au dépôt desdits sédiments organiques et desdits sédiments inorganiques, et SU est un taux de subsidence. 25 Selon une mise en oeuvre de l'invention, on peut déterminer ladite quantité de matière organique terrestre accumulée à partir d'un taux d'accumulation de ladite matière organique terrestre TAM déterminé selon une formule du type : TAM = Su - SRin - SL 30 où SL est un taux de variations eustatiques, SRin est un taux de sédimentation relatif au dépôt desdits sédiments inorganiques, et SU est un taux de subsidence. 7 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 exécuté 5 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 Figures annexées et décrites ci-après. 10 Présentation succincte des figures - La Figure 1 présente un exemple de représentation maillée issue de la mise en oeuvre 15 d'une simulation stratigraphique pour un bassin sédimentaire donné et pour un pas de temps donné. - La Figure 2 présente une zone saturée et une zone insaturée identifiées au moyen du modèle hydrologique selon l'invention pour le bassin sédimentaire et le pas de temps considérés en Figure 1. 20 Description détaillée du procédé De façon générale, l'un des objets de l'invention concerne un procédé pour exploiter les 25 hydrocarbures présents au sein d'un bassin sédimentaire. De manière générale, un bassin sédimentaire est formé par des dépôts de sédiments successifs dans le temps. Les sédiments peuvent être inorganiques (tels que des dépôts clastiques, carbonatés etc) ou bien organique (dépôt de matière organique d'origine lacustre, marine et/ou terrestre). La présente invention s'applique dans le cas où le bassin sédimentaire étudié résulte du dépôt de sédiments 30 inorganiques et de sédiments organiques au moins d'origine terrestre. En particulier, la présente invention vise à quantifier de manière fiable la distribution de la matière organique d'origine terrestre s'étant déposée au sein d'un bassin sédimentaire. Par matière organique terrestre, on entend une matière organique qui n'est au moins ni d'origine marine ni d'origine lacustre. La présente invention est particulièrement mais non exclusivement appropriée dans le 35 cas du dépôt de matière organique terrestre dans des environnements de type tourbière. 8 Le procédé selon l'invention repose sur au moins trois modèles couplés : un modèle stratigraphique permettant de modéliser l'évolution stratigraphique d'un bassin sédimentaire, un modèle hydrologique délimitant notamment la zone insaturée du bassin étudié, et un modèle d'accumulation de la matière organique d'origine terrestre. Selon une mise en oeuvre de 5 l'invention, les simulations numériques reposant sur chacun de ces modèles sont exécutées sur ordinateur de manière séquentielle, ou bien, alternativement, les modèles stratigraphiques, hydrologiques et d'accumulation selon l'invention sont résolus simultanément. La présente invention requiert de disposer de : 10 - de mesures relatives aux dépôts de sédiments inorganiques du bassin étudié : il s'agit de mesures réalisées in situ, permettant notamment de déterminer l'architecture du bassin à l'actuel, de qualifier les différents dépôts sédimentaires inorganiques (composition minéralogique ou au moins du type de dépôt, épaisseurs, âges, conditions de dépôts etc) ou encore les événements géologiques subis par le bassin au cours des temps géologiques 15 (érosion, subsidence, eustatistme etc). Ces mesures peuvent consister en des études d'affleurements et/ou d'analogue, en des mesures diagraphiques réalisées le long de puits forés dans le bassin étudié, en des échantillons de roche prélevés par exemple par carottage, en des images sismiques obtenues suite à des campagnes d'acquisition sismique. A partir des mesures ainsi réalisées, on peut définir les paramètres d'entrée d'une 20 simulation stratigraphique, tels que la morphologie du bassin (notamment la topographie, mais aussi les limites de couches), les apports en sédiments inorganiques, leur transport pour chaque pas de temps. Par exemple, la morphologie pour un pas temps donné peut être obtenue à partir de la morphologie mesurée à l'actuel (sur une image sismique par exemple) et au moyen d'une restauration structurale, technique bien connue du spécialiste, 25 qui peut par exemple être réalisée au moyen du logiciel KronosFlow (IFP Energies nouvelles, France). De même, les apports en sédiments pour un pas de temps donné peuvent être déterminés à partir du volume de la couche sédimentaire observée à l'actuel, divisé par la durée de la sédimentation, et, le cas échéant, en tenant compte d'un taux de production de matière organique terrestre. Enfin, les paramètres de transport pour un pas 30 de temps donné peuvent être estimés à partir de profils de dépôts établis à partir de l'observation du bassin au temps actuel. Ces mesures peuvent également permettre de définir des paramètres d'entrée du modèle hydrologique et du modèle d'accumulation de matière organique terrestre selon l'invention ; - de mesures relatives au dépôt de matière organique d'origine terrestre : il s'agit de mesures 35 réalisées in situ, par exemple dans des tourbières existants à l'actuel, dans le bassin étudié ou bien dans un analogue, à partir desquelles on peut déduire notamment un taux de production de matière organique terrestre. Par exemple, on peut mesurer in situ un taux de 9 production de tourbières observables actuellement, puis au moyen d'une méthode de mise à l'échelle (méthode d'« upscaling »), on transcrit ce taux de production déterminé sur une échelle de temps courte (1 à 10 mm/an) en un taux à une échelle de temps plus grande (de l'ordre du million d'années). De manière classique, les taux de production de matière 5 organique terrestre observés aux échelles géologiques sont compris entre 1 et 200 m/Ma. On connaît par exemple le document (Bohacs et Suter, 1997) qui reporte une gamme des taux de production de matière organique en fonction des climats. Ces mesures peuvent permettent de définir des paramètres d'entrée du modèle stratigraphique et du modèle d'accumulation de matière organique terrestre selon l'invention ; 10 de mesures relatives à l'hydrologie du bassin : il s'agit de mesures réalisées in situ, permettant notamment de qualifier des paramètres hydrologiques du bassin, à l'actuel et de manière préférée pour des temps géologiques antérieurs. A l'actuel, il peut s'agir de mesures de la pluviométrie, de l'évaporation, des flux des aquifères du bassin étudié. Pour des temps géologiques antérieurs, on peut par exemple extrapoler les mesures actuelles à 15 partir de la connaissance des climats antérieurs déduite de l'observation du bassin à l'actuel, ou on peut supposer que les paramètres de transport sont équivalents à ceux des fleuves actuels, et que la pluviométrie et l'évaporation sont définies en se référant aux zones climatiques actuelles. Ces mesures permettent notamment de définir les paramètres de transport nécessaires en entrée d'une simulation stratigraphique, mais également de 20 définir des paramètres du modèle hydrologique et du modèle d'accumulation de matière organique terrestre selon l'invention ; un simulateur stratigraphique : il s'agit d'un logiciel exécuté sur ordinateur visant à reconstituer les processus sédimentaires ayant affecté le bassin depuis un temps géologique t jusqu'au temps actuel. Une simulation stratigraphique numérique est de 25 manière générale implémentée de manière discrète dans le temps, c'est-à-dire qu'une simulation stratigraphique simule l'état stratigraphique du bassin pour une succession de pas de temps. Un pas de temps d'un simulateur stratigraphique correspond à une durée géologique au cours de laquelle des dépôts ou des érosions sédimentaires ont été enregistrées. La simulation de l'histoire du remplissage d'un bassin sédimentaire est 30 réalisée à partir des paramètres d'entrée représentatifs de l'histoire sédimentaire du bassin étudié. Selon une mise en oeuvre de l'invention, les paramètres d'entrée d'une simulation stratigraphique sont au moins (1) l'espace disponible pour la sédimentation, lié à des mouvements tectoniques et/ou eustatiques et à la compaction mécanique des sédiments (ou tassement des sédiments sous l'effet du poids des couches sus-jacentes) (2) l'apport en 35 sédiments dans le bassin, soit par les frontières, soit par l'intermédiaire d'une production ou d'une précipitation in situ, (3) le transport de ces sédiments (capacité de transport estimée à partir des caractéristiques du sédiment comme la taille des grains ou encore leur densité, du flux d'eau s'écoulant à la surface du sol et de la pente locale du bassin) dans l'espace 10 disponible créé. Le système d'équations décrivant ces processus peut par exemple être résolu par une discrétisation spatiale en volume fini, et un schéma explicite en volumes finis. De manière classique, à l'issue d'une simulation stratigraphique, on dispose au moins de l'évolution de la topographie du bassin au cours du temps, de la distribution des sédiments 5 dans le temps et l'espace, de la distribution des environnements de dépôt (marin, lacustre ou terrestre), de la localisation dans le temps et l'espace de la côte, et optionnellement du flux d'eau de surface (tel que le ruissèlement). Selon une mise en oeuvre de l'invention, la distribution des sédiments au sein du bassin pour un pas de temps donné peut se présenter par une représentation maillée pour laquelle chaque maille est au moins remplie en faciès 10 (par exemple on distingue un grés d'un calcaire, mais on peut également distinguer un calcaire d'un autre calcaire par son pourcentage de marne). De manière analogue, la distribution des environnements de dépôt pour un pas de temps donné peut se présenter sous la forme d'une propriété catégorielle représentative du type d'environnement (marin, lacustre, terrestre) ; les mailles voisines ayant la même valeur de propriété catégorielle 15 appartiennent alors à un même environnement de dépôt. Le modèle stratigraphique selon l'invention ne comprend pas une modélisation de la production de matière organique d'origine terrestre. On peut par exemple trouver une description d'un tel simulateur stratigraphique dans le document (Granjeon et Joseph, 1999). Un exemple d'un tel simulateur stratigraphique est le logiciel DionisosFlowTM (IFP Energies nouvelles, France). 20 Selon une mise en oeuvre de l'invention, le procédé selon l'invention comporte au moins les étapes 1) à 4) suivantes, l'enchainement des étapes 1) à 3) étant appliqué pour au moins un pas de temps de la simulation stratigraphique. Alternativement, les étapes 1) à 3) sont réalisées 25 simultanément, ce qui signifie que les modèles sur lesquels reposent ces étapes sont résolus simultanément. Cette alternative permet une plus grande efficacité en termes de temps de calcul et/ou d'éviter que la solution de l'un des modèles soit fonction de la solution d'un modèle résolu à l'étape précédente. 30 1. Modélisation stratigraphique 2. Modélisation hydrologique de la zone insaturée 3. Modélisation de l'accumulation de la matière organique terrestre 4. Exploitation des hydrocarbures de la formation 11 De manière préférée, les étapes 1) à 3) sont répétées pour chaque pas de temps de la simulation stratigraphique. Avantageusement, les étapes 1) à 3) sont au moins répétées pour chaque pas de temps de la simulation stratigraphique pour lesquels une modélisation des dépôts de matière organique d'origine terrestre est souhaitée. 5 1. Modélisation stratigraphique Il s'agit au cours de cette étape de modéliser les dépôts sédimentaires au moins inorganiques au sein du bassin étudié pour le pas de temps considéré. Cette étape est réalisée au moyen d'un simulateur stratigraphique tel que décrit ci-dessus, à partir de paramètres de 10 simulation stratigraphique tels que décrits ci-dessus, ces paramètres étant déterminés au moins à partir de mesures relatives au dépôt de sédiments inorganiques et à l'hydrologie du bassin étudié telles que décrites ci-dessus. Le simulateur stratigraphique selon l'invention peut toutefois permettre la modélisation des dépôts sédimentaires organiques d'origine lacustre et/ou marine. 15 Selon un mode de réalisation de l'invention, on met en oeuvre cette étape au moyen d'un simulateur stratigraphique tel que celui décrit dans le document (Granjeon et Joseph, 1999). Un exemple d'un tel simulateur stratigraphique est le logiciel DionisosFlowTM (IFP Energies nouvelles, France). Selon l'invention, à l'issue de cette étape, on obtient une représentation maillée 20 représentative au moins de la distribution des sédiments inorganiques au sein dudit bassin pour le pas de temps considéré. Selon une mise en oeuvre de l'invention, chacune des mailles de cette représentation est renseignée au moins en termes de faciès lithologique (grès, carbonates etc) A titre purement illustratif, la Figure 1 présente une section verticale au sein d'une 25 représentation maillée issue d'une simulation stratigraphique pour un pas de temps donné. On peut observer sur cette figure que l'architecture du bassin étudié pour le pas de temps considéré est en forme de cuvette dans laquelle se sont déposés des sédiments successifs, une partie de ces sédiments ayant été apportés dans le bassin par un flux d'eau de surface F, au moins une partie du bassin étant situé sous le niveau de la mer M. On peut également 30 distinguer sur cette figure la limite de la côte C, qui délimite le domaine marin du domaine terrestre. La simulation stratigraphique pour le pas de temps considéré a permis de simuler le dépôt sédimentaire D (couche de mailles la plus en surface du bassin). Il s'agit pour l'exemple illustré d'un dépôt de matière inorganique, déposée dans le domaine terrestre (ou zone continentale). 35 12 2. Modélisation hydrologique de la zone insaturée Au cours de cette étape, il s'agit de délimiter la zone insaturée au sein du bassin étudié pour le pas de temps considéré, au moyen d'un modèle hydrologique et de paramètres d'entrée du modèle hydrologique. 5 De manière générale, un bassin sédimentaire ayant connu un dépôt de matière organique terrestre a, au moins à temps géologique donné, présenté une zone continentale (ou domaine terrestre). On peut notamment déterminer la zone continentale d'un bassin donné à partir des environnements de dépôt tels que déterminés par exemple au moyen d'une simulation stratigraphique. De manière générale, une zone continentale d'un bassin sédimentaire est 10 soumise à un flux de fluide météorique (c'est-à-dire d'origine non marine) et comprend : o une zone insaturée ou zone vadose : il s'agit de la partie superficielle de la zone continentale du bassin sédimentaire, non saturée en eau, formée d'eau douce, soumise à des flux hydrologiques liés à la percolation des précipitations au travers du milieu poreux ; 15 o une zone saturée ou phréatique : située sous la zone vadose, il s'agit d'une zone saturée en eau douce, soumise à des flux hydrologiques liés au flux des nappes phréatiques souterraines. Autrement dit, le toit de la zone phréatique correspond à la base de la zone vadose et réciproquement, et le toit de la zone vadose correspond à la topographie du bassin au pas de 20 temps considéré. Selon une mise en oeuvre de l'invention, on détermine le toit de la zone phréatique au moyen dudit modèle hydrologique, à partir duquel on déduit la base de la zone vadose. De manière générale, les paramètres d'entrée du modèle hydrologique selon l'invention sont déterminés à partir d'au moins une partie des mesures relatives aux dépôts de sédiments 25 et à l'hydrologie du bassin telles que décrites précédemment, et d'au moins une partie des informations présentes dans la représentation maillée obtenue à l'issue de l'étape précédente. Selon une mise en oeuvre de l'invention, le modèle hydrologique selon l'invention peut requérir d'une propriété représentative de l'aptitude aux écoulements du milieu poreux contenu dans chaque maille de la représentation maillée, telle que la perméabilité ou la conductivité 30 hydraulique, pour le pas de temps considéré. Selon une mise en oeuvre de l'invention, à partir de la représentation maillée remplie en faciès telle qu'obtenue à l'issue de l'étape précédente et à partir d'une table de correspondance entre type de faciès et valeur de conductivité hydraulique, on remplit chaque maille de la représentation maillée en conductivité hydraulique. De telles tables de correspondance sont bien connues du spécialiste. De préférence, on peut 35 avoir recours à des tables de correspondance établies sur la base d'échantillons de roches prélevés au sein du bassin étudié. Selon une autre mise en oeuvre de l'invention, la conductivité 13 hydraulique peut être déterminée à partir d'un type de faciès, d'une composition minéralogique et d'une porosité (ces deux dernières propriétés pouvant être obtenues par exemple au moyen d'une simulation stratigraphique ou à partir de tables de correspondance), et au moyen d'une méthode d'homogénéisation telle que décrite dans la demande de brevet EP 3104199. 5 Selon une mise en oeuvre de l'invention, en plus des conductivités hydrauliques, le modèle hydrologique selon l'invention peut en outre requérir des propriétés telles que : - un taux de précipitations : il s'agit d'un pourcentage représentatif des précipitations survenues au pas de temps considéré. Ce taux de précipitations est prédéfini par le spécialiste en fonction de ses connaissances sur le bassin au pas de temps considéré 10 (par exemple à partir du climat moyen au pas de temps considéré), et notamment à partir de mesures relatives à l'hydrologie et/ou aux dépôts sédimentaires dans le bassin sédimentaire étudié telles que décrites ci-dessus ; - un coefficient d'infiltration, pondérant l'apport d'eau de la surface (lié aux précipitations) vers le système de nappe phréatique : en effet, seule une partie du flux d'eau lié aux 15 précipitations de surface percole jusqu'à la nappe phréatique. Ce coefficient d'infiltration est par exemple compris entre 1 et 100 (3/0, et est prédéfini par le spécialiste en fonction de ses connaissances sur le bassin au pas de temps considéré, et notamment à partir de mesures relatives à l'hydrologie et/ou aux dépôts sédimentaires dans le bassin étudié telles que décrites ci-dessus ; 20 - un flux d'eau de surface : il s'agit du flux d'eau se déplaçant le long de la topographie du bassin sédimentaire, lié par exemple aux ruissèlements ou aux cours d'eau présents à la surface du bassin au pas de temps considéré. Le flux d'eau de surface est par exemple compris entre 10 et 1000 m3/s. Il peut être prédéfini par le spécialiste en fonction de ses connaissances sur le bassin au pas de temps considéré. Par exemple, à 25 partir d'une image sismique, le spécialiste peut détecter la présence de chenaux dans une couche géologique, et à partir notamment de la géométrie de ces chenaux, définir un ordre de grandeur du flux d'eau au moment où cette couche était à la surface du bassin. Le flux d'eau de surface au pas de temps considéré peut aussi être issu de la simulation stratigraphique exécuté pour le pas de temps considéré à l'étape précédente. 30 Le simulateur stratigraphique DionisosFlowTM (IFP Energies nouvelles, France) permet notamment une estimation du flux d'eau de surface à chaque pas de temps ; - la localisation de la côte du bassin sédimentaire : il s'agit de la limite entre le domaine marin et le domaine terrestre, qui se détermine par exemple à partir des environnements de dépôt déterminés à l'issue de la simulation stratigraphique, et de manière générale à 35 partir de la topographie du bassin, de l'eustatisme et de la subsidence au pas de temps étudié. 14 - la topographie du bassin sédimentaire : il s'agit de la topographie en entrée de la simulation stratigraphique pour le pas de temps considéré. La topographie du bassin sédimentaire évolue au cours du pas de temps, en raison du dépôt de sédiments pour le pas de temps courant et de la subsidence subie par le bassin sédimentaire au cours de 5 ce même pas de temps. De manière préférée, on peut délimiter au moins le toit de la zone saturée (et donc la base de la zone insaturée) à l'aide d'une équation hydrologique dans laquelle les écoulements dans les nappes phréatiques libres sont conditionnés par des gradients de pression dans le 10 sous-sol. Plus précisément, on peut définir un modèle hydrologique reposant sur une équation du type : ah = avec : h - la charge hydraulique [L], 15 - le temps [T], - le tenseur symétrique de conductivité hydraulique (ou perméabilité) [L.T-1], - le coefficient d'emmagasinement spécifique [L-1], - un terme source d'extraction (pompage) [V]. 20 Le terme source Q d'extraction correspond à la recharge en eau (ou au pompage) de la nappe par les précipitations diminuées de l'évapotranspiration, ainsi qu'aux apports en eau aux frontières de la zone simulée (apport des zones amont et export vers les zones aval). Selon cette mise en oeuvre de l'invention, le substratum du bassin sédimentaire est supposé imperméable. La résolution de cette équation permet de déterminer les propriétés hydrauliques 25 (charge hydraulique et vitesse d'écoulement) en tout point du bassin. L'élévation du toit de la nappe phréatique, qui sépare la zone phréatique de la zone vadose, et les gradients de pressions sont déduits des charges calculées. Une telle réalisation du modèle hydrologique permet notamment la détermination du niveau des nappes phréatiques à l'état stationnaire. On peut alors en déduire la localisation de la base de la zone vadose, celle-ci correspondant au toit 30 de la zone phréatique. 15 A titre purement illustratif, la Figure 2 présente un exemple de mise en oeuvre d'un modèle hydrologique pour le même bassin et pour le même pas de temps que celui considéré à la Figure 1. Ainsi le bassin étudié, pour ce pas de temps, comporte une zone vadose météorique ZV et une zone phréatique météorique ZP, le toit de la zone phréatique 5 correspondant à la base de la zone vadose. Sur cette Figure 2 sont également présentés, à titre illustratif, sous la forme de flèches, les flux hydrologiques subis par le bassin étudié : un flux de direction verticale dans la zone vadose ZV, lié à la percolation des précipitations, un flux de direction latérale dans la zone phréatique ZP, lié au déplacement du flux d'eau souterrain sous l'effet des gradients de pression dans le bassin, ainsi que le flux d'eau de surface F lié à la 10 présence d'un fleuve à la surface du bassin. Cette figure présente également l'emplacement de deux colonnes de mailles COL1 et COL2 qui traversent à la fois la zone vadose ZV et la zone phréatique ZP. On peut notamment observer que l'épaisseur de la zone vadose est différente entre ces deux colonnes. 15 3. Localisation et quantification de l'accumulation de la matière organique terrestre Au cours de cette étape, il s'agit de localiser et quantifier l'accumulation de la matière organique d'origine terrestre au cours du pas de temps considéré. Selon un mode de réalisation de l'invention, cette étape peut être mise en oeuvre selon trois sous-étapes : 20 - à partir de la délimitation de la zone insaturée telle que déterminée à l'étape précédente, on détermine des zones d'accumulation potentielles de la matière organique terrestre ; - parmi les zones d'accumulation potentielles de la matière organique terrestre, on détermine les zones d'accumulation effectives de la matière organique terrestre ; - on détermine un taux d'accumulation de la matière organique terrestre dans les zones 25 d'accumulation effectives et on en déduit la quantité de matière organique accumulée. 3.1. Détermination des zones d'accumulation potentielles de la matière organique terrestre Au cours de cette sous-étape, on détermine des zones potentielles d'accumulation de la matière organique terrestre en considérant deux critères basés sur : 30 l'épaisseur de la zone insaturée : en effet, lorsque la zone vadose est de faible épaisseur (par exemple inférieure ou égale à 1m), on peut considérer qu'il s'agit d'une zone « humide », dans laquelle la production de la matière organique terrestre est possible et est même préservée par la présence d'eau qui réduit l'oxydation de la matière organique. A l'opposé, lorsque la zone vadose est épaisse (par exemple 35 supérieure à 1 m), on considère qu'il s'agit d'une région « sèche » non propice à la 16 production de matière organique d'origine terrestre. Selon l'invention, l'épaisseur de la zone insaturée est directement déterminée à partir de la délimitation de la zone insaturée déterminée au moyen du modèle hydrologique selon l'invention, et tel que décrit à l'étape précédente ; et 5 - le flux d'eau de surface (lié au ruissèlement, aux cours d'eau de surface etc) : en effet, si le flux d'eau est trop important (par exemple au-delà de 10 m3/s), aucun dépôt de matière organique terrestre n'est réalisé, l'environnement étant considéré comme trop perturbé hydrologiquement pour permettre une vie biologique suffisante à la production d'une quantité non négligeable de matière organique terrestre. Le flux d'eau de surface 10 peut être un paramètre issu de la simulation stratigraphique ou bien être un paramètre prédéfini par le spécialiste. Selon l'invention, lorsque ces deux critères sont satisfaits pour une maille de la représentation maillée, cette maille correspond à une zone d'accumulation potentielle de la matière organique terrestre. 15 3.2. Détermination des zones d'accumulation effectives de la matière organique terrestre Au cours de cette sous-étape, on détermine si les zones potentielles d'accumulation de la matière organique terrestre sont des zones d'accumulation effectives de la matière organique terrestre en considérant un critère basé sur la stabilité de l'épaisseur de la zone insaturée au 20 cours du pas de temps considéré. En effet, selon l'invention, on considère que si, pour une zone s'accumulation potentielle, l'épaisseur de la zone insaturée : - diminue au cours du pas de temps, alors l'environnement du bassin sédimentaire étudié est propice à des inondations régulières, ce qui a pour conséquence d'empêcher toute production de matière organique terrestre ; 25 - augmente au cours du pas de temps, alors l'environnement du bassin sédimentaire étudié s'assèche, ce qui inhibe la production de matière organique et surtout augmente fortement le taux de dégradation de cette matière organique ; - reste stable au cours du pas de temps, alors les conditions hydriques au sein du bassin sédimentaire permettent la production et la préservation de la matière organique durant 30 le pas de temps considéré. Selon l'invention, on détermine un taux de stabilité de la zone vadose au cours du pas de temps à partir de la variation de la topographie et du taux de sédimentation dans le bassin étudié au cours de ce même pas de temps. Selon une variante de réalisation de l'invention, on évalue un taux de stabilité de la zone vadose durant un pas de temps en fonction d'un taux de variations 35 eustatiques (SL), d'un taux de sédimentation relatif au dépôt de sédiments organiques et 17 inorganiques (SR), et d'un taux de subsidence (Su). Par exemple, selon une mise en oeuvre de l'invention, on considère que la zone vadose est stable si l'égalité suivante est vérifiée : Su = SR + SL Autrement dit, selon ce mode de mise en oeuvre de l'invention, on considère que si le 5 changement topographique associé à la subsidence est compensé par le dépôt de sédiments (organiques et inorganiques) et par les variations eustatiques, alors la zone vadose est stable. Dans le cas où cette égalité n'est pas vérifiée, cela signifie soit que l'épaisseur de la zone vadose augmente, créant ainsi une zone d'inondation non favorable à une production de matière organique, soit que l'épaisseur de la zone vadose diminue, créant une zone 10 d'asséchement, également non favorable à une production de matière organique. Le taux de sédimentation de matière inorganique fait partie des paramètres d'entrée de la simulation stratigraphique et peut être défini à partir des mesures relatives à la stratigraphie telles que définies ci-dessus. Le taux de sédimentation organique peut provenir de mesures relatives aux sédiments organiques telles que décrites ci-dessus, ou bien provenir de documents référençant 15 de tels taux, comme par exemple le document (Bohacs et Suter, 1997). 3.3. Quantification de la matière organique terrestre déposée Au cours de cette sous-étape, il s'agit de déterminer la quantité de matière organique accumulée dans les zones d'accumulation effectives telles que déterminées à la sous-étape 3.2 20 décrites précédemment. Selon l'invention, on détermine cette quantité en fonction de l'estimation d'un espace disponible dans chacune des mailles desdites zones d'accumulation effectives. Selon une mise en oeuvre de l'invention, on détermine dans un premier temps un taux d'accumulation de matière organique terrestre au cours du pas de temps à partir de la variation 25 de la topographie, du taux de sédimentation organique terrestre et du taux de sédimentation inorganique dans le bassin étudié au cours de ce même pas de temps. En particulier, selon une variante de réalisation de l'invention, on évalue le taux d'accumulation de matière organique (TAM) durant un pas de temps en fonction des variations eustatiques (SL), de la sédimentation associée au dépôt de sédiments inorganiques (SRin) et organiques (X), et de la subsidence 30 (Su). Le taux de sédimentation de matière organique X est par exemple estimé en se référant au document (Bohacs et Suter, 1997), qui définit des taux de sédimentation organique compris entre 1 et 180 m/Ma en fonction des climats. Par exemple, selon une mise en oeuvre de l'invention, si Su - SRin - SL X, alors le taux d'accumulation de la matière organique TAM est égale à TAM = Su - SRin - SL. Puis, dans un deuxième temps, une fois un taux d'accumulation 35 de matière organique terrestre au cours du pas de temps déterminé, on détermine la quantité de matière organique terrestre accumulée dans chacune des mailles des zones d'accumulation 18 effectives précédemment identifiées, en multipliant le taux d'accumulation de matière organique terrestre ainsi déterminé par la durée du pas de temps. Selon une mise en oeuvre de l'invention au cours de laquelle les étapes 1 à 3 sont 5 répétées pour une succession de pas de temps de la simulation stratigraphique, la quantité de matière organique terrestre accumulée dans une maille donnée à un pas de temps donné peut être transportée vers d'autres mailles de la représentation maillée par l'exécution de la modélisation stratigraphique au pas de temps suivant. 10 Selon une variante de mise en oeuvre de l'invention, lorsque la matière organique terrestre ne peut être totalement accumulée dans la zone de production initiale, à cause d'un espace disponible insuffisant, celle-ci peut être exportée vers une autre zone de dépôt (par exemple un environnement lacustre ou marin). Selon cette variante, à l'issue des sous-étapes 3.1, 3.2 et 3.3, on détermine une quantité de matière organique exportée en faisant la 15 différence entre la quantité de matière organique terrestre produite (égale au taux de production de matière organique retenu pour la mise en oeuvre de l'étape 3, multiplié par la durée du pas de temps considéré) et la quantité de matière organique terrestre accumulée in situ (égale au taux d'accumulation de matière organique terrestre déterminé à l'issue de l'étape 3.3, multiplié par la durée du pas de temps). La quantité de matière organique exportée ainsi déterminée 20 peut, au pas de temps suivant, être transportée vers une autre zone de dépôt du bassin sédimentaire étudié, via la modélisation du transport de sédiments de la simulation stratigraphique selon l'invention. 25 4- Exploitation des hydrocarbures du bassin sédimentaire A partir de la localisation et de la quantification de la matière organique terrestre déposée au sein du bassin sédimentaire pour au moins un pas de temps et de manière préférentielle pour tous les pas de temps de la simulation, on estime le potentiel en hydrocarbures du bassin sédimentaire étudié. 30 En effet, les hydrocarbures sont issus de la matière organique présente dans le bassin. Ainsi, le potentiel en hydrocarbures peut notamment être estimé à partir de la localisation de la matière organique dans le bassin, de son abondance (on parle de richesse en matière organique), et aussi de son origine (marine, terrestre, lacustre). La présente invention permet notamment de quantifier de manière plus fiable le potentiel en hydrocarbures relatif à de la 35 matière organique d'origine terrestre. 19 Pour cette quantification, on peut par exemple déterminer le TOC (la teneur en carbone organique), les indices d'hydrogène et oxygène, et la minéralogie de la matière organique au sein des couches du bassin sédimentaire. Par exemple, l'indice d'hydrogène final dans les matières organiques préservées est 5 calculé en fonction des proportions de chaque type de matière organique présent. Selon une mise en oeuvre de l'invention, on peut en outre déterminer un indice d'hydrogène propre à la matière organique terrestre et dont la valeur variable dans le temps et dans l'espace, dépend des conditions de préservation de cette matière organique. Ainsi, on peut par exemple définir la qualité de la matière organique à travers une loi du type, HI = HIO. a 1, où HIO est le qualité 10 de référence de la matière organique terrestre, et a et b deux paramètres de calibration permettant de décrire l'impact de la zone vadose sur la qualité de la matière organique terrestre produite dans ces conditions. Ainsi, cette qualité de matière organique terrestre est directement liée aux conditions hydriques de la zone. De telles propriétés peuvent être par exemple utilisées en entrée d'une simulation de 15 bassin (un exemple d'une telle simulation de bassin est le logiciel TemisFlowTm (IFP Energies nouvelles, France) qui permet notamment de quantifier le potentiel pétrolier d'un bassin et de déterminer l'emplacement du réservoir pétrolier au sein de ce bassin. Puis, on définit ensuite un schéma d'exploitation du réservoir identifié au sein de ce 20 bassin. Par exemple, la détermination d'un schéma d'exploitation comprend la définition d'un nombre, d'une géométrie et d'une implantation (position et espacement) de puits injecteurs et producteurs, la détermination d'un type de récupération assistée (par injection d'eau, de tensioactifs, etc.), etc. Un schéma d'exploitation d'un réservoir d'hydrocarbures doit par exemple permettre un fort taux de récupération des hydrocarbures piégés dans le réservoir géologique 25 identifié, sur une longue durée d'exploitation, nécessitant un nombre de puits et/ou d'infrastructures limité. De manière classique, la détermination d'un schéma d'exploitation des hydrocarbures est réalisée à l'aide d'une simulation d'écoulement. Un exemple d'une telle simulation d'écoulement est le logiciel PumaFlowe (IFP Energies nouvelles, France). Puis, une fois un schéma d'exploitation défini, les hydrocarbures piégés dans le réservoir 30 sont exploités en fonction de ce schéma d'exploitation, notamment en forant les puits injecteurs et producteurs du schéma d'exploitation ainsi déterminé, et en installant les infrastructures de production nécessaires au développement du gisement. 35 20 Ainsi, la présente invention repose sur un couplage entre une simulation stratigraphique, un modèle hydrologique et un modèle d'accumulation de la matière organique terrestre permettant la prise en compte de conditions hydriques favorables au dépôt de cette matière organique terrestre et à sa préservation. Plus précisément, la présente invention permet une 5 modélisation de l'évolution des flux hydriques proches de la surface d'un bassin sédimentaire, notamment l'évolution de la zone insaturée de ce bassin, ce qui permet de modéliser de façon plus réaliste la localisation des zones d'accumulation de la matière organique terrestre ainsi que la quantité de matière organique terrestre ainsi accumulée. De cette manière, la présente invention contribue à une meilleure connaissance du potentiel pétrolier du bassin sédimentaire 10 étudié, et permet ainsi de définir des conditions permettant une exploitation optimale de ce bassin. Produit programme d'ordinateur 15 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, REVENDICATIONS1) Procédé pour exploiter des hydrocarbures au sein d'un bassin sédimentaire, ledit bassin résultant notamment du dépôt de sédiments inorganiques et de matière organique terrestre, ledit procédé étant mis en oeuvre à partir de mesures de propriétés relatives auxdits dépôts et à l'hydrologie dudit bassin, au moyen d'un simulateur stratigraphique exécutant sur ordinateur une modélisation stratigraphique pour une succession de pas de temps, caractérisé en ce que, pour au moins un pas de temps : A. au moyen dudit simulateur et de paramètres dudit simulateur pour ledit pas de temps, lesdits paramètres dudit simulateur étant construits à partir desdites mesures, on détermine une représentation maillée représentative au moins de la distribution desdits sédiments inorganiques au sein dudit bassin pour ledit pas de temps ; B. au moyen d'un modèle hydrologique et de paramètres dudit modèle hydrologique, lesdits paramètres dudit modèle hydrologique étant construits à partir desdites mesures et de ladite représentation maillée, on délimite au moins une zone insaturée au sein dudit bassin ; C. au moyen d'un modèle d'accumulation de ladite matière organique et de paramètres dudit modèle d'accumulation, lesdits paramètres dudit modèle d'accumulation comprenant l'épaisseur de ladite zone insaturée audit pas de temps, la stabilité de ladite épaisseur de ladite zone insaturée au cours dudit pas de temps, et le flux d'eau de surface audit pas de temps, on détermine au moins une zone d'accumulation de ladite matière organique terrestre au sein dudit bassin et une quantité de matière organique terrestre accumulée dans ladite zone d'accumulation ; et en ce que, à partir au moins de ladite quantité de matière organique terrestre accumulée et de ladite représentation maillée, on quantifie le potentiel pétrolier dudit bassin, on définit un schéma d'exploitation dudit bassin, et on exploite ledit bassin en fonction dudit schéma. 2) Procédé selon la revendication 1, dans lequel chacune des mailles de ladite représentation maillée déterminée au moyen dudit simulateur stratigraphique comprend au moins une propriété relative au faciès lithologique desdits sédiments inorganiques de ladite maille. 3) Procédé selon la revendication 2, dans lequel lesdits paramètres dudit modèle hydrologique comprennent la conductivité hydraulique en chaque maille de ladite représentation maillée, 22 ladite conductivité hydraulique en une desdites mailles étant déterminée au moins à partir de ladite propriété relative au faciès lithologique desdits sédiments organiques en ladite maille. 5 4) Procédé selon l'une des revendications précédentes, dans lequel ledit modèle hydrologique délimite une zone saturée dudit bassin en tenant compte de gradients de pression dans ledit bassin, et définit la base de ladite zone insaturée par le toit de ladite zone saturée. 5) Procédé selon l'une des revendications précédentes, dans lequel une maille de ladite 10 représentation maillée est comprise dans une zone d'accumulation potentielle de ladite matière organique terrestre si : - ladite épaisseur de ladite zone insaturée audit pas de temps est au plus égale à lm ; et - ledit flux d'eau de surface est au plus égal à 10 m3/s. 15 6) Procédé selon la revendication 5, dans lequel une maille de ladite représentation maillée située dans ladite zone d'accumulation potentielle fait partie de ladite zone d'accumulation de ladite matière organique terrestre si ladite épaisseur de ladite zone insaturée dans ladite zone d'accumulation potentielle est invariante au cours dudit pas de temps. 20 7) Procédé selon la revendication 6, dans lequel ladite zone insaturée est stable si une égalité de la forme suivante est vérifiée : Su = SR + SL où SL est un taux de variations eustatiques, SR est un taux de sédimentation relatif au dépôt desdits sédiments organiques et desdits sédiments inorganiques, et SU est un taux de 25 subsidence. 8) Procédé selon l'une des revendications précédentes, dans lequel on détermine ladite quantité de matière organique terrestre accumulée à partir d'un taux d'accumulation de ladite matière organique terrestre TAM déterminé selon une formule du type : 30 TAM = Su - SRin - SL où SL est un taux de variations eustatiques, SRin est un taux de sédimentation relatif au dépôt desdits sédiments inorganiques, et SU est un taux de subsidence. 23 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 oeuvre du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un 5 ordinateur.