Brevet FR2884636A1 — publié le 20061020

qui représentent généralement des formes énergétiques et continues de l'image. Par définition, on appellera événement cohérent dans une image 3D toute surface cohérente suivant des critères de continuité et d'énergie. Pour une image 2D, un événement cohérent correspondra à une ligne non nécessairement droite. Les procédés habituels de détermination d'évènements cohérents dans des images sismiques sont généralement basés sur des algorithmes dits de propagation. On place des germes sur l'image, c'est-à-dire des points d'ancrage au niveau des évènements cohérents que l'on veut déterminer. Ces germes sont généralement le résultat d'une interprétation humaine. Les algorithmes de propagation déterminent l'intégralité de chacun des évènements cohérents choisis en recherchant dans l'image les chemins de plus grande corrélation spatiale à partir des germes relatifs à chaque évènement. Ils propagent de point en point, pixel par pixel, les évènements cohérents de l'image. Ce type d'algorithme peut se révéler instable, particulièrement en environnement bruité, et un mauvais chemin de plus grande corrélation spatiale est parfois vite emprunté. De plus, cette approche par propagation n'est pas optimale pour une interprétation volumique 3D et en temps réel des évènements cohérents présents dans des cubes sismiques. Une application directe de l'invention concerne le contrôle qualité d'un pointé d'horizons (évènements cohérents correspondant à des interfaces géologiques) sur un cube sismique. Le pointé d'évènements cohérents appelés horizons sur des cubes sismiques est réalisé de façon semi-automatique. Un interpréteur crée des germes qui servent de points d'ancrage pour un algorithme de propagation destiné à générer tout ou partie de l'horizon pointé . Mais, ce type d'algorithme peut se révéler instable, particulièrement en environnement bruité. Cette instabilité se traduit par des décalages du pointé par rapport à 2884636 -3- l'extremum local correspondant à l'événement sismique ciblé, ou bien parfois par des erreurs plus radicales correspondant à des sauts de phase. Ces erreurs peuvent avoir un impact très important sur la valeur de l'amplitude ou de n'importe quel autre attribut extraits par exemple dans le cadre d'une étude orientée réservoir . Or, un contrôle fiable et rapide de la qualité du pointé de ces horizons se révèle indispensable. En effet, une étude orientée réservoir est une étude d'imagerie d'un réservoir et d'estimation des réserves qu'il contient et les enjeux économiques induits par de telles études étant considérables, la précision des résultats qu'elles fournissent se doit d'être optimale. De petites erreurs peuvent avoir de grandes conséquences économiques. Par ailleurs, l'examen visuel est mal approprié pour repérer rapidement des anomalies spatiales s'étendant dans les trois dimensions de l'espace. Une autre application possible concerne la détermination d'évènements cohérents sur différents cubes sismiques d'une même zone géographique, ces cubes étant générés avec des paramètres de traitement différents. Dans les solutions actuellement proposées, il arrive couramment en traitement et en interprétation sismique que deux (ou plusieurs) cubes sismiques soient définis sur une même zone géographique. Dépendant du contexte opérationnel, l'analyse des différences entre les deux cubes est alors généralement instructive. Elle peut mettre en évidence des différences liées: - à la production d'hydrocarbures dans le temps (problématique réservoir 4D c'est à dire défini selon les trois dimensions spatiales et une dimension temporelle) ; - à des natures de variable étudiée différentes comme dans une problématique 4C où on propage deux types d'ondes sismiques dans le sous-sol: des ondes de compression P et des ondes de cisaillement S afin de réaliser deux cubes, un cube en amplitudes P et un cube en amplitudes S; - à des paramètres de traitement différents. - 4 La caractérisation de ces différences (localisation, quantification) n'est jamais aisée, la simple différence entre les deux cubes sismiques se révélant souvent insuffisante, voire stérile. En effet, des décalages géographiques mêmes faibles des évènements cohérents apparentés sur chacun des deux cubes sismiques diminuent notablement la valeur de l'information du cube différence. En vue de minimiser ces inconvénients, l'invention propose de fournir un outil d'aide à la détermination d'évènements cohérents dans une image sismique en associant à chaque pixel de l'image un couple d'indices EII Event Importance Index , ECI Event Confidence Index calculés en appliquant, sur des fenêtres se recouvrant , une segmentation hiérarchique. Par fenêtre se recouvrant , on entend qu'au lieu de réaliser une segmentation sur toute l'image comme cela se fait couramment, on réalise des segmentations successives sur des parties de l'image définies par une fenêtre glissante qui se déplace sur l'image à la manière d'un scanner. Le calcul du couple d'indices EII et ECI résout le problème de la 20 détermination en temps réel des évènements cohérents d'une image sismique. Avantageusement, la détermination des évènements cohérents d'une image sismique se fait grâce au procédé selon l'invention qui comprend les phases 25 opératoires suivantes: une phase de choix: o d'un critère de segmentation fonction d'une variable représentée sous forme d'image sismique tel que l'amplitude, le 30 contraste, le gradient, o d'une zone de l'image à qualifier, o d'au moins une fenêtre glissante sur ladite zone et de ses caractéristiques, une phase de segmentation hiérarchique, pour des positions de la fenêtre glissante se recouvrant, comportant une segmentation de ladite zone en n régions et pour chaque pixel repéré au moins une fois par une frontière de segmentation: o la sauvegarde du nombre de fois où le pixel est repéré par une frontière de segmentation pour chaque fenêtre glissante à laquelle il appartient et le calcul du rapport de ce nombre par le nombre théorique maximal de fois où ce même pixel pourrait être repéré, ce rapport correspondant à l'indice ECI; o l'attribution pour chaque frontière de segmentation d'un niveau hiérarchique basé sur un calcul fonction d'un critère déterminé tel que la surface de la frontière, de manière à obtenir l'indice EII; - une phase d'attribution aux pixels non repérés d'une valeur d'indices El et EIC correspondant à une non-valeur caractéristique numérique ou alphanumérique. - une phase de détermination des évènements cohérents de l'image par seuillage des indices EII et ECI c'est-à-dire en ne sélectionnant que des pixels repérés de la zone correspondant à des valeurs inférieures ou supérieures à un seuil fixé. La variable pourra être un attribut sismique tel que l'amplitude, la réflectivité, 25 l'impédance ou tout autre attribut calculé sur un cube sismique. Ce seuillage pourra dépendre notamment de l'application métier poursuivie. Les caractéristiques de la fenêtre glissante pourront comprendre ses dimensions, le taux de recouvrement d'une position i de la fenêtre par rapport à une position i-r, le taux de recouvrement étant défini par 1/r, r étant le module du vecteur déplacement dans une des directions principales de l'image sismique multiple d'un pas de grille (définition une image est définie sur une grille avec comme paramètre le pas de grille dans chacune des directions principales de l'image (généralement axes euclidiens)) dans cette direction. Le nombre n de régions de segmentation de la zone pourra être dépendant 5 de la zone segmentée. Les frontières de la segmentation pour chaque position de la fenêtre glissante définissent les pixels repérés par le processus de segmentation et pour lesquels les indices EII et ECI peuvent être calculés, les pixels de faible EII et de fort ECI définissant les évènements les plus cohérents de l'image. Avantageusement, la phase de segmentation pourra s'effectuer: par l'algorithme de segmentation hiérarchique de Ligne de Partage des Eaux ; il s'agit d'un algorithme de morphologie mathématique qui permet de réaliser de la segmentation hiérarchique d'images; - sur une image sismique préalablement filtrée c'est-à-dire une image transformée par un filtre, par exemple un filtre permettant d'améliorer la continuité de l'image et d'éliminer du bruit aléatoire haute fréquence Le recouvrement des positions de la fenêtre glissante pourra être défini par un paramètre de taux de recouvrement exprimé par exemple en pixels. La valeur de ce paramètre de taux de recouvrement pourra être diminuée pour alléger les temps de calcul. Bien entendu, dans ce cas, la qualité des résultats ne devra pas se trouver 25 dégradée. Cette approche, conforme au procédé selon l'invention, pourra être mise en oeuvre dans le cadre de tout contexte opérationnel mettant en jeu un cube sismique ou des cubes sismiques d'une même zone géographique. Les pixels repérés par le processus de segmentation de la zone guide entourant l'horizon pointé et caractérisés par leurs indices EII et EIC pourront permettre l'élaboration de cartes d'attributs de qualification d'un horizon afin de permettre d'effectuer un contrôle qualité dudit horizon. Des modes d'exécution de l'invention seront décrits ci-après, à titre d'exemples non limitatifs, avec référence aux dessins annexés dans lesquels: La figure 1 est une représentation en perspective d'un cube sismique et d'un horizon à qualifier; Les figures 2a et 2b sont des sections verticales de la figure 1 10 représentant respectivement l'horizon à qualifier et l'horizon à qualifier entouré d'une zone guide; La figure 3 est une représentation d'une section de la zone guide comportant une fenêtre glissante; La figure 4 est une représentation du déplacement de la fenêtre de la figure 3; Les figures 5a, 5b et 5c sont des représentations d'une fenêtre 20 glissante segmentée respectivement en deux, trois et trente régions; La figure 6 est une représentation du déplacement de la fenêtre de la figure 3 faisant apparaître les frontières suivant leur hiérarchie; Les figures 7a, 7b et 7c sont des représentations de cartes d'attributs de qualification. Les figures 8a et 8b sont des représentations sur des sections verticales d'anomalies, respectivement des écarts par rapport à un extremum local et un saut de phase de l'horizon, mises en évidences grâce aux cartes d'attributs de qualification des figures 7a, 7b et 7c. Le premier exemple est illustré sur les figures 1 à 8b. Cet exemple concerne le contrôle qualité d'un pointé d' horizons (évènements cohérents correspondant à des interfaces géologiques) et en particulier d'un horizon pointé 1 sur un cube sismique 2 (Figures 1 et 2a) . Un horizon est très étendu dans l'espace et représente une interface géologique de sédimentation entre couches géologiques, il résulte généralement d'une interprétation par une géophysicien ou un géologue. La détermination des indices EH et ECI, conformément au procédé selon l'invention, contribue à définir les parties d'horizons jugées anomaliques qu'il convient d'examiner pour garantir une cohérence géologique des horizons pointés. La variable ou attribut sismique utilisé est l'amplitude du cube. Les premières étapes du procédé comprennent le choix: d'un critère de segmentation hiérarchique de la variable tel que l'amplitude, - d'une zone du cube sismique comprenant l'horizon pointé à qualifier c'est-à- dire sur lequel doit être effectué un contrôle de la qualité, - d'une fenêtre glissante 3 sur ladite zone et de ses caractéristiques (dimensions, taux de recouvrement...) (Figure 3), Le critère de segmentation hiérarchique choisi est l'amplitude correspondant aux valeurs de la variable elle-même, c'est-à-dire que l'on prend un critère dynamique lié aux valeurs de la variable. Ainsi, par exemple, plus la valeur moyenne de la variable associée à une frontière de segmentation sera importante plus les pixels associés à cette frontière se trouveront dans les premiers niveaux de hiérarchie. La zone du cube sismique à qualifier est limitée dans l'espace dans une zone guide 4 du cube sismique 2 comprenant et entourant l'horizon pointé à qualifier 1 (Figure 2b). Cette zone 4 sert de masque pour sélectionner les pixels du cube sismique sur lesquels s'effectueront les étapes suivantes du procédé selon l'invention. La largeur de la zone 4 est définie selon une section verticale (Figure 2b) du cube sismique par un paramètre vertical d'écart en pixels de part et d'autre de l'horizon pointé 1. Il est à noter que l'on pourra faire varier spatialement les dimensions de la zone guide entourant l'horizon. Les dimensions de la fenêtre glissante doivent être fixées. Pour cela, on pourra avantageusement réaliser au préalable une analyse géostatistique pour définir des critères objectifs d'aide au choix des dimensions de la fenêtre glissante. Ces critères peuvent comprendre les longueurs de corrélations spatiales par exemple déterminées grâce à l'identification des portées d'un variogramme. Le variogramme est une fonction statistique permettant d'analyser les corrélations spatiales à l'intérieur d'un jeu de données spatial, les portées observées sur le variogramme renseignant sur les dimensions moyennes des structures (évènements) dans une image. Si nécessaire, la valeur du paramètre du taux de recouvrement des différentes positions de la fenêtre glissante 3 pourra être diminuée pour alléger les temps de calcul. Bien entendu, cette valeur sera choisie de façon à ne pas dégrader la qualité 25 des résultats. Sur la figure 4, ledit taux de recouvrement est choisi égal à un pixel et un repère fixe R est placé. La position 6 de la fenêtre est décalée d'un pixel vers la droite par rapport à la position 5 de la fenêtre et la position 7 de la fenêtre est décalée d'un pixel vers la droite par rapport à la position 6 de la fenêtre. Les pixels repérés sont représentés par des points P. Le critère de segmentation choisi, les dimensions de la fenêtre glissante et sa valeur de taux de recouvrement fixés et la zone guide 4 définie, tous les pixels du cube sismique de la zone guide sont segmentés par la fenêtre glissante. Dans cet exemple, la phase de segmentation hiérarchique, pour des 5 positions de la fenêtre glissante se recouvrant,comprend la segmentation de l'image en trente régions et pour les pixels repérés au moins une fois par une frontière de segmentation: - la sauvegarde du nombre de fois où le pixel est repéré par une frontière de segmentation pour toutes les positions de la fenêtre glissante auxquelles il appartient et le calcul du rapport de ce nombre par le nombre théorique maximal de fois où ce même pixel pourrait être repéré, ce rapport correspondant à l'indice ECI (Figure 4) ; l'attribution pour chaque frontière de segmentation d'un niveau hiérarchique selon un barème B basé sur un calcul fonction d'un critère déterminé tel que la surface de la frontière, de manière à obtenir l'indice EII (Figure 6), les différents niveaux pouvant être matérialisés par des couleurs. Il est à noter que le nombre de régions de segmentation varie d'une fenêtre glissante à l'autre. En effet, par exemple, si on place une fenêtre glissante sur la France, la segmentation fera apparaître 22 régions, ce ne sera pas la même chose pour une fenêtre glissante tombant sur l'Allemagne. Plus particulièrement, lorsque l'image est segmentée en deux régions 8, 9, 25 une frontière 10 apparaît (Figure 5a). Cette frontière correspond au premier niveau de hiérarchie et sa valeur est fixée à un. Puis, l'image est segmentée en trois régions 8, 11, 12, une frontière 13 apparaît (Figure 5b). Cette frontière correspond au deuxième niveau de hiérarchie et sa valeur est fixée à deux. La segmentation se poursuit ainsi jusqu'à obtenir trente régions et les frontières correspondantes (Figures 5c et 6). Puis le procédé comporte les étapes suivantes: - une phase d'attribution aux pixels non repérés d'une valeur d'indices El et EIC correspondant à une non-valeur caractéristique numérique; - la détermination des évènements cohérents de la zone par seuillage des indices EII et ECI, l'indice EII devant être inférieur à une valeur de seuil et l'indice ECI devant être supérieur à une valeur de seuil, Tous les pixels repérés par le processus de segmentation de la zone guide entourant l'horizon pointé étant caractérisés par leurs indices EH et EIC, plusieurs cartes d'attributs de qualification de l'horizon pointé peuvent être calculées afin de permettre d'effectuer le contrôle de l'horizon. Dans cette zone, un pixel est défini par trois dimensions x, y et z. Une carte d'attributs de qualification est une représentation selon deux dimensions x et y. Un attribut défini par un couple de coordonnées x et y peut correspondre à une valeur moyenne par verticale (ligne axée selon l'axe vertical z de profondeur du cube sismique) pour l'ensemble des pixels repérés de la zone guide entourant l'horizon (par exemple ECI ou EII moyen par verticale) ou alors à la sélection d'un pixel sur la verticale, par exemple celui de plus fort ECI ou de plus faible EII, et à une valeur associée à ce pixel. La figure 7a représente une carte de l'horizon à qualifier. Un code couleur représente la profondeur Z de l'horizon. La figure 7b représente une carte de qualification de position) des pixels retenus parce qu'appartenant à l'évènement local le plus cohérent. La figure 7c est une carte de l'attribut de qualification du pointer de l'horizon 30 à qualifier, cet attribut étant nommé attribut de cohérence LPE. Pour chaque verticale de la zone guide définie par un couple de coordonnées x et y a été sélectionné le pixel appartenant à l'événement local le plus cohérent défini comme étant le pixel de plus faible EII et de plus fort ECI. La différence entre la position verticale de l'horizon pointé et la position 2884636 -12- verticale des pixels retenus permet de mettre en évidence des anomalies de position de l'horizon pointé. Ces cartes permettent de mettre en évidence des anomalies et de les repérer dans le plan XY. Ces anomalies peuvent être alors retrouvées, à partir de l'interprétation des anomalies sur les cartes 2D de qualification, sur le cube de données, en section par exemple comme illustré sur les figures 8a et 8b). Par ailleurs, il se trouve que dans le cas particulier de l'attribut de cohérence LPE, on peut également représenter cet attribut dans le cube. Ces anomalies peuvent consister, pour un pixel donné, à un écart par rapport à l'extremum local c'est-à-dire un écart entre l'évènement le plus cohérent c'est-à-dire le pixel repéré et l'horizon pointé qui est supérieur à l'écart maximum admis ou inférieur à un écart minimum admis. Sur l'exemple de la figure 8a, l'horizon pointé est représenté par une ligne noire 14 et les pixels repérés qualifiés, c'est-à-dire les pixels de plus fort ECI et de plus faible EII, sont représentés par des ronds blancs 15. Le cadre 16 contient une portion d'horizon pointé situé plus haut que les pixels repérés constituant ainsi. Ces anomalies peuvent également consister, en des sauts de phase de l'horizon correspondant à des erreurs de pointé qui représentent le passage non autorisé d'un horizon à un autre ou d'un évènement à un autre: Sur l'exemple de la figure 8b, l'horizon pointé est représenté par une ligne noire 17 et les pixels repérés qualifiés, c'est-à-dire les pixels de plus fort ECI et de plus faible EN sont représentés par une ligne blanche 18. Le cadre 19 contient un saut de phase de l'horizon pointé. Le deuxième exemple concerne la détermination d'évènements cohérents sur différents cubes sismiques d'une même zone géographique, ces cubes étant générés avec des paramètres de traitement différents. En effet, dans les solutions actuellement proposées, il arrive couramment en traitement et en interprétation sismique que deux ou plusieurs cubes sismiques soient définis sur une même zone géographique. Dépendant du contexte opérationnel, l'analyse des différences entre les deux cubes est alors généralement instructive. Elle peut mettre en évidence des différences liées: à la production d'hydrocarbures dans le temps (problématique réservoir 4D) ; - à des natures de variable étudiée différentes (problématique 4C par exemple) ; - à des paramètres de traitement différents. Or, la caractérisation de ces différences (localisation, quantification) n'est jamais aisée, la simple différence entre les deux cubes sismiques se révélant souvent insuffisante, voire stérile. En effet, des décalages géographiques mêmes faibles des évènements cohérents apparentés sur chacun des deux cubes sismiques diminuent notablement la valeur de l'information du cube différence. Ainsi, la simple différence entre les deux cubes sismiques ne permet pas de faire la distinction entre des différences liées d'un part aux décalages géographiques des évènements sismiques cohérents et d'autre part à des variations en énergie (en valeur) de ces évènements. Or, la quantification de l'impact des paramètres de traitement sur les images en amplitude de cubes sismiques est primordiale dans une étude pétrolière orientée réservoir, cet impact influant in fine sur le calcul des réserves en hydrocarbures. La détermination des indices EII et ECI sur chacun des deux cubes d'une même zone géographique, conformément au procédé selon l'invention, permet de s'affranchir des décalages géographiques spatiaux des évènements cohérents de même origine géologique et ainsi de permettre la mise en oeuvre d'outils quantitatifs complémentaires d'analyse des différences entre cubes sismiques d'une même zone géographique. - 14 - Le procédé comporte alors les étapes suivantes: - la détermination, pour chaque cube sismique, des indices Ell et ECI conformément au procédé selon l'invention développé dans l'exemple 1 en choisissant: o l'amplitude de la variable amplitude sismique comme critère de segmentation, o les dimensions de la fenêtre glissante de façon à ce que, plus on veut retenir d'évènements cohérents de petite taille, plus les dimensions de la fenêtre glissante sont restreintes, o une valeur de seuil importante pour effectuer un seuillage sur l'indice ECI (sélection des valeurs ECI hautes) et l'indice Ell (sélection des valeurs Ell hautes) afin de déterminer un nombre d'évènements cohérents important sur les deux cubes sismiques, une fois les évènements cohérents identifiés sur chacun des cubes sismiques, un appareillement de ces évènements consistant à définir qu'un événement géologique donné est repéré sur chacun des deux cubes et correspond à des évènements sismiques donnés sur les cubes; ainsi appareillés, ces événements servent de support pour le calcul d'attributs de caractérisation des différences dans les cubes sismiques. Le critère d'appareillement pourra par exemple être un critère de proximité spatiale: deux évènements suffisamment proches c'est à dire appartenant à un voisinage commun sont considérés comme correspondant au même évènement géologique. Selon une variante, on pourra se servir des résultats de segmentation de l'un des cubes sismiques pour réaliser sous contrainte la segmentation de l'autre cube. En particulier l'algorithme de Ligne de Partage des Eaux permet de construire des marqueurs pour la segmentation d'un cube sismique, ces marqueurs étant obtenus à partir des résultats de segmentation de l'autre cube. Selon une autre variante, on pourra diminuer la valeur du paramètre taux de recouvrement des fenêtres glissantes pour alléger les temps de calcul de telle façon que la qualité des résultats ne soit pas dégradée. L'appareillement des évènements cohérents sur chacune des deux images peut être réalisé par la recherche, avec une certaine tolérance verticale, des pixels définissant des évènements cohérents d'une des deux images sismiques les plus proches des pixels définissant des évènements cohérents de l'autre image sismique. La sélection des évènements cohérents appareillés permet de s'affranchir des effets de décalage géographique des événements cohérents et par conséquent de mieux caractériser les différences entre cubes sismiques par le calcul par exemple de la simple différence des valeurs, du décalage géographique des évènements cohérents sur chacun des deux cubes ou de tout autre attribut pertinent. Cette approche, conforme au procédé selon l'invention, pourra être mise en 20 oeuvre dans le cadre de tout contexte opérationnel mettant en jeu des cubes sismiques différents mais d'une même zone géographique (4D, 4C, multi-3D, etc.). -16-Revendications 1. Procédé de détermination d'évènements cohérents d'une image sismique, caractérisé en ce qu'il comprend les phases opératoires suivantes: une phase de choix: o d'un critère de segmentation fonction d'une variable représentée sous forme d'image sismique tel que l'amplitude, le contraste, le gradient d'une zone de l'image à qualifier, o d'au moins une fenêtre glissante sur ladite zone et de ses caractéristiques (dimensions, taux de recouvrement...), - une phase de segmentation hiérarchique, pour des positions de la fenêtre glissante se recouvrant, comportant une segmentation de ladite zone en n régions et pour chaque pixel repéré au moins une fois par une frontière de segmentation: o la sauvegarde du nombre de fois où le pixel est repéré par une frontière de segmentation pour chaque fenêtre glissante à laquelle il appartient et le calcul du rapport de ce nombre par le nombre théorique maximal de fois où ce même pixel pourrait être repéré, ce rapport correspondant à l'indice ECI; o l'attribution pour chaque frontière de segmentation d'un niveau hiérarchique basé sur un calcul fonction d'un critère déterminé tel que la surface de la frontière, de manière à obtenir l'indice EII; - une phase d'attribution aux pixels non repérés d'une valeur d'indices El et EIC correspondant à une non-valeur caractéristique numérique ou alphanumérique; - une phase de détermination des évènements cohérents de l'image par seuillage des indices EH et ECI c'est-à-dire en ne sélectionnant que 30 des pixels repérés de la zone correspondant à des valeurs inférieures ou supérieures à un seuil fixé. 2. Procédé selon la revendication 1, caractérisé en ce que la variable est un attribut sismique tel que l'amplitude, la réflectivité, l'impédance 3. Procédé selon la revendication 1, caractérisé en ce que les caractéristiques de la fenêtre glissante comprennent ses dimensions, le taux de recouvrement d'une position i de la fenêtre par rapport à une position i-r, le taux de recouvrement étant défini par 1/r, r étant le module du vecteur déplacement dans une des directions principales de l'image sismique multiple d'un pas de grille dans cette direction. 4. Procédé selon la revendication 1, caractérisé en ce que le nombre de n régions de segmentation de la zone est 15 dépendant de la zone segmentée. 5. Procédé selon la revendication 1, caractérisé en ce que la deuxième phase s'effectue par l'algorithme de segmentation hiérarchique de Ligne de Partage des Eaux . 6. Procédé selon la revendication 1, caractérisé en ce que la deuxième phase s'effectue sur une image sismique préalablement filtrée c'est-àdire image transformée par un filtre. 7. Procédé selon la revendication 6, caractérisé en ce que ledit filtre améliore la continuité de l'image et élimine du bruit aléatoire haute fréquence. 8. Procédé selon la revendication 1, caractérisé en ce que le recouvrement des fenêtres glissantes est défini par un paramètre de taux de recouvrement exprimé par exemple en pixels. 2884636 -18- 9. Utilisation du procédé selon la revendication 1 pour le contrôle qualité d'un pointé d'horizons sur un cube sismique. 10. Utilisation du procédé selon la revendication 1 pour la détermination d'évènements cohérents sur différents cubes sismiques d'une même zone géographique, ces cubes étant générés avec des paramètres de traitement différents.