FR3032532A1 - - Google Patents

Abstract

Procédé pour modéliser une formation souterraine d'un champ. Le procédé comprend : obtenir un volume sismique comprenant une pluralité de traces sismiques de la formation souterraine du champ, calculer, sur la base du volume sismique, une valeur déduite sismiquement d'un attribut structurel représentant une caractéristique structurelle de la formation souterraine, calculer, sur la base d'un modèle structurel, une valeur déduite structurellement de l'attribut structurel, le modèle structurel comprenant une pluralité de couches structurelles de la formation souterraine, comparer la valeur déduite sismiquement et la valeur déduite structurellement pour générer une valeur de différence représentant un écart de modélisation de l'attribut structurel à une position correspondante dans la formation souterraine, et générer un résultat d'interprétation sismique sur la base de la valeur de différence et de la position correspondante.

Description

1 CONTEXTE [0001] Certaines opérations, telles que la prospection, le forage, le test au câble, les complétions, la production, la planification et l'analyse de champs, peuvent être effectuées pour localiser et de recueillir des fluides précieux de fond de puits. Des prospections sont souvent effectuées en utilisant des méthodologies d'acquisition, telles que des analyseurs ou des dispositifs de prospection sismiques pour générer des cartes de formations souterraines. Ces formations sont souvent analysées pour déterminer la présence de gisements souterrains, tels que des fluides ou des minéraux précieux, ou pour déterminer si les formations présentent des caractéristiques appropriées pour le stockage de fluides. Les gisements souterrains ne sont pas limités à un hydrocarbure tel que le pétrole, et dans l'ensemble du présent document, les termes de "champ pétrolifère" et d' "opération de champ pétrolifère" peuvent être utilisés indifféremment avec les termes de "champ" et d' "opération de champ" pour désigner un champ comportant des types quelconques de fluides ou de minéraux précieux et des opérations de champ concernant n'importe lesquels de ces gisements souterrains.

RESUME [0002] De manière générale, selon l'un de ses aspects, l'invention concerne un procédé pour modéliser une formation souterraine d'un champ. Le procédé comprend : obtenir un volume sismique comprenant une pluralité de traces sismiques de la formation souterraine du champ, calculer, sur la base du volume sismique, une valeur déduite sismiquement d'un attribut structurel représentant une caractéristique structurelle de la formation souterraine, calculer, sur la base d'un modèle structurel, une valeur déduite structurellement de l'attribut structurel, le modèle structurel comprenant une pluralité de couches structurelles de la formation souterraine, comparer la valeur déduite sismiquement et la valeur déduite structurellement pour générer une valeur de différence représentant un écart de modélisation de l'attribut structurel à une position correspondante dans la formation souterraine, et générer un résultat d'interprétation sismique sur la base de la valeur de différence et de la position correspondante. [0003] Selon un premier mode de réalisation, le procédé pour modéliser une formation souterraine d'un champ comprend en outre: afficher, dans une fenêtre graphique sur un écran d'affichage physique, la valeur de différence et la position correspondante ; et effectuer, en réponse à l'affichage, une interprétation sismique du volume sismique pour générer un résultat d'interprétation sismique, dans lequel l'interprétation sismique est ajustée sur la base au moins 3032532 2 de la valeur de différence et la position correspondante. Dans ce premier mode de réalisation, le procédé comprend en outre : recevoir, en réponse à la génération du résultat d'interprétation sismique, un ajustement du modèle structurel pour générer le modèle structurel révisé, dans lequel l'ajustement est fondé sur la valeur de différence et la position 5 correspondante ; et effectuer une opération de champ sur la base du modèle structurel révisé. Dans ce premier mode de réalisation, le procédé comprend en outre : calculer, sur la base du modèle structurel révisé, une valeur déduite structurellement ajustée de l'attribut structurel ; comparer la valeur déduite sismiquement et la valeur déduite structurellement ajustée pour générer une valeur de différence ajustée représentant un écart de modélisation de l'attribut 10 structurel sur la base du modèle structurel révisé ; et afficher, dans la fenêtre graphique sur l'écran d'affichage physique, la valeur de différence ajustée, dans lequel l'interprétation sismique fait partie d'une séquence d'interprétations sismiques qui sont réalisées de manière itérative jusqu'à ce que la valeur de différence ajustée respecte un critère prédéterminé. [0004] Selon un deuxième mode de réalisation du procédé pour modéliser une 15 formation souterraine d'un champ, le volume sismique comprend une pluralité de cellules sismiques tridimensionnelles (3D) qui sont associées à une pluralité de valeurs déduites sismiquement de l'attribut structurel, dans lequel la pluralité de valeurs déduites sismiquement comprennent la valeur déduite sismiquement ; le modèle structurel comprend une pluralité de cellules structurelles 3D qui sont associées à une pluralité de valeurs déduites structurellement 20 de l'attribut structurel, dans lequel la pluralité de valeurs déduites structurellement comprennent la valeur déduite structurellement ; la valeur de différence fait partie d'une pluralité de valeurs de différence qui sont générées en comparant la pluralité de valeurs déduites sismiquement et la pluralité de valeurs déduites structurellement ; la comparaison de la pluralité de valeurs déduites sismiquement et de la pluralité de valeurs déduites 25 structurellement comprend une mise en correspondance de la pluralité de cellules sismiques et de la pluralité de cellules structurelles pour générer une pluralité de cellules rééchantillonnées, la pluralité de valeurs de différence sont associées à la pluralité de cellules rééchantillonnées ; le procédé comprenant en outre un affichage, dans une fenêtre graphique sur l'écran d'affichage physique, la pluralité de valeurs de différence relativement à la pluralité de 30 cellules rééchantillonnées. Selon ce deuxième mode de réalisation, le modèle structurel comprend en outre une pluralité de failles, le procédé comprenant en outre : déterminer une pluralité de zones de voisinage de failles à au plus une distance prédéterminée de la pluralité de failles ; et qualifier un sous-ensemble de la pluralité de valeurs de différence sur la base 3032532 3 d'au moins un élément sélectionné dans un groupe constitué (i) d'un sous-ensemble correspondant de la pluralité de cellules rééchantillonnées se situant en dehors de la pluralité de zones de voisinage de failles et (ii) d'une mesure de qualité sismique du sous-ensemble correspondant de la pluralité de cellules rééchantillonnées, l'affichage étant effectué en 5 réponse à la qualification. BREVE DESCRIPTION DES DESSINS [0005] Les dessins annexés illustrent plusieurs modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume et ne doivent pas être considérés comme limitant son étendue, étant entendu que des attributs 10 sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être utilisés dans d'autres modes de réalisation tout aussi efficaces. [0006] La figure 1 est une vue schématique, partiellement en coupe transversale, d'un champ comportant une pluralité d'outils d'acquisition de données positionnés en diverses positions le long du champ pour recueillir des données en provenance de la formation 15 souterraine, dans lequel des modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre. [0007] La figure 2 représente un exemple de système dans lequel un ou plusieurs modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre. 20 [0008] La figure 3 représente un exemple de procédé destiné à des attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume conformément à un ou plusieurs modes de réalisation. [0009] Les figures 4 et 5 représentent des exemples de copies d'écrans représentant des attributs déduits du signal sismique et des attributs déduits de la propriété d'âge 25 géologique relatif d'un modèle basé sur le volume conformément à un ou plusieurs modes de réalisation. [0010] La figure 6 représente un système informatique dans lequel un ou plusieurs modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre.

3032532 4 DESCRIPTION DETAILLEE [0011] Des modes de réalisation sont représentés sur les dessins identifiés ci-dessus et sont décrits ci-dessous. Dans la description des modes de réalisation, des références numériques semblables ou identiques sont utilisées pour identifier des éléments communs ou 5 semblables. Les dessins ne sont pas nécessairement à l'échelle et certains détails et certaines vues des dessins peuvent être représentés selon une échelle exagérée ou de manière schématique pour plus de clarté et de concision. [0012] Les modes de réalisation décrits ici comprennent un processus destiné à calculer des attributs structurels de la formation souterraine. Des valeurs déduites 10 sismiquement et des valeurs déduites structurellement des attributs structurels sont comparées pour générer des différences de comparaison. Dans un ou plusieurs modes de réalisation, les valeurs déduites structurellement sont déduites à partir de la propriété d'âge géologique relatif d'un modèle structurel. Les différences de comparaison correspondent aux différences entre le modèle basé sur le volume et les données sismiques dont a été déduit le modèle basé sur le 15 volume. Par conséquent, l'analyse des données sismiques et/ou le modèle basé sur le volume sont ajustés pour minimiser les différences et d'améliorer la précision des attributs sismiques et des attributs structurels calculés. [0013] La figure 1 illustre une vue schématique, partiellement en coupe transversale, d'un champ (100) dans lequel un ou plusieurs modes de réalisation d'attributs sismiques 20 déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre. Dans un ou plusieurs modes de réalisation, un ou plusieurs des modules et des éléments représentés sur la figure 1 peuvent être omis, répétés et/ou substitués. Par conséquent, les modes de réalisation de contrôle de qualité d'une auto-poursuite d'horizon tridimensionnel (3D) dans un volume sismique ne doivent pas être considérés comme étant 25 limités aux agencements de modules spécifiques représentés sur la figure 1. [0014] Comme illustré sur la figure 1, le champ (100) comporte la formation souterraine (104), des outils d'acquisition de données (102-1), (102-2), (102-3) et (102-4), un système d'emplacement de forage A (204-1), un système d'emplacement de forage B (204-2), un système d'emplacement de forage C (204-3), une unité de surface (202), et un système 30 (208) informatique d'exploration et de production (E&P). La formation souterraine (104) comporte plusieurs structures géologiques, telles qu'une couche de grès (106-1), une couche de calcaire (106-2), une couche de schiste (106-3), une couche de sable (106-4), et une ligne de faille (107). 3032532 5 [0015] Dans un ou plusieurs modes de réalisation, des outils d'acquisition de données (102-1), (102-2), (102-3) et (102-4) sont positionnés en diverses positions le long du champ (100) pour recueillir des données de la formation souterraine (104), cela étant désigné sous le nom d'opérations de prospection. Ces outils d'acquisition de données sont notamment aptes à 5 mesurer la formation souterraine (104) et à détecter les caractéristiques des structures géologiques de la formation souterraine (104). A titre d'exemple, des tracés de données (1081), (108-2), (108-3) et (108-4) sont illustrés le long du champ (100) pour mettre en évidence les données générées par ces outils d'acquisition de données. Plus précisément, le tracé de données statique (108-1) est un temps de réponse sismique bidirectionnel. Le tracé statique 10 (108-2) représente des données de carottage (mesurées à partir d'une carotte de la formation (104)) ou d'autres données de diagraphie pétrophysiques. Le tracé de données statique (108-3) est une trace de diagraphie, également nommée diagraphie de forage. La courbe ou le graphe de décroissance de la production (108-4) est un tracé de données dynamique du débit de fluide au cours du temps. D'autres données peuvent également être recueillies, telles que des 15 données historiques, des entrées d'utilisateurs, des informations économiques, et/ou d'autres données de mesure et d'autres paramètres d'intérêt. [0016] En vue d'acquérir le temps de réponse sismique bidirectionnel sur le tracé de données statique (108-1), les outils d'acquisition de données (102-1) peuvent être un véhicule de prospection qui est apte à mesurer des propriétés de la formation souterraine sur la base de 20 vibrations sonores. Une vibration sonore de ce type (par exemple 186, 188, 190) générée par une source (170) se réfléchit sur une pluralité d'horizons (par exemple 172, 174, 176) dans la formation souterraine (104). Chacune des vibrations sonores (par exemple 186, 188, 190) est reçue par un ou plusieurs capteurs (par exemple 180, 182, 184), tels que des récepteurs à géophones, situés sur la surface de la Terre. Les géophones produisent des signaux de sortie 25 électriques, qui peuvent par exemple être transmis en tant que données d'entrée à un ordinateur (192) sur le véhicule de prospection (102-1). En réponse aux données d'entrée, l'ordinateur (192) peut générer une sortie de données sismiques, telle que le temps de réponse sismique bidirectionnel sous la forme d'une trace sismique. [0017] Comme il est par ailleurs illustré sur la figure 1, chacun du système 30 d'emplacement de forage A (204-1), du système d'emplacement de forage B (204-2) et du système d'emplacement de forage C (204-3) est associé à une installation de forage, un puits de forage, et d'autres équipements d'emplacement de forage configurés en vue d'effectuer des opérations de puits de forage, telles que la diagraphie, le forage, la fracturation, la production, 3032532 6 ou d'autres opérations appropriées. A titre d'exemple, les systèmes d'emplacements de forage (204-1), (204-2), (204-3) sont associés à une installation de forage (101), à un puits de forage (103), et à des équipements de forage en vue d'effectuer une opération de forage. De même, le système d'emplacement de forage B (204-2) et le système d'emplacement de forage C (204-3) 5 sont associés à des installations de forage, à des puits de forage, et à d'autres équipements d'emplacement de forage respectifs, tels que des équipements de production et des équipements de diagraphie en vue d'effectuer respectivement l'opération de production et l'opération de diagraphie. De manière générale, les opérations de prospection et les opérations de puits de forage sont désignées sous le nom d'opérations de champ concernant le champ 10 (100). De plus, les outils d'acquisition de données et les équipements d'emplacement de forage sont désignés sous le nom d'équipements d'opération de champ. Ces opérations de champ peuvent être effectuées de manière commandée par une unité de surface (202). A titre d'exemple, l'équipement d'opération de champ peut être commandé par un signal de commande d'opération de champ envoyé par l'unité de surface (202). 15 [0018] Dans un ou plusieurs modes de réalisation, l'unité de surface (202) est fonctionnellement reliée aux outils d'acquisition de données (102-1), (102-2), (102-3), (1024), et/ou aux systèmes d'emplacements de forage (204-1), (204-2), (204-3). En particulier, l'unité de surface (202) est configurée en vue d'envoyer des commandes aux outils d'acquisition de données (102-1), (102-2), (102-3), (102-4), et/ou aux systèmes 20 d'emplacements de forage (204-1), (204-2), (204-3) et d'en recevoir des données. Dans un ou plusieurs modes de réalisation, l'unité de surface (202) peut être située au niveau des systèmes d'emplacements de forage (204-1), (204-2), (204-3) et/ou en des positions distantes. L'unité de surface (202) peut être pourvue d'installations informatiques destinées à recevoir, stocker, traiter et/ou analyser les données provenant des outils d'acquisition de données (102-1), (102- 25 2), (102-3), (102-4), des systèmes d'emplacements de forage (204-1), (204-2), (204-3), et/ou d'une autre partie du champ (100). L'unité de surface (202) peut également être pourvue de mécanismes d'actionnement ou être fonctionnellement prévue pour lesdits mécanismes au niveau du champ (100). L'unité de surface (202) peut ensuite envoyer des signaux de commande au champ (100) en réponse aux données reçues, par exemple en vue de 30 commander et/ou d'optimiser diverses opérations de champ décrites ci-dessus. [0019] Dans un ou plusieurs modes de réalisation, l'unité de surface (202) est reliée de manière à communiquer avec un système informatique E&P (208). Dans un ou plusieurs modes de réalisation, les données reçues par l'unité de surface (202) peuvent être envoyées au 3032532 7 système informatique E&P (208) à des fins d'analyse plus poussée. De manière générale, le système informatique E&P (208) est configuré en vue d'analyser, de modéliser, de commander, d'optimiser ou d'effectuer des tâches de gestion des opérations de champ mentionnées ci-dessus sur la base des données fournies par l'unité de surface (202). Dans un 5 ou plusieurs modes de réalisation, le système informatique E&P (208) est doté d'une fonctionnalité destinée à manipuler et analyser les données, comme le fait d'effectuer une interprétation sismique ou une interprétation de diagraphies d'images de résistivité de puits de forage pour identifier des surfaces géologiques dans la formation souterraine (104) ou effectuer une simulation, une planification et une optimisation des opérations de production 10 des systèmes d'emplacements de forage (204-1), (204-2), (204-3). Dans un ou plusieurs modes de réalisation, le résultat généré par le système informatique E&P (208) peut être affiché afin d'être visualisé par l'utilisateur au moyen d'un écran d'affichage bidimensionnel (2D), d'un écran d'affichage tridimensionnel (3D), ou d'autres écrans d'affichage appropriés. Bien que l'unité de surface (202) soit représentée sur la figure 1 comme étant séparée du 15 système informatique E&P (208), dans d'autres exemples, l'unité de surface (202) et le système informatique E&P (208) peuvent également être combinés. [0020] La figure 2 représente de manière plus détaillée le système informatique E&P (208) dans lequel un ou plusieurs modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre.

20 Dans un ou plusieurs modes de réalisation, un ou plusieurs des modules et des éléments représentés sur la figure 2 peuvent être omis, répétés et/ou substitués. Par conséquent, les modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume ne doivent pas être considérés comme étant limités aux agencements de modules spécifiques représentés sur la figure 2. 25 [0021] Comme illustré sur la figure 2, le système informatique E&P (208) comporte un outil E&P (230), un référentiel de données (235) pour le stockage de données intermédiaires et de sorties obtenues de l'outil E&P (230), et un moteur d'exécution de tâche de champ (231) destiné à effectuer diverses tâches de l'opération de champ. Dans un ou plusieurs modes de réalisation, le référentiel de données (235) peut comprendre un dispositif 30 de stockage à lecteur de disque, un dispositif de stockage à semi-conducteur, d'autres dispositifs de stockage de données informatiques appropriés, ou des combinaisons de ceux-ci. Dans un ou plusieurs modes de réalisation, un contenu stocké dans le référentiel de données (235) peut être un fichier de données, une liste chaînée, une séquence de données, une base de 3032532 8 données, une représentation graphique, ou toute autre structure de données appropriée. Le référentiel de données (235) comporte une fonctionnalité permettant de stocker des données telles qu'un volume sismique (227), des données de puits (232), un modèle basé sur le volume (228), et des écarts de modélisation (229). Le référentiel de données (235) peut' comprendre 5 une fonctionnalité permettant de stocker des données supplémentaires non représentées sur la figure 2. Divers types de données stockées dans le référentiel de données sont présentés ci-dessous. [0022] Pendant les opérations de champ illustrées sur la figure 1 ci-dessus, des données peuvent être recueillies à des fins d'analyse et/ou de surveillance des opérations. Ces 10 données peuvent par exemple comprendre des informations concernant des formations souterraines, des équipements, et des données historiques et/ou autres. Les données concernant la formation souterraine peuvent être recueillies au moyen de diverses sources. Ces données de formation peuvent être statiques ou dynamiques. Les données statiques sont des données qui ne varient généralement pas, sauf au cours d'événements géologiques majeurs 15 et de périodes géologiques étendues. Les données statiques concernent par exemple la structure d'une formation et la stratigraphie géologique qui définissent des structures géologiques de la formation souterraine. Les données dynamiques sont des données qui peuvent varier au cours du temps. Les données dynamiques concernent par exemple des fluides s'écoulant à travers les structures géologiques de la formation souterraine au cours du 20 temps. Ces données statiques et/ou dynamiques peuvent être recueillies en vue d'en apprendre davantage sur les formations et sur les gisements précieux qu'elles contiennent. [0023] Les données peuvent être utilisées en vue de prédire les conditions de fond de puits et de prendre des décisions concernant les opérations de champ. Ces décisions peuvent mettre en jeu la planification des puits, le ciblage des puits, les complétions de puits, les 25 niveaux d'exploitation, les débits de production, et d'autres paramètres d'opérations et/ou d'exploitation. Un grand nombre de variables et de grandes quantités de données à prendre en compte peuvent être présents lors de l'analyse d'opérations de champ. En raison du grand nombre de variables et des grandes quantités de données, une modélisation du comportement de l'opération de champ pour déterminer la marche à suivre souhaitée peut être utile. Pendant 30 les opérations en cours, les paramètres d'exploitation peuvent être ajustés au fur et à mesure que les conditions de champ varient et que de nouvelles informations sont reçues. Divers aspects des opérations de champ, tels que des structures géologiques, des gisements de fond 3032532 9 de puits, des puits de forage, des installations de surface, ainsi que d'autres parties de l'opération de champ, peuvent être modélisés. [0024] La prospection sismique peut être effectuée en injectant de l'énergie dans le sol en une ou plusieurs positions de sources, à titre d'exemple, au moyen d'une explosion 5 commandée, d'une entrée mécanique, ou d'une autre technique. L'énergie renvoyée est ensuite mesurée sous la forme de signaux sismiques aux positions de récepteurs de surface à des distances et des azimuts variables par rapport à la position de la source. Le temps de propagation de l'énergie de la source au récepteur, mesuré au moyen des réflexions et les réfractions aux interfaces des strates souterraines, indique la profondeur et l'orientation de ces 10 strates. Les signaux sismiques peuvent être mesurés sous la forme de traces sismiques, qui sont des courbes mathématiques représentant les amplitudes de signaux (c'est-à-dire des amplitudes sismiques) mesurées en fonction de la profondeur ou en fonction du temps de propagation sismique. [0025] Une collection de traces sismiques au sein d'un volume d'intérêt peut être 15 désignée sous le nom de volume sismique (227). En d'autres termes, un volume sismique (227) est la collection de données sismiques (c'est-à-dire de traces sismiques) correspondant à une région souterraine particulière de la Terre, qui est désignée sous le nom de volume d'intérêt. Un volume sismique (227) peut être affiché sous la forme d'images sismiques sur la base de différentes résolutions d'échantillonnage et orientations de visualisation, et en utilisant 20 diverses techniques de traitement d'amplitude sismique différentes en vue d'accentuer ou de mettre en évidence des profils de réflexion sismiques. [0026] Dans un ou plusieurs modes de réalisation, le volume sismique (227) peut être partitionné en cellules. En d'autres termes, chaque cellule représente une sous-région particulière du volume sismique (227). Les sous-régions du volume sismique (227) peuvent 25 avoir des tailles et/ou des formes homogènes ou hétérogènes. A titre d'exemple, les cellules peuvent être définies par un maillage du volume sismique (227), de telle sorte que chaque point d'intersection de lignes de grille au sein de la grille définisse un sommet d'une cellule, et de telle sorte qu'une grille régulière ou qu'une grille adaptive puisse être utilisée. [0027] Dans un ou plusieurs modes de réalisation, les traces sismiques illustrées sur la 30 figure 1 ci-dessus sont fournies au système informatique E&P (208) et stockées dans le référentiel de données (235) en tant que volume sismique (227). Dans un ou plusieurs modes de réalisation, le volume sismique (227) peut être présenté sous la forme d'un volume tridimensionnel (3D) à un utilisateur effectuant une interprétation sismique, lequel est désigné 3032532 10 sous le nom d'utilisateur d'interprétation sismique. La partie supérieure du volume sismique 3D affiché représente la position en surface des traces sismiques individuelles. Au sein du volume sismique 3D, les traces sismiques peuvent être représentées sous la forme de lignes verticales de l'amplitude sismique en fonction du temps ou de la distance le long de l'axe z du 5 volume 3D. Plus précisément, chaque trace individuelle est une représentation de l'amplitude en fonction du temps d'une réflexion acoustique provenant de structures géologiques dans la formation souterraine. D'autres représentations de traces sismiques peuvent être utilisées sans s'écarter du cadre d'un ou de plusieurs modes de réalisation. Une séquence de traces x en fonction du temps est désignée par "ligne" ou "en ligne" en interprétation sismique. Une 10 séquence de traces y en fonction du temps est désignée par "ligne transversale". Une "tranche d'horizon" est une tranche prise dans le volume 3D qui est identifiée par l'utilisateur d'interprétation sismique comme correspondant à un horizon (par exemple l'un des horizons (172, 174, 176)) dans la formation souterraine (104) illustrée sur la figure 1 ci-dessus. [0028] Dans l'ensemble de la présente divulgation, le terme "attributs structurels" 15 désigne des caractéristiques structurelles de la formation souterraine. Dans un ou plusieurs modes de réalisation, les attributs structurels peuvent être estimés en appliquant une combinaison de techniques de traitement du signal et de géométrie différentielle au volume sismique. Les valeurs des attributs structurels directement calculés à partir du volume sismique (227) sont nommées valeurs déduites sismiquement ou attributs sismiques. Plus 20 précisément, les valeurs déduites sismiquement ou les attributs sismiques sont des valeurs des attributs structurels qui sont directement calculées à partir du volume sismique (227) sans aucun autre processus intermédiaire de modélisation basée sur le volume. En d'autres termes, le volume sismique (227) est la seule structure de données basée sur le volume utilisée pour calculer les valeurs déduites sismiquement ou les attributs sismiques. Certains attributs 25 sismiques sont des attributs à une seule trace calculés sur la base d'une trace sismique unique indépendamment les uns des autres, alors que les attributs multi-traces sont calculés sur la base de la relation entre des traces sismiques voisines. De nombreux attributs structurels sont des attributs multi-traces. A titre d'exemple, le pendage estimé des couches géologiques est un attribut multi-trace déterminé à partir de l'anisotropie des traces sismiques ou à partir de la 30 corrélation entre des traces sismiques voisines. Dans un autre exemple, la courbure estimée de couches géologiques est un autre attribut multi-trace qui met en évidence la localisation des surfaces de failles. 3032532 11 [0029] Un type particulier d'attributs structurels concerne les caractéristiques sédimentaires (c'est-à-dire stratigraphiques) de la formation souterraine. Des attributs structurels de ce type sont notamment nommés attributs stratigraphiques. A titre d'exemple, les attributs stratigraphiques comprennent des informations concernant les profils 5 sédimentaires, l'organisation des couches stratigraphiques, ou l'environnement sédimentaire associé au gisement souterrain. La densité de convergence est un exemple de calcul des valeurs déduites sismiquement d'un attribut stratigraphique. [0030] Des images sismiques peuvent également être utilisées en vue d'estimer indirectement les attributs structurels de la formation souterraine par l'intermédiaire d'un 10 processus de modélisation afm de créer un modèle basé sur le volume (228). La variabilité spatiale et/ou temporelle de profils d'empilement, ou de séquences, observée dans des images sismiques se rapporte aux environnements sédimentaires et à des processus post-sédimentaires, comme l'érosion et l'activité tectonique. En d'autres termes, les profils de réflexion contenus dans les images sismiques relient les environnements sédimentaires et 15 l'ordre d'empilement vertical à la séquence de sédimentation dans la formation souterraine. Pendant l'interprétation sismique, la datation relative des profils de réflexion des images sismiques permet de déchiffrer l'historique géologique du sous-sol et conduit à l'estimation de caractéristiques sédimentaires probables. Les résultats d'interprétation et d'analyse sismique peuvent être organisés sous la forme d'un modèle basé sur le volume (228). De manière 20 générale, un modèle basé sur le volume est un volume de données tridimensionnel (3D) qui modélise des caractéristiques géologiques au sein d'un volume d'intérêt. Dans un ou plusieurs modes de réalisation, le volume d'intérêt peut être partitionné en de multiples unités telles que des couches géologiques, des blocs faillés, des corps sédimentaires, etc. La création d'un modèle basé sur le volume peut mettre en jeu l'interpolation ou l'extrapolation de données 25 sismiques interprétées ou l'interpolation ou l'extrapolation de données de puits interprétées. Des informations provenant de données de puits (c'est-à-dire de diagraphies de carottes ou pétrophysiques) et des informations provenant de données sismiques peuvent être combinées afin d'obtenir un modèle basé sur le volume unique. A titre d'exemple, le modèle basé sur le volume (228) peut comprendre des informations qui décrivent des couches stratigraphiques et 30 des blocs faillés dans la formation souterraine, et peut également être nommé modèle structurel. En d'autres termes, le terme de "modèle structurel" désigne un modèle basé sur le volume qui décrit des structures (par exemple des couches stratigraphiques et des blocs faillés) de la formation souterraine. Dans un ou plusieurs modes de réalisation, le modèle basé 3032532 12 sur le volume (228) est un modèle structurel et comporte en outre des informations décrivant l'âge géologique relatif des formations souterraines en chaque point du volume d'intérêt. Des exemples du modèle basé sur le volume (228) sont décrits en référence aux figures 4 et 5 ci-dessous. 5 [0031] Dans un ou plusieurs modes de réalisation, comme pour le volume sismique (227), le modèle basé sur le volume (228) peut être partitionné en cellules. En d'autres termes, chaque cellule représente une sous-région particulière du modèle basé sur le volume (228). Les sous-régions région du modèle basé sur le volume (228) peuvent avoir des tailles et/ou des formes homogènes ou hétérogènes. A titre d'exemple, les cellules peuvent être définies 10 par un maillage du modèle basé sur le volume (228), de telle sorte que chacun des points d'intersection de lignes de grille au sein de la grille définisse un sommet d'une cellule, et de telle sorte qu'une grille régulière ou qu'une grille adaptive puisse être utilisée. Le volume sismique (227) et le modèle basé sur le volume (228) sont des types particuliers de structures de données basées sur le volume contenant respectivement des données sismiques mesurées et 15 des données de description structurelle modélisées. Par ailleurs, dans certains modes de réalisation, une correspondance biunivoque peut exister entre des cellules du modèle basé sur le volume (228) et des cellules contenues dans le volume sismique (227). Dans d'autres modes de réalisation, les cellules contenues dans le modèle basé sur le volume (228) ne sont pas en correspondance avec les cellules contenues dans le volume sismique (227). 20 [0032] Le modèle structurel est utilisé dans un processus de modélisation en vue d'estimer les attributs structurels de la formation souterraine. Les valeurs des attributs structurels estimés en utilisant le modèle structurel sont nommées valeurs déduites structurellement. Comme cela a été noté ci-dessus, le modèle structurel (qui est un type de modèle basé sur le volume) peut être fondé sur des données de puits, des données sismiques, 25 ou une combinaison de données de puits et de données sismiques. Par conséquent, les valeurs déduites structurellement sont des valeurs des attributs structurels qui sont calculées en utilisant le modèle structurel, soit d'une manière qui n'est pas fondée sur des données sismiques, soit de manière indirectement fondée sur le volume sismique (227). De manière générale, le volume sismique (227) est la seule structure de données basée sur le volume 30 utilisée pour calculer des valeurs déduites sismiquement ou des attributs sismiques. En comparaison, les valeurs déduites structurellement ou les attributs structurels sont calculés en utilisant au moins une structure de données basée sur le volume (par exemple le modèle basé sur le volume (228)) distincte ou en plus du volume sismique (227). 3032532 13 [0033] Dans un ou plusieurs modes de réalisation, en plus du volume sismique (227), les données de carottage ou d'autres données de diagraphie pétrophysiques (par exemple le tracé statique (108-2)) illustrées sur la figure 1 ci-dessus sont également fournies au système informatique E&P (208) et stockées dans le référentiel de données (235) en tant que données 5 de puits (232). Dans un ou plusieurs modes de réalisation, les parties supérieures de formation dans la formation souterraine peuvent être estimées sur la base des modifications de lithologie analysées à partir de la carotte ou des données de diagraphie pétrophysiques. [0034] Dans un ou plusieurs modes de réalisation, l'outil E&P (230) est configuré pour faciliter une interprétation sismique et une modélisation de la formation. Comme illustré 10 sur la figure 2, l'outil E&P (230) comporte l'analyseur de volume sismique (222), l'analyseur de données de puits (225), le générateur de modèle basé sur le volume (223), le moteur de modélisation (224), et le moteur d'affichage graphique (226). Chacun de ces composants est décrit ci-dessous. [0035] Dans un ou plusieurs modes de réalisation, l'analyseur de volume sismique 15 (222) est configuré en vue d'analyser le volume sismique (227) afin de générer une valeur déduite sismiquement d'un attribut structurel en une ou plusieurs positions dans la formation souterraine. A titre d'exemple, l'analyseur de volume sismique (222) peut générer la valeur déduite sismiquement sur la base des exemples décrits dans les diverses références dont la liste est fournie ci-dessus. Des exemples de la valeur déduite sismiquement d'un attribut 20 structurel sont décrits en référence à la figure 5 ci-dessous. [0036] Dans un ou plusieurs modes de réalisation, le générateur de modèle basé sur le volume (223) est configuré en vue d'analyser le volume sismique (227) afin de générer le modèle basé sur le volume (228) (par exemple un modèle structurel) de la formation souterraine. Dans un ou plusieurs modes de réalisation, l'analyseur de données de puits (225) 25 est configuré en vue d'estimer les parties de formation supérieures dans la formation souterraine sur la base des modifications de lithologie analysées à partir des données de carottage ou de diagraphie pétrophysiques contenues dans les données de puits (232). Dans un ou plusieurs modes de réalisation, le générateur de modèle basé sur le volume (223) est configuré pour générer le modèle basé sur le volume (228) (par exemple un modèle structurel) 30 de la formation souterraine sur la base des résultats de l'analyseur de données de puits (225) sans recourir à aucune donnée sismique. Dans un ou plusieurs modes de réalisation, le générateur de modèle basé sur le volume (223) est configuré pour générer le modèle basé sur le volume (228) (par exemple un modèle structurel) de la formation souterraine sur la base 3032532 14 d'une combinaison du volume sismique (227) et des résultats de l'analyseur de données de puits (225). [0037] Dans un ou plusieurs modes de réalisation, le moteur de modélisation (224) est configuré pour calculer, sur la base du modèle basé sur le volume (228), une valeur déduite 5 structurellement de l'attribut structurel en une ou plusieurs positions dans la formation souterraine. De plus, le moteur de modélisation (224) est en outre configuré pour comparer la valeur déduite sismiquement et la valeur déduite structurellement afin de générer une valeur de différence représentant un écart (parmi les écarts de modélisation (229)) de modélisation de l'attribut structurel à une position correspondante dans la formation souterraine. Dans un ou 10 plusieurs modes de réalisation, les écarts de modélisation (229) sont minimisés ou réduits d'une autre manière en ce que l'analyseur de volume sismique (222) ajuste l'analyse des données sismiques et/ou que le moteur de modélisation (224) ajuste le modèle basé sur le volume (228). Dans un ou plusieurs modes de réalisation, les écarts de modélisation (229) sont présentés à un utilisateur interpréteur sismique ou à un utilisateur modéliseur structurel 15 par l'intermédiaire du moteur d'affichage graphique (226). Par conséquent, l'utilisateur interpréteur sismique ou l'utilisateur modéliseur structurel peut ajuster manuellement l'analyse des données sismiques et/ou le modèle basé sur le volume (228). Dans un ou plusieurs modes de réalisation, l'outil E&P (230) minimise ou réduit d'une autre manière les écarts de modélisation (229) en utilisant le procédé décrit par référence à la figure 3 ci-dessous. Un 20 exemple de façon de minimiser ou de réduire d'une autre manière les écarts de modélisation (229) sur la base des différences de comparaison est décrit en référence aux figures 4 et 5 ci-dessous. Dans un ou plusieurs modes de réalisation, le moteur d'affichage graphique (226) est configuré pour générer un affichage graphique représentant les résultats de l'outil E&P (230). A titre d'exemple, l'affichage graphique peut être présenté à un utilisateur sur un dispositif 25 d'affichage physique 2D ou 3D et correspond à diverses copies d'écrans illustrées sur les figures 4 et 5 ci-dessous. [0038] Dans un ou plusieurs modes de réalisation, le moteur d'exécution de tâche de champ (231) est configuré pour générer un signal de commande d'opération de champ sur la base au moins de l'horizon interprété (229). Comme cela a été noté ci-dessus, l'équipement 30 d'opération de champ illustré sur la figure 1 ci-dessus peut être commandé par le signal de commande d'opération de champ. Par exemple, le signal de commande d'opération de champ peut être utilisé pour commander des équipements de forage, un actionneur, une vanne de 3032532 15 fluide, ou d'autres dispositifs électriques et/ou mécaniques disposés autour du champ (100) illustré sur la figure 1 ci-dessus. [0039] Le système informatique E&P (208) peut comprendre un ou plusieurs ordinateurs de système, qui peuvent être mis en oeuvre sous la forme d'un serveur ou d'un 5 système informatique classique quelconque. Cependant, l'homme du métier, se fondant sur la présente divulgation, notera que les mises en oeuvre de diverses technologies décrites ici peuvent être effectuées selon d'autres configurations de systèmes informatiques, parmi lesquelles des serveurs à protocole de transfert hypertexte (HTTP), des dispositifs portatifs, des systèmes multiprocesseurs, des dispositifs électroniques grand public à base de 10 microprocesseurs ou programmables, des ordinateurs personnels en réseau, des miniordinateurs, des ordinateurs centraux, etc. [0040] Bien que des composants particuliers soient illustrés et/ou décrits pour une utilisation dans des unités et/ou des modules du système informatique E&P (208) et de l'outil E&P (230), une grande variété de composants présentant diverses fonctions peuvent être 15 utilisés en vue d'assurer les fonctions de formatage, de traitement, d'utilitaires et de coordination pour le système informatique E&P (208) et l'outil E&P (230). Les composants peuvent avoir des fonctionnalités combinées et peuvent être mis en oeuvre sous la forme de logiciels, de matériels, de microcodes, ou de combinaisons de ceux-ci. [0041] La figure 3 représente un exemple de procédé destiné à des attributs sismiques 20 déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume conformément à un ou plusieurs modes de réalisation. A titre d'exemple, le procédé représenté sur la figure 3 peut être mis en oeuvre en utilisant le système informatique E&P (208) et l'outil E&P (230) décrits en référence à la figure 2 ci-dessus pour le champ (100) décrit en référence à la figure 1 ci-dessus. Dans un ou plusieurs modes de réalisation, un ou plusieurs des éléments 25 représentés sur la figure 3 peuvent être omis, répétés, et/ou mis en oeuvre dans un ordre différent. Par conséquent, les modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume ne doivent pas être considérés comme étant limités aux agencements particuliers d'éléments représentés sur la figure 3. [0042] En premier lieu, dans l'Elément 301, une valeur déduite sismiquement d'un 30 attribut structurel est calculée sur la base d'un volume sismique. Dans un ou plusieurs modes de réalisation, chaque cellule du volume sismique est associée à une valeur déduite sismiquement d'un ou plusieurs attributs structurels. A titre d'exemple, le pendage et l'azimut de pendage, la courbure minimum, maximum et Gaussienne, ou la variation latérale de 3032532 16 l'épaisseur le long d'une couche sont des exemples d'attributs structurels pour lesquels des valeurs déduites sismiquement peuvent être déduites directement à partir de l'analyse de traces sismiques voisines dans le volume sismique. [0043] Dans l'Elément 302, une interprétation sismique basée sur le volume sismique 5 est tout d'abord effectuée ou ajustée de manière itérative. En d'autres termes, des informations préalablement connues concernant la vitesse de signaux à travers divers types de roches peuvent être utilisées pour vérifier, à partir des traces sismiques dans le volume sismique, des informations concernant des formations sous-terraines estimées, comme le type de roche, la porosité, la perméabilité, etc. Des informations supplémentaires collectées à partir du champ, 10 comme des carottages, peuvent être utilisées en vue d'effectuer l'interprétation sismique. Dans l'Elément 303, un modèle structurel de la formation souterraine est généré ou ajusté de manière itérative à partir de résultats de l'interprétation sismique. En particulier, à partir de l'interprétation sismique correspondant à chacune des traces sismiques, un volume structurel est généré. 15 [0044] Dans l'Elément 304, une valeur déduite structurellement de l'attribut structurel est calculée sur la base du modèle structurel. Dans un ou plusieurs modes de réalisation, sur la base d'une séquence de couches structurelles décrites dans le modèle structurel, les âges géologiques relatifs des couches structurelles sont déterminés. Les âges géologiques relatifs des couches structurelles sont quant à eux utilisés pour déterminer la valeur déduite 20 structurellement de l'attribut structurel. A titre d'exemple, un attribut de pendage de couche, un attribut de densité de convergence, un attribut de courbure de couche, et/ou un attribut d'épaisseur de couche peuvent être déterminés sur la base des âges géologiques relatifs des couches structurelles. [0045] Dans l'Elément 305, la valeur déduite sismiquement et la valeur déduite 25 structurellement sont comparées pour générer une valeur de différence pour chaque cellule du volume sismique et/ou du modèle structurel. En particulier, chaque valeur de différence représente un écart de modélisation de l'attribut structurel à une position correspondante dans la formation souterraine. Dans un ou plusieurs modes de réalisation, des cellules sismiques (c'est-à-dire des cellules contenues dans le volume sismique) et les cellules structurelles (c'est- 30 à-dire les cellules contenues dans le modèle structurel) sont comparées pour générer des cellules rééchantillonnées. Plus précisément, les cellules rééchantillonnées sont définies par une grille rééchantillonnée du volume d'intérêt associé au volume sismique et au modèle structurel. Par conséquent, les valeurs de différence sont générées pour chacune des cellules 3032532 17 rééchantillonnées par interpolation et extrapolation à partir de grilles initiales associées au volume sismique et au modèle structurel pour obtenir la grille rééchantillonnée. [0046] Dans un ou plusieurs modes de réalisation, diverses cellules peuvent ne pas être soumises à une analyse plus poussée. A titre d'exemple, des zones de voisinage de failles 5 sont identifiées à moins d'une distance prédéterminée des failles décrites dans le modèle structurel. Par conséquent, les valeurs de différence des cellules qui se trouvent au sein d'une zone de voisinage de faille prédéfinie ne sont pas soumises à une analyse plus poussée. En d'autres termes, ces cellules ne sont pas utilisées pour déterminer si des mises à jour doivent être effectuées et pour effectuer d'éventuelles mises à jour. A titre d'autre exemple, les valeurs 10 de différence d'autres cellules peuvent également être rejetées sur la base de la valeur d'un attribut sismique concernant le niveau de qualité de l'image sismique (par exemple le rapport signal à bruit, le coefficient de corrélation entre des traces sismiques voisines, etc.) ou le niveau de certitude de l'attribut sismique considéré. Le niveau de qualité de l'image sismique ou le niveau de certitude de l'attribut sismique considéré est nommé mesure de qualité 15 sismique. Dans un ou plusieurs modes de réalisation, après le rejet des valeurs de différence mentionnées ci-dessus, les valeurs de différence qualifiées restantes sont affichées dans une fenêtre graphique sur un écran d'affichage physique. [0047] Dans un ou plusieurs modes de réalisation, l'attribut structurel comporte un attribut de pendage de couche. Dans ces modes de réalisation, les valeurs de différence 20 peuvent être affichées pendant l'interprétation sismique, cela permettant à un utilisateur interpréteur sismique d'identifier une faille supplémentaire du volume sismique n'ayant pas été incluse dans le modèle structurel. Par conséquent, un modèle structurel révisé peut être généré en ajoutant la faille supplémentaire aux failles déjà présentes dans le modèle structurel. [0048] Dans l'Elément 306, il est déterminé si la valeur de différence dépasse un seuil 25 prédéterminé. Dans un ou plusieurs modes de réalisation, un paramètre statistique (par exemple le maximum, la médiane, la moyenne, etc.) des valeurs de différence correspondant aux cellules du volume sismique et/ou du modèle structurel sont utilisées pour cette comparaison. Si la détermination est positive, c'est-à-dire si la valeur de différence (par exemple le maximum, la médiane, la moyenne, etc. des valeurs de différence) dépasse le seuil 30 prédéterminé, le procédé passe à l'Elément 308. Si la détermination est négative, c'est-à-dire si la valeur de différence (par exemple le maximum, la médiane, la moyenne, etc. des valeurs de différence) ne dépasse pas le seuil prédéterminé, le procédé passe à l'Elément 307, dans lequel une opération de champ est effectuée sur la base du modèle structurel. 3032532 18 [0049] Dans l'Elément 308, il est déterminé si les valeurs déduites sismiquement doivent être ajustées pour réduire les valeurs de différence. Si la détermination est positive, c'est-à-dire si les valeurs déduites sismiquement doivent être ajustées, le procédé retourne à l'Elément 301, dans lequel les valeurs déduites sismiquement sont ajustées de manière itérative. A titre d'exemple, les paramètres utilisés pour générer les valeurs déduites sismiquement peuvent être ajustés automatiquement ou sur la base d'une entrée de l'utilisateur. [0050] Si la détermination effectuée dans l'Elément 308 est négative, c'est-à-dire si les valeurs déduites sismiquement ne doivent pas être ajustées, le procédé passe à l'Elément 309, 10 dans lequel une autre détermination est effectuée pour savoir si l'interprétation sismique doit être ajustée pour réduire les valeurs de différence. Si la détermination est positive, c'est-à-dire si l'interprétation sismique doit être ajustée, le procédé retourne à l'Elément 302, dans lequel l'interprétation sismique est ajustée de manière itérative. A titre d'exemple, les valeurs de différence peuvent être présentées à un utilisateur interpréteur sismique à titre de référence, de 15 manière à ce que l'interprétation sismique puisse être ajustée sur la base d'une entrée de l'utilisateur interpréteur sismique. Si la détermination est négative, c'est-à-dire si l'interprétation sismique ne doit pas être ajustée, le procédé retourne à l'Elément 303, dans lequel le modèle structurel est ajusté de manière itérative. A titre d'exemple, une faille supplémentaire identifiée sur la base des valeurs de différence peut être automatiquement 20 ajoutée au modèle structurel pour générer un modèle structurel révisé. [0051] Les figures 4 et 5 représentent des exemples de copies d'écran pour des attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume conformément à un ou plusieurs modes de réalisation. Plus précisément, la copie d'écran A (311), la copie d'écran B (312), la copie d'écran C (313), la copie d'écran D (314), la 25 copie d'écran E (301), et la copie d'écran F (302) sont générées en utilisant le PROCESSUS 1 présenté ci-dessous en vue d'identifier des failles peu marquées sur la base de l'attribut structurel "pendage". En particulier, la procédure de processus déroulée dans le PROCESSUS 1 peut être fondée sur l'outil E&P décrit en référence à la figure 2 et le procédé décrit en référence à la figure 3 ci-dessus.

30 PROCESSUS 1 (A) Calculer le même attribut structurel sur la base (1) des traces sismiques dans un volume sismique et (2) d'un attribut d'âge géologique relatif dans un modèle structurel 3D déduit à partir de ce volume sismique. En particulier, une valeur 3032532 19 déduite sismiquement de l'attribut structurel est calculée pour chaque cellule dans le volume sismique, et une valeur déduite structurellement de l'attribut structurel est calculée pour chaque cellule dans le modèle structurel 3D. Plus précisément, les valeurs déduites sismiquement sont calculées sans utiliser de modèle basé sur le 5 volume, comme le modèle structurel 3D et l'attribut d'âge géologique relatif associé. (B) Rééchantilloner les valeurs déduites structurellement calculées en utilisant l'attribut d'âge géologique relatif pour obtenir les cellules du volume sismique OU augmenter l'échelle des valeurs déduites sismiquement pour obtenir les cellules du modèle structurel 3D. 10 (C) Facultativement, normaliser ou modifier l'échelle des valeurs déduites sismiquement et des valeurs déduites structurellement afin qu'elles se situent dans la même gamme, ou qu'elles aient des distributions semblables. (D) Calculer soit la différence soit le rapport entre les valeurs déduites sismiquement et les valeurs déduites structurellement au niveau de chaque cellule du 15 modèle. (E) Facultativement, soumettre à un post-traitement la différence en appliquant un seuil ou un algorithme de rehaussement de contours. (F) Facultativement, affecter à l'attribut calculé un attribut de "fiabilité" qualifiant le niveau de robustesse de l'attribut calculé en (A). 20 (G) On s'attend à ce que la différence mette en évidence des écarts entre les données sismiques et le modèle structurel 3D extrait des données sismiques. Dans un ou plusieurs modes de réalisation, des techniques de courbes de niveau ou de rendu de volumes peuvent être utilisées pour visualiser les résultats et mettre en relation les zones de fort écart avec des caractéristiques particulières de l'interprétation sismique 25 ou du modèle basé sur le volume. [0052] En plus d'identifier des failles peu marquées sur la base de l'attribut de "pendage", la procédure de processus déroulée dans le PROCESSUS 1 peut être appliquée à d'autres scénarios, comme l'utilisation des attributs de courbure pour mettre en évidence des plis ou des corps géologiques peu marqués dans le modèle structurel, en utilisant les 30 modifications d'épaisseur des couches structurelles (c'est-à-dire la divergence du vecteur de pendage) afin d'identifier des profils sédimentaires, etc. 3032532 20 [0053] Comme illustré sur la figure 4, la copie d'écran A (311) illustre une coupe transversale d'un volume sismique contenant des données sismiques (311-1). Les signaux sismiques sont généralement très bruyants au voisinage de surfaces de failles, qui correspondent à une discontinuité à fort pendage dans la formation souterraine. Cette 5 discontinuité à fort pendage crée de forts contrastes dans la copie d'écran A (311). Il en résulte que les attributs "pendage" et "azimut de pendage" calculés à partir de ces zones dans le volume sismique peuvent être bruités et/ou peu fiables. La copie d'écran B (312) et la copie d'écran C (313) représentent respectivement des valeurs déduites sismiquement de l'attribut "pendage" (c'est-à-dire des valeurs de pendage déduites sismiquement (312-1)) et de l'attribut 10 "azimut de pendage" (c'est-à-dire des valeurs d'azimut de pendage déduites sismiquement (313-1)), sur la base du volume sismique illustré sur la copie d'écran A (311). [0054] La copie d'écran D (314) représente des valeurs de pendage déduites structurellement (314-1) qui sont générées sur la base d'un modèle structurel déduit du volume sismique illustré sur la copie d'écran A (311). En particulier, la copie d'écran D (314) 15 représente des failles (par exemple la faille (314-1)) superposées aux valeurs de pendage déduites structurellement (314-1). Contrairement aux valeurs déduites sismiquement représentées sur la copie d'écran B (312) et la copie d'écran C (313), les valeurs de pendage déduites structurellement (314-1) sont discontinues au franchissement des failles (par exemple la faille (314-1)), mais régulières en tout point éloigné des failles. Par conséquent, les plus 20 grandes différences de comparaison entre les valeurs de pendage déduites sismiquement (3121) et les valeurs de pendage déduites structurellement (314-1) correspondent à des failles et à des voisinages de failles de médiocre qualité sismique. A titre d'exemple, la comparaison peut être effectuée en utilisant la procédure de processus déroulée dans le PROCESSUS 2 ci-dessous.

25 PROCESSUS 2 (A) Calculer le vecteur normal à la couche structurelle sur la base des valeurs de pendage déduites sismiquement (312-1) et des valeurs d'azimut de pendage déduites sismiquement (313-1). (B) Calculer le gradient de l'âge géologique relatif dans le modèle structurel. 30 (C) Rééchantilloner ce gradient pour obtenir le volume sismique. (D) Calculer l'angle entre le vecteur normal et le gradient rééchantillonné. L'angle calculé est l'attribut d'écart de pendage "brut". 3032532 21 [0055] Comme illustré sur la figure 5, la copie d'écran E (301) et la copie d'écran F (302) illustrent divers aspects du modèle structurel déduit du volume sismique illustré sur la copie d'écran A (311). Plus précisément, la copie d'écran E (301) représente un plan de faille (301-7) coupant des couches structurelles (301-1) à (301-6). La copie d'écran F (302) 5 représente des valeurs d'une fonction implicite représentant les âges géologiques relatifs de couches structurelles (302-1) à (302-6), correspondant aux couches structurelles (301-1) à (301-6) illustrées sur la copie d'écran E (301). Par ailleurs, le plan de faille (302-7) illustré sur la copie d'écran F (302) correspond au plan de faille (301-7) illustré sur la copie d'écran E (301). A titre d'exemple, la fonction implicite représentée sur la copie d'écran F (302) peut 10 être déduite des données sismiques (311-1). Par conséquent, la fonction implicite peut être utilisée pour calculer le gradient de l'âge géologique relatif. [0056] Cet attribut d'écart "brut" peut être soumis à un post-traitement au moyen de la procédure de processus déroulée dans le PROCESSUS 3 ci-dessous en vue d'identifier les failles qui sont présentes dans les données sismiques mais qui ont été prises en compte dans le 15 modèle structurel. PROCESSUS 3 (A) Appliquer un seuil à l'attribut d'écart de pendage "brut" pour ne retenir que les valeurs d'écart élevées (par exemple rejeter tous les angles inférieurs à 5 degrés). (B) Filtrer les zones situées dans le voisinage immédiat de failles identifiées 20 dans le modèle structurel : (B.1) Calculer un attribut de "voisinage de faille" sur la base de la distance à la faille la plus proche, par exemple sous la forme d'une Carte de Distance Euclidienne calculée en utilisant l'algorithme de SaitoToriwaki, au voisinage de failles du cadre structurel. 25 (B.2) Définir une région dans le modèle structurel qui regroupe les points proches des surfaces de failles de plus d'un seuil de distance prédéterminé. (B.3) Utiliser cette région de "voisinage de faille" pour annuler des parties sélectionnées de l'attribut d'écart. 30 (C) Appliquer une technique de filtrage du bruit et de rehaussement de contours : 3032532 22 (C.1) Facultativement éliminer des voxels isolés "à fort écart" en utilisant des opérateurs de morphologie mathématiques. (C.2) Appliquer des techniques de rehaussement de contours telles que la technique "ant-tracking" en vue d'identifier des linéaments reliant les zones d'attributs de fort écart. [0057] La sortie du processus déroulé dans le PROCESSUS 3 est un attribut sismique mettant en évidence les failles peu marquées qui sont identifiées dans le volume sismique mais qui n'ont pas été repérées dans le modèle structurel. Plus précisément, les failles peu marquées n'ont pas été repérées dans le modèle structurel lorsque le modèle structurel a été 10 déduit du volume sismique. [0058] Les modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peuvent être mis en oeuvre sur un système informatique. Une combinaison quelconque de matériels mobiles, de bureau, serveurs, intégrés ou d'autres types de matériels peut être utilisée. A titre d'exemple, le système 15 informatique peut être constitué par un ou plusieurs dispositifs mobiles (par exemple un ordinateur portable, un téléphone intelligent, un assistant numérique personnel, un ordinateur tablette, ou un autre dispositif mobile), des ordinateurs de bureau, des serveurs, des cartes dans un châssis de serveur, ou tout autre type de dispositif ou de dispositifs informatique(s) qui comporte/comportent au moins la puissance de traitement, la mémoire, et le(s) 20 dispositif(s) d'entrée et de sortie minimum pour mettre en oeuvre un ou plusieurs modes de réalisation. A titre d'exemple, comme illustré sur la figure 6, le système informatique (600) peut comprendre un ou plusieurs processeur(s) d'ordinateur (602), la mémoire associée (604) (par exemple une mémoire vive (RAM), une mémoire cache, une mémoire flash, etc.), un ou plusieurs dispositif(s) de stockage (606) (par exemple un disque dur, un lecteur optique tel 25 qu'un lecteur de disque compact (CD) ou un lecteur de disque numérique polyvalent (DVD), une clé de mémoire flash, etc.), et de nombreux autres éléments et fonctionnalités. Le/les processeur(s) d'ordinateur (602) peut/peuvent être un circuit intégré destiné à traiter des instructions. A titre d'exemple, le/les processeur(s) d'ordinateur peut/peuvent être constitué(s) par un ou plusieurs coeurs, ou micro-coeurs d'un processeur. 30 [0059] Le système informatique (600) peut également comprendre un ou plusieurs dispositif(s) d'entrée (610), comme un écran tactile, un clavier, une souris, un microphone, un pavé tactile, un stylet électronique, ou tout autre type de dispositif d'entrée. Par ailleurs, le système informatique (600) peut comprendre un ou plusieurs dispositif(s) de sortie (608), 3032532 23 comme un écran (par exemple un affichage à cristaux liquides (LCD), un affichage à plasma, un écran tactile, un écran à tube cathodique (CRT), un projecteur, ou un autre dispositif d'affichage), une imprimante, une unité de stockage externe, ou tout autre dispositif de sortie. Un ou plusieurs dispositif(s) de sortie peut/peuvent être identique(s) au dispositif d'entrée ou 5 différent(s) de celui-ci. Le système informatique (600) peut être connecté à un réseau (612) (par exemple un réseau local (LAN), un réseau étendu (WAN) tel que l'Internet, un réseau mobile, ou tout autre type de réseau) via une connexion par une interface réseau (non représentée). Le(s) dispositif(s) d'entrée et de sortie peut/peuvent être connecté(s) localement ou à distance (par exemple via le réseau (612)) au/aux processeur(s) d'ordinateur (602), à la 10 mémoire (604), et au/aux dispositif(s) de stockage (606). Il existe de nombreux types différents de systèmes informatiques, et le(s) dispositif(s) d'entrée et de sortie mentionné(s) ci-dessus peut/peuvent se présenter sous d'autres formes. [0060] Des instructions logicielles sous la forme de code de programme lisible par ordinateur pour mettre en oeuvre des modes de réalisation peuvent être stockées, en totalité ou 15 en partie, de manière temporaire ou permanente, sur un support non transitoire lisible par ordinateur tel qu'un CD, un DVD, un dispositif de stockage, une disquette, une bande, une mémoire flash, une mémoire physique, ou tout autre support de stockage lisible par ordinateur. Plus précisément, les instructions logicielles peuvent correspondre à un code de programme lisible par ordinateur qui, lorsqu'il est exécuté par un/des processeur(s), est 20 configuré pour mettre en oeuvre des modes de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume. [0061] Par ailleurs, un ou plusieurs éléments du système informatique mentionné ci- dessus (600) peut/peuvent être situé(s) en un lieu distant et être connecté(s) aux autres éléments par l'intermédiaire d'un réseau (612). Par ailleurs, les modes de réalisation peuvent 25 être mis en oeuvre sur un système distribué comportant une pluralité de noeuds, dans lequel chaque partie des attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume peut être située sur un noeud différent au sein du système distribué. Dans un mode de réalisation d'attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume, le noeud correspond à un dispositif informatique 30 distinct. Le noeud peut correspondre à un processeur d'ordinateur auquel est associée une mémoire physique. Le noeud peut correspondre à un processeur d'ordinateur ou à un micro-coeur d'un processeur d'ordinateur comportant une mémoire et/ou des ressources partagées. 3032532 24 [0062] Les systèmes et procédés présentés se rapportent à l'acquisition d'hydrocarbures à partir d'un champ pétrolifère. Il est à noter que des systèmes et procédés identiques peuvent être utilisés pour effectuer des opérations souterraines, telles que l'exploitation minière, l'extraction d'eau, et l'acquisition d'autres fluides souterrains ou d'autres géomateriaux à partir d'autres champs. Par ailleurs, des parties des systèmes et procédés peuvent être mises en oeuvre sous la forme de logiciels, de matériels, de microcodes, ou de combinaisons de ceux-ci. [0063] Bien que l'objet des attributs sismiques déduits de la propriété d'âge géologique relatif d'un modèle basé sur le volume ait été décrit en rapport avec un nombre limité de 10 modes de réalisation, l'homme du métier, se fondant sur la présente divulgation, notera qu'il est possible de concevoir d'autres modes de réalisation qui ne s'écartent pas du cadre divulgué ici. Par conséquent, l'étendue de l'invention doit être considérée comme n'étant limitée que par les revendications annexées.

Claims (15)

1. REVENDICATIONS1. Procédé pour modéliser une formation souterraine (104) d'un champ (100), comprenant: obtenir un volume sismique (227) comprenant une pluralité de traces sismiques de la formation souterraine (104) du champ (100) ; calculer, sur la base du volume sismique (227), une valeur déduite sismiquement d'un attribut structurel représentant une caractéristique structurelle de la formation souterraine (104) ; calculer, sur la base d'un modèle structurel, une valeur déduite structurellement de l'attribut structurel, le modèle structurel comprenant une pluralité de couches structurelles de la formation souterraine (104) ; comparer la valeur déduite sismiquement et la valeur déduite structurellement pour générer une valeur de différence représentant un écart (229) de modélisation de l'attribut structurel à une position correspondante dans la formation souterraine (104) ; et générer un résultat d'interprétation sismique sur la base de la valeur de différence et de la position correspondante.
2. 2. Procédé selon la revendication 1, comprenant en outre: afficher, dans une fenêtre graphique sur un écran d'affichage physique, la valeur de différence et la position correspondante ; et effectuer, en réponse à l'affichage, une interprétation sismique du volume sismique (227) pour générer un résultat d'interprétation sismique, dans lequel l'interprétation sismique est ajustée sur la base au moins de la valeur de différence et la position correspondante.
3. 3. Procédé selon la revendication 2, comprenant en outre: recevoir, en réponse à la génération du résultat d'interprétation sismique, un ajustement du modèle structurel pour générer le modèle structurel révisé, dans lequel l'ajustement est fondé sur la valeur de différence et la position correspondante ; et effectuer une opération de champ sur la base du modèle structurel révisé.
4. 4. Procédé selon la revendication 3, comprenant en outre: 3032532 26 calculer, sur la base du modèle structurel révisé, une valeur déduite structurellement ajustée de l'attribut structurel ; comparer la valeur déduite sismiquement et la valeur déduite structurellement ajustée pour générer une valeur de différence ajustée représentant un écart (229) de 5 modélisation de l'attribut structurel sur la base du modèle structurel révisé ; et afficher, dans la fenêtre graphique sur l'écran d'affichage physique, la valeur de différence ajustée, dans lequel l'interprétation sismique fait partie d'une séquence d'interprétations sismiques qui sont réalisées de manière itérative jusqu'à ce que la valeur de 10 différence ajustée respecte un critère prédéterminé.
5. 5. Procédé selon la revendication 3, dans lequel l'attribut structurel comprend un attribut de pendage de couche, dans lequel le modèle structurel comprend en outre une pluralité de failles, le procédé comprenant en outre: détecter, pendant l'interprétation sismique, une faille supplémentaire du volume 15 sismique (227) sur la base d'au moins la valeur de différence de l'attribut de pendage de couche, dans lequel la faille supplémentaire n'a pas été introduite dans le modèle structurel ; et générer le modèle structurel révisé en ajoutant la faille supplémentaire à la pluralité de failles. 20
6. 6. Procédé selon la revendication 1, comprenant en outre: déterminer, sur la base d'une séquence de la pluralité de couches structurelles, une pluralité d'âges géologiques relatifs de la pluralité de couches structurelles, dans lequel la valeur déduite structurellement est déterminée sur la base de la pluralité d'âges géologiques relatifs. 25
7. 7. Procédé selon la revendication 1, dans lequel l'attribut structurel comprend au moins un élément sélectionné dans un groupe constitué d'un attribut de pendage de couche, d'un attribut de densité de convergence, d'un attribut de courbure de couche, et d'un attribut d'épaisseur de couche.
8. 8. Procédé selon la revendication 1, 30 dans lequel le volume sismique (227) comprend une pluralité de cellules sismiques tridimensionnelles (3D) qui sont associées à une pluralité de valeurs déduites 3032532 27 sismiquement de l'attribut structurel, dans lequel la pluralité de valeurs déduites sismiquement comprennent la valeur déduite sismiquement, dans lequel le modèle structurel comprend une pluralité de cellules structurelles 3D qui sont associées à une pluralité de valeurs déduites structurellement de 5 l'attribut structurel, dans lequel la pluralité de valeurs déduites structurellement comprennent la valeur déduite structurellement, dans lequel la valeur de différence fait partie d'une pluralité de valeurs de différence qui sont générées en comparant la pluralité de valeurs déduites sismiquement et la pluralité de valeurs déduites structurellement, 10 dans lequel la comparaison de la pluralité de valeurs déduites sismiquement et de la pluralité de valeurs déduites structurellement comprend une mise en correspondance de la pluralité de cellules sismiques et de la pluralité de cellules structurelles pour générer une pluralité de cellules rééchantillonnées, dans lequel la pluralité de valeurs de différence sont associées à la pluralité de cellules 15 rééchantillonnées, le procédé comprenant en outre: afficher, dans une fenêtre graphique sur l'écran d'affichage physique, la pluralité de valeurs de différence relativement à la pluralité de cellules rééchantillonnées.
9. 9. Procédé selon la revendication 8, dans lequel le modèle structurel comprend en outre une pluralité de failles, le procédé comprenant en outre: 20 déterminer une pluralité de zones de voisinage de failles à au plus une distance prédéterminée de la pluralité de failles ; et qualifier un sous-ensemble de la pluralité de valeurs de différence sur la base d'au moins un élément sélectionné dans un groupe constitué (i) d'un sous-ensemble correspondant de la pluralité de cellules rééchantillonnées se situant en dehors 25 de la pluralité de zones de voisinage de failles et (ii) d'une mesure de qualité sismique du sous-ensemble correspondant de la pluralité de cellules rééchantillonnées, dans lequel l'affichage est effectué en réponse à la qualification.
10. 10. Système pour modéliser une formation souterraine (104) d'un champ (100), comprenant : 30 un analyseur de volume sismique (227) s'exécutant sur un processeur d'ordinateur et configuré pour : calculer, sur la base d'un volume sismique (227) comprenant un ensemble de traces sismiques d'une formation souterraine (104) du champ (100), une 3032532 28 valeur déduite sismiquement d'un attribut structurel représentant une caractéristique structurelle de la formation souterraine (104) ; un moteur de modélisation s'exécutant sur le processeur d'ordinateur et configuré pour : 5 calculer, sur la base d'un modèle structurel, une valeur déduite structurellement de l'attribut structurel, le modèle structurel comprenant une pluralité de couches structurelles de la formation souterraine (104) ; comparer la valeur déduite sismiquement et la valeur déduite structurellement pour générer une valeur de différence représentant un écart (229) de 10 modélisation de l'attribut structurel à une position correspondante dans la formation souterraine (104) ; et générer un résultat d'interprétation sismique sur la base de la valeur de différence et de la position correspondante ; et un référentiel de données relié au processeur d'ordinateur et configuré pour stocker le 15 volume sismique (227), le modèle structurel, et la valeur de différence avec la position correspondante.
11. 11. Système selon la revendication 10, comprenant en outre : un moteur d'affichage graphique configuré en vue d'afficher, dans une fenêtre graphique sur un écran d'affichage physique, la valeur de différence et la 20 position correspondante, dans lequel l'analyseur de volume sismique (227) est en outre configuré pour effectuer, en réponse à l'affichage, une interprétation sismique du volume sismique (227) pour générer un résultat d'interprétation sismique, dans lequel l'interprétation sismique est ajustée sur la base au moins de la valeur de 25 différence et la position correspondante.
12. 12. Système selon la revendication 11, dans lequel le moteur de modélisation est en outre configuré pour : recevoir, en réponse à la génération du résultat d'interprétation sismique, un ajustement du modèle structurel pour générer le modèle structurel révisé, dans 30 lequel l'ajustement est fondé sur la valeur de différence et la position correspondante, dans lequel le système comprend en outre : 3032532 29 un générateur de modèle structurel configuré en vue d'analyser le volume sismique (227) afin de générer le modèle structurel de la formation souterraine (104) ; et un moteur d'exécution de tâche de champ (100) configuré en vue d'effectuer 5 une opération de champ sur la base du modèle structurel révisé.
13. 13. Système selon la revendication 12, dans lequel le moteur de modélisation est en outre configuré pour : calculer, sur la base du modèle structurel révisé, une valeur déduite structurellement ajustée de l'attribut structurel ; 10 comparer la valeur déduite sismiquement et la valeur déduite structurellement ajustée pour générer une valeur de différence ajustée représentant un écart (229) de modélisation de l'attribut structurel sur la base du modèle structurel révisé ; et afficher, dans la fenêtre graphique sur l'écran d'affichage physique, la valeur de différence ajustée, 15 dans lequel l'interprétation sismique fait partie d'une séquence d'interprétations sismiques qui sont réalisées de manière itérative jusqu'à ce que la valeur de différence ajustée respecte un critère prédéterminé.
14. 14. Système selon la revendication 12, dans lequel l'attribut structurel comprend un attribut de pendage de couche, 20 dans lequel le modèle structurel comprend en outre une pluralité de failles, dans lequel l'analyseur de volume sismique (227) est en outre configuré pour : détecter, pendant l'interprétation sismique, une faille supplémentaire du volume sismique (227) sur la base d'au moins la valeur de différence de l'attribut de pendage de couche, 25 dans lequel la faille supplémentaire n'a pas été introduite dans le modèle structurel, et dans lequel le générateur de modèle structurel est en outre configuré pour générer le modèle structurel révisé en ajoutant la faille supplémentaire à la pluralité de failles. 30
15. 15. Programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes d'un procédé selon l'une quelconque des revendications 1-9 lorsque ledit programme est exécuté sur un ordinateur.