Brevet FR2575298A1 — publié le 19860627

L'invent.i.on concezne 1une manière générale un procédé utile pour corriger les effets de la couche alté- rée dans le traitement de données sismiques. L'invention concerne plus particulièrement un procédé destiné à utiliser les arrivées premières réfractées, à partir de plusieurs jeux de données sismiques de déports constants pour calculer en conformité avec la surface, des corrections de temps unidi- rectionnelies de la couche altérée pour chaque récep- teur et chaque emplacement de tir dans une prospection sismique. L'exploration sismique est l'une des techniques les plus puissantes pour étudier la configuration des couches rocheuses au-dessusde la surface du sol. Le produit final classique d'un levé sismique est une carte, appelée "coupe de profondeur sismique", indiquant l'épaisseur et l'orientation des diverses couches s'étendant au-dessous de la partie de la surface du sol à partir de laquelle la prospection a été menée. En mettant en corrélation la coupe de profondeur sismique avec d'autres informations géologiques, telles que des données concernant les affleurements de diver- ses couches, des carottages et des diagraphies de puits, et des coupes sismiques précédentes, on peut élaborer une infor- mation étonnanment détaillée concernant, sur plusieurs kilo- mètres, la croûte terrestre extérieure. L'exploration sis- mique est surtout utilisée dans la recherche de structures souterraines favorables à la présence de gisements de pétrole et de gaz. Les levés de sismique réflexion, le type le plus courant de levé sismique,sont réalisés par le déclen- chement d'une onde de choc à la surface de la terre et l'écoute, en plusieurs lieux de la surface, des réflexions de cette perturbation sur les formations souterraines sous- jacentes. Ces réflexions proviennent de zones o le terrain présente une variation de l'impédance acoustique, en général l'interface entre des couches adjacentes. Les dispositifs utilisés pour l'écoute des réflexions sont appelés géophones. Le signal enregistré par chaque géophone représente, en fonction du temps, l'amplitude des réflexions détectées par ce géophone. Pour une bonne approximation, les réflexions détectées par chaque géophone proviennent d'un point de chaque surface réfléchissante situé sur une verticale passant par le point se trouvant à mi-distance entre la source et le géophone. Ainsi, pour chaque perturbation sismique ("tir"), chaque géophone enregistre un signal ("trace"} qui représente des caractéristiques des formations situées verticalement au-dessous d'un point connu de la surface de la terre. Dans la réalisation d'une prospection ou d'un levé sismique, un grand nombre de géophones, habituellement entre 48 et 1024, sont positionnés le long de la ligne du levé. En conséquence, pour chaque tir, on obtient de nombreu- ses traces. Chacune des traces résultant d'un seul tir repré- sente les réflexions en provenance des interfaces suivant une ligne verticale unique traversant les formations souterraines. A l'instant o chaque trace est enregistrée, elle est désignée d'une façon qui lui est propre d'après les positions de la source et du détecteur. De cette manière, chaque trace est identifiée de façon particulière par rapport à toutes les autres traces. Cette information est utilisée par la suite pour la correction et l'affichage des traces. Après chaque tir, la source est déplacée le long de la ligne du levé et un second tir est réalisé, donnant un nouveau jeu de traces. En général, les géophones sont espacés régulièrement et le mouvement de la source par rapport aux géophones est établi avec soin de façon que, à l'exception de l'un des deux géophones extrêmes, chaque point central source-détecteur à partir du premier tir corresponde à un point central source-détecteur à partir du second tir. Dans le type le plus simple de levé sismique, ceci est obtenu par déplacement de la source et des détecteurs d'un seul bloc, sur une distance égale à l'écartement des géo- phones, après chaque tir. En continuant le levé de cette manière, les nombreuses traces correspondent chacune à l'un d'un nombre plus petit de points centraux source-détecteur. Par conséquent, chaque point central est représenté par plusieurs traces. Des traces ayant des points centraux communs représentent des réflexions provenant des interfaces situés sur une ligne commune descendant verticalement à - partir du point central. La conduite du levé sismique de manière que chaque trace provenant d'un premier tir corres- ponde à une trace provenant du tir suivant est appelée pros- pection en couverture multiple. Les traces obtenues dans l'exécution du levé ou de la prospection doivent être corrigées avant l'affichage final et l'analyse afin de compenser divers facteurs qui empêchent une comparaison directe des traces d'origine. L'une des plus ennuyeuses de ces corrections consiste à compenser les effets, sur les traces, de la couche la plus haute du terrain, généralement d'une épaisseur de 10 à 100 mètres, appelée "couche de faible vitesse" ou "couche altérée". La vitesse des ondes sismiques de compression (ondes p) à travers la couche altérée est généralement de l'ordre de 500-1000 mètres/seconde, alors que les vitesses des ondes p dans les couches situées au-dessous de la couche altérée sont géné- ralement supérieures à 1500 mètres/seconde. Etant donné que la couche altérée diffère souvent notablement en épaisseur sur des distances horizontales relativement courtes, le temps de transit d'une onde sismique à travers la couche altérée peut varier notablement sur la ligne de levé sismique. Si elle n'est pas corrigée, cette variation peut altérer notablement la configuration et la profondeur observées des couches sous-jacentes. Par exemple, en supposant une strati- fication horizontale des couches sous-jacentes à la couche altérée, une zone mince d'une couche altérée ayant par ailleurs une épaisseur constante peut faire apparaître comme convexe une partie des couches horizontales. Etant donné que des variations même faibles de l'orientation calculée des couches rocheuses peuvent avoir un effet important sur des décisions concernant la probabilité que l'on a de trouver du pétrole et du gaz dans certains emplacements souterrains, il est important que des aberrations' dues à la couche altérée soient déterminées avec la plus grande précision possible. Un procédé antérieur de correction des effets de la couche altérée est décrit dans le brevet des Etats- Unis d'Amérique N 2 276 306. Dans cette technique, on utilise de la dynamite pour déclencher la perturbation sismique. Chaque charge de dynamite est placée au fond d'un trou foré à travers la couche altérée. Un géophone placé à proximité de l'embouchure du trou de la source enregistre le temps de transit vertical à travers la couche altérée, en même temps qu'un groupe de géophones espacés le long de la ligne de prospection enregistrent les réflexions à partir des couches sousjacentes. A partir du géophone placé à l'embouchure de chaque trou de source, on connaît le temps de transit vertical en fonction de la position le long de la ligne du levé sismique. Ceci permet aux diverses traces enregistrées pendant la prospection d'être soumises à une correction simple et précise portant sur les effets de la couche altérée. Dans l'exploration sismique moderne, ce procédé est rarement utilisable, car les sources placées en surface ont largement remplacé l'utilisation des sources souterraines en raison du co t relativement élevé et des difficultés d'environnement associés à ces dernières. Une autre technique bien connue pour établir l'effet de la couche altérée sur des données sismiques réfléchies est le procédé de réfraction à ordonnée à l'ori- gine'. Dans ce procédé, pour chaque tir, les traces résul- tantes sont examinées afin de déterminer le temps nécessaire à l'onde sismique pour parcourir un trajet allant de la sources passant t -.e a couc e. y ' 1y l i' neC=sé;ece situe aei bos d5e te-g --'- at.,g._qe "? _11sSec z.e touche, et susin!aque3l!e elle est rfractée- jusqu - ce qu'ell- arrive ae- récepteurs. Ce trajet de réfraction est illustré sur la figure 1 des dessins annexés et décrits ci-après. Etant donné que ce trajet sismique est en général le plus rapide pour aller de la source au récepteur, il s'agit du premier signal reçu par ce dernier. Pour chaque position de source, le temps des arrivées premières réfractées est tracé en fonction du déport sourcerécepteur. Un calage par moindres carrés est appliqué aux données pour que l'on obtienne une ligne droite dont la pente représente la vitesse de réfraction V2, et ayant un temps d'ordonnée à l'origine spécifique tr, corres- pondant à une distance de déport nulle. Par l'application de la loi de Snell, on démontre que le temps de transit vertical unidirectionnel à travers la couche altérée, tLVL, est donné par l'équation: t_ tLVL = r2 (V1/2 o V1 = vitesse de l'onde P sismique dans la couche altérée. Cependant, étant donné que le graphique sur lequel le procédé de réfraction à ordonnée à l'origine est basé ne fait pas la distinction entre le temps de transit dans la couche altérée à la source et au récepteur, le temps de transit unidirectionnel vertical calculé par ce procédé est une moyenne des deux valeurs. Dans les méthodes de traitement sismique modernes, l'incertitude introduite par ce calcul de moyennes représente, dans de nombreux cas, une fraction notable de l'erreur totale d'une trace corrigée, obtenue par sommation. Il est souhaitable de prévoir un procédé pour établir des corrections statiques pour la couche altérée, procédé dans lequel la correction effectuée à l'emplacement de chaque source est indépendante de la couche altérée à l'emplacement d'un récepteur quelconque, et la correction effectuée à l'emplacement de chaque récepteur est indépen- dante de la couche altérée à l'emplacement d'une source quelconque. Il est en outre souhaitable que ce procédé de correction statique portant sur la couche altérée n'exige aucun processus ni équipement spécial sur place, dépendant uniquement des données normalement acquises dans le cours d'une prospection sismique par ondes de compression, en couverture multiple. Il est en outre souhaitable que ce procédé donne une correction absolue plutôt que relative portant sur la couche altérée. L'invention concerne donc un procédé qui con- vient à la correction sismique des effets de la couche altérée dans des levés ou prospections sismiques par réflexion et réfraction d'ondes p et d'ondes s. Selon un aspect de l'invention, on établit trois groupements possé- dant chacun des traces d'un déport commun choisi. Les traces de chaque groupement sont affichées en fonction d'un emplacement de couverture multiple ou position miroir. A partir de chaque groupement, on détermine la première rupture pour les traces, la première rupture représentant le temps de transit sur le trajet de parcours de l'onde p, de la source au récepteur, le long de l'interface de réfraction au bas de la couche altérée. Pour chaque trace présentant l'un des trois déports communs choisis, on développe une équation mettant en relation le temps de transit connu avec la somme de trois facteurs: 1) la dis- tance source-récepteur divisée par la vitesse de l'onde p à l'interface de réfraction; 2) une participation de la source au temps de parcours; et 3) une participation du récepteur au temps de parcours. Chacune de ces équations possède deux inconnues (facteurs 2 et 3); cependant, au moins une inconnue de chaque équation est partagée par deux autres équations, correspondant aux autres distances de déport. Par conséquent, les inconnues sont surdéterminées et les équations peuvent être résolues par des méthodes classiques de résolution simultanée par les moindres carrés, appliquées par ordinateur. Dans une autre forme de réali- sation de l'invention, on obtient une autre amélioration dans la détermination de la correction statique concernant la couche altérée en utilisant au moins un groupement de déport supplémentaire et en procédant à une résolution portant sur les variations latérales de la vitesse de réfraction. L'invention sera décrite plus en détail en regard du dessin annexé à titre d'exemple nullement limitatif et sur lequel: - la figure 1 est un schéma illustrant la trans- mission de signaux sismiques à partir d'une paire de sources vers un certain nombre d'emplacements de récepteur, par réfraction sur la limite inférieure de la couche altérée et réflexion sur une interface réfléchissante située au-dessous de la couche altérée; et - la figure 2 est un schéma illustrant le trajet des rayons d'un signal sismique le long de la limite inférieure de la couche altérée. L'invention concerne un procédé de traitement de donnéessismiques pour établir une correction statique, conforme à la surface, portant sur les effets de la couche altérée. Dans la forme préférée de réalisation, ce procédé est appliqué à des données sismiques multiples, obtenues en couverture multiple à partir d'un levé ou d'une prospec- tion par réflexion en progression, réalisé suivant un trajet de levé sensiblement linéaire. Cependant, il est évident à l'homme de l'art que la présente invention peut être appliquée au traitement de données provenant de nombreux types différents de levés sismiques. Le fait que la descrip- tion suivante porte de façon particulière sur des levés par réflexion en progression est illustratif plutôt que limi- tatif. La figure 1 est une coupe transversale de la surface de la terre illustrant la transmission d'un signal sismique dans un levé ou une prospection classique de sismique réflexion. Une perturbation sismique est introduite dans le sol à partir de chacun de plusieurs emplacements de source sismique S1; S2...Sn disposés en une rangée sensiblement linéaire le long de la surface 10 de la terre. Immédiatement au-dessous de la surface 10 se trouve une couche altérée (couche de faible vitesse) 12, et, au-dessous de la couche altérée 12, une couche non altérée (couche à grande vitesse) 14. La limite entre les couches altérée, et non altérée 12, 14 définit une interface 16 de réfraction. La figure 1 montre qu'à un angle critique d'incidence, une onde sismique descendant à travers la couche altérée 12 est réfractée le long de l'interface de réfraction 16 plutôt que transmise dans la couche non altérée 14. Ceci est illustré sur la figure 1 par les deux trajectoires 24 et 26 issues de la source S2. Une première trajectoire 26 arrive à l'inter- face 16 de réfraction au-dessous de l'angle critique d'inci- dence et est donc transmise dans la couche non altérée 14 afin d'être réfléchie par les diverses interfaces plus profondes, y compris l'interface 20. Une seconde trajectoire 24 arrive à l'interface 16 de réfraction sous l'angle critique d'incidence et est donc réfractée le long de l'interface 16 de réfraction au lieu d'être transmise dans la couche non altérée 14. L'angle critique, mesuré par rapport à une perpendiculaire à l'interface 16 de réfraction au point d'incidence, est donné par la loi de Snell: ecritique = cos-1 équation 1 o V1 = vitesse de l'onde sismique p dans la couche altérée V2 = vitesse de l'onde sismique p dans la couche non altérée LYJ... c t- par la trajectoire, quelque peu plus longue, suivant l'inter- face de réfraction 16. Ceci est dû au fait que la vitesse dans la couche altérée est généralement de 500-1000 mètres/ seconde, tandis que celle dans la couche non altérée est en général supérieure à 1500 mètres/seconde. Par conséquent, hormis dans les cas d'espacements source-détecteur relative- ment courts, l'allongement de la trajectoire représenté par le transit du signal vers et à partir de l'interface de réfraction 16 est plus que compensé par la plus grande rapidité de transmission de l'onde p de l'interface de réfraction 16. En conséquence, la première partie du signal sismique à détecter (la "première rupture") par chaque récepteur suit généralement une trajectoire passant le long de l'interface 16 de réfraction, comme montré sur la figure 2. Dans une forme préférée de réalisation de l'invention, la correction statique portant sur la couche altérée est obtenue à partir d'une résolution d'un groupe surdéterminé d'équations simultanées basées sur le temps de transit de l'onde p le long de l'interface de réfraction 16. Ces équations sont de la forme: T.=S. + Rj +X. Tij = Si + R j équation 2 V2 o T.. est le temps de transit de l'onde sismique p sur un parcours passant par l'interface de réfraction, à partir d'une source placée en i jusqu'à un récepteur placé en j Si est la participation de la source au temps de parcours Rj est la participation du récoepteur au temps de parcours J Xij est la distance entre les points i et j V2 est la vitesse de réfraction, qui est également la vitesse de l'onde sismique p dans les couches situées immédiatement au-dessous de l'interface de réfraction. Tij et Xij sont évidemment connus pour chaque paire de récepteurs. Dans la présente forme de réalisation, on suppose que V2 est constante en tous points le long de la ligne du levé. V2 peut être mesurée au moyen d'un levé par réfraction ou selon d'autres manières bien connues de l'homme de l'art. Dans la mise en oeuvre préférée du procédé de l'invention, on choisit trois distances de déport source- récepteur à partir du levé sismique en couverture multiple. Ainsi qu'il apparaîtra ci-après, il est commode, mais non nécessaire à la mise en oeuvre du présent procédé, de se baser sur un déport commun dans le choix des traces à utiliser dans le présent procédé. Pour chacune des distances de déport choisies, on établit et on affiche un groupement de toutes les traces de récepteurs en fonction de l'empla- cement du point central source-détecteur. A partir de ce groupement, on peut déterminer visuellement les arrivées premières réfractées, représentant Tij pour chaque trace. Ensuite, on établit un groupe d'équations simultanées de la forme de l'équation 2 pour chaque trace comprise dans l'un quelconque des trois groupes de traces représentant les trois distances de déport choisies. Chacune des équations ne contient que deux inconnues, Si et Rj. Le groupe d'équations simultanées contient des termes communs. Rj est identique pour toutes les traces représentant des paires source-récepteur ayant un emplacement de récepteur commun. Ceci est dû au fait que, comme montré sur la figure 1, seule la partie initiale de la trajectoire du signal sismique diffère pour des signaux sismiques allant de sources différentes vers un récepteur commun. De façon similaire, Si est identique pour toutes les traces représentant des paires source-récepteur ayant un emplacement de source commun. En outre, lorsque le levé sismique est effectué de façon que les emplacements des sources coïncident avec les emplacements des récepteurs, pour obtenir une bonne approximation, Ri sera égal à Sj i J lorsque R. et Sj représentent un emplacement de surface commun. Le nombre d'inconnues dans le groupe d'équations simultanées développé à partir de toutes les traces repré- sentant les trois distances de déport choisies est inférieur au nombre d'équations indépendantes. En conséquence, le système d'équationsest surdéterminé et toutes les valeurs de Si et Rj peuvent être établies par une solution simultanée par moindres carrés. La solution simultanée par moindres carres pour toutes les valeurs de Si et Rj dans les équations peut être calculée par des méthodes mises en oeuvre par ordinateur, bien connues de l'homme de l'art. Une meilleure précision dans la détermination des diverses valeurs de Si et Rj peut être obtenue par un accroissement du nombre de valeurs de déport représentées dans le système d'équations simultané. Dan.s le cas d'une vitesse de réfraction cons- tante, il est apparu que pour des données sismiques de qualité moyenne, le meilleur équilibre entre la précision du résultat et l'efficacité du temps d'utilisation de l'ordinateur est obtenu avec l'utilisation de six distances de déport. Après l'obtention de la correction du temps de parcours pour chaque emplacement particulier de surface représenté dans le levé, on peut obtenir le temps de transit vertical unidirectionnel à travers la couchealtérée, à l'emplacement de chaque source ou de chaque récepteur, par une équation de la forme: s. iL_ 1/2 équation 3 (v2j La valeur, pour chaque emplacement de surface, est utilisée pour la correction de temps statique absolu pour des coupes-temps. Il est simple d'effectuer une conver- sion donnant la correction de profondeur statique pour des coupesprofondeur. Un avantage important de la mise en oeuvre de la présente invention est que les temps de transit calculés sont en conformité avec la surface. Autrement dit, le temps de transit dans la couche altérée, calculé pour l'emplacement de chaque récepteur, est indépendant du temps de transit dans la couche altérée pour les emplacements de toutes les sources servant d'origine à tous les signaux sismiques reçus à cet emplacement de récepteur. En outre, le présent procédé donne une correction statique absolue plutôt que relative concernant la couche altérée. Ceci facilite l'établissement de liens de temps entre des lignes sismiques sécantes. Un autre avantage est qu'il n'est pas nécessaire de procéder à des opérations spéciales dans l'exécution du levé sismique. Ce procédé de réalisation de la correction statique portant sur la couche altérée n'exige aucune autre donnée que celle obtenue dans un levé sismique normal. Un perfectionnement de la présente invention permet de prendre en considération les variations latérales de la vitesse de réfraction sur la longueur du levé. Il est important de corriger ces variations lorsqu'elles sont présentes. Si la vitesse de réfraction latérale est supposée à tort comme étant constante, les variations de la vitesse latérale seront traduites incorrectement en variations d'épaisseur dans la couche altérée. Pour tenir compte de ces variations, on modifie l'équation 2 afin de prendre en 1 3 compte une comp sante tl_ de vitrsable deg 4= À=. + +_ % X.. Tij = isl Vk équation 4 o Vk est une inconnue représentant la vitesse moyenne de réfraction à un emplacement k0 pour toutes les paires sourcerécepteur ayant un emplacement k en tant que point central. Ainsi, pour un groupe unique de traces de déport commun, chaque trace correspond à un emplacement unique k et, si Vk n'est pas constante, à une vitesse unique de réfraction Vk. En présence de groupes multiples de traces, chaque groupe ayant un déport commun, chaque emplacement k et chaque valeur Vk pour un groupe de traces à déport commun correspondent à un emplacement k et à une valeur Vk dans chaque autre groupe de traces à déport commun. Pour sur- déterminer le système d'équations de la forme de l'équation 4, il faut quatre groupes de traces à déport commun. Cepen- dant, dans cette situation à vitesse de réfraction variable, il est apparu que le meilleur équilibre entre la précision du résultat et l'efficacité du temps d'utilisation de l'ordinateur est obtenu avec l'utilisation de dix distances de déport. En mettant en oeuvre ce procédé, il est avantageux que les déports choisis soient symétriques de part et d'autre d'un profil à cheval. Ceci minimise les effets du pendage sur la vitesse de réfraction Vk. On a, après la solution par moindre carré: si tLVL i= 1/2 i 1/ Il va de soi que de nombreuses modifications peuvent être apportées au procédé décrit et représenté sans sortir du cadre de l'invention. REVENDICATIONS 1. Procédé pour corriger des effets de la couche altérée (12) des données sismiques multiples, carac- térisé en ce qu'il consiste: - (a) à choisir au moins trois groupes de traces sismiques à partir desdites données sismiques, les traces de chaque groupe ayant un déport commun; (b) à déterminer pour chaque trace desdits groupes de traces sismiques le temps de transit de la source (S1, S2...) au récepteur (R1... Rn) sur la trajectoire du signal sismique correspondant à une réfraction le long de la limite inférieure (16) de la couche altérée; (c) à établir pour chaque trace d'un groupe de trace sismique une équation de la forme: T. S +R Tij i Rj T1 o Tij = temps de transit source-récepteur le long de l'interface de réfraction à la limite inférieure de la couche altérée Si = part; dans le temps,de la source à l'emplace- ment i R. = part, dans le temps,du récepteur à l'emplace- ment j T1 = temps de transit d'une onde sismique parcourant une distance préalablement choisie le long de la limite inférieure de la couche altérée; et (d) à résoudre les équations établies à l'étape c. 2. Procédé selon la revendication 1, caractérisé en ce que T X.. V2 o Xij = distance de la source au récepteur V2 = vitesse de réfraction uniforme. 3. Procédé selon la revendication 1, caractérisé en ce que l'étape d consiste à effectuer une résolution simultanée par moindres carrés pour les valeurs de S. et R. à chacun des emplacements i et j d'une source Et d'un J récepteur. 4. Procédé selon la revendication 1, caractérisé en ce qu'il consiste en outre (e) à convertir la participation au temps dans la couche altérée, due à l'emplacement de chaque source et de chaque récepteur, en une épaisseur de la couche altérée a cet emplacement. 5. Procédé selon la revendication 1, caractérisé en ce qu'au moins quatre groupes de traces sismiques sont utilisés, chaque groupe possédant une valeur de déport commune différente, et dans lequel T1 = Xij Vk o Vk est une inconnue représentant la vitesse moyenne de réfraction à l'emplacement k pour toutes les paires source- récepteur englobées dans lesdits groupes choisis de traces ayant un emplacement k comme point central. 6. Procédé selon la revendication 5, caractérisé en ce que les groupes choisis de traces sismiques sont disposés symétriquement de part et d'autre d'un profil à cheval. 7. Procédé pour établir la correction statique de couche altérée (12) sur des données sismiques multiples, caractérisé en ce qu'il consiste: a) à choisir parmi les données résultant d'un levé sismique à un groupe de traces représentant.une série d'emplacements de couverture multiple le long de la ligne du levé sismique, le groupe choisi de traces représentant un nombre total d'emplacements de sources individuelles (S1, S2...) et de récepteurs (R1.. Rn) qui est inférieur au nombre total de traces choisies b) à déterminer pour chaque trace faisant partie du groupe de traces choisi le temps de transit source- récepteur pour la trajectoire du signal d'onde p corres- pondant à une réfraction le long de la limite inférieure (16) de la couche altérée; c) à établir pour chaque trace du groupe de traces choisies une équation de la forme T.. = S. + R. + X.. V2 o Tij = temps de transit source-récepteur le long de l'interface de réfraction à la limite inférieure de la couche altérée Si = part, dans le temps, de la source à l'emplace- ment i R. = part, dans le temps,du récepteur à l'emplace- J ment j Xi. = distance de la source au récepteur V2 = vitesse uniforme de réfraction le long de la surface inférieure de la couche altérée; et d) à résoudre le système d'équations simultanées établi dans l'étape c pour chacun des emplacements parti- culiers i et j du récepteur et de la source. 8. Procédé selon la revendication 7, caractérisé en ce que le groupe choisi de traces représente toutes les traces provenant du levé sismique, ayant un déport égal à l'une d'au moins trois valeurs de déport choisies. 9. Procédé selon la revendication 7, caractérisé en ce que V2 = Vk, Vk étant une inconnue représentant la vitesse moyenne de réfraction à l'emplacement k pour toutes les paires source-récepteur comprises dans lesdits groupes choisis de traces, ayant un emplacement k comme point central. 10. Procédé selon la revendication 7, caractérisé en ce qu'il consiste en outre e) à convertir la participa- tion au temps dans la couche altérée de chaque emplacement de source et de récepteur en un temps de transit de l'onde p verticale dans la couche altérée pour chaque emplacement particulier de source et de récepteur représenté dans ledit groupe úl t.Z:ECfS t-ii 11. Procedé pour calculer le temps de parcours d'une onde de compression à travers la couche altérée (14) d'une formation souterraine en chacun de plusieurs points espacés le long de la surface de la formation, caractérisé en ce qu'il consiste: à déclencher plusieurs perturbations sismiques en des emplacements de sources (S1, S2.o o) espacés le long de la surface de la formation; à effectuer un enregistrement de chaque pertur- bation sismique en plusieurs emplacements de récepteurs (R1... Rn)espacés le long de la surface de la formation; à déterminer pour chaque trace correspondant à plusieurs paires choisies source-récepteur le temps de transit de la source au récepteur pour un signal d'onde p réfracté le long de l'interface inférieure (16) de la couche altérée; et à résoudre l'équation Tij = S + R + pour Xij our V2 lesdites paires choisies sourcerécepteur o Tij = temps de transit source-récepteur le long de l'interface de réfraction à la limite inférieure de la couche altérée Si = part, dans le temps,de la source à l'emplacement i Rj = part, dans le temps,du récepteur à l'emplacement j Xij = distance de la source au récepteur V2 = vitesse de réfraction le long de l'inter- face à la limite inférieure de la couche altérée o le nombre total de paires choisies source-récepteur, chacune étant la base d'une équation unique de la forme donnée ci-dessus, est plus grand que le nombre total d'emplacements indépendants de source et de récepteur représentés par lesdites paires source-récepteur. 12. Procédé selon la revendication 11, caractérisé en ce que les paires choisies source-récepteur correspondent chacune à l'une de plusieurs distances choisies de déport. 13. Procédé selon la revendication 11, caracté- risé en ce qu'au moins certains des emplacements de source et des emplacement de récepteur représentent une position commune le long de la surface de la formation. 14. Procédé selon la revendication 11, caracté- risé en ce que la vitesse de réfraction V2 est supposée constante. 15. Procédé selon la revendication 14, caracté- risé en ce que la vitesse de réfraction V2 est égale à Vk, qui est une inconnue supplémentaire représentant la vitesse moyenne de réfraction à l'emplacement k pour toutes les paires source-récepteur ayant un emplacement k comme point central, et dans lequel le nombre total de paires choisies source-récepteur, chacune étant la base d'une équation unique de la forme donnée dans la revendication 11, est plus grand que la somme du nombre total d'emplacements de source et de récepteur indépendants et du nombre total de valeurs indépendantes de Vk.