Brevet FR2808336A1 — publié le 20011102

DESCRIPTION La présente invention concerne une méthode d'interprétation chronostratigraphique d'une section ou d'un bloc sismique, c'est-à-dire un habillage géologique d'une section ou d'un bloc sismique. La présente invention se situe dans le contexte de l'exploration pétrolière et permet un passage du domaine géophysique vers le domaine géologique. Etat de l'art La méthode selon l'invention, s'applique sur des sections sismiques ou des blocs sismiques. Une section sismique est formée par la juxtaposition dans un plan de signaux monodimensionnels, échantillonnés et appelés traces sismiques. De la même façon, un bloc sismique est formé par la juxtaposition de traces sismiques dans un volume. On appelle "coupe sismique" soit une section sismique, soit une tranche de bloc sismique. Une coupe sismique offre une vue de la juxtaposition des traces sismiques contenues dans le plan de coupe. Ces vues sont des images sismiques, qui seront appelées coupe-images sismiques dans l'exposé de la mise en oeuvre de la méthode. Dans une image sismique, I'intensité lumineuse d'un pixel est proportionnelle à la grandeur sismique représentée par les signaux monodimensionnels. L'interprétation chrono- stratigraphique de sections sismiques ou de blocs sismiques passe par la synthèse d'horizons sismiques dans la section ou le bloc. Plusieurs méthodes ont été élaborées pour assurer des synthèses d'horizons. Leurs résultats sont plus ou moins bons; ils dépendent en fait de l'environnement géologique dont la coupe sismique offre l'image. Ainsi, dans les régions o les couches géologiques sont à dominante monoclinale, la synthèse des horizons par mesure de similarité entre traces voisines donne de bons résultats. En revanche, dans des zones o la géologie est plus perturbée, il est préférable de calculer d'abord les vecteurs gradients de l'intensité lumineuse entre pixels voisins puis de mettre en oeuvre une synthèse d'horizons par intégration du champ d'orientation des vecteurs gradients calculés. La thèse de Marc Donias, soutenue le 28/01/99 à l'Université Bordeaux I et intitulée "Caractérisation de champs d'orientation par analyse en composantes principales et 2 2808336 estimation de la courbure. Application aux images sismiques", décrit en détail les modes de réalisation de synthèse d'horizons évoqués ci-dessus. La synthèse d'horizons appliquée à chaque pixel de la coupe-image sismique crée autant d'horizons qu'il y a de pixels dans l'image. Les horizons sismiques, en tant que marqueurs de la géologie locale, ne peuvent pas se croiser. En revanche, ils peuvent converger, fusionner localement et se confondre, voire diverger. La fusion d'horizons conduit à la notion d'accumulation de synthèses. Pour réaliser une accumulation des synthèses d'horizons, on définit une matrice identique en taille à l'image sismique. Chaque élément de la matrice est associé à un pixel de l'image et est initialement affecté d'une valeur nulle. Pour chaque pixel de l'image, on calcule une courbe de continuité qui correspond à la synthèse de l'horizon passant par ledit pixel. Toutes les courbes de continuité sont transverses à la dimension verticale de l'image. Lors du calcul des courbes de continuité, on incrémente d'une unité un élément de la matrice chaque fois que le pixel auquel il est associé dans l'image est traversé par une courbe de continuité. Le calcul d'une courbe de continuité transverse à la dimension verticale de la coupe- image en un pixel donné consiste à calculer les gradients d'intensité lumineuse pour tous les pixels inclus dans un voisinage du pixel choisi, puis à calculer un gradient local à partir des mesures du gradient obtenues sur le voisinage et à affecter le gradient local au pixel choisi. La courbe de continuité est ensuite développée en cheminant transversalement de pixel en pixel à partir du pixel choisi jusqu'aux limites latérales verticales de l'image, dans les deux directions indiquées par le gradient local et son opposé, en itérant les deux étapes précédentes. Lorsque tous les pixels de la coupe-image sismique ont été balayés, la matrice réalisant les accumulations de synthèses est représentée sous la forme d'une nouvelle image dans laquelle chaque pixel possède une intensité lumineuse proportionnelle au nombre cumulé dans l'élément de la matrice correspondant, lequel est au moins égal à un. Les limites observées sur cette image donnent une bonne idée de l'organisation des couches géologiques dans le sous-sol. Apport de l'invention La méthode selon l'invention, exploite l'image obtenue par accumulation des synthèses pour déterminer les dépôts géologiques tels qu'ils se sont déposés et non 3 2808336 pas tels qu'ils sont observés aujourd'hui sous la forme de couches. Pour ce faire, la méthode définit une transformation de l'échelle verticale de la coupe sismique mesurée en temps sismiques en une échelle verticale géologique mesurée en temps géologiques. La méthode permet ainsi de définir les vitesses de sédimentations qui ont présidé aux dépôts des couches géologiques. En particulier, elle fait apparaître les hiatus géologiques, c'est-à-dire les érosions et les lacunes. Définition de l'invention La présente invention a pour objet une méthode d'interprétation chrono- stratigraphique d'une coupe-image sismique S, comportant une dimension horizontale ou largeur et une dimension verticale ou hauteur en direction du sous-sol et constituée de colonnes de pixels, qui consiste à: - définir une matrice M identique en taille à la coupe-image S et constituée d'éléments dont chacun est associé à un pixel de la coupe- image S et est initialement affecté d'une valeur nulle, - pour chaque pixel i de la coupe-image S, calculer une courbe de continuité Ci passant par ledit pixel et transverse à la dimension verticale de la coupe-image S, incrémenter d'une unité un élément de la matrice M chaque fois que le pixel auquel il est associé dans la coupe-image S est traversé par une courbe Ci, ladite méthode se caractérisant en ce qu'elle consiste en outre à: - construire, pour chaque colonne c de la matrice M, un histogramme Hc constitué d'un nombre de classes égal au nombre d'éléments de ladite colonne c, chaque classe correspondant à l'un des éléments de la colonne c et contenant un nombre d'échantillons égal à la valeur cumulée dans l'élément considéré de la matrice M, laquelle valeur est égale au nombre de courbes passant par le pixel associé audit élément, le nombre total d'échantillons distribués dans l'histogramme construit pour chaque colonne étant égal au nombre total de pixels dans la coupe-image S, - égaliser chaque histogramme Hc pour produire un histogramme égalisé Hc', - définir une coupe-image S' vide, de largeur en pixels identique à la largeur en pixels de la coupe-image S et de hauteur en pixels égale au nombre de classes des histogrammes Hc', - affecter à chaque colonne c' de la coupe-image S' la distribution définie par l'histogramme égalisé Hc' en attribuant à chaque pixel de la colonne c' la cardinalité du contenu de la classe associée de Hc', - délimiter dans la coupe-image S' les groupes de pixels contigus contenant des échantillons et étiqueter chacun desdits groupes, - attribuer à chaque pixel de la coupe-image S l'étiquette donnée au groupe de pixels auquel il a été affecté dans la coupe-image S' et visualiser la coupe-image S étiquetée. Selon une autre caractéristique, le calcul de la courbe de continuité Ci transverse à la dimension verticale de la coupe-image S en un pixel i donné consiste à: - calculer les gradients d'intensité lumineuse pour tous les pixels inclus dans un voisinage Vi du pixel i, - calculer un gradient local Gi à partir des mesures du gradient obtenues sur le voisinage Vi et affecter le gradient Gi au pixel i, - cheminer transversalement de pixel en pixel à partir du pixel i jusqu'aux limites latérales verticales de la coupe- image S dans les deux directions indiquées par le gradient Gi et son opposé -Gi, en itérant les deux étapes précédentes. La présente invention peut aussi s'appliquer sur des bloc-images sismiques en trois dimensions. La présente invention a donc également pour objet une méthode d'interprétation chrono-stratigraphique d'un bloc- image sismique B comportant deux dimensions horizontales, à savoir largeur et profondeur, et une dimension verticale ou hauteur en direction du sous-sol et constitué de colonnes de pixels, qui consiste à - définir un bloc N identique en taille au bloc-image B et constitué d'éléments dont chacun est associé à un pixel du bloc-image B et est initialement affecté d'une valeur nulle, - pour chaque pixel i du bloc-image B, calculer une surface de continuité Si passant par ledit pixel et transverse à la dimension verticale du bloc-image B, - incrémenter d'une unité un élément du bloc N chaque fois que le pixel auquel il est associé dans le bloc-image B est traversé par une surface Si, ladite méthode se caractérisant en ce qu'elle consiste en outre à: - construire, pour chaque colonne c du bloc N, un histogramme Hc constitué d'un nombre de classes égal au nombre d'éléments de ladite colonne c, chaque classe correspondant à l'un des éléments de la colonne c et contenant un nombre d'échantillons égal à la valeur cumulée dans l'élément considéré du bloc N, laquelle valeur est égale au nombre de surfaces passant par le pixel associé audit élément, le nombre total d'échantillons distribués dans l'histogramme construit pour chaque colonne étant égal au nombre total de pixels dans le bloc-image B, - égaliser chaque histogramme Hc pour produire un histogramme égalisé Hc', - définir un bloc-image B' vide, de largeur en pixels et de profondeur en pixels identiques à la largeur en pixels et à la profondeur en pixels du bloc-image B, et de hauteur en pixels égale au nombre de classes des histogrammes Hc', - affecter à chaque colonne c' du bloc-image B' la distribution définie par l'histogramme égalisé Hc' en attribuant à chaque élément de la colonne c' le contenu de la classe associée de Hc', délimiter dans le bloc-image B' les groupes de pixels contigus contenant des échantillons et étiqueter chacun desdits groupes, - attribuer à chaque pixel du bloc-image B l'étiquette donnée au groupe de pixels auquel il a été affecté dans le bloc-image B' et visualiser le bloc-image B étiqueté. Selon une autre caractéristique, le calcul de la surface de continuité Si transverse à la dimension verticale du bloc-image B en un pixel i donné consiste à: - calculer les gradients d'intensité lumineuse pour tous les pixels inclus dans un voisinage de chaque pixel i d'une colonne de pixels Ki, - pour chaque pixel i de la colonne Ki, effectuer une analyse en composantes principales sur les gradients calculés dans le voisinage Vi du pixel i afin de déterminer un couple de vecteurs directeurs du plan tangent à la surface Si au pixel i de la colonne Ki, - cheminer concentriquement de colonne en colonne à partir de la colonne Ki jusqu'aux limites latérales verticales du bloc-image B, en itérant les deux étapes précédentes. Figures Les figures la et lb représentent, à des échelles différentes, la même coupe-image sismique avant l'application de la méthode selon l'invention. La figure 2a représente une courbe de continuité associée à un pixel de la coupe- image sismique et la figure 2b représente une vue de toutes les courbes de continuité. La figure 3a représente un histogramme associé à une verticale de la figure 2b et la figure 3b montre le même histogramme égalisé. La figure 4 représente une coupe-image égalisée, appelée coupe-image chrono- stratigraphique. La figure 5a représente une coupe-image égalisée binarisée et la figure 5b montre la même coupe-image étiquetée. La figure 6 représente une coupe-image sismique étiquetée après application de la méthode selon l'invention. Description complète La méthode selon l'invention est une méthode d'interprétation chronostratigraphique automatique d'une coupe-image sismique. En se référant aux figures, un mode de mise en oeuvre de cette méthode est donné ci-après. La figure la représente une coupe-image sismique S. La coupe-image S est qualifiée de sismique car elle représente une vision du sous-sol issue d'une campagne d'exploration sismique. La coupe-image S comporte deux dimensions; elle est définie par une extension horizontale selon un axe horizontal et par une extension verticale selon un axe vertical en direction du sous-sol. La coupe-image S est composée de pixels régulièrement répartis selon un pas horizontal sur l'axe horizontal et un pas vertical sur l'axe vertical. La coupe-image S contient en particulier un nombre de colonnes de pixels égal au quotient de l'extension horizontale par le pas horizontal et un nombre de pixels par colonne égal au quotient de l'extension verticale par le pas vertical. En particulier, le trait noir vertical 10 sur la figure la représente une colonne de pixels qui est le support des figures suivantes dans la description de la méthode. 7 2808336 Pour une mise en oeuvre de la méthode selon l'invention, on définit une matrice M identique en taille à la coupe-image S. Le nombre de lignes de la matrice M est égal au nombre de pixels dans une colonne de la coupeimage S et le nombre de colonnes de la matrice M, égal au nombre de colonnes de la coupe-image S, est égal au nombre de pixels dans une ligne de cette coupe-image. La matrice M est ainsi constituée d'autant d'éléments qu'il y a de pixels dans la coupe-image S et chaque élément est associé à un pixel de la coupe-image S. Tous les éléments de la matrice M sont des entiers ayant initialement une valeur nulle. Pour chaque pixel i de la coupe-image S, on calcule une courbe de continuité Ci passant par ledit pixel et transverse à la dimension verticale de la coupe-image S. Le calcul de cette courbe Ci passe par le calcul d'un gradient local Gi de l'intensité lumineuse au pixel i. On définit l'intensité lumineuse associée à un pixel comme la représentation dans une palette, par exemple une palette de niveaux de gris, d'un attribut sismique, par exemple l'amplitude du signal sismique; dans les conditions prises dans l'exemple, une amplitude forte serait manifestée par un pixel clair et une amplitude faible par un pixel sombre. On peut alors calculer un gradient d'intensité lumineuse Gi entre un pixel i et ses voisins, lequel gradient Gi est en fait le gradient de l'attribut sismique considéré. Le gradient Gi est affecté au pixel i. Le gradient Gi comporte une composante horizontale et une composante verticale et il est déterminé notamment par une analyse en composantes principales appliquée sur tous les gradients calculés sur tous les pixels inclus dans un voisinage Vi du pixel i. Le voisinage Vi est défini par une fenêtre centrée sur le pixel i, de préférence une fenêtre d'une taille égale à 7 pixels par 7 pixels. La direction d'allongement principale du nuage des mesures de gradients dans le voisinage Vi, donnée par la direction du premier axe d'inertie, permet de déterminer une direction locale du gradient Gi, laquelle est attribuée au pixel i. La figure 2a montre une courbe de continuité Ci associée à un pixel i. La courbe Ci est qualifiée de courbe de continuité car elle s'étend dans l'image à partir du pixel i vers des pixels présentant des caractéristiques semblables à celles du pixel i. La courbe Ci est développée transversalement à la dimension verticale de la coupe- image S, laquelle est montrée sur la figure lb à la même échelle que la figure 2a et sur laquelle on peut déceler l'horizon synthétisé par la courbe Ci. La courbe Ci est 8 2808336 obtenue en cheminant de pixel en pixel dans les directions successives déterminées par les gradients successifs Gi calculés en continu et leurs opposés -Gi. En conséquence, la valeur de la dérivée de la courbe de continuité Ci en chacun des point-pixels qui la constituent est la valeur des gradients locaux calculés en ces mêmes pixels. Chaque fois qu'une courbe de continuité Ci passe par un pixel i, on incrémente d'une unité l'élément correspondant au pixel i de la matrice M. Les éléments de la matrice M sont en quelque sorte des compteurs associés aux pixels de la coupe-image S. La figure 2b montre une vue de toutes les courbes de continuité transverses. Plus le tracé des courbes transverses apparaît gras, plus le nombre de courbes transverses superposées est élevé. Plus concrètement, si la coupe-image S matérialise, par exemple, la valeur de l'amplitude du signal sismique, la courbe Ci, démarrée au pixel i, recherchera une continuité d'amplitude au travers de la coupe-image S sur la base d'une similarité avec l'amplitude observée au pixel i. Si la valeur de l'amplitude au pixel i est très forte, une courbe développée depuis le pixel i dans l'image tendra à suive un marqueur géologique. Une fois toutes les courbes de continuité calculées, les éléments de la matrice M contiennent les résultats des comptages des courbes de continuité passant par chaque pixel de la coupe-image S, selon la stricte correspondance élément à pixel. C'est précisément le contenu de la matrice M qui est visualisé sur la figure 2b. Pour chaque colonne de la matrice M, colonne qui correspond à une colonne de pixels prise dans l'image représentée sur la figure 2b, on construit un histogramme Hc comme celui qui est représenté sur la figure 3a et qui est associé à la colonne c indiquée par le trait noir vertical sur la figure 2b. Cet histogramme est constitué d'un nombre de classes égal au nombre de pixels dans ladite colonne c. Chaque classe de l'histogramme est associée à un élément dans la colonne de la matrice M, lequel élément correspond à un pixel dans la colonne de pixels correspondante de la coupe-image S; plus le nombre de courbes de continuité passant par ce pixel est élevé, plus la classe correspondante de Hc est remplie. De fait, I'histogramme Hc de la figure 3a présente la distribution des courbes de continuité sur la verticale 10 tracée sur la coupe-image S, distribution qui est précisément stockée dans la colonne correspondante de la matrice M associée à la coupe-image S. Le nombre d'échantillons contenus dans les classes des différents histogrammes, associés aux différentes colonnes, varie entre un minimum et un maximum, tous deux étant des entiers positifs. Le minimum est un car il existe toujours au moins une courbe de continuité passant par chaque pixel. Le maximum est variable mais il est borné par le nombre total de courbes de continuité, qui n'est autre que le nombre total de pixels de la coupe-image S. A chaque histogramme Hc, on applique alors un algorithme d'égalisation pour produire un histogramme égalisé correspondant Hc'. Pour ce faire, on pourra par exemple utiliser les algorithmes décrits dans l'ouvrage de MM R. Gonzalez et R. Woods intitulé "Digital Image Processing" et édité en 1992 par Addison-Wesley, ou bien dans l'ouvrage de MM J-P Coquerez et S. Philipp intitulé "Analyse d'images: filtrage et segmentation" et publié en 1995 chez Masson. L'application de l'algorithme d'égalisation d'histogramme sur l'histogramme de la figure 3a redistribue la population classée selon une nouvelle distribution représentée par l'histogramme Hc' de la figure 3b. L'algorithme d'égalisation d'histogramme, dont le but premier est d'accroître le contraste d'une image, a ici pour but de redistribuer la population classée en l'éclatant dans un plus grand nombre de classes et de façon plus équirépartie. Le nombre de classes de l'histogramme égalisé Hc' est égal au nombre de classes de l'histogramme Hc multiplié par un facteur de dilatation fonction de la résolution souhaitée. Dans la mise en oeuvre de la méthode selon l'invention, on définit une nouvelle coupe-image S' dont la largeur en pixels est identique à celle de la coupe-image S et dont la hauteur en pixels est égale au nombre de classes de l'histogramme Hc'. Initialement, la coupe-image S' est vide, c'est-à-dire que tous ses pixels apparaissent noirs. A partir de tous les histogrammes égalisés, on construit une nouvelle image dans la coupe-image S' en affectant à chaque colonne de la coupeimage S' la distribution définie par l'histogramme égalisé correspondant. Pour ce faire, on utilise les échantillons nouvellement répartis dans les classes des histogrammes égalisés pour calculer une intensité lumineuse pour chaque pixel associé à chaque classe. L'intensité lumineuse d'un pixel sur la coupe-image S', dans une palette de niveaux de gris par exemple, est fonction du nombre d'échantillons contenus dans la classe associée au pixel. Ainsi l'intensité lumineuse d'un pixel sera nulle et le pixel 2808336 apparaîtra noir sur la coupe-image S', s'il n'y a pas d'échantillon dans la classe associée au pixel. De même, I'intensité lumineuse sera non nulle et le pixel apparaîtra gris, plus ou moins clair, quand la classe associée contiendra un plus ou moins grand nombre d'échantillons. Il est important de noter que la population classée est la même pour tous les histogrammes, qu'ils soient égalisés ou non. Chaque histogramme en présente une distribution différente mais la somme du contenu de toutes les classes est identique pour tous les histogrammes. C'est précisément cet invariant qui permet de considérer la coupe-image S' comme une coupe-image chrono-stratigraphique, c'est- à-dire que les limites qu'elle fait apparaître sont considérées comme des isochrones au sens géologique. La figure 4 montre la coupe-image chrono-stratigraphique S' construite à partir des histogrammes égalisés. La transformation effectuée ici permet de passer d'une échelle verticale géophysique en profondeur graduée en mètres ou en milli-secondes à une échelle de temps géologique graduée en million d'années. En particulier, on notera le trait vertical 11 qui situe sur la figure 4 la colonne de pixels initialement choisie sur la figure la. On binarise alors la coupe-image chrono-stratigraphique S' en attribuant la valeur un en tout pixel contenant une valeur non-nulle et en conservant la valeur zéro en tout pixel sans information. Sur l'image binarisée, représentée par la figure 5a, apparaissent des zones blanches sur un fond noir. Une zone blanche indique que les pixels de la zone contiennent au moins une information sismique alors que le fond noir indique une absence d'information, c'est-à-dire qu'il n'y a pas d'échantillon dans les classes associées des histogrammes égalisés. D'un point de vue géologique, le fond noir révèle des périodes de lacune ou des périodes d'érosion, car il y a absence de dépôt de sédiments durant les périodes indiquées par tout segment noir observé sur une verticale de la coupe-image chrono-stratigraphique S'. On étiquète l'image binarisée en attribuant des motifs ou des couleurs tous différents aux zones blanches connexes de sorte que deux zones blanches distinctes ne possèdent jamais le même motif ou la même couleur. La figure 5b montre le résultat d'un étiquetage de la coupe-image binarisée 5a. il 2808336 On applique enfin à chaque colonne de la coupe-image 5b la transformation inverse de l'égalisation d'histogramme, c'est-à-dire que l'on remet les échantillons classés par l'histogramme égalisé dans leur classe initiale de l'histogramme originel. On reconstruit alors la section originelle en attribuant à chaque pixel la couleur que l'étiquetage lui a attribuée. Le résultat de la transformation inverse est présenté sur la figure 6. Les zones étiquetées et numérotées de 12 à 15 sur la figure 5b correspondent aux zones semblablement numérotées sur la figure 6. En particulier, sur la figure 5b, on note que la zone 12 se subdivise presque en deux zones 12 et 12' mais comme ces deux zones sont reliées par un petit isthme w, elles héritent du même étiquetage. Sur la figure 6, on peut aussi passer de 12 à 12' puisque leur séparation est discontinue du fait de l'existence de l'isthme w. La figure 6 est très proche d'une véritable coupe géologique. Les zones de même motif sont des zones de dépôt de sédiments relativement régulier et homogène. Deux zones de motifs différents en contact tendent à indiquer un changement de formation géologique, c'est-à-dire un nouveau régime de dépôt comme, par exemple, celui créé par le passage d'un phénomène transgressif à un phénomène régressif. En outre, la reconnaissance de marqueurs géologiques éventuels contenus dans la coupe-image étiquetée présentée sur la figure 6 permet de dater de façon cohérente dans l'espace l'ensemble des événements géologiques observés. L'ensemble de la méthode, selon l'invention, est aussi applicable en trois dimensions dans un bloc-image sismique. Un bloc-image sismique est défini par deux axes horizontaux et un axe vertical en direction du sous- sol. Dans cette mise en oeuvre, des surfaces de continuité transverses à la dimension verticale du bloc-image sont les homologues des courbes de continuité transverses à la dimension verticale de la coupe-image et seule l'étape de calcul de ces surfaces se différencie dans l'application de la méthode au bloc-image. Pour déterminer les surfaces de continuité du bloc-image passant par tous les pixels du bloc-image, on balaye successivement toutes les colonnes du bloc-image. Pour une colonne Ki choisie du bloc-image, on détermine simultanément tous les plans tangents aux surfaces de continuité en chacun des pixels formant la colonne Ki. Pour ce faire, on calcule d'abord un gradient local d'intensité lumineuse Gi en chaque pixel de la colonne Ki, dans un voisinage Vi du pixel considéré; le voisinage Vi est ici un petit cube centré sur le pixel considéré, par exemple un cube de 7x7x7 pixels. Les 12 2808336 gradients sont des vecteurs à trois composantes selon les trois dimensions du bloc- image. Ils sont déterminés en chaque pixel i de la colonne Ki par une analyse en composantes principales appliquée sur tous les gradients calculés sur tous les pixels inclus dans le voisinage Vi du pixel i. L'analyse en composantes principales en trois dimensions sur les gradients calculés dans le voisinage Vi d'un pixel i conduit à la détermination de trois vecteurs orthogonaux, à savoir un vecteur principal Fil associé au premier axe d'inertie définissant la ligne de plus grande pente du plan tangent au pixel i, un vecteur principal Fi2 associé au deuxième axe d'inertie orthogonal au précédent et directeur du plan tangent au pixel i et un vecteur principal Fi3 associé au troisième axe d'inertie orthogonal au plan tangent au pixel i. On obtient ainsi, en chaque pixel i de la colonne Ki, un couple de vecteurs directeurs (Fil, Fi2) du plan tangent à la surface de continuité recherchée et passant par ce pixel. On itère le processus décrit ci-dessus en propageant les analyses en composantes principales concentriquement de colonne en colonne à partir de la colonne Ki jusqu'aux limites latérales verticales du bloc-image. 13 2808336 REVENDICATIONS 1. Méthode d'interprétation chrono-stratigraphique d'une coupe-image sismique S qui comporte une dimension horizontale ou largeur et une dimension verticale ou hauteur en direction du sous-sol et qui est constituée de colonnes de pixels, consistant à: - définir une matrice M identique en taille à la coupe-image S et constituée d'éléments dont chacun est associé à un pixel de la coupe-image S et est initialement affecté d'une valeur nulle, - pour chaque pixel i de la coupe-image S, calculer une courbe de continuité Ci passant par ledit pixel et transverse à la dimension verticale de la coupe-image S, - incrémenter d'une unité un élément de la matrice M chaque fois que le pixel auquel il est associé dans la coupe-image S est traversé par une courbe Ci, et caractérisée en ce qu'elle consiste en outre à: - construire, pour chaque colonne c de la matrice M, un histogramme Hc constitué d'un nombre de classes égal au nombre d'éléments de ladite colonne c, chaque classe correspondant à l'un des éléments de la colonne c et contenant un nombre d'échantillons égal à la valeur cumulée dans l'élément considéré de la matrice M, laquelle valeur est égale au nombre de courbes passant par le pixel associé audit élément, le nombre total d'échantillons distribués dans l'histogramme construit pour chaque colonne étant égal au nombre total de pixels dans la coupe-image S, - égaliser chaque histogramme Hc pour produire un histogramme égalisé Hc', - définir une coupe-image S' vide, de largeur en pixels identique à la largeur en pixels de la coupe- image S et de hauteur en pixels égale au nombre de classes des histogrammes Hc', - affecter à chaque colonne c' de la coupe-image S' la distribution définie par l'histogramme égalisé Hc' en attribuant à chaque pixel de la colonne c' la cardinalité du contenu de la classe associée de Hc', - délimiter dans la coupe-image S' les groupes de pixels contigus contenant des échantillons et étiqueter chacun desdits groupes, 14 2808336 - attribuer à chaque pixel de la coupe-image S l'étiquette donnée au groupe de pixels auquel il a été affecté dans la coupe-image S' et visualiser la coupe-image S étiquetée. 2. Méthode selon la revendication 1, caractérisée en ce que le calcul de la courbe de continuité Ci transverse à la dimension verticale de la coupe-image S en un pixel i donné consiste à: - calculer les gradients d'intensité lumineuse pour tous les pixels inclus dans un voisinage Vi du pixel i, - calculer un gradient local Gi à partir des mesures du gradient obtenues sur le voisinage Vi et affecter le gradient Gi au pixel i, cheminer transversalement de pixel en pixel à partir du pixel i jusqu'aux limites latérales verticales de la coupe-image S dans les deux directions indiquées par le gradient Gi et son opposé -Gi, en itérant les deux étapes précédentes. 3. Méthode d'interprétation chrono-stratigraphique d'un bloc-image sismique B comportant deux dimensions horizontales, à savoir largeur et profondeur, et une dimension verticale ou hauteur en direction du soussol et constitué de colonnes de pixels, consistant à: - définir un bloc N identique en taille au bloc-image B et constitué d'éléments dont chacun est associé à un pixel du bloc-image B et est initialement affecté d'une valeur nulle, - pour chaque pixel i du bloc-image B, calculer une surface de continuité Si passant par ledit pixel et transverse à la dimension verticale du bloc-image B, - incrémenter d'une unité un élément du bloc N chaque fois que le pixel auquel il est associé dans le bloc-image B est traversé par une surface Si, et caractérisée en ce qu'elle consiste en outre a: - construire, pour chaque colonne c du bloc N, un histogramme Hc constitué d'un nombre de classes égal au nombre d'éléments de ladite colonne c, chaque classe correspondant à l'un des éléments de la colonne c et contenant un nombre d'échantillons égal à la valeur cumulée dans l'élément considéré du bloc N, laquelle valeur est égale au nombre de surfaces passant par le pixel associé audit élément, le nombre total d'échantillons distribués dans l'histogramme 2808336 construit pour chaque colonne étant égal au nombre total de pixels dans le bloc- image B, - égaliser chaque histogramme Hc pour produire un histogramme égalisé Hc', - définir un bloc-image B' vide, de largeur en pixels et de profondeur en pixels identiques à la largeur en pixels et à la profondeur en pixels du bloc-image B, et de hauteur en pixels égale au nombre de classes des histogrammes Hc', - affecter à chaque colonne c' du blocimage B' la distribution définie par l'histogramme égalisé Hc' en attribuant à chaque élément de la colonne c' le contenu de la classe associée de Hc', - délimiter dans le bloc-image B' les groupes de pixels contigus contenant des échantillons et étiqueter chacun desdits groupes, attribuer à chaque pixel du bloc-image B l'étiquette donnée au groupe de pixels auquel il a été affecté dans le bloc-image B' et visualiser le bloc-image B étiqueté. 4. Méthode selon la revendication 3, caractérisée en ce que le calcul de la surface de continuité Si transverse à la dimension verticale du bloc-image B en un pixel i donné consiste à: calculer les gradients d'intensité lumineuse pour tous les pixels inclus dans un voisinage de chaque pixel i d'une colonne de pixels Ki, - pour chaque pixel i de la colonne Ki, effectuer une analyse en composantes principales sur les gradients calculés dans le voisinage Vi du pixel i afin de déterminer un couple de vecteurs directeurs du plan tangent à la surface Si au pixel i de la colonne Ki, - cheminer concentriquement de colonne en colonne à partir de la colonne Ki jusqu'aux limites latérales verticales du bloc-image B, en itérant les deux étapes précédentes.