FR 2773618 A1 19990716 FR 9900797 A 19990114 Atténuation de réflexions multiples de données du fond marin à composantes multiples La présente invention concerne des méthodes de réduction ou d'atténuation de réflexions multiples dans des données du fond marin à composantes multiples. ARRIÈRE-PLAN DE L'INVENTION Les levés sismiques marins sont habituellement réalisés par remorquage d'une source d'énergie et de détecteurs sismiques derrière un navire. La source communique une onde acoustique à l'eau, en créant un champ d'ondes qui se propage de façon cohérente dans le substratum. Lorsque le champ d'ondes frappe les interfaces entre les formations terrestres, ou strates, il est réfléchi à travers la terre et l'eau jusqu'aux détecteurs, o il est converti en signaux électriques et enregistré. L'analyse de ces signaux permet de déterminer la forme, la position et la lithologie des formations sous-jacentes. Dans d'autres méthodes de levés sismiques, les détecteurs et/ou sources sont placés au niveau ou à proximité du fond marin ou dans des puits. Un problème que l'on recontre lors de levés sismiques - et lors de la réalisation de profils sismiques verticaux inverses - est celui de la réverbération dans une colonne d'eau. Le problème, qui résulte du pouvoir de réflexion intrinsèque de la surface et du fond de l'eau, peut être expliqué comme suit. Une onde sismique générée dans (réfléchie par) les strates terrestres passe dans l'eau dans une direction généralement ascendante. Cette onde, que l'on appelle la "primaire", se propage dans l'eau et passe par le détecteur sismique qui enregistre sa présence. Le champ d'ondes continue sa progression vers la surface de l'eau, o il est réfléchi vers le bas. Ce champ d'ondes réfléchi ou "fantôme" se propage également dans l'eau et passe par le ou les détecteur(s) o il est à nouveau enregistré. Selon la nature des matériaux terrestres au fond de l'eau, le champ d'ondes fantôme peut être lui-même réfléchi vers le haut dans l'eau, en donnant ainsi naissance à une série d'une ou de plusieurs réflexions fantômes ou réflexions multiples ultérieures. Les réflexions multiples de surface libre peuvent être classées suivant leur ordre, qui est égal au nombre de réflexions depuis la surface libre. Ainsi, -2- les réflexions de surface libre du premier ordre comprennent de l'énergie se propogeant initialement vers le bas en partant des sources (par opposition aux réflexions fantômes dans lesquelles l'énergie se propage vers le haut et est réfléchie depuis la surface libre), réfléchie vers le haut depuis le fond marin ou une limite sous le fond marin, puis réfléchie vers le bas depuis la surface libre jusqu'aux hydrophones. Les réflexions multiples de surface libre du deuxième ordre font l'objet de deux réflexions vers le bas depuis la surface de l'océan avant d'être détectées par les hydrophones, et ainsi de suite. Cette réverbération du champ d'ondes sismique dans l'eau obscurcit les données sismiques, en amplifiant certaines fréquences et en en atténuant d'autres, ce qui rend difficile l'analyse des formations terrestres sous- jacentes. Dans les cas o les matériaux terrestres au fond de l'eau sont particulièrement durs, un excès d'énergie acoustique ou de bruit généré par la source sismique peut également se trouver piégé dans la colonne d'eau, en se réverbérant de la même manière que les ondes sismiques réfléchies elles-mêmes. Ce bruit présente souvent une grande amplitude et a de ce fait tendance à couvrir les signaux plus faibles des réflexions sismiques que l'on recherche pour l'étude. Dans l'art, Ruehie divulgue, dans le brevet américain n 4 486 865, une technique destinée à réduire les réflexions fantômes, dans laquelle un capteur de pression et un détecteur de la vitesse acoustique sont positionnes dans l'eau, à proximité immédiate l'un de l'autre. La sortie de l'un au moins des détecteurs est réglée en gain et filtrée, en utilisant une opération de déconvolution présentant une quantité prédéterminée de bruit blanc par rapport au décalage nul de la fonction d'autocorrélation. Le brevet suggère qu'en ajoutant ce signal déconvolué/réglé en gain à la sortie de l'autre détecteur, il est possible d'annuler les réflexions fantômes. Dragoset décrit, dans le brevet américain n 5 365 492, une méthode itérative de détermination d'un facteur d'échelle utilisé pour sommer des mesures simultanées de pression et de vitesse. Haggerty divulgue, dans le brevet américain n 2 757 356, un système de levé sismique marin par réflexion dans lequel deux dispositifs -3sismographiques sont disposés à deux profondeurs distinctes dans l'eau de façon à ce que les réverbérations dans la colonne d'eau qu'ils reçoivent soient déphasées de 180 degrés. En combinant la sortie des détecteurs, le brevet suggère que les réverbérations s'annuleront. White, dans 'Seismic Wave Radiation - Transmission and Attenuation', McGraw-Hill, 1965, puis Barr, dans le brevet américain n 4 979 150, proposent tous deux l'utilisation d'une paire hydrophone/géophone pour séparer les ondes se propageant vers le haut et vers le bas. Tous deux utilisent une décomposition en ondes planes de l'équation d'onde afin de déterminer un facteur d'échelle permettant d'additionner des mesures issues des différents types de détecteurs. La solution proposée fait toutefois l'hypothèse d'une incidence normale des ondes. Ikelle et coll. utilisent, dans la demande de brevet britannique GB-A- 9710435, une technique par série de Born en vue d'éliminer des réflexions multiples des levés du fond marin. Berni propose, dans le brevet américain n 4 345 473, I'utilisation d'un accéléromètre à composante verticale combiné à un hydrophone pour annuler le bruit réfléchi par la surface dans les opérations sismiques marines. Gal'perin propose, dans 'Vertical Seismic Profiling' Publication Spéciale n 12 de la Society of Exploration Geophysicists, I'utilisation d'un détecteur sismique qui combine la sortie d'un capteur de pression et de vitesse en vue d'une utilisation dans des levés en profils sismiques verticaux. Moldoveanu utilise, dans le brevet américain n 5 621 700, un câble de fond marin sensible à la pression et à la vitesse. Les réverbérations sont atténuées par un processus de moyennage mettant en jeu les étapes consistant à ajouter le produit des données de pression par la valeur absolue des données de vitesse et.le produit des données de vitesse par la valeur absolue des données de pression. L. Amundsen et A. Reitan ont publié dans Geophysics, vol. 60, 2, 1995, pages 563-572, une méthode de décomposition de données du fond marin à composantes multiples utilisant des filtres de décomposition déterminés par analyse en ondes planes. À- 4 - K.M. Schalkwijk et coll. ont publié dans SEG, Expanded Abstracts, 1997, pages 8-11, une méthode de décomposition simultanée du champ d'ondes enregistré en ondes ascendantes et descendantes et de compression (P) et de cisaillement (S) par une décomposition en une étape. La méthode décrite n'a pas donné de résultats satisfaisants lorsqu'on l'a appliquée à des données réelles. A la lumière des références susmentionnées, on considère qu'un objet de la présente invention est de procurer des méthodes d'élimination des réflexions fantômes ou multiples de levés marins, en particulier de mesures sur le fond marin. Un objet spécifique de l'invention est de procurer une telle méthode sans faire aucune hypothèse quant à la direction d'incidence des ondes acoustiques. RESUME DE L'INVENTION Conformément à la présente invention, on procure une méthode de réduction de réflexions multiples dans des données enregistrées au cours d'un levé sismique marin, comprenant les étapes consistant à (a) fournir des données de champ d'ondes sismique enregistrées par des récepteurs sismiques de composantes multiples, lesdites données comprenant des mesures liées à des composantes de pression et de vitesse dudit champ d'ondes; (b) choisir une composante du champ d'ondes; (c) sélectionner deux composantes ou plus des données de champ d'ondes fournies, nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (d) filtrer l'une ou plusieurs des composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (e) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à réflexions multiples réduites de ladite composante choisie. Contrairement aux tentatives connues de mise à l'échelle des différentes composantes du champ d'ondes préalablement à la sommation, telles que décrites par exemple dans les brevets américains no 4 486 865 et 5 365 492, le filtre de la présente invention est un filtre spatial utilisant une opération de convolution qui combine plusieurs traces pour générer une trace en sortie. Une trace est la mesure d'une composante du champ - 5 - d'ondes en une position. Dans son mode de réalisation le plus général, le filtre varie également en fréquence. Il est donc possible de considérer des contributions du champ d'ondes avec un angle d'incidence non nul, de manière unifiée, sans avoir à calculer séparément de facteurs d'échelle dépendant de l'angle. Relativement à la méthode décrite par K. M. Schalkwijk et coll., dans SEG, Expanded Abstracts, 1997, pages 8-1 1, la présente invention sépare I'étape consistant à supprimer des réflexions multiples descendantes issues de la décomposition en ondes de compression (P) et de cisaillement (S). Cette technique conduit à de meilleurs résultats lorsqu'on l'applique à des données réelles. Le filtre de suppression des ondes de réflexion dérive de préférence d'une solution de l'équation d'onde él&stodynamique grâce au théorème de représentation. Le filtre de suppression des ondes de réflexion tente de reconstruire le champ d'ondes descendant à partir de mesures de composantes multiples. En soustrayant le champ d'ondes descendant reconstruit du champ d'ondes total mesuré, on peut en déduire le champ d'ondes ascendant. Le filtre est de préférence un filtre spatial utilisant une opération de convolution qui combine plusieurs traces (mesures de récepteurs uniques) pour générer une trace en sortie. Dans son mode de réalisation le plus général, le filtre varie également en fréquence. Puisque le filtre de suppression de réflexions multiples dérive du théorème de représentation plutôt que d'une décomposition en ondes planes, le fond marin peut être arbitrairement irrégulier. Le filtre élimine les réverbérations de l'eau à l'image des méthodes décrites dans l'art antérieur. Il atténue toutefois également d'autres réflexions multiples qui ne sont pas supprimées lorsqu'on utilise une méthode connue. Par ailleurs, la nouvelle méthode englobe naturellement les signaux à décalage non nul sans qu'on ait à introduire de facteur directionnel séparé. En pratique, les mesures qui combinent des signaux issus de deux récepteurs différents, comme des géophones et des hydrophones, nécessitent un calibrage supplémentaire en vue d'une adaptation de leur 6 - sortie et, en particulier dans le cas de mesures au fond de l'océan, en vue d'adapter des différences de couplage du récepteur et du fond marin. Les hommes de l'art apprécieront et comprendront ces caractéristiques de I'invention, ainsi que d'autres, des modes de réalisations préférés et des variantes de celle-ci, et des avantages supplémentaires de l'invention grâce à la description détaillée et aux dessins ci-dessous. BREVE DESCRIPTION DES DESSINS La figure 1 illustre schématiquement des événements sismiques marins tels qu'enregistrés dans un câble du fond marin; les figures 2A, B illustrent schématiquement un champ d'ondes ascendant et descendant tel qu'enregistré dans un câble du fond marin; la figure 3 représente graphiquement les conditions aux limites pour un levé sismique marin utilisant un câble de fond marin; la figure 4 illustre un organigramme représentant des étapes élémentaires d'une méthode conforme à l'invention. MODE(S) DE MISE EN OEUVRE DE L'INVENTION Si l'on se réfère à la figure 1, il y est illustré des événements sismiques typiques du fond marin. Un cercle plein noir désigne une source. Un triangle désigne un récepteur. Dans les exemples, les sources sont situées à proximité de la surface 10 de la mer. Les récepteurs sont situés sur le fond marin 11. Sous le fond marin, il est illustré un réflecteur unique 12 à une profondeur arbitraire. Des trajets suivis par les signaux sont illustrés pour divers événements sismiques ayant occasionné des réflexions multiples dans les données enregistrées. En fonction de la direction d'incidence au niveau du site du récepteur, l'événement peut être séparé en événements ascendants (figure 2A) et descendants (2B). Comme illustré par la figure 2A, le champ d'ondes descendant contient essentiellement des réflexions multiples, à savoir les réverbérations 21 dans la couche d'eau et les réflexions multiples 22 aux jalons des récepteurs. Il comporte également la primaire 20 du fond marin. Le champ 7 - d'ondes ascendant, illustré à la figure 2B, contient toutes les primaires 23 à l'exception de la primaire du fond marin et de certains événements de réflexions multiples, en particulier les réflexions multiples 24 aux jalons des sources qui ne sont pas incluses dans le champ d'ondes descendant. La nouvelle méthode de suppression de réflexions multiples telle que décrite ci-dessous est basée sur le théorème de représentation élastodynamique. Soit vj la j-ième composante de la vitesse acoustique et p la pression. Les champs correspondants débarrassés de réflexions multiples, juste en dessous du fond marin, sont désignés par un indice supérieur - (tilde). On peut montrer que la valeur de v; débarrassée de réflexions multiples peut s'écrire sous la forme: Vj((o, X) = Vj(O, X) -- [1] J[ dS (x') [ (iop (co, x') G3j (x', x, ko, k) - Vi(o), X') lI "lij (X', X, ka, kg)] en un point x dans un domaine D donné limité par une surface S, et o Gij est le tenseur de Green du champde déplacement élastodynamique, Epqk est le tenseur des tensions de Green et kox = Moa et kg = --g sont les nombres d'onde des ondes P et S, respectivement, et co est la fréquence circulaire. En utilisant la loi de Hooke, on peut écrire une mesure de pression débarrassée de réflexions multiples sous la forme: [2] (), x) = i((3 p(x))-1{ (02(x) - 2 2(x)) [ 2] 2] [aiVi(o, x) -+ 2v2(o, x)] + oc (x) a3v3(', x)} o les vélocités v débarrassées de réflexions multiples sont données par l'équation [1], et le di désigne des dérivées partielles, p(x) donnant la densité au point x. Il est important de noter que les équations [1] et [2] sont des équations générales de suppression de réflexions multiples pour des données de fond océanique, et elles sont valides pour un fond marin de forme arbitraire et des paramètres de milieu à variations latérales. Pour évaluer l'équation, on peut utiliser le traçage dynamique de rayons pour déduire -8- une expression numérique des éléments du tenseur de Green. Le profil S du fond marin est une limite absorbante pour le tenseur de Green. On peut déduire une solution analytique des équations [1] et [2] en faisant des hypothèses simplificatrices. Dans l'exemple ci-dessous, on suppose par exemple que les paramètres du fond marin sont constants. Dans le cadre de l'évaluation, la notation utilisée dans les équations [1] et [2] est modifiée. Les points sur la surface S sont à présent désignés par, ' et ceux dans le volume D sont désignés par x, x'. Les limites utilisées pour évaluer l'équation intégrale [1] sont illustrées à la figure 3, le domaine D étant limité par le fond marin 31 et un autre niveau de référence 32 juste en dessous du fond marin. Les limites de droite et de gauche du domaine D s'étendent à l'infini. En amenant le niveau de référence très près du fond marin, on peut réduire le calcul du champ d'ondes descendant à l'évaluation de l'intégrale de représentation élastodynamique le long du fond marin. En supposant également que les paramètres du milieu juste en dessous du fond marin sont constants, on peut utiliser les tenseurs de Green dans l'espace libre pour Gij(x, o), x') et Ipqk(x, co, x'). Grâce à ces approximations, on peut évaluer analytiquement l'intégrale de représentation élastodynamique. Soient {p,vx,vy,vz} les données sismiques à quatre composantes enregistrées directement au niveau du fond marin; p désigne le champ de pression, vx et vy désignent les composantes horizontales de la vitesse acoustique, et vz la composante verticale de la vitesse acoustique. Conformément à la décomposition haut et bas standard en un lieu particulier, par exemple le fond marin, {p,vx,vy,Vz} est décomposé en un champ d'ondes descendant {pD,vxD,vyDvzD} et en un champ d'ondes ascendant {pUvxuvyUvzu}l de sorte que [3] {P,Vx, vy, vz} = {pDVD, VY, V zD}+ {pU VXvvz} En utilisant le théorème de représentation comme décrit précédemment, on trouve que le champ d'ondes ascendant {pU,vxU,vyUvzU} peut être déterminé en utilisant 99- p4(xS,;,)) = 2 [p(x, t, o)) - FP;vz (xS, t, Cl) * v (x, t, o))} vIU(x4, t;] C)) = 2 [v(x6, c;,), c) * v2(xs, t;, vyU(xs,,), ) - rvY;Vz(xS, wc) *v(x, t, vzu(xs, t; cl) = 2 [v.(xS, t, c) - rvz;vx(x, t, () * v"(x5, t, ) - vJU(x., {, [V X,) :;(S XXI rvz;VY(x t, ") * vY(xS, t; cl) - rvz;P( t, l) * p(x,, t;Co, o Xs =(xs,ys,Zs) est la position des sources et,:=(1,2,) est la position des récepteurs. Le symbole * indique une convolution spatiale par rapport à. Les coefficients F sont connus analytiquement. Par exemple: rP;Vz (t, o) - icop: + 4 + a+) + [51 4 L;) (81 + 22j gp(o, t,)) - iwp [4@5)(a2 + a'] gs(0,;, Co) o 9p(Oi E cl)) =4 [C exp - i, [6] 1 gs(0, t, C) - exp Co I} 4-n Ik- e xp, et o les constantes p, cc et D sont respectivement la densité, la vitesse des ondes P et la vitesse des ondes S juste en dessous du fond marin. Comme les récepteurs sont positionnes sur le fond marin, les dérivées dans [5] sont effectuées le long de la surface du fond marin (la variable de dérivation est). Pour les applications pratiques, le filtre F peut être approximé de façon à limiter sa longueur et donc augmenter la vitesse à laquelle la méthode peut être mise en oeuvre. Parmi les diverses méthodes connues permettant de réduire la longueur d'un filtre tout en préservant son efficacité, on décrit ci-dessous deux méthodes: - 10- En utilisant un développement de Taylor sur le filtre, dans le domaine des nombres d'onde, on peut écrire le filtre sous la forme d'une série infinie dans le nombre d'onde horizontal. Si l'on ne conserve qu'un certain nombre de termes, ceux-ci peuvent subir une transformation de Fourier inverse vers un domaine spatial pour donner un filtre spatial approximatif de suppression de réflexions multiples. Ce filtre s'exprime en fonction de puissances d'opérateurs de dérivée seconde, qui peuvent être ici encore approximés au degré de précision requis. Cette technique équivaut à développer des équations d'ondes unidirectionnelles pour l'extrapolation/la migration de champs d'ondes pour donner les équations bien connues dites à 5, 15, 30, 45 etc. degrés - pour manipuler correctement des ondes se propageant avec un angle allant jusqu'à 5, 15, 30, 45 degrés par rapport à l'axe vertical. En utilisant des filtres optimisés numériquement pour déterminer par optimisation un filtre spatial de longueur prédéfinie dont l'expression en nombres d'onde correspond, d'une certaine façon, au vrai filtre à nombres d'onde sur un intervalle de nombres d'onde prédefini. Ces filtres peuvent être très courts par comparaison avec celui donné par l'équation [4]. Cette technique est similaire à celle utilisée pour optimiser des opérateurs de dérivée dans la dérivation numérique. On peut par ailleurs accroître la précision et la portée de la méthode décrite ci-dessus en utilisant diverses extensions: D'après l'équation [5], les coefficients F sont fonction de la géométrie du fond marin et de paramètres élastiques. Lorsque ces informations ne sont pas disponibles ou ne sont pas suffisamment précises, une optimisation similaire à celle décrite dans la demande de brevet international WO96/20417 pour déterminer ces coefficients numériquement. Du fait que la séparation du champ d'ondes utilisée ici est basée sur le théorème de représentation au lieu d'une décomposition en ondes planes, le fond marin peut être arbitrairement irrégulier dans le calcul du champ d'ondes ascendant décrit dans l'équation [4]. En appliquant une technique par série de diffusion, on peut éliminer davantage de réflexions multiples du champ d'ondes ascendant. - 11 - Bien que l'équation [4] permettant de calculer le champ {pU,vx, vyUvzu} élimine un nombre significatif de réflexions multiples, y compris toutes les réverbérations dans la masse d'eau et toutes les réflexions multiples au niveau des jalons des récepteurs, comme décrit à la figure 2A, on peut étendre le processus: en utilisant une opération de type diffusion inverse, comme décrit par exemple dans la demande de brevet britannique GB-A- 9710435, on peut éliminer les réflexions au niveau des jalons des sources (figure 2B). Les premiers termes d'une série de diffusion pour le champ de pression ascendant sont: [7] DU(x, co) =pU(x, C) + Al (Co)p U(x,, o) + A2 (C..)p2U (X,, +..., ou O O [8] Pnu(xs ' 0) = dxp U(x, 0b) dx'gp(x, (o x')pn1 u(x',, ) Si l'on se réfère à présent à la figure 4, il y est décrit les étapes élémentaires pour appliquer l'atténuation de réflexions multiples décrite ci- dessus à des collectes de tirs. Suivant le champ d'ondes duquel doivent être éliminées des réflexions multiples, on sélectionne un ensemble de deux composantes, ou plus, des champs d'ondes. L'une des combinaisons suivantes doit être sélectionnée: {p, vz}, {vx, vz}, {vy, vz}, {p,vx,vy, vz}. Toute combinaison de la composante horizontale peut être formée avec le champ de pression ou la vitesse verticale. On décrit les principales étapes pour l'application: 1. Collecter des données sismiques pré-superposées à composantes multiples. 2. Sélectionner la composante devant être débarrassée de réflexions multiples et sélectionner les autres composantes nécessaires à la soustraction. Par exemple, pour débarrasser la pression p de réflexions multiples, sélectionner vz. 3. Déterminer le filtre de suppression de réflexions multiples par optimisation ou par utilisation directe de l'équation [5], ou les deux. 4. Effectuer la somme dans [4] pour obtenir la pU' débarrassée de réflexions multiples. -12 - Une extension de la présente invention comprend l'étape consistant à décomposer le champ d'ondes mesuré en ondes de compression (P) et de cisaillement (S) en utilisant, par exemple, la méthode décrite ci-dessus par Amundsen et Reitan. Il est important de noter, pour la portée de l'invention, que le filtre r' des équations [4,5] est déduit dans le domaine x-co. En utilisant des méthodes de transformation standard, le filtre peut être transformé dans d'autres domaines, tels que f-k ou t-p. -13- REVENDICATIONS 1.- Méthode de réduction de réflexions multiples dans des données enregistrées au cours d'un levé sismique marin, comprenant les étapes consistant à (a) fournir des données de champ d'ondes sismique enregistrées par des récepteurs sismiques de composantes multiples, lesdites données comprenant des mesures liées à des composantes de pression et de vitesse dudit champ d'ondes; (b) choisir une composante du champ d'ondes; (c) sélectionner deux composantes ou plus des données de champ d'ondes foumrnies, nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (d) filtrer l'une ou plusieurs desdites composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (e) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à réflexions multiples réduites de ladite composante choisie. 2.- Méthode selon la revendication 1, dans laquelle le filtre comporte des contributions du champ d'ondes avec des angles d'incidence nuls et non nuls à une position de récepteur donnée. 3.- Méthode selon la revendication 1, comprenant en outre l'étape consistant à décomposer le champ d'ondes en ondes P et S en utilisant un filtre séparé. 4.- Méthode selon la revendication 1, dans laquelle le filtre est fonction de la fréquence. 5.- Méthode selon la revendication 1, dans laquelle le filtre est fonction de l'espace et de la fréquence. 6.- Méthode selon la revendication 1, dans laquelle les composantes choisies sont choisies parmi un groupe comprenant la pression p et trois composantes de vitesse orthogonales vx, vy, vz. -14- 7.- Méthode selon la revendication 1, dans laquelle les composantes choisies sont choisies parmi un groupe comprenant la pression p et trois composantes de vitesse orthogonales vx, vy, vz, et les composantes nécessaires à l'atténuation de réflexions multiples dans lesdites composantes définies sont sélectionnées selon les étapes suivantes: si la pression p est choisie, sélectionner alors la pression p et la composante de vitesse verticale vz pour réduire les réflexions multiples; si la composante verticale vz est choisie, sélectionner alors vz, la pression p et les deux composantes de vitesse horizontales vx et vy; et si l'une des deux composantes de vitesse horizontales vx et vy est choisie, sélectionner alors ladite composante horizontale et la composante verticale vz. 8.- Méthode selon la revendication 1, dans laquelle le filtre est basé sur une solution du théorème de représentation élastodynamique. 9.- Méthode selon la revendication 1, dans laquelle le filtre reconstruit le champ d'ondes descendant à partir de mesures à composantes multiples du champ d'ondes total. 10.- Méthode selon la revendication 1, comportant en outre l'étape consistant à utiliser une correction pour réaliser l'adaptation de la sortie de différents types de récepteurs et des différences de couplage avec le fond marin. 11.- Méthode selon la revendication 1, dans laquelle le traitement de la variante à réflexions multiples réduite se poursuit en utilisant des méthodes connues de traitement de données sismiques. 12.- Méthode selon la revendication 1, dans laquelle les données de champ d'ondes sont enregistrées par des récepteurs situés au niveau du fond marin. 13.- Levé sismique marin, comprenant l'étape consistant à poser un support de récepteurs sur le fond marin, ledit support de récepteurs comprenant une pluralité de récepteurs de composantes multiples capables de recevoir au moins deux composantes indépendantes de mesures de vitesse et/ou de pression; générer de l'énergie acoustique pour qu'elle se propage à travers des formations terrestres en dessous -15- dudit support de récepteurs; utiliser lesdits récepteurs pour mesurer des données de champ d'ondes à composantes multiples de ladite énergie acoustique en des positions des récepteurs; et les étapes supplémentaires consistant à: (a) choisir une composante du champ d'ondes; (b) sélectionner deux composantes ou plus des données de champ d'ondes nécessaires à l'atténuation de réflexions multiples dans ladite composante choisie; (c) filtrer l'une ou plusieurs desdites composantes sélectionnées en utilisant un filtre spatial qui combine des composantes de données mesurées à des positions différentes des récepteurs et opère indifféremment sur des ondes P et S; et (d) sommer les composantes sélectionnées et sélectionnées filtrées pour générer une variante à réflexions multiples réduites de ladite composante choisie.