FR 2769718 A1 19990416 FR 9712622 A 19971009 1 2769718 Procédé de séparation d'ondes en sismique de puits pour des acquisitions de type à déport croissant La présente invention concerne une méthode pour la séparation des ondes élastiques reçues et enregistrées après propagation dans un milieu et, plus particulièrement, le traitement des signaux enregistrés lors d'une acquisition effectuée par émission d'ondes élastiques en plusieurs positions d'un ou plusieurs émetteurs disposés en surface du milieu, le ou les récepteurs étant fixe(s) dans un puits foré dans ledit milieu pendant toute la série d'émission, le ou les récepteurs étant ensuite déplacé(s) à une autre position dans le puits pour une autre série d'émissions. Une telle acquisition est usuellement désignée par le terme anglais "walkaway", c'est-à-dire que l'émission est effectuée à déport croissant d'une émission à la suivante. En sismique de puits terrestre et/ou marine, on utilise plusieurs émetteurs séparés par une distance constante ou non et disposés en surface du milieu à explorer le long d'une ligne d'acquisition. Un outil de réception est descendu dans le puits dont la direction peut être verticale ou déviée. L'outil comprend généralement plusieurs récepteurs séparés par une distance constante ou non qui est généralement différente de la distance constante ou non séparant deux émetteurs ou sources consécutifs. Chaque récepteur comprend trois capteurs ou géophones qui sont orientés suivant trois directions différentes généralement orthogonales, une des directions étant déterminée pour disposer d'un axe de référence. L'axe de référence peut être non vertical ou vertical si on monte ces capteurs sur un système à cardan. Dans les deux cas, l'axe de référence V des capteurs est perpendiculaire au plan défmini par les deux autres directions H1 et H2 dont les positions ne sont pas connues généralement. Les axes V, Hl, H2 constituent un premier système d'axes. Les signaux représentatifs des ondes se propageant dans le milieu sont enregistrés en fonction du temps sous forme de traces. Les ondes qui se propagent dans le milieu sont principalement les ondes P montantes ou descendantes (ondes directes) et les ondes S. 2 2769718 Les traces enregistrées selon les axes V, Hl, H2 peuvent être regroupées en collections de traces à émetteur commun par exemple. Un traitement actuel des traces consiste à défminir un deuxième système d'axes R, N, T, dans lequel l'axe R est situé dans le plan d'énergie maximum. Une telle transformation est effectuée à l'aide de ce qu'on appelle usuellement une double rotation, une rotation d'angle 1H étant effectuée autour de l'axe V et la deuxième rotation d'angle (>v étant effectuée autour d'une direction dénommée Hmi qui est perpendiculaire à un vecteur Hm, qui est la projection de l'onde directe descendante dans le l0 plan défmi par H1 et H2, l'angle d)H étant l'angle que font entre eux H1 et Hmax. L'angle (v est l'angle formé par les axes V et R. Cette transformation du premier système d'axes V, H1, H2 en un deuxième système d'axes R, N, T est effectuée notamment à l'aide d'un logiciel bien connu des spécialistes et dénommé SEISLINK, commercialisé par la société WESTERN. On peut produire des traces dans le système d'axes R, N, T en projetant les composantes des traces repérées suivant les axes V, Hl, H2 et constituer des collections de traces à émetteur commun repérées en fonction de la profondeur z et du temps t, à partir d'au moins les traces suivant l'axe R (ou traces R). Les collections de traces, par exemple collections V, R ou N, sont ensuite traitées pour séparer les différents types d'ondes, surtout les ondes P descendantes et montantes et les ondes S montantes et descendantes. Une technique de séparation desdites ondes peut être mise en oeuvre en utilisant le logiciel SEPAR de la Compagnie Générale de Géophysique (CGG). Toutefois, cela nécessite de connaître les vitesses apparentes des types d'ondes recherchés. Généralement, ces vitesses apparentes sont déterminées à l'aide d'un pointé manuel sur les traces enregistrées. Les inconvénients de cette technique résident dans l'imprécision inhérente aux pointés manuels et parfois à l'impossibilité de discerner les types d'ondes directement sur les traces enregistrées. De plus, le coût du traitement est élevé pour des résultats peu précis. 3 2769718 Une autre technique préconisée par la société SCHLUMBERGER consiste à utiliser un logiciel qui est capable d'analyser les collections de traces à émetteur commun et d'en déduire les vitesses apparentes des ondes à séparer lorsqu'elles existent. Un inconvénient de cette technique réside dans son instabilité et surtout dans le fait que la séparation des ondes est peu efficace en raison du faible nombre de capteurs utilisé. Ces techniques de séparation d'ondes sont décrites respectivement dans les articles suivants: l0 B. SEEMAN and L. HOROWICZ: "Vertical seismic profiling: Separation of upgoing and downgoing acoustic waves in a stratified medium. GEOPHYSICS, Vol. 48, No. 5 (May 1983), p. 555- 568; C. ESMERSOY: "Inversion of P and SV waves from multicomponent offset vertical seismic profiles", GEOPHYSICS, Vol. 55, No. 1, (January 1990), p. 39-50. La présente invention a pour objet une nouvelle méthode de séparation des ondes élastiques, dans laquelle on met en oeuvre un traitement particulier des collections de traces qui permet de remédier aux inconvénients des méthodes antérieures telles que rappelées succinctement ci-dessus. La méthode selon l'invention pour la séparation des ondes se propageant dans un milieu est du type consistant à: - générer des ondes dans le milieu successivement à partir de plusieurs positions d'au moins un émetteur et situées à la surface dudit milieu, - recevoir lesdites ondes sur au moins un récepteur susceptible d'occuper successivement au moins deux positions dans un puits foré dans le milieu, chaque récepteur comprenant trois capteurs orientés dans un premier système d'axes (V, H1, H2), un desdits axes (V) étant orienté suivant une direction déterminée, chaque capteur délivrant un signal qui est enregistré en fonction du temps sous forme d'une trace qui est repérée par les positions du couple émetteur-récepteur ayant produit ladite trace, - constituer des collections C de traces à émetteur commun (z, t) à partir d'au moins un ensemble de traces V, H1, H2, 4 2769718 et elle est caractérisée en ce qu'elle consiste, en outre, à a) extraire les premières arrivées C signal d'au moins une des collections de traces C, b) réaliser à partir de ladite collection de traces C une collection de traces C résiduelles dans laquelle chaque trace résiduelle est la différence entre la trace C et l'estimation C signal, c) sélectionner dans la collection de traces C résiduelles une trace pivot située à une profondeur z, d) faire passer une droite de pente p variable en chaque point échantillonné de la trace pivot, e) sommer le long de ladite droite, pour chaque valeur de p, les amplitudes des événements sismiques situés aux intersections des traces avec ladite droite, f) reporter les sommes amplitude pour les temps et pente correspondants dans un domaine (z, p) dans lequel p représente les valeurs de ladite pente variant entre deux limites extrêmes et -r le temps, g) réaliser un sismogramme synthétique à partir des données du domaine (T, p), h) comparer le sismogramme synthétique à la collection émetteur commun de traces C résiduelles, i) procéder à des itérations jusqu'à ce que le sismogramme synthétique soit acceptable par rapport aux données initiales de ladite collection émetteur commun des traces C résiduelles. Selon une caractéristique avantageuse la méthode consiste, en outre, à défminir pour chaque couple de positions émetteur-récepteur un deuxième système d'axes (R, N, T) qui est calculé, par des moyens connus, à partir des composantes suivant les trois axes V, H1, H2, du vecteur déplacement correspondant à la première arrivée de l'onde d'intérêt et de la projection Hma, dudit vecteur déplacement dans le plan Hl, H2, l'axe R étant situé dans le plan V, Hmax, l'axe N étant perpendiculaire à l'axe R et situé dans le plan V, Hmax et l'axe T étant perpendiculaire au plan R, N, - produire des traces dans le système d'axes R, N, T en projetant les composantes des traces repérées suivant les axes V, H1, H2, 2769718 et à constituer les collections de traces (C) à émetteur commun (z, t) à partir des traces R, - extraire les premières arrivées R signal de la collection de traces R, - réaliser à partir de la collection de traces R une collection de traces R s résiduelles dans laquelle chaque trace résiduelle est la différence entre la trace R et l'estimation R signal, sélectionner dans la collection de traces R résiduelles une trace pivot située à une profondeur z, - effectuer les étapes d) à i) sur les traces R résiduelles, Selon une autre caractéristique, la méthode consiste, en outre, à: - définir un troisième système d'axes (X, Y, Z), fixe géographiquement, - déterminer dans le système d'axes (X, Y, Z) la déviation (W) et l'azimut (O) du puits, - aligner l'axe V du premier système d'axes (V, H1, H2) avec la direction du puits, - définir un vecteur hl perpendiculaire aux axes V et Z et contenu dans le plan Hl, H2, - déterminer l'angle ô entre le vecteur hl et l'axe H1. Selon une autre caractéristique, les composantes C ou R résiduelles, ainsi que les composantes H1 et H2 ou N et T sont projetées dans le système fixe (X, Y, Z) pour constituer des composantes résiduelles Xr et Yr, Zr, sur les axes X, Y, Z, la décomposition (t, p) étant réalisée à partir des composantes des traces Xr et Yr, Zr Selon une autre caractéristique, on calcule, par des moyens connus, pour chaque couple émetteur-récepteur et dans les domaines (t, p), un attribut angulaire polarisation et un attribut angulaire azimut Y. Lorsque l'onde d'intérêt est l'onde P, on utilise les attributs polarisation q, azimut y, angle de déviation y et azimut 0 du puits ainsi que la pente p pour calculer l'attribut vitesse Vp des ondes P dans le milieu. 6 2769718 A partir des attributs polarisation et vitesse Vp, on calcule l'attribut lenteur apparente pz par rapport à l'axe Z d'un événement sismique correspondant à des valeurs -c et p dans le domaine (t, p). Les attributs pz et Vp peuvent être utilisés pour opérer une sélection (mute) des événements à conserver dans le domaine (t, p), lesdits événements étant notamment ceux liés aux ondes P montantes. D'autres avantages et caractéristiques apparaîtront plus clairement à la lecture de la description de la méthode selon l'invention, ainsi que des dessins annexés sur lesquels: - la figure 1 est une représentation schématique du dispositif d'émission et de réception dans un milieu, - la figure 2 est une représentation schématique des directions de réception V, H1, H2 des capteurs d'un récepteur du dispositif de réception, - la figure 3 est une représentation schématique du repérage du puits de la figure 1 dans un système d'axes X, Y, Z fixe géographiquement, - les figures 4 et 5 sont des représentations des angles tH et 7 2769718 - la figure 12 représente schématiquement, dans le système d'axes X, Y, Z, la synthèse de la trace pivot de la figure 10, représentative de l'estimation du champ d'onde P montante. Selon la méthode de la présente invention, on utilise (figure 1) une ou plusieurs sources sismiques ou émetteurs 1 qui sont disposés à la surface 2 d'un milieu 3 dans lequel un puits 4 a été foré. Le puits 4 peut être vertical ou dévié comme sur la figure 1. Dans le puits dévié 4 est disposé un dispositif de réception 5 qui peut occuper au moins deux positions dans le puits. Le dispositif de réception peut comprendre un ou o0 plusieurs récepteurs, par exemple quatre récepteurs 6 constitués chacun par trois capteurs dont les directions sont suivant un trièdre rectangle formé par trois axes, V, HI, H2, l'axe V étant dirigé selon l'axe 7 du puits 4, les autres axes Hl et H2 étant situés dans un plan perpendiculaire à l'axe V (figure 2). La distance séparant deux positions consécutives de l'émetteur 1 est généralement constante et par exemple égale à 25 m, alors que la distance séparant deux récepteurs consécutifs 6 est généralement constante et par exemple égale à 10 m. Selon la présente invention, on détermine la déviation 'y et l'azimut 0 de l'axe 7 du puits 4 dans un système d'axes orthogonaux X, Y, Z, fixe géographiquement, l'axe Z étant vertical (figure 3). Les coordonnées du vecteur V sont données en fonction de ' et 0 par les formules suivantes, dans le système d'axes X, Y, Z fcos 0 sin y V sin 0 sin W [cos y les angles y et 0 étant compris respectivement dans les intervalles [0, rf/2] et [0, 2r]. Un nouveau vecteur hl est défini pour être perpendiculaire à la fois au vecteur V et à l'axe Z, et dont les coordonnées sont données par les formules suivantes, dans le système d'axes X, Y, Z sin0 hl -cosO O 8 2769718 Le vecteur h1 est alors contenu dans le plan H1, H2. De la même manière, on définit un vecteur h2 situé dans le plan Hl, H2 et perpendiculaire à h1. Les coordonnées du vecteur h2 dans le système d'axes X, Y, Z sont données par les formules suivantes [cos cos O h2 cos sin 0 sinm Les positions des axes H1 et H2, par rapport aux vecteurs h1 et 1o h2 sont données par les équations suivantes sin 0 cos ô + cos y cos 0 sin 6 H1 -cos 0 cos 6+cos ysin 0 sin - s sin sin et - sin 0 sin 8 + cos ucos 0 cos 8 H2 cos 0 sin 8 + cos y sin 0 cos 8 - sin vy cos 6 Dans ces formules, 6 représente l'angle formé par l'axe H1 et le vecteur h1 (figure 4). L'émetteur 1 occupe successivement plusieurs positions sur la surface 2, lesdites positions de l'émetteur 1 étant situées suivant des configurations géométriques données, rectangulaires ou circulaires concentriques autour du puits foré 4. Le déport entre l'émetteur 1 et le dispositif de réception 5 croît d'une émission à l'autre, de manière à effectuer une acquisition du type walkaway. L'émetteur 1 génère des ondes élastiques qui se propagent dans le milieu 3 et qui sont reçues et enregistrées sur chacun des capteurs de chaque récepteur 6 du dispositif de réception 5, les ondes étant enregistrées en fonction du temps t sous la forme de traces dans les directions V, H1 et H2. Ainsi, à chaque émission, on obtient un ensemble de traces qui sont repérées par la profondeur z du récepteur 6 et la position de l'émetteur 1. Lorsque l'enregistrement correspondant à une émission est terminé, on procède à une autre émission pour une autre 9 2769718 position de 1 'émetteur et ainsi de suite pour toutes les positions prédéterminées des points d'émission. Après quoi, on déplace le dispositif de réception 5 à une autre profondeur dans le puits 4, ainsi que cela est représenté sur la figure 1 sur laquelle le récepteur inférieur 6 vient occuper la position 6'. Une autre série d'émissions et d'enregistrements est alors effectuée pour cette autre position du dispositif de réception 5. Chaque trace enregistrée sur chaque capteur, c'est-à-dire dans les directions V, Hl, H2, est repérée par les positions du couple émetteur o0 1-récepteur 6 ayant produit ladite trace. Dans une autre étape, on définit pour chaque couple de positions émetteur-récepteur, un deuxième système d'axes orthogonaux R, N, T qui est calculé à partir des composantes suivant les trois axes V, Hl, H12 du vecteur déplacement correspondant à la première arrivée de l'onde 1s d'intérêt, par exemple l'onde P, et de la projection Hmax du vecteur déplacement dans le plan Hl, H2. La majeure partie de l'énergie associée à l'onde d'intérêt P ou onde directe se situe sur cet axe R qui correspond audit vecteur déplacement. L'axe N est perpendiculaire à l'axe R et il est situé dans le plan V, Hmax et l'axe T est perpendiculaire au plan R, N (figure 5). Le passage du système d'axes V, Hl, H2 au système d'axes R, N, T est effectué par une rotation d'angle tiH autour de l'axe V et une autre rotation d'angle Iv autour d'un vecteur Hmin qui est perpendiculaire au vecteur Hmax,. Les angles Hi et 1v sont calculés à partir des trois composantes V, Hl, H2 présentes dans une fenêtre temps située autour de l'arrivée directe. Les angles D et tDv peuvent être calculés en utilisant par exemple le logiciel SEISLINK (version 1996) commercialisé par la société WESTERN. Pour chaque couple émetteur-récepteur, les données initiales V, Hl, H2 fournies par les capteurs sont projetées dans le deuxième système d'axes R, N, T. Les composantes selon les axes R, N, T s'expriment à l'aide des équations suivantes: R= V cos (v +(H1 cos>H + H2 sin (pI) Sin (V N = -vcos(v T = - H1 sin (H + H2 Cos (DH 2769718 L'angle 6 est calculé pour chaque position de réception, à l'aide d'un processus d'inversion qui utilise les angles (I)H, les positions relatives de l'émetteur par rapport au récepteur, la déviation y et l'azimut 0 du puits. La valeur de l'angle 6 est calculée par minimisation de la fonction suivante: F = sin() cos( + cIH) sin(O - a) + sin(8 + H)[cos(Occ)sinCIcosP - cos( sin T] sources 1o dans laquelle ca représente l'angle formé par le plan contenant le couple émetteur-récepteur et l'axe X. Dans une autre étape et après avoir constitué des collections émetteur commun des composantes des traces suivant les axes R, N, T, (figure 7), on traite chaque collection émetteur commun de traces R pour extraire le champ d'onde P descendant. Pour cette opération, menée indépendamment sur chaque collection émetteur commun des traces R, on utilise la technique de filtrage matriciel (CEPHAG-IFP). Après horizontalisation des données à l'aide du pointé de la première arrivée, on applique cette technique sur toute la plage temps disponible ou sur une portion seulement (on peut ainsi éviter de faire intervenir un autre évènement descendant d'amplitude importante qui fausserait l'estimation - onde descendante S par exemple). La technique repose sur une décomposition des matrices interspectrales après moyennage fréquenciel et/ou spatial; les paramètres de filtrage sont choisis de façon à permettre l'extraction du champ P descendant dans les meilleures conditions. Le moyennage fréquenciel est réalisé de façon à ne conserver que les premiers termes des fonctions d'intercorrélation, aucun moyennage spatial n'étant effectué. Pour chaque collection initiale R, la présente invention consiste à fournir, en sortie, une collection émetteur commun dénomnmée R signal, ces signaux constituant l'estimation du champ P descendant pour la position source considérée. Une collection émetteur commun R résidu ou R résiduelle est obtenue par différence entre la collection initiale R et l'estimation R signal (figure 8). Les collections R signal ainsi estimées peuvent être utilisées ultérieurement pour le calcul des opérateurs anti-multiple. l l 2769718 Les composantes R résiduelles, N et T sont projetées dans le repère fixe X, Y, Z, en utilisant les angles (I)H, I)v et ô ainsi que la déviation y et l'azimut O du puits. Les composantes ainsi obtenues sont dénommées X résidu, Y résidu et Z résidu (figure 9). Pour chaque collection de traces émetteur commun, d'une part, et chaque composante X résidu, Y résidu et Z résidu, d'autre part, on procède à une décomposition dans un domaine (-c, p) en utilisant une technique d'inversion. La figure 10 représente schématiquement une collection de traces résiduelles émetteur commun. Dans cette collection de 1o traces résiduelles, on sélectionne une trace pivot 10 située à une profondeur z. On fait passer une droite 11 de pente p variable en chaque échantillon t de la trace pivot. On somme le long de ladite droite 11, pour chaque valeur de p, les amplitudes des évènements sismiques 12 situés aux intersections des traces 10 et 13 à 16 avec ladite droite 11. On reporte les sommes amplitudes A pour les temps et pentes correspondants dans le domaine (T, p) dans lequel p représente les valeurs de la pente et ' le temps correspondant à t. A partir des données de chacun des X résidu, Y résidu, Z résidu dans le domaine (T, p), on réalise un sismogramme synthétique que l'on compare à la collection correspondante X résidu, Y résidu ou Z résidu et on procède à des itérations jusqu'à ce que le sismogramme synthétique soit acceptable par rapport aux données de la collection correspondante X résidu, Y résidu ou Z résidu. Pour éviter l'étalement des évènements dans le domaine (-T, p), on applique dans le processus d'inversion une régularisation en norme Ll. La figure 11 représente la décomposition dans le domaine (-c, p) des collections de traces X résidu, Y résidu et Z résidu, de la figure 9. A partir de la décomposition dans le domaine (T, p) des composantes X résidu, Y résidu et Z résidu, on calcule l'attribut polarisation et l'attribut azimut y pour chaque valeur de p et Tc et ce, sur une fenêtre glissante, en utilisant les éléments propres de la matrice variance et covariance incluse dans le logiciel précité SEISLINK. Ces attributs q ou polar(c, p) et y ou azi(T, p) sont utilisés pour calculer l'attribut Vp (vitesse de l'onde P dans le milieu considéré), à l'aide de la formule suivante 12 2769718 Vp(TP) (cos[polar(T,p)]cos'' +sin[polar(T,p)]sinPcos[azi(T,p)-0]) La lenteur apparente pz(t, p) par rapport à l'axe Z, de 1'évènement considéré peut être également évaluée par la formule cos[polar(t, p)] Pz('Py) = Vp (T,P) L'axe Z étant vertical, la lenteur pz(T, p) est utilisée pour discriminer les évènements montants d'avec les évènements descendants. La connaissance des attributs Vp(T, p) et pz(t, p) permet de définir une sélection (mute) des évènements à conserver dans le domaine (t, p). Il peut être utile à cette occasion d'utiliser la connaissance que l'on peut avoir de la vitesse géologique Vp du milieu dans lequel est positionné le dispositif de réception. On peut par exemple appliquer un mute (sélection) dans le domaine (T, p) en utilisant les différents attributs précédemment décrits de façon à ne conserver que les évènements P montants. Ces évènements sont ensuite synthétisés pour construire un sismogramme dans le domaine (X, Y, Z, t) représentatif du champ d'onde P montant. Sur la figure 12, on a représenté l'estimation des composantes X, Y et Z du champ d'onde P montant correspondant à la trace pivot 10 à laquelle on a appliqué la décomposition dans le domaine (T, p). Un traitement identique peut être effectué pour extraire les champs d'onde S descendant et montant. Dans le cas o l'axe 7 du puits 4 est vertical (J = 0), on fixe arbitrairement la valeur 0 à zéro. La méthode selon l'invention est plus souple que les méthodes antérieures de filtrage optimal ou d'inversion paramétrique. Elle ne nécessite pas de pointé en-dehors de celui de la première arrivée et il n'est pas nécessaire de définir le nombre d'ondes à extraire; de plus, il n'est pas nécessaire de déterminer les vitesses apparentes. La méthode selon l'invention est par conséquent très simple d'utilisation et elle peut être utilisée telle quelle pour des acquisitions 3D. 13 2769718 De plus, elle fournit une possibilité d'analyser les données et de contrôler le résultat de la séparation. Bien entendu, il est possible de procéder à la décomposition (t, p) des données représentées dans le système d'axes R, N, T. Il est également possible de calculer l'attribut polarisation et l'attribut azimut y par rapport à ce repère R, N, T, ainsi que l'attribut Vp dans le domaine (T p). Les collections de traces R, N, T peuvent être remplacées par les collections C obtenues à partir des composantes V, H1, H2, les autres étapes de la méthode selon l'invention étant effectuées sur lesdites collections C de manière à obtenir les collections de traces C résiduelles nécessaires à la décomposition dans le domaine (t, p). De même, il est possible de projeter dans le système d'axes (X, Y, Z) les collections de traces C résiduelles et leurs collections de traces associées pour obtenir les collections de traces X résidu, Y résidu et Z résidu qui sont ensuite soumises à la décomposition dans le domaine (, p). La présente invention peut être appliquée pour l'imagerie fine des réservoirs ainsi que pour la calibration des mesures AVO (mesure du coefficient de réflexion versus angle d'incidence). En effet, à partir de données walkaway, il est possible de mesurer les champs incident et réfléchi et par déconvolution, de faire une estimation possible du coefficient de réflexion, cette mesure étant d'autant plus exacte que le capteur est situé près de l'objectif étudié. Par ailleurs, les analyses de polarisation conduisent à la connaissance des angles d'incidence. On a de ce fait, tous les éléments pour évaluer sans biais le comportement du coefficient de réflexion vis-à-vis de l'angle d'incidence. Enfm, la présente invention peut être utilisée pour la mesure de l'anisotropie au niveau de la ligne du dispositif de réception, ladite anisotropie pouvant être obtenue par une analyse conjointe de la vitesse apparente et de la polarisation de 1' arrivée directe. 14 2769718 REVENDICATIONS 1. Méthode pour la séparation des ondes se propageant dans un milieu, du type consistant à: - générer des ondes dans le milieu successivement à partir de plusieurs positions d'au moins un émetteur et situées à la surface dudit milieu, - recevoir lesdites ondes sur au moins un récepteur susceptible d'occuper successivement au moins deux de positions dans un puits foré dans le milieu, chaque récepteur comprenant trois capteurs orientés dans un premier système d'axes (V, H1, H2), un desdits axes (V) étant orienté suivant une direction déterminée, chaque capteur délivrant un signal qui est enregistré en fonction du temps sous forme d'une trace qui est repérée par les positions du couple émetteur- récepteur ayant produit ladite trace, - constituer des collections de traces (C) à émetteur commun (z, t) à partir d'au moins un ensemble de traces V, Hl, H2, caractérisée en ce qu'elle consiste, en outre, à a) extraire les premières arrivées C signal d'au moins une des collections de traces C, b) réaliser à partir de ladite collection de traces C une collection de traces C résiduelles dans laquelle chaque trace résiduelle est la différence entre la trace C et l'estimation C signal, c) sélectionner dans la collection de traces C résiduelles une trace pivot située à une profondeur z, d) faire passer une droite de pente p variable en chaque point échantillonné de la trace pivot, e) sommer le long de ladite droite, pour chaque valeur de p, les amplitudes des événements sismiques situés aux intersections des traces avec ladite droite, f) reporter les sommes amplitude pour les temps et pente correspondants dans un domaine (r, p) dans lequel p représente les valeurs de ladite pente variant entre deux limites extrêmes et r le temps, 2769718 g) réaliser un sismogramme synthétique à partir des données du domaine (X, p), h) comparer le sismogranmme synthétique à la collection émetteur commun de traces C résiduelles, i) procéder à des itérations jusqu'à ce que le sismogramme synthétique soit acceptable par rapport aux données initiales de ladite collection émetteur commun des traces C résiduelles. 2. Méthode selon la revendication 1, caractérisée en ce qu'elle consiste en outre à: 0 - définmir pour chaque couple de positions émetteur-récepteur un deuxième système d'axes (R, N, T) qui est calculé, par des moyens connus, à partir des composantes suivant les trois axes V, HI, H2, du vecteur déplacement correspondant à la première arrivée de l'onde d'intérêt et de la projection Hma, dudit vecteur déplacement dans le plan H1, H2, l'axe R étant situé dans le plan V, Hmax, l'axe N étant perpendiculaire à l'axe R et situé dans le plan V, Hmax et l'axe T étant perpendiculaire au plan R, N. - produire des traces dans le système d'axes R, N, T en projetant les composantes des traces repérées suivant les axes V, H, H2, et à constituer les collections de traces (C) à émetteur commun (z, t) à partir d'au moins les traces R, - extraire les premières arrivées R signal de la collection de traces R, - réaliser à partir de la collection de traces R une collection de traces R résiduelles dans laquelle chaque trace résiduelle est la différence entre la trace R et l'estimation R signal, - sélectionner dans la collection de traces R résiduelles une trace pivot située à une profondeur z, - effectuer les étapes d) à i) sur les traces R résiduelles. 3. Méthode selon la revendication 1 ou 2, caractérisée en ce qu'elle consiste en outre à: - définir un troisième système d'axes (X, Y, Z), fixe géographiquement, - déterminer dans le système d'axes (X, Y, Z) la déviation Y et l'azimut 0 du puits, 16 2769718 - aligner l'axe V du premier système d'axes (V, Hl, H2) avec la direction du puits, - définir un vecteur hl perpendiculaire aux axes V et Z et contenu dans le plan H1, H2, - déterminer l'angle 6 entre le vecteur hl et l'axe Hl. 4. Méthode selon la revendication 3, caractérisée en ce qu'elle consiste en outre à projeter dans le système fixe (X, Y, Z) la composante résiduelle (C ou R) et les composantes associées H1 et H2 ou N et T pour constituer des composantes résiduelles Xr et Yr, Zr, sur les axes X, Y, Z io la décomposition dans le domaine (t, p) étant réalisée à partir des composantes de traces Xr et Yr, Zr. 5. Méthode selon l'une des revendications 1 à 4, caractérisée en ce qu'elle consiste à calculer, par des moyens connus, pour chaque couple émetteur-récepteur, et dans le domaine (T,p), un attribut angulaire polarisation q et un attribut angulaire azimut y. 6. Méthode selon la revendication 5, caractérisée en ce que l'onde d'intérêt est l'onde P et en ce qu'on utilise les attributs polarisation , azimut y, angle de déviation y et azimut 0 du puits ainsi que la pente p pour calculer l'attribut (Vp) de la vitesse des ondes P dans le milieu. 7. Méthode selon la revendication 6, caractérisée en ce que l'attri- but lenteur apparente Pz par rapport à l'axe Z d'un événement sismique correspondant à des valeurs -c et p est calculé dans le domaine (T, p) à partir des attributs polarisation q et Vp. 8. Méthode selon la revendication 7, caractérisée en ce que les attributs pz et Vp sont utilisés pour opérer une sélection (mute) des événements à conserver dans le domaine (Tr, p). 9. Méthode selon la revendication 8, caractérisée en ce que les événements conservés sont ceux liés aux ondes P montantes. 10. Méthode selon la revendication 1, caractérisée en ce qu'on construit une section synthétique, à partir des évènements conservés, représentative du champ d'onde P montant. 11. Méthode selon les revendications 1 à 5, caractérisée en ce que l'onde d'intérêt est l'onde S 12. Méthode selon la reveMicaion 3, cAmaisée en ce que lorsque l'axe du puits est vertical on fixe l'aziu dudit puits à une valeur mlle.