FR 3076626 A1 20190712 FR 1860364 A 20181109 PRÉDICTION DE PROPRIÉTÉS DE ROCHES SISMIQUES EN TEMPS AVANCÉ SUR LA BASE D'UNE ANALYSE SISMIQUE 4D DOMAINE DE L'INVENTION [0001] La présente invention concerne de manière générale l'analyse sismique de formations rocheuses souterraines, et en particulier des techniques d'analyse sismique permettant de modéliser des propriétés variables dans le temps de formations rocheuses souterraines. CONTEXTE [0002] Dans l'exploration et la production d'hydrocarbures, les techniques d'imagerie sismique peuvent être utilisées pour mieux comprendre la profondeur et les structures de formations géologiques souterraines. Diverses sources sismiques, telles que la dynamite, les camions « chute de poids », les canons à air comprimé et d'autres sources de bruit situées à la surface d'une formation contenant des hydrocarbures, peuvent être utilisées pour propager des ondes d'énergie sismique à travers une formation souterraine. Les ondes propagées sont réfléchies à travers la formation et acquises à l'aide de divers dispositifs récepteurs de signaux sismiques, par exemple, des géophones, des hydrophones, et similaires. Un enregistrement des réflexions des ondes sismiques acquises à différents moments dans la formation souterraine peut ensuite être utilisé pour estimer les propriétés de roches de différentes zones de formation susceptibles de présenter un intérêt pour l'exploration et/ou la production d'hydrocarbures. Par exemple, les données sismiques peuvent être utilisées pour générer des modèles tridimensionnels (3D) de structures géologiques souterraines, y compris des failles et d'autres caractéristiques stratigraphiques susceptibles de piéger des gisements d'hydrocarbures et de minéraux, à l'intérieur d'une zone d'intérêt. Ces modèles 3D peuvent être utilisés pour simuler l'écoulement de fluides à l'intérieur de la formation souterraine. Les résultats d'une telle simulation peuvent ensuite être utilisés pour effectuer diverses activités liées à la production et à l'évaluation économique d'hydrocarbures, y compris, mais sans s'y limiter, la prédiction de la production future, la mise en place de nouveaux puits et l'évaluation de scénarios de gestion de réservoirs alternatifs. [0003] Diverses techniques de modélisation et de simulation sismiques peuvent être utilisées pour permettre l'intégration quantitative de données de terrain dans un modèle de réservoir pouvant être mis à jour. Ces techniques peuvent inclure, par exemple, la pétrophysique et la physique de roche intégrées pour déterminer la gamme des lithotypes et des propriétés de roches, l'inversion géostatistique pour déterminer un ensemble de modèles plausibles de propriétés de roches d'origine sismique à résolution et hétérogénéité verticales suffisantes pour la simulation d'écoulement, le 2 transfert de grille stratigraphique pour déplacer avec précision des données d'origine sismique sur le modèle géologique, et la simulation d'écoulement pour la validation et le classement de modèles afin de déterminer le modèle qui correspond le mieux à toutes les données. La simulation sismique réalisée à l'aide de telles techniques peut être considérée comme réussie si le modèle sous-jacent reflète avec précision les diagraphies de puits d'origine, les données sismiques et l'historique de production. Cependant, l'intégration de données sismiques du terrain dans un modèle dans le cadre d'un processus de simulation ne réussit pas toujours, car les données sismiques peuvent être inexactes, incomplètes ou totalement indisponibles. BRÈVE DESCRIPTION DES DESSINS [0004] [Fig. 1] La figure 1 est un diagramme schématique en coupe transversale d'un environnement de puits de forage pour l'acquisition de données sismiques à partir d'un relevé sismique terrestre d'une formation souterraine. [0005] [Fig. 2] La figure 2 est un diagramme schématique en coupe transversale d'un environnement marin pour l'acquisition de données sismiques à partir d'un relevé sismique marin d'une formation souterraine. [0006] [Fig. 3] La figure 3 est un schéma fonctionnel d'un système d'analyse et de modélisation sismiques d'une propriété de roche variable dans le temps d'une formation souterraine en temps avancé. [0007] [Fig. 4A, 4B, 4C] Les figures 4A à 4C sont des schémas d'un profil illustratif d'une propriété de roche variable dans le temps lors de sa transition sur des intervalles de temps successifs. [0008] [Fig. 5] La figure 5 est un graphique d'un profil illustratif pour une propriété de roche variable dans le temps au niveau d'un emplacement fixe à l'intérieur d'une formation souterraine. [0009] [Fig. 6] La figure 6 est un graphique d'un autre profil illustratif pour une propriété de roche variable dans le temps au niveau d'un emplacement fixe à l'intérieur d'une formation souterraine. [0010] [Fig. 7] La figure 7 est un organigramme d'un processus illustratif d'analyse et de modélisation sismiques d'une propriété de roche variable dans le temps d'une formation souterraine en temps avancé. [0011] [Fig. 8] La figure 8 est un schéma fonctionnel d'un exemple de système informatique dans lequel des modes de réalisation de la présente invention peuvent être mis en oeuvre. 3 DESCRIPTION DE MODES DE RÉALISATION ILLUSTRATIFS [0012] Des modes de réalisation de la présente invention concernent l'analyse et la modélisation sismiques de propriétés de roches variables dans le temps d'une formation souterraine en temps avancé. Bien que la présente invention soit décrite ici par référence à des modes de réalisation illustratifs pour des applications spécifiques, il faut comprendre que les modes de réalisation ne sont pas limités à ceux mentionnés ici. D'autres modes de réalisation sont possibles, et des modifications peuvent être apportées aux modes de réalisation conformément à l'esprit et à la portée des enseignements de la présente invention, ainsi que dans des domaines supplémentaires dans lesquels les modes de réalisation seraient d'une utilité significative. [0013] Il apparaîtrait également évident pour l'homme du métier que les modes de réalisation, tels que décrits ici, peuvent être mis en oeuvre dans de nombreux modes de réalisation différents de logiciels, de matériels, de micrologiciels et/ou d'entités illustrés dans les figures. Tout code logiciel réel conçu pour commander spécifiquement le matériel nécessaire à la mise en oeuvre des modes de réalisation ne limite pas la description détaillée. Ainsi, le comportement fonctionnel des modes de réalisation sera décrit étant entendu que des modifications et variations des modes de réalisation sont possibles, compte tenu du niveau de détail présenté ici. [0014] Dans la présente description détaillée, des références à « l'un des modes de réalisation », « un mode de réalisation », « un exemple de mode de réalisation », etc., indiquent que le mode de réalisation décrit peut inclure un aspect, une structure ou une caractéristique particulière, mais chaque mode de réalisation n'inclut pas nécessairement cet aspect, cette structure ou cette caractéristique particulière. En outre, de telles expressions ne se réfèrent pas nécessairement au même mode de réalisation. De plus, lorsqu'un aspect, une structure ou une caractéristique particulière est décrite en relation avec un mode de réalisation, il est estimé que la mise en oeuvre dudit aspect, de ladite structure ou caractéristique en relation avec d'autres modes de réalisation est connue de l'homme du métier, que cette mise en oeuvre soit ou non expressément décrite. [0015] Bien que des modes de réalisation de la présente invention puissent être décrits dans le contexte d'un puits de forage horizontal ou d'un puits de forage vertical, sauf indication contraire, il convient de noter que les modes de réalisation décrits conviennent également à une utilisation dans des puits de forage ayant d'autres orientations, y compris des puits de forage verticaux, des puits de forage inclinés, des puits de forage multilatéraux ou similaires. De même, sauf indication contraire, même si des modes de réalisation peuvent être décrits dans le contexte d'une opération en mer, il convient de noter que l'appareil selon la présente invention convient également pour une utilisation dans des opérations à terre et inversement. En outre, sauf indication contraire, même si 4 des modes de réalisation peuvent être décrits dans le contexte d'un trou tubé, il convient de noter que l'appareil selon la présente invention convient également à une utilisation dans des opérations à trous ouverts. [0016] Comme cela sera décrit plus en détail ci-dessous, des modes de réalisation de la présente invention peuvent être utilisés pour analyser et estimer des propriétés de roches variables dans le temps d'une formation souterraine. Dans un ou plusieurs modes de réalisation, des ensembles de données d'inversion sismique représentant des propriétés de roches variables dans le temps de la formation peuvent être acquis sur une pluralité d'intervalles de temps. Par exemple, de telles données peuvent inclure des mesures sismiques acquises périodiquement à différents moments afin de créer une vue à intervalles de temps d'une ou de plusieurs propriétés de roches de la formation. Ces données à intervalles de temps peuvent être appelées données sismiques «4D », où le temps est considéré comme une quatrième dimension dans l'acquisition, le traitement et l'analyse de relevés sismiques répétés pour évaluer des changements dans les propriétés de roches de la formation au fil du temps. Des exemples de telles propriétés de roches incluent, mais sans s'y limiter, l'impédance, le coefficient de Poisson, le rapport entre la vitesse de compression et la vitesse de cisaillement (Vp/Vs), la porosité, la perméabilité, la saturation en fluide, la pression et d'autres propriétés géomécaniques ou géophysiques de la formation. [0017] Dans un ou plusieurs modes de réalisation, les ensembles de données à intervalles de temps acquis peuvent être utilisés pour estimer ou prédire les propriétés de roches variables dans le temps de la formation à un moment futur. Par exemple, une analyse des ensembles de données acquis peut être utilisée pour déterminer les tendances de données existantes, pouvant être utilisées pour produire des réseaux de propriétés sismiques estimées dans le futur et/ou de multiples réalisations équiprobables d'une ou de plusieurs propriétés de roches sismiques inversées à un instant donné auquel les mesures sismiques sont manquantes ou doivent encore être acquises. Les propriétés de roches variables dans le temps estimées peuvent ensuite être co-analysées quantitativement et co-visualisées qualitativement avec les résultats d'une simulation d'écoulement afin d'évaluer comment les propriétés de roches variables dans le temps peuvent varier dans le temps en raison de déplacements de fluide et de changements de propriétés de roches élastiques. Les propriétés estimées peuvent également être utilisées pour prédire des changements dans des propriétés de roches variables dans le temps, dérivées de la modélisation pétro-élastique basée sur le temps, qui peuvent à leur tour être dérivées d'une simulation d'écoulement de fluide dans le milieu poreux de la formation. Le calcul et l'analyse des données de propriétés de roches estimées peuvent être effectués pour tout scénario futur lorsqu'il existe au moins deux études ou relevés sismiques inversés. Cependant, il convient de noter que les estimations des propriétés de roches sismiques en temps avancé peuvent être mieux limitées par des données de relevés à intervalles de temps sensiblement plus longs, par exemple, lorsque le nombre de relevés sismiques à intervalles de temps inversés est de loin supérieur à deux. [0018] Des modes de réalisation illustratifs de la présente invention et leurs méthodologies associées sont décrits ci-dessous par référence aux figures 1 à 8 car ils pourraient être utilisés, par exemple, dans un système informatique d'analyse et de modélisation sismiques de propriétés de roches variables dans le temps d'une formation souterraine. D'autres caractéristiques et avantages des modes de réalisation décrits seront évidents ou apparaîtront évidents à l'homme du métier lors de l'examen des figures et de la description détaillée suivantes. Il est prévu que toutes ces caractéristiques et tous ces avantages supplémentaires soient inclus dans la portée des modes de réalisation décrits. De plus, les figures illustrées ne représentent que des exemples et ne sont pas destinées à établir ou à suggérer toute limitation par rapport à l'environnement, l'architecture, la conception ou le processus dans lesquels différents modes de réalisation peuvent être mis en oeuvre. [0019] La figure 1 est un diagramme schématique en coupe transversale d'un environnement de puits de forage 100 pour l'acquisition de données sismiques à partir d'un relevé sismique terrestre d'une formation souterraine. L'environnement de puits de forage 100 inclut un derrick 102 positionné sur une surface 104 de la formation. Le derrick 102 peut supporter un train de forage 106 et d'autres composants de l'environnement de puits de forage 100. Le train de forage 106 peut inclure des tuyaux segmentés qui se prolongent au-dessous de la surface 104 dans un puits de forage 108. Le train de forage 106 peut transmettre le fluide de forage (ou la boue de forage) nécessaire pour faire fonctionner un trépan de forage 110 positionné à l'extrémité du train de forage 106. La boue transmise par le train de forage 106 peut fournir le couple nécessaire pour faire fonctionner le trépan de forage 110. Le poids du train de forage 106 peut exercer sur le trépan de forage 110 une force axiale qui, conjointement avec la rotation du trépan de forage 110, peut aider à forer le puits de forage 108 depuis la surface 104 à travers diverses couches de formation 112. [0020] Le train de forage 106 inclut un ensemble de fond de trou 114 positionné sur le train de forage 106 près du trépan de forage 110. L'ensemble de fond de trou 114 inclut une combinaison de divers composants y compris, mais sans s'y limiter, un ou plusieurs colliers de forage 116, un outil sismique 118 et un ensemble de moteur de fond de puits 120 logeant un moteur pour le trépan de forage 110. L'outil sismique 118 peut inclure un certain nombre de dispositifs de mesure pour mesurer diverses propriétés de la formation et d'autres données relatives à l'opération de fond de trou susceptible d'être en cours. De tels dispositifs de mesure peuvent inclure, par exemple, un 6 réseau de capteurs sismiques 122, par exemple, des géophones ou d'autres types de dispositifs de mesure sismique en fond de trou, pour détecter et enregistrer de l'énergie sismique en fond de trou. [0021] L'énergie sismique peut être générée par une source sismique 126 située à la surface 104 de la formation. La source sismique 126 peut être l'une quelconque des différentes sources d'énergie sismique qui peuvent être utilisées pour effectuer un relevé sismique de la formation. Des exemples de telles sources sismiques incluent, mais sans s'y limiter, un canon à air comprimé, une source sonore à effet plasma, un camion chute de poids, un ou plusieurs dispositifs explosifs, une source d'énergie à impulsions électromagnétiques (« EMP ») et un vibrateur sismique. Comme représenté dans la figure 1, le relevé sismique dans cet exemple peut être effectué en propageant de l'énergie sismique sous la forme d'ondes sismiques 124 depuis la source sismique 126 à travers les couches de formation 112 entourant le puits de forage 108. Certaines des ondes sismiques 124 générées par la source sismique 126 peuvent être réfléchies ou réfractées par les couches de formation 112 et ensuite échantillonnées par les capteurs sismiques 122 de l'outil sismique 118 à l'intérieur de l'ensemble de fond de trou 114. [0022] Les données sismiques échantillonnées peuvent ensuite être transférées vers une unité d'acquisition de données 128 située à la surface 104, par exemple, par l'intermédiaire d'une connexion câblée ou sans fil entre l'unité d'acquisition de données 128 et l'outil sismique 118 (ou un dispositif de télémétrie de celle-ci). Les données sismiques reçues par l'unité d'acquisition de données 128 depuis l'outil sismique 118 peuvent être stockées à l'intérieur d'un dispositif de stockage interne de l'unité d'acquisition de données 128. Dans certaines mises en oeuvre, l'outil sismique 118 peut également inclure une mémoire interne ou un dispositif de stockage dans lequel les données sismiques échantillonnées peuvent être stockées en fond de trou avant d'être transmises à l'unité d'acquisition de données 128. Il convient de noter que l'unité d'acquisition de données 128 peut être mise en oeuvre à l'aide de n'importe quel type de dispositif informatique ayant au moins un processeur et une mémoire. Un tel dispositif informatique peut également inclure une infrastructure de communication capable de recevoir et de transmettre des données sur un réseau, ainsi qu'une interface d'entrée/de sortie pour recevoir une entrée d'un utilisateur (par exemple, par l'intermédiaire d'une souris, d'un clavier QWERTY, d'un écran tactile, d'une tablette graphique ou d'un microphone) et présenter des informations à l'utilisateur, par exemple, par l'intermédiaire d'un écran ou d'un autre périphérique de sortie. Un tel dispositif informatique peut être, par exemple et sans limitation, un téléphone mobile, un assistant numérique personnel (PDA), une tablette, un ordinateur portable, un ordinateur de bureau, une station de travail, un serveur, un groupe de serveurs ou un autre type similaire de dispositif capable de traiter des instructions et de recevoir et de transmettre des données depuis et vers des personnes et d'autres dispositifs informatiques. 7 [0023] Dans un ou plusieurs modes de réalisation, les capteurs sismiques 122 peuvent échantillonner les ondes sismiques 124 en fond de trou à des intervalles de temps prédéterminés. En plus ou en variante, la source sismique 126 peut générer et transmettre les ondes sismiques 124 depuis la surface 104 à des intervalles prédéterminés. En conséquence, les données sismiques transmises à la surface 104 peuvent être acquises périodiquement par l'unité d'acquisition de données 128, par exemple, une fois par mois, une fois tous les trois mois, une fois tous les six mois, une fois par an, etc. Dans certaines mises en oeuvre, de telles données sismiques à intervalles de temps acquises par l'unité d'acquisition de données 128 peuvent être stockées localement et traitées afin d'analyser et d'estimer les propriétés de roches variables dans le temps d'une formation souterraine, comme cela sera décrit plus en détail ci-dessous. En variante, les données sismiques acquises par l'unité d'acquisition de données 128 peuvent être transmises par l'intermédiaire d'un réseau de communication à un dispositif informatique à distance, par exemple, le dispositif informatique 306 de la figure 3, comme cela sera décrit plus en détail ci-dessous. Le dispositif informatique à distance peut être utilisé pour mettre en oeuvre les techniques d'analyse et d'estimation sismiques de propriétés de roches variables dans le temps décrites ici. Dans certaines mises en oeuvre, l'unité d'acquisition de données 128 peut transférer les données sismiques par l'intermédiaire du réseau vers une base de données à distance (non représentée) pour le stockage. La base de données peut également être accessible au dispositif informatique à distance pour récupérer les données sismiques stockées par l'intermédiaire du réseau. [0024] Alors que seule l'unité d'acquisition de données 128 est représentée dans la figure 1, il convient de noter que l'environnement de puits de forage 100 peut inclure un nombre quelconque de ces unités ou dispositifs de traitement de données pour acquérir et traiter des informations provenant de l'outil sismique 118. De même, alors que seule la source sismique 126 est représentée, il convient de noter que n'importe quel nombre de sources sismiques peut être utilisé comme souhaité pour une mise en oeuvre particulière. [0025] La figure 2 est un diagramme schématique en coupe transversale d'un environnement marin pour l'acquisition de données sismiques à partir d'un relevé sismique marin d'une formation souterraine. Comme représenté dans la figure 2, un navire sismique 200 est positionné sur une surface 202 de l'océan. Le navire sismique 200 peut remorquer une ou plusieurs sources sismiques 204, telles qu'une source d'impulsion ou une source vibratoire. Les sources sismiques 204 peuvent transmettre des ondes sismiques 206 à travers le fond océanique 208. Les ondes sismiques 206 peuvent être réfléchies ou réfractées par des formations souterraines 210 situées sous le fond océanique 208 et reçues par un réseau de capteurs sismiques 212, tels que des hydrophones, traînant derrière le navire sismique 200 sur une ou plusieurs flûtes 214. Dans certains aspects, les 8 flûtes 214 peuvent inclure un câblage électrique ou à fibres optiques pour connecter le réseau de capteurs 212 à un équipement sismique sur le navire sismique 200, y compris une unité d'acquisition de données 216. Les capteurs 212 peuvent mesurer les réflexions des ondes sismiques 124 et transmettre les mesures à travers les flûtes 214 pour le stockage dans l'unité d'acquisition de données 216. [0026] La figure 3 est un schéma fonctionnel d'un système 300 d'analyse et d'estimation sismiques de propriétés de roches variables dans le temps d'une formation souterraine. À des fins de discussion et d'explication, le système 300 sera décrit par référence aux composants de l'environnement de puits de forage 100 de la figure 1, comme décrit ci-dessus. Par exemple, le système 300 tel que représenté dans la figure 3 inclut l'unité d'acquisition de données 128 de la figure 1. Cependant, il convient de noter que les modes de réalisation de la présente invention ne sont pas destinés à se limiter à ceux-ci et que le système 300 peut être adapté pour inclure d'autres composants comme souhaité pour une mise en oeuvre particulière. Par exemple, le système 300 peut inclure l'unité d'acquisition de données 216 et un outil sismique incluant le réseau de capteurs sismiques 212 de la figure 2 pour mettre en oeuvre les techniques d'analyse et d'estimation sismiques décrites ici dans un environnement de puits de forage marin. [0027] Comme décrit ci-dessus, l'unité d'acquisition de données 128 peut être couplée en communication à l'outil sismique 118 de la figure 1 par l'intermédiaire d'une connexion câblée ou sans fil. L'outil sismique 118 peut inclure un ou plusieurs des capteurs sismiques 122 pour détecter et échantillonner des ondes d'énergie sismique générées par une source sismique (par exemple, la source sismique 126 de la figure 1) et réfléchies ou réfractées par des couches de formation souterraine adjacentes à un puits de forage (par exemple, les couches de formation 112 adjacentes au puits de forage 108 de la figure 1). Les informations sismiques échantillonnées peuvent être transmises par l'outil sismique 118 sur la connexion câblée ou sans fil à l'unité d'acquisition de données 128 à la surface. L'unité d'acquisition de données 128 peut stocker les informations sismiques qu'elle reçoit depuis l'outil sismique 118 à l'intérieur d'un dispositif de stockage 302 en tant que données sismiques 304. Le dispositif de stockage 302 peut être n'importe quel type de support d'enregistrement couplé à un circuit intégré qui commande l'accès au support d'enregistrement. Le support d'enregistrement peut être, par exemple et sans limitation, une mémoire à semi-conducteur, un disque dur, ou un autre type similaire de mémoire ou de dispositif de stockage. [0028] Dans un ou plusieurs modes de réalisation, les données sismiques 304 peuvent être transmises de l'unité d'acquisition de données 128 à un dispositif informatique 306 par l'intermédiaire d'un réseau 308. Le réseau 308 peut être tout type de réseau ou de combinaison de 9 réseaux utilisé(e) pour communiquer des informations entre différents dispositifs informatiques. Le réseau 308 peut inclure, mais sans s'y limiter, un réseau câblé (par exemple, Ethernet) ou sans fil (par exemple, Wi-Fi et 3G). De plus, le réseau 308 peut inclure, mais sans s'y limiter, un réseau local, un réseau moyen et/ou un réseau étendu tel qu'Internet. L'unité d'acquisition de données 128 et le dispositif informatique 306 peuvent être couplés à, ou inclure, des dispositifs de communication respectifs 310A, 310B. Les dispositifs de communication 310A et 310B incluent ou sont couplés aux antennes 312A et 312B, respectivement, pour émettre et recevoir des informations par l'intermédiaire du réseau 308. Cependant, il convient de noter que les modes de réalisation décrits ne sont pas limités à ceux-ci et que tout moyen approprié peut être utilisé pour transmettre des informations entre l'unité d'acquisition de données 128 et le dispositif informatique 306. Dans un ou plusieurs modes de réalisation, les données sismiques 304 peuvent inclure des informations brutes provenant des capteurs 122 de l'outil sismique 118, qui peuvent avoir besoin d'être traitées et formatées de manière appropriée par l'unité d'acquisition de données 128 avant la transmission des données sismiques 304 au dispositif informatique 306. [0029] Le dispositif informatique 306 peut inclure un dispositif de traitement 314, un bus 316 et un dispositif de mémoire 318. Le dispositif de traitement 314 peut exécuter une pluralité de fonctions ou d'opérations pour estimer une valeur d'au moins une propriété de roche variable dans le temps de la formation souterraine à un moment futur sur la base des données sismiques 304 reçues de l'unité d'acquisition de données 128 et stockées à l'intérieur du dispositif de mémoire 318. Le dispositif de traitement 314 peut exécuter de telles fonctions ou opérations à la suite de l'exécution d'instructions 320 stockées dans le dispositif de mémoire 318. Le dispositif de traitement 314 peut être un processeur et le dispositif de mémoire 318 peut être n'importe quel type de support de stockage lisible par un processeur, y compris, mais sans s'y limiter, une mémoire à semi-conducteur, un disque dur ou un autre type similaire de mémoire ou de dispositif de stockage. Les instructions 320 peuvent inclure des instructions spécifiques au processeur générées par un compilateur ou un interpréteur à partir de code écrit dans tout langage de programmation informatique approprié. Les instructions 320 peuvent inclure des instructions exécutables par le dispositif de traitement 314 pour effectuer diverses opérations relatives aux techniques d'analyse et de modélisation sismique décrites ici. De telles opérations peuvent inclure, par exemple, des opérations d'estimation d'une valeur de la propriété de roche variable dans le temps à un instant ou à un intervalle de temps futurs sur la base de divers schémas d'approximation qui tiennent compte d'une ou de plusieurs valeurs observées de la propriété de roche depuis un intervalle de temps antérieur, comme cela sera décrit plus en détail ci-dessous. [0030] Dans un ou plusieurs modes de réalisation, les fonctions d'estimation de propriétés de roches variables dans le temps susmentionnées, exécutées par le dispositif de traitement 314, peuvent inclure des fonctions pour analyser les données sismiques 304 et déterminer des tendances de données existantes. Une visualisation de telles tendances de données, par exemple, sous la forme d'un graphique linéaire ou d'un autre type de visualisation graphique, peut également être affichée pour un utilisateur par l'intermédiaire d'une unité d'affichage 322. Comme décrit ci-dessus, les données sismiques 304 peuvent inclure des mesures d'ondes sismiques réfléchies et réfractées acquises par l'outil sismique 118 et l'unité d'acquisition de données 128 sur une série d'intervalles de temps. Ces mesures peuvent être acquises dans le cadre d'un relevé sismique bidimensionnel (2D) ou tridimensionnel (3D) de la formation, qui est périodiquement effectué à différents moments, produisant ainsi un certain nombre d'études de données sismiques. [0031] Dans un ou plusieurs modes de réalisation, l'estimation de valeurs de propriétés de roches variables dans le temps peut être effectuée pour un intervalle de temps futur sur la base des mesures sismiques acquises pour au moins deux intervalles de temps précédents, c'est-à-dire, lorsqu'au moins deux études sismiques sont disponibles. Le tableau suivant représente un exemple des mesures sismiques acquises pour une propriété de roche variable dans le temps d'une formation. -ème Intervalle de temps pour la roche variable dans temps i mesure 0 to 1 t1 2 t2 Q tQ [0032] Comme représenté dans le tableau ci-dessus, des mesures sismiques de la propriété de roche variable dans le temps peuvent être acquises sur une pluralité d'intervalles de temps t,, ti+i, t,+2, tQ, où un indice «I» représente une base de temps pour des mesures de la propriété de roche à chaque intervalle et la valeur de i est un entier compris entre 0 et un certain nombre Q. L'état de la propriété de roche variable dans le temps à un instant donné « t » peut dépendre de la fréquence à laquelle des relevés sismiques de la formation sont effectués et à laquelle des mesures sismiques des relevés sont acquises. La fréquence d'acquisition sismique dans cet exemple peut être basée sur la taille ou la durée de chaque intervalle de temps t,, par exemple, mesuré en jours, mois, années, etc. 11 [0033] Dans un ou plusieurs modes de réalisation, des relevés sismiques de la formation souterraine peuvent être effectués sur une pluralité d'intervalles de temps afin de surveiller les changements dans une ou plusieurs propriétés de roches variables dans le temps de la formation. De tels changements peuvent être dus, par exemple, à l'élimination de pétrole et/ou de gaz de la formation pendant les opérations de production d'hydrocarbures. Comme cela sera décrit plus en détail ci-dessous, les études sismiques ou les ensembles de données à intervalles de temps obtenus à partir des relevés sismiques répétés peuvent être analysés quantitativement et qualitativement à l'aide d'un ou de plusieurs schémas d'approximation afin d'évaluer dans quelle mesure les propriétés de roches variables dans le temps peuvent varier dans le temps en conséquence de déplacements de fluide et de changements de propriétés de roches élastiques à l'intérieur de la formation. Un exemple des changements qui peuvent survenir pour une propriété de roche d'une formation souterraine dans le temps et dans l'espace est représenté dans les figures 4A à 4C. [0034] Les figures 4A, 4B et 4C représentent différents instants ou états d'un profil illustratif d'une propriété de roche variable dans le temps sur des intervalles de temps successifs to, ti et t2, respectivement. La propriété de roche variable dans le temps peut être, par exemple et sans limitation, l'impédance de la roche à l'intérieur d'une zone de la formation, par exemple, une zone ou une surface de formation qui a été ciblée pour l'exploration et/ou la production d'hydrocarbures. Le profil de la propriété de roche à chaque intervalle de temps peut être une image de la propriété de roche variable dans le temps sous la forme d'une carte 2D spatialement variable de la propriété de roche. Une telle carte peut représenter, par exemple, une distribution spatiale de la propriété de roche variable dans le temps à l'intérieur d'une zone sismique 2D de la formation souterraine. Cependant, il convient de noter que les techniques décrites ne sont pas limitées à celles-ci et que ces techniques peuvent également être appliquées à des volumes et à des ensembles de données sismiques 3D. [0035] Comme représenté dans les figures 4A à 4C, la propriété d'impédance variable dans le temps de la roche varie de valeurs élevées à des valeurs inférieures lorsque la distribution de la roche variable dans le temps (zone ombrée) augmente de la partie gauche de la carte à la partie droite au fil du temps. L'augmentation peut être due à un pourcentage plus élevé de la roche variable dans le temps déplacée, par exemple, à mesure que davantage d'hydrocarbures sont produits à partir de la formation d'un intervalle de temps à l'autre. [0036] Dans un ou plusieurs modes de réalisation, des profils de propriétés de roches variables dans le temps pour des intervalles de temps ou des instants futurs après l'heure actuelle ou indiquée peuvent être estimés sur la base des propriétés de roches variables dans le temps mesurées ou calculées à partir des intervalles de temps actuels et précédents. Comme décrit ci-dessus, les 12 valeurs d'une ou de plusieurs propriétés de roches variables dans le temps peuvent être calculées sur la base de données sismiques acquises pour une zone ciblée de la formation à chacun d'une pluralité d'intervalles de temps. Le calcul peut inclure une inversion des données sismiques acquises pour générer un ou plusieurs ensembles de données d'inversion sismique représentant les propriétés de roches sur la pluralité d'intervalles de temps et le calcul des valeurs de propriétés de roches sur la base des ensembles de données d'inversion sismique générés. [0037] Dans un ou plusieurs modes de réalisation, les valeurs calculées peuvent être analysées pour déterminer des tendances dans les données existantes sur la pluralité d'intervalles de temps. L'analyse peut impliquer, par exemple, l'utilisation des ensembles de données d'inversion sismique pour simuler l'écoulement de fluide à l'intérieur du support de propriétés de roches souterraines sur la pluralité d'intervalles de temps. Les résultats de la simulation peuvent ensuite être utilisés pour générer un modèle pétro-élastique de la formation souterraine pour prédire les propriétés de roches variables dans le temps à un intervalle de temps futur, par exemple, le long d'une série en temps avancé. Dans un ou plusieurs modes de réalisation, le modèle pétro-élastique peut être défini par une grille stratigraphique 2D ou 3D de cellules dans laquelle les propriétés de roches de la formation peuvent être réparties. Par exemple, une épaisseur verticale de chaque cellule dans la grille peut représenter des valeurs de la propriété de roche dans le domaine temporel et une largeur horizontale de chaque cellule dans la grille peut représenter des valeurs de la propriété de roche dans un domaine spatial. [0038] À des fins de continuité entre l'espace et le temps, une ou plusieurs techniques de modélisation numérique peuvent être utilisées pour effectuer les simulations d'écoulement de manière à ce qu'un domaine de dépendance physique existe à l'intérieur du domaine de dépendance numérique. Cela peut fournir une synergie entre les ensembles de données sismiques inversés et un modèle numérique utilisé pour la simulation d'écoulement. Cela peut également permettre de co-analyser et d'étalonner l'écoulement de fluide de simulation numérique à l'intérieur de la zone de formation ciblée avec les données sismiques acquises sur les mêmes intervalles de temps. Pour garantir que la condition de dépendance ci-dessus entre les domaines physique et numérique est respectée, la précision de la prédiction de série en temps avancé peut être vérifiée à l'aide de l'équation (1) : z lu At' REVENDICATIONS1. Procédé de prédiction de propriétés de roches variables dans le temps mise en oeuvre par ordinateur et de forage d'un puits de forage, le procédé comprenant : l'acquisition de données sismiques pour une formation souterraine sur une pluralité d'intervalles de temps ; le calcul d'une valeur d'au moins une propriété de roche de la formation souterraine pour chacun de la pluralité d'intervalles de temps, sur la base des données sismiques correspondantes acquises pour cet intervalle de temps ; la détermination d'au moins l'une d'une tendance ou d'une relation spatio-temporelle dans les données sismiques sur la base de la valeur de l'au moins une propriété de roche calculée pour chacun de la pluralité d'intervalles de temps ; l'estimation d'une valeur de l'au moins une propriété de roche pour un intervalle de temps futur, sur la base de la détermination ; la sélection d'un emplacement pour un puits de forage à forer à l'intérieur de la formation souterraine, sur la base de la valeur estimée de l'au moins une propriété de roche ; et le forage du puits de forage à l'emplacement sélectionné. 2. Procédé selon la revendication 1, dans lequel le calcul comprend : la mise en oeuvre d'une inversion des données sismiques acquises pour générer un ou plusieurs ensembles de données d'inversion sismique représentant l'au moins une propriété de roche sur la pluralité d'intervalles de temps ; et le calcul des valeurs de l'au moins une propriété de roche, sur la base des ensembles de données d'inversion sismique générés. 3. Procédé selon la revendication 2, dans lequel la détermination comprend : la simulation de l'écoulement de fluide à l'intérieur de la formation souterraine sur la pluralité d'intervalles de temps, sur la base des un ou plusieurs ensembles de données d'inversion sismique ; la génération d'un modèle pétro-élastique de la formation souterraine, sur la base de la simulation ; et la détermination d'au moins l'une de la tendance ou de la relation spatio-temporelle dans les valeurs de l'au moins une propriété de roche, sur la base du modèle pétro-élastique généré. 30 4. Procédé selon la revendication 3, dans lequel l'estimation comprend : la génération d'un ensemble de réalisations équiprobables de l'au moins une propriété de roche pour l'intervalle de temps futur ; la sélection d'au moins une réalisation parmi l'ensemble de réalisations équiprobables, sur la base du modèle pétro-élastique de la formation souterraine ; et l'estimation de la valeur de l'au moins une propriété de roche pour un intervalle de temps futur, sur la base de la réalisation sélectionnée. 5. Procédé selon la revendication 4, dans lequel la génération de l'ensemble de réalisations équiprobables comprend en outre : l'ajout d'un bruit blanc gaussien aux réalisations équiprobables afin de tenir compte d'une ou de plusieurs sources d'incertitude. 6. Procédé selon la revendication 4, dans lequel le modèle pétro-élastique de la formation souterraine est défini par une grille tridimensionnelle (3D) de cellules, et la valeur de l'au moins une propriété de roche est estimée pour chaque cellule de la grille 3D de cellules correspondant à un volume 3D de la formation souterraine. 7. Procédé selon la revendication 3, comprenant en outre : l'obtention de mesures d'une valeur réelle de l'au moins une propriété de roche à l'intervalle de temps futur ; la comparaison de la valeur estimée de l'au moins une propriété de roche avec la valeur réelle ; et la mise à jour du modèle pétro-élastique de la formation souterraine, sur la base de la comparaison ; éventuellement, dans lequel la relation spatio-temporelle est déterminée à l'aide d'un variogramme spatio-temporel basé sur une fonction aléatoire représentant des changements de valeurs de l'au moins une propriété de roche pour un ou plusieurs emplacements dans l'espace et le temps ; et éventuellement, dans lequel la comparaison comprend le calcul d'une covariance entre la valeur estimée de l'au moins une propriété de roche et la valeur réelle pour chacun des un ou plusieurs emplacements dans l'espace et le temps. 8. Système comprenant : au moins un processeur ; et 31 une mémoire couplée au processeur dans laquelle sont stockées des instructions qui, lorsqu'elles sont exécutées par le processeur, amènent le processeur à exécuter une pluralité de fonctions, y compris des fonctions pour : acquérir des données sismiques pour une formation souterraine sur une pluralité d'intervalles de temps ; calculer une valeur d'au moins une propriété de roche de la formation souterraine pour chacun de la pluralité d'intervalles de temps, sur la base des données sismiques correspondantes acquises pour cet intervalle de temps ; déterminer au moins l'une d'une tendance ou d'une relation spatio-temporelle dans les données sismiques sur la base de la valeur de l'au moins une propriété de roche calculée pour chacun de la pluralité d'intervalles de temps ; estimer une valeur de l'au moins une propriété de roche pour un intervalle de temps futur, sur la base de la détermination ; et sélectionner un emplacement pour un puits de forage à forer à l'intérieur de la formation souterraine, sur la base de la valeur estimée de l'au moins une propriété de roche, dans lequel le puits de forage est foré à l'emplacement sélectionné. 9. Système selon la revendication 8, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour : exécuter une inversion des données sismiques acquises pour générer un ou plusieurs ensembles de données d'inversion sismique représentant l'au moins une propriété de roche sur la pluralité d'intervalles de temps ; et calculer les valeurs de l'au moins une propriété de roche, sur la base des ensembles de données d'inversion sismique générés. 10. Système selon la revendication 9, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour : simuler un écoulement de fluide à l'intérieur de la formation souterraine sur la pluralité d'intervalles de temps, sur la base des un ou plusieurs ensembles de données d'inversion sismique ; générer un modèle pétro-élastique de la formation souterraine, sur la base de la simulation ; et déterminer au moins l'une de la tendance ou de la relation spatio-temporelle dans les valeurs de l'au moins une propriété de roche, sur la base du modèle pétro-élastique généré. 32 11. Système selon la revendication 10, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour : générer un ensemble de réalisations équiprobables de l'au moins une propriété de roche pour l'intervalle de temps futur ; sélectionner au moins une réalisation parmi l'ensemble de réalisations équiprobables, sur la base du modèle pétro-élastique de la formation souterraine ; et estimer la valeur de l'au moins une propriété de roche pour l'intervalle de temps futur, sur la base de la réalisation sélectionnée. 12. Système selon la revendication 11, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour ajouter un bruit blanc gaussien aux réalisations équiprobables afin de tenir compte d'une ou de plusieurs sources d'incertitude. 13. Système selon la revendication 11, dans lequel le modèle pétro-élastique de la formation souterraine est défini par une grille tridimensionnelle (3D) de cellules, et la valeur de l'au moins une propriété de roche est estimée pour chaque cellule de la grille 3D de cellules correspondant à un volume 3D de la formation souterraine. 14. Système selon la revendication 10, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour : obtenir des mesures d'une valeur réelle de l'au moins une propriété de roche à l'intervalle de temps futur ; comparer la valeur estimée de l'au moins une propriété de roche avec la valeur réelle ; et mettre à jour le modèle pétro-élastique de la formation souterraine, sur la base de la comparaison ; éventuellement, dans lequel la relation spatio-temporelle est déterminée à l'aide d'un variogramme spatio-temporel basé sur une fonction aléatoire représentant des changements de valeurs de l'au moins une propriété de roche pour un ou plusieurs emplacements dans l'espace et le temps ; et éventuellement, dans lequel les fonctions exécutées par le processeur incluent des fonctions pour calculer une covariance entre la valeur estimée de l'au moins une propriété de roche et la valeur réelle pour chacun des un ou plusieurs emplacements dans l'espace et le temps. 33 1 5 . Support de stockage lisible par ordinateur dans lequel sont stockées des instructions qui, lorsqu'elles sont exécutées par un processeur, amènent le processeur à exécuter une pluralité de fonctions, y compris des fonctions pour : acquérir des données sismiques pour une formation souterraine sur une pluralité d'intervalles de temps ; calculer une valeur d'au moins une propriété de roche de la formation souterraine pour chacun de la pluralité d'intervalles de temps, sur la base des données sismiques correspondantes acquises pour cet intervalle de temps ; déterminer au moins l'une d'une tendance ou d'une relation spatio-temporelle dans les données sismiques sur la base de la valeur de l'au moins une propriété de roche calculée pour chacun de la pluralité d'intervalles de temps ; estimer une valeur de l'au moins une propriété de roche pour un intervalle de temps futur, sur la base de la détermination ; et sélectionner un emplacement pour un puits de forage à forer à l'intérieur de la formation souterraine, sur la base de la valeur estimée de l'au moins une propriété de roche, dans laquelle le puits de forage est foré à l'emplacement sélectionné.