FR 2629603 A1 19891006 FR 8904199 A 19890321 METHODL Lú TRAITEIMNT DE DONNERS SISMIQUES LP présente invention est relative à une méthode pour le traitement ae données sismiques, et se rapporte plus particulièrement, mais non exclusivement, à une méthode perfectionnee pour la trznsformation d'une série de trrces sismioues de dérivation non-zéro distribuées spstilement en une nouvelle serie de traces sismiques de sortie non- ambulatoires, de dérivation zéro, resolutives en dépression, representatives de lI configuration tri-dimensionnelle (3-D) d'un volume de la terre. C'est chose bien connue que le recueil de Point Médian Commun (CMP), (parfois appelé également Point ae Profondeur Commun ou Point de Réflexion Commun, CDP ou CRP respective- ment) de données de cheamp sismioue. Voir par exemple les brevets U.S. Mendenhall et all, 3.217.828, et Mayne, 2.732.906, oui sont Incorporés dans le présent mémoire des- criptif a titre d'indication de la situation technicue en ce oui concerne le CMP. Le brevet U.S. 3.217.828 divulgue un traitement de données bi-dimensionnel (2-D) dàns leouel - des couches de terre dépressives sont projetées dans un plan bidimensionnel, perpendiculaire à la surface de la terre, le long d'une limne de profil dasignee. Il con'ilent aussi ae signaler notamment à ce sujet le brevet U.S. 4.742.4'7. Le terre est un objet à trois oimensions. C'est la con- figuration d'un volume de la terre oui présente de l'inté- rêt pour certographier des structures tri-dimensionnelles susceptibles ce contenir oes ressources minérales ovant de la valeur. Dans le traitement 2-D (bi-dimensionnel, c'est seule-ent le composante dépressive parallèle à la ligné de levé de plan oui est cartopraphiée. Les composantes dépres- sives de couche de terre se situent en dehors du plan 2-D sont déformees. Les limitations des calculs 2-D traditionnels se pré- sentent sous trois aspects: Premièrement, en ce qui con- cerne un réflecteur de dépression, le point de réflexion véritable pour 1l réflexion a'une déviation ou dérivation non-7éro se situe en haut de dépression à partir du point médian entre une source sismique et un récepteur de signaux sismicues. Lorscue les données provenant des différentes dérivations sont recueillies, il se présente une maculature du point de réflexion. Deuxièmement, les vitesses de recueil des données sont en fonction des dépressions, et en conséquence, lorsque des incidents avec dépressions différentes se croisent sur une section sismique en traitement, aucune fonctionAde vi- tesse unitaire n'alignera correctement les différents in- cidents sur le rassembleur CMP (Point Médian Commun) à cet emplacement précis. Cette situation a pour résultat Oue le mode de recueil traditionnel recueille de façon erronée les incidents pour lesquels la vitesse choisie n'était pas ap- propriée. Troisièmement, les vitesses de recueil sont en fonction de l'azimut, ce nui veut dire oue la vitesse ce recueil dépend de l'orientation de la dépression maximum par rap- port au plan vertical entre la source sismiQue et le récep- teur. Etrnt donné que l'exploration 2-D (bi-dimensionnelle) se limite à de seules lignes de profil, la résolution 3-D (tri- dimensionnelle) des détails d'un volume de la terre, basée sur l'exploration 2-D, est souvent plus ou moins déformée. La couverture Tonale 3-D se présente frécuemment sous la forme d'une grille de lignes sismicues se composant essen- tiellement de deux séries de lignes disposées ort.hoonale- ment l'une par rapport à l'aeutre. Chaoue ligne inclue une pluralité de récepteurs de signaux sismioues parsemée d'une pluralité de sources sismiques qui sont disposées sur des emplacements désirés à l'intérieur de la grille. Toutefois, la grille peut être rectangulaire, circulaire, hexagonale, en forme d'étoile, ou de toute autre configuration géomé- trique désirée. La grille de sources et de récepteurs est habituellement disposée horizontalement sur la surface de 1P terre, mais, tout au moins dans le plan conceptif, il n'existe aucune raison pour oue des dispositions verticales soient exclues. Avec des adaptations appropriées en ce qui concerne l'équipement, des opérations 3-D peuvent être ef- fectuées sur terre ou sur mer. Dens les limites d'un rayon raisonnable de deux milles (3 km 5) ou à peu près, l'amorçage de toute source cuelcon- que isonifiera un système entier de capteurs ou récepteurs sismioues dans une aire de grande étendue. En conséquence, toute paire de récepteurs de sources, ou toute combinaison de récepteurs de sources, le long de tout azimut désiré entre ces récepteurs, peut être désignée par l'expression "ligne de profil". -4- L';ire d'explor-tion peut être aivisée en cellules ie en cellules ne forme géométrioue hppropriÉe telle ci)'une forme rectangulsire. Les idmensions de 1 cellule sont en fonction ae la resolution fondamentale désirée, mais des aimensions telles cue 25 par 50 mètres ou moins ne sont pas inhabituelles. Lors du traitement, pour toute paire donnee de récepteurs de sources, un point médian peut être désirné; le point médikn est affecté à une cellule appro- priee, pour cartographier et dans d'autres buts. Les brevets U.S. 4.241.429; 4.330.873; et 4.404.312 sont des divulgations typicues ies systèmes d'exploration sismioue 3-D (tri-dimensionnels) de la situation technioue antérieure. F.K. Levin, dans l'article "Apparent Velocity from DippinF Interface Reflectionsl" (Vitesse Appprente à partir de Réflexions d'InterfFce de' Dépression), GeoDhv- sics, 1971, pages 510 - 516, examine la vitesse moyenne en fonction de la dépression et de l'azimut. Ordinirement, le géophysicien n'a pas une connaisscnce antérieure de lF dépression locale ou de l'azimut local d'une ligne de profil par rapport à le direction de 1P dé- pression maximum. En conséquence, il existe toujours le pro- blème de l sélection de la fonction de vitesse de migra- tion appropriée pour la migration dépressive. Or ce pro- blème est d'un intérêt essentiel en cFs de l'existence de dépressions Pbruptes, sur un angle substantiellement plus important Que 10 degrés environ. Ainsi que Levin l'a fait remurquer, et comme le ouestion sera examinée plus loin, si l1 vitesse de migration n'est pas corrigée Fussi bien pour 1P dépression que pour l'azimut, lorsqu'on recueille des données de points de réflexion à partir de différentes déviations et de différents azimuts, les données sont macu- lées ou brouillées. -5 - La méthode l1 plus directe pour obvier à le muculature au point de réflexion consiste à assurer 1 mi-ration avant le recueil. Cette technioue convertit les ionnées de éviation non-?éro directement sur l'image de migra- tion finale, s ns aucune restriction en ce cui concerne le dépression, la déviation ou l'azimut. Toutefois, la migra- tion totale avent le recueil est une opération très cod- teuse au point de vue informaticue, étant donné oue pour une exploration 3-D (tri-dimensionnelle), jusqu'à 2 à 10 millions de traces, incluant chacune de 1500 à 2000 échan- tillons, doivent être traitées individuellement. En outre, il est nécessaire de connaitre de façon extrêmement pré- cise le champ des vitesses de migration, avant d'avoir lE possibilité d'une migration de pré-recueil. Cette condi- tion de rigueur implicue d'onéreuses analyses additionnel- les de vitesse. Il est évident, par conséquent, nue la fa- culté de compression des données offerte par le recueil avant la migration est fortement désirable. Une variante séduisante de la migration totale avant le recueil dans l'exploration 2-D (bi-dimensionnelle) est l'utilisation de la "Pre-stack Partial Migration" (Migra- tion Partielle de Pré-recueil), connue également sous la désignation "Dip Moveout Correction" (Correction de Sor- tie de Dépression) (DM0). Ce mode de traitement s'attaoue à la tache relativement modeste consistent à convertir sur une déviation zéro réelle des données de déviation finie ayant été soumises A une correction de sortie normple conventionnelle (NMO) appliquée avec utilisation de vites- ses de migration. La correction DM0 supprime l'effet de la maculature de point de réflexion, et permet de recueillir les circonstances d'une manière cohérente, indépendamment de la dépression ou de l'azimut. Les données peuvent en- suite etre recueillies et mises en migratiorn o façon tradi- tionnelle, mais avec un aegre considérable de compression le ces données. Un avantage supplementrire consiste en ce eue la correction DIi0 est une correction relativement petite cui a tendance à être insensible Pux erreurs d ns lR vitesse es- timée pour la correction NIO0. Hale a proposé un algorithme Pour la réduction de données 2-D (bi-dimensionnelle) dfns sa thèse de Doctorat "Dip Move- out by Fourier Transform" (Sortie de Dépression par Trans- formation Fourier), présentée au Département de Géophysicue de l'Université de Stanford, IMii 1983. Dans les pages 67 - 71 de cette thèse, il propose d'etenore son FnPlyse 2-D de 2-D à 3-D (tri-dimensionnelle), rais il échoue dFns sa ten- tative pour arriver à une solution analytique pratique au point de vue informatique. Un inconvénient majeur de la transformation de Hrle en 3-D est que celleci ne peut pas être refornulée directement comme une transformation Fourier. Sa méthode exige l'effort de calcul informatique équivalent de celui d'une trPnsforme- tion discontinue, et non un effort de Fast Fourier Transform (Transformation Fourier accelerée). Dans cette façon de pro- céder, la quantité de calcul informatique est la même pour toutes les dérivations (à l'exception de le dérivation zéro), et pour toutes les durées de translation, Elors qu'il est bien connu cue l'opérateur DM0 exerce peu d'effet pour les petites dérivations et les durées de translation croissantes. L transforwemation de Hale n'est pas en mesure d'exvloiter cette propriété. Un problème insoluble dans le traitement 3-D DMO0 consiste en ce eue la dérivation réceptrice de source est un vecteur au lieu d'être une quFntite échelle comme ceci est le ceps dans le traitement 2-D. En outre, les données sont brouil- ? - lEes ou maculées, sur une surface ellipsoale. eu lieu ie s'étendre le long d'une ligne. D-ns une exploration 3-D (tri-nirensionnelle) t.picue, o les amplitules ces ocrivations et des a7imuts peuvent varier très consicéra- 5.lement, le nombre des dérivations séparées qui aoivent être traitées en utilisant oes mores àe résolution trpdi- tionnels peut se rapprocher ou nombre des traces orns l'exploraiion. Le but de la présente invention consiste à Pssurer 18 mise en oeuvre efficace, et économique sur le plan infor- mptioue, avec assistance par ordinateur, d'une mèthode pour la conversion d'un ensemble ce traces de champ sismicue dis- tribuées en séries en un nouveau volume tri-dimensionnel clairement résolu de traces sismicues représentatives de la terre en présence de couches oepressives le long de derivF- tions Ayant des azimuts variables. Selon un rspect de la présente invention, une plurlité de trPces sismioues sont enrendrées sur une plurslité de récepteurs de signaux sismioues oui sont isonifiés par une pluralité de sources sismiaues disposées en déploiements ronux dans une région à explorer. Une ligne est établie entre une source sélectionnée et un récepteur, et une plu- ralité de cellules mi sont définies le long de cette li- mne. La trFce sismicue dérivée du récepteur sélectionné est corrigée en sortie normale et pondérée en raison in- verse de la dérivation entre un point médian m et la o cellule respective mi et proportionnellement à la racine carrée du temps de trace sismioue pour la création d'un ensemble ce traces pondérées. Les traces sont modifiées par application d'un coefficient de sortie. Les traces ro- difiées sont recueillies esans les cellules i correspon- azntes. Le processus est répète pour tous les erplacerents récepteurs de sources susceptibles de présenter ie l'inté- rêt. Les traces recueillies à partir de toutes les cellules peuvent être affichées en tant que représentation structurale tridimensionnelle de dérivation zéro, sens migration, a'un volume de la terre. Les buts et les methoces ce la présente invention se cor- prennront mieux en se reportant aux dessins, Fns leseuels: la Figure 1 définit les pistes de rayonnement associées avec une source acoustique et un récepteur en présence d'une couche ce terre non dépressive; la Figure 2 montre les effets a'une maculature ce point de réflexion en présence a'une couche de terre dépressive; la Figure 3 illustre un exemple d'un déploiement sismi- cue de récepteurs de sources; la Figure 4 est une vue isométrioue d'un volume de la terre destinée à illustrer la signification. de différentes quantités symbolicues; la Figure 5 montre une série de réflexions synthétiaues planes avec manifestations sismioues dépressives; la Figure 6 montre le même série de réflexions synthéti- cues que celle représentée à la Figure 5, après application d'un algorithme anti-adultérant; et la Figure 7 est un schéma séquentiel a'une vise en oeu- vre d'ora ina teur pour l'application de la méthode de la pré- sente invention. On admettra cue les spécialistes de la partie sort suffi- samment au courant des brevets '828 et '906 ainsi eue des autres éléments de documentation de référence mentionnés. précédemment, pour cu'il ne soit pas nécessaire d'examiner davantage ici les profils CIP (Point Médian Commun). - 9 - Les paragraphes oui suivent son: essentiellement oe na- ture didactique. En ce oui concerne maintenant la Figure 1, nous Définirons certaines cuantités cui seront examinées plus loin. Une partie de la surface de la terre est montrée en 10. Une source sonore S peut être dérivée à partir d'un récepteur R pour une distance de dérivation 2k, avec un point médian en mo. Un champ a'onde engendré en S, peut courir le long du chemin de rayonnement 12 jusou'à un point de réflexion 14 sur le réflecteur 16, en étant réfléchi en retour le long du chemin ou piste de rayonnement 18 jusou' au récepteur R. L'amplitude du champ d'onde en fonction du temps, f(t). peut être enregistrée sur un enregistrement à graduation chronologioue, mentionné ci-après sous 1E dé- signation "trace'. Les données sismiques sont généralement ouantifiees sous forme d'échantillons digitaux. lans les processus devant être examinés ci-après, on peut aeir in- dividuellement sur chaque échantillon de trace sismioue. Dans un but àe concision, l'emploi ce l'expression "trace sismioue" en liaison avec un processus signifiera que chaque échantillon de donnée de cette trace s été soumis individuellement au processus considéré. En raison de la distance de dérivation 2k, le temps d'arrivée t d'une ondulette particulière le long du che- min ou piste S-14-R est plus important aue le temps de cheminement tn d'une ondulette oui aurait cheminé le long d'une piste directe ou chemin direct mo-14-mo. La diffé- rence de temps est appelée sortie normale ou écart de che- minement normal (NM0). La quantité ft) meut être corrigée pour NM0 par la relation - 10 - n=t- (4K / V) (1) o V représente la moyenne oe racine carrée (RMS) au média par lequel se déplaçait l'ondulette. La relation ci-avant est exacte pour des composants dépression- réflecteurs d'approximativement cix degrés ou moins. Il convient de faire observer que pour aes dépressions insignifiantes et une dérivation zéro, les points de réflexion pour f(tn) tels cue les points 14 se situent le long c'un demi-cercle tel sue le demi-cercle 19 cen- tre autour de m.. En raison de la dispersion géométri- aue, l'amplitude o'ondulette varie en fonction de la loi des carrés inverses. Lorsou'il y a présence de dépression importante, la formulation pure et simple de l'Equation (1) n'est plus valide. Pour en venir maintenant à la Figure 2, une on- dulette peut se propager à partir ae S le long de la piste ou chemin de rayonnement 20 jusqu'au point de ré- flexion 22 sur le réflecteur dépressif 24 et de là jus- ou'à R le long ae la piste ou chemin 26, R étant dérivé à partir de S pour une distance 2k' Les points de ré- flexion pour une dérivation finie ne se situeront plus alors le long d'un demi-cercle, comme ceci est le cas pour des pistes ou chemins de déplacement de dérivation zéro; ces points de réflexion se situeront le long d'une ellipse 28 ayant ses foyers en S et en R4 Ppr conséquent, pour la piste ou chemin de rayonnement de dérivation fi- nie, le point de réflexion 22 véritable se situe en haut de dépression à partir du point de réflexion de dériva- tion réro 30. Le temps de parcours de dérivation finie - il - corripé pour NO0 (sarti= normale ou ecart ce cheminement normal) résulte enr une qurntité F(tn).u psint 30,- proje- tée le lonp d'une perpenriiculaire par r-pport à mi, qui est effectivement le point de profondeur s-ns migration pour une piste ou cr; erin de parcours re derivation ?ero. Le point de profondeur re migration réelle f(to) en 22 se o situe en aesous de ricui est séparé ce mj ppr l- distance Ax sec. En conséquence, si on essaie ae recueillir aes échantillons de traces ayant le même point médiat.Tais des dérivations différentes, les données sont brouillées ou ma- culees. Dans le cas o l'on aurait une connaissance préala- ble de sla dépression locale, on pourrait appliquer le fac- teur NiO en utilisant une vitesse dépencante ce la dépres- sion, mais pour des raisons qui seront expliquées, cette mise au point n'est pas toujours possible habituellement, et même ouand elle est possible, elle n'est correcte oue pour la seule dépression considérée. Nos remarcues introductoires ont concerné 1l géométrie bi-dimensionnelle le long d'un plan parallèle à la direc- tion de notre dépression. Toutefois, le but que nous-re- cherchons est la mise en image d'un volume tri-dimensionnel de ls terre. La Figure 3 est une vue en plan c'une disposition 34 de récepteurs (représentés par des cercles) mise en place selon un ces nombreux schémas oui peuvent être désirées. à la surface de la terre ou à proximité de cette surface. Le disposition peut être adoptée pour un grand nombre d'opérations sur terre ou en mer. Uniquement dans le but du présent examen, on admet cue le sens de progression le long d'une ligne de profil s'étp- blit de l'ouest à l'est comme par exemple le long de la ligne 36, parallèlement A la direction de dépression mpxi- - 12 - mum, comme montré pr le symbole ce tracé de dépression 38, ayant un angle de dépression %f. Une source S, à le position initiale représentée par le triangle renversé , est allumee de manière à produire aes traces sismi- ques à graduation chronoloFique pour chaoue récepteur dpns la disposition situé dOns la portée oe la source S. La source est ensuite déplacée sur une nouvelle position telle que la position de récepteur Ri, o la source est de nouveau allumée. L'exploration se poursuit au fur et à mesure oue la source est déplacée le long de la ligne de profil 36, puis réciproquement d'avant ern arrière le long de chacue rangée de récepteurs. Les lignes de ran- gée peuvent être explorées le long de 1P direction Nord- Sud comme on le désire. Chaoue fois cu'une source est al- lumée, une nouvelle trace sismique est engendrée pour chaque récepteur de la disposition. Le groupe de données réception par choaque récepteur peut être assigné à une cellule telle oue la cellule 42 associée avec les récep- teurs. Pour le traitement ultérieur des données, les cellules respectives peuvent correspondre à des emplace- ments appropriés dans une mémoire d'ordinateur ou dans une mémoire à disque, dans laauelle les traces sismicues peuvent être mémorisées sous forme d'échantillons digi- taux. Le long de toute ligne donnée telle cue la ligne 36, une image 2-D (bi-dimensionnelle), d'une portion de la terre, peut être produite comme il a été décrit pour la Figure 2. Les points de rayonnement incident pour les dé- rivations non-zéro se situent le long d'une ellipse. Pour toute ligne 44 se situant sur un angle arbitraire (azi- mut) par rapport à la direction ce dépression mpximum - 13 - dans une disposition 3-D (tri-dimensionnelle), les points de rayonnement incident s'établissent sur la surface d'un ellipsoïde dont le grpno axe est aligné le long de l'angle d'azimut e. Les quantités de Fraduation dans une projection 2-D (bi-dimensionnelle) deviennent des vecteurs en 3- D. Levin, oui a été cité précédemment, montre cue le rFpport entre une vitesse dépendant de la dépression Vd et 1l vi- tesse moyenne V est donnee par l'expression Vd/V = (1 - sin2 e cos t)}T (2) Pour des dépressions importantes et de petits Pzimuts, le rapport devient très important. Par conséquent, en ce oui concerne le traitement 3-D, la vitesse de recueil est aussi bien en fonction de la dépression Qu'en fonction de l'azimut. Nous ne connaissons pas habituellement la valeur de la dépression locale et du tracé local, bien oue nous puissions avoir une certaine connaissance de la dépression régionale du secteur considéré dans son ensemble.. De ce fait, nous ne connaissons pas la vitesse de recueil cor- recte. Or le but de la présente invention consiste à pré- voir une méthode de traitement de données oui soit indé- pendante de la dépression et de l'azimut, et qui réduise ainsi îa nécessité de se baser sur une connaissance précise de la fonction de vitesse. Reportons-nous maintenant à la Figure 4. Cette Figure représente la surface de 1a terre en coordonnées x, y. La Figure 4 est une vue isométrique éclatée de la partie de la disposition 34 qui inclue S, R1, R2 et Ri. La ligne 44 décrit un Fngle d'azimut de 8 par rapport à la direc- tion de dépression maximum. Une partie d'interface réflé- - 14 - chisspnte 46 est montrée en dessous ie la surface -vec une depression de y, comme dans les Figures 2 et 3. En dessous du segment de ligne 48 projeté verticalement sous la ligne 44, le réflecteur 46 F une composante de réflexion de &. La source S, ayant des coordonnées (xS, yS, est dérivée à partir du récepteur Ri, ayant des coordonnées (XR, R) sur une distence L = 2K et ay. nt un point méedian mo. La source S engendre un champ d'onde f(t) qui est réflécni à partir de l'interface dépressive 46 et qui est réceptionné par le rece- veur Ri. Nous définissons maintenant un plan ?f le long de la ligne 44, perpendiculaire à la surface et ayant des coor- données d'extrémité (us YS), (xR RYR)' Nous choisissons maintenant une distance iDf, et nous définissons des points mi sur la ligne 44 de façon que (mo, mi) = d(mo, mi) = iDf. (3) - Il est possible ou'une pluralité de points mi soit pro- duite de manière oue l'index mobile i se situe a-ns les li- mites -n oui est inférieur ou égal à i lequel est inférieur ou égal à n, nDf étant inférieur ou egal à k. Les légendes * mi désignent les centres de cellules telles que 42 à la Fi- i gure 3, ou des subdivisions plus petites de ces cellules. Les dimensions des cellules dépendent de la résolution dis- ponible à partir de la disposition 3-D (tri-dimensionelle=, et de celle désirée pour l'affichage final 3-D des données. En première étape, nous appliquons NM0 (sortie normple ou écart de cheminement normal) à f(t) afin de trouver f(tn). Pour cette première étape, la correction NM0 est ba- sée sur la moyenne des racines carrées de vitesse. L'étape suivante consiste à définir une ouantité Ci pour chaque va- leur i supérieure ou égale à -n, mais inférieure ou égale à n. telle qell - 15 - C i=- (iDf/k)2]. 4) On calcule ensuite un coefficient ponréré d'ampli-uie W, ve C wi(iDf,tn,2k) = Df/[(1+4k(2s-s2)) x[ (1-(iDf)2/(k2))33k (5) o s est égal à tn(s/f),.et F est la fréquence dominPnte. La fonction poncérée est asymptotiouement correcte comme approximation de éaerivation et de temps soit pour zéro, soit pour l'infini. Elle fait apparaitre un carect4re similaire à celui de la fonction pondérée (w) divulguée dans i0 la Demande de Brevet U.S. 06/779.176, mentionnée ci-avant, dans le zone dans leouelle w était acceptable. Les résultats analytioues laissent pronostiquer une bonne concordance avec l'opérateur de DMO (correction de sortie ou d'écart de dé- pression), f-k, pour toutes les valeurs de temps, de déri- vation et d'écerte7ent de traces. Ceci est confirmé par les résultats expérimentaux. Ensuite, les trpces Donderées f (tn) deviennent égales par définition wif(tn), modifié à le valeur mi respective 2. n i par application d'un coefficient de sortie ou d'écart de cheminement Ci qui est une fonction du rapport entre mo et la valeur de sortie respectivemi,etune moitié de la distance de dérivation (k), selon l'expression to = Ci x tn (6) et recueillent le résultat dans les traces sismiques conte- nues dans les cellules centrées à chacune des valeurs mi dans le plan p f. la quantité recueillie (totalisée) cor- respond A une expérimentation dérivation zéro dans l'espace tridimensionnel. - 16 - Les carictéristioues au coefficient ae sortie ou d'écart de cheminement Ci ne mpnquent pas o'interât. upnd iDF est zéro, c'est à dire quand mi. mo' C. est l'unité et f(to) f(tn), le temos de parcours à deux voies corrigé à partir de mo, en migration sur mj, par correction YIO (sortie nor- male ou écart de cheminement normal). Quand iDf se rapproche de k, c se rapproche de zéro, et f(to) se rpproche de zéro. Ce résult-t n'est naturellement pas plausible. En conséquence, ouand iDf se rapproche de k, les temps d'extrémités de l'el- lipse de point incident de réflexion sont modifiés de façon eue le grpdient ae temps soit limité en proportion de l'ex- pression: dt/dx peut être admise dans la zone d'exploration, et x est la dis- tance entre m. adjacents. i Après l'application, selon l'Ecuation 6, de la correction de sortie ou d'écart de dépression (DM0) à la trace corrigée NMO (écart de cheminement normal), de manière à obtenir la trace composante corrigée DM0, la teneur en fréouence locale de la trace doit être limitée afin d'éviter des mpnifesta- tions factices adultérentes. La fréquence permissible maxi- mum se détermine de préférence au moyen de l'expression: fma (to,iDf,Df(k 2-(iDf)2)/(2iDf x Df x to) (8) o Df est l'espacement des cellules le long de la ligne de source à récepteur S-mo-R. Reportons-nous maintenant à la Figure 5 qui illustre une section sismique synthétique. Les effets bénéfinues de la - 17 - caractéristique anti-'dultérante sont démontrés par une série de manifestations de réflexion planes 50 avec res manifestations de dépression 52 traversantes. Ceci est une section de dérivation commune à laouelle a été appli- ouee une DMO (correction de sortie ou d'écart de dépres- sion) trrditionnelle. On notera eue les manifestations planes et dépressives 50 et 52 ne sont pas une phase zéro, des précurseurs significatifs 54 apparoissnt en rvpnt de l'ondulette 56. L- Figure 6 représente la même section synthétioue après l'utilisation de l'algorithme anti- adultérant décrit ci-avant. Les précurseurs 54 ont été effacés dans une large mesure, et sont à peine discernables, ce oui donne une onoulette de phase zéro consistante 56. Aux fins d'explication, l'exemple de la Figure 4 8 été limité à un segment de ligne 44 s'étendant de S à Ri. Mais le traitement peut ne pas être borné à ce segment de ligne de courte longueur; il peut être développé de S à Rj et au-delà, comme montré à l1 Figure 3. En outre, d'autres paires de récepteurs de source peu- vent être choisies, comme par exemple à partir d'une source Si pr l'intermédiaire d'un récepteur Rk le long de la li- gne tiretée 45 dans la Figure 3. Certaines des cellules mi le long de Si - Rk coïncideraient alors avec d'autres cel- lules mj le long de la ligne S - Rj. Des trPces sismiques communes à toute cellule ouelconque, ouels oue soient les emplacements de récepteurs de source, seraient alors re- cueilliesensemble, ces traces sismiques représentant le même point de réflexion de sous-surface. Dans une explorr- tion de champ réelle, il peut y avoir autant d'emplacements de sources qu'il y a d'emplacements de récepteurs. Le pro- cessus exposé ci-avant s'exécute pour chaque combinaison - 18 - récepteur-source présent-nt de l'intérêt d.-ns l'rire n'e.- ploration. Usuellement, les valeurs recueillies finales d'ns les cellules respectives sont normalisses afin de compenser le nombre vpriable de traces totalisées. Les traces sismiques recueillies dnr.s chpcue cellule peuvent être affichées en volume tri-dimensionrel de la terre sous la forme, par exemple, d'un cube isométricue, ou bien les données peuvent être affichees en groupes de panneaux bi-dimensionnels, au gré de l'utilisateur, en employant une technioue de traitement traaitionnelle. De façon avantageuse, 18 procédure evite l'utilispion te vitesses dépendant de la dépression et dépendent de l'azimut. La seule vitesse implicuée est l'estimation la meilleure es- timation de la moyenne des carrés de vitesse employée pour le calcul du NM0 (écart de cheminement normal). En raison du nombre massif de calculs nécessaire, le processus est de préférence mis en oeuvre par calcul sur or- dinateur, selon le schema séquentiel de la Figure 7. La présente invention a été décrite avec un certain de- gré de spécificité à titre d'exemple uniquement. D'putres méthodes et techniques peuvent venir à l'esprit des spécip- listes de la partie, tout en restint dans le cadre de la matière divulguée par la présente description - 19 - REV.NDICATIONS 1. Une méthode oe recueil de donnees sismicues tri- dimensionnelles de manière à former une section sisricue, les donnees consistant en une pluralité oe traces corres- ponda.nt à des paires source-à-récepteur individuelles, re- cueillies dans une pluralité de cellules définies géométri- cuement sur une représentation de la surface de la terre, le recueil des données étant basé sur une moyenne de carré de vitesse sélectionnee, la dite méthode comprenant les phases suivantes: (a) correction d'écart de cheminement normal de cheoue trace préalablement 8u recueil ce la trace, la correction étant basée sur la distance de source à récepteur pour cha- oue trace; (b) la pondération de chaque trace corrigée par correc- tion d'écart ce cheminement normal de façorn à compenser la dispersion géométrioue entre les paires de source-à- récepteur, selon l'expression: wi(iDf,tn,2k) = Df/f (1+4k(2s-s2)) x / (l-(iLf)2/(k))3]. dns laouelle: wi est un coefficient de pondération d'amplitude; iDf est une moitié de la dimension de la cellule; tn est un temps corrigé par correction d'écart de chemi- nement normal; k est une moitié de 1P distance de dérivation; Df est l'espacement de la cellule entre la dite pluralité de cellules le long d'une ligne incluent les paires de source à récepteur; s est égal à tn(2/F); et F est égal à la fréquence dominante; - 20 - (c) 1F Molis4ication oes trHces corrigées par correction d'écart normal et pondérées, par un coefficient ae correc- tion d'écprt oui est en fonction A'une distance relative entre la dite pluralité ae cellules d.ns laouelle chacue trace doit être recueillie et le point médian entre la source et le receveur correspondent à cette trFce, l> modi- fication des traces étant effectuée selon les.expreseions: t = Ci x t; et l o i n C. = (l-(iiDf/K chaque trace corrigée par correction d'écart normal de ma- nibre à éviter des mnnifestations factices adultérantes dans les données sismiques recueillies, selon l'expres- sion: FMX(to,iDf,Df) = (k' - (iDf)2)/(2iDf x Df x to), dans laquelle: Fmx est la fréquence permissible maximum. 2. Dans une méthode de recueil de données sismioues tri- dimensionnelles, les données consistant en une pluralité de traces correspondant à une source individuelle pour des paires de récepteurs et ayant une gamme de fréouence pré- déterminée, les dites données étent recueillies dans une pluralité de cellules géométriquement définies sur une re- présentation de la surface de la terre, et le recueil des données étant basé sur une moyenne de carré de vitesse sélectionnée, le perfectionnement comprenant; (a.) préalablement au recueil des traces, le correction - 21 - d'écart normial ne chaque trace b&sée sur 1l distance source à récepteur pour chaque tr-ce; (b) la pondération de cbaoue trace corrigée par eorrec- tion d'écart norrmal Pfin de compenser la oispersion réoré- trioue entre la source et les paires ae recepteurs; (c) la modification des traces corrigees par correction d'écart nor-al et pondérées obtenues pzr un coefficient d'écart oui est en fonction d'une distFnce relative entre la cellule cans laquelle chaQue trace doit être recueillie et le point médian entre la source et le récepteur corres- pondnnt à cette trace; et (d) 1 limitation de le teneur en fréauence de chacue trace de manière à éviter des manifestations artificielles adultérentes, selon l'expression: F ax(to,iDf,Df) = Sk - (iDf)2/(kiDf Df x t), dans laouelle: FmaX est la fréquence maximum permissible; t est le temps de parcours modifié du dignPl sisrioue; iDf est une moitié de la dimension de l1 cellule; Df est l'espacement de la cellule entre la source et les paires de récepteurs; et k est une moitié de la distance de dérivation. 3. La méthode définie selon la revendication 1, caractérisée en ce qu'elle inclut en outre la phase consistant à définir un gradient de temps entre des cellules adjacentes en propor- tion d'une dépression géologique maximum à prévoir. 4. La méthode définie selon la revendication 2, caractérisée en ce qu'elle inclut en outre la phase consistant à définir un gradient de temps entre des cellules adjacentes en propor- tion d'une dépression géologique maximum à prévoir.