Brevet FR2283449A1 — publié le 19760326

Dans l'exploration sismique, on utilise des ondes acoustiques pour sonder la croûte du globe terréstre afin de déterminer la nature et iremplacement -de formations souterraines. La croûte terrestre peut être considérée comme un milieu de transmission ou un filtre dont les caractéristiques sont déterminées en y faisant passer des ondes acoustiques. Dans le procédé sismique par réflexion, des ondes acoustiques sont émises en un point d'émission à la surface du sol ou prés de celle-ci et les ondes acoustiques réfléchies par les limites réfléchissantes des surfaces souterraines sont reçues en un point déterminé. Les ondes reçues sont détectées et enregistrées sous une forme permettant leur exploitation. Des spécialistes qualifiés peuvent déduire d'une telle exploitation la forme et la profondeur des limites des surfaces réfléchissantes souterraines et la probabilité de trouver des accumulations de matières minérales telles que le pétrole et le gaz naturel. Une des premières réalisations dans le domaine de l'ex- ploration sismique émettait de l'énergie sismique en un point, sous forme d'impulsions provenant en général de l'explosion, ou tir, d'une charge de dynamite. Une grande quantité d'énergie était ainsi injectée dans le sol pendant un court instant. Dans un autre mode d'exploration sismique, décrit dans le brevet des Etats-Unis d'Amérique NO 2 688 124, de l'énergie sismique peut être injectée dans le sol pendant une période de temps prolonge et les ondes reçues en un point déterminé sont sommées. Par conséquent, on injecte, par des procédés du type de celui dénommé "Vibroseis" de énergie sismique dans le sol sous forme d'une onde continue de fréquence variable. tes ondes résultantes reçues en un point de détection sont mises en autocorrélation avec un signal analogue provenant de l'onde injectée pour comparer les ondes émises et les ondes reçues afin de réaliser un enregistrement représentatif des caractéristiques acoustiques du sol. Lorsqu'on utilise le procédé dit "Vibroseis" d'exploration sismique, le minutage, l'amplitude et es relations de phase des signaux émis et enregistrés ont de 1'importance. Pour toutes les fréquences d'émission d'un signal d'entrée continuement variable, la reproduction dans l'appareil d'enregistrement doit être fidèle, sinon il faut procéder à une correction ultérieure pour rétablir la phase et l'amplitude correctes. Il faut pour cela que l'équipement d'enregistrement soit exempt de déphasages impossibles à compenser dans la totalité de la bande de fréquences du signal émis. Aucun filtrage en fréquences du bruit ambiant dans la gamme entière des fréquences émises n'est toléré.Par ailleurs, si une onde entretenue de fréquence variable est utilisée comme signal d'entrée, ni les ensembles émetteurs ni les ensembles récepteurs ne peuvent être accordés sur une fréquence particulière pour supprimer le bruit indésirable Dans un autre procédé d'exploration sismique décrit dans le brevet des Etats-Unis d'Amérique NO 7 182 743, une série discontinue de trains d'ondes, chacun à fréquence constante, mais dont les valeurs des fréquences forment une progression arithnetique, sont émis en direction du sol. Les ondes reçues en un point déterminé sont sommées avec des phases appropriées pour produire un seul signal de sortie avec une forme d'onde qui permet de reconnaitre des phénomènes tels que des réflexions ou réfractions sismiques. Avec une solution telle que celle du brevet NO 3 182 743 précité, le brevet des Etats-Unis d'Amérique NO 3 259 878 décrit émission d'ondes entretenues sous forme de trains d'énergie sismique, chaque train ayant une fréquence différente, et la dé termination-à partir des ondes sismiques émises et reçues- de l1amplit;ude et de la phase relative de ces ondes, ensuite la création de signaux équivalents ayant une amplitude et un déphasage appropriés, qui sont additionnés de manière à produire un signal réfléchi complexe. Le procédé du brevet des Etats-Unis d'Amérique NO 3 259 878 précité est semblable à celui décrit dans les brevets des Etats-Unis d'Amérique NO 3 221 297 et NO 3 568 142.Le brevet des Etats-Unis d'Amérique NO 3 332 512 décrit encore un autre procédé d'exploration par fréquences séparées dans lequel chaque train dténergie est produit par addition, par la méthode de Fourier de plusieurs ondes sinusoîdales de fréquences différentes pour produire un train d'ondes à enveloppes courbes. Il est connu que 11 énergie sismique émise à travers le sol subit une. atténuation dont le taux varie en fonction de la fréquence de l'onde. En général, les fréquences élevées sont plus atténuées que les fréquences basses et l'atténuation totale est fonction de la distance parcourue. Cette atténuation relative des ondes sismiques aux différentes fréquences influe beaucoup sur la forme du signal qui peut être produit par sommation d'ondes de fréquences différentes émises en un point déterminé. La présente invention concerne un procédé perfectionné d'exploration sismique grâce auquel une correction appropriée de l'atténuation et des variations de vitesse(fonctions de la fréquence, des ondes sismiques se propageant dans le sol est réalisée. Selon l'invention, une série de trains d'ondes distincts à fréquence constante sont émis'en direction du sol, chacun de ces trains ayant une fréquence médiane différente, mais faisant partie de préférence d'une série de fréquences dont les valeurs sont en progression arithmétique ou géométrique. les signaux résultants reçus en un point déterminé après réflexion et réfraction par les formations souterraines sont sommés avec une pondération, variable dans le temps, de 1'amplitude de chaque signal pour compenser les variations fonction de la fréquence de l'at ténuation dans le sol. La sommation des ondes reçues en un point déterminé conduit à la création d'un signal "idéal" qui correspond à l'irjection d'une impulsion au point d'émission. melon une forme de réalisation de l'invention, chaque train d'ondes est modulé en amplitude pour produire un train d'ondes tronqué. Cette modulation en amplitude donne le signal le plus parfait par sommation, avec le nombre minimal de trains d'ondes dans une série, lorsque les fréquences médianes et les largeurs de bande respectives de ces trains d'ondes sont convenablement choisies. L'invention permet de se rapprocher d'un signal de recherche idéal à tous les instants de ltenregistrement et de conserver une résolution optimale pour toutes les profondeurs intéressantes. Par ailleurs étant donné que tous les trains d'ondes séparés ont une fréquence constante, tous les enregistrements séparés peuvent être effectués avec un filtre à bande étroite adapté à la fréquence du train d'ondes émis de manière à atténuer considérablement les parasites ambiants dus au ol. De plus, les ensembles de géophones et de sources peuvent être réglés sur la fréquence du signal de manière à donneur une atténuation optimale du bruit engendré par la source pour chaque fréquence cmEsadérée. L'invention sera décrite plus en détail en regard des destins annexés à titre d'exemple nullement limitatif et sur lesq2tB la figure 5 est un tracé dtun signal émis pour chacune des fréquences choisies ; la figure 2 est un tracé d'un signal obtenu par sommation pour un groupe particulier de fréquences émises les figures A et 3B sont des tracés de signaux de. som- mation pour des groupes additionnels de fréquences d'entrée destinés à illustrer la variation du signal de sommation en fonction de la largeur de bande à l'émission ;; la figure A est une courbe en coordonnées bi-logarithmiques dans laquelle la largeur de bande est portée en abscisses et 11 amplitude relative est portée en ordonnées et représente deux lobes lat8raux se rapportant au signal de sommation selon l'invention et à celui obtenu par le procédé dit "Vibroseis" les figures SA à 5D sont des tracés de signaux représentant. par un exemple fictif la mise en oellvre du procédé selon l'invention en vue de corriger l'atténuation, fonction de la fré qAenee, par le sol. Par ailleurs, le tableau I est un organigramme représentant le traitement selon un mode d'exécution de l'invention ; et le tableau II est un organigramme représentant le traitement selon un autre mode d'exécution de l'invention. Le tableau I représente une série, ou groupe, de trains d'ondes à fréquence constante, de fréquences différentes F1, F2... Fn, qui sont émis en direction du sol comme l'indiquent les opérations 10 à 12. Chaque train d'ondes du groupe peut être considéré comme un train d'énergie distinct et séparé du type représenté en général par le signal indiqué sur la figure 1. Lés fréquences médianes de ces trains d'ondes font, de préférence, partie d'une progression arithmétique ou géométrique. Pour une progression arithmétique des fréquences, la différence entre- chaque fréquence de la série et la précédente est constante ; par-exemple, les fréquences F1, 2....?n peuvent être respectivement 20, 25, 30, 35 et 40 Hz.Pour une progression géométrique des fréquences, le rapport entre chaque fréquence de la série et la fréquence précédente est constant ; par exemple les fréquences S 2....Fn peuvent être égales à 15,0, 18,9, 23,8, 30,0, 37,8, 47,6 et 60,0 Hz. Ces trains d'ondes sont de préférence modulés en amplitude comme on l'explique plus en détail ci-après, mais ce n'est pas indispensable. Lorsqu'fun équipement approprié est disponible, tous les trains d'ondes de fréquences médianes différentes peuvent être émis en même temps et les signaux de sommation formés directement en temps réel lorsque les signaux réfléchis sont reçus en un point déterminé. Ceci nécessite un certain nombre de sources pour émettre simultanément les trains d'ondes séparés. Toutefois, dans un mode préféré de mise en oeuvre, tous les trains d'ondes sont engendrés et mis successivement en direction du sol par la même source rayonnante et on réalise des enregistrements, comme l'indiquent les opérations 13 à 15, des ondes reçues en un ou plusieurs points de détection pour chaque fréquence. Il peut être avantageux de répétsr plusieurs fois la production et l'émission des trains d'ondes pour certaines de ces fréquences, si bien que ltenregistrement correspondant à chaque fréquence présente un rapport signal/bruit ambiant satisfaisant.On fait ensuite passer les signaux enregistrés à travers des filtres à bande étroite représentés en 16, 17 et 18, chaque filtre étant destiné à laisser passer principalement les ondes dont la fréquence correspond au canal d'émission associé et à éliminer toutes les autres fréquences y compris le bruit ambiant. Tous les signaux filtrés produits par les opérations 16 à 18 sont ensuite soumis à une pondération d'amplitude pour compenser les pertes par atténuation par transmission dans le -sol, qui varient en fonction de la fréquence. Cette pondération d'amplitude est de-préférence variable avec le temps selon une loi connue ou à déterminer, et destinée à compenser la variation des pertes par atténuation. On fait ensuite passer tous les signaux, après filtrage et pondération d'amplitude, par une jonction sommatrice 22 qui produit un signal de sommation de sortie lequel est enregistré en 23. Il va de soi que chaque signal peut être enregistré séparément etque les opérations de pondération et de sommation peuvent être effectuées postérieurement à 1'enregistrement sur le terrain. Le schéma fonctionnel du tableau I représente la série d'opérations d'émission d'enregistrement et de traitement des signaux pour une seule station à géophone. Dans la pratique normale de l'exploration sismique, on utilise plusieurs stations à géophones et plusieurs canaux d'enregistrement sont reliés à ces stations à géophones. L'organigramme du tableau II ci-agrès est destiné à représenter les opérations d'enregistrement et de traitement pour un tel canal d'enregistrement qui est utilisez pour enregistrer successiemant des signaux pour un nombre de trains d'ondes séparés égal au nombre de trains ondes émis. Par conséquent, selon l'invention, on enregistre un groupe de trains d'ondes émis dont chacun contient une petite quantité d'énergie sismique et dont le total une fois sommé à la sortie produit un signal de grande amplitude, comparé au bruît. Le signal de sortie obtenu par sommation, de la manière décrite ci-dessus, comprend des représentations des trains d'ondes émis pour chaque fréquence particulière, modifiés par les trajets de propagation dans le sol, plus l'effet du bruit. Le signal effectivement transmis pour le résultat combiné peut être déterminé par sommation des signaux correspondants à chacun des trains d'ondes. Par exemple, la figure 2 représente le signal obtenu par sommation des signaux du type représenté sur la figure 1 avec des fréquences de 20, 25, 70, 35.et 40 Hz. On voit que, à la suite d'une interférence destructive, la plupart des pulsations latérales de la forme d'onde résultante représentée sur la figure 2 ont subi une réduction d'amplitude et on observe une seule pointe principale 26. Le signal représenté sur la figure 2 est donc très proche d'un signal idéal à introduire dans le trajet de propagation dans le sol. Comme on l'a indiqué ci-dessus, chaque train d'ondes émis selon l'invention est de préférence modulé en amplitude. Cette modulation d'amplitude supprime la répétition des impulsions de chaque coté de la pointe principale (par exemple 26 sur la figure 2) sans qutun très grand nombre d'ondes séparées soit nécessaire. Pour les fréquences dont 1es valeurs sont en progression arithmétique, l'impulsion principale 26 est répétée à des - intervalles égaux à.la période correspondant à l'inverse de l'augmentation de fréquence entre deux fréquences émises successives. Par exemple, pour la série de fréquences 20, 25, 30, 35 et 40 Hz, l'impulsion principale sera émise tous les 1/5 s.Par conséquent, une période de récurrence des impulsions de 5 s, avec une gamme de fréquences de 20 à 40 Hz exigera 5 x 20 ou 100 fréquences distinctes à intervalles réguliers de 0,2 Hz entre 20 et 40 Hz à moins que tous les trains d'ondes ne soient modulés en amplitude. Pour une progression arithmétique, la fréquence de modu lation est de préférence telle que l'amplitude instantanée de chaque train d'ondes varie d'un maximum à zéro en un temps égal à l'inverse de la différence de fréquences entre les trains d'ondes contigus de cette progression. En d'autres termes, la fréquence de modulation doit être au moins égale à la moitié de la différence de fréquences entre deux trains d'ondes consécutifs de ladite progression. La fréquence de modulation la plus élevée admissible en pratique dépend des caractéristiques du vibrateur ou émetteur d'ondes acoustiques et du coefficient d'amortissement électrique des filtres d'enregistrement utilisés. La forme de l'enveloppe de l'onde modulatrice est de préférence sinusoïdale. Un grand nombre de types d'enveloppes d'ondes modulatrices conviennent. avec des conditions appropriées imposées à leur fréquence de modulation. La forme actuellement préférée d'onde de modulation a une enveloppe qui a la forme d'une demi-période d'une fonction cosinus carré qui peut être définie mathématiquement comme suit A(t) = cos2 2s M(t-t0) (1) pour t # t0 tM/2 > to ~ TMx2 A(t) = o (2) pour t H7 to + - to TM/2 formules dans lesquelles : N = l'inverse de la durée de l'onde to = l'instant de référence correspondant à l'axe de symétrie de l'onde. On peut montrer que la résolution en amplitude de l'impulsion effective de sommation augmente en même temps que la largeur de bande de la série de trains d'ondes émis. Une largeur de bande de l'ordre de 2 à 1 (un octave) est considérée comme la largeur minimale qui assure une résolution utilisable. ta figure 3A représente l'impulsion résultant d'une sommation en phase de trains d'ondes modulés en amplitude de fréquences 15, 20, 25, 30, 35, 40 et 45 Hz (largeur de bande =3 La figure 3B représente l'impulsion résultant d'une sommation en phase de trains d'ondes de fréquences 10, 15, 20, 25, 30, 35, 40, 45 et 50 Hz (largeur de bande = 5/1). Dans tous ces exemples, chaque train d'ondes est modulé en amplitude ave une seule alternance d'une onde sinusoïdale de fréquence 2,5 Hz.Bien que-, selon les figures 3A et 33, la largeur de bande ait été augmentée en augmentant le nombre de termes dans les séries, un nombre moindre de trains d'ondes produirait des résultats semblables si la modulation d'amplitude devait être modifiée pour éviter la répétition des signaux pulsés. La figure 4 représente à titre de comparaison l'amplitude des lobes secondaires du signal de sommation par rapport à l'amplitude de crête des odelettes pour la présente invention et pour le système connu sous le nom de "Vibroseis" pour différentes largeurs de bande. La courbe 30 représente l'amplitude relative du lobe secondaire le plus important de même polarité que lton- delette de phase zéro, ou ondelette de crête rapportée à l'amplitude de cette dernière ondelette pour le procédé selon la présente invention ; la même grandeur correspondant au système "Vibroseis1, est représentée par la courbe 31.Le lobe secondaire le plus important de même polarité que l'ondelette maximale est représenté en 26A sur la figure 2. la courbe 33 représente l'amplitude relative du premier lobe secondaire, de polarité opposée à l'ondelette de crête, rapportée à l'amplitude de ladite onuelatte pour le procédé selon l'insention ; et la courbe 34 représenta le même type de grandeur poule système '1Vibroseis". Le premier lobe secondaire avec une polarité opposée à celle de ltondelette de crête est représenté en 263 sur la figure 2. On voit ainsi que comparée au système "Vibroseis", la présente invention procure une impulsion de sommation plus nette et plus pointue avec des lobes secondaires plus fortement amortis qutavec le système "Vibroseis", pour la même largeur de bande. Il va de soi qu'une caractéristique propre de l'appareil émetteur d'ondes acoustiques, du type radiateur ou vibrateur utilisé pour produire des trains d'ondes selon la présente invention est la production de signaux d'amplitude plus faible aux fréquences les plus élevées, qui sont soumises à l'atténuation maximale au cours de leur propagation dans le sol. Par conséquent, il peut être avantageux de réaliser et combiner un nombre plus important d'enregistrements successifs à meure que la fréquence augmente pour obtenir un rapport signal/bruit ambiant suffisant pour la fréquence maximale et la zone intéressante la plus profonde.Ce renforcement sélectif des fréquences les plus élevées n'est pas possible avec le système 'Vibroseis" décrit dans le brevet des Etats-Unis d'Amérique NO 2 688 124 précité. Si l'on revient au tableau I, un traitement selon le procédé décrit et illustré peut être facilement effectué avec un équipement connu et classique par les spécialistes. Il est donc inutile de donner une description détaillée de types particuliers d'équipement ou de s'y référer. Les opérations de traitement peuvent être réalisées en temps réel ou en temps non réel avec un retard, après mémorisation dans une mémoire appropriée. Les opérations de filtrage 16 à 18 peuvent être réalisées avant les opérations d'enregistrement 13 à 15 par un équipement d'enregistrement sur le terrain, de manière analogique ou numérique. En variante, les opérations de filtrage 16 à 18 peuvent être réalisées après les opérations d'enregistrement 13 à 15 sur place ou dans un centre de traitement de ltinformation. Ires opérations de pondération de l'amplitude 19 à 20 peuvent être accomplies en temps réel à l'instant de l'enregistrement sur le terrain ou ultérieurement dans un centre de traitement de l'information, de manière analogique ou numérique. Par exemple, cette pondération. peut être réalisée, dans le cas d'une pondération constante dans le temps, simplement par des amplificateurs dans chaque canal dtenregistrement, si une amplification est nécessaire, ou par des atténuateurs si une diminution d'amplitude est nécessaire. En variante, si une pondération en amplitude variable dans le temps convient pour compenser les pertes par atténuation par propagation dans le sol, desWcircuits de régulation du gain peuvent être utilisés. Il est également évident que diverses opérations de pondération et de sommation peuvent être exécutées plusieurs fois jusqu'à ce que le résultat optimal désiré soit atteint. Un appareillage convenable pour l'émission des trains d'ondes selon la présente invention peut être choisi parmi les appareillages utilisables actuellement connus ou disponibles sur le marché. Par exemple, un vibrateur utilisable pour des opérations sur la terre ferme est vendu par la firme Westinghouse Air Brake Company. Evidemment, des types de vibrateurs différents doivent être utilisés pour les opérations sous l'eau ou à la mer. L'affaiblissement des diverses fréquences au cours de la propagation dans le sol peut être déterminé approximativement à partir de nombreuses mesures effectuées sur des sols types. Dans des cas particuliers, il serait possible de préciser cette information par des mesures dans des sondages réels pour la profondeur de pénétration choisie. Toutefois, il va de soi que l'aspect des courbes d'information elles-mêmes donnera des indications valables concernant les coefficients de pondération optimaux. Divers programmes de pondération peuvent être appliqués successivement pour optimaliser empiriquement le résultat final. les figures 5A à 5D montrent avec un modèle artificiel réel l'efficaeité de la présente invention. La figure 5A représente un exemple de pointes concernant la caractéristique de transmission des ondes acoustiques du sol. Chaque maximum sur la figure 5 représente le coefficient de réflexion d'une discontinuité acoustique dont ltemplaeemer,t correspond à la durée de la propagaticn dans les deux sens dans le sol, en fonction des valeurs représentées. Un signal positif pour une -pointe indique un coefficient de réflexion positif correspondant à une absence d'inversion de phase et un signe négatif pour une pointe indique une inversion de phase pour la discontinuité acoustique considérée. Pour déterminer les effets de ltatténuation, tous les trains d'ondes à fréquence constante faisant partie d'une série avec des fréquences différentes, chaque fréquence ayant une valeur définie par une progression arithmétique, sont combinés avec le modèle de pointe représenté sur la figure 5A. Les fré quences utilisées sont 15, 20, 25, 30, 35, 40-et 45 Hz. Il est tenu compte des effets de l'atténuation par le sol par l'utilisa tion d'un coefficient a d'atténuation donné par l'équation ci après Atténuation en dB = 8000 a . (FN - FR) . T dans laquelle F N est la fréquence du canal en Hz FR est la fréquence de référence en Hz T est le temps nécessaire à l'aller et retour, en secondes On admet que la vitesse est constante et égale à 2400 m/s dans l'équation ci-dessus. La figure 5B représente le signal de sommation obtenu pour un facteur a d'atténuation égal à 0,12, valeur plutôt éle vée. Une comparaison de la figure 5B avec la figure 5A montre qu'il est difficile de choisir les discontinuités correspondant aux maxima ou pointes de coefficient de réflexion. La figure 5C représente le signal de sommation obtenu quand le coefficient a d'atténuation est égal à 0,012, une valeur très faible. La figure 5D représente le signal de sommation produit selon la présente invention quand les amplitudes des signaux ré sultants produits par composition de chaque train d'ondes avec le modèle de pointe de la figure 5A sont ajustés à des valeurs d'amplitude comparables avant sommation. Il est évident, par la comparaison de la figure 5D avec les figures 5B et 5C, que la présente invention conduit à des résultats améliorés. Sur la figure SD, les ondelettes pointues correspondent plus directement et d'une manière plus identifiable aux pointes ou maxima du mo dèle de sol selon la figure 5A. La description ci-dessus est basée sur l'hypothèse que la vitesse de propagation des ondes sismiques dans le sol est exactement la même à toutes les fréquences, et cela peut être con sidéré comme vrai dans'la plupart des cas. Toutefois, il existe des cas où il est connu que la vitesse de propagation varie avec la fréquence. Dans un tel cas, la présente invention permet aussi d'appliquer une correction temporelle variable avec le temps aux signaux des diverses fréquences avant la sommation et ceci conserve le signal idéal en assurant une coïncidence exacte des phases des signaux de toutes les fréquences. Il est nécessaire dans ce cas d'appliquer une correction légèrement différente, variant avec le temps, à chaque enregistrement pour une fréquence parti coudière avant toute sommation.Par exemple, l'enregistrement correspondant à la fréquence pour laquelle la vitesse est maxi- male peut être choisi comme référence. Ensuite, on applique un decalage temporel variant avec le temps pour tous les autres enregistrements. Le décalage nécessaire augmentera progressivement à mesure que le temps de propagation augmentera et le taux d'augmentation sera fonction de la différence des vitesses de propagation de la fréquence de référence et de la fréquence en question. Les lois régissant ce phénomène ne sont pas bien définies ; en fait, ce phénomène peut être négligé dans la plupart des cas. Toutefois, dans tous les cas,des programmes appropriés peuvent être établis à partir de mesures spéciales dans un sondage ou déduits des caractéristiques des informations recueillies au cours de diverses expériences. Il va de soi que la présente invention n'a été décrite qu'à titre indicatif, mais nullement limitatif, et qutelle est susceptible de diverses variantes sans sortir de son cadre. TABLEAU I Fréquence 1 | Fréquence Fréquence émise -10 émise J il émise 412 Fi F2 F3 fréquence Fréquence Fréquence enregistrée 13 enregistrée 14 F3 l 15 Fi F2 F3 Filtre L W 16 | Filtre W 17 I Filtre | 18 LT |Pondération 19 | |Pondération | .'amplitude L-n > 22 Enre stre- 23 TABLEAU Il Fréquence Fréquence Fréquence émise ~ 10 émise il emise 12 Fi F2 F3 Fréquence Fréquence, réquence enregistrée 13 enregistrée r 14 tenregistrée r t5 Fi F2 F Filtre 16 Filtre )7 Filtre t8 ondération - - pondération -M Pondération 'amplitude r 19 d'amplitude r 20 | d'amplituderl Correction | | Correction Correction -tO temporeUe r 40 | temporelle empOB I 22 | Enregistremen + 23 RE\TEI I CATI ONS 1. Procédé d'exploration sismique dans lequel un groupe de trains d'ondes de fréquence constante est émis en direction du sol et les signaux résultants sont détectés en un point espacé du point d'émission et sont combinés, caractérisé en ce qu'il comprend une opération d'application d'une pondération, variable avec le temps, à l'amplitude des signaux détectés avant cette combinaison, pour tenir compte des variations des caractéristi- ques d'atténuation du sol en fonction de la fréquence. 2. Procédé selon la revendication 1, caractérisé en ce que les fréquences médianes des trains d'ondes forment une progression arithmétique. 3. Procédé selon la revendication 1, caractérisé en ce que les fréquences médianes des trains d'ondes forment une progression géométrique. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que les trains d'ondes émis en direction du sol sont modulés en amplitude. 5. Procédé selon l'une- quelconque des revendications 1 à 4, caractérisé en ce que tous les signaux détectés passent par un filtre à bande étroite destiné à laisser passer principalement la fréquence du train d'ondes émis correspondant et à éliminer les fréquences de bruit. J 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que le nombre de signaux combinés pour chaque fréquence varie en fonction-de la fréquence, de manière à corriger -les variations en fonction de la fréquence de l'atténuation due au sol et de l'intensité des trains d'ondes émis. 7.- Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que la largeur de bande du groupe de trains d'ondes est au moins égale à un octave environ. 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que les signaux détectés subissent, avant leur combinaison, une variation de phase qui est fonction du temps et de la fréquence, de manière à compenser les variations de la vitesse de propagation dans le sol en fonction de la fréquence, de manière que tous les signaux enregistrés soient maintenus en phase, quel que soit le temps de propagation. 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les trains d'ondes sont modulés en amplitude par une onde équivalant sensiblement à une demi-période d'une fonction cosinus carré.