Brevet FR2941055A1 — publié le 20100716

PROCEDE D'ACQUISITION DE DONNEES VIBROSISMIQUES CONCERNANT UNE ZONE DU SOUS-SOL, ET PROCEDE D'EXPLORATION SISMIQUE INCLUANT UN TEL PROCEDE L'invention concerne l'exploration sismique d'une zone du sous-sol terrestre. L'utilisation des ondes sismiques générées artificiellement pour la reconnaissance des couches géologiques est utilisée depuis plus de 50 ans. Les ondes les plus utilisées sont de loin les ondes réfléchies et plus précisément les ondes de compression réfléchies. Lors d'opérations de prospection sismique, un équipement vibrateur génère un signal vibrosismique qui se propage notamment sous la forme d'une onde qui se réfléchit sur des interfaces des couches géologiques. Ces ondes sont reçues par des géophones qui transforment le mouvement du sol en un signal électrique enregistré au moyen d'un équipement d'enregistrement. L'analyse des temps d'arrivée et des amplitudes de ces ondes permet de construire une représentation des couches géologiques sur lesquelles les ondes se sont réfléchies. L'une des complications de la méthode sismique terrestre est que les ondes réfléchies ne sont pas les seules à être enregistrées. En particulier, l'onde associée à la plus grande proportion d'énergie est l'onde de surface qui se propage le long de la surface du sol. Cette onde est similaire à celle qui se propage à la surface de l'eau quand on y fait tomber par exemple un caillou. Elle est communément appelée �ground roll �, mot qui désigne (en anglais) le mouvement de rotation des particules à la surface du sol. L'onde de surface ou ground roll interfère avec les ondes réfléchies et en brouille l'image. Il y a d'autres ondes parasites générées mais elles sont en général moins gênantes. Dans la suite, les termes � ondes de surface � et � ondes réfléchies � seront utilisés pour désigner les parties d'un sismogramme représentant les ondes respectives correspondantes. Il existe actuellement deux techniques principales pour séparer les ondes de surface dans les données sismiques enregistrées. La technique la plus répandue utilise un filtrage en vitesse apparente qui permet de séparer les ondes de surface lentes des ondes 5 réfléchies plus rapides. La seconde technique est celle du filtre de polarisation. Les ondes de surface s'accompagnent d'un mouvement circulaire de la surface du sol alors que les ondes réfléchies s'accompagnent d'un mouvement linéaire. C'est cette différence que met à profit le filtre de polarisation. 10 Il est par ailleurs connu ( Martin, J.E., Jack, I.G. "The behaviour of a seismic vibrator using different phase control methods and drive levels" First Break, November 1990) que la variation du niveau de la force générée par un vibrateur se traduit par une variation des temps de propagation des ondes réfléchies. Cette variation de temps est en général interprétée 15 comme une variation de vitesses de propagation dans la zone compressée sous la plaque du vibrateur. Il a été observé que cette variation de temps n'était pas identique pour tous les types d'ondes. En particulier, pour les ondes réfléchies utilisées pour obtenir la quasi-totalité des images en sismique pétrolière, 20 cette variation est souvent très faible alors que pour les ondes de surface elle peut être beaucoup plus importante. La présente invention tire parti de cette observation. Selon un premier aspect, l'invention a pour objet un procédé d'acquisition de données vibrosismiques concernant une zone du sous-sol 25 fournissant les informations permettant de séparer les composantes d'ondes de surface contenues dans les données vibrosismiques. Il est prévu selon ce premier aspect un procédé d'acquisition de données vibrosismiques concernant une zone du sous-sol, comprenant les étapes suivantes : 30 - activer une source vibrosismique pour qu'elle émette un premier signal vibrosismique présentant une première amplitude, - activer la source pour qu'elle émette un deuxième signal vibrosismique présentant une deuxième amplitude, différente de la première amplitude, - enregistrer des premières données correspondant au premier 5 signal après propagation dans le milieu à explorer, - enregistrer des deuxièmes données correspondant au deuxième signal après propagation dans le milieu à explorer. Le procédé proposé tire parti de la différence de temps de propagation liée à la variation du niveau de force générée par la source. 10 Grâce à cette différence, une estimation des ondes de surface peut être obtenue à partir des premières et deuxièmes données en vue de leur filtrage. Le procédé peut également être utilisé pour séparer les ondes de compression converties en ondes de cisaillement dans le cas où celles-ci 15 sont enregistrées au moyen de géophones à 3 composantes. En effet, le phénomène est associé à la source et affecte indifféremment les ondes de compression qui se réfléchissent et celles qui se convertissent en ondes de cisaillement. De même, le procédé peut être appliqué aux ondes de cisaillement 20 émises par une source de cisaillement. Enfin, le procédé proposé s'applique directement aux ondes de surface diffractées. Selon une première possibilité, les premier et deuxième signaux émis par la source ne diffèrent que par leur amplitude. 25 Selon une deuxième possibilité, les premier et deuxième signaux présentent une forme et une durée différentes. Dans tous les cas, les données enregistrées à partir des premier et deuxième signaux contiennent chacune des ondes réfléchies et des ondes de surface. Celles-ci ne diffèrent que par leurs amplitudes et par leur temps 30 de propagation. Selon un deuxième aspect, l'invention a pour objet un procédé d'exploration sismique d'une zone du sous-sol permettant de séparer les composantes d'ondes de surface contenues dans des données vibrosismiques. Il est prévu selon ce deuxième aspect un procédé d'exploration sismique d'une zone du sous-sol, comprenant les étapes suivantes : - activer une source vibrosismique pour qu'elle émette un premier signal vibrosismique présentant une première amplitude, - activer la source pour qu'elle émette un deuxième signal vibrosismique présentant une deuxième amplitude, différente de la première amplitude, - enregistrer des premières données correspondant au premier signal après propagation dans le milieu à explorer, - enregistrer des deuxièmes données correspondant au deuxième signal après propagation dans le milieu à explorer, - obtenir une estimation des ondes de surface contenues dans ces 15 données à partir des premières et deuxièmes données enregistrées, et - soustraire cette estimation aux premières données. De façon appropriée, l'étape d'obtention d'une estimation des ondes de surfaces comprend les sous-étapes suivantes : 20 - comparer les premières et deuxièmes données enregistrées. - transformer les deuxièmes données pour égaliser la composante d'ondes réfléchies contenue dans les deuxièmes données avec la composante d'ondes réfléchies contenue dans les premières données, - soustraire les deuxièmes données transformées aux premières 25 données de manière à éliminer la composante d'ondes réfléchies afin de ne conserver que les composantes d'ondes de surface, - estimer les ondes de surface à partir de la soustraction ci-dessus. La sous-étape de transformation peut comprendre 30 - multiplier les deuxièmes données par un coefficient pour égaliser l'amplitude des composantes d'ondes réfléchies contenues dans les premières et deuxièmes données, - décaler les deuxièmes données d'un temps qui ramène les composantes d'ondes réfléchies contenues dans les premières et deuxièmes données à une référence temporelle commune. Dans ce cas, la sous-étape de comparaison consiste de façon 5 appropriée à mesurer le temps de décalage entre les deuxièmes et premières données. La sous-étape d'estimation des ondes de surface consiste de façon appropriée en une simple intégration. 10 Alternativement, la sous-étape de transformation peut comprendre : - appliquer aux deuxièmes données un opérateur de passage. Dans ce cas, la sous-étape de comparaison consiste à déterminer un opérateur de passage entre les premières et deuxièmes données. Et la sous-étape d'estimation de la composante d'ondes de 15 surfaces est réalisée par un calcul et une application d'un filtre inverse. L'opérateur de passage peut être déterminé à partir de zones des premières et deuxièmes données dans lesquelles les ondes réfléchies et les ondes de surface n'interfèrent sensiblement pas entre elles. L'opérateur de passage prend en compte de façon appropriée 20 plusieurs paramètres incluant la fréquence et des nombres d'onde dans des directions horizontales. PRESENTATION DES FIGURES - La figure 1 représente de manière schématique un dispositif 25 d'émission et de réception d'ondes vibrosismiques, - la figure 2 représente de manière schématique un premier signal (signal faible) et un deuxième signal (signal fort), tels qu'ils sont émis par le dispositif de la figure 1, - la figure 3 est un diagramme représentant de manière 30 schématique des enregistrements obtenus par le dispositif de la figure 1 après propagation du deuxième signal dans le sous-sol, - la figure 4 est un diagramme représentant de manière schématique des enregistrements obtenus par le dispositif de la figure 1 après propagation du premier signal dans le sous-sol et égalisation, - la figure 5 est un diagramme obtenu par soustraction des diagrammes des figures 3 et 4, - la figure 6 est un diagramme obtenu après intégration du diagramme de la figure 5, correspondant à une estimation des signaux de surface, - la figure 7 est un diagramme obtenu après soustraction des 10 diagrammes des figures 3 et 6, correspondant à un enregistrement dans lequel les signaux de surface ont été éliminés, - la figure 8 représente de manière schématique des étapes d'un procédé d'exploration sismique conforme à un premier mode de réalisation de l'invention, 15 - la figure 9 représente de manière schématique des étapes d'un procédé d'exploration sismique conforme à un deuxième mode de réalisation de l'invention. DESCRIPTION DETAILLEE 20 La figure 1 représente de manière schématique un dispositif d'émission et de réception d'ondes vibrosismiques destiné à l'exploration sismique en milieu terrestre. Le dispositif comprend une source constituée d'un vibrateur 1 apte à générer un signal sismique, un ensemble de récepteurs 2 (ou géophones) 25 aptes à recevoir un signal sismique et à le convertir en un signal électrique et un enregistreur 3 apte à enregistrer les signaux électriques générés par les récepteurs. La source 1, les récepteurs 2 et l'enregistreur 3 sont positionnés à la surface du sol 5. La figure 1 représente un seul vibrateur mais il est entendu que la source peut être composée de plusieurs 30 vibrateurs, comme il est bien connu de l'homme de métier. En fonctionnement, la source 1 est commandée pour générer un signal vibrosismique. Celui-ci se propage d'une part à la surface du sol, sous la forme d'ondes 4 de surface, et d'autre part dans le sous-sol, sous la forme d'ondes 6 qui se réfléchissent sur une interface 7 entre deux couches géologiques. Chaque récepteur 2 reçoit à la fois une onde 4 de surface et une onde 6 réfléchie et les convertit en un signal électrique dans laquelle sont superposées la composante correspondant à l'onde réfléchie et celle qui correspond à l'onde de surface, qui est indésirable et qu'il s'agit de filtrer. Premier mode de réalisation : modèle simplifié Un modèle simple consiste à décrire un sismogramme SR(t,d) issu d'un signal émis par la source vibrosismique S avec une force d et observé par le récepteur R en fonction du temps t comme la somme de deux composantes associées respectivement aux ondes réfléchies P(t,d) et aux ondes de surface G(t, d) . SR(t,d)=P(t,d)+G(t,d) [1] On utilise la variable d pour désigner une force de façon à éviter une confusion avec la fréquence. La variable d désigne en réalité le � drive �, c'est-à-dire le pourcentage de la force maximale disponible pour un type de vibrateur donné. C'est le drive qui est utilisé par les géophysiciens de terrain pour modifier la force appliquée sur le sol. Dans ce modèle, la variation de la force de b (bas) à h (haut) se traduit par une variation d'amplitude identique et une variation de temps différente pour chacune des deux composantes : SR(t, h) = P(t, h) + G(t, h) [2] Avec : SR(t, b) = P(t, b) + G(t, b) [3] PU, b) = h ù.PU +t p, h) [4] G(t,b) = h .G(t+tiG,h) [5] où tp désigne le retard de la composante basse amplitude par rapport à la composante haute amplitude, pour les composantes associées aux ondes 30 réfléchies,25 et tiG désigne le retard de la composante basse amplitude par rapport à la composante haute amplitude, pour les composantes associées aux ondes de surface. Par substitution, la relation [3] devient : SR(t, b) = h •[P(t+tip,h)+G(t+tiG,b)] Donc : SR(t ùtip, b) = h • [P(t, h)+G(t+tiG ùtip, b)] Par différence avec la relation [2], on obtient : SR(t ùtip, b)ùhSR(t, h)=-•[G(t+tiGùtip, b)ùG(t,h)] [8] Ce qui peut également s'écrire : SR(t ùtip,b)ù h SR(t,h) _ (tiG ùtip). dG(t,h) [9] Et par intégration permet de déterminer l'onde de surface : G(t, h) = 1 -tip) J [SR(t ùti p, b) ù b SR(t, h) dt [10] (tiG h La figure 8 représente de manière schématique des étapes d'un procédé d'exploration sismique comprenant un procédé d'acquisition vibrosismique conforme à l'invention, composé des étapes 101-104, et un procédé de filtrage des ondes de surface, selon un premier mode de réalisation (étapes 105-109). Dans une première étape 101, la source 1 constituée d'un vibrateur 20 est commandée pour générer un premier signal vibrosismique S~ (signal fort) présentant une première amplitude. Cette première amplitude est de façon appropriée l'amplitude optimale du point de vue de l'exploitation des signaux sismiques pour obtenir une image du sous-sol. A cet effet, le vibrateur est réglé pour appliquer sur le sol des vibrations avec un premier 25 niveau de force h (niveau haut), par exemple avec un drive d de 80%. Dans une deuxième étape 102, l'enregistreur 3 enregistre des premières données. [6] [7] Les premières données enregistrées sont représentées sur la figure 3. Ces premières données comprennent un ensemble de signaux (ou traces) enregistrées en fonction du temps et résultant de la propagation du premier signal sismique, chaque signal ayant été généré par un récepteur 2. Les premières données comprennent des données correspondant à des ondes de surface 4 et des données correspondant à des ondes réfléchies 6. Dans une troisième étape 103, le vibrateur 1 est commandé pour générer un deuxième signal vibrosismique S2 (signal faible) présentant une deuxième amplitude. A cet effet, le vibrateur est réglé pour appliquer sur le sol des vibrations avec un deuxième niveau de force b (niveau bas), par exemple avec un drive d de 20%. Le deuxième signal S2 est identique au premier signal SI , excepté qu'il présente une amplitude inférieure à celle du signal SI , comme cela est représenté sur la figure 2. Il est à noter que les termes � premier � et � deuxième � ne sont utilisés ici que pour distinguer d'une part les signaux émis au niveau haut et d'autre part les signaux émis au niveau bas, ainsi que les données respectivement enregistrées à la suite de ces émissions. Les termes � premier � et � deuxième � n'indiquent en rien un ordre chronologique. L'ordre dans lequel ont lieu les deux émissions et par suite les enregistrements est indifférent pour atteindre les objectifs de l'invention. L'émission d'un signal haut peut précéder celle d'un signal faible, ou l'inverse. Dans une quatrième étape 104, l'enregistreur 3 enregistre des 25 deuxièmes données. Dans une cinquième étape 105, les retards tp et tiG entre les composantes P(t,h) et P(t,b) et entre les composantes G(t,h) et G(t,b) sont déterminés à partir des premières et deuxièmes données. Dans une sixième étape 106, les deuxièmes données sont traitées 30 de la manière suivante : les amplitudes des traces sont multipliées par un facteur b et les traces sont décalées dans le temps de tiG ûtip . La figure 4 représente les deuxièmes données ainsi traitées. Cette sixième étape a pour effet de transformer les deuxièmes données de sorte que les composantes P d'onde réfléchie contenue dans les deuxièmes données soient superposables aux composantes P d'onde réfléchie contenues dans les premières données. Autrement dit, les composantes P d'onde réfléchie sont ramenées à une référence commune dans les premières et deuxièmes données. Dans une septième étape 107, les deuxièmes données ainsi 10 traitées ou transformées (figure 4) sont soustraites aux premières données (figure 3). Cette septième étape a pour effet d'éliminer les composantes P d'onde réfléchie. La figure 5 représente de manière schématique les données soustraites obtenues par l'étape de soustraction. On constate sur 15 cette figure que les ondes réfléchies ont disparu. Dans une huitième étape 108, les données soustraites sont intégrées par rapport au temps de manière à obtenir des données d'estimation GE des composantes d'onde de surface (ground roll). La figure 6 représente de manière schématique les données 20 d'estimation GE. Dans une neuvième étape 109, les données d'estimation du ground roll (figure 6) sont soustraites aux premières données (figure 3), de manière à éliminer les composantes d'onde de surface. Cette neuvième étape permet d'obtenir une estimation PE des 25 données correspondant aux ondes réfléchies. Le ground roll a sensiblement disparu. Le modèle correspondant au premier mode de réalisation décrit ci- dessus, basé sur une description de l'effet d'une variation de la force appliquée par les équations [4] et [5], est un modèle simplifié, qui ne prend 30 pas en compte d'autres types d'ondes présents dans les enregistrements, en particulier le bruit ambiant (non généré par la source), et certaines variables, telles que la fréquence, l'azimut et l'inclinaison des rayons qui peuvent avoir une incidence significative sur cet effet. Pour prendre en compte ces différentes variables, un deuxième mode de réalisation peut être utilisé, basé sur un modèle plus complexe. 5 Deuxième mode de réalisation : modèle complexe Pour prendre en compte l'ensemble des variables, il est possible de décrire le phénomène dans le domaine de Fourier {f, kx, ky}où f représente la fréquence, et kx et ky les nombres d'onde dans des directions 10 horizontales X et Y. De plus, dans ce domaine, l'effet du bruit ambiant est facile à contrôler et à réduire. Pour simplifier les notations on se contentera d'indiquer la variable f à la place de f, kx, ky . 15 En prenant des notations identiques, les équations [2] et [3] deviennent : SR(f, h) = P(f, h)+ G(f, h) [12] SR(f, b) = P(f, b)+ G(f, b) [13] A ce stade, on peut calculer un opérateur de passage 01, par 20 exemple par division spectrale : O1(f) = P(f' h) [14] P(f, b) L'opérateur le plus général dépend de f, kx, ky , le plus simple est un scalaire. L'opérateur 01 est appliqué aux données SR(f,b) enregistrées 25 avec une force basse b. SR2(f, b) = SR(f, b) *01(f) [15] En soustrayant la relation [14] de la relation [11], on obtient une estimation filtrée du ground roll : SR(f, h) ù SR2(f, b) = GE(f, h) * (1ù O1(f )) [16] L'estimation finale du ground roll est obtenue par application d'un filtre inverse. GE(f, h) = SR(f' h) ù SR2(f, h) [17] 1ù01(f) Il reste à soustraire cette estimation GE(f,h) du ground roll aux 5 données SR(f,h) pour obtenir une estimation PE(f,h) des ondes réfléchies. PE(f, h) = SR(f, h) ù GE(f, h) [18] La figure 9 représente de manière schématique des étapes d'un procédé d'exploration sismique comprenant un procédé d'acquisition 10 vibrosismique conforme à l'invention, composé des étapes 201-204, et un procédé de filtrage des ondes de surface selon un deuxième mode de réalisation (étapes 205-209). Selon une première étape 201, la source 1 est commandée pour générer un premier signal sismique SI (signal fort) présentant une première 15 amplitude. A cet effet, le vibrateur est réglé pour appliquer sur le sol des vibrations avec un premier niveau de force h (niveau haut), par exemple avec un drive d de 80%. Selon une deuxième étape 202, l'enregistreur 3 enregistre des premières données. 20 Les premières données enregistrées sont représentées sur la figure 3. Ces premières données comprennent un ensemble de signaux (ou traces) enregistrées en fonction du temps et résultant de la propagation du premier signal sismique, chaque signal ayant été généré par un récepteur 2. Les premières données comprennent des données correspondant à 25 des ondes de surface 4 et des données correspondant à des ondes réfléchies 6. Selon une troisième étape 203, la source 1 est commandée pour générer un deuxième signal S2 (signal faible) présentant une deuxième amplitude. A cet effet, le vibrateur est réglé pour appliquer sur le sol des 30 vibrations avec un deuxième niveau de force b (niveau bas), par exemple avec un drive d de 20%. Le deuxième signal S2 est identique au premier signal SI , excepté qu'il présente une amplitude inférieure à celle du signal S~ , comme cela est représenté sur la figure 2. Comme dans le cas du premier mode de réalisation, il est à noter que les termes � premier � et � deuxième � ne sont utilisés ici que pour identifier de façon pratique d'une part les signaux émis au niveau haut et d'autre part les signaux émis au niveau bas, ainsi que les données respectivement enregistrées à la suite de ces émissions. Les termes � premier � et � deuxième � n'indiquent en rien un ordre chronologique. L'ordre dans lequel ont lieu les deux émissions, et par suite les enregistrements de données, est indifférent pour atteindre les objectifs de l'invention. L'émission d'un signal haut peut précéder celle d'un signal faible, ou l'inverse. Selon une quatrième étape 204, l'enregistreur 3 enregistre des deuxièmes données. Selon une cinquième étape 205, un opérateur de passage O1(f) est déterminé à partir des première et deuxième données (équation [14]). La détermination de l'opérateur de passage repose sur le fait qu'il existe généralement des zones dans lesquelles les ondes réfléchies P et les ondes de surface G n'interfèrent pas ou sensiblement pas et peuvent être séparées . Le fait d'utiliser plusieurs sismogrammes associés à différents récepteurs permet de réduire l'effet d'un non-respect strict de la condition de non-interférence. Au besoin, le calcul de l'opérateur de passage peut faire l'objet 25 d'une itération, d'une manière bien connue de l'homme du métier. Selon une sixième étape 206, l'opérateur de passage O1(f) est appliqué aux deuxièmes données (équation [15]). Selon une septième étape 207, les deuxièmes données ainsi traitées ou transformées sont soustraites aux premières données. 30 Cette septième étape a pour effet d'éliminer les composantes P d'onde de pression réfléchie. La figure 5 représente de manière schématique la transformée de Fourier inverse des données soustraites. On constate sur cette figure que les ondes réfléchies ont sensiblement disparu. Selon une huitième étape 208, l'opérateur inverse 1 est 1ù01(f) appliqué aux données soustraites de manière à obtenir des données d'estimation GE des composantes d'onde de surface (ground roll). La figure 6 représente de manière schématique la transformée de Fourier inverse des données d'estimation GE des composantes d'onde de surface (ground roll). Selon une neuvième étape 209, les données d'estimation du ground 10 roll (figure 6) sont soustraites aux premières données (figure 3), de manière à éliminer les composantes d'onde de surface. Cette étape permet d'obtenir une estimation PE des données correspondant aux ondes réfléchies. Dans cette estimation, le ground roll a sensiblement disparu. 15 On a noté plus haut que dans les deux modes de réalisation décrits, l'ordre chronologique est indifférent : l'émission du signal de niveau haut peut précéder celle du signal de niveau faible, ou l'inverse. En ce qui concerne l'intervalle de temps séparant les deux émissions, plusieurs options sont possibles. Le plus simple du point de vue de la séparation des 20 deux enregistrements est de réaliser les deux émissions et enregistrements séquentiellement, sans recouvrement. On peut néanmoins aussi envisager de réduire l'intervalle entre les émissions, pour des raisons de productivité des opérations, et avoir des émissions et enregistrements qui se recouvrent. Dans ce dernier cas, il sera nécessaire de réduire la distorsion 25 due au bruit harmonique. Cela pourra être obtenu de façon appropriée par la méthode décrite dans le document FR 2845483. REVENDICATIONS1. Procédé d'acquisition de données vibrosismiques concernant une zone du sous-sol, comprenant les étapes suivantes : - activer une source vibrosismique pour qu'elle émette un premier signal vibrosismique présentant une première amplitude (101, 201), - activer la source pour qu'elle émette un deuxième signal vibrosismique présentant une deuxième amplitude, différente de la première amplitude (103, 203), - enregistrer des premières données correspondant au premier signal après propagation dans le milieu à explorer (102, 202), - enregistrer des deuxièmes données correspondant au deuxième signal après propagation dans le milieu à explorer (104, 204). 2. Procédé selon la revendication 1, dans lequel les premier et deuxième signaux (S(t,h),S(t,b)) émis par la source ne diffèrent que par leur amplitude. 3. Procédé selon la revendication 1, dans lequel les premier et 20 deuxième signaux (S(t, h), S(t, b)) présentent une forme et une durée différentes. 4. Procédé selon la revendication 1, dans lequel les premier et deuxième signaux sont émis de façon séquentielle. 5. Procédé selon la revendication 1, dans lequel les premier et deuxième signaux sont émis de façon à présenter un recouvrement. 6. Procédé d'exploration sismique d'une zone du sous-sol, 30 comprenant des étapes de : - (101, 201) activer une source vibrosismique pour qu'elle émette un premier signal vibrosismique présentant une première amplitude, 25 - (103, 203) activer la source pour qu'elle émette un deuxième signal vibrosismique présentant une deuxième amplitude, différente de la première amplitude, - (102, 202, 104, 204) enregistrer des premières données (SR(t, h) et des deuxièmes données (SR(t, b) correspondant respectivement au premier signal et au deuxième signal après leur propagation dans le milieu à explorer, - (105-109, 205-209) obtenir une estimation des ondes de surface contenues dans ces données à partir des premières et deuxièmes données 10 enregistrées, et soustraire cette estimation aux premières données. 7. Procédé selon la revendication 6, dans lequel les premières et deuxièmes données enregistrées contiennent des composantes d'ondes réfléchies (P(t,h), P(t,b)), résultant d'une réflexion des signaux dans le 15 milieu et des ondes de surface (G(t,h)) résultant d'une propagation des signaux en surface du milieu, le procédé comprenant des étapes de : - (106, 206) transformer les deuxièmes données pour égaliser la composante d'ondes réfléchies (P(t,b)) contenue dans les deuxièmes données avec la composante d'ondes réfléchies (P(t,h)) contenue dans les 20 premières données, - (107, 207) soustraire les deuxièmes données transformées (SR2(t,b)) aux premières données de manière à éliminer la composante d'ondes réfléchies afin de ne conserver que les composantes d'ondes de surface. 25 8. Procédé selon la revendication 7, dans lequel l'étape de transformation (106 comprend : - multiplier les deuxièmes données (SR(t,b)) par un coefficient (h/b) pour égaliser l'amplitude des composantes d'ondes réfléchies 30 contenues dans les premières et deuxièmes données, 5- décaler les deuxièmes données (SR(t, b)) d'un temps de décalage (tp -TG) pour ramener les composantes d'ondes réfléchies contenues dans les premières et deuxièmes données à une référence temporelle commune. 9. Procédé selon la revendication 8, comprenant une étape préalable de : - (105) mesurer le temps de décalage (tip ùTG) entre les deuxièmes et premières données. 10 15 20 25 10. Procédé selon la revendication 7, dans lequel l'étape de transformation comprend : - (206) appliquer aux deuxièmes données (SR(f, b)) un opérateur de passage (O1(f) ). 11. Procédé selon la revendication 10, comprenant une étape préalable de : - (205) déterminer un opérateur de passage (O1(f)) entre les premières et deuxièmes données. 12. Procédé selon la revendication 11, dans lequel l'opérateur de passage (O1(f)) est déterminé à partir de zones des premières et deuxièmes données, dans lesquelles les ondes réfléchies et les ondes de surface n'interfèrent sensiblement pas entre elles. 13. Procédé selon l'une des revendications 10 à 12, dans lequel l'opérateur de passage (O1(f)) est une fonction dépendant de plusieurs paramètres incluant la fréquence (f) et les nombres d'ondes (kx, ky) dans des directions horizontales. 30