FR 2154707 A1 19730511 FR 7234419 A 19720928 La présente invention concerne un procédé destiné à accentuer les amplitudes des diffractions élémentaires ou événements sismiques et/ou à les éliminer en revenant à leur source. Dans la pratique oourante de sommation des informations sismiques en provenanoe d'une oote de référence (PCC) à la suite d'un tirs les informations sont recherchées et sommées le long d'hyperboles, oe qui permet de déterminer au mieux la vitesse à des instants de référence échelonnés le long d'un groupe donné de traces sismiques. Ce procédé est basé sur l'typo thèse que des réflexions provenant d'un miroir commun doivent apparaftre sur une hyperbole particulière b travers les traces. Cependant, on a observé que stil y a une faille ou un plissement le long des intersurfaces, la relation temps-distanoe entre les événements sismiques et les traces pour un PCC ntest pas hyperbolique. Les franges de diffraction provenant d'une telle struoture sont particulièrement intéressantes. Etant donné que les événements engendrant d'une autre manière ces franges ne se trouvent pas sur des hyperboles ooupant les traces pour un PCC, une reoherohe et une sommation le long de oes hyperboles provoquent la sommation d'un grand nombre de ces éirénements avec un décalage de phase et ils sont ainsi "détruits" ou éliminés.Par conséquent, la structure recherchée, qui est å l'origine des franges de diffraction, peut ne pas être reconnaissable. La présente invention concerne un procédé de sommation des traces de PCC ou des traces d'un profil à alignement vertical, qui est mis en oeuvre dans deux oas sur des diffractions élémentaires. Ce procédé est mis en oeuvre aveo un équipement de calcul automatique sur chaque trace sismique à sommer et détermine le long de chaque trace les instants où l'enregistrement d'ondes sismiques est attendu et les instants d'émis- sien a partir de points source bypothétiques a des profondeurs au-dessous d'un point de référence en surface, écarté de la trace sommée.Des éohantillons de trace sont, à ces instants, sommés de manière a former un ensemble sommé à dea instants qui sont fonction du temps de propagation sol-miroirol entre les points source hypothétiques et le point en surface correspondant à ceux de la trace sommée. Ces opérations sont mises en oeuvre pour chaque trace d'un groupe de traces sismiques différentes. Dans le mode d'exécution où les traces à sommer sont des traces de PCCI le procédé est mis en oeuvre en augmentant de quantités finies différents temps correspondant a des points au-dessous du point de réfé- renoe en surface. Pour un groupe donné de traces de PCC, on caloule le temps (Le propagation sol-miroir-sol des ondes sismiques à partir de chaque instant échelonné en direction du point de surface de référence d'un groupe donné. On calcule aussi les temps de propagation des ondes sismiques entre chaque couple "point de tir-récepteur de l'ensemble" sur le terrain et chaque instant échelonné.Des échantillons de chaque trace sont ohoisis à des instants correspondant au temps de propagation point de tir-réoopteur calculé. Ces échantillons de trace sont sommés pour former un ensemble de temps correspondant au temps de propagation associé sol-miroir-sol associé. Ces opérations sont mises en oeuvre aveo plusieurs groupes de traces en se basant sur des calculs à partir de temps éohelonnéXF au-dessous dudit point de référence en surface. Quand la sommation est effectuée sur des traces sommées ante- rieurementl il existe deux modes d'exécution. Dans l'un, le temps est augmenté de quantitéa tnies(d onnement) le long d'une trace de référence correspondant au point de référence en surface et la vitesse en sous-sol pour chaque instant échelonné sert à calculer le temps de propagation sol-miroir-sol à partir de chaque temps échelonné jusqulå la surface correspondant å la trace sommée. Les temps de propagation sol-miroir-sol calculés sont employés pour choisir des échantillons le long des traces. Ces échantillons de trace sont sommés pour former un ensemble sommé aux instants correspondant aux temps échelonnés au-dessous du point de référenoe en surfaoe. Dans un autre mode d'etéoutions le temps est augmenté de quants tés finies le long de chaque trace à additionner et la vitesse au-dessous de la surface pour chaque temps échelonné le long de oes traces sert à calculer les temps de propagation tertioaus sol-miroir-sol au-dessous des points de référence en surface qui, par hypothèse, sont les sources d'énergie diffractée. Les échantillons de trace pour chaque temps 60he- lonné le long de chaque trace à sommer sont sommés pour former un ensemble sommé à des instants correspondant au temps de propagation vertioal solmiroir-sol calculé autour du point de référence en surface qui > par hypothèse, est la source de l'énergie diffractée. On fait remarquer que l'utilisation nouvelle d'un système de type connu de traitement de l'informations utilisation définie par les programmes exprimant en langage symbolique les divemsopérationa que comprend le procédé de l'invention > a pour effet technique de réaliser concrètement une machine nouvelle à caractère industriel certain. D'autres objets et avantages de la présente invention seront mieux compris à la leoture de la desoription qui va suivre de plusieurs exemples de réalisation et en se référant aux dessins annexés dans lesquelss les figures 1A et 13 représentent une teohnique pour obtenir sur le terrain des informations oonoernant une oote de référence (PCC) ou un point en surface de référence (PSC); la figure 2 représente une technique de recherche par hyperbole mise en oeuvre le long d'un groupe de traces sismiques pour des PCC; la figure 3 est une courbe représentant une grandeur obtenue par une teohnique de recherche par hyperboles, mise en oeuvre à un instant T0 donné;; la figure 4 est une courbe représentant une grandeur obtenue par un procédé de recherche par hyperboles ainsi que de détection, mis en oeuvre sur l'enregistrement complet d'un groupe donné de traces de PCC; la figure 5 représente l'ensemble des grandeurs obtenues par une technique classique de recherche par réflexion et sommation; la figure 6 représente l'ensemble des grandeurs obtenues par une teohnique de recherche par diffraction et sommation selon l'invention;; la figure 7 est un profil obtenu par addition du profil de la figure 5 au profil de la figure 63 la figure 8 représente la manière dont les calculs sont effeotué. pour mettre en oeuvre l'organigramme nO 13 les figures 9 et 10 représentent le mode de sommation de 1 'in- formation lors de la mise en oeuvre du procédé de l'organigramme nO 1 et de la figure 8g la figure 11 représente les calculs à partir d'une droite sismique jusqu'a une ligne de faille connue;; les figures 12 et 13 représentent un mode d'exécution pour la "compression" des diffractions élémentaires par traitement d'un profil somié; les figures 14 et 15 représentent un autre procédé de compres- sion des diffractions élémentaires;; la figure 16 représente les extrémités accentuées d'une diffraotion élémentaire qui apparaissent quand on oomprime les diffractions élémentaires en remontant jusqu'à leur source par le procédé des figures 12 b 153 la figure lÂ représente un procédé d'obtention d'une série de traces sismiques qui est destiné à mettre en oeuvre un tir aveo oote de référence, parfois dénommé aussi "tir avec point de surface de référence". La couverture représentée est une couverture d'ordre deux très simplifiée, autrement dit il y a deux traces sismiques pour chaque PCC. Les ondes sismiques sont émises aux points de tir S1 et S2 et l'énergie réfléchie par une intersurface réfléchissante en un point à cote de référenoe(PCC)est détectée par des récepteurs en E et R2 pour la produotion d'une série de traces sismiques. Ces traces sont de préférence enregistrées sous forme numérique sur bande magnétique.Les emplaoements de tir et d'enregistrement sont symétriques par rapport au point de surface de référenoe(PS; Dans un exemple, la couverture de la figure 1A est obtenue en utilisant une source S et un groupe de trois récepteurs R en un endroit et en déplaçant ensuite la source et l'ensemble de récepteurs vers la gauche pour opérer de la manière représentée sur la figure 13. On obtient N groupes de traces de PCC.On se reportera à l'article " Common Refleotion Point Horizontal Data Starking Techniques" par W,Harry Marne GS0PSICS, Vol. tXVII No. 6, Il (Décembre, 1962) pp.927-938, 8 figures, pour une description d'autres appareillages et techniques d'obtention des ootes de référence. Si les informations sismiques étaient obtenues par réflexion sur des miroirs plans souterrains superposés, les réflexions pour un groupe de traces sismiques de PCC, à partir d'un miroir oommun, tendraient å se placer sur une hyperbole particulière. Dans le mode de traitement classi- que, les traces pour des PCC sont recherchées pour des temps échelonnés (les T0 de la figure 2) à travers un groupe d'hyperboles (E sur la figure 2) ou de vitesses pour ohoisir à chaque instant l'hyperbole ou la vitesse pour laquelle les réflexions élémentaires d'un miroir commun sont le mieux alignées.Les informations des traces provenant de chaque groupe de PCC sont sommées le long des hyperboles optimales pour former des ensembles sommés destinés à l'établissement d'un profil à alignement vertical. L'hyperbole définissant les réflexions élémentaires sur les traces de PCC est définie par l'équation: dans laquelle: T est le temps de propagation source-réoepteur des ondes sismiques passant par une intersurfaoe souterraineS T est le temps, avec un déport horizontal nul, de propagation sol miroir-sol entre la surface et l'intersurface pour une distance point de tir-réoepteur nulle, X est la distance horizontale souroe-réoepteur Va est la vitesse efficace ou vitesse moyenne quadratique. La figure 3 représente l'amplitude ou énergie d'un signal détecté pour chaque vitesse admise pour un T0 donné obtenue par la technique de recherche par hyperboles. Le maximum de la trace représente la détermination optimale de vitesse pour le temps T01. Des traces du type représenté sur la figure 4 sont produites par le procédé de recherche et de deteo- tion par hyperboles mis en oeuvre le long de l'enregistrement complet d'une série de traces pour des PCC. Les traces horizontales de la figure 4 représentent les amplitudes des signaux détectés en fonction de la vitesse pour chaque T0 échelonné tandis que la courbe 21 est le lieu des valeurs maximales ou "pios" de chacune des courbes d'amplitude de signaux.La oourbe 23 est la courbe des temps de propagation pour les réflexions primaires. Elle représente la meilleure détermination de vitesse des réflexions primaires le long de l'enregistrement en utilisant la teohnique de recherche par hyperboles. Les valeurs oorrespondantes de la vitesse sont destinées à la oorreotion d'obliquité et à la sommation des traces du groupe, ou des groupes de PCC. L'enregistrement de profil de la figure 5 est un profil å alignement vertioal obtenu par sommation de 160 groupes de résultats d'ordre 6 obtenus sur le terrain pour des PCC en utilisant la teohnique de reoherohe et sommation par hyperboles. On se reportera i ce sujet à l'article "Velo- oi4 Speotra - Digital Computer Derivation of Velooity Functions" de X.Turhan Taner et Fulton Koehler, présents au 47e congrès annuel de la "Society of Exploration Geophysicists" Oklahoma City, Etats-Unis d 'Âmérique, 2 novembre 1967, pp. 1 -16; appendice pp.la-7a, (18 figures), pour une étude plus poussée de oe type de traitement. Comme indiquez ci-dessus > on suppose dans le type de traitement précédent que les miroirs sont des miroirs plans superposés et par oonse- quent que les réflexions sont alignées le long d'hyperboles pour des enre gietrements de PCC. Cependant, cette hypothèse n'est pas valable dans les cas où des irrégularités marquées telles que des bords de failles existent dans les oouohes suuterraines. Comme décrit dans la littérature, ces bords ou d'autres irrégularités peuvent agir comme des sources ponctuelles le long de la droite sismique d'exploration et produire des phénomènes ou des franges de diffraction. Voir, par exemples 1 'ouvrage SEISMIC PROSPECTING FOR OIL, C. Hevitt Dix, Harper & Brotherss Noir York (1952). Ces franges sont intéressantes étant donné qu'elles facilitent l'identification des struotures souterraines; cependant, elles sont souvent difficiles à observer sur les informations traitées par des procédés classiques. Par exemple, pour les informations rassemblées pour des PCC ou des PSC les évXnements engendrant les franges de diffraction ne sont pas obligatoirement alignés sur des hyperboles. Par oonséquentw la reoherche et la sommation le long d'hyperboles a pour oonséquenoe qu'un grand nombre d'év8nements qui donneraient lieu dans le cas oontraire à des franges de diffraction sont sommés aveo un décalage de phase et disparaissent. La présente invention concerne selon une de ses caractéristiques, un procédé qui utilise les traces de PCC pour calculer et définir les temps d propagations ou les temps d'arrivée, pour des événements de diffraction provenant de points ou régions intéressants souterrains. Ces temps définissent des échantillons de traces qui sont sommés le long des courbes de temps de propagation calculés pour former des ensembles sommés en vue d'accentuer les évenementss de diffraction. Pour mettre en oeuvre le procédé selon l'invention on procède à une détermination initiale de la stratification et de la structure souterraines. Ces déterminations peuvent englober le profil de la figure 5 obtenu par les procédés classiques de recherche et de sommation par hyperboles. À ce sujet, les informations rassemblées sur place pour des POC ou PSC sont obtenues à des temps T0 échelonnés pour une gamme déterminée de vitesses en vue d'obtenir la détermination optimale de vitesse & BR Organigramne nO 1 - (voir page 24) 34) recherche de l'ajustement optimal de la vitesse. 35) correction dtobliquit et sommation. 36) interprétation du profil à alignement vertical pour détecter les sources ponctuelles. 37) calcul des courbes de temps de propagation pour l'énergie diffractée le long des traces pour des PCC et provenant de sources d'énergie diffractée, à partir des traces de référence. 38) sommation de plusieurs échantillons de traces pour chaque groupe de PCC le long de leur courbe de temps de propagation pour former des ensembles sommés afin d'augmenter l'amplitude des événements de diffraction provenant des sources spécifiées. Après le choix des traces de référence on calcule les courbes de tempe de propagation de I 'énergie diffractée le long des traces pour des PCC provenant de sources d'énergie diffraotée, à partir des traces de référence0 Ces calculs figurent dans le rectangle 37 de L'organigramme n 1. Les échantillons de traces pour chaque groupe de PCC sont sommés le long de leur courbe de temps de propagation pour former des ensembles sommés en vue d'aooentuer les évenements de diffraction a partir des sources spécifiées. Cette dernière opération est représentée en 38. Les opérations 36 37 et 38 sont décrites plus en détail en référeno aux figures 5 et 8. Lors de l'étude de l'enregistrement de profil de la figure 5, un événement intéressant est défini en 40. Cependant, on ignore si cet éveflement est dû b des diffractions, étant donné qu'il n'y a pas de franges de diffraction bien définies à partir des bords 40a et 40b. Ces bords sont, par conséquent, considérés comme des sources ou zones souterraines intéressantes. Les groupes de PCC employés pour obtenir le profil de la figure 5 sont traités et sommés pour augmenter l'amplitude des franges de diffraction provenant des sources 40a et 40b. Le profil å alignement vertioal aveo les franges de diffraction résultantes est reprd- senté sur la figure 6. La technique de traitement et de sommation employée pour obtenir le profil de la figure 6 est déorite oi-aprè. Sur la figure 5, les traces sommés 32 et 92 passent respectivement par les sources ponctuelles 40a et 40b, Ces traces servent de traces de référence en vue des calculs et des sommations pour obtenir le profil de la figure 6. Ils correspondent respectivement aux groupes de PCC n 32 et 92. Comme l'indique la figure 8, chaque groupe de traces de référence ou de PCC à partir duquel des calculs doivent être mis en oeuvre est désigné par LFT. Le groupe de PCC sur lequel on effectue des calculs et qui est sommé est désigné par NCDP. Sur la figure 8, un ECDP est représenté aveo deur sources ou points de tir S1 et S2 et deux récepteurs ou empla oements de réoepteurs Rî et R2. Cependant, il va de soi qu'il y aura, en pratique davantage de trajets point de tir-récepteur, par exemple 6 ou 12. Un point de tir S1 et un emplacement de récepteur RI forment, par défini tion,un couple point de tir-récepteur, S2 et R2 sont un autre couple point de tir-récepteur, eto. Des calculs effectués à partir d'un LFT donné pour obtenir un groupe donné de PCC afin de définir des échantillons de trace pour des sommations sont décrits oi-après. D'après la figure 8, pour un LFT, le temps TSUM est identifié, au départ, à zéro. Ce temps est augmenté, à chaque fois, le long du profil, d'une quantité delta représente par DELT, de manière que TSUM = TSUM + DELT pour chaque nouvel accroissement. Le nombre d'augmentations de TSUM est égal à NSA, le nombre d'éohantillons de trace. A chaque instant TSUM, pour un LFT donnés on caloule TA, TB et TKF pour une trace donnée. TA est le temps de propagation des ondes sismiques entre ltemplacement défini par TSUM et un emplacement donnés par exemple R1, de récepteur.TB est le temps de propagation de l'énergie sismique montre TSUM et S1, le point de tir associé à R1. . TFK est le temps de parool1rs sol-mircir-sol de l'énergie sismique entre PSC et TSUM > PSC représentant le point de surface de référence correspondant au PCC sommé. Lorsque TSUM est augmenté de DELT, la source ponctuelle 40a est inoluse sous forme d'un emplacement à partir duquel on effectue les calculs de TA, TB et TSK. La somme de TA et TB est identifiée à TSK qui est le temps de propagation entre S1, TSUM et E . Par conséquent: TSK r TA + TB (2) Sur la tracas enregistrée par un récepteur donné, du groupe de PCC sommé TSK localise l'échantillon de trace qui représente l'énergie reçue et enregistrée par le récepteur et se propageant entre le point de tir et le réoepteur associé en passant par le point souterrain TSUM correspondant. Voir les figures 9 et 10 sur lesquelles X1, x2 et X3 représentent des traces sismiques.Si l'on somme cet échantillon et donne un ensemble STKA à l'instant correspondant TKF, STKA correspond avec le point de surface de référence correspondant au groupe de PCC sommé. Pour chaque instant échelonné TKF, on calcule TSK pour une trace donnée d'un groupe donné de PCC pour définir un échantillon de trace sur cette trace et l'instant de sommation associé de l'ensemble sommé. On fait oela pour chaque trace du groupe de PCC sommé et les échantillons de trace définis par chaque TSK calculé sont sommés en donnant l'ensemble sommé à l'instant TKF approprié. En d'autres termes, les temps TSK calculés à partir d'un groupe de PCC récepteurs,en partant d'un TSUM donné,définissent une courbe de temps de propagation de l'énergie diffractée le long des traces sismiques correspondantes de PCC. Des courbes de ce genre sont représentées sur les figures 9 et 10 en TX1, TX2 et TX3.Des éohantillons de trace le long de ces courbes sont sommés pour former l'ensemble sommé à l'instant TKF approprié pour acoentuer les diffractions élémentaires. Après qu'un groupe de PCC a été sommé pour former l'ensemble sommé TSKA, il est enregistré et un autre groupe de PCC est sommé et enregistré en se basant sur des calculs à partir du même LFT. Ceoi eet fait pour plusieurs groupes de PCC. Ensuite, les calculs oommencent à partir. du seoond LFT et plusieurs groupes de PCC sont sommés et enregistrés pour aocentuer les phénomènes de diffraotion provenant de sources ponctuelles pour le seoond LFT. Pour l'obtention du profil de la figure 6, les groupes de PCC 1 à 62 ont été sommés et enregistrés en se basant sur des calculs å partir de LFT = 32. Les groupes de PCC 63 à 160 ont été ensuite sommés et enregistrés en se basant sur des calculs à partir de LFT . 92. Sur la figure 6, des franges de diffraction 40a' sont nettement visibles et partent de l'emplacement de la source ponctuelle 40a. Ceoi indique que ladite source 43 a été correctement ohoisie. Les franges de diffraction 40b' correspondant à la source ponctuelle 40 b ne sont pas très visibles, ce qui indique que la source ponctuelle 40b n'a peut-être pas été ohoisie oorreotement. La figure 7, qui représente l'addition des profils à alignement vertioal des figures 5 et 6 représente nettement les franges de diffraction 40a' partant de la source ponctuelle 40a. On peut en oonolure que l'événement 40 est dû à des diffractions. I1 va de soi que, au lieu de ohoisir des traces de référence déterminées à partir desquelles on effectue les calculs, ces traces de référence peuvent être échelonnées dans le profil. A partir de chaque trace de référence on peut effectuer des calculs pour déterminer chaque groupe de PCO, et chaque groupe de PCC est sommé d'après ces calculs. Ces traces sommées peuvent être mémorisées et la trace de référence dé placée et on reoommenoe oes opérations jusqu'à oe que la trace de réfê- rence ait balayé tout le profil. Toutes les traces mémorisées oorreepondant à la trace n 1 peuvent être sommées et imprimées; de momie, toutes les traces mémorisées correspondant à la trace n 2 peuvent être sommées et imprimées,etc. On peut établir de cette manière un profil somma par diffraction. On décrit oi-aprbs plus en détail en se reportant b nouveau å la figure 8, le mode de calcul de TÀ, TB et TKF. Les distances source récepteur de chaque groupe de PCC sont oonnues à partir des opérations sur le terrain. La distance source-réoepteur pour la trace traitée est désignée par XÂ. La distance horizontale entre LST et PSC est désignée par IF. On voit sur la figure 11 que la droite d'exploration sismique fait un angle déterminé par rapport a une faille connue. La perpendiculaire abaissée de la ligne de faille sur chaque PSC est employée comme IF, et l'on a IF = XF' sin #, p étant l'angle entre la ligne de faille LF et la ligne sismique. La distance entre les points de surface de référence est oonnue par des opérations sur le terrain. Ces distances sont désignées par XCDP. Le LFT sur laquelle on effeotue des sommations est oonnu- tout comme les PCC sommés.Par conséquent, on peut calculer IF par l'équation oi-aprèss IF - (DPT-NoDP)(XCDP) (3) Comme on l'a indiqué oi-dessus, le temps TSUM est augmenté de DELT pour vaque calcul de TA et TB. Dans un mode d'exécution, chaque augmentation du temps DELT est égale à 4 ms. Par conséquent, l'instant d'un TSU > déterminé est oonnu ou peut être calculé. Une courbe de temps de propagation ou une fonction de vitesse du type représenté en 23 sur la figure 4 est déterminée ou oonnue par des explorations et des traitements antérieure. La valeur de la vitesse est connue ≈partir de la courbe pour chaque TSUM.Par conséquent, la distance verticale Y entre un TSUM et la surface est déterminée par la relation Y = TSUM x V/2. V/2. (4) La distance horizontale entre un TSUM et S1 et un TSUM et HI peut Autre déterminée de la manière représentée sur la figure 8. Le temps de propagation TA entre HI et un TSUN est déterminé par l'équation oi-aprees Le temps de propagation entre la source SI et un TSUM est donné par l'équation ci-aprèss Le tempe de propagation sol-miroir-sol TKS entre un PCC et un TSUM peut Autre déterminé par l'équation ci-après: Les organigrammes 2, 3 et 4 décrivent des procédés d'obtention des profils des figures 5, 6 et 7.Ces organigrammes sont suffisants pour permettre å un spécialiste qualifié d'écrire un programme de calculatrice en langage FOR?RLX en vue de donner des instructions à une calculatrice numérique d'usage universel, par exemple du type CDC 6400, pour effectuer les opérations. Les organigrammes 2 et 3 concernent un programme principal intitulé DIFSTK pour traiter les informations provenant de PCC afin de mettre en oeuvre le procédé. Ils appellent un sous-programme GENDLT (organigramme n 4) qui, associé au programme principal, traite les informations de PCC pour produire le profil de la figure 6.Les informations traitées sont appliquées å un enregistreur de profil à surface variable pour obtenir les profils des figures 5, 6 et 7. Ci-après quelques définitions concernant les organigrammes: STK est l'ensemble sommé de réflexions alors que STKA est l'ensemble sommé de diffractions. NCD est un compteur tandis que NFOLD est l'ordre N de couverture des PCGs par exemple 6 ou 12. NC est un compteur d'échantillons de trace tandis que TSUbI est un compteur de temps. DELT est l'intervalle temporel entre les échantillons de traces. NSA est le nombre d'échantillons par trace tandis que NTROUT est le nombre de traces à obtenir à la sortie. Organigrme n 2 - (voir page 25) Programme DIFSTK 51) leoture des informations à l'entrée. 53) NCD = O 55) amener au zéro l'ensemble sommé par réflexions STK ainsi que 1 'ensem- ble sommé par diffractions STKÂ. 57) lire le temps, les informations de vitesse, Si nécessaire, au point de commande. 61) DO 996 M = 1, NFOLD 63) inscrire la trace à sommer. 65) localiser la distance point de tir-géophone pour l'ensemble des dis tances. 67) NO-O TSUM - O 69)69 TSUM = TSUM + DELT 71) NC " NC + 1 75) identifier et ajouter l'échantillon de trace défini par T pour obte nir l'ensemble sommé TSK. 77) IF (NC.EQ.NSA) 79,69 Sur l'organigramme2 ci-dessus, les informations à l'entrée sont lues comme indiqué en 51. Ces informations comprennent la distance point de tir-récepteur, l'espacement des PCC ou l'espacement entre les points en surface commun des groupes de PCC adjacents. Une fonction de la vitesse ou courbe de temps de propagation, qui est une fonction de la vitesse et du temps est également lue. Un compteur NCD est réglé sur zéro en 53. En 55 (instruction 55) des ensembles STK et STKA sont amenés au zéro. Les informations de vitesse sont lues à la demande au point de commande en 57. En 61, les instructions DO LOOP, DO 996, de 1 à NFOLD (ordre N) commandent la leoture successive de N traces; on détermine les instants de sommation correcte pour les ensembles sommés de réflexion et de diffraction et on somme les informations. La trace à sommer est inscrite en 63 et la dis- tance point de tir -géophone de l'ensemble des distances est définie en 6T, Les compteurs NC et TSUkI sont ramenés au zéro en 67, tandis que TSUM est augmenté de DELT en 69. L'indication du compteur NC est également augmentée de 1 en 71. L'équation de correotion du temps de propagation normal pour la sommation des réflexions est établie en 73.En 75, l'échantillon de trace défini par T est identifié et sommé pour former l'ensemble STK sommé. Er. - NC est vérifié pour déterminer s'il est égal à NSA. Si c'est Te cas, @@@@ les échantillons qui ont été lus de la trace traitée ont été sommés et a commande est transférée à l'instruction 79. Sinon, la commande est retransférée à l'inst:action 69 et un autre échantillon de la trace est éshartilloriné. Organigramme n0 3 - (voir page 26) 79) 79TSUM = O 81) appeler GENDLT. 83) effeotuer 83 J r SA TSUM n TSUM + DELT A partir de TSUM, localiser l'échantillon de trace approprié défini part et calculé par le sous-programme GENDLT et ajouter cet éohan tillon de trace à l'ensemble de diffraction sommé STKA en TSUM. 85) instruction continue 996. 89) enregistrer les traces sommées de réflexions sur le disque 5. 91) enregistrer les traces sommées de diffractions sur le disque 6. 93) NCD = NCD + 1. 95) IF (NcD. N. NTRGUT) 97) DO 103 N = 1 NCD 99) extraire une trace sommée de réflexions du disque 5 et mémoriser. 101) extraire des traces sommées de diffraotiom du du disque 6 et mémoriser. 103) ajouter les traces sommées de réflexions aux traces sommées de diffrac tions. 105) enregistrer la somme sur l'unité 3. 107) enregistrer les traces sommées de diffraotiom sur l'unité 3. 109) enregistrer les traces sommées de réflexions sur l'unité 3. Sur l'organigramme n 3 TSUM est identifié à zéro en 79 et le sous-programme GENDLT est appelé en 81, de manière à mettre en oeuvre la sommation des diffractions à l'aide d'un ensemble T calculé par le sousprogramme GENDLT. Dans ce cas, le sous-programme GENDLT calcule un inter valle # T à partir de chaque TKF pour la trace traitée. Chaque # T définit l'échantillon de trace approprié relatif à son TKF. Les accroissements A T sont mémorisés dans l'ensemble T. La sommation des diffractions pour une trace unique est mise en oeuvre par l'instruction DO LOOFA DO 83 figurant dans le reotangle 83.D'abord TSUM est augmenté de #T et ltéohan- tillon de trace approprié défini par #T et calculé par GENDLT est défini et ajouté à l'ensemble sommé des diffraotions STKB en TSUM. On fait cela pour tous les échantillons de trace NSÂ de la trace traitée. Après que tous les échantillons de la trace ont été correctement sommés, la oommande est retransférée à l'instruction DO 996 par l'instruction continue du rec- tangle 85. Une autre trace est captée et sommée. Cela est fait N fois (EFOLD fois) pour obtenir une trace de sortie pour un groupe donné de PCC. En 89, une traoe de réflexions sommées est enregistrée sur le disque 5 tandis qu'en 91 une trace de diffraction est enregistrée sur le disque 6. NCD est augmenté de 1 en 93 et on vérifie si le nombre désiré de traces NTROUT a été enregistré. Sinon, la oommande est retransferée à l'instruction 53 de manière que les ensembles sommés soient ramenés i zéro et len sommations de réflexions et de diffractions soient mises en oeuvre aveo un autre groupe de PCC pour produire une autre trace de sortie. En 97, PO LOOP, l'instruction DO 103 est mise en oeuvre de 1 i NCD foin. En 99, une trace sommée de réflexions est extraite du disque 5 et mémorisée tandis qu'une trace sommée de diffraction 101 est extraite du disque 6 et mémorisée. En 103, la trace sommée de réflexions est ajoutée b la trace sommée de diffractions et la oommande est retransférée au début de l'instruction DO LOOP, DO 103 de manière que les opérations des reotangles 99, 101 et 103 soient recommencées NCD fois. En 105, les sommes des traces sommées de réflexions et de diffractions sont enregistrées sur l'unité 3. En 107 et 109, les traces sommées de diffractions et de réflexions sont aussi respectivement enregistrées sur l'unité 3. On décrit oi-après l'organigramme 4 du sous-programme GENDLT. Organigramme n 4 - (voir page 27) Sous-pro gramme GENDLT 121) insorire les numéros des PCC (LFT à partir desquels on exécute les oalouls). 123) ramener à zéro ensemble T où les intervalles b T sont mémorisés. 125) attribuer un numéro à un ensemble fiotif de PCC au-delà de l'ensemble de PCC le plus à gauohe du profil. 127) mettre en place les points centraux entre les PCC correspondant aux numéros à partir desquels les calculs doivent être effectués. 129) si le numéro du PCC qui est sommé correspond à un point au-delà du milieu, commencer les calculs à partir du numéro suivant de PO, au dessous duquel l'interprétateur suppose la présence de source de phénomènes de diffraction. 131) TSUM n DELT NC n 133) DO 133, J " 1, NSA calculer TKF et TSK à l'aide des équations 2 à 7 T(NC)= TSK - TKF TSUM = TSUM + DELT 133 NC n NO + 1 133A) instructions de retour. Les numéros du PCC ou du LFT à partir desquels les calculs sont effectués sont inscrite en 121. Pour obtenir le profil de la figure 6, on choisit les numéros de PCC correspondant aux traces 32 et 92. En 123, l'ensemble T où sont mémorisées les valeurs calculées des # T est ramené à zéro. Des instructions figurent en 125, 127 et 129 pour vérifier les groupes de PCC sommés lorsque les bouoles de programme passent par ces groupes.En faits ces instructions permettent de sommer lee traces de diffraction des groupes 1 à 62 de PCC pour former les ensembles sommés basés sur les calculs effectués à partir du premier groupe de POC de ref4- renoue (n 32) et permettent ensuite de sommer les traces de diffraction des groupes 63 à 160 de PCC pour former des ensembles sommés basés sur les calculs, à partir du groupe suivant de PCC de référence (n 92). En 131, le TSUK est identifié à DELT et NC est identifié à 1.Les instructions DO LOOP, DO 133 dans le reotangle 133 servent à calculer l'ensemble d'intervalles AT de sortie employé pour la sommation des traces de diffraotions. Les calculs sont effectués en acoord aveo les équations 2 à 7 oi- dessus. L'intervalle #T pour le TSUM considéré est calculé par la relation T(NC)= TSK -TKF. NC est un indice d'échantillonnage qui détermine l'emploi cément de l'intervalle #T calculé dans l'ensemble des T. TSUM et NC sont ensuite augmentés et les oalouls sont effectués sur un autre échantillon. Le sous-programme GENDLT est appelé par le programme principal pour chaque trace. Ce sous-programme calcule ensuite les intervalles #T pour cette trace, puis on revient au programme principal. Dans le mode d'exécution oi-dessuss la sommation a été mise en oeuvre sur seulement deux groupes de PCC n0 32 et 92 qui sont choisis par un interprétateur. Comme indiqué ci-dessus, tous les groupes de PCC peuvent être sommés par diffraction en se basant sur des calculs à partir d'une trace de référence déplacée par échelons le long du profil. Les programmes des organigrammes 2 à 4 peuvent être modifiés pour mettre en oeuvre cette variante par introduction d'une boucle extérieure enoadrant l'instruction DO LOOP, DO 996 et en remplaçant par l'organigramme n 5 les instructions DO LOOP, DO 103 ainsi que les instructions 105, 107 et 109 de l'organigramme n 3. Organigramme n 5 - (voir page 28) 201) exécuter 137, J = 1,NCD. 202) exécuter 139, N - 1NCD 203) @@@@@@@@ la trace sommée de diffraction de l'unité N. 204) @@@@@@@ cette trace pour former un ensemble STKB. 205) instruction continue 139. 206) enregistrer la somme obtenue sur l'unité ;CD + 1 207) instruction continue 137. La boucle extérieure peut être fermée en intercalant les instructions DO 1000, L - 1, NSETS entre les instructions 57 et 61 de l'or ganigramme 2 et une instruction continue, 1000 entra les instruc- tions 85 et 89 de l'organigramme n 3. NSETS indique le nombre de soupes de PCC. I1 va de soi que la sommation des traces de réflexions ne sera pas mise en oeuvre dans les instructions de transfert. Dans le sous-programme GENDLT (orgarigramme n 4), les instructions 125, 127 et 1;9 servaient annulées et les numéros des PGC lus successivement en 121 comprendraient tous ceux formant le profil à traiter.Chaque trace sommée de diffraction serait enrg-.istrée en 91 sur un disque N différent. Avec cette variante, le premier groupe de PCC de référence serait choisi et les calculs seraient effectué jusqu'au premier groupe de PCC à sommer. Il serait sommé par diffraotion et enregistré. On effectuerait ensuite les calculs jusqu'au groupe suivant le PCC. I1 faudrait sommer par diffraction, enregistrer et recommencer les opérations jusqu'à ce que tous les groupes de PCC soient sommés par diffraction en se basant sur des calculs effectués à partir de la premiere référence. Ceci formerait un premier profil. Le groupe de PCC de référence suivant serait choisi et les opérations recommencées pour former un second profil.De cette manière, tous les groupes de PCC seraient sommés par diffraction en se basant sur tous les groupes de PCC le référence choisis. Dans ltorganigramme n 5, NCD est égal au nombre de traces de chaque profil et au nombre de profils produits. -La boucle extérieure appelle ensuite un numéro de trace tandis que la boucle intérieure ajoute les nombres de traces correspondants de chaque profil. Par conséquentS toutes les traces n 1 de tous les profils sont sommées et enregistrées; toutes les traces n0 2 de tous les profils sont sommées et enregistrées jusqu'à ce qu'un profil résultant sommé par diffraction soit obtenu. On décrit ci-apres un autre procédé ou mode d'exécution dans lequel les événements de diffraction sont sommés ou éliminés en retournant à leur source. En fait, oe procédé "comprime" les diffractions en retournant à leur sourie de telle maniere que, sur le profil traité, on obtient des événements de grande amplitude aux sources ponctuelles, ce qui accentue les sctr%-er d'énergie diffractée. Ce procédé peut traiter un profil sommé par réf3:xion ou par diffraction, par exemple ceux des figures 5 et 6. Lors de la mise en oeuvre du procédé selon l'invention, les traces sommées à l'intérieur d'une fenêtre de largeur prédéterminée sont sommées pour former une trace de référence qui est la trace centrale de cette snêtre. Cette fenêtre est déplacée le long d'une trace et les traces sommées à l'intérieur de cette feutre sont sommées tour à tour pour former la trace centrale. Le déplacement et la sommation recommencent jusqu'à obtention d'un nombre prédéterminé de traces de sortie. Dans un mode d'exé- cution, la fenêtre comportait 31 traces préalablement sommées; cependant, il va de soi que la largeur de la fenêtre peut entre augmentée ou diminuée suivant le pouvoir résolvant désiré. La figure 12 représente une première variante de ce mode d'est cution destinée à comprimer les évenements de diffraction. Sur la figure 12, la trace centrale ou de référence de la fenêtre dans laquelle les traces du profil à alignemont vertical doivent être sommées est représentee en 153. Une trace sommée à traiter est représentée en 155. Le long des traces de référence, le temps est augmenté d'une quantité finie et à chaque instant ainsi échelonné, un échantillon de trace faisant partie de la trace 155 representant un événement de diffraction enregistré stlr cette trace, est mis en place et sommé pour former un ensemble sommé correspon- dant à la trace de référence au temps échelonné correspondant. On suppose que ce temps est augmenté à 153a sur la trace 153. A partir de oe point, ltévenement de diffraction correspondant sur la trace 155 est en place pour la sommation afin de former l'ensemble sommé à l'instant 153. On suppose que cet événement de diffraction est celui représente en l=5a sur la trace 155. Les événements de diffraction sur le profil à alignement vertical, qu'il s'agisse d'un profil sommé par réflexion ou par diffraction, se trouvent sur des hyperboles. Pour sommer ces évenements de diffraction se rapportant à un point d'un profil à alignement vertical, on emploie une équation de sommation hyperbolique.Sur la figure 12, l'événement de diffraotion 155a ee trouve sur l'hyperbole 157 qui part de l'instant 153a sur la trace de référence. L'événement de diffraction 155a peut être placé sur la trace 155 en caloulant le temps de propagation TP sol-miroir-sol à partir de 153~ jusqu'à un point en surface de la trace 155 représentée en 155 S. Ce point en surface correspond au point de surface de référence du groupe de PCC à partir duquel on a établi la trace 155.En faisant pivoter la courbe de temps de propagation TP sol-miroir-sol le long de l'arc 159, l'échantillon de trace recherchée 155a sur la trace 155 et le long de l'hyperbole 157 est mis en place. Quand Cet échantillon de trace est mis en place, il est sommé pour former l'ensemble sommé correspondant b la trace de référence 153 pour le temps de propagation vertical sol-miroir-sol oorreepondant à l'instant 153a. Par conséquent, l'évènement de diffraction 155a est sommé en direction de la source ponotuelle le long de la trace de référence 153. Lorsqu'on met en oeuvre le procédé de la figure 12, la vitesse pour chaque temps échelonné le long de la trace de référence 153 est employée dans les calculs pour localiser ltéohantillon dede trace approprié sur chaque trace sommée. Les calculs destinés à localiser les événements de diffraction le long de la trace sommée pour les sommer en direction de la source en vue de former la trace de référence sont effectués tout d'abord en teppprr, puis en distance ou profondeur, et iinalement en temps. Â oe sujet, le temps est augmenté le long de la trace de référence 153 en augmentant les nombres NO sur ltéohantillon de trace et en multipliant chaque numéro d'échantillon de trace par DELT, l'intervalle temporel entre des échantillons de trace. Par exemple, oeoi détermine le temps de propagation vertical sol-miroir-sol à l'instant 153a. La profondeur à cet instant est déterminée en multipliant le temps calculé par la vitesse attribuée à l'échantillon de trace considéré le long de la trace de référenoe 153. Cette vitesse est désignée par V.NC.a. Etant donné que le temps est le tempe de propagation sol-miroil-sols on divise par deux oe produit pour obtenir la profondeur en ce point. On a par conséquent: DV n T. V.NO.a (8) 2 La figure 13 représente la prcfondeur calculée suivant une ver tical DV à l'instant 153a' qui oorrespond à l'instant 153a. Les calculs sont maintenant effectués en distance ou en profondeur.La distance de la trace de référence à la trace sommée est désignée par XH et est oonnue ou peut être calculée par des opérations sur le terrain. C'est la distance entre les points de surfaoe de référence des traces de PCC à partir desquelles les traces 153 et 155 ont été établies. La distance DD en diagonale entre les points 153a' en profondeur et 155S est calculée par l'équa- tion oi-aprèsss Cette distance sert ensuite à déterminer le temps de propagation TF solmiroir-sol entre le point en profondeur 153a' et le point 155 s a la surface.Ce temps IF de propagation se calcule par 1'équation oi-apres Après avoir déterminé par le calcul l'emplacement de ltéohantil- lon de trace au point 155a correspondant au point source 153a et après avoir sommé cet échantillon de trace pour former l'ensemble sommé à un instant correspondant à 153a, on augmente le temps le long de la trace de référence 153 e-t oaloule et effeotue une sommation sur un autre éohantillon de traoe. On effeotue les opérations jusqu'à ce que tous les éohantillons de trace soient déterminés par le calcul et sommés.On choisit une autre trace et la somme pour former l'ensemble sommé correspondant å la trame de référence 153. On effeotue ces opérations jusqu'à oe que toute. les traces à l'intérieur de la fenêtre aient été sommées pour former l'ensemble sommé correspondant à la traoe oontrale 153 de la fenOtro, On déplace ensuite cette entre le long d'une trace et exécute des calculs et une sommation semblables. On opère ainsi jusqu' l'obtention du nombre prédéterminé de traces de sortie. On se reportera oi-apres aux figures 14 et 15 pour étudier l'autre variante en vue de traiter un profil à alignement vertical afin de sommer les événements de diffraction en retournant i leur source. Dans cette variante, toutes les traces à l'intérieur de la fenêtre prédéterminée sont également sommées pour former une trace oentrale, ou de référence, 153; puis on déplaoe la fenêtre le long d'une trace et recommence l'opéra- tion jusqu'à obtention du nombre prédéterminé de traces à la sortie. Cependant, pour cette variante, le temps augmente le long de la trace à sommer au lieu d'augmenter le long de la trace de référence. On oaloule à partir de chaque échantillon de la trace à sommer l'emplacement de la source sur la trace de référence 153 dloù provient l'énergie diffractée. La vitesse pour chaque échantillon le long de la trace sommée sert à localiser la source ponctuelle correspondante le long de la trace de re- férence. Une fois cette source ponctuelle localisée, on somme ensuite l'échantillon de trace pour former l'ensemble sommé correspondant à la trace de référence pour le temps de propagation vertical sol-miroirsol 'afin de déterminer la profondeur de cette source ponctuelle. Dans cette variante, les calculs sont également exécutés en temps, puis en distance et finalement en temps. On peut supposer, en référence à la figure 14, que la trace centrale, ou de référence, est celle représentée en 153 et que les éohantillons le long de la trace 155 doivent être sommés pour former un ensemble sommé correspondant à la trace de référence 153. On augmente par échelons le temps le long de la trace 155 en augmentant par échelons le numéro 10 de l'échantillon sur la trace et en multipliant chaque numéro d'échantillon par DELT pour mettre en place l'échantillon de trace au point 155b par exemple.Ce temps de propagation sol-miroir-sol est désigné par TFb. I1 est intéressant de connaître l'emplacement sur la trace 153 de la source ponctuelle 153b d'où provient ltévenement de diffraction correspondant à l'échantillon 155b de la trace. L'emplacement de la source ponctuelle 153 sur la trace 153 peut entre défini comme étant l'intersection de la droite de temps de propagation EFb avec la trace 153 (voir figure 14).Le temps de propagation vertical sol-miroir-sol jusqu'à la source ponctuelle l53b paut être calculé en calculant d'abord la profondeur du point échantillon 155b et en utilisant cette valeur pour déterminer la profondeur et par conséquent le temps de propagation sol-miroir 801 vertioal au point 153b. La vitesse en chaque point échantillon le long de la trace sommée est utilisée dans les calcules, comme indiqué oidessus. Cette variante permet d'effeotuer des calouls, compte tenu dee variations de vitesse le long de la trace sommée et présente des avantages pour le traitement d'informations comportant des variations de vitesse très importantes. Après avoir calculé le temps de propagation sol-miroir-sol au point 155b de la trace 155, on caloule la profondeur DV' au point souterrain correspondant 155bu (figure 15). Pour cela, on multiplie la vitesse attribuée au numéro de l'échantillon en ce point par le temps calculé puis on divise le produit par deux. Cette vitesse est désignée par V.NC.b. Par conséquent: DVI T . V0NC.b (11) 2 On voit sur la figure 15 que cette profondeur est égale à la distance oblique entre le point en surface 155S et le point souterrain 153' correspondant à l'instant 153b.Si l'on connatt DV' et XE on peut calculer la profondeur en 153b' par la relation ci-aprèst Cette profondeur Z estoenvertie en temps de propagation verticale solmiroir-sol par l'équation ci-après; TFP n .2 (13) V,NG,b Ensuite, l'échantillon de trace 155b est sommé pour former l'ensemble sommé pour le temps de propagation sol-miroir-sol calculé TFF. Comme on l'a indiqué oi-dessus, les procédés des figures 12 à 15 peuvent traiter un profil sommé par réflexion ou un profil sommé par diffraction. Dien que leprosédé olassique de sommation des réflexions puisse faire disparaftre des évenements de diffraction, il ne les élimine pas tous. Ceux qui sont présents sur un profil sommé par réflexion sont sommés ou bien comprimés en direction de leur source. Lorsqu'on traite des informations concernant des réflexions, les procédés des figures 12 à 15 présentent la caractéristique additionnelle de déplacer l'iniorma- tion pour supprimer les chevauchements ou interférences du profil en temps. La figure 16 représente la manière dont on traiterait par lee opérations des figures 12 à 15 l'évenement 40 en opérant sur le profil sommé par diffraction de la figure 6. Les franges de diffraction sont comprimées et des événements de grande amplitude apparaissent aux limites de ltévenement 40 et aooentuent les limites des sources d'énergie diffrao tée. Si l'on revient à la figure 10, on observe que lorsquton met en oeuvre les procédés des figures 12 à 15 sur un profil sommé par diffraotion, le temps IKP pour lequel les événements de diffraction intéressant un groupe de PCC sont sommés en vue d'une accentuation d'amplitude permet aussi de localiser la source ponctuelle sur la trace de référence pour laquelle l'événement de diffraction accentué est sommé ou comprimé. Les organigrammes nO 6 et 7 décrivent la mise en oeuvre des procédés des figures 12 à 15 en vue de oomprimer les événements de diffraction en revenant à leur source. Ces organigrammes traitent un profil sommé par réflexion, mais ils peuvent entre modifiés pour traiter un profil sommé par diffraction. Les organigrammes 6 et 7 sont suffisamment détaillés pour permettre à un spécialiste qualifié d'écrire un programme en langage FORTRAN pour mettre en oeuvre le procédé, par exemple avec une calculatrice numérique d'usage universel CDC 6400. Si l'on étudie d'abord l'organigramme 6, l'unité 5 citée dans ltorganigramme 3 contenant le profil sommé par réflexion est ramonée au point de départ comme indiqué en 161. Le sous-pro gramme STKD est ensuite appelé comme indiqué en 163. Ce sous-programme est représenté dans l'organigramme 7.Dans oe sous-programme, on donne les définitions oi-aprèss NCTR est un compteur employé pour établir si le nombre prédéterminé de traces de sortie, NTROUT^a été obtenu. 3 est l'ensemble dans lequel les traces doivent entre sommées tandis que NTRIN est le nombre de traces a l'intérieur de la entre. Les autres variables et noms ont été décrits ci-dessus. Organigramme n 6 - (voir page 29) 161) ramener au point de départ l'unité 5 contenant l'enregistrement du profil sommé par réflexion. 163) appeler STKD. 163a) fin Organigramme n 7 - (voir page 29) Sous-pro gramme STKD 171) NOTE n O 173) B(J) n 0 ramener l'ensemble au zéro pour y faire figurer les traces de sortie sommées. 175) DO 191 I = 1,NTRIN 181) pour déterminer la distance X de la trace, calculer le numéro d'échantillon correspondant à TF. 183) sommer dans l'ensemble B l'échantillon dont le numéro correspond au oentre d'une fenêtre. 185) NC = NO + 1 187) NC > NSA ? 189) passer à 179 191) instruction continue 191 193) enregistrer l'ensemble sommé B sur un disque ou une bande magnétique, à la sortie, ou tracer une courbe. 195) NCTR n NCTR + 1 déplacer la fenêtre d'une unité 197) NCTR # NTROUT ? 197a) instruction de retour NOTE et B sont ramenés à zéro en 171 et 173. La sommation est mise en oeuvre par l'instruction DO LOOP, DO 191 qui est représentée débutant dans le reotangle 175. On suppose dans ce procédé que l'ensemble B sommé est au centre de la fenêtre qui correspond à la trace n 16 de cette fenêtre, si l'on suppose que oelle-ai comprend 31 traces. Pour traiter la première trace, on part du centre de la fenêtres on se déplace de 15 traces et on multiplie par 15 l'intervalle entre les groupes de PCC. Si l'on traite la seconde trace, on se déplace de 14 traces à partir du centre de la fenêtre et multiplie ce nombre par l'intervalle entre les groupes de PCC. On traite de cette manière chaque trace qui se trouve à l'intérieur de la entre. En 177, on identifie NO à 1 et en 179, on résout l'équation 10 pour déterminer le temps de propagation sol-miroir-sol de la manière décrite à propos des figures 12 et 13.En 181, pour la distance de la trace correspondant à celle du centre de la fenêtre à la trace traitée, l'échantillon, dont le numéro correspond au temps de propagation TF sol-miroir-sol calculé, est placé et sommé en 163 pour former l'ensemble 3 b un instant oorre- pondant b l'instant échelonné T de la figure 12. On augmente d'une unité en 185 la valeur affichée par le compteur Na et on procède à un essai en 187 pour déterminer ni BC est supérieur à NSA, le nombre d'échantillons de traces dans une trace donnée.Si la réponse est "non" alors l'instruc- tion Upasser à 179" du rectangle 189 retransfère la oommande à l'instruc- tion 179 selon laquelle le tempe de propagation sol-miroir-sol et oalculé pour l'échantillon de trace suivant le long de la trace de reiérenoe. Si NO est supérieur à SA la oommande est retransférée & l'instruction 191 continue au début de l'instruction DO LOOP, DO 191 et des opérations sont effectuées sur la trace suivante à l'intérieur de la fenêtre pour mettre en place et sommer ces échantillons pour former l'ensemble sommé 3. Après que toutes les traces à l'intérieur de la entre ont été traitées et sommées pour former l'ensemble B, la trace sommée est enregistrée comme indiqué en 193. Ensuite, comme indiqué en 195, la entre est déplacée d'une trace et on compare en 197 NOTER à NTROUT pour établir s'il est égal ou supérieur à oe dernier. Si la réponse est"non' la oommande est retransférée A 173, ou ensemble sommé 3 est ramené au zéro et la sommation est mise en oeuvre dans la fenêtre suivante. Si NCTR est supérieur ou égal à NTROUT, le sous-programme revient au programme principal. Si le procédé des figures 14 et 15 doit être mis en oeuvre, alors on remplace dans le reotangle 179, l'équation 10 par l'équation 13 ou une équation équivalente. Lors de la mise en oeuvre du mode d'exécution des figures 12 & 15 il va de soi que des fonotions vitesse-temps appropriées doivent être insorites dans le système de calcul pour mettre en oeuvre les calculs. Dans un mode d'exécution, la même fonotion vitesse-temps sert pour tous les calculs correspondant à chaque position de la entre Bien entendu, diverses modifications peuvent titre apportées par l'homme de l'art aux dispositifs ou procédés qui viennent d'être décrits uniquement à titre d'exemples non limitatifs, sans sortir du cadre de l'invention. Organigramme n 1 34 recherche de I'ajustement optimal de la34 vitesse 35 correction d'obliquité et sommati I interpretation du profil à alignement 36 vertical pour détecter les sources ponctuelles calcul des courbes de temps de propagation 37 pour l'énergie diffractée le long des traces pour des PCC et provenant de sources d'énergie diffractée, à partir des traces de référence sommation de plusieurs échantillons de traces pour chaque groupe de PCC le long de leur courbe de temps de propagation pour former des ensembles sommés afin d'augmenter l'amplitude des événements de diffraction provenant des sources spécifiées. r Programme D I FSTK I |inforlatounes 'entré 51 Organigramme nO 2 I I NcD=0 c 53 7 55 amener au zéro I'ensemble somné par réflexions STK ainsi que I'ensemble | lire le temps, leE informations de P lire le tezs, les informations de itss sinlce - au tée cotmnande DO 996 - 1, NFOLD 61 inscrire la trace à so,mer 63 localiser la distance point de tir- 65 géophone pour I'ensemble des distimces NC=O 6769 TSUM --O 69 TSUM: TSUM +DELT I . 71 NC=~C+I r 7/ 2 + X 2/ V 2 T: TSUM 75H identifier et ajouter 1'dch ntillon de trace défini par T pour obtenir 1'ensemble sommé TSK s IF (NC. EQ. NSAJ 79, 69 F I 79TSUM:O 79 Organigraimne 3 ; 3 ss appeler GENDLT 83 JI 83 effectuer 83 J = 1,NSA TSUM = TSUM + DELT A partir de TSUM, localiser l'échantillon de trace approprié défini par T et calculé par le sous programme CENDLT et ajouter cet échantillon de trace à l'ensemble de diffraction sommé STKA en TS T 996 lncsonrtuincuelon 85 89 1 inue 89 enregistrer les traces sommées de réflexions sur le disque 5. 91 enregistrer les traces sommées de diffract sur le disque 6. NCO=NCD+I 93 95 N 95 55 IF(NCD. NE. NTROUT) 55 I DO 103 iN= I, NCD97 extraire une trace sonrmée de réflexions f- 99 I du disque 5 et mémoriser extraire des trace J. I !s extraire des traces sommées de diffrac- loi tions tions du disque 6 et mémoriser. 103 ajouter les traces sommées de réflexions aux traces sommées de diffractions I enregistrer a somme sur 105 t n; to SOIllae t enregistrer es tra WsunslggemSes de diffraction sur P 107 I unSe" enregistrer les traces so-nes e re réflexions sur 109 urganigramme n0 4 Suus-programme GENDLT V l inscrire les numéros des PCC (LFT à partir 121 desquels on exécute les calculs) I ramener à zéro 11 ensemble T où les intervalles 123 A T sont mémorisés attribuer un numéro à un ensemble fictif de PCC au-delà 125 de I'ensemble de PCC le plus à gauche du profil P mettre en place les points centraux entre les - 127 PCC correspondant aux numéros à partir desquels les calculs doivent être effectués 7 à partir du numéro suivant de PCC, au-desous 129 duquel I'interprétateur suppose la présence de source de phénomènes de diffraction 7 NC= l | DO 133 J= I, NSA calculer TKF et TSK à I'aide des 133 équations 2 à 7 T (NC = TSK - TKF TSUM = TSUM+DELT 133 NC = NC+ I instructions Wo133A de de retour Organigramme nO 5 exécuter 137, J = 1,NCD | 201 1 exécuter 139, N ~ 1,NCD 202 extraire la trace sonrmée de diffraction ,203 de I'unit8 N de l'unité N ajouter cette trace pour former un ensemble 204 STKB > Z instruction continue 139 ",205 t r enregistrer la somme obtenue sur l'unité ,206 NCD + 1 instruction continue 137 207 &verbar; Sous-programme STKD &verbar; Organigramme n" 7 AICTR=O 171 ramener I'ensemble a zéro PtkUCHs faire figurer les traces de sortle sommées f75 DO 191 I- I, HTRIN 1 DO 191 ;;;TRI11 177 179 ' I(NC)a. NC. DELT )2 179 TF--2 L pour pour déterminer la distance X de la trace, calculer le numéro d'échantillon corres r > ondant à TF sommer dans l'ensemble B l'échantil 181 lon dont le numéro correspond au centre d'une entre ::NC+I 163 Organigra-ee,n06 NC = nu I I les ramener au point de NC NSA? 187 départ l'unité 5 conte nant 1 enregistremeDt 169 protil somBF par retExBa F 189 161 |appeler STKD I 191 161 ap , , r 163 Ins truc lnstruction C10 191 l63A fin | | B sur sommé B sur un disaue.ou une bande magnétiqpe, a la s rtle. ou tracer e bande une muroe I NÇTR = TI?+I I 195 Deplacer la fenêtre d'une unité r NCTR N TROUT F r ns ruc ion 197 dP In+mir 197A REVENDICATIONS 1 - Système de traitement de l'information utilisant une calculatrice numérique, caraotérisé en ce qu'il est programmé conformément å un programme particulier permettant d'additionner plusieurs traces sismiques produites à partir d'informations sismiques obtenues par plusieurs récepteurs placés à des distances "point de tir-récepteur" connues, croissantes, de part et d'autre , de plusieurs points de surface de référenoe espaoés le long d'une ligne d'erploration sismique, oe programme partioulier étant tel que défini ci-dessous. a) ohoix d'une trace sismique espacée d'un point de référence en surface, à sommer. b) détermination d'un instant auquel on prévoit l'enregistrement d'ondes sismiques le long de ladite trace choisie, ladite énergie provenant d'une source ponctuelle hypothétique au-dessous dudit point de référence en surface. c) choir dtun échantillon de trace sur ladite trace choisie, éohantillon correspondant audit instant déterminé auquel on prévoit l'enre- gistrement d'énergie diffractée. d) sommation desdits éohantillons de trace pour former un ensemble sommé à un instant qui est fonotion du temps de propagation ssol-miroir- sol des ondes sismiques, entre ladite source ponctuelle hypothétique et le point en surface correspondant à ladite trace choisie. e) répétition des opérations a) & d) pour plusieurs sources ponotuelles hypothétiques différentes placées à différentes profondeurs au-dessous dudit point de référence en surface. f) répétition des opérations a) à e) pour plusieur. traces différen- tes choisies. 2 - Système selon la revendication 1 dans lequel lesdites traces sismiques à sommer sont des traces sommées antérieurement, produites chacune par sommation d'informations enregistrées par des récepteurs placés à des distances "point de tir-réoepteur" variables de part et d'autre d'un point de surface de référence, caractérisé en oe qu'il com- prend les opérations ci-après, détermination le long de chaque trace choisie des instants auxquels on prévoit l'enregistrement d'ondes sismi- ques par augmentation de quantités finies de différents temps oorrespondant aux instants au-dessous dudit point de référence en surface et emploi de la vitesse de propagation pour la oouohe considérée à ohàque instant échelonné pour calculer le temps de propagation sol-miroir-sol pour chaque temps ainsi échelonné jusqu'à la surface, correspondant à ladite trace choisies utilisation des temps de propagation sol-miroir-sol calculés pour chaque trace choisie en vue de choisir des échantillons de trace respectivement le long de chacune desdites traces choisies et sommation des échantillons de trace de chaque trace choisie pour former ledit ensemble sommé, à des instants correspondant audits temps échelonnés k partir desquels lesdits temps de propagation sol-miroir-sol sont calculés. 3 - Système selon la revendioation 1 caractérisé en ce que les traces sismiques à sommer sont des traces sommées antérieurement produite chacune par sommation d'informations enregistrées par des réoepteurs placés à des distances "point de tir-récepteur" croissantes, de part et d'autre d'un point de surface de référence, chaque trace sommée antérieurement comportant un point de surface de référence, oarac- térisé en oe qu'il comprend les opérations oi-aprèss détermination sur chaque trace ohoisie des instants auxquels on prévoit l'enregistrement d'ondes sismiques par augmentation de quantités finies de temps diffe- rents le long de chaque trace choisie, emploi de la vitesse de propagation dans une couche souterraine pour chaque temps ainsi échelonné le long de ladite trace ohoisie pour calculer les temps de propagation solmiroir-sol correspondant au-dessous dudit point de référence en surface dtoù l'on suppose que les ondes diffractées proviennent, emploi desdits temps échelonnée le long de chaque trace choisie pour choisir des éohan- tillons de trace respeotivement le long de chacune desdites traces ohoisies et sommation des échantillons de trace sur chaque trace choisie pour former ledit ensemble sommé à des instants correspondant audit temps de propagation sol-miroir-sol calculé au-dessous dudit point de référence en surface, d'où lton suppose que l'énergie diffraotée provient. 4 - Système selon la revendioation 1, caractérisé en oe que lesdites traces sismiques à sommer oomprennent un groupe de traces sismiques produites par des récepteurs placés à des distances croissantes du point de tir, de part et d'autre d'un point de surface de référence, le long d'une ligne d'exploration, chaque distance "point de tir-réoepteur" variable de chaque groupe définissant un groupe "point de tir-réoepteur" caractérisé en'ce qu'il comprend:: l'augmentation de quantités finies des temps différents correspondant aux sources ponctuelles hypothétiques audessous dudit point de référence en surfaces calcul du temps de propagation sol-miroir-sol des ondes sismiques à partir de chaque temps échelons né en direction du point de surface de référence dudit groupe, calcul du temps de propagation des ondes sismiques entre chaque paire "point de tir-réoepteur" et chaque temps échelonné, choix d'échantillons de trace à partir de chaque trace dudit groupe de traces correspondant audit temps de propagation "point de tir-réoepteur" calculé et sommation des échantillons de trace choisis dans chaque trace dudit groupe de traces donné pour former ledit ensemble sommé, à des instants correspondant au temps de propagation sol-miroir-sol associé calculs 5 - Système de mise en oeuvre de la sommation par diffraotion de plusieurs traces sismiques produites, à partir d'informations sismiques obtenues par plusieurs récepteurs placés à des distances oroissantes "point - tir-réoepteurt' oonnues de part et d'autre de plusieurs points de surface de référence espacés le long d'une ligne d'exploration sismique, ledit système utilisant un équipement de oaloul automatique, caractérisé en ce qu'il comprend les opérations oi-aprèss a) choix d'une trace à sommer pour former un ensemble sommé, dont la position correspond à oelle d'une trace de référence espacée horizontalement de ladite trace ohoisie, b) échelonnement jusqu'à un point sur une desdites traces faisant l'objet de ltopération 4, o) calcul d'une distance qui définit un aro dont le oentre est sur la surface correspondant à celle desdites traces ohoisies et qui ooupe ladite trace auxdits points échelonnés et l'autre trace en un point intéressant, d) choix de l'échantillon de ladite trace choisie à un instant correspondant au point d'intersection dudit aro aveo ladite trace choisie, e) sommation dudit échantillon de trace pour former ledit ensemble sommé pour le temps de propagation, dans la direction verticale, solmiroir-sol correspondant au point d'interseotion dudit aro aveo ladite trace de référence, f) répétition plusieurs fois des opérations b) à e), g) répétition des opérations a) à f) pour chaque trace d'un groupe de traces différentes à sommer pour former ledit ensemble sommé lesdits échantillons de trace ayant pour origine plusieurs traces comportant différents points de surface de référence. 6 - Système selon la revendication 5 caractérisé en ce que l'opération b) est mise en oeuvre par éohelonnement jusqu'à un point le long de ladite trace de référence et l'opération o) est mise en oeuvre en calculant la distance qui définit un aro ayant son oentre sur la surface correspondant à ladite trace choisie et qui coupe ladite trace de réfé- rence audit point obtenu par échelonnement et ladite trace ohoisie audit point intéressant. 7 - Système selon la revendication 5, caractérisé en ce que l'opération b) est mise en oeuvre par échelonnement jusqu'à un point le long de ladite trace choisie et l'opération o) est mise en oeuvre en calculant la distance qui définit un aro dont le oentre est sur la sur faoe correspondant à ladite trace ohoisie et qui coupe ladite trace ohoisie audit point obtenu par échelonnement et ladite trace de référence audit point intéressant. 8 - Système de mise en oeuvre de la sommation par diffraction de plusieur. traces sismiques produites à partir d'informations sismiques obtenues å partir de plusieurs récepteurs placés à des distances oroissantes "point de tir-réoopteur" oonnuess de part et d'autre de plusieurs points de surface de référence espacés le long d'une droite d'exploration sismique, ledit système comportant un équipement de calcul automatique, caractérisé en oe qu'il comprend les opérations oi-aprèss a) choix d'une trace à sommer pour former un ensemble sommé dont l'emplacement oorresponià une trace de référence espaoée horizontalement de ladite trace choisie, b) échelonnement du temps jusqu'à un point le long de ladite trace de référence o) calcul du temps de propagation sol-miroir-sol d'ondes sismiques entre ledit point le long de ladite trace de référence et la surface correspondant à ladite trace ohoisie, d) choix,eur ladite trace ohoisie,de l'échantillon de trace qui correspond audit temps de propagation sol-miroir-sol calculé, e) sommation dudit échantillon de trace choisi pour former ledit ensemble sommé pour le temps de propagation dans la direction vertioale sol-miroir-sol en vue d'obtenir ledit point déterminé par échelonnement le long de ladite trace de référence, f) répétition, plusieurs fois, des opérations b) à e), g) répétition des opérations a) à f) pour chaque trace d'un groupe de traces différentes en vue de les sommer pour former ledit ensemble sommé, les échantillons de trace provenant de plusieurs traces aveo des points de surface de référence différents. 9 - Système de mise en oeuvre de la sommation par diffraction de plusieurs traces sismiques produites à partir d'informations sismiques obtenues par plusieurs récepteurs placés à des distances croissantes "point de tir-réoepteur" connues, de part et d'autre de plusieurs points de surface de référence le long d'une droite d'exploration sismique, ledit système comportant un équipement de calcul automatique, caractérisé en ce qu'il comprend les opérations ci-après:: a) choix d'une trace à sommer pour former un ensemble sommé dont l'emplacement correspontà une trace de référence espacée horizontalement de ladite trace choisie, b) échelonnement du temps jusqu'à un point sur ladite trace choisie, o) calcul de la distance dans la direction verticale audit point déterminé par échelonnement de temps et qui définit un aro dont le oentre est sur la surface correspondant à ladite trace choisie, et qui ooupe ladite traie de référence en un point intéressant, d) choix de l'échantillon de trace sur ladite trace choisie audit point déterminé par échelonnement du temps, e) sommation desdits échantillons de trace pour former ledit ensemble sommé pour le temps de propagation sol-miroir-sol, dans la direction verticale, correspondant à la profondeur du point d'intersection dudit aro avec ladite trace de référence i) répétition, plusieurs fois des opérations b) à e), g) répétition des opérations a) b i) pour chaque trace d'un groupe de traces différentes à sommer pour former ledit ensemble à sommera les échantillons de trace provenant de plusieurs traces comportant des points de surface de référence différents. 10 - Système d'accentuation d'éwdenements de diffraction sur un profil sismique formé à partir de plusieurs groupes de traces sismiques, les traces de chaque groupe étant produites par des récepteurs placés b des distances "point de tir-récepteur" connues et croissantess de part et d'autre d'un point de surface de référence, le long d'une droite d'exploration sismique, chaque intervalle croissant n point de tir récepteur d'un groupe donné définissant un couple "point de tir-réoepteur le point de surface de référence de ohaoun desdits groupes oi-dessus se trouvant à une distance croissante du point de surface de référence d'un groupe prédéterminé, ledit système comportant un équipement de calcul automatique, caractérisé en ce qutil comprend les opérations ci-apress a) choix d'un point de référence en surface basé sur lthypothese qu'une source d'énergie diffractée se trouve au-dessous dudit point de référence en surface, b) établissement de toutes les diverses fonctions de temps de propagation "point de tir-récepteur" , représentant le temps de propagation des ondes sismiques respectivement entre chaque paire "point de tir récepteur" d'un groupe donné, grace à un échelonnement des temps pour des points au-dessous dudit point de référence en surface, chaque fonction de temps de propagation "point de tir-réoepteur" définissant l'instant d'apparition d'un échantillon de trace sur une trace enregistrée par un récepteur d'une paire "point de tir-récepteur" donnée, o) établissement de fonctions de temps de propagation sol-miroir-sol représentant le temps de propagation sol-miroir-sol entre le point en surface commun dudit groupe et lesdits points correspondant auxdits temps échelonnés, d) mise en place sur une tracte dudit groupe donné de l'échantillon de trace défini par sa fonction de temps de propagation "point de tir-récep- teur" oorrespondante, e) audition dudit échantillon de trace mis en place par l'opération d) pour former un ensemble sommé pour un temps déterminé par sa fonction de temps de propagation sol-miroir-sol correspondante, f) répétition des opérations d) et e) pour chaque échantillon de chaque trace dudit groupe de traces, g) répétition des opérations b) à f) pour chaque groupe desdites tra oes pour produire un profil comprenant un grand nombre de traces sommées dans lesquelles les événements de diffraction sont accentués. 11 - Système selon la revendication 10, caraotérisé en ce qu'il comprend: le choix d'un point de référence en surface différent et la répétition des opérations b) à g). 12 - Système de mise en oeuvre de la sommation par diffraction de plusieurs traces sismiques sommées antérieurement produites à partir d'informations sismiques obtenues par plusieurs réoepteurs plaoés à des distances "point de tir-réoepteur" connues et croissantes, de part et d'autre d'un groupe de points de surface de référence espacés le long d'une droite d'exploration sismique,ledit système comportant un équipement de calcul automatique, caractérisé en ce qu'il comprend les opérations oisprèss S a) choix d'une fenêtre oomprenant un nombre prédéterminé desdites traces sommées antérieurement, b) choix d'une trace sismique à l'intérieur de ladite fenêtres o) détermination, le long de ladite trace sismique choisies des temps auxquels on prévoit ltenregi trement d'énergie diffractée et provenant de sources ponctuelles hqzothétinues à diverses profondeurs, le long d'une trace sismique prédéterminée, à l'intérieur de ladite fenêtres d) ohoix d'échantillons faisant partie de ladite trace choisie et correspondant auxdits temps déterminés en o) pour lesquels on prévoit 1 'enregistrement d'énergie diffractée, e) sommation desdits échantillons de trace pour former un ensemble sommé pour des temps correspondant respectivement aux temps déterminés en o) pour lesdites sources ponctuelles hypothétiques, f) répétition des opérations b) à e) pour chaque trace sismique à l'intérieur de ladite fen8tre afin de sommer, en vue de former ledit ensemble sommé, des échantillons de trace provenant de plusieurs traces sismiques comportant des points de surface de référence différents, g) déplacement de ladite fenêtre jusqu'à ce qu'elle occupe différentes :positions le long desdites traces précédemment sommées, h) mise en oeuvre des opérations b) à f) pour chaque position pour former un profil, comportant un grand nombre de traces sommées dans lequel les évenements de diffraction provenant desdites traces sommées antérieurement sont comprimés.