FR3043213A1

FR3043213A1 -

Info

Publication number: FR3043213A1
Application number: FR1560583A
Authority: FR
Inventors: Jimmy Klinger
Original assignee: Services Petroliers Schlumberger SA
Current assignee: Services Petroliers Schlumberger SA
Priority date: 2015-11-04
Filing date: 2015-11-04
Publication date: 2017-05-05
Also published as: EP3371627B1; WO2017076951A1; US11385369B2; EP3371627A1; US20180313966A1

Abstract

L'invention concerne des procédés, systèmes et supports lisibles par ordinateur pour le traitement de données sismiques. Le procédé comporte la réception (104) d'une pluralité de traces sismiques représentant un domaine souterrain, et la réception (102) d'un modèle stratigraphique implicite d'au moins une partie du domaine souterrain. Le procédé comporte également la sélection (108) d'une isovaleur dans le modèle stratigraphique implicite, et la définition (110), au moyen d'un processeur, d'un intervalle géologiquement cohérent dans le modèle stratigraphique implicite sur la base au moins en partie d'une position de l'isovaleur dans le modèle stratigraphique implicite. Le procédé comporte en outre le calcul (124) d'un ou plusieurs attributs de la pluralité de traces sismiques dans l'intervalle.

Description

ATTRIBUTS SISMIQUES TRIDIMENSIONNELS STRATIGRAPfflQUEMENT

COHÉRENTS

Contexte de l'invention [0001] L'interprétation sismique est un processus par lequel des traces sismiques enregistrées sont utilisées pour déterminer les caractéristiques d'un volume (ou d'une zone) souterrain. Des traces sismiques peuvent être acquises à l'aide de géophones positionnés en des emplacements différents autour d'une source sismique, qui peut être une explosion. L'explosion peut générer des ondes sismiques qui se propagent à travers le sol. Certaines des ondes peuvent se propager vers le bas, jusqu'à ce qu'elles soient réfléchies par un réflecteur tel qu'une interface entre deux types de roches. Les ondes réfléchies peuvent ensuite se propager en sens inverse vers la surface, par exemple à une certaine distance horizontale de la source. Les géophones peuvent mesurer l'arrivée des ondes sismiques, et permettent de produire une ou plusieurs traces sismiques à partir de celles-ci.

[0002] Les traces sismiques peuvent être agencées sous la forme d'un cube sismique (par exemple en trois dimensions), qui peut représenter le volume souterrain, y compris des faciès, des strates, etc. dans la roche à travers laquelle les ondes sismiques se propagent. Une fonction implicite peut être développée pour le domaine souterrain, sur la base des données sismiques et de toutes autres données disponibles (par exemple des diagraphies de puits, des carottes, etc.). Diverses techniques de modélisation, par exemple des modèles d'écoulement de fluide, des modèles de bassins, etc., peuvent ensuite être construites à partir des données sismiques, ainsi que de toutes autres données disponibles pour le volume (par exemple des diagraphies de puits, des carottes, etc.).

[0003] Dernièrement, des progrès ont été réalisés dans les flux de travaux d'interprétation sismique, en ayant pour cible la génération d'un modèle interprété automatiquement du domaine souterrain. Cela permet de faciliter l'isolation d'une surface concordante correspondant à un certain âge géologique en sélectionnant une isovaleur à l'intérieur d'un volume de fonction stratigraphique. Ces progrès peuvent de manière générale viser à générer un cube de Temps Géologique Relatif (RGT) à l'aide des informations contenues dans le cube sismique. A cet effet, le volume peut avoir un niveau de détail élevé, de manière à ce qu'un interprète puisse extraire des tranches d'isovaleurs pour l'analyse topographique. Résumé de l'invention [0004] Des modes de réalisation de l'invention permettent de fournir des procédés, des systèmes, et des supports lisibles par ordinateur pour le traitement de données sismiques. A titre d'exemple, le procédé peut comprendre la réception d'une pluralité de traces sismiques représentant un domaine souterrain, et la réception d'un modèle stratigraphique implicite d'au moins une partie du domaine souterrain. Le procédé peut également comprendre la sélection d'une isovaleur dans le modèle stratigraphique implicite, et la définition, au moyen d'un processeur, d'un intervalle géologiquement cohérent dans le modèle stratigraphique implicite sur la base au moins en partie d'une position de l'isovaleur dans le modèle stratigraphique implicite. Le procédé peut en outre comprendre le calcul d'un ou plusieurs attributs de la pluralité de traces sismiques dans l'intervalle.

Brève description des dessins [0005] Les dessins annexés, qui font partie intégrante de la présente demande, illustrent des modes de réalisation des présents enseignements et, en association avec la description, servent à expliquer les principes des présents enseignements. Sur les figures : [0006] La figure 1 illustre un organigramme d'un procédé de traitement de données sismiques, selon un mode de réalisation.

[0007] La figure 2 illustre une vue d'un modèle d'un volume souterrain, comportant un cube sismique et une fonction stratigraphique implicite/modèle du volume, selon un mode de réalisation.

[0008] La figure 3 illustre une vue agrandie du modèle, représentant une construction d'une trace stratigraphique en combinant deux traces sismiques de part et d'autre d'une faille, selon un mode de réalisation.

[0009] La figure 4 illustre des vues de discordances qui sont prises en compte lors de la construction des traces stratigraphiques, selon un mode de réalisation.

[0010] La figure S illustre un affichage d'un attribut calculé à l'aide d'un mode de réalisation du procédé, selon un mode de réalisation.

[0011] La figure 6 illustre une vue schématique d'un système informatique, selon un mode de réalisation.

Description détaillée [0012] La description détaillée présentée ci-après fait référence aux dessins annexés. Chaque fois que cela est approprié, les mêmes références numériques sont utilisées dans les dessins et la description qui suit pour désigner des parties identiques ou semblables. Bien que plusieurs modes de réalisation et caractéristiques de la présente invention soient décrits ici, des variantes, des adaptations, et d'autres formes de réalisation sont possibles sans que l'on s'écarte de l'esprit et de la portée de la présente invention.

[0013] La figure 1 illustre un organigramme d'un procédé 100 de traitement de données sismiques, selon un mode de réalisation. Dans un mode de réalisation spécifique, le procédé 100 peut faciliter une interprétation des données sismiques, par exemple en tant que partie du processus de détermination de divers attributs sismiques. A titre d'exemple, le procédé 100 peut faciliter l'extraction d'attributs géométriques, de volume et/ou de surface stratigraphiquement cohérents des données sismiques. Dans certains modes de réalisation, les attributs sismiques peuvent être calculés sur la base d'intervalles géologiquement cohérents autour de surfaces d'isovaleurs d'un modèle implicite, comme cela sera décrit plus en détail ci-après. Bien que cela soit décrit de manière générale en termes de données/modèles tridimensionnelles (les isovaleurs étant par exemple des surfaces), il est à noter que des modes de réalisation du procédé 100 peuvent facilement s'appliquer à des données bidimensionnelles, ou de tout autre nombre de dimensions. Par conséquent, les termes “domaine” et “isovaleur”, tels qu'ils sont utilisés ici, ne désignent pas en eux-mêmes un nombre particulier quelconque de dimensions.

[0014] Les attributs peuvent être calculés à l'aide d'un ou plusieurs processeurs d'ordinateurs exécutant des instructions contenues dans un ou plusieurs supports non volatils lisibles par ordinateur. En outre, le procédé 100 peut transformer des données sismiques acquises d'un domaine physique en des attributs qui peuvent être visualisés sur un dispositif d'affichage. Ces exemples, et d'autres encore, d'aspects du procédé 100, selon certains modes de réalisation, ressortiront à la lecture de la description qui suit.

[0015] Se référant de nouveau au mode de réalisation illustré sur la figure 1, le procédé 100 peut comporter la réception des entrées sous la forme d'un modèle stratigraphique implicite, comme indiqué en 102, et des données sismiques représentant un volume ou une section, par exemple un domaine souterrain (ayant un nombre quelconque de dimensions), comme indiqué en 104. Les données sismiques peuvent comprendre des traces sismiques, enregistrées à partir d'ondes sismiques qui se sont propagées à travers le domaine souterrain. Par ailleurs, dans certains cas, les données sismiques peuvent avoir une grille ou un maillage qui leur est associé, par exemple afin de discrétiser les données selon le temps/la profondeur et/ou la position dans le domaine souterrain, pour divers processus d'interprétation, de simulation et/ou de modélisation. Dans d'autres cas, les données sismiques peuvent ne pas être maillées.

[0016] La fonction stratigraphique implicite peut être construite de n'importe quelle manière appropriée et, par exemple, des valeurs déduites à l'aide de la fonction peuvent remplit les cellules appropriées de la grille ou du maillage du cube sismique. A titre d'exemple, on peut utiliser une technique de modélisation implicite qui peut comporter l'utilisation d'une ou plusieurs fonctions implicites. A titre d'exemple, une telle technique peut comporter la représentation des horizons géologiques en trois dimensions à l'aide de surfaces d'isovaleurs spécifiques d'un champ scalaire de propriété (par exemple une fonction implicite) défini sur un maillage d'arrière-plan tridimensionnel. Dans d'autres exemples, le maillage d'arrière-plan peut avoir deux dimensions ou tout autre nombre de dimensions. Dans cet exemple, la continuité de la propriété du champ scalaire peut être régie par la continuité du maillage d'arrière-plan.

[0017] La technique de modélisation implicite peut comporter l'élaboration (ou la réception) d'un maillage d'arrière-plan approprié pour l'interpolation d'une fonction implicite, l'identification d'un ensemble de séquences concordantes du type géologique d'horizons stratigraphiques, et la correction du maillage d'arrière-plan sur lequel l'interpolation est effectuée afin de traiter une première séquence concordante ou entre le traitement de deux séquences concordantes successives. Cette correction peut comprendre la création de sous-volumes dans le maillage d'arrière-plan en le subdivisant par des frontières entre séquences précédemment interpolées, l'identification des sous-volumes correspondant à une séquence concordante "courante", et la limitation d'autres processus d'interpolation et d'extraction de surface d'isovaleurs aux sous-volumes identifiés. En outre, la correction peut comprendre, par exemple, la gestion de l'activité des failles dans un ou plusieurs des sous-volumes identifiés, par exemple, en introduisant et/ou en supprimant une ou plusieurs discontinuités internes dans le maillage d'arrière-plan.

[0018] En outre, la technique de modélisation stratigraphique implicite peut comporter l'interpolation d'une ou plusieurs fonctions implicites "séquence concordante par séquence concordante", par exemple, de façon facultative, une séquence concordante à la fois. Une analyse non limitative supplémentaire de séquences concordantes est fournie ci-après. Dans d'autres modes de réalisation, le domaine sismique peut être traité sur toute autre base, et/ou deux séquences concordantes ou plus peuvent être considérées en parallèle.

[0019] Le processus d'inteipolation peut se dérouler de n'importe quelle manière appropriée, par exemple de manière itérative, comme en tenant compte des effets d'éléments aux limites, etc. Le procédé de modélisation implicite peut fournir en sortie un maillage (par exemple, ou des maillages). A titre d'exemple, un maillage (par exemple, ou des maillages) peut être considéré comme étant un modèle d'un environnement géologique.

[0020] Il est à noter que l'exemple ci-dessus d'une technique de modélisation stratigraphique implicite ne doit pas être considéré comme étant limitatif mais qu'il ne constitue qu'un exemple parmi un grand nombre des techniques de modélisation proposées ici.

[0021] Revenant à la figure 1, le procédé 100 peut comporter la détermination d'une taille d'intervalle, comme indiqué en 106. La taille d'intervalle peut être un intervalle géologiquement cohérent, par exemple de temps, dans le modèle stratigraphique implicite. Plus précisément, dans un mode de réalisation, l'intervalle peut s'étendre au-dessus et en dessous d'une surface d'isovaleurs sélectionnée sur une distance de AS qui peut être une valeur constante (ou sensiblement constante) dans les unités du modèle stratigraphique implicite. A titre d'exemple, les unités AS peuvent être en âge géologique, qu'il soit relatif ou absolu, par exemple en fonction du traitement/marquage du domaine souterrain qui est modélisé.

[0022] La taille de l'intervalle peut être reçue en tant qu'entrée d'un utilisateur, dans au moins un mode de réalisation. Dans d'autres modes de réalisation, la taille de l'intervalle peut être prédéterminée, sélectionnée selon divers attributs globaux des données sismiques (résolution, variance, etc.). Dans encore d'autres modes de réalisation, le procédé 100 peut être itératif, comme cela sera expliqué plus en détail ci-après, et ainsi la taille de l'intervalle peut être sélectionnée sur la base au moins en partie de la taille de l'intervalle lors d'une itération précédente ou suivante.

[0023] Avant, pendant, ou après avoir déterminé la taille d'intervalle, le procédé 100 peut également comporter la définition (par exemple sélection) d'une surface d'isovaleurs du modèle stratigraphique, comme indiqué en 108. La surface d'isovaleurs peut représenter une topographie de surface, telle qu'elle est représentée dans le modèle implicite, ayant une valeur géologique globalement uniforme ou temporelle relative dans le modèle implicite. La surface d'isovaleurs peut être une surface bidimensionnelle, capturant un, quelques, ou la totalité des points positionnés au temps géologique (absolu ou relatif) dans la fonction implicite représentant le domaine souterrain.

[0024] La figure 2 illustre un exemple d'un modèle 200 représentant un cube sismique sur lequel est superposée une fonction stratigraphique implicite (ou du moins chevauchant celle-ci), selon un mode de réalisation. Comme illustré, le modèle 200 comprend des surfaces d'isovaleurs 202, parmi lesquelles une surface d'isovaleurs 204 est sélectionnée en 108. Comme cela pourra être facilement apprécié, le modèle 200 peut être une visualisation affichée sur un dispositif d'affichage informatique.

[0025] Le procédé 100 peut comporter la définition d'un intervalle stratigraphique géologiquement cohérent dans le modèle, sur la base de la position de la surface d'isovaleurs, comme indiqué en 110. A titre d'exemple, l'intervalle peut s'étendre au-dessus et en dessous de la surface d'isovaleurs sélectionnée 204 sur une distance (temps) de grandeur AS. Se référant de nouveau à la figure 2, un intervalle 206 peut être défini sur la base de la position de la surface d'isovaleurs sélectionnée 204, par exemple au-dessus et en dessous de la demi distance de l'intervalle AS. Par conséquent, l'intervalle 206 peut avoir une hauteur de 2AS. Des frontières d'intervalle 208, 210 peuvent également être des surfaces d'isovaleurs, fixées à un AS globalement constant dans la direction verticale de la surface d'isovaleurs sélectionnée 204. Les frontières d'intervalle 208,210 peuvent donc définir entre elles l'intervalle géologiquement cohérent 206.

[0026] Le procédé 100 peut ensuite procéder à la sélection de traces sismiques contenues dans l'intervalle stratigraphique, comme indiqué en 112. Dans certains modes de réalisation, les données de l'intervalle peuvent comprendre des données autres que des traces sismiques, ou s'ajoutant à celles-ci, telles que des attributs sismiques (par exemple l'amplitude, la fréquence, le temps, l'atténuation, la variance, le maximum, le minimum, ou l'amplitude moyenne, etc.). Par conséquent, bien que l'exemple illustré fasse référence au fait de sélectionner des traces sismiques, il est à noter que des modes de réalisation dans lesquels d'autres types de données, qui peuvent le cas échéant être de type sismique, peuvent être envisagés conformément à la présente invention.

[0027] Le procédé 100 peut également comporter la définition des “traces stratigraphiques” sur la base des traces sismiques, comme indiqué en 114. A titre d'exemple, lesdites une ou plusieurs traces sismiques du cube sismique peuvent être sélectionnées par déplacement perpendiculaire à, ou le long d'un gradient de la fonction stratigraphique dans l'intervalle (par exemple l'intervalle 206). Dans divers modes de réalisation, cela peut comporter l'échantillonnage d'uns ou plusieurs (par exemple plusieurs en trois dimensions) traces sismiques dans le cube sismique, et la connexion des parties des traces sismiques les unes aux autres.

[0028] Dans certains modes de réalisation, par exemple en tant que partie de la définition des traces stratigraphiques, le procédé 100 peut comporter une opération pour déterminer si une ou plusieurs des traces sismiques passent par une discontinuité (par exemple une faille), comme indiqué en 116. Cette détermination peut être effectuée selon une mise en correspondance de failles, dont les résultats peuvent être fournis en tant que partie du modèle stratigraphique implicite. En réponse à une détermination positive en 116, le procédé 100 peut, en 118, comporter le pontage des deux côtés de la faille. A titre d'exemple, le pontage en 118 peut comporter la sélection d'une trace sismique provenant des deux côtés de la faille. Par conséquent, à titre d'exemple, lors du calcul d'une racine quadratique moyenne d'amplitude pour une trace stratigraphique, la trace stratigraphique peut prendre en compte les deux parties de traces sismiques provenant des deux côtés de la faille.

[0029] La sélection des parties de traces sismiques à assembler et ainsi former la trace stratigraphique à travers une discontinuité (par exemple le pontage indiqué en 118) peut se dérouler de plusieurs manières différentes. A titre d'exemple, des paléo-coordonnées (u, v) peuvent être utilisées, un rejet de faille le plus petit peut être trouvé, et/ou des orientations définies par l'utilisateur peuvent être utilisées.

[0030] Plus précisément, des paléo-coordonnées peuvent par exemple être calculées sous la forme d'un processus distinct et peuvent permettre de parcourir une surface isostratigraphique, en tenant compte du déplacement de faille (par exemple, l'espace sédimentaire peut générer une représentation chronostratigraphique pseudo-3D d'un modèle structurel, dans lequel une méthode des éléments finis peut faciliter le rétro-calcul du déplacement de faille se produisant au cours du temps). En outre, le rejet le plus petit pour un point peut se situer au niveau du toit de faille/paroi inférieure du contact de faille d'horizon, qui peut être la plus faible distance entre ce point et le toit de faille/paroi inférieure correspondant du même contact de faille d'horizon. De plus, un utilisateur peut spécifier un “azimut” unique pour la totalité du cube, qui peut représenter l'axe principal d'une formation de faille (en extension ou en compression). De ce fait, une variété de facteurs de ce type peuvent être utilisés, seuls ou en combinaison, de manière à sélectionner les parties de traces sismiques à assembler lors de la formation de la trace stratigraphique.

[0031] La figure 3 illustre une vue agrandie d'un partie 300 du modèle 200, illustrant de manière conceptuelle un exemple de pontage d'une faille 302 qui est rejointe par deux traces sismiques 304, 306, selon un mode de réalisation. Comme illustré, la faille 302 passe par l'intervalle 206, cela conduisant aux deux traces sismiques 304, 306 ayant des trajectoires différentes dans l'intervalle 206. Lorsque les deux traces 304, 306 sont sélectionnées pour l'opération de pontage, par exemple sur la base des facteurs analysés plus haut, les traces 304, 306 peuvent être considérées ensemble comme étant une trace stratigraphique unique 308, par exemple à des fins de détermination d'attributs sismiques (ou autres) dans l'intervalle 206.

[0032] Se référant de nouveau à la figure 1, le procédé 100 peut également comporter une opération pour déterminer si les traces sismiques dans l'intervalle 206 franchissent une limite de séquence ou une "discordance”, comme indiqué en 120. Si cela est le cas, le procédé 100 peut répondre en ajustant l'intervalle 206 par rapport à la trace sismique passant par la discordance, de manière à ne plus prendre en compte au moins une partie de la trace sismique lors de la construction de la trace stratigraphique. En d'autres termes, une partie de la trace sismique peut être ignorée si elle se trouve du "mauvais" côté de la discordance (à titre d'exemple, dans le cas de la discordance entre la surface d'isovaleurs sélectionnée et la partie de la trace sismique). Ainsi, une cohérence géologique peut être préservée dans les calculs d'attributs sismiques, selon certains modes de réalisation, étant donné que les roches se trouvant sur des côtés différents d'une discordance peuvent présenter une cohérence géologique faible ou nulle les unes par rapport aux autres.

[0033] La figure 4 illustre quatre types de séquences pouvant être prises en compte, selon un mode de réalisation. Plus précisément, celle-ci illustre des exemples de formations qui comprennent une ou plusieurs séquences, par exemple, des séquences de structures sédimentaires (par exemple, des strates, des horizons, etc.) se produisant dans des roches sédimentaires. Comme illustré, une formation 410 peut comprendre une séquence unique, des formations 420 et 430 peuvent chacune comprendre deux séquences, et une formation 440 peut comprendre trois séquences, la séquence médiane étant réduite à une surface de discontinuité unique.

[0034] A titre d'exemple, un horizon concordant peut être un horizon entre un horizon inférieur et un horizon supérieur où les horizons ont été soumis à une histoire géologique relativement commune, par exemple, en étant déposés successivement (par exemple, en continu dans le temps). Se référant à la formation 410, les horizons ne se coupent pas les uns les autres, et chacun des horizons peut être considéré comme étant concordant avec des horizons adjacents (par exemple, inférieur et supérieur ou plus ancien et plus jeune).

[0035] A titre d'exemple, l'érosion peut avoir pour effet de dénuder la roche, par exemple, du fait de la décomposition et/ou du transport physique, chimique et/ou biologique. L'érosion peut par exemple se produire lorsque du matériau (par exemple altéré par la roche, etc.) est transporté par des fluides, des solides (par exemple le vent, l'eau ou la glace) ou par mouvement de masse (par exemple, comme dans le cas de chutes de pierres et de glissements de terrain). Se référant à la formation 420, parmi les deux séquences représentées, la séquence inférieure peut avoir été érodée et la séquence supérieure avoir été déposée sur le dessus de la séquence inférieure érodée. Dans cet exemple, la frontière entre les deux séquences peut être désignée sous le nom d'érosion ; en notant qu'elle est concordante avec la séquence supérieure, plus jeune. A titre d'exemple, l'érosion peut avoir pour effet de “tronquer” une séquence d'horizons et former une surface sur laquelle les matériaux suivants peuvent se déposer (par exemple, le cas échéant, de manière concordante).

[0036] A titre d'exemple, un recouvrement de base peut être un type de caractéristique dans une formation, comme par exemple un biseau de progradation ou un biseau d'aggradation. A titre d'exemple, un biseau de progradation peut être une terminaison de strates susjacentes plongeant plus profondément contre une surface ou des strates sous-jacentes qui ont des pendages apparents inférieurs. A titre d'exemple, un biseau de progradation peut être observé à la base de clinoformes progradantes et peut représenter la progradation d'une lisière de bassin. Dans le cas d'un biseau d'aggradation, il peut par exemple s'agir d'une terminaison de strates plus jeunes à pendage faiblement profond contre des strates plus anciennes à pendage plus profond (par exemple, une stratigraphie de séquences qui peut se produire pendant des périodes de transgression). Se référant à la formation 430, la direction indiquée par “z” étant la profondeur, le type de recouvrement de base représenté peut être considéré comme étant un biseau de progradation (par exemple, des strates inférieures ayant des pendages apparents inférieurs). Dans cet exemple, la frontière du recouvrement de base est concordante avec des horizons immédiatement plus anciens (séquence inférieure).

[0037] Quant à la formation 440, celle-ci comprend trois séquences et peut être désignée comme une discontinuité, étant donné que la frontière n'est concordante ni avec des horizons plus anciens ni avec des horizons plus jeunes. Dans les exemples de la figure 4, les érosions, les recouvrements de base et les discontinuités peuvent être désignés comme des discordances ou des horizons (ou par exemple de surfaces, de couches, etc.) non concordants.

[0038] A titre d'exemple, lorsqu'il est déterminé, en tant que partie du procédé 100, qu'une trace stratigraphique (ou qu'une trace sismique qui doit être reconstruite dans une partie d'une trace stratigraphique) passe par une érosion, le procédé 100 peut contraindre l'intervalle 206 à la position de la surface d'isovaleurs 204 plus AS, par exemple en divisant en deux l'intervalle 206 pour la trace stratigraphique. Dans un autre mode de réalisation, le procédé 100 peut comporter le décalage de l'intervalle 206 vers le haut. En outre, la taille et/ou la position de l'intervalle 206 peut être ajustée par toute autre méthode de manière à éviter la prise en compte de deux côtés de l'érosion 402.

[0039] Lorsque le procédé 100 inclut la détermination de la trace stratigraphique passant par un recouvrement de base, le procédé 100 peut de façon similaire contraindre l'intervalle 206, mais, dans ce cas, à la position de la surface d’isovaleurs 202 moins AS (ou un certain multiple de AS). Cela peut être frit pour des raisons opposées au décalage (et/ou sinon ajustement) de l'intervalle 206 lorsque la trace stratigraphique passe par l'érosion 402. Plus précisément, la roche située au-dessus du recouvrement de base peut ne pas être cohérente ou pertinente par rapport à la roche située au-dessus.

[0040] Lorsque le procédé 100 détermine que la trace stratigraphique passe par une discontinuité à trois séquences, le procédé 100 peut omettre la trace stratigraphique de toute autre prise en compte d'attributs sismiques volumétriques dans l'intervalle 206. La trace stratigraphique située au-dessus de la roche peut être peu liée ou ne pas être liée à la roche située en-dessous, et inversement, de telle sorte que les attributs volumétriques concernant la discontinuité à trois séquences peut ne pas être calculée.

[0041] Pour l'intervalle stratigraphique défini selon un intervalle géologiquement cohérent 206, par exemple centré ou positionné d'une autre manière sur la surface d'isovaleurs sélectionnée 204, et les traces stratigraphiques se situant dans l'intervalle substitué à des traces sismiques verticales ou orientées suivant le pendage, le procédé 100 peut procéder au calcul d'un ou plusieurs attributs (par exemple sismiques, géométriques, superficiels, etc.) dans l'intervalle 206 sur la base des traces stratigraphiques, par exemple en trois dimensions, comme indiqué en 124. De plus, le procédé 100 peut comporter l'affichage des attributs calculés, par exemple à l'aide d'un écran d'affichage (par exemple un moniteur, un écran tactile, un projecteur, etc.). Les attributs affichés peuvent représenter une conclusion différente qui peut être tirée du procédé 100 lorsqu'on utilise des intervalles géologiquement cohérents pour déterminer les attributs. La figure 5 illustre un exemple d'un attribut calculé sur la base d'intervalles de temps géologiquement cohérents.

[0042] Le procédé 100 peut en outre être itératif. Par conséquent, le procédé 100 peut déterminer si des intervalles supplémentaires doivent être pris en compte, comme indiqué en 126. Si des intervalles supplémentaires doivent être déterminés, le procédé 100 peut revenir à la sélection d'une isovaleur en 108, qui peut être différente de la surface d'isovaleurs précédemment définie. Ainsi, le procédé 100 peut parcourir les intervalles du modèle 200 jusqu'à ce que la totalité ou un volume d'intérêt de ceux-ci soit pris en compte. En outre, lors du calcul de certains attributs, les attributs calculés pour des intervalles voisins peuvent affecter le calcul d'au moins une partie des attributs des autres. Par conséquent, le procédé 100 peut effectuer une itération pour calculer une ou plusieurs fois l'attribut en 124 en utilisant n'importe lesquels ou la totalité des intervalles du modèle 200 (ou une partie quelconque de ceux-ci).

[0043] Le procédé 100 peut en outre comporter un ajustement de la taille d'intervalle, comme indiqué en 128, ayant pour effet de modifier la taille de la fenêtre pour les attributs sismiques fenêtrés. Cela peut autoriser des déterminations d'attributs plus précises ou plus généralisées.

[0044] Par ailleurs, lorsqu'on utilise le volume de fonction stratigraphique, des modes de réalisation du procédé 100 peuvent permettre un accès direct aux informations topographiques. Λ titre d'exemple, une trace stratigraphique échantillonnée quelconque dans le volume peut avoir une valeur attribuée dans le cube stratigraphique, permettant une extraction de la trace stratigraphique dans le cube avec la même valeur afin de recréer une surface d'âge géologique commun. Celle-ci peut être corrélée à un événement sismique. Par conséquent, le procédé 100 peut comporter le calcul des “attributs de surface”, directement en utilisant le volume sismique visé de manière combinée avec la fonction stratigraphique. L'“attribut de surface” peut ensuite être stocké sous la forme d'un nouveau domaine sismique (par exemple, un volume). Les attributs de surface peuvent ensuite être extraits et à nouveau stockés comme un enfant d'une surface, ou soumis à un rendu combiné sur une isovaleur de la fonction stratigraphique.

[0045] Le procédé 100 peut contenir de l'énergie à l'intérieur d'un événement stratigraphique, évitant d'étaler l’énergie à travers les discontinuités (failles ou discordances stratigraphiques). Ces considérations peuvent faciliter les processus de lissage et/ou d'annulation du bruit.

[0046] Dans certains modes de réalisation, le procédé 100 (et/ou des processus quelconques de celui-ci) peut être exécuté par un système informatique. La figure 6 illustre un exemple d'un tel système informatique 600, conformément à certains modes de réalisation. Le système informatique 600 peut comprendre un ordinateur ou un système d'ordinateur 601 A, qui peut être un système d'ordinateur individuel 601A ou un agencement de systèmes d'ordinateurs répartis. Le système d'ordinateur 601A comprend un ou plusieurs modules d'analyse 602 qui sont configurés pour exécuter diverses tâches selon certains modes de réalisation, telles qu'un ou plusieurs des procédés décrits ici (par exemple, les procédés 100, 500, et/ou des combinaisons et/ou des variantes de ceux-ci). Pour exécuter ces diverses tâches, le module d'analyse 602 s'exécute indépendamment sur, ou en coordination avec, un ou plusieurs processeurs 604, qui est (ou sont) connecté(s) à un ou plusieurs supports de stockage 606A. Le(s) processeurs) 604 est (ou sont) également connecté(s) à une interface de réseau 607 pour permettre au système d'ordinateur 601A de communiquer par l'intermédiaire d'un réseau de données 609 avec un ou plusieurs systèmes d'ordinateurs et/ou systèmes informatiques supplémentaires, tels que les systèmes informatiques 601 B, 60IC, et/ou 601D (on notera que les systèmes d'ordinateurs 601B, 601C, et/ou 601D peuvent ou non partager la même architecture que le système d'ordinateur 601A, et peuvent être localisés à des emplacements physiques différents, les systèmes d'ordinateurs 601A et 601B pouvant par exemple être localisés dans une installation de traitement, tout en étant en communication avec un ou plusieurs systèmes d'ordinateurs tels que 601C et/ou 601D qui sont localisés dans un ou plusieurs centres de données, et/ou localisés dans divers pays sur différents continents).

[0047] Un processeur peut comprendre un microprocesseur, un microcontrôleur, un module ou un sous-système processeur, un circuit intégré programmable, une matrice prédifiusée programmable, ou un autre dispositif de commande ou informatique.

[0048] Les supports de stockage 606A peuvent être mis en œuvre sous la forme d'un ou plusieurs supports de stockage lisibles par ordinateur ou lisibles par machine. On notera que, bien que, dans l'exemple de mode de réalisation de la figure 6, le support de stockage 606A est représenté comme étant au sein du système d'ordinateur 601A, dans certains modes de réalisation, le support de stockage 606A peut être réparti à l'intérieur et/ou à travers de multiples enceintes internes et/ou externes du système informatique 601A et/ou de systèmes informatiques supplémentaires. Le support de stockage 606A peut comprendre une ou plusieurs formes de mémoire différentes, incluant des dispositifs de mémoire à semi-conducteur tels que des mémoires vives dynamiques ou statiques (DRAM ou SRAM), des mémoires mortes effaçables et programmables (EPROM), des mémoires mortes effaçables et programmables électriquement (EEPROM) et des mémoires flash, des disques magnétiques tels que des disques fixes, souples et amovibles, d'autres supports magnétiques, incluant des supports optiques à bande tels que des disques compacts (CD) ou des disques vidéo numériques (DVD), des disques BLUERAY®, ou d'autres types de stockage optique, ou d'autres types de dispositifs de stockage. On notera que les instructions discutées plus haut peuvent être fournies sur un support de stockage lisible par ordinateur ou lisible par machine, ou en variante, peuvent être fournies sur de multiples supports de stockage lisibles par ordinateur ou lisibles par machine répartis dans un large système pouvant comporter une pluralité de nœuds. Ce ou ces support(s) de stockage lisible(s) par ordinateur ou lisible(s) par machine est (sont) considéré(s) comme faisant partie d'un article (ou d'un article de fabrication). Un article ou article de fabrication peut être n'importe quel(s) composant unique ou composants multiples fabriqués. Le ou les supports) de stockage peut ou peuvent être localisé(s) soit dans la machine exécutant les instructions lisibles par machine, soit être localisé(s) sur un site distant depuis lequel des instructions lisibles par machine peuvent être téléchargées par l'intermédiaire d'un réseau à des fins d'exécution.

[0049] Dans certains modes de réalisation, le système informatique 600 contient un ou plusieurs module(s) de sélection de modèle 608. Dans l'exemple de système informatique 600, le système d'ordinateur 60 IA comprend un module de modélisation 608. Dans certains modes de réalisation, un module de modélisation unique peut être utilisé pour exécuter certains ou la totalité des aspects d'un ou plusieurs des modes de réalisation du procédé 100. Dans d'autres modes de réalisation, une pluralité de modules de modélisation peuvent être utilisés pour exécuter certains ou la totalité des aspects du procédé 100.

[0050] Il est à noter que le système informatique 600 n'est qu'un exemple d'un système informatique, et que le système informatique 600 peut avoir plus ou moins de composants que ceux représentés, qu'il peut combiner des composants supplémentaires non représentés dans l'exemple de mode de réalisation de la figure 6, et/ou que le système informatique 600 peut avoir une configuration différente ou un agencement différent des composants représentés sur la figure 6. Les divers composants représentés sur la figure 6 peuvent être mis en œuvre sous forme matérielle, logicielle, ou d'une combinaison de des deux formes matérielle et logicielle, y compris un ou plusieurs circuits intégrés de traitement du signal et/ou spécifiques d'applications.

[0051] En outre, les divers aspects des procédés de traitement décrits ici peuvent être mis en œuvre par exécution d'un ou plusieurs modules fonctionnels dans des appareils de traitement d'informations tels que des processeurs polyvalents ou des puces spécifiques d'applications, tels que des ASIC, des FPGA, des PLD, ou d'autres dispositifs appropriés. Ces modules, des combinaisons de ces modules, et/ou leur combinaison avec des matériels généraux entrent tous dans la portée de protection de l'invention.

[0052] Les divers aspects des procédés de traitement décrits ici peuvent être mis en œuvre par exécution d'un ou plusieurs modules fonctionnels dans des appareils de traitement d'informations tels que des processeurs polyvalents ou des puces spécifiques d'applications, tels que des ASIC, des FPGA, des PLD, ou d'autres dispositifs appropriés. Ces modules, des combinaisons de ces modules, et/ou leur combinaison avec des matériels généraux, entrent tous dans la portée de protection de l'invention.

[0053] Par ailleurs, les interprétations géologiques, les modèles et/ou d'autres aides à l'interprétation peuvent être affinés de façon itérative ; ce concept peut s'appliquer au procédé 100, 500, tel qu'il est analysé ici. Cela peut comporter l'utilisation de boucles de rétroaction exécutées de manière algorithmique, comme un dispositif informatique (par exemple le système informatique 600, figure 6), et/ou par commande manuelle par un utilisateur qui peut effectuer des déterminations pour savoir si une étape donnée, une action donnée, un gabarit donné, un modèle donné, ou un ensemble de courbes donné est devenu suffisamment précis pour l'évaluation de la formation géologique tridimensionnelle souterraine considérée.

[0054] La description présentée ci-dessus à titre d'explication a été fournie en référence à des modes de réalisation. Cependant, les analyses présentées ci-dessus à titre d'illustration ne doivent pas être considérées comme étant exhaustives ou comme limitant l'invention aux formes précises mentionnées. De nombreuses modifications et variantes sont possibles au vu des enseignements présentés plus haut. A titre d'exemple, l'ordre dans lequel les aspects du procédé 100 sont mis en œuvre peut être défini à l'avance, sans que l'on s'écarte de la portée de la présente invention. Les modes de réalisation ont été choisis et décrits afin d'expliquer au mieux les principes de l'invention et ses applications pratiques, pour ainsi permettre à l'homme du métier d'exploiter au . mieux l'invention et ses divers modes de réalisation avec diverses variantes qui sont appropriées pour l'utilisation particulière envisagée.

Claims

REVENDICATIONS

1. Procédé de traitement de données sismiques, comportant les opérations pour : recevoir (104) une pluralité de traces sismiques représentant un domaine souterrain ; recevoir (102) un modèle stratigraphique implicite d'au moins une partie du domaine souterrain ; sélectionner (108) une isovaleur dans le modèle stratigraphique implicite ; définir (110), à l'aide d'un processeur, un intervalle géologiquement cohérent dans le modèle stratigraphique implicite sur la base au moins en partie d'une position de l'isovaleur dans le modèle stratigraphique implicite ; et calculer (124) un ou plusieurs attributs de la pluralité de traces sismiques dans l'intervalle.
2. Procédé selon la revendication 1, dans lequel la définition de l'intervalle géologiquement cohérent comporte : un réglage d'une première frontière décalée par rapport à l'isovaleur d'un premier temps géologique ; et un réglage d'une seconde frontière décalée par rapport à l'isovaleur d'un second temps géologique, dans lequel l'intervalle est défini entre les première et seconde frontières.
3. Procédé selon la revendication 2, dans lequel une grandeur du premier temps géologique est constante à travers la première frontière et est égale à une grandeur du second temps géologique.
4. Procédé selon l'une quelconque des revendications précédentes, comportant en outre une construction d'une ou plusieurs traces stratigraphiques orientées perpendiculairement à ou le long d'un gradient d'une fonction stratigraphique sur laquelle est fondé le modèle stratigraphique implicite, dans lequel lesdites une ou plusieurs traces stratigraphiques sont construites sur la base d'une ou plusieurs de la pluralité de traces sismiques, et dans lequel lesdits un ou plusieurs attributs sont calculés sur la base desdites une ou plusieurs traces stratigraphiques.
5. Procédé selon la revendication 4, dans lequel la construction desdites une ou plusieurs traces stratigraphiques comporte des opérations pour: déterminer qu'une première trace sismique de la pluralité de traces sismiques et qu'une seconde trace sismique de la pluralité de traces sismiques passent par une discontinuité ; et en réponse, ponter la discontinuité en construisant lesdites une ou plusieurs traces stratigraphiques à partir d'une combinaison des première et seconde traces sismiques, dans lequel les première et seconde traces sismiques ont des trajectoires différentes.
6. Procédé selon la revendication 5, dans lequel le pontage de la discontinuité comporte une sélection des première et seconde traces sismiques pour les combiner en lesdites une ou plusieurs traces stratigraphiques sur la base de paléo-coordonnées, d'un rejet de faille, d'une orientation définie par l'utilisateur, ou d'une combinaison de ceux-ci.
7. Procédé selon la revendication 4, dans lequel la construction desdites une ou plusieurs traces stratigraphiques comporte des opérations pour : déterminer qu'au moins une de la pluralité de traces sismiques passe par une discordance dans l'intervalle ; et ajuster l'intervalle par rapport à la au moins une de la pluralité de traces sismiques, de manière à ce que la au moins une de la pluralité de traces sismiques soit seulement considérée sur un côté de la discordance.
8. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre un affichage d'un ou plusieurs attributs au moyen d'un dispositif d'affichage.
9. Procédé selon l'tme quelconque des revendications précédentes, comprenant en outre des opérations pour : stocker les un ou plusieurs attributs sous la forme d'un nouveau domaine sismique ; et extraire les un ou plusieurs attributs comme un enfant de l'isovaleur, ou afficher les un ou plusieurs attributs sur une isovaleur du modèle stratigraphique implicite, ou les deux.
10. Système informatique, comprenant : un ou plusieurs processeurs ; et un système de mémoire comprenant un ou plusieurs supports non volatils lisibles par ordinateur stockant des instructions qui, lorsqu'elles sont exécutées par au moins l'un des un ou plusieurs processeurs, amènent le système informatique à exécuter des opérations, les opérations comprenant : une réception (104) d'une pluralité de traces sismiques représentant un domaine souterrain ; une réception (102) d'un modèle stratigraphique implicite d'au moins une partie du domaine souterrain ; une sélection (108) d'une isovaleur dans le modèle stratigraphique implicite ; une définition (110) d'un intervalle géologiquement cohérent dans le modèle stratigraphique implicite sur la base au moins en partie d*une position de l'isovaleur dans le modèle stratigraphique implicite ; et un calcul (124) d'un ou plusieurs attributs de la pluralité de traces sismiques dans l'intervalle.
11. Système selon la revendication 10, dans lequel la définition de l'intervalle géologiquement cohérent comporte : Un réglage d'une première frontière décalée par rapport à l'isovaleur d'un premier temps géologique ; et un réglage d'une seconde frontière décalée par rapport à l'isovaleur d'un second temps géologique, dans lequel l'intervalle est défini entre les première et seconde frontières.
12. Système selon la revendication 11, dans lequel une grandeur du premier temps géologique est constante à travers la première frontière et est égale à une valeur du second temps géologique.
13. Système selon l'une quelconque des revendications 10-12, dans lequel les opérations comprennent en outre une construction d'une ou plusieurs traces stratigraphiques orientées perpendiculairement à ou le long d'un gradient d'une fonction stratigraphique sur laquelle est basé le modèle stratigraphique implicite, dans lequel lesdites une ou plusieurs traces stratigraphiques sont construites sur la base d'une ou plusieurs de la pluralité de traces sismiques, et dans lequel les un ou plusieurs attributs sont calculés sur la base des une ou plusieurs traces stratigraphiques.
14. Système selon la revendication 13, dans lequel la construction desdites une ou plusieurs traces stratigraphiques comporte : une détermination d'une première trace sismique de la pluralité de traces sismiques et d'une seconde trace sismique de la pluralité de traces sismiques qui passent par une discontinuité ; et en réponse, un pontage de la discontinuité en construisant les une ou plusieurs traces stratigraphiques à partir d'une combinaison des première et seconde traces sismiques, dans lequel les première et seconde traces sismiques ont des trajectoires différentes.
15. Système selon la revendication 13, dans lequel la construction des une ou plusieurs traces stratigraphiques comprennent : une détermination d'au moins une de la pluralité de traces sismiques qui passe par une discordance dans l'intervalle ; et un ajustement de l'intervalle par rapport à la au moins une de la pluralité de traces sismiques, de manière à ce que la au moins une de la pluralité de traces sismiques soit seulement considérée sur un côté de la discordance.
16. Support non volatil lisible par ordinateur stockant des instructions qui, lorsqu'elles sont exécutées par un ou plusieurs processeurs d'un système informatique, amènent le système informatique à exécuter des opérations, les opérations comprenant : une réception (104) d'une pluralité de traces sismiques représentant un domaine souterrain ; une réception (102) d'un modèle stratigraphique implicite d'au moins une partie du domaine souterrain ; une sélection (108) d'une isovaleur dans le modèle stratigraphique implicite ; une définition (110) d'un intervalle géologiquement cohérent dans le modèle stratigraphique implicite sur la base au moins en partie d'une position de l'isovaleur dans le modèle stratigraphique implicite ; et un calcul (124) d'un ou plusieurs attributs de la pluralité de traces sismiques dans l'intervalle.
17. Support selon la revendication 16, dans lequel la définition de l'intervalle géologiquement cohérent comporte : un réglage d'une première frontière décalée par rapport à l'isovaleur d'un premier temps géologique ; et un réglage d'une seconde frontière décalée par rapport à l'isovaleur d'un second temps géologique, dans lequel l'intervalle est défini entre les première et seconde frontières.
18. Support selon la revendication 16, dans lequel les opérations comprennent en outre une construction d'une ou plusieurs traces stratigraphiques orientées perpendiculairement à ou le long d'un gradient d'une fonction stratigraphique sur laquelle est basé le modèle stratigraphique implicite, dans lequel les une ou plusieurs traces stratigraphiques sont construites sur la base d'une ou plusieurs de la pluralité de traces sismiques, et dans lequel les un ou plusieurs attributs sont calculés sur la base des une ou plusieurs traces stratigraphiques.
19. Support selon la revendication 18, dans lequel la construction des une ou plusieurs traces stratigraphiques comporte : une détermination d'une première trace sismique de la pluralité de traces sismiques et d'une seconde trace sismique de la pluralité de traces sismiques passent par une discontinuité ; et en réponse, un pontage de la discontinuité en construisant les une ou plusieurs traces stratigraphiques à partir d'une combinaison des première et seconde traces sismiques, dans lequel les première et seconde traces sismiques ont des trajectoires différentes.
20. Support selon la revendication 18, dans lequel la construction desdites une ou plusieurs traces stratigraphiques comporte : une détermination d'au moins une de la pluralité de traces sismiques qui passe par une discordance dans l'intervalle ; et un ajustement de l'intervalle par rapport à la au moins une de la pluralité de traces sismiques, de manière à ce que la au moins une de la pluralité de traces sismiques soit seulement considérée sur un côté de la discordance.