FR 2939520 A1 20100611 FR 0858439 A 20081210 -1 � Procédé de modélisation géologique de données sismiques par corrélation de traces � Domaine technique La présente invention concerne un procédé pour élaborer de façon automatique un modèle géologique à partir de données sismiques. Elle trouve des applications particulièrement pertinentes dans le domaine de l'analyse de données sismiques, notamment lors de la recherche de nouveaux gisements pétroliers. Etat de la technique antérieure D'une façon générale, les données sismiques représentent un sous sol et permettent de ce fait d'appréhender la structure géologique de la zone recouverte. Un signal sismique représente typiquement la dérivée de l'impédance acoustique des couches traversées par des impulsions sonores émises au cours de l'acquisition desdites données sismiques. De ce fait, la sismique est une représentation indirecte et dégradée de la géologie sous jacente. Indirecte car la valeur du pixel n'est pas représentative de la nature même du sous-sol en ce point, et dégradée car la donnée finale est le fruit d'une chaîne de traitements (déconvolution, � stacking process �) qui introduisent des artéfacts et des approximations. Par conséquent, la reconstruction de la géologie à partir de la sismique passe en général par des moyens d'analyse complexes qui cherchent à identifier les couches géologiques une par une et indépendamment les unes des autres. Les méthodes d'analyse existantes utilisent comme postulat qu'une même couche géologique présente des signatures sismiques similaires. Les surfaces sont donc extraites par le biais de propagateurs qui démarrent à partir d'un germe, ou par des classifieurs (Kmean) qui regroupent les signatures sismiques de même calibre. La principale limitation de ces techniques découle du postulat de base lui-même qui n'est jamais vraiment respecté. De plus, les algorithmes employés présentent de nombreux inconvénients: influence trop grande du germe, dérives des propagations, surfaces invalides qui connectent différents étages géologiques... Plus récemment sont apparues des approches globales qui visent à calculer au cours d'un même processus la totalité du modèle géologique directement à partir des données sismiques. Parmi ces méthodes deux 2939520 -2- grandes familles se distinguent : celles qui utilisent le pendage local des couches et celles qui utilisent des primitives. Dans la famille des méthodes utilisant le pendage local des couches, on connaît notamment le document US 6,771,800 de Keskes et al. (2002) qui 5 repose sur la propagation de lignes 2D suivant la direction du pendage. Cette technique ne fonctionne qu'en 2D et ne gère pas les ruptures de continuités de la géologie. On connaît également le document de Lomask et al., � Flattening without picking �, Geophysics vol. 11 no. 4, 2006, pp.13-20, qui propose une 10 méthode de compensation des pendages par un processus itératif de minimisation. Au cours de ce processus, les pendages, obtenus par calcul de gradient 3D sur chaque point sismique, sont peu à peu annulés jusqu'à obtenir un bloc sismique plat, sans déformations globales, qui permet de mettre en évidence les taux de sédimentation par analyse stratigraphique. 15 Cette technique fonctionne en 3D mais la gestion des ruptures du modèle nécessite la connaissance a priori des discontinuités. Dans la famille de méthodes utilisant des primitives comme éléments de base du modèle géologique, ces éléments sont, dans un premier temps, extraits de l'image sismique et sont ensuite agencés les uns par rapport aux 20 autres pour construire le modèle. Dans la méthode décrite par Monsen et al. dans le document US 2008/0140319 notamment, les primitives sont des surfaces géologiques de tailles variables obtenues automatiquement par propagation sur les extrema de la sismique. Ces primitives sont ensuite triées par ordre topologique sur la 25 base d'un graphe acyclique (� Directed Acyclic Graph �, DAG) qui représente les relations de superposition entre les primitives. L'ordre topologique affecté à chaque surface permet de construire le modèle géologique. Cette méthode fonctionne en 3D, elle gère les discontinuités de type érosion mais ne permet pas la gestion automatique des failles. 30 Une variante de cette approche est proposée par Verney et al. dans le document � An Approach of Seismic Interpretation Based on Cognitive Vision �, 70th EAGE Conference & Exhibition - Rome, Italy, 9 - 12 June 2008. Dans ce document, les relations entre mini surfaces sont hiérarchisées en fonction de leur distance verticale moyenne et que les critères de fusion entre 35 primitives reposent sur des principes provenant de la science cognitive. 2939520 -3- Dans le contexte actuel d'un marché pétrolier de plus en plus concurrentiel et cher, la recherche de nouveaux gisements requiert l'usage d'outils d'analyse automatiques plus précis et rapides. La présente invention a pour but un procédé travaillant à l'échelle du point sismique pour élaborer de 5 manière automatique un modèle géologique à partir de données sismiques, pouvant gérer automatiquement les discontinuités géologiques telles que les failles, et permettant l'implémentation de stratégies multi-échelles. Exposé de l'invention Cet objectif est atteint avec un procédé d'élaboration d'un modèle 10 géologique à partir de données sismiques préalablement collectées, caractérisé en ce qu'il comprend des étapes de : - échantillonnage desdites données sismiques selon au moins un ensemble de traces d'inclinaison prédéterminée constituées chacune de points d'échantillonnage, 15 - analyse de similitudes entre les données sismiques autour desdits points d'échantillonnage, - détermination de liens entre points d'échantillonnage appartenant à des traces différentes à partir de ladite analyse, - constitution du modèle géologique en attribuant à chaque point 20 d'échantillonnage un âge géologique relatif qui est calculé au moins à partir des liens associés audit point d'échantillonnage et de l'âge géologique relatif des autres points d'échantillonnage voisins de la trace incluant ledit point d'échantillonnage. Le procédé selon l'invention permet ainsi d'élaborer un modèle 25 géologique à partir d'une grille d'échantillonnage arbitraire des données sismiques, en exploitant les similitudes locales des données. En particulier, il ne repose pas sur la détection d'extrema et la propagation de surfaces à partir de ces extrema comme les procédés de l'art antérieur, mais il permet une modélisation plus fine de l'ensemble des données. 30 Le procédé peut en outre avantageusement comprendre une étape de remplissage du modèle géologique par propagation d'une mini-surface par similarité dans le voisinage de chaque point d'échantillonnage comprenant un lien, en prenant pour germe ledit point. -4- Ces mini-surfaces ne doivent toutefois pas être confondues avec les surfaces propagées comme primitives dans des procédés de l'art antérieur. En effet : - elles ne sont utilisées que pour compléter le modèle géologique, et en 5 particulier les relations topologiques existant entre elles ne sont pas exploitées, - la forme de leur contour et leur taille sont essentiellement déterminées par la grille d'échantillonnage, - du fait de leur surface très réduite, les problèmes de dérives liés aux 10 algorithmes de propagation sont négligeables. Les données d'un volume sismique peuvent être échantillonnées dans un plan sensiblement horizontal selon un maillage sensiblement de l'une quelconque des formes suivantes : rectangulaire, hexagonal, et les traces peuvent être d'orientation sensiblement verticale et comprendre des points 15 d'échantillonnage sensiblement équidistants. Avantageusement, la détermination de liens entre points d'échantillonnage appartenant à des traces différentes peut comprendre une étape de calcul d'images de corrélation, - chaque image de corrélation correspondant à un couple de traces, 20 - chaque trace étant représentée selon une dimension de l'image de corrélation, et - la valeur de chaque point de l'image représentant la corrélation locale des données sismiques au voisinage des deux points respectifs de chaque trace. 25 Le procédé selon l'invention peut comprendre en outre des étapes de : - construction de segments dans l'image de corrélation par propagation parmi des points correspondant à des valeurs de corrélation supérieures à au moins un seuil prédéterminé, lesquels segments correspondant à des ensembles de liens ou peignes de corrélations entre points d'échantillonnage 30 des traces ; - sélection des segments de l'image de corrélation, laquelle sélection ne laissant au plus qu'un segment par ligne et par colonne de l'image de corrélation, correspondant à un seul lien par point d'échantillonnage, et les segments étant sélectionnés sur la base de la valeur des points les 35 constituants, et de leur longueur ; 2939520 -5- - validation des liens par l'analyse des chemins constitués par les successions de points reliés par lesdits liens, une condition pour que les liens soient validés étant, dans un volume sismique, que les chemins soient bouclés. Avantageusement, les étapes de sélection de segments dans l'image de corrélation et de validation des liens peuvent être répétées de manière itérative sur les segments restants dans l'image de corrélation, en utilisant comme contrainte les liens validés aux itérations précédentes, jusqu'à ce qu'il n'y ait plus de segments exploitables. Avantageusement, l'attribution d'un âge géologique relatif aux points d'échantillonnage peut être un processus itératif exécuté tant qu'il reste des points d'échantillonnage non datés comprenant au moins un lien, qui comprend des étapes de : - sélection d'un point d'échantillonnage non daté comprenant au moins un lien, - attribution d'un âge relatif à ce point en tenant compte des points déjà datés de sa trace, de telle sorte que l'âge le long d'une trace varie toujours dans le même sens, et - attribution du même âge relatif à tous les points reliés à ce point par des liens. Le procédé selon l'invention peut comprendre en outre une étape de normalisation des âges relatifs, l'incrément d'âge entre points d'âge relatif consécutifs étant fixé à une valeur prédéfinie. Suivant un autre aspect du procédé selon l'invention, un modèle géologique multi-échelles peut être élaboré en fusionnant des ensembles de liens entre points d'échantillonnage, lesquels ensembles de liens étant obtenus avec des grilles d'échantillonnage différentes. Avantageusement, le procédé selon l'invention peut être appliqué de manière automatisée, notamment : - à l'extraction directe des surfaces de dépôt sur la base d'iso-valeurs, - à l'horizontalisation de la sismique en cartographiant sur une ligne ou un plan le signal sismique correspondant à un étage géologique, - à l'analyse de failles, et - à l'analyse de chenaux, 2939520 -6- Suivant un autre aspect de l'invention, il est proposé un programme d'ordinateur mettant en oeuvre le procédé suivant l'invention. Description des figures et modes de réalisation D'autres avantages et particularités de l'invention apparaîtront à la 5 lecture de la description détaillée de mises en oeuvre et de modes de réalisation nullement limitatifs, et des dessins annexés suivants : - la figure 1 illustre la construction de l'image de corrélation correspondant à deux traces de l'ensemble des données sismiques. La figure la illustre l'échantillonnage vertical des données en traces sismiques tandis que la figure 10 lb illustre la construction de l'image de corrélation ; - la figure 2b illustre le peigne de corrélation correspondant à une géologie dite iso pack avec des variations de pendage telle que illustrée à la figure 2a ; - la figure 3b illustre le peigne de corrélation correspondant à une géologie présentant à la fois des variations de pendage et d'épaisseur telle que 15 illustrée à la figure 3a ; - la figure 4 illustre l'établissement des peignes de corrélation entre points de part et d'autres d'une faille présentant un rejet normal (figure 4a) ou inverse (figure 4b) ; - la figure 5 illustre un processus d'extraction des segments de droite de 20 forte intensité dans l'image de corrélation ; - la figure 6 illustre des exemples de stratégies d'échantillonnage du cube sismique ; - la figure 7 illustre la validation des peignes de corrélation par analyse des chemins ; 25 - la figure 8 illustre le principe de l'insertion itérative de peignes de corrélation pour remplir l'espace des données ; - la figure 9 illustre un exemple d'application de l'algorithme de calcul des positions de référence ; - La figure 10 illustre un exemple de mini-surfaces obtenues dans une 30 grille d'échantillonnage carrée ; - la figure 11 illustre le remplissage du modèle géologique ; Le procédé selon l'invention exploite un ensemble de données sismiques pour en déduire un modèle géologique, et ceci en travaillant à l'échelle des points sismiques. -7- Les données telles que décrites dans le mode de mise en oeuvre préférentiel font référence à un cube ou volume 3D de données sismiques. Le procédé selon l'invention peut toutefois également être mis en oeuvre sur un ensemble de données sismiques bidimensionnel, tels que des coupes. Le volume de données sismiques peut être caractérisé par 3 axes, soit respectivement des axes Inline et Xline sensiblement parallèles à la surface du sol et un axe TWT orienté dans le sens des profondeurs croissantes. Le procédé selon l'invention peut être appliqué à la totalité des données du volume sismique, mais dans le cas général il est suffisant de se limiter à un sous ensemble de points échantillonnés pour limiter notamment le temps de calcul. En référence à la figure 1, une trace sismique ou trace 3, 4 comprend un ensemble de points 5 sensiblement alignés et équidistants, échantillonnés aux profondeurs ti, t2, t3, ... dans la direction de l'axe TWT. La totalité du volume sismique peut être échantillonné par un ensemble de traces disposées dans le plan Inline-Xline selon un maillage par exemple carré, hexagonal, voire arbitraire, tel qu'illustré à la figure 6. En référence à la figure 1, une première étape du procédé selon l'invention consiste à calculer pour un ensemble des couples de traces (3, 4) d'intérêt extrait des données sismiques 1 une image de corrélation I(X1, X2), X1 et X2 étant les positions respectives des traces sismiques 3, 4. Cette image de corrélation 2 représente l'ensemble des coefficients de corrélations calculés entre les Ni points provenant de X1 et les N2 points provenant de X2. Dans cette image de dimensions (N1*N2) pixels, si X1 est en ordonnée et X2 en abscisse, la valeur du point P(i, j) situé à la colonne i et à la ligne j est le coefficient de corrélation Cf(i, j) calculé entre le ième point de la trace X2 et le jème point de la trace X1. Les coefficients de corrélation Cf(i, j) peuvent être calculés par exemple à partir de la corrélation d'une mini trace Vi de taille S, centrée autour d'un point X2(i) et d'une mini trace sismique Vj de taille S centrée autour d'un point X1(j). Le mode de calcul des coefficients Cf(i, j) peut également et de manière non limitative comprendre tout type de calcul de distance ou de différences entre les deux vecteurs Vi et Vj. Dans la figure 1 par exemple, la valeur du point 8 de l'image de corrélation 2 est le résultat du calcul de la corrélation entre la mini trace 6 autour du point 5, de position t3 dans la trace X2 avec la mini trace autour 2939520 -8- du point de position ti de la trace X1. Par une convention non limitative, le mode de calcul de la valeur des points de l'image de corrélation 2 est défini de telle sorte que cette valeur augmente avec la similitude des vecteurs Vi et Vj. Chaque point de l'image de corrélation I(X1, X2) tel que 8 correspond 5 donc à un lien, par exemple 7, entre un point de la trace X1 et un point de la trace X2. Un point qui présente une forte valeur dans I(X1, X2) représente un lien de forte probabilité entre X1 et X2. L'image de corrélation met donc en évidence tous les liens de forte probabilité entre deux traces sismiques. En référence à la figure 2, lorsque les points de forte probabilité sont contigus 10 dans l'image de corrélation 2, leur pertinence augmente car ils ne sont plus isolés dans le contexte géologique mais forment un ensemble consistant 10 qui lie deux étages géologiques. Ces ensembles contigus correspondent à des peignes de corrélation 11 dans l'espace des données 1, lesquels peignes sont définis par les ensembles E de points P(i, j) de l'image I(X1, X2) tels 15 qu'aucun point P(i j) de E ne présente la même abscisse ni la même ordonnée qu'un des autres points de E. Cette propriété traduit le fait que dans un peigne donné, un point sismique n'est connecté qu'à un seul lien. Dans le cas d'une géologie complètement laminaire à pendage nul, les traces X1 et X2 sont similaires et le segment de plus forte corrélation se situe 20 sur la diagonale de l'image 2, avec un coefficient directeur sensiblement égal à 1. Le peigne résultant comprend alors des liens horizontaux. Dans le cas d'une géologie dite iso pack avec des variations de pendage telle qu'illustrée à la figure 2, les deux traces X1 et X2 présentent des signatures similaires mais sont décalées. Le segment de plus forte corrélation 25 10 est alors parallèle à la diagonale, avec un coefficient directeur sensiblement égal à 1, et la distance 12 entre la diagonale et le segment représente la valeur du rejet, c'est-à-dire du décalage suivant l'axe TWT, entre X1 et X2. De même, dans le cas plus général où la géologie présente à la fois des 30 variations de pendage et d'épaisseur, tel qu'illustré à la figure 3, le segment de plus forte corrélation 10 possède un coefficient de direction différent de 1 car le rejet varie le long du peigne 11. Avantageusement, le procédé selon l'invention n'est pas sensible à la présence de failles, contrairement aux méthodes de propagation utilisées 35 dans l'art antérieur pour obtenir des surfaces sismiques et des primitives. En 2939520 -9- utilisant la méthode des peignes de corrélation, des relations ou liens 7 entre les points peuvent être établies de part et d'autre de la faille 20, quelque soit la valeur et le sens (normal ou inverse) du rejet, tel qu'illustré à la figure 4. Les peignes 11 peuvent être détectés automatiquement en extrayant des 5 segments de droite 10 de forte intensité de l'image de corrélation I(X1, X2). Cette détection peut avantageusement être réalisée par exemple sur la base d'un algorithme comprenant les étapes suivantes, illustrées aux figures 5a à 5f, respectivement : - classification des points de l'image I(X1, X2) en trois familles : 10 INVALID, VALID et INCERTAIN en fonction de leur valeur de corrélation (figure 5a) ; - utilisation de tous les points classés VALID comme germe pour propager des segments. Les points candidats pour les propagations sont uniquement les points classés VALID et INCERTAIN (figure 5b) ; 15 - recherche de tous les points de I(X1, X2) qui présentent la valeur maximale sur une ligne ou sur une colonne. Ces points sont marqués MAXP (figure 5c) ; - élimination des segments qui ne contiennent aucun point marqués MAXP (figure 5d) ; 20 - élimination des segments dont la longueur est inférieure à la longueur médiane des segments obtenus dans l'étape précédente (figure 5e) ; - présentation de chaque segment en partant du plus long au plus court. Chaque segment présenté élimine tous les segments qui ont un point sur la même ligne ou sur la même colonne de l'image de corrélation (figure 5f). On 25 obtient ainsi le plus long segment 10 correspondant au plus long peigne 11 en s'assurant que les points sismiques du peigne 11 ne sont associés qu'à un seul lien. Pour détecter les peignes 11 sur l'ensemble des données du bloc sismique, il faut choisir des positions de couples de traces X1 et X2 à partir 30 desquels une image de corrélation est calculée. Ces traces peuvent être positionnées n'importe où dans le bloc sismique. Si les capacités machines le permettent, on peut utiliser toutes les traces du cube sismique et donc tous les points. Cependant, en règle générale, un échantillonnage du cube sismique est nécessaire. 2939520 -10- Avantageusement, on peut opter pour une stratégie d'échantillonnage régulière dans laquelle on prend une trace sur N traces possibles et on calcule toutes les images I(Xi, Xj) tel que les positions Xi et Xj soient distantes de N, suivant par exemple un maillage sensiblement carré ou hexagonal (figures 5 6a et 6b). On peut aussi mettre en oeuvre une stratégie aléatoire dans laquelle tous les couples (Xi, Xj) sont tirés aléatoirement (figure 6c). On peut également avantageusement combiner plusieurs stratégies, la seule finalité étant d'obtenir un nombre suffisant de peignes 11 pour que l'essentiel du bloc sismique soit caractérisé. 10 Cette souplesse d'usage confère au procédé selon l'invention un avantage considérable car elle permet de mettre en oeuvre facilement des stratégies multi-échelles, ou multi-résolutions, en combinant par exemple des peignes distants de N échantillons avec des peignes distants de 2*N échantillons. L'utilisation des peignes à plusieurs échelles ou niveaux de 15 résolution permet de mieux traiter par exemple des facteurs géologiques à grande échelle tels que des zones de failles ou d'érosion. Avantageusement, le procédé selon l'invention peut comprendre une étape de validation des liens, sur la base d'une analyse des chemins constitués par les suites de points sismiques reliés entre eux par ces liens. En 20 effet, dans la mesure où un lien relie par définition des points qui se situent sur une même surface géologique, en parcourant un chemin de lien en lien on doit normalement rester sur la même surface. Si un chemin ne se reboucle pas sur lui même il est rejeté car il n'y a aucune certitude au regard de sa pertinence. Si un chemin comprend deux points différents appartenant à une 25 même trace 31, tel qu'illustré à la figure 7b, il n'est également pas validé car il traverse différents étages géologiques lors de son parcours. Si, au contraire, le chemin revient exactement au point de départ 30 sans recouvrement (figure 7a), il est validé car la probabilité qu'il soit sur la même surface est forte. Lorsqu'un chemin est validé, tous des liens qui le composent sont 30 également validés. De même, lorsqu'un certain nombre de liens est validé dans un peigne, la totalité du peigne est validé car on considère que les peignes sont des ensembles indissociables et pertinents. Bien que ce critère de validation des liens soit exigeant, il permet de remplir l'espace sismique par des peignes ayant une consistance géologique. -11- Avec cette technique un lien valide peut parcourir des surfaces géologiques présentant de fortes variations de pendage et traversées par des failles. Il peut être avantageusement mis en oeuvre dans le procédé selon l'invention un processus itératif pour remplir au maximum le volume sismique étudié, de telle sorte à bâtir un modèle final le plus complet possible avec le maximum de liens. Chaque itération de ce processus comprend une détection puis une validation d'un ensemble de peignes tels que cela est décrit précédemment. Chaque nouvelle itération est contrainte par les itérations précédentes. L'espace sismique se remplissant peu à peu, le nombre de possibilités d'insertion de nouveaux peignes diminue. Soit un couple (X1, X2) de deux traces. Si au cours de la première itération, un peigne a été validé entre X1 et X2, ce dernier servira de contrainte pour l'itération suivante car il limite le champ de solutions disponibles pour un second peigne. Ce principe est illustré à la figure 8. Dans cette figure, les segments correspondant à trois peignes pertinents sont présents dans l'image de corrélation. Le segment 40 correspondant au peigne 43, le plus évident, est validé lors de la première itération. Pour la seconde itération, toute la partie grisée de l'image n'est plus disponible conformément au principe selon lequel il ne peut y avoir qu'un lien par couple de points. Par conséquent, le segment 42, dont une trop grande partie des liens recouvre dans la zone grisée, est éliminé de la liste des candidats. Seul le segment 41 correspondant au peigne 44 reste disponible. Les premiers peignes valides se positionnent sur la sismique de bonne qualité, qui présente les meilleures corrélations. Ces peignes vont ensuite servir de contrainte pour la deuxième vagues de peignes, lors de la deuxième itération, qui va utiliser l'espace disponible avec une sismique de moins bonne qualité. Ainsi, la contrainte géologique apportée par les liens existants va s'enrichir peu à peu de façon à compenser la perte graduelle de qualité de la sismique qui reste non traitée. Avantageusement, l'espace est peu à peu recouvert même dans les zones présentant un faible rapport signal/bruit. Une condition d'arrêt du processus itératif peut être par exemple l'absence de nouveaux peignes détectés. Le procédé selon l'invention permet ainsi de remplir l'espace des données de peignes de corrélations reliant les points des couples de traces (Xi, Xj). Soit Ep, l'ensemble des points sismiques qui se situent sur tous les 2939520 -12- couples (Xi, Yj). Avantageusement, il est possible d'affecter à chaque point P de Ep une position de référence calculée avec les liens validés et qui traduit la position de P dans le contexte géologique. Cette position de référence permet d'établir une hiérarchie des points P, du plus jeune au plus ancien. Ce calcul 5 peut être effectué par exemple sur la base de l'algorithme suivant : ValRef(P) = valeur de référence pour un point P Step est l'incrément des valeurs de position, par exemple Step = 1 10 (Initialisation) Pour tous les points P faire P = UNMARKED insérer P dans GlobalStack. VaIRef(P) = 0. 15 Fin Faire (Calcul des positions) Dépiler GlobalStack dans nextP Tant Que nextP existe Faire 20 P = nextP nextP = UNKNOW insérer P dans LocalStacki et dans LocalStack2 minVaiRef = UNDEFINED et maxValRef = UNDEFINED Tant Que LocalStacki n'est pas vide Faire 25 Pop LocalStacki dans K Si nextP =UNKNOW Alors Chercher PAbove le premier point UNMARKED situé au dessus de K Chercher PBelow le premier point UNMARKED situé au dessous 30 de K Si PAbove existe alors nextP = PAbove Si PBelow existe alors nextP = PBelow Fin Si Chercher PAboveMark le premier point MARKED situé au dessus 35 de K -13- Chercher PBelowMark le premier point MARKED situé au dessous de K Si PAboveMark existe et si ValRef(PAbove)>minValRef alors minValRef = ValRef(PAbove) Si PBelowMark existe et si ValRef(PBelow) Dans le cas du procédé selon l'invention, l'ensemble E peut représenter tous les points de la surface du bloc sismique et S les positions des traces Xi, utilisées pour calculer les images de corrélation. Avantageusement, tous les points d'une mini surface propagée à partir d'un germe G se voient attribuer la position de référence associée à G. Cette technique permet de remplir tout l'espace en respectant les variations de pendage local. Il est à noter que, contrairement aux procédés de l'art antérieur, la méthode ne se limite pas aux extrema de la sismique mais une mini-surface est créée pour chaque point comprennent un lien valide. De même, par comparaison avec les procédés de l'art antérieur qui utilisent des -15- surfaces propagées comme primitives, les mini-surfaces calculées dans le procédé selon l'invention présentent des caractéristiques très particulières : - la forme de leur contour et leur taille sont essentiellement déterminées par la grille d'échantillonnage, - du fait de leur surface très réduite, les problèmes de dérives liés aux propagations sont négligeables, - elles ne sont utilisées que pour compléter le modèle, et en particulier les relations topologiques existant entre elles ne sont pas exploitées. La figure 10 montre un exemple de mini-surfaces obtenues dans une grille d'échantillonnage carrée. Le remplissage du modèle géologique est illustré à la figure 11, (a) : image sismique, (b) : ensemble des liens validés, (c) : mini surfaces propagées à partir des points connectés aux liens validés, (d) : positions de référence obtenues pour chaque point de l'image (a). Selon des caractéristiques avantageuses du procédé selon l'invention, - les modèles obtenus sont très précis car les primitives de référence sont les points sismiques. Dans les procédés de l'art antérieur qui utilisent des surfaces comme primitives, la modélisation réelle n'est possible que sur les extrema, le reste de l'information étant obtenu par interpolation verticale entre les extrema. Dans le procédé selon l'invention aucune interpolation n'est requise, ce qui permet moins d'approximation et donc plus de précision ; - dans les approches de l'art antérieur basées sur des primitives, une surface est indivisible, c'est un ensemble de points indissociables qui ont toujours la même valeur dans le modèle géologique. Cette rigidité écarte toute possibilité de contrôle du modèle au sein d'une même surface. La qualité du modèle est donc tributaire de la qualité d'extraction des primitives. Dans le procédé selon l'invention le modèle peut être contraint en tout point de l'espace sismique en ajoutant simplement un lien. Le modèle conserve toute sa souplesse et le contrôle peut se faire sans limitation ; - il peut être intéressant de combiner plusieurs échelles d'analyse dans l'interprétation des données sismiques. Par exemple, on peut mixer une analyse locale qui s'intéresse aux petites variations du signal sismique avec une analyse plus globale qui regarde les grands événements tel que les forts rejets de failles ou les variations d'épaisseur. Avec les méthodes de modélisation de l'art antérieur, l'usage de plusieurs échelles au sein d'un seul 2939520 -16 et même calcul est très difficile à mettre en oeuvre, contrairement au procédé selon l'invention dans lequel il suffit de faire varier les distances entre les traces sur lesquelles les images de corrélation sont calculées ; - les peignes de corrélations utilisés dans le procédé selon l'invention 5 peuvent trouver des liens entre les couches géologiques situées de part et d'autre des failles et ce quelque soit l'amplitude et le sens du rejet. Le procédé selon l'invention permet ainsi de mieux gérer les discontinuités géologiques que les autres méthodes de l'art antérieur tributaires d'une analyse locale telle que la propagation ou le calcul local du pendage. 10 Le procédé selon l'invention peut avantageusement être appliqué à des analyses automatiques de paramètres géologiques tels que : - l'extraction directe des surfaces de dépôts sur la base d'iso-valeur. Le procédé selon l'invention permet un gain de temps important car il permet de s'affranchir du processus d'extraction manuelle par pointé qui est 15 extrêmement long et pénible. Les surfaces obtenues sont une source d'information très riche, car elles constituent l'élément fondamental et premier pour de nombreux systèmes et méthodes d'analyse de la géologie ; - l'horizontalisation de la sismique qui consiste à cartographier sur une ligne ou sur un plan horizontal le signal sismique correspondant à un étage 20 géologique. En traitant ainsi plusieurs étages successifs on arrive à reconstruire une sismique dite redressée dans laquelle tous les reliefs des dépôts géologiques sont aplanis. Le procédé selon l'invention est avantageusement mis en oeuvre pour ce genre d'applications, très utilisées dans la recherche d'hydrocarbures notamment, car le relief de chaque surface 25 est alors connu avec de surcroît une très bonne précision et une grande fiabilité ; - l'analyse des failles qui sont des glissements des couches géologiques le long d'une surface de fracturation. Cette fracturation engendre une rupture de continuité des dépôts géologiques qui, en sismique, se traduit par une 30 perturbation du signal. Ainsi, dans l'image de modélisation géologique résultante la présence d'une faille se traduira par des variations brutales des valeurs topologiques qui se déclinent le long de surfaces. Ces deux critères rendent possible l'extraction des failles sous forme de surfaces vectorisées et aussi de cartographie l'ampleur du rejet le long des failles ; 2939520 -17- - l'analyse des chenaux qui sont des anciens lits de rivières terrestres ou sous-marines comblés par des sédiments. En exploration pétrolière, ces chenaux constituent des cibles privilégiées car présentent toutes les caractéristiques des pièges à hydrocarbures. Les phénomènes de creusement 5 et de remplissage des lits d'écoulement sont rapides à l'échelle géologique et se traduisent, dans l'image de modélisation géologique obtenue par l'application du procédé selon l'invention, par des variations d'épaisseur brutales des étages géologiques sous forme de lentilles. Ces caractéristiques morphologiques sont assez discriminantes pour envisager une cartographie 10 tridimensionnelle automatique des chenaux. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. REVENDICATIONS1. Procédé d'élaboration d'un modèle géologique à partir de données sismiques préalablement collectées, caractérisé en ce qu'il comprend des 5 étapes de : - échantillonnage desdites données sismiques (1) selon au moins un ensemble de traces (3, 4) d'inclinaison prédéterminée constituées chacune de points d'échantillonnage (5), - analyse de similitudes entre les données sismiques (6) autour desdits 10 points d'échantillonnage (5), - détermination de liens (7) entre points d'échantillonnage (5) appartenant à des traces différentes à partir de ladite analyse, - constitution du modèle géologique en attribuant à chaque point d'échantillonnage (5) un âge géologique relatif qui est calculé au moins à 15 partir des liens (7) associés audit point d'échantillonnage (5) et de l'âge géologique relatif des autres points d'échantillonnage voisins de la trace (4) incluant ledit point d'échantillonnage. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend en 20 outre une étape de remplissage du modèle géologique par propagation d'une mini surface par similarité dans le voisinage de chaque point d'échantillonnage (5) comprenant un lien (7), en prenant pour germe ledit point (5). 25 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que : - Les données d'un volume sismique (1) sont échantillonnées dans un plan sensiblement horizontal selon un maillage sensiblement de l'une quelconque des formes suivantes : rectangulaire, hexagonal, et - les traces (3, 4) sont d'orientation sensiblement verticale et 30 comprennent des points d'échantillonnage (5) sensiblement équidistants. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la détermination de liens (7) entre points d'échantillonnage appartenant à des traces différentes comprend une étape 35 de calcul d'images de corrélation (2), - chaque image de corrélation correspondant à un couple de traces (3, 4), 2939520 -19 - chaque trace étant représentée selon une dimension de l'image de corrélation (2), et - la valeur de chaque point de l'image (2) représentant la corrélation locale des données sismiques (6) au voisinage des deux points (5) respectifs 5 de chaque trace (3, 4). 5. Procédé selon la revendication 4, caractérisé en ce qu'il comprend en outre une étape de construction de segments (10) dans l'image de corrélation (2) par propagation parmi des points (8) correspondant à des valeurs de 10 corrélation supérieures à au moins un seuil prédéterminé, lesquels segments (10) correspondant à des ensembles de liens ou peignes de corrélations (11) entre points d'échantillonnage (5) des traces (3, 4), 6. Procédé selon la revendication 5, caractérisé en ce qu'il comprend en 15 outre une étape de sélection des segments (10) de l'image de corrélation (2), - laquelle sélection ne laissant au plus qu'un segment (10) par ligne et par colonne de l'image de corrélation (2), correspondant à un seul lien (7) par point d'échantillonnage (5) et - les segments (10) étant sélectionnés sur la base de la valeur des points 20 (8) les constituants, et de leur longueur. 9. Procédé selon la revendication 6, caractérisé en ce qu'il comprend en outre une étape de validation des liens (7) par l'analyse des chemins (30, 31) constitués par les successions de points (5) reliés par lesdits liens (7), une 25 condition pour que les liens (7) soient validés étant, dans un volume sismique, que les chemins soient bouclés (30). 10. Procédé selon la revendication 7, caractérisé en ce que les étapes de sélection de segments (10) dans l'image de corrélation (2) et de validation 30 des liens (7) sont répétées de manière itérative sur les segments restants (41, 42) dans l'image de corrélation (2), en utilisant comme contrainte les liens validés (40) aux itérations précédentes, jusqu'à ce qu'il n'y ait plus de segments exploitables. 35 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'attribution d'un âge géologique relatif aux points d'échantillonnage (5) est un processus itératif exécuté tant qu'il reste des -20- point d'échantillonnage non daté (5) comprenant au moins un lien (7), qui comprend des étapes de : - sélection d'un point d'échantillonnage non daté (5) comprenant au moins un lien (7), - attribution d'un âge relatif à ce point (5) en tenant compte des points déjà datés de sa trace (4), de telle sorte que l'âge le long d'une trace varie toujours dans le même sens, et - attribution du même âge relatif à tous les points reliés à ce point par des liens (7). 10. Procédé selon la revendication 9, caractérisé en ce qu'il comprend en outre une étape de normalisation des âges relatifs, l'incrément d'âge entre points (5) d'âges relatifs consécutifs étant fixé à une valeur prédéfinie. 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un modèle géologique multi-échelles est élaboré en fusionnant des ensembles de liens entre points d'échantillonnage, lesquels ensembles de liens étant obtenus avec des grilles d'échantillonnage différentes. 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il est appliqué à l'extraction directe des surfaces de dépôt sur la base d'iso-valeurs. 13. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce qu'il est appliqué à l'horizontalisation de la sismique en cartographiant sur une ligne ou un plan le signal sismique correspondant à un étage géologique. 14. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce qu'il est appliqué à l'analyse de failles. 15. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce qu'il est appliqué à l'analyse de chenaux. 16. Programme d'ordinateur mettant en oeuvre le procédé selon l'une quelconque des revendications précédentes.