FR2502794A1

FR2502794A1 - METHOD OF SEISMIC EXPLORATION OF A MEDIUM, ESPECIALLY OF GEOPHYSICAL PROSPECTION BY SEISMIC WAVES

Info

Publication number: FR2502794A1
Application number: FR8106082A
Authority: FR
Inventors: Maurice Barbier
Original assignee: GULF INTERSTATE GEOPHYSICAL
Current assignee: GULF INTERSTATE GEOPHYSICAL
Priority date: 1981-03-26
Filing date: 1981-03-26
Publication date: 1982-10-01
Also published as: GB2095835A; ES8303713A1; NL8200976A; IT1155671B; BR8201711A; IT8267384A0; DE3210581A1; JPS57173776A; AU8103682A; ES510784A0

Abstract

CE PROCEDE CONSISTE A EMETTRE UNE PREMIERE SEQUENCE DE N IMPULSIONS MECANIQUES SEPAREES PAR DES TEMPS EGAUX T SUIVIE PAR UN TEMPS D'ECOUTE AU MOINS EGAL AU TEMPS ALLER-RETOUR D'UNE IMPULSION DU SOL A LA COUCHE LA PLUS PROFONDE, A CORRELER LES SIGNAUX RECUS, A EMETTRE UNE DEUXIEME SEQUENCE DE N IMPULSIONS SEPAREES PAR DES TEMPS EGAUX T DIFFERENTS DE T, A PROCEDER COMME PRECEDEMMENT, A EMETTRE UN NOMBRE SUFFISANT DE SEQUENCES NI, TI, A SOMMER LES RESULTATS DES DIFFERENTES CORRELATIONS, LES NI, TI, ETANT TELS QUE LA SOMME DES FONCTIONS D'AUTOCORRELATION SE PRESENTE SOUS LA FORME D'UN PIC CENTRAL D'AMPLITUDE N N N... NI ENCADRE DE PICS SECONDAIRES AVEC UN RAPPORT SUPERIEUR AU RAPPORT DES AMPLITUDES DES SIGNAUX RECUS DANS DES INTERVALLES DE TEMPS CORRESPONDANT AUX INTERVALLES SEPARANT LE PIC PRINCIPAL DE CHACUN DES PICS SECONDAIRES.THIS PROCESS CONSISTS OF SENDING A FIRST SEQUENCE OF N MECHANICAL PULSES SEPARATED BY EQUAL TIMES FOLLOWED BY A LISTENING TIME AT LEAST EQUAL TO THE BACK AND RETURN TIME OF A PULSE FROM THE SOIL TO THE DEEP LAYER, CORRECTING THE SIGNALS RECEIVED, TO EMIT A SECOND SEQUENCE OF N PULSES SEPARATED BY TIME EQUAL T DIFFERENT FROM T, TO PROCEED AS PREVIOUSLY, TO EMIT A SUFFICIENT NUMBER OF SEQUENCES NI, TI, TO SUM THE RESULTS OF THE VARIOUS CORRELATIONS, THE NI, TI, SUCH AS THE SUM OF THE AUTOCORRELATION FUNCTIONS IS IN THE FORM OF A CENTRAL PIC OF AMPLITUDE NN N ... NOR FRAMED BY SECONDARY PEAKS WITH A RATIO GREATER THAN RATIO OF SIGNAL AMPLITUDES RECEIVED IN TIME INTERVALS CORRESPONDING TO INTERVALS SEPARATING THE MAIN PEAK FROM EACH OF THE SECONDARY PEAKS.

Description

1. La présente invention concerne un procédé d'exploration sismique d'un1. The present invention relates to a method of seismic exploration of a

milieu et plus particulièrementmiddle and more particularly

un procédé de prospection géophysique par ondes sismiques. a geophysical prospecting process using seismic waves.

Un procédé connu consiste à émettre une suite -d'impulsions d'amplitude sensiblement constante dans laquelle l'intervalle de temps entre deux impulsions successives est inférieur au temps de parcours T aller et retour d'une impulsion entre la surface du sol et l'interface géologique le plus profond auquel on s'intéresse, et dont la fonction d'autocorrélation se présente sous la forme d'un pic central encadré par des pics secondaires tel que le rapport entre ce pic central et les pics secondaires soit le plus grand possible. La valeur de ce rapport est très importante car c'est elle qui limite la profondeur de l'interface qu'il sera possible de détecter lorsque le procédé est mis en oeuvre. A known method consists in emitting a series of pulses of substantially constant amplitude in which the time interval between two successive pulses is less than the travel time T back and forth of a pulse between the ground surface and the deepest geological interface in which we are interested, and whose autocorrelation function is in the form of a central peak framed by secondary peaks such that the ratio between this central peak and the secondary peaks is as large as possible . The value of this ratio is very important because it is this which limits the depth of the interface that it will be possible to detect when the method is implemented.

Une façon d'obtenir un grand rapport entre l'am- One way to get a great relationship between am-

plitude du pic central et celles des pics secondaires consis- the central peak and the secondary peaks consist of

te à utiliser des intervalles de temps inégaux entre les -impulsions successives comme il est décrit dans le brevet français nO 1583239 (voir figure 1 annexée au présent mémoire descriptif). Suivant ce brevet, la suite d'impulsions émise a une durée tn supérieure à T et peut être représentée par une fonction du temps f (t) qui prend à certains instants -la valeur 1 et garde la valeur 0 le reste du temps, comme te to use unequal time intervals between successive pulses as described in French Patent No. 1583239 (see Figure 1 appended to this specification). According to this patent, the sequence of pulses emitted has a duration tn greater than T and can be represented by a function of time f (t) which takes at certain times the value 1 and keeps the value 0 the rest of the time, as

illustré sur cette figure 1, avec (ti + 1 - ti) (T. illustrated in this figure 1, with (ti + 1 - ti) (T.

Si on utilise des intervalles de temps t1, t2, t3... ti tn, on trouvera des pics secondaires à des temps qui sont donnés par la valeur de chacun des intervalles, la -somme de deux intervalles successifs, la somme de trois intervalles successifs, etc..., l'amplitude du pic secondaire à chacun des temps étant égale au nombre de fois que le même intervalle a été utilisé dans la séquence ou le nombre de fois que la somme d'intervalles successifs donnera la même valeur (voir figure 2 qui illustre la fonction d'auto-corrélation de la suite d'impulsions émise avec un maximum d'amplitude n s'il y a "n" impulsions et des maxima secondaires aux temps: + t1; + t2***. + tn* En théorie, on peut donc obtenir de très grands If we use time intervals t1, t2, t3 ... ti tn, we will find secondary peaks at times which are given by the value of each of the intervals, the sum of two successive intervals, the sum of three intervals successive, etc ..., the amplitude of the secondary peak at each time being equal to the number of times that the same interval was used in the sequence or the number of times that the sum of successive intervals will give the same value ( see figure 2 which illustrates the auto-correlation function of the sequence of pulses emitted with a maximum amplitude n if there are "n" pulses and secondary maxima at times: + t1; + t2 *** . + tn * In theory, we can therefore obtain very large

rapports en utilisant un grand nombre d'intervalles diffé- reports using a large number of different intervals

rents mais en pratique on est très vite limité à un nombre réduit d'intervalles de temps possibles par l'emploi de la source sismique qui doit transmettre les impulsions sismiques dans le sol. En effet, une source est construite pour fonctionner à une certaine fréquence de répétition pour laquelle elle fournira sa puissance nominale. Si on augmente cette fréquence de répétition, l'amplitude des impulsions sismiques produites va diminuer et si on diminue cette fréquence, la source est alors mal employée et les résultats sismiques qu'on obtiendra seront moins bons puisqu'on aura dans le temps imparti à la transmission envoyé moins d'énergie dans le sol. On peut donner comme exemple de sources sismiques un canon à air comprimé pour lequel la cadence de répétition dépendra de la taille du compresseur qui l'alimente ou un marteau hydraulique pour lequel la cadence de répétition dépendra des caractéristiques de la pompe à huile. A signaler que la méthode actuelle qui oblige quelquefois à passer d'un intervalle de temps très court entre deux émissions à un intervalle de temps long n'est pas satisfaisant du point de vue de l'utilisation de la source, lui faisant subir des à-coups auxquels l'inertie propre à tout système physique l'empoche de répondre avec but in practice one is very quickly limited to a reduced number of possible time intervals by the use of the seismic source which must transmit the seismic pulses into the ground. Indeed, a source is built to operate at a certain repetition frequency for which it will provide its nominal power. If we increase this repetition frequency, the amplitude of the seismic pulses produced will decrease and if we decrease this frequency, the source is then misused and the seismic results that we will obtain will be less good since we will have in the time allowed to the transmission sent less energy into the ground. An example of a seismic source is an air cannon for which the repetition rate will depend on the size of the compressor that powers it, or a hydraulic hammer for which the repetition rate will depend on the characteristics of the oil pump. Note that the current method which sometimes requires passing from a very short time interval between two emissions to a long time interval is not satisfactory from the point of view of the use of the source, causing it to undergo - blows to which the inertia proper to any physical system prevents it from responding with

une pleine efficacité.full efficiency.

Pour une source donnée, on est donc limité à un nombre d'intervalles de temps compris entre deux valeurs t1 et t2, t1 étant l'intervalle de temps minimum et t, l'intervalle de temps maximum. Entre t1 et t2, le nombre de valeurs distinctes dont on dispose est fini, car -l'enregistrement des signaux s'effectue avec un pas 3. d'échantillonnage et il n'y a aucune différence physique entre deux intervalles dont la différence est inférieure à la valeur du pas d'échantillonnage utilisé. Pour mieux faire comprendre la différence entre l'état de la technique actuelle et le procédé selon l'invention, on donnera For a given source, we are therefore limited to a number of time intervals between two values t1 and t2, t1 being the minimum time interval and t, the maximum time interval. Between t1 and t2, the number of distinct values available is finite, because - the recording of the signals is carried out with a sampling step 3. and there is no physical difference between two intervals whose difference is lower than the value of the sampling step used. To better understand the difference between the state of the art and the process according to the invention, we will give

ci-après un exemple chiffré.below is an encrypted example.

On supposera que l'on dispose d'une source dont on peut faire varier efficacement la cadence d'émission entre valeurs distinctes caractérisées par 10 intervalles de temps entre émissions successives t1 à atIO Il s'agit d'un nombre très réduit qui a été choisi pour faciliter les explications mais le principe resterait le même si ce We will assume that we have a source whose emission rate can be effectively varied between distinct values characterized by 10 time intervals between successive emissions t1 to atIO. This is a very small number which has been chosen to facilitate explanations but the principle would remain the same if it

nombre était plus élevé.number was higher.

Dans l'état de la technique actuelle, avec ces 10 valeurs on peut construire une suite d'impulsions telle que la fonction d'autocorrélation se présente sous la forme d'un pic central d'amplitude 10 et des pics secondaires d'amplitude 1 à des temps égaux à t1, t2... t10, puis (t1 + t2) etc... (voir figure 3). Si ce nombre de 10 impulsions n'est pas suffisant pour émettre assez d'énergie dans le sol, il faut émettre d'autres impulsions en utilisant à nouveau les mêmes intervalles t1 à t10. L'aihplitude des pics secondaires va alors augmenter. Par exemple, si on a besoin d'émettre 100 impulsions dans le sol, le mieux qu'on puisse faire est d'utiliser 10 fois chacun des intervalles de temps, pour une amplitude du pic central de 100 et une amplitude de chacun des pics secondaires de 9 (voir figure 4). Du point de vue du rapport de l'amplitude du pic central à l'amplitude des pics secondaires, il n'y a - eu pratiquement aucune amélioration. La technique actuelle est donc limitée par l'amplitude de ces pics secondaires In the current state of the art, with these 10 values, it is possible to construct a sequence of pulses such that the autocorrelation function is in the form of a central peak of amplitude 10 and secondary peaks of amplitude 1 at times equal to t1, t2 ... t10, then (t1 + t2) etc ... (see Figure 3). If this number of 10 pulses is not sufficient to emit enough energy into the ground, other impulses must be emitted using the same intervals t1 to t10 again. The amplitude of the secondary peaks will then increase. For example, if you need to emit 100 pulses into the ground, the best you can do is use 10 times each of the time intervals, for an amplitude of the central peak of 100 and an amplitude of each of the peaks secondary 9 (see Figure 4). In terms of the ratio of the amplitude of the central peak to the amplitude of the secondary peaks, there has been - practically no improvement. The current technique is therefore limited by the amplitude of these secondary peaks

qu'on appelle "bruit de corrélation". called "correlation noise".

Pour décrire le procédé selon l'invention on reprendra les chiffres de l'exemple précédent, à savoir 100 impulsions à émettre avec 10 intervalles de temps différents, et on se référera aux figures 5 à 9 qui sont 4. des diagrammes des signaux entrant en jeu et de leurs phases To describe the process according to the invention, the figures of the previous example will be used, namely 100 pulses to be transmitted with 10 different time intervals, and reference will be made to FIGS. 5 to 9 which are 4. diagrams of the signals entering game and their phases

de traitements.of treatments.

Ce procédé consiste alors à émettre d'abord une suite de 10 impulsions équidistantes l'intervalle de temps entre les impulsions étant t1, à enregistrer les signaux reçus par les détecteurs pendant une durée égale au temps séparant la première impulsion de la dernière augmentée du temps de trajet aller et retour T de la dernière impulsion entre la surface et l'interface le plus profond auquel on s'intéresse, à corréler les signaux reçus avec la suite représentant les temps d'émission des 10 impulsions, à transmettre une nouvelle suite de 10 impulsions équidistantes, l'intervalle de temps étant t2, puis à procéder comme pour la première séquence, à continuer à émettre de nouvelles suites d'impulsions caractérisées par le même nombre 10 mais des intervalles de temps t3 puis t4 etc... jusqu'à This method then consists in first emitting a series of 10 equidistant pulses, the time interval between the pulses being t1, in recording the signals received by the detectors for a duration equal to the time separating the first pulse from the last one increased by time. back and forth path T of the last pulse between the surface and the deepest interface in which we are interested, to correlate the signals received with the sequence representing the transmission times of the 10 pulses, to transmit a new sequence of 10 equidistant pulses, the time interval being t2, then proceed as for the first sequence, to continue to emit new sequences of pulses characterized by the same number 10 but time intervals t3 then t4 etc ... until 'at

t10, enfin à sommer les résultats des 10 opérations succes- t10, finally summing up the results of the 10 successful operations

sives ainsi réalisées.sives thus produced.

Le détail de chacune des opérations ci-dessus est donné-ci-après: The details of each of the above operations are given below:

La figure 5 représente l'une des 10 suites d'impul- Figure 5 shows one of the 10 pulse sequences.

sions que l'on désignera par Ci (t) [illustré sur le schéma (a)] ainsi que les signaux reçus [schéma (b)] pendant un temps égal à la durée de la suite émise augmentée du -temps T précédemment défini. Si chaque impulsion de la suite Ci (t) a produit dans le sol un signal émis s(t) dont la forme dépend à la fois de la source choisie et de la nature de la surface du sol, on peut écrire que le signal sions that will be designated by Ci (t) [illustrated in diagram (a)] as well as the signals received [diagram (b)] for a time equal to the duration of the transmitted sequence increased by the time T previously defined. If each pulse of the sequence Ci (t) has produced in the soil an emitted signal s (t) whose shape depends on both the chosen source and the nature of the soil surface, we can write that the signal

transmis s'exprime par l'expression Ci(t) a s (t) (1). transmitted is expressed by the expression Ci (t) a s (t) (1).

Par ailleurs, si on appelle R.C.log(t) le signal reçu par les détecteurs de surface lorsqu'une seule impulsion est transmise, le signal reçu lorsqu'on a transmis la suite Furthermore, if the signal received by the surface detectors when a single pulse is transmitted is called R.C.log (t), the signal received when the sequence is transmitted

Ci(t) s'exprime par: Ci(t)ys(t)xR.C.log(t) (2). Ci (t) is expressed by: Ci (t) ys (t) xR.C.log (t) (2).

Il a été dit que le signal reçu était alors corrélé avec la suite émise, c'est-à-dire avec Ci (t). Le résultat de la i 5. corrélation peut alors s'écrire AFC Ci(t)xs(t)XR.C.log(t) (3), AFC Ci (t) It has been said that the signal received was then correlated with the sequence sent, that is to say with Ci (t). The result of the i 5. correlation can then be written AFC Ci (t) xs (t) XR.C.log (t) (3), AFC Ci (t)

étant la fonction d'autocorrélation de la séquence C. (t). being the autocorrelation function of the sequence C. (t).

La figure 6 représente cette fonction d'autocorré- Figure 6 represents this function of autocorré-

lation de Ci(t) (AFC) dans sa totalité, donc une durée égale lation of Ci (t) (AFC) in its entirety, therefore an equal duration

à 9 t. ainsi que la durée de la fonction R.C.log(t) recher- at 9 t. as well as the duration of the R.C.log (t) function sought-

i chée dont on sait seulement que la durée T est supérieure 3 ti. L'origine des temps à prendre en considération est le maximum de la fonction d'autocorrélation qui doit être calé i chée which we only know that the duration T is greater than 3 ti. The origin of the times to be taken into account is the maximum of the autocorrelation function which must be calibrated

sur le temps zéro de la fonction R.C.log(t). on time zero of the R.C.log (t) function.

La figure 7 représente les 10 fonctions d'auto- Figure 7 shows the 10 auto functions

corrélation de C1(t) à C1 (t) calées sur leur maximum ainsi correlation of C1 (t) to C1 (t) set on their maximum as well

que leur somme qui doit être effectuée selon la description that their sum to be made as described

de la nouvelle méthode. On remarque que le résultat final of the new method. We notice that the final result

est un pic central d'amplitude 100 encadré de pics secon- is a central peak of amplitude 100 framed by second peaks

daires d'amplitude 9. A ce stade le résultat est identique amplitude daries 9. At this stage the result is identical

à celui obtenu par la technique précédente qui a été illus- to that obtained by the previous technique which was illus-

tré par la figure 4.shown in Figure 4.

Le premier avantage du procédé selon l'invention réside dans une utilisation plus efficace de la source The first advantage of the process according to the invention lies in a more efficient use of the source

sismique qui fonctionne à plusieurs cadences fixes successi- seismic which operates at several successive fixed rates

ves plutôt que d'une façon désordonnée. Un deuxième avantage est aussi la possibilité d'atténuer les bruits de surface rather than in a haphazard way. A second advantage is also the possibility of attenuating surface noise

lorsqu'on utilise au moins deux sources. when using at least two sources.

On supposera que l'on utilise deux sources synchrones S1et S2 espacées d'une distance di le long de la ligne reliant les détecteurs G et situées toutes les deux du même côté de ces détecteurs. On supposera également que l'émission de l'une de ces sources à la surface du sol produise une onde de surface se propageant horizontalement à la vitesse V (figure 8'o S1 et S2 sont les sources d'émission, V est la vitesse de propagation de l'onde de surface, ti est l'intervalle de temps entre les émissions des deux sources synchrones, di est la distance entre les deux sources It will be assumed that two synchronous sources S1 and S2 are used, spaced apart by a distance di along the line connecting the detectors G and both located on the same side of these detectors. We will also suppose that the emission of one of these sources on the surface of the ground produces a surface wave propagating horizontally at the speed V (figure 8'o S1 and S2 are the sources of emission, V is the speed propagation of the surface wave, ti is the time interval between the emissions of the two synchronous sources, di is the distance between the two sources

lorsqu'elles émettent à la cadence ti). when they transmit at the rate ti).

6. Si on choisit la distance di entre les deux sources telle que ti l'intervalle de temps de la séquence émise soit égal à di/V, on constate que la deuxième émission de la source S2 se produira au moment o la première de S1 passera par -S2, il en résultera un renforcement de cette onde. Par contre, si on choisit d. tel que (k+1)d./V = 3t./2, il en résultera une atténuation du fait que des émissions de S1 et S2 se trouvent en opposition de phase. En ajustant la distance entre les sources pour chacun des intervalles de -temps utilisés l'effet se produira pour chacune des séquences d'émission alors que ceci n'est pas possible avec la technique actuelle qui utilise des intervalles de temps 6. If one chooses the distance di between the two sources such that ti the time interval of the transmitted sequence is equal to di / V, one notes that the second emission of the source S2 will occur at the moment when the first of S1 will go through -S2, it will result in a strengthening of this wave. On the other hand, if we choose d. such that (k + 1) d./V = 3t./2, it will result in an attenuation of the fact that emissions of S1 and S2 are in phase opposition. By adjusting the distance between the sources for each of the time intervals used the effect will occur for each of the emission sequences while this is not possible with the current technique which uses time intervals

différents pour la même distance entre les sources. different for the same distance between the sources.

Un troisième avantage du nouveau procédé est -de pouvoir obtenir avec un nombre limité d'intervalles de temps différent un rapport de l'amplitude du pic central A third advantage of the new process is to be able to obtain, with a limited number of different time intervals, a ratio of the amplitude of the central peak

sur le bruit de corrélation aussi grand qu'on le désire. on the correlation noise as large as desired.

On notera une différence essentielle entre le procédé selon l'invention et la technique connue après une émission de -10 impulsions. Dans le cas de l'invention le bruit de corrélation est constitué en un certain nombre de pics secondaires de relativement forte amplitude mais apparaissant à des temps bien définis ti, 2ti, 3 ti, etc... alors que dans la technique actuelle l'amplitude des pics secondaires 25.est limitée à 1 mais ils apparaissent à des temps très différents t1, t2, t3.... t10, (t1 + t2), Puis (t2 + 3 etc... On peut tirer parti de cette différence en appliquant An essential difference will be noted between the method according to the invention and the known technique after an emission of -10 pulses. In the case of the invention, the correlation noise consists of a number of secondary peaks of relatively high amplitude but appearing at well defined times ti, 2ti, 3 ti, etc ... whereas in the current technique the amplitude of secondary peaks 25.is limited to 1 but they appear at very different times t1, t2, t3 .... t10, (t1 + t2), Then (t2 + 3 etc ... We can take advantage of this difference in applying

après émission de chaque série de 10 impulsions et corréla- after emission of each series of 10 pulses and correlation

tion, un opérateur de filtrage à 3 points constitué d'un pic -d'amplitude -1 suivi d'un pic d'amplitude 2 au temps ti et un pic d'amplitude -1 au temps 2t.. La figure 9 représente de gauche à droite: la fonction AFC Ci (t), l'opérateur de filtrage à trois points précédemment défini et le résultat du filtrage. On voit que l'amplitude des pics secondaires -a été annulée alors que celle du pic central est encore égale à 2. En théorie, on pourrait donc se contenter d'une seule séquence, tous les pics secondaires ayant 7. disparu. En pratique, l'atténuation pourrait être moins efficace s'il se produit de légères -variations des ti lors de l'émission d'une même séquence et on préférera répéter la tion, a 3-point filtering operator consisting of a peak-of amplitude -1 followed by a peak of amplitude 2 at time ti and a peak of amplitude -1 at time 2t. Figure 9 represents left to right: the AFC Ci (t) function, the previously defined three-point filtering operator and the result of the filtering. We see that the amplitude of the secondary peaks has been canceled while that of the central peak is still equal to 2. In theory, we could therefore be satisfied with a single sequence, all the secondary peaks having 7. disappeared. In practice, the attenuation could be less effective if slight ti -variations occur during the emission of the same sequence and it is preferable to repeat the

même opération avec plusieurs ti de valeurs différentes. same operation with several ti of different values.

Il va de soi que la présente invention n'a été décrite qu'à titre purement explicatif et nullement limitatif et que toutes modifications utiles pourront y être apportées sans sortir du cadre de l'invention, tel que défini par les It goes without saying that the present invention has been described for purely explanatory and in no way limitative and that any useful modifications may be made thereto without departing from the scope of the invention, as defined by the

revendications ci-après.claims below.

8.8.

Claims

1. - Method for exploring a medium by coded emission of mechanical pulses, characterized in that

that it consists in emitting a first sequence of impulse NI

mechanical ions separated by equal time intervals t1 followed by a listening time at least equal to the round trip time of one of the pulses from the ground surface to the deepest layer that we want to study, at

correlate signals received by mechanical wave detectors

pics with the sequence transmitted, to issue a second

sequence of N2 mechanical pulses separated by inter-

equal time t2 but different from those of the first sequence, to proceed then as for the first sequence, to continue to transmit then a sufficient number of sequences themselves characterized by numbers

of impulses Ni and at different intervals ti, at som-

mer the results of the various correlations, the Ni and ti having been chosen in such a way that the sum of the autocorrelation functions of each of the sequences emitted is in the form of a central peak of amplitude N1 + N2 + N3 +. .. Nor framed by secondary peaks such that the ratio between this central peak and the secondary peaks is greater than the ratio of the amplitudes of the reflected signals received in time intervals corresponding to the given time intervals separating the main peak from each of the secondary peaks .

2. - Method according to claim 1, characterized in that after the correlation of the signals received with each sequence, the result of this correlation is convolved with a 3-point operator consisting of a peak of amplitude -1 followed by a amplitude peak 2 at a time equal to the time interval between two successive transmissions and an amplitude peak -1 at a time equal to twice the interval

of time between two successive transmissions.

9.

3. - Method according to claim 1 or 2, characterized in that the pulses are emitted by at least two sources whose distance between them is variable

from one broadcast sequence to another.

4. - Method according to claims 1, 2 or 3,

characterized in that the distance between the emission sources is chosen in each sequence according to the

speed of propagation of surface waves.