FR2646520A1 - Procede de pointe automatique et d'aide a l'interpretation, en particulier de section sismique, par des techniques d'analyse d'images - Google Patents

Abstract

Procédé d'analyse d'image composée de linéations appliqué en particulier à l'étude structurale et stratigraphique d'une section sismique basé sur un pointé par " contouring " des linéations et une décomposition des lignes obtenues après le pointé en segments suffisamment courts pour être de caractéristiques homogènes. Des paramètres caractérisant chacun des segments sont calculés. Des alignements de caractéristiques homogènes vis-à-vis d'un certain paramètre sont ensuite reconstitués à partir de la juxtaposition de segments élémentaires dont le paramètre associé satisfait un même critère.

Description

La présente invention concerne un procédé de pointé particulièrement

adapté aux images o une direction est privilégiée, comme par exemple l'image d'une section sismique. Appliqué à la sismique, le procédé selon l'invention permet notamment d'effectuer une sélection automatique des portions d'horizons sismiques de la section sismique présentant certaines caractéristiques et d'aider à l'interprétation structurale et stratigraphique de cette section. De façon très générale, on peut considérer que le sous-sol est constitué d'un empilement de couches géologiques de caractéristiques différentes qui s'organisent

dans l'espace suivant une certaine géométrie.

Pour définir la géométrie du sous-sol, les spécialistes de l'exploration du sous-sol, en particulier les géophysiciens en exploration pétrolière, utilisent une technique particulière appelée "sismique réflexion" qui consiste à émettre à la surface du sol des signaux acoustiques, et à les enregistrer après que lesdits signaux se soient propagés dans le sous-sol et se soient réfléchis sur les frontières acoustiques particulières que constituent les limites entre couches géologiques différentes, et que

l'on appelle réflecteurs.

Les signaux enregistrés sont traités pour obtenir des images particulières du sous-sol appelées- "sections sismiques", que l'on peut assimiler à des plans de coupe verticale du sous-sol, et sur lesquelles les réflecteurs apparaîtront comme des linéations superposées les unes aux autres, plus ou moins marquées, plus ou moins continues, subhorizontales ou déformées, à l'image des couches du sous-sol, et appelées par les spécialistes "horizons sismiques". Une section sismique est constituée par une succession, de gauche à droite ou de droite à gauche, d'enregistrements verticaux appelés "traces sismiques" constitués de signaux de plus ou moins forte amplitude, composés d'arches positives et négatives, chacun de ces signaux étant l'image d'un point d'un réflecteur. Les signaux images d'un même réflecteur corrèlent d'une trace à

l'autre faisant apparaître à l'oeil des lingations.

Certaines de ces linéations peuvent être très courtes et avoir néanmoins une signification: c'est le cas par exemple des "doublets". On dit qu'on a présence de doublets sur une section sismique quand entre deux linéations majeures bien marquées, associées à des arches d'amplitude de même signe et séparées sur une grande longueur par une seule arche de signal de signe contraire, vient s'intercaler sur quelques traces une arche de signal de fréquence plus élevée et de même signe que les arches

composant les linéations majeures.

Les sections sismiques élaborées par les géophysiciens sont ensuite transmises à des spécialistes -géophysiciens ou géologues- qui les interpréteront pour y

déceler par exemple les éventuels réservoirs.

L'interprétation consiste à pointer sur plusieurs sections sismiques décrivant une même région géographique et pouvant se recouper les unes les autres, un certain nombre

d'horizons sismiques déterminés, et à les cartographier.

L'opération de pointage s'effectue classiquement à la main.

Elle consiste à relier par corrélation optique entre eux, à l'aide d'un crayon de couleur, les signaux que l'on

considère comme appartenant à un même horizon sismique.

Usuellement, on pointera différents horizons avec différentes couleurs pour mieux les différencier. Les signaux images d'un même réflecteur présenteront des caractéristiques -forme, contenu fréquentiel, nombre

d'arches, amplitude, etc...- similaires en général.

Un tel pointage représente, on s'en doute, un travail long et fastidieux, donc coûteux. Il n'est pas rare en effet qu'une interprétation occupe plusieurs interprétateurs durant plusieurs semaines, quand un grand nombre d'horizons doit être pointé sur un grand nombre de sections. C'est le cas par exemple quand on effectue une

exploration de reconnaissance sur une très large zone-

géographique ou quand on effectue une "sismique de gisement"

sur zone restreinte mais avec des sections très rapprochées.

Ce travail manuel peut de plus et conséquemment conduire à des erreurs de pointé qui seront d'autant plus faciles à commettre que les horizons à pointer seront mal individualisés et que la géologie de la zone étudiée est agitée. L'existence de faille notamment peut interrompre un réflecteur et embarrasser l'interprétateur qui pourra corréler, de part et d'autre de la faille, deux linéations

de signaux qui ne seront pas associées au même réflecteur.

L'interprétateur de documents sismiques est donc confronté à un important volume de données figées sur une représentation dont il doit extraire visuellement des informations. Une grande partie de ce travail est longue, fastidieuse et reste très subjective, dépendant de la

qualité de la représentation comme de l'interprétateur.

La vision par ordinateurs a réalisé des' progrès tels, ces dernières années, qu'il est possible maintenant d'utiliser l'ordinateur non seulement comme moyen de visualisation et de stockage mais aussi comme moyen

d'analyse d'image.

Une section sismique peut être considérée comme une image 2D à deux dimensions I(k,1) ou l'indice k représente en ordonnée le temps ou la profondeur et l'indice 1 représente en abscisse le numéro de la trace sismique, I étant l'amplitude fonction de la valeur du coefficient de

réflexion au niveau du réflecteur.

A supposer que l'on représente I, l'amplitude, par différents niveaux de gris, la détection d'horizon sismique pourrait théoriquement se ramener à un problème de détection de bord ou "contouring". Cette approche se heurte cependant à un problème global: les images sismiques sont extrêmement bruitées et les méthodes classiques de traitement d'image s'avèrent très inadaptées à leur analyse. Ainsi les résultats obtenus après application de méthodes classiques de détection de contours ou de "pointés d'extrema" sont très décevants. Cette dernière technique, que l'on utilise classiquement en traitement de données géophysiques, par exemple pour la détermination des pendages, consiste, à partir du pointé d'un maximum d'amplitude sur une trace caractérisée par un temps To, à chercher sur la trace voisine le maximum associé au temps le plus proche de To, et ainsi de suite de trace en trace. Cette méthode, très sensible au bruit, nécessite de plus la donnée de nombreux

seuils par l'interprétateur.

La prise en compte de certaines caractéristiques particulières à l'image d'une section sismique, comme son bruit ambiant ou encore l'unidirectionnalité de l'information portée par elle, a permis l'élaboration d'un procédé d'analyse d'image à deux dimensions représentant des valeurs prises par un paramètre suivant deux axes, lesdites, valeurs s'organisant suivant des linéations, particulièrement adapté aux images sur lesquelles les

linéations sont préférentiellement distribuées suivant un-

des axes, et permettant de déterminer les zones sur

lesquelles les linéations sont de caractéristiques voisines.

Le procédé selon l'invention est caractérisé en ce que: - on effectue tout d'abord un pointé des linéations portées par l'image, - on décompose ensuite chacune de ces lingations en segments curvilignes jointifs tels que la distance la plus grande entre un point de ce segment et le segment de droite joignant ses deux extrémités ne dépasse pas une distance limite donnée, - on attribue ensuite à chacun de ces segments au moins une grandeur caractérisant l'ensemble des valeurs du paramètre associées au segment, - on visualise enfin l'ensemble des segments pour lesquels au moins une des grandeurs appartient à un même intervalle

prédéterminé.

Selon une variante, on calcule et on attribue à chacun des alignements constitués de la jonction de segments jointifs associés à au moins une grandeur appartenant à un même intervalle, des paramètres le caractérisant, afin de visualiser l'ensemble des alignements pour lesquels au moins un des paramètres associés à l'alignement appartient à un même domaine défini par l'utilisateur. Ce procédé trouve évidemment une application particulière à l'analyse structurale et stratigraphique d'une section ou image sismique, pour laquelle l'information -en l'occurrence l'image des horizons sismiques- est préférentiellement distribuée suivant l'horizontale, et permet notamment d'obtenir automatiquement le pointé de tous les horizons sismiques d'une section et de délimiter les ensembles géologiques de caractéristiques voisines. Le procédé selon l'invention, adapté à l'application à la sismique est caractérisé - d'une part en ce que, pour effectuer le pointé des horizons, on applique un programme de "suivi de bord" à l'image de la section sismique après que celle-ci ait été binarisée, la binarisation consistant à remplacer chacune des traces par une fonction binaire ne pouvant prendre que deux valeurs, la première de ces valeurs étant prise si la pente de l'amplitude de la trace est positive, la seconde valeur si la pente de l'amplitude est négative, ce programme permettant de déterminer des contours de zones

associées à l'une ou l'autre valeur.

On décompose ensuite chacun des contours en une suite de chaînes limitées par des points de rebroussement qui sont les points du contour se situant au lieu de changement de direction gauche à droite ou droite à gauche quand on décrit le contour dans un sens donné, lesquelles chaînes

sont considérées comme l'image des horizons sismiques.

- d'autre part, en ce que l'on transforme, avant la décomposition en chaînes, les contours présentant des replis caractérisés, chacun, par deux points de rebroussement rapprochés, en sus des deux points de rebroussement présents sur tout contour, en un contour simple avec deux points de rebroussement seulement et un nombre, multiple de deux, de chaînes isolées, ceci afin de prendre en compte les anomalies qui, sur les contours, peuvent être considérées comme l'image de doublets sismique. Enfin, pour l'application à la sismique, les paramètres caractéristiques attribués aux segments et alignements sont au moins: un pendage, correspondant par exemple à l'angle formé par la ligne droite joignant les deux extrémités du segment ou de l'alignement avec un axe de référence vertical, une amplitude, correspondant à la moyenne des amplitudes optimum le long de la portion de réflecteur associée au segment ou à l'alignement, une variance en amplitude, qui correspond à la variance sur l'ensemble des amplitudes optimum le long de la portion de réflecteur, la longueur, qui correspond au nombre de points du segment ou de l'alignement, et une pseudo-période. Selon des variantes, l'image obtenue après "contouring" peut être filtrée afin d'en éliminer les éléments de trop petite surface, par exemple ne s'étendant

pas sur plus de deux traces.

Le procédé donne bien entendu de bons résultats sur les images sismiques, mais il a également été appliqué avec succès à d'autres types d'images, par exemple à des images VLSI (Very Large Scale Integration), des scènes naturelles telles des images de satellites, SPOT, LANDSAT,

et à la reconnaissance de pièces industrielles.

Un avantage de la présente invention, en particulier dans son application à la sismique, réside dans le fait, outre des gains de temps et de coût évidents, qu'elle permet d'extraire de l'image sismique une

information fiable et robuste liée à l'information sismique.

A titre d'exemple, concernant le gain de temps, une interprétation manuelle de deux mois peut être ramenée, avec ce procédé, à un traitement de quelques jours, voire

quelques heures.

D'autres avantages et caractéristiques

apparaîtront à la lecture de la description d'un mode de

réalisation de la présente invention, se référant aux dessins annexes sur lesquels: - la figure 1 représente l'image d'une section sismique; - la figure 2 représente l'image de la figure 1 binarisée; - la figure 3 représente l'image de la figure 2 après application d'un programme de "contouring"; - la figure 4 représente la superposition du résultat du pointé automatique des extrema locaux et de l'image initiale de la figure 1; - la figure 5 représente la sélection des alignements associés à un pendage négatif; - le schéma de la figure 6 permet de comprendre sur un zoom de section fictive ce que recouvrent certains termes employés; - les figures 7A, 7B montrent sur un exemple fictif la

démarche suivie pour le traitement des doublets.

La section sismique représentée sur la figure 1 a été obtenue par la technique d'exploration dite "sismique réflexion" dont le principe a été donné. Une importante littérature existe sur ce sujet et on pourra s'y reporter pour savoir comment une section sismique est réalisée. Le plus souvent, comme c'est le cas sur la figure 1, et comme on l'a détaillé sur la figure 6, une section sismique est constituée par la juxtaposition d'enregistrements discrets

d'amplitude variable, appelés traces (1).

La représentation analogique la plus courante de ces traces, telle celle schématisée sur la trace la plus à droite de la figure 6, consiste a représenter en noir les arches du signal d'amplitude positive, et en blanc les arches d'amplitude négative, l'axe de référence (2) des amplitudes nulles pouvant ne pas être représenté. Les traces appartenant à la section sismique sont référencées dans un plan dont l'axe des abscisses représente des distances X et l'axe des ordonnées des temps T ou des profondeurs P. L'image de section sismique de la figure 1 est composée de colonnes de pixels représentant chacune une trace de la section, un codage en niveaux de gris permettant d'apprécier l'amplitude de la trace, le noir étant associé aux amplitudes fortes et le blanc aux amplitudes faibles. Ce que l'oeil discerne dans une section sismique, ce sont les ensembles d'arches cohérents qui constituent les horizons sismiques. Ce sont les cohérences que l'on recherche. Chacune de ces arches peut être approximée par son extremum local, qui donne une bonne localisation en temps de l'arche et de son amplitude. Ces extrema locaux peuvent être reliés sous certaines contraintes choisies par l'interprétateur pour constituer des chaînes susceptibles de

représenter une portion d'horizon sismique.

La figure 6 montre à titre d'exemple trois chaînes de maxima locaux (4), (5), (6). L'extremum local peut être considéré également comme le lieu de changement du signe de la pente de l'amplitude du signal. La première étape du procédé selon l'invention appliquée à la sismique consiste à remplacer la trace par une fonction en escalier telle la fonction (7) sur la figure 6, prenant la valeur 1 sur les plages temps pour lesquelles la pente d'amplitude est positive et la valeur 0 sur les plages o la pente est négative. On "binarise" ainsi la section, trace par trace, pour aboutir à une représentation similaire à celle représentée sur la figure 2, mettant à l'évidence des "objets" blancs, l'objet étant un ensemble de valeurs 1 entouré par des 0, et de "trous" noirs, le trou étant un ensemble de valeurs 0 entouré par des 1. Dans le cas particulier présenté, les objets sont associés aux pentes positives et les trous aux pentes négatives: Si l'on suit le sens croissant des ordonnées, chaque passage blanc-noir correspond à un maximum local et chaque passage noir-blanc à

un minimum local.

Il peut être avantageux d'appliquer sur cette image un filtre morphologique pour éliminer les surfaces homogènes de très petitedimension, s'étendant par exemple sur une ou deux traces, et qui ne peuvent être que la

traduction de la présence de bruit.

On applique ensuite à l'image binarisée, un programme classique en traitement d'images de "suivi de bord", appelé encore programme de "contouring", tels les algorithmes de "contouring" de Pavlidis ou Rosendfeld ou encore la méthode "couple-bord" développée par N. KESKES, l'inventeur, dans sa thèse, permettant de délimiter les

limites ou "contour" des zones associées à une même valeur.

Un tel contour (8) est schématisé sur la figure 6.

On peut sur chacun des contours fermés ainsi définis, déterminer un certain nombre de points particuliers, au moins au nombre de deux pour les contours simples et non coupes par un bord de section. Ces points se situent aux lieux de changement de direction, -de gauche à droite ou de droite à gauche dans notre exemple- quand on décrit le contour dans un sens de rotation donné, tels les points (9) et (10) de la figure 7A, et sont appelés points

de rebroussement.

A partir d'un même point de rebroussement, le suivi du contour dans un sens de rotation particulier (sens des aiguilles d'une montre ou sens inverse des aiguilles d'une montre) permet de décrire l'objet tandis que le suivi en sens contraire permet de décrire le trou. Le sens de rotation et la direction du mouvement (droite à gauche ou gauche à droite) permettent de dissocier sans ambiguité si la frontière suivie est associée à un maximum ou à un

minimum d'amplitude.

Certains contours tel le contour schématisé sur la

figure 7A présentent plus de deux points de rebroussements.

Ils seront, dans l'application sismique du procédé selon l'invention, traités de façon particulière car ces points de rebroussement supplémentaires traduisent la présence de replis (13) du contour qui peuvent être significatifs de ce

qu'en sismique réflexion on appelle des "doublets".

Ces doublets (14) peuvent être liés è des phénomènes importants du point de vue de l'interprétation: par exemple présence de banc carbonaté, d'huile ou de gaz et

il est important de les traiter avec soin.

Dans le procédé selon l'invention, on caractérise chacun des replis présents sur la section binarisge par un certain nombre de paramètres qui sont: - les distances dx1 et dx2 qui existent respectivement entre chacun des 2 points de rebroussement appartenant au repli et le point le plus proche du contour auquel le repli appartient et qui est également situé sur la verticale passant par le point de rebroussement; - la distance horizontale d qui sépare les deux verticales passant par les deux points de rebroussement associés au même repli, et qui définit la taille du repli; - éventuellement la distance D (non représente sur les figures) qui sépare le repli du repli suivant rencontré en décrivant le contour dans un sens de rotation et une

direction particuliers, s'il y en a un.

La prise en compte ou la non prise en compte du

repli comme image de doublet est effectuée lors de la phase.

de décomposition en "chaînes" de chacun des contours. On entend par "chaîne" l'ensemble des maxima ou minima locaux associés à la portion de contour qui est comprise entre deux points de rebroussement. Si le contour étudié est celui d'un objet sans repli (avec les conventions de signe et valeurs données jusqu'à présent en exemple), comprenant seulement deux points de rebroussement, il sera constitué de deux chaînes, une chaîne supérieure associée à des minimas locaux

et une chaîne inférieure associée à des maximas locaux.

Le contour comportant un repli sera, avant traitement des doublets, constitué de quatre chaînes, telles les chaînes (15), (16), (17) et (18) de la figure 7A. Le traitement des doublets consiste à valider ou non les points de rebroussement du repli en tant qu'extrémité de chaîne. Si le repli n'est pas considéré comme l'image d'un doublet la chaîne est interrompue au niveau du repli et deux chaînes telles les chaînes (18) et (16) de la figure 7A seront étudiées indépendamment pour la suite du procédé. Si le repli est considéré comme l'image d'un doublet, la chaîne n'est pas interrompue au niveau du repli mais son tracé est cependant modifié; ainsi qu'on l'a représenté sur la figure 7B, on peut définir pour chaque repli trois portions de contour qui se superposent entre les verticales passant par les deux points de rebroussement le caractérisant, à savoir la portion'de contour (19) supérieure définie par les temps les plus faibles, et appartenant à la chaîne supérieure(18), la portion de contour (21) inférieure définie par les temps les plus élevés et appartenant à la chaîne inférieure (16) du contour et la portion de contour milieu (20), comprise entre les deux précédentes et qui se superpose avec une

chaîne entière.

En présence de repli, on considère qu'il existe d'une part une chaîne majeure (22), qui est le résultat de la réunion de la chaîne supérieure (18), à laquelle appartient la portion de contour supérieure (19), avec la partie de la chaîne inférieure qui ne comprend pas la portion de contour inférieure, et d'autre part deux chaînes simples constituées l'une par la portion de contour (20)

milieu, l'autre par la portion de contour inférieure (21).

Les critères de sélection de doublet porte sur les paramètres dx1 dx2, d, et D. Dans l'exemple de réalisation décrit, un repli caractérisé par ces quatre paramètres est considéré comme l'image d'un doublet si l'une au moins des distances dx1 et dx2 est inférieure à un seuil fixé par l'interprétateur fonction du contenu fréquentiel de la section sismique par exemple de l'ordre de 16 ms pour une sismique acquise entre 10 et 80 Hz- et si la longueur d est inférieure à une distance adaptative fonction de la distance D qui sépare le doublet traité du doublet suivant, par exemple la moitié de la distance D. Tous les replis supérieurs ou inférieurs d'un même contour sont traités successivement, par ordre de taille croissante, la distance D étant recalculée après chaque

traitement d'un doublet.

Tous les contours étant décomposés en chaînes, chaque chaîne est ensuite décomposée en segments, un segment étant une portion de chaîne telle qu'aucun point du segment n'est distant de plus d'une valeur h de la droite reliant ses deux extrémités, h étant un seuil fixé par

l'interprétateur, de l'ordre de quelques millisecondes.

La raison de cette décomposition réside dans le fait qu'il est impossible de garantir qu'une chaîne appartienne à un seul horizon sismique et qu'il importe de définir les portions homogènes d'horizon sachant que les points de rupture d'homogénéité se sitent sur les points de

forte courbure des chaînes.

Il existe, pour effectuer ce traitement que l'on peut considérer comme une approximation polygonale, des algorithmes de traitement classiques tel l'algorithme de Pavlidis, fournissant un résultat indépendant du sens dans

lequel on décrit la chaîne.

De façon avantageuse, dans le procédé selon l'invention, en particulier dans son application à la sismique, la sélection et la visualisation de la sélection est effectué automatiquement à l'aide d'un micro-ordinateur et d'un écran couleur sur lequel on visualise la section et les résultats des différentes étapes du procédé, de façon interactive. Chaque segment est référencé en mémoire machine et un certain nombre de paramètres caractérisant les arches de la portion d'horizon sismique qui lui est associé sont calculés et mémorisés. Ces paramètres sont, par exemple; une longueur, correspondant au nombre de traces interceptées par le segment, - un pendage local, correspondant à l'angle formé par l'axe vertical de la section et le segment de droite joignant les deux extrémités du segment, - une amplitude: l'amplitude qui est un des paramètres les plus perceptibles sur une section sismique est prise en compte sous la forme d'une moyenne des amplitudes optimum des arches le long du segment, et de l'écart type associé, - une pseudo période moyenne le long du segment et

l'écart-type associé.

En deux dimensions la seule existence des segments permet à l'interprétateur de visualiser les éléments sismiques existants dans son image. Le procédé selon l'invention permet à l'interprétateur, après lecture d'une bande magnétique, d'observer après quelques minutes de calcul l'influence et la répartition des paramètres tels

l'amplitude, la continuité, le pendage, la fréquence etc...

Il peut visualiser l'intégralité des chaînes, pour obtenir une image similaire à celle de la figure 3, ou les faire apparaître en surimpression de. la sismique, pour

obtenir une image similaire à celle de la figure 4.

La décomposition en segments, et le calcul de paramètres caractérisant chacun. des segments, permet à l'interprétateur de sélectionner, dans son image sismique, l'ensemble des segments qui sont de caractéristiques voisines et donc de visualiser les ensembles sismiques homogènes. La sélection s'effectue en fonction de critères de distance portant sur les paramètres calculés, nécessitant la donnée de seuils et limites par l'interprétateur. La figure 5 montre à titre d'exemple la sélection des segments de la section de la figure 1 caractérisés par un pendage positif. Certains des alignements visibles sur cette figure 5 sont le résultat de la juxtaposition de différents segments adjacents satisfaisant ce même critère. Ces alignements constituent des ensembles homogènes d'extrêma

locaux appartenant à un même horizon.

Il est possible de recalculer des paramètres caractéristiques tels que la longueur, l'amplitude moyenne, la variance en amplitude, la fréquence apparente ou la pseudo période moyenne sur les alignements qui sont déjà le

résultat d'une sélection de premier ordre.

Les critères de sélection peuvent être plus sophistiqués et porter sur plusieurs paramètres simultanément. L'interprétateur peut ainsi sélectionner les ensembles de segments (donc de portions d'horizons sismiques) caractérisés par une amplitude appartenant à une certaine tranche de valeurs données et par des variations de pendage, entre segments adjacents, ne dépassant pas une

autre valeur donnée.

Le codage en couleurs permet d'améliorer la perception de la stratigraphie du milieu géologique étudié et de ses caractéristiques. Il permet par exemple la représentation de l'intégralité des chaînes sismiques et de différentes classes d'amplitudes définies par l'interprétateur, chacun des ensembles des alignements appartenant à une classe donnée apparaissant dans une couleur particulière, une chaîne pouvant être constituée de

plusieurs alignements appartenant à des classes différentes.

La couleur peut aider également à la perception d'un paramètre supplémentaire à celui surlequel a porté la sélection: la représentation d'alignements codés en couleurs différentes en fonction de leur longueur peut par exemple être utile. Dans ce cas à chaque alignement est

associé une couleur particulière et une seule.

En pratique, les contours, dès la première phase du procédé, sont très allongés, certains traversant toute la largeur de la section, et sont assez souvent, du moins pour les horizons bien marqués de forte amplitude, composés de la réunion de seulement deux ou trois alignements, en fonction

du ou des critères choisis.

La continuité des pointés d'horizon peut encore être améliorée par la jonction des alignements situés, dans un contexte spatial d'horizons continus, dans les prolongements l'un de l'autre et dont les extrémités ne sont

séparés que par deux ou trois traces par exemple.

La description détaillée ci-dessus qui porte sur

l'application du procédé selon l'invention à l'analyse de section sismique-ne limite pas l'invention à cette seule application. De nombreuses variantes sont accessibles à l'homme de l'art, suivant les applications envisagées, sans

pour cela sortir du cadre de l'invention.

Le procédé selon l'invention a par exemple été adopté avec succès à d'autres domaines que la géophysique, et a été utilisé pour le comptage de cellules, la détection

simple de contours, l'étude du réseau de peau etc...

Si l'information est distribuée suivant une direction plus proche de l'horizontale que de la verticale, on travaillera colonne par colonne et si l'information est distribuée suivant une direction plus proche de la verticale

que de l'horizontale, on travaillera ligne par ligne.

Claims (3)

REVENDICATIONS

1 - Procédé d'analyse d'image è deux dimensions représentant des valeurs prises par un paramètre suivant deux axes, lesdites valeurs s'organisant suivant des lingations, particulièrement adapté aux images sur lesquelles les linéations sont préférentiellement distribuées suivant un des axes, permettant de déterminer les zones sur lesquelles les lingations sont de caractéristiques voisines, caractérisé en ce que: - on effectue un pointé des linéations portées par l"image, - on décompose chacune de ces linéations en segments curvilignes jointifs tels que la distance la plus grande entre un point de ce segment et le segment de droite joignant ses deux extrémités ne dépasse pas une distance limite donnée, - on attribue ' chacun de ces segments au moins une

grandeur caractérisant l'ensemble des valeurs du.

paramètre associées au segment, - on visualise l'ensemble des segments pour lesquels au moins une des grandeurs appartient à un même

intervalle prédéterminé.

2 - Procédé selon la revendication 1 caractérisé en ce qu'on calcule et on attribue à chacun des alignements constitués de la jonction de segments jointifs associés à au moins une grandeur appartenant au même intervalle, des paramètres le caractérisant, afin de visualiser l'ensemble des alignements pour lesquels au moins un des paramètres associés à l'alignement appartient à un même

domaine défini par l'utilisateur.

3 - Application du procédé selon la revendication 1 ou 2 à l'analyse structurale et stratigraphique d'une section sismique dont l'image est composée de colonnes représentant chacune des traces de la section, permettant notamment d'obtenir automatiquement le pointe de tous les horizons sismiques d'une section et de délimiter les ensembles géologiques de caractéristiques voisines, caractérisé en ce que d'une part, pour effectuer le pointé des horizons, - on applique un programme de suivi de bord (ou "contouring") à l'image de la section sismique après que celle ci ait été binarisée, la binarisation consistant 3 remplacer chacune des traces par une fonction binaire ne pouvant prendre que deux valeurs, la première de ces valeurs étant prise si la pente de l'amplitude de la trace est positive, la seconde valeur étant prise si la pente de l'amplitude de la trace est négative, ce programme permettant de déterminer les contours des zones associées à l'une ou l'autre valeur sur la section binarisée, on décompose chacun des contours en une suite de chaînes limitées par les points de rebroussement qui sont les points du contours se situant au lieu de changement de direction gauche à droite ou droite à gauche quand on décrit le contour dans un sens de rotation donné, lesquelles chaînes sont considérées comme l'image des horizons sismiques, d'autre part en ce que l'on transforme, avant la décomposition en chaînes, les contours présentant des replis caractérisés, chacun, par deux points de rebroussement rapprochés, en sus des deux points de rebroussement présents sur tout contour, en un contour simple avec deux points de rebroussement seulement et un nombre, multiple de deux, de chaînes isolées, et en ce que les paramètres caractéristiques attribués aux segments et alignements sont au moins: un pendage, une amplitude, une variance en amplitude, une longueur, une pseudo-période.