Brevet FR2963833A1 — publié le 20120217

Procédé de détermination de la correction en temps pour un détecteur disposé au fond de la mer La présente invention concerne les techniques d'exploration du sous-sol, en particulier un procédé de détermination pour un détecteur disposé sous la mer notamment sur la surface du fond marin, du temps de propagation verticale to et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission. Il est connu, notamment dans l'exploration pétrolière, de produire des images sismiques à partir de séries de mesures géophysiques effectuées depuis la surface du sous-sol. Dans la technique de la sismique, ces mesures impliquent l'émission dans le sous-sol d'une onde et la mesure d'un signal comportant diverses réflexions de l'onde sur les structures géologiques rencontrées. Ces structures sont typiquement des surfaces séparant des strates géologiques distinctes ou des failles. Les images sismiques sont des représentations du sous-sol à deux ou trois dimensions. On appelle la dimension verticale une de ces dimensions correspondant soit aux temps de propagation des ondes sismiques, soit aux profondeurs. Elles sont obtenues par des techniques connues sous l'appellation de "migration" qui utilisent un modèle de vitesse estimée fournissant une cartographie de la vitesse de propagation des ondes sismiques dans les roches constituant la zone explorée. Ce modèle de vitesse est utilisé pour estimer les positions des réflecteurs du sous- sol à partir des enregistrements sismiques. Bien entendu, les images sismiques ainsi produites, ainsi que les modèles 2 de vitesse sous-jacents, présentent certaines distorsions car ce ne sont que des estimations dérivées d'un nombre nécessairement limité de mesures. Dans le cas d'exploration de sous-sol sous marin, des détecteurs peuvent être placés au fond de la mer sur le sous-sol à explorer. Des ondes sismiques sont émises à partir de points situés près de la surface de la mer. Ces ondes se propagent dans l'eau et pénètrent le sous-sol. Les détecteurs placés au fond de l'eau sur la surface du sous- sol vont détecter l'arrivée de l'onde sismique directe ainsi que les ondes réfléchies par le sous-sol. Afin de suivre l'évolution d'un réservoir pétrolier d'un sous-sol, il est possible de réaliser une première image sismique du sous-sol à un temps donné puis de réaliser une deuxième image sismique du même sous-sol après un certain temps. En particulier, pour suivre les changements en teneur en hydrocarbure d'un réservoir exploité en production, il peut être intéressant de suivre l'évolution de l'image sismique du sous-sol dans le temps. Afin de pouvoir comparer deux images sismiques d'un même sous-sol, il est important de savoir positionner le plus précisément possible chaque détecteur placé sur la surface dudit sous-sol. La demande ER 10 52600 décrit un procédé permettant de déterminer de façon précise la position des détecteurs placés sur la surface du sous-sol. Une particularité des détecteurs utilisés est d'être des détecteurs autonomes. Ces détecteurs sont munis d'une horloge interne qui doit être synchronisée avant déploiement et re-synchronisée après récupération. La durée de collecte des données pouvant atteindre plusieurs mois, il peut être observé des dérives temporelles. -3 La durée d'acquisition conduit à observer des dérives résiduelles des horloges internes des détecteurs, des erreurs de synchronisation entre détecteurs, en particulier des variations du temps origine, des décalages liés aux possibles variations de la vitesse de propagation des ondes sismiques dans la tranche d'eau. Ceci se superpose aux effets de marée, aux possibles instabilités du signal de la source sismique, à un positionnement imparfait des détecteurs et points sources, à l'imprécision de la vitesse de propagation des ondes dans la tranche d'eau et à l'incertitude sur les données de bathymétrie. La plupart des méthodes utilisées pour compenser ces effets sont sensibles à la phase du signal source. En général, ces méthodes ne peuvent être réalisées totalement en une seule passe de traitement, et elles ne fournissent pas de mesure de la qualité de la compensation autre que le contrôle d'une bonne mise en phase des évènements sismiques. Ainsi, il existe un besoin pour une méthode permettant de déterminer de manière précise l'ensemble des éléments permettant les corrections de temps sur des détecteurs placés sur la surface du fond marin. L'invention propose ainsi un procédé de détermination pour un détecteur disposé sur la surface du fond marin du temps de propagation verticale (to) et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission, le procédé comprenant les étapes suivantes pour au moins un des N points d'émission: - émission d'au moins une onde à partir dudit point d'émission, enregistrement de l'onde reçue par le détecteur, détermination du temps de propagation verticale to au moyen de la relation suivante : z Ao -tnnd -tdir -8tz o , avec tdir le temps de propagation de l'onde directe entre le point d'émission et le détecteur et tmn1 le temps de propagation de la 5 première onde multiple entre le point d'émission et le détecteur, détermination de la vitesse V de propagation de l'onde au moyen de la relation suivante : 2 2 (2~X~ A =9ttr -tmtf =(2 V avec X la distance entre le 10 point d'émission et le point à la surface de l'eau à la verticale du détecteur. Selon un mode de réalisation de l'invention, les quantités Al et A2 peuvent être mesurées sur les images 15 sismiques obtenues lors de l'enregistrement de l'onde reçue par le détecteur. Avantageusement, le procédé selon l'invention permet de déterminer le temps de propagation vertical et la vitesse de propagation de l'onde dans la tranche d'eau 20 situé au dessus du détecteur. Le temps de propagation vertical et la vitesse de propagation permettent de déterminer les corrections de temps pour un détecteur disposé sur la surface du fond marin. Un procédé selon l'invention peut en outre 25 comporter une ou plusieurs des caractéristiques optionnelles ci-dessous, considérées individuellement ou selon toutes les combinaisons possibles : - le temps de propagation vertical to est déterminé par les étapes suivantes . 30 ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=t2 auto-corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné de l'onde reçue par le détecteur détermination de la quantité to à partir de la valeur de T correspondant au pic principal du signal auto-corrélé, - la vitesse de propagation est déterminée par les étapes suivantes : ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=9t2, corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=9t2 avec le signal s'(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=t2, détermination de la quantité O1 à partir de la valeur de T correspondant au pic principal du signal corrélé, - la quantité A est déterminée au moyen de la relation suivante : Al = 8td;,2 -Do tain étant déterminé au moyen du temps d'arrivée directe Tdir pointé sur l'enregistrement de l'onde reçue par le détecteur, ledit temps d'arrivée directe pointé Tdir corrigé du décalage temporel dTo=To-to avec To le temps de propagation vertical estimé et to le temps de propagation vertical déterminé au moyen de dos - à partir des oo déterminés pour chacun des N points d'émission, on modélise les ,Ao obtenus pour chacun des points d'émission au moyen de l'équation suivante ~o~ = D+Bxx+Byy avec x et y les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur et avec D, BX et By des paramètres de modélisation dont les valeurs sont déterminées de sorte que Sod s'ajuste le mieux aux Do obtenus, par exemple au sens des moindres carrés, le temps de propagation vertical (ta) est déterminé D au moyen de la relation suivante: t° J8D-V2(BS +B2) ' à partir des L1 déterminés pour chacun des N points d'émission, on modélise les Al obtenus pour chacun des points d'émission au moyen de l'équation suivante Am0d=A(x2+y2)+Fxx+Fyy+E avec x et y les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur, et avec A, FX, Fy et E des paramètres de modélisation dont les valeurs sont déterminées de sorte que 0'I 0d s'ajuste le mieux aux Al obtenus, par exemple au sens des moindres carrés, la vitesse de propagation (V) est déterminée au moyen 2N/i de la relation suivante: V = - l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle rl entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : =tan-I , y= tan-' (VtoVBz +B2 D - l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle p entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : (F Y Fr1 q=tan' , y = sin JF2 + F2 AVto 1 , et - après l'enregistrement de l'onde reçue par le détecteur, les données en champs montant et descendant sont séparées et seules les données en champs descendant sont utilisées pour déterminer le temps de propagation vertical to et la vitesse de propagation V. L'invention se rapporte également à un procédé de correction de l'origine des temps des signaux enregistrés par un détecteur disposé sur la surface d'un fond marin, lesdits signaux correspondant à des ondes émises à partir d'au moins un point d'émission parmi N points d'émission, dans lequel les signaux reçus par le détecteur sont corrigés au moyen d'un procédé de type move-out et en tenant compte de l'origine des temps avec le temps to et la vitesse de propagation V déterminés au moyen d'un procédé selon l'invention. Avantageusement, le procédé selon l'invention n'est pas sensible aux effets dits de stretch ou d'étirement qui peuvent être créés par les méthodes usuellement utilisées. - 8 Enfin, il peut être mis en oeuvre tôt dans la séquence de traitement, tandis que, dans des méthodes plus classiques, les corrections propres à chaque détecteur sont réalisées après migration, ce qui oblige à suivre toujours le même traitement, ne permettant pas de mettre en oeuvre des techniques plus élaborées, par exemple l'Offset Vector Tiling. L'invention se rapporte également à un produit de programme d'ordinateur, comprenant une série d'instructions qui lorsque chargée dans un ordinateur entraîne l'exécution par ledit ordinateur des étapes du procédé selon l'invention. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : - la figure 1 est une représentation schématique d'une disposition d'un détecteur et d'un point d'émission selon un premier mode de réalisation, - les figures 2a et 2b sont des représentations schématiques d'une disposition d'un détecteur et d'un point d'émission selon un deuxième mode de réalisation, - la figure 3a représente le décalage temporel déterminé selon l'invention pour un ensemble de détecteurs placés sur la surface d'un sous-sol marin, - les figures 3b et 3c représentent les vitesses non- corrigées et corrigées du décalage temporel obtenues par un procédé selon l'invention, et - la figure 4 représente les enregistrements des ondes émises par une source et reçues par différents détecteurs placés sur la surface du sous-sol marin sans et avec correction des origines des temps. Pour des raisons de clarté, les différents éléments représentés sur les figures ne sont pas nécessairement à l'échelle. Le procédé selon l'invention exploite les temps d'arrivée de l'onde directe et du multiple du fond de l'eau pour estimer le temps de propagation verticale à la verticale du détecteur ainsi que la vitesse de propagation dans la tranche d'eau située entre le détecteur et le point d'émission. Au sens de l'invention, on entend par « temps de propagation verticale » le temps que mettrait une onde directe émise à partir d'un point d'émission placé directement à la verticale du détecteur pour parcourir la distance entre ledit point d'émission et ledit détecteur. Au sens de l'invention, on entend par « méthode de type move-out » des méthodes permettant de corriger le décalage de temps dû à la position dans un plan horizontal du point d'émission. Des méthodes de correction de type linear move-out ou normal move-out sont bien connues de l'homme du métier. La figure 1 représente une projection selon un plan vertical défini par le détecteur 10 disposé sur la surface 11 d'un fond marin et le point d'émission 12 disposé à la surface de la mer. Comme illustré sur la figure 1, lorsqu'une onde, en particulier une onde sismique, est émise à partir du point d'émission 12, ladite onde se propage jusqu'au détecteur 10. L'onde peut se propager directement du point 30 d'émission 12 au détecteur, on parle alors d'onde directe 14. -10- Comme il apparaît sur la figure 1, le temps de propagation de l'onde directe entre le point d'émission et le détecteur vérifie l'équation suivante : tir = - + - , avec Z la profondeur à laquelle est ~V, placé le détecteur sur le fond marin, X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur et V la vitesse moyenne de propagation de l'onde dans la tranche d'eau. L'onde émise au point d'émission 12 peut également se propager du point d'émission 12 au détecteur 10 après des réflexions sur le fond marin et sur la surface de la mer, on parle alors d'onde multiple. La première onde multiple 16 correspond à une réflexion sur le fond marin et une réflexion sur la surface 15 de l'eau avant d'arriver au détecteur. Comme il apparaît sur la figure 1, lorsque le fond marin 11 peut être considéré comme un plan horizontal le temps de propagation de la première onde multiple entre le point d'émission et le détecteur vérifie l'équation 20 suivante : _ 3z 2 (X\z t2 =1 T + ti ~ \Vi, avec Z la profondeur à laquelle est placé le détecteur sur le fond marin, X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur et V la vitesse moyenne de 25 propagation de l'onde dans la tranche d'eau. D'après les équations vérifiées par tain et tmul, il est possible d'établir les relations suivantes : Ao = tmzul - =8to (équation 1) et 01 = 9t z1ir - tmul 2~X `2 V ~ (équation 2). -11- Selon un mode de réalisation de l'invention, l'équation (1) permet de déterminer le temps de propagation vertical à l'aplomb du détecteur 10. Avantageusement, cette approche ne nécessite aucune connaissance a priori de la vitesse de propagation dans l'eau. L'équation (2) conduit à la détermination de la vitesse de propagation de l'onde dans l'eau. Selon un mode de réalisation de l'invention, pour estimer la quantité Ao, le procédé selon l'invention peut comprendre les étapes suivantes : les traces sismiques s(t) sont ré-échantillonnées selon la variable T = t2, puis la quantité Lo est obtenue en pointant le pic principal de l'auto-corrélation du signal s(T). Avantageusement, le procédé selon l'invention est indépendant de la phase du signal correspondant à l'enregistrement de l'onde reçue par le détecteur. Selon un mode de réalisation de l'invention, la quantité ni peut être réalisée en suivant un procédé analogue avec un ré-échantillonnage selon la variable T = 9t2 et une cross-corrélation avec le signal s(T=t2). Selon un mode de réalisation de l'invention, la quantité Al peut être obtenue au moyen de la relation =8(Tdfr -To+to)z-Ao, avec Tdir le temps de propagation de l'onde directe mesuré sur l'enregistrement de l'onde reçue par le détecteur, to le temps de propagation vertical déterminé au moyen de no et To le temps de propagation vertical estimé. En effet, d'après les définitions de no et Ol, il apparaît que A'=81lr-Do L'analyse de no conduit à déterminer sans erreur le temps de propagation vertical to. On effectue ensuite un pointé classique du temps d'arrivée directe pointé Tdir sur l'enregistrement de l'onde reçue par le détecteur s(t). De 2963833 -12- préférence, ce pointé est cohérent entre les traces enregistrées par un même détecteur et par différents détecteurs (même technique de pointé). Ce temps d'arrivée directe pointé Tdir est connu à 5 une constante près en fonction du détecteur. En suivant la procédure décrite dans le document FR 10 52600 ce temps peut être utilisé pour déterminer avec précision la position du détecteur, mais également pour fournir une estimation du temps à offset nul To. 10 En général, ce temps est biaisé car l'origine temporelle des données est inexacte, et le pointé dépend de la phase du signal source. En comparant To et to on obtient le décalage temps dTo= To - to du détecteur, ce qui permet d'en déduire de 15 façon non biaisée la quantité 41 en prenant tdir = Tdir -dTo. La figure 2a représente une projection selon un plan vertical défini par le détecteur 10 disposé sur la surface 11 du fond marin et du point d'émission 12 disposé à la surface de la mer. Comme représenté sur la figure 2a, 20 la bathymétrie peut être plus complexe qu'un fond marin plan. La figure 2b représente une vue du dessus du détecteur 10 et du point d'émission 12. Les inventeurs ont observé que pour un détecteur 25 donné les variations de 4o et 41 avec la distance horizontale entre le détecteur et le point d'émission sont respectivement décrites par les modèles plan et quasi-parabolique suivants Ao°d =D+Bax+Byy (équation 3), 30 AT0d = A(x2 + y2)+FFx+Fyy+E (équation 4 ) , Selon un mode de réalisation, pour un détecteur donné, on recherche donc les surfaces paramétriques définies par les équations 3 et 4. A partir des paramètres 2963833 - 13 - définissant ces surfaces, il est possible d'obtenir une estimation : du temps vertical to, du plan de réflexion défini en particulier par les angles y, q et la cote dz, et de la vitesse moyenne V dans la tranche d'eau au moyen des 5 relations suivantes . D B ~ Vto.~Bx + By v = tan .\/8D-V2(BX +Bÿ) X70 -tan ~B I Yo D , (équation 5) = 2fi V = tan-' ~~ F ~F 2 = sin 2E 3AGio (équation 6) -NT -À AVto Le plan de réflexion est défini comme étant le plan 10 tangent à la surface du fond marin à l'endroit où est placé le détecteur. Comme illustré par la figure 2b, lorsque le fond marin n'est plus modélisé comme un plan horizontal, l'onde multiple peut se propager en dehors du plan vertical 15 comprenant le trajet de l'onde directe. Selon un mode de réalisation de l'invention, pour un détecteur placé sur le fond marin on commence par déterminer pour chaque point d'émission les 8o, par exemple par auto-corrélation. 20 Il est possible de déterminer les paramètres D, Bx et By de l'équation 3 de sorte que Lô°d se rapproche le plus de l'ensemble des à0. L'homme du métier peut utiliser toute méthode d'inversion connue permettant d'ajuster les paramètres D, Bx et By de sorte à obtenir la meilleure 25 correspondance possible avec les Ao déterminés pour chaque couple détecteur point d'émission. A partir des valeurs des paramètres D, Bx et By, il est possible de déterminer la valeur du temps de -14- propagation vertical au moyen de l'équation : D t0 _ . j8D-V2(B2 +Bv) (équation 5). Une estimation de la vitesse de propagation V est suffisante pour déterminer to car, dans l'équation 5, la 5 vitesse est du second ordre. Une première évaluation de temps to permet de calculer sans erreur Al et d'en déduire une estimation de la vitesse au moyen de l'équation 6. La méthode peut également être utilisée en boucle pour 10 réévaluer to si la valeur a priori de V était éloignée de la valeur déterminée a posteriori. L'exploitation des équations 5 et 6 permet deux estimations indépendantes de l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est 15 placé le détecteur et la direction verticale et de l'angle i entre la projection de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe des abscisses dans le repère centré sur le détecteur. 20 La cohérence des estimations peut être utilisée comme critère de qualité du procédé selon l'invention. Comme illustré par la figure 3, lorsqu'une pluralité de détecteurs est disposée sur la surface des fonds marins, le procédé selon l'invention permet de déterminer un 25 décalage temporel dTo pour chaque détecteur placé sur la surface des fonds marins. Il apparaît également sur la figure 3 que la correction des vitesses moyennes par la méthode selon l'invention permet un gain dans la cohérence des résultats 30 et permet également de détecter les détecteurs présentant une anomalie de fonctionnement 20. -15- La figure 4 illustre le gain en cohérence de la correction de l'origine des temps obtenu grâce à un procédé selon l'invention. Les enregistrements sont corrigés au moyen d'une 5 méthode de type move-out (en anglais « Move-out x2 +y2 +z2 correction ») t- ~ avec x, y, z les coordonnées du point d'émission de l'onde dans un repère centré sur chaque détecteur et l'origine des temps est déterminé au moyen du temps de propagation vertical to. Le temps to et la vitesse 10 de propagation V sont déterminés au moyen d'un procédé selon l'invention. L'invention ne se limite pas aux modes de réalisation décrits et doit être interprétée de façon non limitative, en englobant tout mode de réalisation équivalent. 15 REVENDICATIONS1. Procédé de détermination pour un détecteur disposé sur la surface du fond marin du temps de propagation verticale (to) et de la vitesse de propagation V dans l'eau d'au moins une onde émise à partir d'au moins un point d'émission parmi N points d'émission, le procédé comprenant les étapes suivantes pour au moins un des N points d'émission: - émission d'au moins une onde à partir dudit point d'émission, enregistrement de l'onde reçue par le détecteur, détermination du temps de propagation verticale to au moyen de la relation suivante : Ao = tmud -td2;r = Sto , avec td, r le temps de propagation de l'onde directe entre le point d'émission et le détecteur et tmul le temps de propagation de la première onde multiple entre le point d'émission et le détecteur, - détermination de la vitesse de propagation V de l'onde au moyen de la relation suivante : _ ~ Al 9tdir = - tmulavec X la distance entre le point d'émission et le point à la surface de l'eau à la verticale du détecteur. 2. Procédé selon la revendication 1, dans lequel le temps de propagation vertical to est déterminé par les étapes suivantes : ré-échantillonnage de l'enregistrement de l'onde 30 reçue par le détecteur selon le changement de variable T=t2 auto-corrélation du signal s(T) correspondant à-17- l'enregistrement ré-échantillonné de l'onde reçue par le détecteur détermination de la quantité ,o à partir de la valeur de T correspondant au pic principal du signal auto- corrélé. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel la vitesse de propagation est déterminée par les étapes suivantes . ré-échantillonnage de l'enregistrement de l'onde reçue par le détecteur selon le changement de variable T=9t2, corrélation du signal s(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=9t2 avec le signal s'(T) correspondant à l'enregistrement ré-échantillonné selon le changement de variable T=t2, détermination de la quantité Al à partir de la valeur de T correspondant au pic principal du signal corrélé. 4. Procédé selon l'une des revendications 1 ou 2, dans lequel la quantité Al est déterminée au moyen de la relation suivante =8td;,' -Ao tdir étant déterminé au moyen du temps d'arrivée directe Tdir pointé sur l'enregistrement de l'onde reçue par le détecteur, ledit temps d'arrivée directe pointé Tdjr corrigé du décalage temporel dTo=To-to avec To le temps de propagation vertical estimé et to le temps de propagation vertical déterminé au moyen de Lo. 5. Procédé selon l'une quelconque des revendications 2963833 -18- précédentes dans lequel, à partir des Do déterminés pour chacun des N points d'émission, on modélise les Do obtenus pour chacun des points d'émission au moyen de l'équation suivante : flood = D+Bxx+Byy avec x et y 5 les coordonnées horizontales de chaque point d'émission dans un repère centré sur le détecteur et avec D, BX et By des paramètres de modélisation dont les valeurs sont déterminées de sorte que A0d s'ajuste le mieux aux Do obtenus, par exemple au sens 10 des moindres carrés, le temps de propagation vertical (to) est déterminé au moyen de la relation suivante: to v18D_v20x 2 +By) 6. Procédé selon l'une quelconque des revendications 15 précédentes, dans lequel à partir des Al déterminés pour chacun des N points d'émission, on modélise les Ai obtenus pour chacun des points d'émission au moyen de l'équation suivante : An' =A(x' +y2)+Fxx+Fvy+E avec x et y les coordonnées horizontales de chaque point 20 d'émission dans un repère centré sur le détecteur, et avec A, Fx, Fy et E des paramètres de modélisation dont les valeurs sont déterminées de sorte que A7d s'ajuste le mieux aux A obtenus, par exemple au sens des moindres carrés, 25 la vitesse de propagation (V) est déterminée au moyen de la relation suivante: V =2N/i 7. Procédé selon la revendication 5 dans lequel, l'angle y entre la normale au plan tangent à la surface du 30 fond marin à l'endroit où est placé le détecteur et D-19- la direction verticale et/ou l'angle n entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : rB y Bx 1 Ii =tan-1 , y= tan 7Vto.\/Bx + Bÿ D 8. Procédé selon la revendication 6 dans lequel, l'angle y entre la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et la direction verticale et/ou l'angle n entre la projection sur le plan horizontal (X, Y) de la normale au plan tangent à la surface du fond marin à l'endroit où est placé le détecteur et l'axe (X) dans le repère centré sur le détecteur, est déterminé par : Il =tan '(Fy Y Fx, y=sin' iVF2 +Fy AVto J 9. Procédé de correction de l'origine des temps des signaux enregistrés par un détecteur disposé sur la surface d'un fond marin, lendit signaux correspondant à des ondes émises à partir d'au moins un point d'émission parmi N points d'émission, dans lequel les signaux reçus par le détecteur sont corrigés au moyen d'un procédé de type move-out et en tenant compte de l'origine des temps avec le temps to et la vitesse de propagation V déterminés au moyen d'un procédé selon l'une quelconques des revendications précédentes.30-20- 10. Produit de programme d'ordinateur comprenant une série d'instructions qui lorsque chargée dans un ordinateur entraîne l'exécution par ledit ordinateur des étapes du procédé selon l'une quelconque des 5 revendications précédentes.