Brevet FR3162529A1 — publié le 20251128

La présente invention concerne un procédé pour le calcul de l’effet de couverture dans le modelage de systèmes de pétrole, comprenant : l’obtention de données lithosphériques et sédimentaires d’entrée ; la discrétisation du modèle différemment dans le domaine sédimentaire et dans le domaine lithosphérique ; le démarrage du calcul de l’effet de couverture ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire est supérieur au temps présent et si l’âge du bassin est supérieure au temps de déposition correspondant à chaque couche sédimentaire, en s’abstenant de réaliser le procédé ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire est inférieur ou égal au temps présent et si l’âge du bassin est inférieure ou égale au temps de déposition correspondant à chaque couche sédimentaire et la vérification de si elle se trouve dans la période d’effondrement, dans lequel : si elle se trouve dans la période d’effondrement, le calcul de l’effondrement ; si elle ne se trouve pas dans la période d’effondrement ou si elle se trouve après l’étape de calcul de l’effondrement, le calcul de la sédimentation ; dans la période d’effondrement, le calcul de la vitesse de déplacement vertical advective à cause de l’amincissement lithosphérique pour chaque coordonnée XY ; le calcul du flux de fluide ; l’utilisation d’éléments finis pour réaliser des calculs à chaque instant dans le temps ; la réalisation de l’addition de l’effet thermique des domaines lithosphérique et sédimentaire à chaque étape de temps. PROCÉDÉ POUR LE CALCUL DE L’EFFET DE COUVERTURE DANS LA MODÉLISATION DE SYSTÈMES DE PÉTROLE ET SUPPORT DE STOCKAGE LISIBLE PAR ORDINATEUR Domaine de l’invention La présente invention concerne le domaine technique de la modélisation de systèmes de pétrole. La présente invention concerne notamment un procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole et support de stockage lisible par ordinateur. Backgrounds de l’invention Dans l’étude de bassins sédimentaires et de systèmes de pétrole, dans le but de mieux évaluer les risques géologies dans l’exploration de nouvelles zones, il est de plus en plus nécessaire d’obtenir de résultats prédictifs plus précis de température et de flux thermique. Ces résultats prédictifs permettent l’analyse de nouvelles zones potentielles, spécialement dans des régions éloignées, où des données obtenues à partir de puits deviennent de plus en plus rares. Afin d’obtenir des résultats prédictifs plus précis pour la température et le flux thermique, des simulations de processus géologiques simulations sont souvent réalisés, où le backstripping est calculé, ce qui consiste à l’enlèvement de couches sédimentaires et leur décompactage pour chaque événement géologique, en obtenant des estimations paléogéométriques pendant l’évolution géologique de la section sédimentaire. Normalement, un logiciel d’analyse de systèmes de pétrole en 3D calcule le backstripping en utilisant plusieurs équations à 1D, en considérant des bassins géologiques sous l’effet de l’isostasie locale. Dans d’autres applications de logiciel, le modèle implémenté calcule le backstripping d’une façon plus raffinée, par couplage des effets d’isostasie de flexion dans le bassin en 3D et en s’approchant alors des conditions observées dans la nature dues aux effets de flexion de la couche lithosphérique. La réponse de flexion permet une détermination plus précise des paléogeométries du socle et des couches sédimentaires et est directement liée à l’état thermique de la lithosphère, qui, à son tour, contrôle le flux thermique dans le socle du bassin sédimentaire à chaque étape de temps. Pour déterminer le flux thermique à utiliser comme condition de limite inférieur dans des simulations à l’intérieur du bassin sédimentaire, il est nécessaire de quantifier l’évolution thermique résultant du gradient géothermique, de l’étirement lithosphérique (effondrement) et de la chaleur radiogénique présente dans la portion supérieure de la croûte. Néanmoins, l’effet thermique transitoire lié à la déposition et au compactage de sédiments froids, connu sous le nom d’effet de couverture, n’est pas mesuré actuellement, effet de couverture qui peut agir pendant et après la sédimentation, provoquant des altérations telles que : réduction du gradient géothermique de la croûte et, par conséquent, réduction du flux thermique dans le socle et chauffage à la fin de chaque phase de déposition par le processus de diffusion de la chaleur transitoire, pour le rééquilibrage thermique du champ de températures de la lithosphère, ce qui dépend de la conductivité thermique, de l’épaisseur et du taux de sédimentation du arrangement sédimentaire. État de la technique Le document intitulé « Impacto do efeito blanketing na história térmica de bacias sedimentares: modelos sintéticos e estudo de caso na Bacia de Santos » (impact de l’effet de couverture sur l’histoire thermique de bassins sédimentaires : des modèles synthétiques et étude de cas dans le bassin de Santos), par Cleriston Ferreira Silva (accessible dans https://pantheon.ufrj.br/handle/11422/13673?mode=full) enseigne un étude approfondi sur l’effet de la déposition de sédiments froids (couverture) sur les flux thermiques dans la partie supérieure du socle, tenant compte ou pas de l’effet de effondrement (amincissement lithosphérique). Spécifiquement, ce document ne décrit que des tests spécifiques et une portée unidimensionnelle (1D) pour l’effet de couverture. Le document « Thermal evolution of the intracratonic Paris Basin: insights from 3D basin modelling » par Martina Torelli, Renaud Traby, Vanessa Teles et Mathieu Ducros (disponible dans https://www.sciencedirect.com/science/article/abs/pii/S0264817220302701?via%3Dihub) englobe un modèle en 3D d’un bassin, en identifiant des paramètres qui touchent à la distribution de la température au cours du temps, montrant l’évolution dans le temps de bassins sédimentaires. Néanmoins, ce document n’enseigne pas l’application de l’effet de couverture. Brève description de l’invention La présente invention, selon un mode de réalisation préféré de celle-ci, définie un procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole comprenant : l’obtention de données lithosphériques et sédimentaires d’entrée ; la discrétisation du modèle différemment dans le domaine sédimentaire et dans le domaine lithosphérique ; le démarrage du calcul de l’effet de couverture ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire est supérieur au temps présent et si l’âge du bassin est supérieure au temps de déposition correspondant à chaque couche sédimentaire, en s’abstenant de réaliser le procédé ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire est inférieur ou égal au temps présent et si l’âge du bassin est inférieure ou égale au temps de déposition correspondant à chaque couche sédimentaire et la vérification de si elle se trouve dans la période d’effondrement, dans lequel : si elle se trouve dans la période d’effondrement, le calcul de l’effondrement ; si elle ne se trouve pas dans la période d’effondrement ou si elle se trouve après l’étape de calcul de l’effondrement, le calcul de la sédimentation ; le calcul de la vitesse de déplacement vertical advective à cause de l’effondrement pour chaque coordonnée XY ; le calcul du flux de fluide ; l’utilisation d’éléments finis pour réaliser l calcul à chaque instant dans le temps ; la réalisation de l’addition de l’effet thermique des domaines lithosphérique et sédimentaire à chaque étape de temps. Spécifiquement, selon un autre mode de réalisation préféré de la présente invention, les données lithosphériques et sédimentaires comprennent au moins l’un de : dimensions de la zone ; paramètres thérmophysiques, tels que : conductivité de la croûte (W/m/K) ; masse volumique de la croûte (kg/m³) ; chaleur spécifique de la croûte (J/kg/K) ; chaleur radiogénique de la croûte (μW/m³) ; conductivité du manteau lithosphérique (W/m/K) ; masse volumique du manteau lithosphérique (kg/m³) ; chaleur spécifique du manteau lithosphérique (J/kg/K) ; manteau lithosphérique chaleur radiogénique (μW/m³) ; masse volumique de l’eau (kg/m³) ; conductivité de l’eau (W/m/K) ; chaleur spécifique de l’eau (J/kg/K) ; masse volumique de l’asthénosphère (kg/m³) ; coefficient d’élasticité (N/m²) ; coefficient de Poisson ; coefficient de dilatation thermique (1/°C) ; accélération de la gravité (m/s²) ; conditions de limite dans la lithosphère, telles que la température de surface et la température de base ; données lithosphériques, correspondant au nombre de couches et aux propriétés thérmophysiques de la couche, comprenant au moins l’une de : épaisseur (m) ; masse volumique (kg/m³) ; chaleur spécifique (J/kg/K) ; conductivité (W/m/K) ; chaleur radiogénique (W/m³*1e-6) ; cartes d’horizons et lithofaciès pour chaque couche sédimentaire ; dans lequel les cartes de salifères et ignées sont optionnelles, en fonction de la région d’intérêt ; valeurs pour la discrétisation du modèle numérique, dans lequel les valeurs peuvent être définies par l’utilisateur, dans la colonne de « nombre d’étapes » ou « séparation », pour la discrétisation temporelle ; et dans les colonnes « nombre d’éléments » ou « séparation », pour la discrétisation spatiale dans les dimensions x, y et z. Une discrétisation plus importante est suggérée : 1) dans la période d’effondrement (des milliers d’années) et 2) dans la partie sédimentaire (des dizaines ou des centaines de mètres dans l’épaisseur en z, correspondant à la profondeur maximale du bassin) ; pour les autres périodes et zones, une discrétisation d’entre 2-3 mille mètres garantira la convergence et des bons résultats) ; âge du bassin (en millions d’années – mya.) ; temps et période de l’effondrement (en millions d’années – mya.) ; valeurs ou cartes des facteurs d’étirement de la lithosphère (β) ou de la croûte et du manteau (δ et β). En outre, selon un autre mode de réalisation préféré de la présente invention, l’étape de calcul de l’effondrement, si elle se trouve dans la période d’effondrement, comprend le calcul de la chaleur transitoire à cause de l’amincissement lithosphérique, ce qui provoque des variations dans le flux thermique, à partir du modèle thermique advectif : où : v s désigne la vitesse des particules solides ; θ est la porosité ; ρ est la masse volumique ; c est la constante de décroissance ; et λ est la conductivité ; le calcul de la valeur de la croûte δ et la valeur du manteau β au sein de l’effondrement, selon : le calcul des épaisseurs et horizons des couches de sédiment en fonction de la valeur de la croûte δ et la valeur du manteau β ; la quantification de la chaleur radiogénique dans la portion supérieure de la croûte. En outre, selon un autre mode de réalisation préféré de la présente invention, l’étape de calcul de la sédimentation comporte le calcul de l’effet de couverture à cause de la déposition de sédiments, par le biais d’un modèle advectif, comportant : le décompactage des couches de sédiment (backstripping) ; l’empilement et le compactage de couches de sédiment jusqu’au temps de déposition correspondant et selon les paramètres thermiques calculés à nouveau en profondeur ; le nouveau calcul de la profondeur du socle ; le calcul de la chaleur radiogénique dans la portion sédimentaire. En outre, selon un autre mode de réalisation préféré de la présente invention, l’étape de calcul de la vitesse de déplacement vertical advective à cause de l’effondrement pour chaque coordonnée XY correspond au calcul de vbulk = G*(h-z) pour chaque coordonnée z, où vbulk = vitesse moyenne du soulèvement de l’asthénosphère à cause de l’effondrement, variant linéairement en z ; G = magnitude du gradient de vitesse verticale le long de la lithosphère h, indiquant combien se déplace verticalement la vitesse de fluide (asthénosphère), le long de l’épaisseur définie comme « h » dans la formulation ; h = limite de la profondeur de la lithosphère. Une fois réalisé ce calcul, on réalise l’addition des vitesses advectives à la vitesse de sédimentation, dans lequel plus le compactage est important, plus le soulèvement à cause de l’effondrement est lent. En outre, selon un autre mode de réalisation préféré de la présente invention, l’étape de calcul du flux de fluide comporte la définition de propriétés thermiques dans les nœuds du maillage numérique en 3D, comprenant au moins l’une de conductivité, chaleur radiogénique, masse volumique*chaleur spécifique de la partie solide (hc_bulk), masse volumique*chaleur spécifique de la partie fluide (hc_eau). En outre, le procédé de la présente invention comprend en outre l’étape consistant à la génération de sorties pour le temps de déposition correspondant, comportant la génération de cartes, de sections et de profils de puits avec les informations thermique et structurelle du bassin sédimentaire. En outre, selon un autre mode de réalisation préféré de la présente invention, il est défini un support de stockage lisible par ordinateur comprenant, stocké dans celui-ci, un ensemble d’instructions lisibles par ordinateur, qui, lorsqu’elles sont exécutées par un ordinateur, exécutent le procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole de la présente invention. Brève description des figures Afin de compléter la présente description et de mieux comprendre les caractéristiques de la présente invention, et selon un mode de réalisation préféré de celle-ci, un ensemble de figures est présenté ci-joint, où son mode de réalisation préféré est représenté d’une manière représentative, mais pas limitative. FIG. 1 La FIG. 1 illustre un organigramme du procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole. FIG. 2 La FIG. 2 montre un organigramme du procédé de solution numérique du procédé de la présente invention. FIG. 3 La FIG. 3 illustre un écran d’une application pour introduire les dimensions de la zone à l’étude, des paramètres thérmophysiques et des conditions de limite thermiques. FIG. 4 La FIG. 4 montre une interface d’une application pour introduire les caractéristiques de la lithosphère, le nombre de couches et des propriétés thérmophysiques. FIG. 5 La FIG. 5 montre un écran d’une application pour introduire des cartes, des lithofaciès, des couches salifères et ignées, la définition de la bibliothèque de lithofaciès est illustrée. FIG. 6 La FIG. 6 montre un écran d’une application pour introduire des valeurs pour la discrétisation du modèle numérique, l’âge du bassin et le temps, la période d’effondrement, des valeurs ou des cartes des facteurs d’étirement de la lithosphère, ou croûte et manteau. FIG. 7 La FIG. 7 présente des exemples de cartes utilisées comme entrée de données. FIG. 8 La FIG. 8 représente le domaine pour un modèle à 1D simplifié. FIG. 9 La FIG. 9 est un dessin schématique du processus d’effondrement lithosphérique. FIG. 10 La FIG. 10 montre le modèle schématique pour mettre à jour les propriétés dans les nœuds du maillage d’éléments finis. FIG. 11 La FIG. 11 illustre la division du maillage d’éléments finis. FIG. 12 La FIG. 12 montre des cartes de flux thermique pour tous les temps de dépositions sédimentaires. FIG. 12 La FIG. 12 illustre des cartes d’épaisseur élastique effective. FIG. 12 La FIG. 12 présente des paléogéométries du socle et des couches sédimentaires. FIG. 12 La FIG. 12 montre un graphe avec des résultats de flux thermique. FIG. 13 La FIG. 13 montre un graphe illustrant le couplement de l’effet de couverture, lorsque l’effet du taux de sédimentation et le processus de déposition sédimentaire sont pris en compte. FIG. 13 La FIG. 13 présente un graphe qui montre un scénario de bassin avec des sédiments (carte de lithofaciès) et effondrement (carte bêta). FIG. 13 La FIG. 13 présente un graphe qui montre un scénario de bassin sans sédiments. FIG. 13 La FIG. 13 illustre un graphe d’un scénario où la conductivité de sédiments est égale à 2,5 W/m/K et la conductivité de la croûte est égale à 2,5 W/m/K. FIG. 13 La FIG. 13 illustre un graphe d’un scénario où la conductivité de sédiments est égale à 1,25 W/m/K et la conductivité de la croûte est égale à 2,5 W/m/K. FIG. 13 La FIG. 13 illustre un graphe d’un scénario où la conductivité de sédiments est égale à 5 W/m/K et la conductivité de la croûte est égale à 2 W/m/K. Description détaillée de l’invention La FIG. 1 illustre un organigramme du procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole, selon un mode de réalisation de la présente invention, et sa description respective présentée ci-dessous explicite les étapes et les calculs impliqués dans la simulation en 3D de processus thermiques et structurels dans des bassins sédimentaires, en prenant compte de l’implémentation de l’effet de couverture, ce qui aura un influence significative sur l’évolution thermique des zones analysées, au cours du temps géologique. Le procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole implique le couplement d’événements thermiques d’origine sédimentaire et lithosphérique, générant de résultats raffinés de températures, flux thermique et paléogéométries de bassins sédimentaires sous forme de cartes en 3D. Le procédé comprend l’inclusion de l’effet de la déposition de sédiments dans les résultats thermiques de bassins, par le biais de l’accouplement entre les portions lithosphérique et sédimentaire dans le modèle thermique computationnel, ce qui améliore la précision des résultats en comparaison avec des données de puits, spécialement ceux concernant la température et la vitrinite. Les résultats obtenus au moyen du procédé de la présente invention peuvent être utilisés en tant que conditions de limite thermiques et structurelles pour l’intégration avec un autre logiciel de modélisation de systèmes de pétrole et également pour la calibration thermique directe, par comparaison des résultats obtenus à partir du modèle avec des données de puits d’une région ou bassin à l’étude. En particulier, l’effet de couverture est résolu dans la boucle temporelle, en combinant les effets lithosphériques et sédimentaires, respectivement. En particulier, le procédé calcule : a) le processus de diffusion thermique ; b) les effets d’advection à cause de l’effondrement de la lithosphère ; c) les effets d’advection à cause de la déposition de sédiments (effet de couverture) ; d) la génération de chaleur radiogénique. Le procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole comprend les étapes suivantes : l’obtention de données lithosphériques et sédimentaires d’entrée 1, dans lequel les données lithosphériques et sédimentaires comprennent au moins l’un de : dimensions de la zone ; paramètres thérmophysiques, tels que : conductivité de la croûte (W/m/K) ; masse volumique de la croûte kg/m³) ; chaleur spécifique de la croûte (J/kg/K) ; chaleur radiogénique de la croûte (μW/m³) ; conductivité du manteau lithosphérique (W/m/K) ; masse volumique du manteau lithosphérique (kg/m³) ; chaleur spécifique du manteau lithosphérique (J/kg/K) ; chaleur radiogénique du manteau lithosphérique (μW/m³) ; masse volumique de l’eau (kg/m³) ; conductivité de l’eau (W/m/K) ; chaleur spécifique de l’eau (J/kg/K) ; masse volumique de l’asthénosphère (kg/m³) ; coefficient d’élasticité (N/m²) ; coefficient de Poisson ; coefficient de dilatation thermique (1/°C) ; accélération de la gravité (m/s²) ; température isotherme pour Te (°C) ; conditions de limite dans la lithosphère, telles que la température de surface et la température de base ; données lithosphériques, correspondant au nombre de couches et aux propriétés thérmophysiques de la couche, comprenant au moins l’une de : épaisseur (m) ; masse volumique (kg/m³) ; chaleur spécifique (J/kg/K) ; conductivité (W/m/K) ; chaleur radiogénique (μW/m³) ; cartes d’horizons et lithofaciès pour chaque couche sédimentaire ; dans lequel les cartes de salifères et ignées sont optionnelles, en fonction de la région d’intérêt ; valeurs pour la discrétisation du modèle numérique, dans lequel les valeurs peuvent être librement définis par l’utilisateur dans la colonne de « nombre d’étapes » ou « séparation », pour la discrétisation temporelle ; et dans les colonnes de « nombre d’éléments » ou « séparation », pour la discrétisation spatiale dans les dimensions x, y et z (tel qu’illustré dans la FIG. 6 ). Une discrétisation plus importante est suggérée : 1) dans la période d’effondrement (des milliers d’années) et 2) dans la partie sédimentaire (des dizaines ou des centaines de mètres dans l’épaisseur en z correspondant à la profondeur maximale du bassin) ; pour les autres périodes et zones, une discrétisation d’entre 2-3 mille mètres garantira la convergence et des bons résultats) ; âge du bassin (en millions d’années – mya) ; temps et période de l’effondrement (en millions d’années – mya) ; valeurs ou cartes des facteurs d’étirement de la lithosphère (β) ou de la croûte et du manteau (δ et β) ; la discrétisation du modèle numérique, comportant la discrétisation du modèle différemment dans le domaine sédimentaire et dans le domaine lithosphérique 2 ; le démarrage du calcul de l’effet de couverture 3 ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire t est supérieur au temps présent et si l’âge du bassin est supérieure au temps de déposition correspondant à chaque couche sédimentaire t 4 et l’abstention de réaliser le procédé 5 ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire t est inférieur ou égal au temps présent et si l’âge du bassin est inférieure ou égale au temps de déposition correspondant à chaque couche sédimentaire t 4 et la vérification de si elle se trouve dans la période d’effondrement 6, dans lequel : si elle se trouve dans la période d’effondrement, le calcul de la période d’effondrement 7, comportant le calcul de la chaleur transitoire à cause de l’amincissement lithosphérique, ce qui provoque des variations dans le flux thermique, à partir du modèle thermique advectif : (Équation 1) où : v s désigne la vitesse des particules solides ; θ est la porosité ; ρ est la masse volumique ; c est la constante de décroissance ; et λ est la conductivité ; le calcul de la valeur de la croûte δ et la valeur du manteau β au sein de l’effondrement, selon l’équation 2 : (Équation 2) le calcul des épaisseurs et des horizons des couches de sédiment en fonction de la valeur de la croûte δ et la valeur du manteau β ; la quantification de la chaleur radiogénique dans la portion supérieure de la croûte. si elle ne se trouve pas dans la période d’effondrement ou si elle se trouve après l’étape de calcul de l’effondrement 7, le calcul de la sédimentation 8, comportant le calcul de l’effet de couverture à cause de la déposition de sédiments, par le biais d’un modèle advectif, comportant : le décompactage des couches de sédiment (backstripping) ; l’empilement et le compactage de couches de sédiment jusqu’au temps de déposition correspondant et selon les paramètres thermiques calculés à nouveau en profondeur ; le nouveau calcul de la profondeur du socle ; le calcul de la chaleur radiogénique dans la portion sédimentaire ; le calcul de la vitesse de déplacement vertical advective à cause de l’effondrement 9, pour chaque coordonnée XY, en utilisant la formulation vbulk = G*(h-z) pour chaque z coordonnée, où vbulk = vitesse moyenne du soulèvement de l’asthénosphère à cause de l’effondrement, variant linéairement en z ; G = magnitude du gradient de vitesse verticale le long de la lithosphère h, indiquant combien se déplace verticalement la vitesse de fluide (asthénosphère), le long de l’épaisseur définie comme « h » dans la formulation ; h = limite de profondeur de la lithosphère. Une fois réalisé ce calcul, les vitesses advectives sont ajoutées à la vitesse de sédimentation, dans lequel plus le compactage est important, plus le soulèvement à cause de l’effondrement est lent ; le calcul du flux de fluide 10, comportant la définition propriétés thermiques dans les nœuds du maillage numérique en 3D, comportant la conductivité, la chaleur radiogénique, la masse volumique*chaleur spécifique de la partie solide (hc_bulk), la masse volumique*chaleur spécifique de la partie fluide (hc_eau) l’utilisation d’éléments finis pour réaliser des calculs dans un moment donné 11 ; la réalisation de l’addition de l’effet thermique des domaines lithosphérique et sédimentaire à chaque étape de temps 12 ; la génération de sorties pour le temps de déposition correspondant 13, comportant la génération de cartes, de sections et de profils de puits avec des informations thermal thermique et structurelle du bassin sédimentaire ; en particulier, l’étape d’utilisation d’éléments finis pour réaliser des calculs dans un moment donné 11 comporte l’utilisation du procédé de PCG – tolérance par défaut du solveur : 1.e-07 et traitement distribué, en utilisant OpenMP directement dans le code. La FIG. 3 illustre un écran d’une application qui met en œuvre le procédé de la présente invention pour introduire les dimensions de la zone à l’étude, les paramètres thérmophysiques et les conditions de limite thermiques. La FIG. 4 montre une interface d’une application qui met en œuvre le procédé de la présente invention pour introduire les caractéristiques de la lithosphère, le nombre de couches et les propriétés thérmophysiques. La FIG. 5 montre un écran d’une application qui met en œuvre le procédé de la présente invention pour introduire des cartes, des lithofaciès, des couches salifères et ignées ; la définition de la bibliothèque de lithofaciès est illustrée. La FIG. 6 montre un écran d’une application qui met en œuvre le procédé de la présente invention pour introduire des valeurs pour la discrétisation du modèle numérique, l’âge du bassin et le temps, la période d’effondrement, des valeurs ou des cartes des facteurs d’étirement de la lithosphère (β) ou de la croûte et du manteau (δ et β). La bibliothèque de lithofaciès est crée pour définir les lithologies distribuées dans chaque couche sédimentaire. La bibliothèque comporte un ID, utilisé pour identifier les lithofaciès dans les cartes d’horizons, et des propriétés correspondantes, associées avec l’ID, tel qu’indiqué dans le tableau 1 : Tableau 1 : propriétés pour les lithologies composeront chaque couche sédimentaire Masse volumique kg/m³ Porosité de surface % Coefficient de compactage 1/m Diffusivité m²/s Conductivité W/m/K Chaleur radiogénique μW/m³ Concernant les cartes, une pluralité de cartes peuvent être utilisées comme entrée de données : des horizons, des lithofaciès, des conditions de limite, l’épaisseur de salifères, des facteurs d’étirement, des propriétés, entre autres. Les cartes sont fournies sous forme de fichiers, typiquement avec le format . xyz (coord. x., coord. y., valeur z). La FIG. 7 présente des exemples de cartes utilisées comme entrée de données dans le procédé de la présente invention. Plus spécifiquement, ce qui caractérise l’effet de couverture ou l’effet de couverture est l’implémentation où les paramètres de l’équation 1, θ , ρ , c et λ (porosité, masse volumique, constante de décroissance et conductivité, respectivement) varient en fonction de la relation de porosité/profondeur de la matrice sédimentaire. Selon la loi d’Athy (Athy L. F., 1930, « Density, porosity, and compaction of sedimentary rocks », AAPG Bulletin, v. 14, no. 1, pp. 1-24), la variation de porosité selon la profondeur z est donnée par l’équation 3 ci-dessous : (Équation 3), où est la porosité de surface et c est le coefficient de compactage de la matrice sédimentaire. Avec ceci, les autres paramètres de l’équation 1 varieront selon des valeurs de porosité/profondeur données en fonction de la variable θ. La conductivité effective de la matrice sédimentaire est obtenue à partir de la moyenne géométrique entre la conductivité de la phase solide λ_s et la conductivité de l’eau qui remplit la partie poreuse, λ_w, autrement dit : (Équation 4). La chaleur spécifique effective (ρc) est obtenue par la moyenne pondérée entre la partie solide (ρc) s et la partie fluide (ρc) w , c’est-à-dire : (Équation 5). La chaleur radiogénique effective (ρr) est obtenue seulement en pondérant la partie solide, c’est-à-dire : (Équation 6). En outre, les vitesses advectives de la partie solide et la partie fluide correspondant aux sédiments sont obtenues à partir de la loi de conservation de masse appliquée aux couches sédimentaires. La formulation pour le cas en 1D avec seulement une couche sédimentaire est développée pour un domaine, selon la FIG. 8 , qui est la représentation du domaine pour un modèle en 1D simplifié, pris et adapté au français à partir de Hutchison, I., 1985 ; The effects of sedimentation and compaction on oceanic heat flow, Geophysical Journal of the Royal Astronomical Society 82 (3) 439{459. doi :10.1111/j.1365-246X.1985.tb05145.x.3. La vitesse de déposition à z=0 est : (Équation 7). En outre, la vitesse de mouvement du socle à z=B est : (Équation 8). Alors, la vitesse de mouvement du socle pour toute t, toute z, sera donné par : (Équation 9). La méthodologie complète est donnée en ajoutant l’effet de la déposition sédimentaire, l’effet de couverture, au processus d’effondrement dans la lithosphère et la chaleur radiogénique, dans les domaines lithosphérique et sédimentaire, qui sont insérés dans l’équation de diffusion thermique/advection dans un milieu solide en 3D (équation 1). L’effondrement est le processus de rétrécissement/amincissement de la lithosphère et de déplacement de matériau de l’asthénosphère dans la direction verticale, tel que montré dans la FIG. 9 . La FIG. 9 représente schématiquement le processus d’effondrement lithosphérique. Sur la gauche, la lithosphère se trouve dans une configuration normale. Sur la droite, la lithosphère est devenue mince par un paramètre β et il y a un soulèvement conséquent de matériau asthénosphérique. C’est ainsi que l’asthénosphère est soulevée, à partir de la limite de profondeur , afin de remplir l’espace laissé par la sortie de matériau lithosphérique à cause de son amincissement et ceci est représenté dans le terme advectif de l’équation 1, avec le processus de déposition sédimentaire (vitesse verticale). L’effet de couverture de sédiments et d’advection à cause de l’effondrement sont des processus en compétition verticale, de sorte que plus le compactage est important, plus la vitesse de soulèvement à cause de l’amincissement lithosphérique est lente. Solution numérique en 3D du problème de transport de la chaleur transitoire dans des éléments finis : La solution numérique de la méthodologie présentée pour quantifier l’effet de couverture couplé avec d’autres processus thermiques a été développée dans des éléments finis. La solution pour le cas en trois dimensions est obtenue en évaluant le compactage indépendamment à chaque nœud du maillage d’éléments finis. Le modèle numérique utilise des éléments hexaédriques H8, avec 8 nœuds pour la solution dans la direction x-y et des éléments avec 12 nœuds H12, dans la direction z. La déposition de sédiments et l’effondrement causent des modifications dans les propriétés géologiques de chaque élément fini au cours du temps, c’est-à-dire, les éléments du maillage sont soumis à des variations dans les propriétés géologiques de leurs domaines ( FIG. 11 ). Dans l’implémentation, on suppose que les propriétés sont interpolées par les mêmes fonctions de forme de l’élément concerné. La FIG. 10 montre le modèle schématique pour mettre à jour les propriétés dans les nœuds du maillage d’éléments finis. Dans la FIG. 10 , (.) i désigne les valeurs nodales des propriétés géologiques obtenues à partir de la forme exacte du domaine dans l’instant concerné. De manière générale, le code reçoit comme données d’entrée la distribution spatiale (en x, y) des épaisseurs des couches sédimentaires et leurs temps de déposition correspondants. Pendant la simulation, l’épaisseur de chaque couche z_i est évaluée selon l’instant actuel. Une fois évaluées les épaisseurs actuelles, la porosité de chaque nœud du maillage d’éléments finis est calculée à partir de la loi d’Athy et les propriétés thermiques sont ensuite évaluées selon les équations pour la variation des propriétés thermiques. Une fois obtenues les propriétés nodales, le modèle d’éléments finis est défini, par interpolation des propriétés thermiques avec les mêmes fonctions de forme de l’élément sélectionné, pour générer le maillage. Formulation de solution numérique : procédé du gradient conjugué préconditionné (PCG) La solution numérique du procédé de la présente invention utilise le procédé implicite/explicite d’intégration dans le temps, arrivant au système algébrique d’équations linéaires suivant : (Équation 10). où la matrice du côté gauche est une matrice définie positive symétrique. Trois caractéristiques principales de ce système d’équations se distinguent : symétrie, creuse et propriété définie positive. La symétrie, en combinaison avec la propriété définie positive, fait que le système soit résolu de manière générique par le procédé du gradient conjugué (CGM). Étant donné qu’il s’agît d’un procédé itératif (sans modification de la structure du système), il est possible d’utiliser des structures de données creuses pour stocker le système linéaire. Dans l’implémentation du modèle décrit dans le présent procédé, on utilise une variation du CGM, qui est le procédé du gradient conjugué préconditionné (PCG), qui est une stratégie pour optimiser le procédé traditionnel, améliorant le conditionnement de la matrice A du système linéaire (équation 10), par le biais d’une multiplication préalable du système par une matrice qui est proche de , mais qui peut être facilement inversée. Alors, le système équivalent suivant peut être résolu : (Équation 11). La nouvelle matrice présente un conditionnement meilleur que . La matrice inverse est obtenue en inversant chacun des éléments de la diagonale principale de la matrice et cela est connu sous le nom de technique de Jacobi, utilisé dans le simulateur numérique implémenté pour résoudre le procédé proposé. L’approximation numérique du problème thermique est donc résolue en utilisant plusieurs d’itérations à partir d’une approximation initiale, jusqu’à atteindre une convergence en fonction de la tolérance à erreur définie entre l’itération actuelle et l’itération précédente : Approximation initiale : (Équation 12) (Équation 13) Itérations suivantes : (Équation 14) (Équation 15) (Équation 16) (Équation 17) (Équation 18). L’organigramme dans la FIG. 2 illustre le procédé mis en œuvre pour la solution numérique du procédé de la présente invention tel que présenté, en utilisant l’approximation d’éléments finis. Parallélisation du code développé pour implémenter le procédé de la présente invention Pour réduire le temps computationnel du procédé de simulation pour le calcul de l’effet de couverture, les sections les plus critiques du code ont été parallélisées. Le code a été parallélisé selon la bibliothèque d’OpenMP, qui a été intégrée dans la méthodologie présentée. Le parallélisme a été implémenté en deux étapes, pour optimiser les étapes auxquelles est consacrée la plus part du temps de traitement : Tout d’abord, le calcul et l’assemblage des matrices globales du FEM sont parallélisés. Dans la deuxième étape, la solution du système linéaire du problème de transport de la chaleur est parallélisée, plus précisément, la multiplication de la matrice globale dans un format creux au sein du processus itératif du PGC. La parallélisation de l’ensemble et l’évaluation des matrices globales du FEM ont été réalisées en fonction d’une division disjointe du maillage d’éléments finis. Cette stratégie est essentielle pour que la parallélisation du code soit réalisée efficacement, sans problèmes de concurrence critique et sans besoin de créer des verrous de type mutex. À cause de la structure de maillage du maillage de FEM utilisé dans le procédé, la division a été réalisé facilement, en subdivisant tout simplement les couches verticales en couches pairs et impairs, tel que montré dans la FIG. 11 , qui illustre la division du maillage d’éléments finis. Résultats du procédé de la présente invention En utilisant le procédé de la présente invention, on obtient les résultats suivants, qui peuvent être utilisés directement ou comme entrée pour des étapes ultérieures dans la modélisation de systèmes de pétrole, tel que montré dans le tableau 2 ci-dessous. Tableau 2 : résultats en 1D, 2D et 3D obtenus à partir du procédé de la présente invention Température de socle (°C) Température de surface de bassin (°C) Flux thermique de socle (mW/m²) Flux thermique de surface de bassin (mW/m²) Épaisseur sédimentaire (m) Taux de sédimentation (m/Mya.) Épaisseur élastique effective – Te (m) Températures dans chaque couche sédimentaire (°C) Flux thermique dans chaque couche sédimentaire (mW/m²) Flux thermique à Moho (mW/m²) Température à Moho (°C) Cartes de profondeur isothermes (m) Cartes de paléobathymétrie (m) Cartes de profondeur de couches sédimentaires à chaque temps de déposition (m) Cartes de restauration (profondeur) de socle (m) En outre, les résultats conventionnels sont générés sous forme de cartes de la région à l’étude, comme les exemples montrés dans la FIG. 12 , la FIG. 12 , la FIG. 12 et la FIG. 12 . La FIG. 12 montre des cartes de flux thermique pour tous les temps de dépositions sédimentaires. La FIG. 12 illustre des cartes d’épaisseur élastique effective (Te), en fonction des températures. La FIG. 12 présente des paléogéométries du socle et des couches sédimentaires. La FIG. 12 montre un graphe avec des résultats de flux thermique : réduction de flux thermique dans la partie supérieure et la partie inférieure du bassin (Qs et Qb), en utilisant l’effet de couverture, en comparaison avec le modèle de Jarvis&McKenzie traditionnel (Jarvis, G.T., Mckenzie D.P., 1980 ; D. P. ; « Sedimentary Basin Formation with Finite Extension Rates », Earth and Planetary Science Letters, 48, pp. 42-52.). La FIG. 13 , la FIG. 13 , la FIG. 13 , la FIG. 13 , la FIG. 13 et la FIG. 13 montrent certains exemples de résultats de flux thermique obtenus à partir de l’implémentation du procédé de la présente invention. Le calcul thermique sans l’effet de couverture a pour résultat le profil de flux thermique très simplifié de la FIG. 13 , variant seulement selon le processus d’effondrement défini pour le bassin et aucun détail additionnel n’est observé dans les caractéristiques thermiques de la région à l’étude. Lorsque l’effet de couverture est couplé, lorsque l’effet du taux de sédimentation et le processus de déposition sédimentaire sont pris en compte ( FIG. 13 ), plus de détails sont observés dans les profils de flux thermique obtenus ( FIG. 13 , FIG. 13 , FIG. 13 , FIG. 13 ). Dans ces exemples, des variations dans les profils thermiques sont observées à cause de l’utilisation de différentes valeurs de conductivité de sédiments, celle-ci est la caractéristique principale de l’effet de couverture. Avantages et application de la présente invention La présente invention proposes l’intégration de l’effet de couverture dans la quantification de la température et du flux thermique dans des bassins sédimentaires et il est pertinent dans des régions où les épaisseurs sédimentaires sont supérieures à 3 km, les conductivités thermiques sont inférieures à 3,5 W/m/K, et les taux de sédimentation sont supérieurs à 200 m/Mya. En outre, des variations dans les propriétés thérmophysiques résultant de la déposition de sel et du compactage contribuent aux variations dans la structure thermique de la région à l’étude. Dans ce contexte, l’application du procédé proposé, par le biais de l’utilisation d’un modèle en 3D, qui présente la capacité de capturer des différences géologiques dans une zone large, fournira des résultats plus précis avec plus de résolution et de détail à une échelle régionale. Comme exemple, selon des études menés dans la portion proximale du bassin de Santos, il est estimé que le flux thermique dans le socle est réduit de 36% en comparaison avec des modèles sans l’application de l’effet de couverture présenté dans le présent document, étant donné que les taux de sédimentation sont élevés (environ 200 m/Mya.) dans des régions où les épaisseurs sédimentaires pré-salifère et post-salifère sont d’environ 3,5 km et 4 km, respectivement. Dans la portion distale du bassin de Santos, où les épaisseurs des dispositions pré-salifère et post-salifère sont de moins de 1,5 km et l’épaisseur de salifère varie entre 3 et 5 km, l’effet de couverture n’est pas significatif, une différence d’environ 2% étant observée par rapport au modèle d’effondrement étendu traditionnel. Dans la pratique, le procédé décrit dans le présent document a déjà été appliqué à des analyses de bassins tels que Santos/Campos, Pelotas, marge équatoriale, Foz do Amazonas et Namibia, fournissant des résultats excellents pour le flux thermique et les températures. Ces résultats, obtenus sous forme de cartes en 3D, ont été comparés avec des données de puits, par rapport aux paléotempératures et réflectances de la vitrinite, et ont montré une calibration très bonne. Outils utilisés Le procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole de la présente invention, implémenté par le biais d’un programme informatique, comprend les outils suivants dans un environnement de programme informatique pour permettre l’utilisation de celui-ci : un simulateur numérique avec optimisation du temps computationnel ; et une interface graphique pour l’introduction de données, l’exécution du simulateur numérique et l’obtention des résultats. Le simulateur numérique est responsable du couplement des événements lithosphériques avec de processus sédimentaires (effet de couverture) qui n’étaient pas quantifiés auparavant dans le calcul thermique du bassin. Avec ceci, les simulations réalisées dans les applications de simulation commencent à quantifier l’effet thermique transitoire lié à la déposition et au compactage de sédiments froids, connu comme l’effet de couverture, intégrant celui-ci avec les autres processus existant dans le modèle, tels que l’évolution du gradient géothermique, l’étirement lithosphérique (effondrement) et la chaleur radiogénique présente dans la portion supérieure de la croûte. Une interface graphique peut être conçue spécialement pour permettre l’utilisation du procédé de la présente invention, qui est facile et intuitive d’utiliser. Par exemple, une interface graphique peut être développée pour appliquer les algorithmes qui définissent la séquence proposée de processus, permettant l’introduction d’une grande quantité de données et la fourniture de résultats avec un niveau élevé de détail à une échelle régionale, pour les zones à l’étude. Le procédé pour le calcul de l’effet de couverture dans la modélisation de systèmes de pétrole de la présente invention produit des résultats améliorés, avec le couplement d’événements lithosphériques avec des processus sédimentaires (couverture), qui n’étaient pas quantifiés auparavant dans le calcul thermal du bassin. Autrement dit, les simulations réalisées dans les programmes de simulation ont commencé à quantifier l’effet thermique transitoire lié à la déposition et au compactage de sédiments froids, connu comme l’effet de couverture, intégrant celui-ci avec les autres processus existant dans le modèle, tels que l’évolution du gradient géothermique, l’étirement lithosphérique (effondrement) et la chaleur radiogénique présente dans la portion supérieure de la croûte. Ce processus, connu sous le nom d’effet de couverture, a comblé l’une des lacunes importantes qui manquaient pour le raffinement du calcul thermique des bassins sédimentaires. Interface graphique développée en C++/QT Le simulateur numérique et l’interface graphique peuvent être développés dans les langages C++ et QT, respectivement, en plus de l’utilisation de bibliothèques du domaine publique, telles que OpenMP, pour optimiser le temps computationnel, entre autres. Pour commencer le procédé, il est nécessaire de choisir le type de simulation (« thermique » ou « backstripping en 3D + thermique ») dans l’interface graphique. L’interface graphique développée en C++/QT permet l’application du procédé de la présente invention dans des cas réels de bassins et la présentation de résultats sous forme de cartes de flux thermique et de température. L’interface graphique est conçue et implémentée spécialement pour l’implémentation du procédé de la présente invention et permet le développement de modèles thermomécaniques pour la simulation de bassins d’une manière intuitive, présentant des résultats avec un niveau élevé de détail sous forme de fichiers de texte, qui peuvent être importés vers des outils de visualisation. En complément, la présente invention concerne un support de stockage lisible par ordinateur, qui comprend, stocké dans celui-ci, un ensemble d’instructions lisibles par ordinateur, dans lequel, lorsque l’ensemble d’instructions lisibles par ordinateur est exécuté par un ou plusieurs processeurs, l’un ou plusieurs processeurs mettent en œuvre le procédé de la présente invention, tel que décrit ci-dessus. En particulière, le support de stockage lisible par ordinateur peut être une mémoire, dans lequel la mémoire peut être d’un type non volatile, telle qu’une unité de disque dur (HDD) ou une unité à semi-conducteurs (SSD), ou il peut s’agir d’une mémoire volatile, telle qu’une mémoire vive (RAM). En outre, le support de stockage lisible peut être tout autre moyen ou support qui peut transporter ou stocker ou enregistrer le code de programme prévu sous forme d’une instruction ou d’une structure de données ou d’un ensemble de instructions, et auquel peut avoir accès un ou plusieurs ordinateurs ou un ou plusieurs processeurs, mais il n’est pas ainsi limité. Alternativement, le support de stockage lisible peut être un circuit ou tout autre dispositif ou moyen qui peut mettre en œuvre une fonction de stockage ou de transport ou d’enregistrement, tel qu’un signal ou une porteuse. Spécifiquement, l’ensemble d’instructions lisibles par ordinateur représente l’algorithme ou le code de programme informatique ou une structure de données qui réalise le procédé de la présente invention tel que décrit ci-dessus. Le processeur peut être un processeur à usage général, qui peut être un microprocesseur ou tout processeur conventionnel ou similaire. L’homme du métier appréciera les enseignements présentés dans le présent document et sera capable de reproduire l’invention dans les modes de réalisation présentés et dans d’autres variantes, englobées par la portée des revendications ci-jointes. Procédé mis en œuvre par ordinateur pour l’intégration de l’effet de couverture dans la quantification de la température et du flux thermique d’un bassin sédimentaire dans la modélisation de systèmes de pétrole, caractérisé en ce qu’il comprend les étapes consistant à : l’obtention de données lithosphériques et sédimentaires d’entrée (1) ; la discrétisation du modèle différemment dans le domaine sédimentaire et dans le domaine lithosphérique (2) ; le démarrage du calcul de l’effet de couverture (3) ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire (t) est supérieur au temps présent et si l’âge du bassin est supérieure au temps de déposition correspondant à chaque couche sédimentaire (t) (4), en s’abstenant de réaliser le procédé (5) ; la vérification de si le temps de déposition correspondant à chaque couche sédimentaire (t) est inférieur ou égal au temps présent et si l’âge du bassin est inférieure ou égale au temps de déposition correspondant à chaque couche sédimentaire (t) (4) et la vérification de si elle se trouve dans la période d’effondrement (6), dans lequel : si elle se trouve dans la période d’effondrement, le calcul de l’effondrement (7) ; si elle ne se trouve pas dans la période d’effondrement, ou si elle se trouve en séquence par rapport à l’étape de calcul de l’effondrement (7), le calcul de la sédimentation (8) ; le calcul de la vitesse de déplacement vertical advective à cause de l’effondrement (9) pour chaque coordonnée XY ; le calcul du flux de fluide (10) ; l’utilisation d’éléments finis pour réaliser le calcul à chaque instant dans le temps (11) ; la réalisation de l’addition de l’effet thermique des domaines lithosphérique et sédimentaire à chaque étape de temps (12) ; la fourniture de résultats de flux thermique et de température. Procédé selon la revendication 1, caractérisé en ce que les données lithosphériques et sédimentaires comprennent au moins l’un de : dimensions de la zone ; paramètres thérmophysiques, tels que : conductivité de la croûte (W/m/K) ; masse volumique de la croûte (kg/m³) ; chaleur spécifique de la croûte (J/kg/K) ; chaleur radiogénique de la croûte (μW/m³) ; conductivité du manteau lithosphérique (W/m/K) ; masse volumique du manteau lithosphérique (kg/m³) ; chaleur spécifique du manteau lithosphérique (J/kg/K) ; chaleur radiogénique du manteau lithosphérique (μW/m³) ; masse volumique de l’eau (kg/m³) ; conductivité de l’eau (W/m/K) ; chaleur spécifique de l’eau (J/kg/K) ; masse volumique de l’asthénosphère (kg/m³) ; coefficient d’élasticité (N/m²) ; coefficient de Poisson ; coefficient de dilatation thermique (1/°C) ; accélération de la gravité (m/s²) ; température isotherme pour Te (°C) ; conditions de limite dans la lithosphère, telles que la température de surface et la température de base ; données lithosphériques, correspondant au nombre de couches et aux propriétés thérmophysiques de la couche, comprenant au moins l’une de : épaisseur (m) ; masse volumique (kg/m³) ; chaleur spécifique (J/kg/K) ; conductivité (W/m/K) ; chaleur radiogénique (W/m³*1e-6) ; cartes d’horizons et lithofaciès pour chaque couche sédimentaire ; dans lequel les cartes de salifères et ignées sont optionnelles, en fonction de la région d’intérêt ; valeurs pour la discrétisation du modèle numérique ; âge du bassin (en millions d’années – mya.) ; temps et période de l’effondrement (en millions d’années – mya.) ; valeurs ou cartes des facteurs d’étirement de la lithosphère (β) ou de la croûte et du manteau (δ et β). où : v s désigne la vitesse des particules solides ; θ est la porosité ; ρ est la masse volumique ; c est la constante de décroissance ; et λ est la conductivité ; le calcul de la valeur de la croûte δ et la valeur du manteau β au sein de l’effondrement, selon : le calcul des épaisseurs et des horizons des couches de sédiment en fonction de la valeur de la croûte δ et la valeur du manteau β ; la quantification de la chaleur radiogénique dans la portion supérieure de la croûte. Procédé selon la revendication 1, caractérisé en ce que le calcul de la sédimentation (8) comporte le calcul de l’effet de couverture à cause de la déposition de sédiments, par le biais d’un modèle advectif, comportant : le décompactage des couches de sédiment (backstripping) ; l’empilement et le compactage de couches de sédiment jusqu’au temps de déposition correspondant et selon les paramètres thermiques calculés à nouveau en profondeur ; le nouveau calcul de la profondeur du socle ; le calcul de la chaleur radiogénique dans la portion sédimentaire. Procédé selon la revendication 1, caractérisé en ce que le calcul de la vitesse de déplacement vertical advective à cause de l’effondrement (9), pour chaque coordonnée XY, comporte le calcul par le biais de la formulation vbulk = G*(h-z) pour chaque coordonnée verticale z, où vbulk = vitesse moyenne du soulèvement de l’asthénosphère à cause de l’effondrement, variant linéairement en z ; G = magnitude du gradient de vitesse verticale le long de la lithosphère h, indiquant combien se déplace verticalement la vitesse de fluide (asthénosphère), le long de l’épaisseur définie come « h » dans la formulation ; h = limite de profondeur de la lithosphère ; comportant en outre l’addition des vitesses advectives à la vitesse de sédimentation, dans lequel plus le compactage est important, plus le soulèvement à cause de l’effondrement est lent. Procédé selon la revendication 1, caractérisé en ce que le calcul du flux de fluide (10) comporte la définition de propriétés thermiques dans les nœuds du maillage numérique en 3D, comprenant au moins l’une de conductivité, chaleur radiogénique, masse volumique*chaleur spécifique de la partie solide (hc_bulk), masse volumique*chaleur spécifique de la partie fluide (hc_eau). Procédé selon la revendication 1, caractérisé en ce qu’il comprend en outre l’étape consistant à la génération de sorties pour le temps de déposition correspondant (13), comportant la génération de cartes, de sections et de profils de puits avec les informations thermique et structurelle du bassin sédimentaire. Support de stockage lisible par ordinateur, caractérisé en ce qu’il comprend, stocké dans celui-ci, un ensemble d’instructions lisibles par ordinateur, qui, lorsqu’elles sont exécutées par un ordinateur, exécutent le procédé selon l’une quelconque des revendications 1 à 7.