FR 2844362 A1 20040312 FR 0211200 A 20020910 La présente invention concerne une méthode pour compacter et filtrer des données sismiques dans le but d'en extraire au mieux des informations sur la nature du sous-sol. Cette méthodologie s'applique à la compaction et au filtrage des jeux de données sismiques "multicubes", i.e. soit des données sismiques avant sommation issues d'une 5 campagne sismique 3D regroupées en un certain nombre de cubes de données correspondant chacun par exemple à un même déport (cubes de traces sismiques dits iso offsets) ou à un même angle d'incidence (dits cubes iso angles d'incidence), ou à des cubes d'attributs sismiques obtenus par traitement des amplitudes sismiques comme, par exemple, des cubes de propriétés élastiques issus d'une inversion stratigraphique jointe de 10 données sismiques avant sommation, ces cubes étant obtenus par la mise en oeuvre de la méthode décrite par exemple dans le brevet FR 2 800 473 (EP 1 096 270) du demandeur. Il s'agit alors de mettre en forme l'information sismique avant sommation contenue dans des portions de traces sismiques formées par exploration d'une zone souterraine, à l'aide de nouveaux attributs représentatifs qui sont ensuite exploités pour l'interprétation géologique 15 ou l'interprétation d'un réservoir du sous-sol. Elle peut aussi s'appliquer à la compaction et au filtrage de données sismiques 3D avant ou après sommation de campagnes répétées dans le temps (sismique 4D), ceci afin d'interpréter les variations de contenu en fluides d'un réservoir pétrolier en cours de production. Etat de la technique Les mesures sismiques sont classiquement utilisées pour fournir des informations complémentaires, par rapport aux données de puits, sur les variations des propriétés géologiques des formations du sous-sol. L'interprétation sismique est basée sur l'analyse d'attributs afin d'extraire ces propriétés géologiques à partir de l'information sismique, et ce dans une zone d'intérêt d'une formation souterraine (zone réservoir par exemple). Une approche standard pour l'analyse d'attributs implique l'utilisation de méthodes statistiques de reconnaissance de forme et d'estimation afin de calibrer les attributs sismiques pertinents avec l'information disponible obtenue au droit de différents puits forés. Par exemple, par le brevet FR -A 2 768 818 (US 6 051 651 du demandeur), on connaît une méthode permettant de détecter les classes potentielles dans une population d'événements liés aux propriétés physiques d'un milieu complexe tel que le sous-sol, repérés à partir de données obtenues par exploration du sous-sol ou par des mesures in situ, 5 (événements d'ordre géologique, géophysique, événements liés à des données de production de fluides etc.), ces événements étant repérés par des points dans un espace multivariable et définis par un vecteur d'attributs. Un autre exemple est donné par le brevet FR 01/05.675 du demandeur, qui décrit une méthode pour faciliter le suivi au cours du temps de l'évolution d'une zone souterraine 10 par analyse comparée d'un certain nombre n de jeux d'enregistrements sismiques obtenus successivement à l'issue de n campagnes sismiques 3D successives (sismique dite 4D), dans laquelle on utilise une technique de reconnaissance des formes appliquée à l'ensemble des événements sismiques de plusieurs campagnes, considérés et analysés simultanément. Citons aussi le brevet EP -A-671 017 (US 5 638 269) du demandeur, qui décrit une 15 méthode permettant de rapprocher des données géologiques obtenues par carottage ou diagraphies dans des puits et des données sismiques obtenues par une campagne d'exploration sismique qui s'appuie sur une technique de calibrage statistique avec un étalonnage statistique obtenu en rapprochant les données géologiques locales (mesurées dans des puits) avec des attributs sismiques lus sur des traces sismiques obtenues au 20 voisinage immédiat de chaque puits. Les attributs sismiques qui font l'objet de cette interprétation sont calculés dans le domaine après sommation, ou dans le domaine avant sommation. Dans le domaine après sommation, les attributs classiques sont calculés à partir des amplitudes au niveau du réservoir, ou à partir de l'impédance P estimée à l'aide d'une 25 inversion stratigraphique. Dans le domaine avant sommation, le nombre d'attributs pertinents peut considérablement augmenter. Dans ce cas, les attributs sont soit les amplitudes à différents offsets (déports) ou angles d'incidence, soit des paramètres issus d'une inversion stratigraphique jointe avant sommation (Impédances P et S, masse spécifique, produit des 30 paramètres de Lamé avec la masse spécifique, etc.). Une méthode dite d'analyse en composantes principales (ACP) bien connue des gens de l'art, peut être appliquée à ces attributs afin d'analyser les relations existantes entre eux et de réduire le nombre d'attributs significatifs. Les composantes principales extraites de l'ACP sont les nouveaux attributs: ils définissent une base vectorielle orthogonale ou 5 non (rotation) et sont des combinaisons linéaires de variables initiales. Ils peuvent être utilisées dans des algorithmes statistiques de reconnaissance de forme et correspondent à un filtrage multivarié de l'information sismique initiale. Différentes applications de l'ACP sont décrites, par exemple, dans: Dumay, J., Fournier, F., 1988, "Multivariate statistical analyses applied to seismic facies 10 recognition", Geophysics, 53, 1151-1159; Hagen, D., C., 1982, "The application of principal components analysis to seismic data sets": Geoexpl., 20, 93-111; ou encore dans le brevet FR-A- 2 772 138 (US 6 345 108) du demandeur. Mais l'ACP ne prend pas en compte la notion de groupe d'attributs et, par 15 conséquent, il est souvent difficile de donner une interprétation physique aux composantes principales, ou bien de les relier clairement aux attributs initiaux, surtout si ceux-ci sont organisés en groupes de même signification physique. L'analyse dite canonique est aussi une méthode statistique bien connue des gens de l'art, qui permet d'établir les relations pouvant exister entre deux groupes de variables afin 20 de savoir si ces deux groupes décrivent les mêmes propriétés. Cette méthode est décrite par exemple par: - Hotelling, H., 1936, "Relations between two sets of variables", Biometrika, 28, 321377. Un exemple d'application de l'analyse canonique est décrite par exemple par: - Fournier, F., and Derain, J.F., 1995, "A statistical methodology for deriving reservoir properties from seismic data", Geophysics, 60, 1437-1450. 4 2844362 Cette méthode se limite à l'étude de deux groupes de variables et donc ne s'applique pas de façon générale à l'information sismique multicubes. De plus, elle définit des variables synthétiques dans chacun des sousespaces associés aux deux groupes, et non une seule base vectorielle permettant de compacter l'ensemble des variables initiales et de décrire les relations entre les deux groupes. Plusieurs méthodes permettant de généraliser l'analyse dite canonique en étendant l'analyse à plus de deux groupes de variables existent. Différents aspects de l'art antérieur dans le domaine considéré sont décrits, par exemple, dans les publications suivantes: - Horst, P., 1961, "Relations among m sets of measures", Psychometrika, 26, n 2, 12910 149; - Caroll, J.D., 1968, "A generalisation of canonical correlation analysis to three or more sets of variables", Proc. 76th Conv. Amer. Psych. Ass; - Kettenring, J.R., 1971, "Canonical analysis of several sets of variables", 58, 3, 433450; ou - Saporta, G., 1990, "Probabilités, analyse des données et statistiques", Technip, Paris. Ces méthodes présentent certaines limites. En particulier, elles ne permettent pas de décrire chaque groupe séparément, les variables synthétiques ne peuvent être reliées facilement aux différents groupes, et, par conséquent, être physiquement interprétées. De plus, seule la part d'inertie globale qu'elles représentent peut être connue, et non la part 20 d'inertie qu'elles représentent pour chaque groupe. Il n'est par conséquent pas possible de filtrer chaque groupe indépendamment. Une technique d'analyse également connue dite analyse en composantes principales généralisée (ACPG), permet de comparer différents groupes de variables tout en décrivant chacun d'eux. On la trouve mise en oeuvre sur des données à caractère économique par 25 exemple dans: - Casin, Ph, 2001, "A generalization of principal component analysis to K sets of variables", Computational Statistics & Data Analysis, 35, 417-428. 2844362 La méthode selon l'invention La méthode selon l'invention propose une approche o en exploitant une technique statistique multivariable appliquée à plusieurs jeux de traces d'enregistrement sismiques 3D disponibles, obtenus par exploration sismique d'une zone souterraine, on compacte et on filtre les événements sismiques lisibles sur ces jeux de traces. La méthode comporte - la répartition de ces événements sismiques dans un certain nombre p de familles de variables sismiques ayant chacune une signification physique particulière; - la formation par combinaison des dites variables sismiques, de variables synthétiques en nombre très inférieur au nombre des dites variables, que l'on obtient par 10 construction d'une base vectorielle orthogonale dans chacun des p ensembles d'analyse constitués par les données de chacune des dites familles, d'o découle la formation d'une base vectorielle orthonormée décrivant les p ensembles d'analyse; - l'utilisation de cette base vectorielle orthonormée pour filtrer et décrire les dits événements sismiques. La formation des variables synthétiques est obtenue par exemple en combinant une technique d'analyse en composantes principales ou ses techniques dérivées basées sur la diagonalisation d'une norme (recherche de vecteurs propres) et une technique de régression simple. Suivant un mode de mise en oeuvre, on peut soumettre les variables synthétiques à 20 une rotation pour améliorer les corrélations entre variables initiales et variables synthétiques. Les variables synthétiques peuvent être formées à partir de familles comportant par exemple un ou plusieurs cubes de traces sismiques obtenues à l'issue d'une campagne sismique 3D, chacun correspondant à un même déport ou offset. On peut également former des variables synthétiques à partir de familles constituées de cubes de traces sismiques obtenus à l'issue d'une campagne sismique 3D, chacun correspondant à un même angle d'incidence. On peut encore former des variables synthétiques à partir de familles constituées de cubes d'attributs sismiques formés par traitement des amplitudes de traces sismiques obtenues à l'issue d'une campagne sismique 3D. On peut encore former des variables synthétiques à partir de familles constituées de 5 cubes de traces sismiques obtenues à l'issue de campagnes sismiques successives de la zone. On peut encore former des variables synthétiques à partir de familles constituées de cubes d'attributs sismiques formés par traitement des amplitudes de traces sismiques obtenues à l'issue de campagnes sismiques successives de la zone. La méthode est particulièrement avantageuse en ce qu'elle permet, appliquée à de l'information sismique: - de définir de nouveaux attributs sismiques, combinaisons linéaires des attributs sismiques initiaux, résumant l'information sismique initiale et permettant de filtrer cette information, tout en respectant la structure multicubes des données - de montrer si tous les cubes d'attributs sismiques initiaux, ou seulement certains d'entre eux, apportent une information redondante, et de mesurer la part de cette information redondante - d'apporter une interprétation physique aux nouveaux attributs sismiques. Les nouveaux attributs sismiques ainsi formés peuvent ensuite être analysés afin 20 d'extraire les propriétés géologiques de l'information sismique. Deux types de méthode peuvent classiquement être utilisés à cette fin: une première qui utilise des techniques d'estimation statistique et fournit une interprétation quantitative en termes de propriétés géologiques, et une deuxième qui utilise des techniques de reconnaissance de forme pour analyser les faciès sismiques et fournir une interprétation qualitatitive. Des exemples de ces deux types de méthode sont mis en oeuvre respectivement dans les brevets EP-A-671 017 (US 5 638 269) et FR-A-2 768 818 (US 6051651) du demandeur, déjà cités. Présentation succincte des figures - les figures IA à 1C montrent l'organisation des enregistrements sismiques en p groupes: chaque groupe rassemble dans un cube l'ensemble des portions de traces correspondant à un critère physique identique; - les figures 2A à 2C montrent respectivement les cartes de trois premiers nouveaux attributs sismiques Z(l), Z(2) et Z(3); - les figures 3A à 3D montrent respectivement le premier plan temporel vrai (3A) de la fenêtre d'analyse et le même plan reconstitué/filtré avec 1 (3B), 3 (3C)ou 8 (3D) des nouveaux attributs sismiques pour le cube d'impédance P; - les figures 4A à 4D montrent respectivement le premier plan temporel vrai (4A) de la fenêtre d'analyse et le même plan reconstitué/filtré avec 1 (4B), 3 (4C)ou 8 (4D) des nouveaux attributs sismiques pour le cube d'impédance S; - les figures 5A à 5D montrent respectivement le premier plan temporel vrai (5A) et le même plan reconstitué/filtré avec 1 (5B), 3 (5C)ou 8 (5D) des nouveaux attributs 15 sismiques pour le cube de masse spécifique p; - les figures 6 montrent respectivement les sections verticales, passant par le puits W2, vraie (6A), reconstituée/filtrée avec les 8 premiers nouveaux attributs (6B), ainsi que les différences (résidus) (6C) entre ces deux sections pour le cube d'impédance P; - -les figures 7 montrent respectivement les sections verticales, passant par le puits W2, 20 vraie (7A), reconstituée/filtrée avec les 8 premiers nouveaux attributs (7B), ainsi que les différences (résidus) (7C) entre ces deux sections pour le cube d'impédance S; - la figure 8 montre la localisation des traces d'apprentissage utilisées pour calibrer la fonction de classement dans le cadre de la reconnaissance des formes supervisée des 8 premiers nouveaux attributs sismiques; - les figures 9A, 9B montrent respectivement les attributions des événements sismiques et les attributions les plus fiables de ces mêmes événements dans le cadre de la reconnaissance des formes supervisée des 8 premiers nouveaux attributs sismiques; et - les figures 10A, 10B représentent respectivement les attributions des événements sismiques et les attributions de ces mêmes événements sismiques dans le cadre de la reconnaissance des formes non supervisée des 8 premiers nouveaux attributs sismiques. Description détaillée On dispose de plusieurs jeux d'enregistrements sismiques obtenus à des instants différents par exploration sismique d'une zone souterraine étudiée. Les événements sismiques que l'on cherche à compacter/filtrer sont caractérisés par des attributs sismiques souvent organisés en groupes ou familles à signification physique particulière (figures lA à 1C). Ces attributs peuvent être la succession des valeurs des amplitudes de portions de 10 traces sismiques ou de paramètres élastiques issus d'une inversion élastique jointe avant sommation, obtenue par exemple à l'aide de la méthode décrite dans le brevet FR 2 800 473 (EP 1 096 270)déjà cité le long de la fenêtre temporelle étudiée. - 1- Notations Soit X une matrice composée de p tableaux juxtaposés Xi,... Xj... XP: X = {Xl... Xi... Xp}, Xj = [Xj,... Xj,k... Xj,mj] o Xj est un tableau composé de mj variables Xj,k de n individus (au sens statistique du terme). La matrice X rassemble tous les événements sismiques disponibles. Ces événements sont tout d'abord organisés en p groupes ou familles Xj (figure lA). Chaque famille Xj correspond à une propriété physique particulière. Elle se présente sous la forme d'un cube de valeurs correspondant à cette propriété, repérées par leur position spatiale 20 (ligne, cdp, temps/profondeur), (figures lB et 1C). Ces cubes peuvent être formés par exemple par un ensemble des valeurs d'amplitude de n traces ou portions de traces associées à un déport commun ou à un angle commun ou encore la distribution spatiale des valeurs d'un attribut commun dérivé de ces amplitudes telle que, par exemple, l'impédance acoustique en ondes P ou S, la masse spécifique ou des paramètres de Lamé, dans une 25 fenêtre temporelle choisie de taille mj. Dans cet exemple, chaque variable Xj,k rassemble les valeurs de l'attribut sismique considéré (amplitude sismique ou propriété élastique) pour un échantillon temporel à l'intérieur de la fenêtre d'analyse. Dans cet exemple, un cube comprend autant d'individus que de traces disponibles. Les mj variables Xj,k d'une table Xj définissent un espace dénoté Sj et sont supposées centrées. Soit D, la matrice diagonale des poids des individus. Soit 1-I un opérateur de projection orthogonale sur l'espace Sj défini par: FIj = Xj(XjDnXj) XjDn 2- Calcul des nouveaux attributs Le but de l'Analyse en Composantes Principales Généralisée (ACPG) est à la fois de décrire les proximités entre les p familles de variables et de décrire séparément chaque groupe grâce à des variables synthétiques qui sont des combinaisons linéaires des variables initiales. Alors, ce qui est commun à tous les groupes ou familles, ou à certains d'entre eux 10 seulement, et par conséquent ce qui est différent, peut être mis en évidence. L'ACPG permet de réduire le nombre de variables utilisées pour décrire un jeu de données et permet un filtrage multivarié tout en conservant la structure multicubes des données initiales. Afin d'atteindre ce but, la première variable synthétique calculée, que l'on nomme Z(l), calculée doit être telle que: 1. les p projections orthogonales Zj(l) de Z(l) sur les espaces Sj doivent être le plus proche les unes des autres; 2. Dans chaque groupe j la dispersion des n individus projetés sur Zj() est maximum. p 1. max IR2 (Z(1) Z(l)h j=l mj 2. max CoV2 (Z(l),Xj,) j v=l Ceci se traduit mathématiquement par les conditions suivantes: o R2 est le coefficient de corrélation multiple et Zj(l) est la projection de Z(1) sur l'espace 20 Sj normée à 1. Z(1) = rI(1)z() Ensuite, on cherche une nouvelle variable synthétique Z(2) selon les mêmes critères. On ajoute une condition d'orthogonalité entre les projections les variables Z(k) du même groupe. Le processus s'arrête quand une base orthogonale dans chaque espace Sj est trouvée. Le nombre de variables Z(k) n'excède pas la dimension la plus grande des espaces Sj. Des méthodes de rotation bien connues des gens de l'art, peuvent être mises en oeuvre dans chacun des espaces ainsi trouvés, afin d'optimiser les coefficients de corrélations entre les nouvelles variables et les variables initiales. La première variable z(1) est vecteur propre normé de la matrice de covariance 10 XXtDn correspondant à la plus grande valeur propre. Z(1) est ensuite projeté sur les p espaces Sj et, dans chaque espace, la régression des mj variables par Zj(1) est calculée. La matrice des p résidus juxtaposés Xj(2) définit le nouvel espace d'étude. Ce nouvel espace correspond à l'espace initial duquel l'influence de Zj(l) a été retiré dans chaque groupe Sj. Les mj variables de la table Xj(2) définissent un sous espace Sf2) de Sj. La seconde 15 composante Z(2) est le vecteur propre normé de la matrice de covariance X(2)X(2)tDD des résidus correspondant à la valeur propre la plus grande. Z(2) est projetée sur les p espaces Sj(2). Et la régression des mj variables par Zj(2) est calculée dans chaque espace Sf(2). La régression effectuée dans chaque espace Sj assure que la projection de Z(2) sur l'espace Sj(2 est orthogonale à la projection de Z(l) sur l'espace Sj. Le calcul des vecteurs propres est effectué par toute méthode connue des gens de l'art. En pratique, la matrice XXtD. étant d'ordre n, n pouvant être très grand, à une étape k, on calcule le vecteur propre normé V(k) de la matrice X(k)tD 1X(k) associé à la plus grande valeur propre. Le vecteur Z(k) est ensuite simplement calculé par transformation linéaire 25 Z(k)=X(k)V(k) et normé. Chaque groupe peut être reconstitué avec un petit nombre d'attributs. Les variables initiales du groupe j peuvent s'exprimer comme la somme d'une combinaison linéaire des variables Zj(k) et d'un résidu. Si toutes les variables Zj(k) d'un groupe sont utilisées, le résidu est nul. Les coefficients de la combinaison linéaire ne sont autres que les 30 coefficients calculés à l'issu de chaque régression dans les espaces Sj. 1l La méthode selon invention a pour propriété, entre autres, que les variables Z(k) sont orthogonales deux à deux. En vue de l'interprétation des variables synthétiques calculées, on calcule ensuite pour chaque groupe j: 1. le carré de la corrélation entre la variable Z(k) et la variable Zj(k) du groupe j, indiquant le degré de proximité entre ces deux variables; 2. le pourcentage d'inertie du groupe j expliqué par la variable Zj(k). Les corrélations entre les variables Z(k) et les variables initiales Xj,k sont également calculées pour mettre en évidence d'éventuelles relations entre les Xj,k. Les corrélations 10 entre les variables Zj(k) et Xj,k d'un même groupe j. ou les corrélations entre variables Zj(k) pour un même k peuvent également être calculées. Les nouvelles variables synthétiques ainsi créées constituent les nouveaux attributs sismiques recherchés. 3 - Exemple d'application La méthode selon invention a été appliquée à trois cubes de propriétés élastiques (impédance P notée Ip, impédance S notée Is et masse spécifique notée p) obtenus après inversion élastique jointe avant sommation à l'aide de la méthode décrite dans le brevet FR 2 800 473 déjà cité d'enregistrements sismiques regroupés en cubes à angles d'incidence commun. Chaque propriété élastique engendre un groupe de variables ou 20 d'attributs sismiques. La fenêtre temporelle d'analyse se compose de 19 échantillons de signaux prélevés sur les enregistrements sismiques. Nous considérons comme variables les valeurs de la propriété élastique pour un échantillon temporel à l'intérieur de la fenêtre d'analyse. Chaque groupe est alors composé de 19 variables composées d'autant d'individus que de traces disponibles. Comme il y a trois groupes de variables, nous 25 disposons au total de 57 variables. Les trois cubes de paramètres ont été préalablement horizontalisés par rapport à une surface de référence. Les variables synthétiques Z(k) obtenues résument les informations respectives des trois propriétés élastiques, tout en prenant en compte les relations entre elles. 2844362 12 3.1 - Analyse des relations entre les impédances P, S et la masse spécifique p Le tableau 1 par les projections groupes. Il montre 5 l'impédance P, 94,4 ci-dessous donne la proportion de variance en pourcentage expliquée des 8 premières variables synthétiques Z(> à Z8) dans les trois qu'il suffit de 8 variables pour résumer 87,4% de la variance de % de la variance de l'impédance S et 88,4% de la variance de la masse spécifique. Z(L) zPj) Z[) zffl ZM) a Ip 23,4 24,5 6,3 13,1 9,8 3,1 3,8 3,4 87,4 Is 26,6 14,2 15, 2 12,6 5,7 7,2 5,8 7,1 94,4 16,9 14,1 13,6 10,2 8,0 9,4 10,1 6,1 88,4 Le tableau 2 ci-après donne le carré de la corrélation entre les variables synthétiques Z(k) et leurs projections sur les trois groupes. Cette valeur indique le degré de 10 proximité entre la variable et le groupe. Plus la valeur est proche de 1, plus la variable est proche de sa projection et donc du groupe. z(I) Z(2) Z(4) Z(È) Z(7) Ip 0,48 0,79 0,38 0,72 0,71 0,24 0,02 0,30 Is 0, 91 0,43 0,86 0,53 0,36 0,59 0,09 0,73. 0,78 0,49 0,77 0,46 0,65 0,79 0,96 0,46 Par exemple, Z(1) est très proche de Is et proche de p, mais très éloigné de Ip. Alors que Z(2) est proche de Ip, mais très éloigné de Is et de p. Sur les 8 premières variables Z(k), trois sont proches de Ip seul, deux sont proches à la fois de Is et de p, une est proche de Is seul, et deux sont proches de p seul. Ainsi, il apparaît que Ip et Is ne 5 contiennent pas la même information et qu'une part de l'information contenue dans Is et dans p semble commune. 3.2 - Filtrage des impédances P, S et de la masse spécifique Les trois propriétés élastiques peuvent être reconstituées avec un nombre restreint de nouveaux attributs sismiques. Le nombre de variables synthétiques utilisées pour la 10 reconstitution peut être différent d'un groupe à l'autre: il dépend du pourcentage d'inertie que l'on désire reconstituer/filtrer. 3.3 - Reconnaissance des formes supervisée et non supervisée des nouveaux attributs Il est possible d'interpréter géologiquement les nouveaux attributs sismiques au 15 travers d'une analyse de faciès avec un algorithme de reconnaissance de forme supervisé ou non. En reconnaissance des formes supervisée, il est nécessaire de choisir des classes d'apprentissage regroupant des événements sismiques situés aux voisinages de puits dont l'interprétation géologique a été effectuée. Cette méthodologie permet de classer 20 l'ensemble des événements sismiques en fonction des caractères géologiques observés aux puits. Dans notre cas, les connaissances géologiques sur les trois puits disponibles Wl, W2 et W3, ainsi que les connaissances géologiques générales sur la zone d'étude nous amènent à considérer trois classes: une première classe correspondant aux puits W1 et W2 25 représentent des sables turbiditiques, une seconde classe correspondant au puits W3 représente des boues turbiditiques et une troisième classe représente des hémipélagites. Les puits Wl et W2 montrant des caractéristiques géologiques proches, seuls les événements sismiques situés autour du puits W2 sont utilisés, le puits Wl servant de test aveugle. La localisation du puits W3 étant incertaine, les événements sismiques utilisés pour 30 caractériser la classe 2 sont situés autour d'un puits virtuel défini à proximité de W3. Enfin, ne disposant d'aucun puits pour la classe 3, un puits virtuel est défini à partir de nos connaissances géologiques générales. La Fig.7 montre la localisation géographique des trois classes d'apprentissage. Chaque classe est constituée de 121 événements sismiques (llxll). Ces événements servent à calibrer une fonction de classement qui est ensuite s appliquée à l'ensemble des événements sismiques. Les Fig.8A, 8B montrent respectivement les affectations obtenues et les affectations les plus probables, c'est-à-dire ayant une probabilité d'affectation supérieure ou égale à 80%. Le puits Wl, utilisé en test aveugle, est correctement attribué à la même classe que le puits W2. Le puits W3 est localisé en bordure des classes 1 et 2. La classe 2 se distribue le 10 long d'une forme chenalisante. La carte de la figure 8B montre que les résultats sont stables: peu d'événements présentent une mauvaise probabilité d'affectation. En reconnaissance des formes non supervisée, les classes d'apprentissage sont constituées en utilisant les modes d'une fonction de densité de probabilité multivariée calculée à partir de l'ensemble des événements sismiques, dans un espace engendré par les 15 huit premiers nouveauxattributs. Les événements sismiques de plus forte densité de probabilité constituant les pics de cette fonction de densité sont sélectionnés pour réaliser l'apprentissage de la fonction de classement. Cette fonction est ensuite appliquée à l'ensemble des événements sismiques considérés pour les affecter à l'une des classes: les Fig.9A, 9B montrent les affectations obtenues, ainsi que les affectations les plus probables, 20 c'est-à-dire ayant une probabilité d'affectation supérieure ou égale à 80%. Les cartes obtenues sont très proches des cartes des fig. 8A, 8B. La classe 3 de l'analyse supervisée correspond aux classes 3, 4 et 5 de l'analyse non supervisée. Les puits Wl et W2 sont attribués à la même classe et le puits W3 est toujours situé à la frontière entre les classes 3, et surtout 1 et 2. Cette convergence des résultats valide l'interprétation de l'approche 25 supervisée. REVENDICATIONS 1) Méthode pour compacter et filtrer des événements sismiques 3D lisibles sur plusieurs jeux de traces d'enregistrement sismiques disponibles, obtenus par exploration sismique d'une zone souterraine, caractérisée en ce qu'elle comporte: s - la répartition de ces événements sismiques dans un certain nombre p de familles de variables sismiques ayant chacune une signification physique particulière; - la formation par combinaison des dites variables sismiques, de variables synthétiques en nombre très inférieur au nombre des dites variables, que l'on obtient par construction d'une base vectorielle orthogonale dans chacun des p ensembles d'analyse 10 constituées par les données de chacune des dites familles, d'o découle la formation d'une base vectorielle orthonormée décrivant les p ensembles d'analyse; et - l'utilisation de cette base vectorielle orthonormée pour filtrer et décrire les dits événements sismiques. 2) Méthode selon la revendication 1, caractérisé en ce que la formation des 15 variables synthétiques est obtenue en combinant une technique d'analyse en composantes principales ou des méthodes dérivées de diagonalisation d'une norme et une technique de régression simple. 3) Méthode selon la revendication 1 et 2, caractérisée en ce que l'on soumet à une rotation les variables synthétiques pour améliorer les corrélations entre variables initiales et 20 variables synthétiques. 4) Méthode selon la revendication 1 ou 2 ou 3, caractérisée en ce que l'on forme des variables synthétiques à partir de familles constituées de cubes de traces sismiques obtenus à l'issue d'une campagne sismique 3D, chacun correspondant à un même déport ou offset. ) Méthode selon la revendication 1 ou 2 ou 3, caractérisée en ce que l'on forme des 25 variables synthétiques à partir de familles constituées de cubes de traces sismiques obtenus à l'issue d'une campagne sismique 3D, chacun correspondant à un même angle d'incidence. 6) Méthode selon la revendication 1 ou 2 ou 3, caractérisée en ce que l'on forme des variables synthétiques à partir de familles constituées de cubes d'attributs sismiques obtenus par traitement des amplitudes de traces sismiques obtenues à l'issue d'une campagne sismique 3D. 7) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on forme des variables synthétiques à partir de familles constituées de cubes de traces sismiques obtenus par des explorations sismiques successives de la zone. 8) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on forme des variables synthétiques à partir de familles constituées de cubes d'attributs 10 sismiques formés par traitement des amplitudes de traces sismiques obtenues par des explorations sismiques successives de la zone.