FR 2494450 A1 19820521 FR 8024276 A 19801114 PROCEDE D'EXPLORATION SISMIQUE PAR LA TECHNIQUE DU PROFIL SISMIQUE VERTICAL ET INSTALLATION POUR SA MISE EN OEUVRE L'invention concerne un procédé pour l'étude des formations souterraines par la technique dite du profil sismique vertical, dans lequel un détecteur d'ondes acoustiques est placé dans un puits successi- vement à différentes profondeurs, des ondes acoustiques sont émises au moyen d'une source en surface et les signaux produits par le détecteur sont enregistrés. Un objectif essentiel de ces mesures est de déceler les hori- zons réflexifs ou miroirs situés plus bas que le fond du puits, grâce à l'analyse des ondes acoustiques réfléchies par ces miroirs et remon- tant vers le détecteur appelées ci-après ondes remontantes. Toutefois, les ondes détectées comprennent, en plus des ondes remontantes, les ondes qui se sont propagées directement de la source au détecteur, des ondes ayant subi des réflexions multiples et des ondes parasites de diverse nature. En regroupant tous les enregistrements en un document unique, on peut repérer des cohérences entre les différents tracés mais étant donné la multiplicité des composantes qui se superposent dans chaque signal, l'interprétation de ce document est extrêmement délicate. Un filtrage des signaux recueillis apparaît donc nécessaire pour améliorer le rapport signal/bruit et mieux mettre en évidence les ondes remontantes. Il est bien connu, afin d'accentuer les composantes d'ondes remontantes, de décaler les signaux de façon à aligner les ondes remon- tantes, et d'additionner les signaux ainsi décalés. Cette technique re- pose sur le fait que, pour deux signaux donnés si et Sk, le décalage entre les temps d'arrivée des ondes remontantes est de même valeur que le décalage tk relatif aux ondes descendantes et de sens contraire, sous réserve que le miroir ne présente pas un pendage trop accusé. Et les décalages tk, qui représentent le temps de propagation des ondes acous- tiques d'un niveau à l'autre, sont rès facilement accessibles, par exemnle à partir d'une diagraphie acoustique du puits indiquant la vites- se de propagation dans les formations traversées. Cependant, cette tech- nique de sommation ne fournit pas en elle-même une réjection suffisante du bruit. Par ailleurs, en ce qui concerne les niveaux de détection, la technique classique consiste à les choisir équidistants (voir l'ouvrage "Vertical Seismic Profiling" de E.I. Gal'perin, Society of Exploration Geophysicists, publication spéciale No.12, pages 11-12). La demande de brevet français No. 2 432 177 prévoit de traiter deux signaux correspondant à des niveaux consécutifs en les décalant pour aligner les ondes descendantes et en les soustrayant l'un de l'au- tre. Sous réserve que le décalage At ait une valeur déterminée, égale au quart de la pseudo-période des ondes acoustiques, on obtient la sup- pression de la composante d'onde descendante et le renforcement de la composante d'onde remontante. Mais ceci n'est vrai en fait que pour la fréquence correspon- dant à la pseudo-periode ci-dessus. De plus, le fait d'assigner une valeur constante à tous les temps At revient à fixer a priori les niveaux de détection, ce qui en pratique est un sérieux inconvénient, pour la raison suivante. La qualité des signaux recueillis dépend d'une façon détermi- nante du couplage acoustique entre le détecteur et la formation. Le couplage acoustique est lui-même fonction de l'état du puits au niveau considéré. Ainsi, si la mesure a lieu dans un trou découvert, le couplage d9tecteur-formation sera bon si le diamètre du puits est régulier et la formation compacte. Dans un trou tubé, c'est la liaison entre le tubage et la for- mation- c'est-à-dire la qualité de la cimentation- qui sera déterminante. Il est donc important, pour obtenir des signaux de détection de qualité optimale, de pouvoir choisir librement les niveaux de détec- tion et donc de pouvoir s'écarter de la règle énoncée dans la demande de brevet précitée. Dans la même optique, il est souhaitable de pouvoir éliminer les signaux de détection dont la qualité apparaît déficiente. Cette possibilité n'existe pas si le procédé de filtrage utilisé repose sur l'hypothèse que tous les décalages entre signaux de niveaux consécutifs sont égaux. L'invention vise donc à améliorer le rapport signal/bruit des signaux de détection de façon à accentuer très fortement les composantes d'onde remontante et cela en autorisant un libre choix des niveaux de détection et l'élimination éventuelle des signaux de qualité déficiente. Selon l'invention, le filtrage des signaux de détection com- prend les opérations suivantes: - prélever parmi les N signaux à traiter un groupe de n signaux s... sk...sn détectés à des niveaux xl...xk...x successifs; - déterminer les temps de propagation tk des ondes acoustiques du niveau x à chaque niveau xk - décaler chaque signal sk en l'avançant du temps tk par rapport au signal s1 de façon a aligner les composantes d'onde descendante; additionner le signal s et les signaux sk ainsi avances pour obtenir une somme z1; - décaler chaque signal sk en le retardant du temps tk par rapport au signal s1; - additionner le signal s et les signaux sk ainsi retardés pour obtenir une somme y1; - générer un signal ul, défini comme une estimation optimale de la com- posante d'onde remontante pour l'ensemble des signaux s.8, a partir des signaux y1 et z1; - enregistrer ce signal u1; - et répéter ces opérations pour obtenir des signaux u2...u...uNn+1 Le mode de filtrage selon l'invention n'impose aucune condi- tion aux temps tk et par conséquent les niveaux de détection xk peuvent %tre choisis en fonction de l'état du puits de manière à obtenir un cou- plage acoustique optimal entre le détecteur et la formation. Pour la même raison, rien n'interdit d'éliminer certains des signaux recueillis, s'il apparaît que leur qualité est insuffisante. Cette possibilité contribue à améliorer le rapport signal/bruit des signaux obtenus. Pour déterminer les signaux un, l'invention envisage deux so- lutions différentes. La première consiste à passer dans le domaine des fréquences. On obtient u. a partir de sa transformée de Pourier U., que l'on déter- i i mine par une relation simple à partir des transformées de Fourier Y. et i Z. respectivement de Yi et z.. Dans la seconde solution, on opère uni- quement dans le domaine des temps. On obtient alors une valeur ui appro- chée directement à partir de Yi et zio. Le caractère approximatif des résultats est compensé par une grande simplicité du traitement. Une autre caractéristique avantageuse du procedé selon l'in- vention est qu'il permet de déterminer, tout aussi simplement, des esti- mations optimales d.i des composantes d'onde-descendante. La connaissan- ce des valeurs d.i est nécessaire pour le traitement de déconvolution qui est appliqué de façon classique à la suite du filtrage. Au sujet de la déconvolution, on se reportera à l'article de P. Kennett et al "Vertical Seismic Profiles: Their Applications in Exploration Geophy- sics", Geophysical Prospecting, 1980, 28, pages 684 et suivantes. L'invention sera bien comprise à la lecture de la description ci-après, faite en se référant aux dessins annexes. Dans les dessins - la figure 1 représente schématiquement une installation d'exploration sismique pour l'obtention d'un profil sismique vertical; la figure 2 montre de façon très simplifiée le trajet des ondes acoustiques descendantes et remontantes, pour deux positions dif- férentes du détecteur; - la figure 3 représente sous forme de schéma-blocs un disposi- tif selon l'invention pour le filtrage des signaux de détection; - la figure 4 illustre l'effet de filtrage procuré par le dis- positif de la figure 3; - les figures 5 et 6 montrent à titre d'exemple un PSV avant filtrage et après filtrage à l'aide d'un dispositif de la figure 3; - la figure 7 représente une autre forme de réalisation du dis- positif de filtrage. La figure 1 représente une installation d'exploration sismique qui comprend une tour de forage 11 disposée au-dessus d'un puits 12 foré dans des formations terrestres 13. Un appareil de fond 14 comprenant un capteur d'ondes acousti- ques tel qu'un géophone est suspendu dans le puits 12 au moyen d'un câble électrique 15 passant sur des poulies solidaires de la tour de forage. Le câble 15 permet de déplacer l'appareil 14 dans le puits et assure en mmee temps la transmission des signaux de détection produits par le cap- teur vers un équipement de surface 16. L'équipement de surface 16 com- prend de façon classique un treuil pour enrouler le câble et des moyens pour déterminer la profondeur à laquelle se trouve l'appareil de fond 14, et des moyens pour traiter et enregistrer les signaux de détection trars- mis par le câble 15. L'appareil de fond 14 comprend de façon également classique un élément d'ancrage mobile 17 qui peut s'écarter de l'appareil 14 et s'ancrer dans la paroi pour assurer un contact approprié entre le cap- teur et la paroi du puits. Cela ne se produit que lorsque l'appareil 1h est arrivé à une profondeur o une mesure doit être exécutée, l'équipe- ment de surface 16 émettant alors un message pour commander le déploie- ment de l'élément d'ancrage. Une source d'ondes acoustiques 18 est placée à la surface, à une distance déterminée de la tête de puits. Cette source peut être constituée par tout dispositif approprié tel qu'un canon à air. La source 18 est commandée à partir de l'équipement de surface 16 pour effectuer un "tir" lorsque l'appareil de fond 14 est ancré à l'un des niveaux de profondeur prévus. Plusieurs tirs sont effectués pour chaque niveau, et l'ensemble des signaux ainsi produits sont mis en mémoire et additionnés dans l'équipement de surface 16 de manière que pour chaque niveau, on obtienne un signal dans lequel les composantes significatives soient renforcées au détriment des composantes de bruit aléatoires. On conviendra dans la suite d'appeler ce signal le signal de détection. A partir de l'ensemble des signaux de détection obtenus pour les différents niveaux, on produit des enregistrements graphiques ou traces sismiques que l'on regroupe en un document unique appelé profil sismique vertical (PSV) dont la figure 5 montre un exemple. La figure 2 représente de façon très simplifiée le trajet des ondes acoustiques détectées par le capteur à deux niveaux. La référence désigne un horizon réflexif ou "miroir" plus profond que le fond 22 du puits. Le miroir 20 est constitué par l'interface entre deux couches d'imoédances acoustiques nettement différentes. Le capteur placé en au niveau xI reçoit une onde acoustique qui se propage directement depuis la source, appel'eonde descendante, et une onde réfléchie par le miroir 20, dite onde remontante. C2est cette onde remontante qu'il faut mettre en évidence pour pouvoir déterminer la profondeur du miroir et son coefficient de réflexion. Le signal de détection sl(t) obtenu au niveau x1 est donc la somme d'une composante d'onde descendante dl(t) et d'une composante d'onde remontante ul(t) et comprend en outre des composantes de bruit dont certaines proviennent de réflexions multiples sur des miroirs si- tués entre la surface et le fond du puits. Le capteur placé en k au niveau xk recevra de même une onde descendante et une onde remontante, mais l'onde descendante arrivera plus tard et l'onde remontante plus tôt. Le décalage tk entre les temps d'arrivée des ondes descendantes vaut t = k xl k V vek étant la vitesse moyenne de propagation des ondes acoustiques dans la tranche délimitée par les niveaux x et xk. Le décalage t' entre les temps d'arrivée des ondes remontantes est égal à - t. si le miroir 20 est sensiblement horizontal. Un terme correctif doit Ztre ajouté si le miroir 20 présente un pendage important. Nous supposerons toutefois dans un premier temps qu'il n'y a pas lieu d'introduire un tel correctif. Compte tenu de ce qui précède, le signal sk(t) obtenu au ni- k veau xk comprend une composante d'onde descendante dk (t - tk) et une composante d'onde remontante uk (t + tk) Sur les diagrammes de la figure 4, r6alisés à partir de donnes synthétiques, comme sur le profil sismique vertical de la figure 5, résultant de mesures réelles, on observe des pics correspondant aux ondes descendantes et on peut vérifier qu'ils se décalent régulièrement d'une trace à l'autre. En revanche, dans les deux cas, les composantes d'ondes remontantes ont un rapport signal/bruit trop faible pour qu'on puisse les détecter. Un filtrage des signaux recueillis est donc nécessaire pour mettre en évidence l'onde remontante avec une netteté suffisante. Le principe de filtrage propos6 est le suivant. Parmi les N signaux de détection à traiter s1...sn, prélevons une série de n signaux, par exemple sl..Sk...sn (nN). Posons qu'il existe des estimations t et d des composantes d'onde remontante et d'onde descendante qui sont optimales pour l'ensem- ble des n signaux ci-dessus. Compte tenu de ce qui précède, on peut écrire, suivant la formulation des moindres carrés, qu'il faut minimiser l'expression n 12 E sk(t) - u1 (t + t) - d (t - tk) (I) k=l Dans le domaine des fréquences, l'expression (I) ci-dessus devient E IS..cUaDI 2I k=l k kUl o Sk, l et D1 sont les transformées de Fourier de sk(t), u(t) et dl(t) et ak et "k sont respectivement les opérateurs exp(- itkw) et exp(itkw). On a bien entendu tl =0 et a1 = "1 = 1. En annulant les dérivées par rapport à et D de l'express1n (II) on obtient les équations n n E k k U. k - = 0 (III) k=l k=l k n n Di C 0 (IV) E k=l akSk -a-n = k1- k=l On en tire, en résolvant les équations (III) et (IV) en l i n n n s c - e k a Sk Sk d l k=l k=l k k). 2 n 2 n nE ak E n -- Z 7 k2 k=l k=1 Dak k est la transformée de Fo urier du signal s k(t + tk), obtenu en avançant le signal sk (t) du temps t par rapport au signal' s (t). En d6calant ainsi sk par rapport à Sl, on réalise l'alignement des composantes d'onde descendante. La somme. est donc la transformme de Fourier Z d'un signal E (t) obtenu en additionnant le 1 1 signal s1(t) et les signaux sk (t + tk)- De sm.e ú Sk est la transfortue de Fourier du signal sk(t - t) obtenu en retardant s k(t) par rapport sl(t), ce décalage produisant l'alignement des composantes d'onde remontante. n La somme E k S est donc la transformée de Fourier Y, d'un k=l k k 1 signal yl(t) obtenu en additionnant le signal s1(t) et les signaux sk(t - t). L'expression devient alors n 2 nY1 - Z E1 ak 1 1 k=l k U1 = (VI) 2 n 2 n 2 n - E. ak E "k k=l k=l et à partir de L-, on calcule u par transformation de Fourier inverse. 1 ' Les temps tk peuvent 'être déterminés par simple lecture à par- tir d'une diagraphie acoustique du puits. Toutefois, un mode de calcul pré erable consiste à détecter les pies correspondant aux premières arr:i- vées des ondes descendantes et à mesurer le temps qui les sépare. Le procédé ci-dessus fournit donc de façon simple une estima- tion optimale des composantes d'onde remontante. On recommence ensuite ces opérations pour déterminer un signal u2 partir des signaux s2..s n+ et ainsi de suite jusqu'au groupe SN 11..sN. Finalement, à partir de N signaux à traiter, on obtient une série de (N-n+l) signaux u.. On obtient de façon correspondante D1l par l'expression n -2 ÀnZ 1 -- Y1 7 ak D 1 1 k=l (k VII) i 2 n a2 2 n -- a k' kk k=l k=l etenentire la transformée- de Fourier inverse d, estimation optimale des composantes d'onde descendante, puis successivement tous les signaux d.i. Les signaux d. sont nécessaires pour l'application de la déconvo- lution à laquelle on procède de façon classique apres le filtrage. En alignant les enregistrements des signaux u sous la forme d'un PSV, on observera que les pics correspondants sont parfaitement alignés et très nets si le miroir 20 est sensiblement horizontal. Si l'on constate que les pics sont moins nets et sont décalés l'un par rapport à l'autre, il faut supposer que le miroir 20 présente un pendage assez accusé. Dans ce cas, on détermine les écarts ck entre la position des pics et la position qu'ils occuperaient en l'absence de pendage et l'on recommence les opérations précédentes en remplaçant chaque opérateur a k par un opérateur yk = exp (i(tk + Sk)W). Le procédé de filtrage décrit est mis en oeuvre à l'aide du dispositif de la figure 3. Les signaux de détection ont été préalablement enregistrés et examinés pour que les signaux de qualité déficiente soient éliminés. Il reste après cette sélection N signaux sl...sN que l'on ins- crit dans une mémoire 30, La mémoire 30 est commandée pour fournir une série de n signaux correspondant à des niveaux successifs. On envisagera ici le traitement de la première série de signaux sl.. Sk...sn. Les signaux s2...sk... sn sont appliques à un premier circuit de décalage 31 qui les retarde res- pectivement des temps t2...tk...tn par rapport à s1. Le signal sl(t) et les signaux sk(t + tk) obtenus sont appliqués à un addtionneur 32 qui délivre le signal-somme yl (t).Les signaux s2...s sont également appli- qués à un second circuit de décalage 33 qui les avance par rapport à sI respectivement des temps t2...tk...tn. Le signal sl(t) et les signaux sk(t - tk) sont appliqués à un additionneur 34 qui fournit le signal- somme zl(t). Les circuits de décalage 31 et 33 peuvent être constitués de façon connue pour l'homme de l'art, par un ensemble de registres à déca- lage. Les temps tk sont déterminés en mesurant l'intervalle de temps en- tre les premières arrivées des ondes descendantes pour les signaux s1 et sk et sont enregistrés dans une mémoire 35. Des dispositifs 36 et 37 fournissant les transformées de Fourier délivrent des signaux représentatifs des fonctions Y1 et Z 1 partir des sommes Y1 et z1. Un calculateur 38 est programmé pour fournir les grandeurs 1 et D1 à partir des signaux Y et Z et des temps tk conformément 1 1 1 i k aux relations (VI) et (VII) ci-dessus. Un dispositif 39 fournissant la transformée de Fourier inverse délivre un signal u à partir du signal. Le signal * est inscrit dans une mémoire 40. Un dispositif 41 fournit de mme la transformée de Fourier inverse d, qui est inscrite dans une mémoire 42. Ce traitement est ensuite répété en décalant d'un niveau, c'est-à-dire avec la série de signaux s2... sn+l, et ainsi de suite jus- qu'à ce que le dernier signal de la série soit le signal sN. Les signaux dé obtenus seront utilisés dans un traitement de déconvolution ultérieur. il Les signaux u. obtenus sont appliqués à un enregistreur gra- phique 43 qui produit un profil sismique vertical regroupant tous les diagrammes de signaux. La figure 6 montre le profil sismique obtenu à partir des traces de départ constituant le profil sismique de la figure 5. Les diagrammes de la figure 4 illustrent les résultats obtenus à l'aide du procédé de filtrage décrit ci-dessus. Les diagrammes 4-1 à 4-11 représentent des signaux de détec- tion correspondant à Il niveaux successifs réalisés à partir de données synthétiques. Tous les temps de propagation entre niveaux consécutifs sont égaux. On observe sur toutes les traces, à l'exception des traces 4- 2, 4-5 et 4-8, des pics correspondant aux ondes descendantes. Les ondes remontantes en revanche ne peuvent pratiquement pas être distinguées du bruit. Pour le filtrage, on élimine les traces 4-2, 4-5 et h-8 dont la qualité est manifestement déficiente. On obtient par le procédé de filtrage ci-dessus un signal u qui présente un pic (U) très marqué correspondant à l'onde remontante. La composante d'onde descendante (D) est au contraire très atténuée. Cet exemple montre que le procédé de filtrage selon l'inven- tion n'impose aucune condition précise en ce qui concerne les temps de propagation et donc les niveaux de détection. En particulier, il n'est nullement nécessaire que les temps de propagation entre niveaux consécu- tifs soient égaux. Ainsi, sur la figure 4, si l'on considère uniquement les si- gnaux effectivement utilisés, les temps de propagation entre niveaux consécutifs ne sont pas tous égaux, mais on obtient un rapport signal/ bruit tout à fait satisfaisant. Le fait que le procédé selon l'invention ne repose par sur l'hypothèse de temps de propagation égaux entre niveaux consécutifs est un avantage pratique très important. Cela permet en premier lieu de choisir librement les niveaux de détection de façon à optimiser le couplage acoustique du détecteur avec la formation. Si la mesure a lieu en trou découvert, on pourra déterminer les zones les plus favorables à l'aide d'une diagraphie de diamètre, indiquant les zones de diamètre régulier, et d'une diagraphie acoustique indicative de la vitesse de propagation des ondes acoustiques et donc de la compacité des formations traversées. Si l'on opère en trou tube, le facteur déterminant pour le couplage sera la qualité de la liaison du ciment au tubage. On pourra déterminer les zones favorables à l'aide d'une diagraphie telle qu'obte- nue par le dispositif selon le brevet français No. 2 400oo 613. La liberté dans le choix des niveaux de détection n'existe toutefois qu'à l'intérieur d'une certaine plage. Au-delà d'un certain intervalle entre niveaux consécutifs, la corrélation entre traces consécutives deviendrait moins fiable. On admet généralement que cet intervalle ne doit pas dépasser environ le tiers de lailongueur d'onde des ondes acoustiques (voir l'ouvrage précité de E.I. Gal'perin, page 27). Pour.fi'xer les idées, on peut indiquer un écart maximal d'environ 40 mètres. On peut d'autre part considérer comme raisonnable un écart mi- nimal d'environ 10 mètres. En ce qui concerne le nombre n de signaux que l'on utilise pour déterminer chaque signal u., le choix dépend du rapport signal/bruit souhaité. On préfère par commodité choisir un nombre impair, et un bon compromis entre la -qualité du signal produit et la durée du traitement est obtenue avec n = 7, 9 ou 11. Une variante simplifiée du procédé de filtrage consiste à dé- terminer les signaux u. sans passer dans le domaine des fréquences, les valeurs obtenues n'étant alors que des approximations. Si l'on part de l'expression (I) ci-dessus n 2 n Isk(t) - u (t + tk) d(t - tk)2 et qu'on annule sa dérivée par rapport à d*(t), on obtient n k= sk t)-d () i (t + 2tk) = 0 (VIII) E (sk(t + Y k-d (t)-u t +2tk)=0 k=1 c 'est-a-dire n z(t) - nd (t) - ut k=l On obtient de façon correspondante, rapport à u (t) i 1 Z u (t) - nY(t) + E n i=1 (t + 2tk) = 0 (IX) en annulant la dérivée de (I) par d (t - 2ti) = O (x) En exprimant le troisième terme de cette équation à l'aide de l'6quation (IX), et en reportant l'expression obtenue dans l'équation (X), on obtient 1t(t=yt) n 1n n u (t)= - y(t) E z(t-2 -- 1 1n2 12 E ui(t-2ti-tk) k=l k=l i=l (XI) En négligeant le troisième terme, on arrive à la valeur appro- ximative ()=1 y(t)-12 (t-n u (t≥ -y(t)- 2 E z (t - 2t) k=l (XII) Le traitement nécessaire pour déterminer cette variante est extremement simple puisqu'il ne transformation de Fourier. On obtient de façon correspondante, pour la valeur approximative les signaux u. dans i comporte aucune les ondes remontantes, d* (t) 1(t -y(t+ n d (t) = - z (t)- 1 E y(t+2t) n n2 k=1 k (XIII) La figure 7 représente le dispositif pour la mise en oeuvre de cette variante. Le bloc 50 sur la figure 7 regroupe l'ensemble des éléments -34 de la figure 3 et fournit les signaux yl(t) et zl(t). * Un circuit de dScalage 51 recevant le signal zl(t) fournit les signaux décalés z1 (t - 2tk), les temps 2t étant fournis par une mémoi- k k re 52. Les signaux zl(t) et zl(t - 2tk) sont appliqués à un additionneur 53. La somme obtenue est divisée par n2 dans un diviseur 54. Le signal yl(t) est divisé par n dans un diviseur 55, et les sorties des diviseurs 53 et 54 sont reliées à un soustracteur 56 qui fournit le signal (t) recherché, lequel est inscrit dans une mémoire 57. Le signal Y1 (t) est appliqué d'autre part à un circuit de décalage 61 qui fournit les signaux yl(t + 2tk), qui sont appliques avec le signal yl(t) à un additionneur 63 dont le signal de sortie est divisé par n dans un diviseur 64. Le signal z1 (t) est divisé par n dans un diviseur 65, et un soustracteur 66 recevant les sorties des diviseurs 64 et 65 fournit le signal dl *(t) qui est inscrit dans une mémoire 67. REVENDICATIONS 1. Procédé d'exploration sismique dans lequel on filtre, pour en améliorer le rapport signal/bruit, des signaux produits par un détecteur d'ondes acoustiques placé à différentes profondeurs dans un puits, chaque signal étant produit en réponse à une émission d'ondes acoustiques à partir d'une source en surface, lesdits si- gnaux comprenant des composantes d'onde descendante et des composan- tes d'onde remontante, caractérisé par le fait que le filtrage com- prend les opérations suivantes: - prélever parmi les N signaux à traiter un groupe de n signaux s...sk...sn détectés à des niveaux xl...xk... xn successifs; - déterminer les temps de propagation tk des ondes acoustiques du niveau x1 à chaque niveau xk; - décaler chaque signal s en l'avançant du temps tk par rapport au signal s de façon à aligner les composantes d'onde descendante; - additionner le signal s et les signaux sk ainsi avancés pour obtenir une somme z1; - - décaler chaque signal sk en le retardant du temps tk par rapport au signal s; - additionner le signal s et les signaux 5k ainsi retardés pour ob- tenir une somme y1 - générer un signal u1, défini comme une estimation optimale de la composante d'onde remontante pour l'ensemble des signaux si...s à partir des signaux y1 et z1; - enregistrer ce signal ul; - et répéter ces opérations pour obtenir des signaux u..... 2. Procédé selon la revendication 1, caractérisé par le fait que le signal u1 est obtenu par transformation de Fourier inverse de l'expression U1 reliée à y1 et z1 par la relation n 2 nY 1- Z1 À Z ak k=l U1 nn n - r- _ 2 a2 z a2 - k=l k=l o Y. et Z. sont les transformées de Fourier de Yi et Z. et ak et __ 1 11 ak sont respectivement les opérateurs exp(-itkw) et exp(itkw). 3. Procédé selon la revendication 2, caractérisé par le fait que l'on génère en outre un signal dl estimation optimale de la composante d'onde descendante pour l'ensemble des signaux s1...s, a partir des signaux y1 et zl, par transformation de Fourier inverse du signal D obtenu conformément a la relation n n-2 nZ1 Y1 ' ak k=l 1 2 n 2 n n -- Z k ' k k=l k k=l k 4. Procédé selon la revendication 1, caractérisé par le fait qu'une valeur approchée du signal u.i est calculée par la relation: i i 1n Ul = nyl-n--2E z1 (t - 2t) n k=l 1 k 5. Procédé selon la revendication 4, caractérisé par le fait que l'on calcule une valeur approchée du signal d, estimation optimale de la composante d'onde descendante pour l'ensemble des signaux s... sn, à partir des signaux Y1 et Zl, par la relation 1 n t _ 1 d = n Zl- 2 E Y (t + 2t) k=l k 6. Procédé selon l'une des revendications 1 à 5, caractérisé par le fait que l'on fixe les profondeurs de détection, à l'intérieur d'une plage déterminée, de manière à optimiser le couplage acoustique entre le détecteur et la formation. 7. Procédé selon l'une des revendications 1 à 6, caractérisé par le fait que l'on élimine avant le filtrage les signaux de détec- tion de qualité déficiente. 8. Installation d'exploration sismique pour la mise en oeu- vre du procédé selon la revendication 1, comprenant un détecteur d'ondes acoustiques et des moyens pour déplacer ce détecteur dans un puits, une source d'ondes acoustiques, des moyens pour enregistrer les signaux produits par le détecteur, lesdits signaux comprenant des composantes d'onde descendante et des composantes d'onde remontante, et des moyens pour filtrer ces signaux afin d'en améliorer le rapport signal/bruit, caractérisée parle fait que les moyens de filtrage com- prennent - une mémoire (30) contenant les N signaux à traiter et commandée pour fournir une série de n signaux s1.. sk**.s correspondant à des niveaux successifs xl......xn; - une mémoire (35) contenant les temps de propagation tk des ondes acoustiques du niveau x1 à chaque niveau xk; - un premier circuit de décalage (31) pour avancer les signaux sk respectivement du temps tk par rapport au signal si; - un premier additionneur (32) fournissant la somme z1 du signal si et des signaux sk ainsi avancés; - un second circuit de décalage (33) pour retarder les signaux sk respectivement du temps tk par rapport au signal s1 - un second additionneur (34) fournissant la somme y1 du signal si et des signaux sk ainsi retardés; - des moyens de calcul (36-39; 51-55) pour déterminer à partir des signaux y1 et zl un signal u.., défini comme une estimation optima le de la composante d'onde pour l'ensemble des signaux s...s - et des moyens d'enregistrement reliés aux moyens de calcul. 9. Installation selon la revendication 8, caractérisée par le fait que les moyens de calcul comprennent des dispositifs (36,37) reliés aux additionneurs (32, 34) et fournissant les transformées de Fourier Y et Z1 de Y, et Zl, un calculateur (38) fournissant la grandeur 1 par la relation: n 2 nY1 -Z. lXE ck U1 = - k=l 2 n 2 n 2 n - Z ak k=l k=l1 o i et k sont respectivement les opérateurs exp(-itkw) et exp(itkw), et un dispositif (39) fournissant la transformée de Fourier inverse et délivrant le signal u à partir de U1. 10. Installation selon la revendication 9, caract6risée par le fait que le calculateur (38) génère 6galement un signal D1 conforme- ment à la relation n nZl-Y1 Z ak2 D= k=l 2 n n --_ Z ak' Z"k k=l k=l ce signal étant appliqué au dispositif (39) fournissant la fansfor- mee de Fourier inverse qui délivre un signal d1t, estimation optima- le de la composante d'onde descendante pour l'ensemble des traces s.s... s 1 k n 11. Installation selon la revendication 8, caract6riséepar le fait que les moyens de calcul comprennent un dispositif de décalage (51) fournissant à partir du signal zl(t) des signaux décalés zl(t - 2tk), un additionneur (53) fournissant la somme des signaux zl(t) et zl(t - 2tk), un circuit (54) divisant cette somme par n, un circuit (55) divisant par n le signal yl et un soustracteur (56) relié aux deux circuits diviseurs (54, 55) et fournissant un signal u1 approché. 12. Installation selon la revendication 11, caractériséepar le fait que les moyens de calcul comprennent en outre un dispositif de décalage (61) fournissant à partir du signal yl(t) des signaux décalés y1(t + 2tk), un additionneur (63) fournissant la somme des signaux y1(t) et yl(t + 2tk), un circuit (64) divisant cette somme par n, un circuit (56) divisant par n le signal z et un soustrac- teur (66) relié aux deux circuits diviseurs (64, 65) et fournissant un signal dl1 représentant de façon approximative une estimation optimale de la composante d'onde descendante pour l'ensemble des traces s1... sk.. sn