Brevet FR2797501A1 — publié le 20010216

L'invention concerne l'exploration sismique marine, et plus particulièrement la détermination des coefficients de réflexion du fond de l'eau à partir de sismogrammes enregistrés pour une utilisation dans la déréverbération de données sismiques sans endommager de façon nuisible des séquences de réflexion qui ont une périodicité proche de celle du temps d'aller et retour dans l'eau d'ondes sismiques. Dans l'exploration sismique marine, de l'énergie sismique est générée par des sources remorquées derrière un navire. Les ondes de pression acoustique résultantes sont réfléchies depuis les formations souterraines et enregistrées sous forme de sismogrammes. Habituellement, les sismogrammes enregistrés contiennent des bruits qui masquent les caractéristiques des formations souterraines. Ces bruits comprennent des réflexions multiples du fond de l'eau comme illustré sur la figure 2 des dessins annexés et décrits ci-après et des réverbérations comme illustré sur les figures 3 et 4 de ces mêmes dessins. Pendant que l'impulsion de tir initial se déplace vers le fond de l'eau, elle est partiellement transmise et partiellement réfléchie par le fond de l'eau. La partie qui est transmise est réfléchie par des discontinuités situées au-dessous de l'eau pour constituer des réflexions primaires qui sont importantes dans la formation d'images du sous-sol. La partie réfléchie revient par rebond vers la surface o elle est de nouveau réfléchie vers le fond de l'eau. Une partie de cette impulsion passe de nouveau à travers la limite et une partie est réfléchie. Le processus se poursuit indéfiniment. A chaque fois qu'une impulsion rebondit entre la surface et le fond de l'eau, une impulsion émerge de la limite située au fond pour suivre l'impulsion initiale dans la section. Ces impulsions multiples sont retardées par des intervalles de temps égaux les unes par rapport aux autres, et leurs amplitudes sont diminuées successivement par le coefficient de réflexion au fond de l'eau. Un processus similaire se produit avec les réflexions primaires qui rebondissent dans la couche d'eau. Le résultat net en est la production d'un sismogramme constitué d'une réflexion primaire unique suivie d'un train d'impulsions multiples ou de réverbération espacées de façon égale. La trace sismique qui est générée par une impulsion de tir peut être décrite mathématiquement sous la forme (1) X(t) = S(t) * R(t) * H(t) + N(t) o X(t) est le signal enregistré, S(t) est le signal de la source, R(t) est une séquence de coefficients de réflexion caractérisant le sous-sol, H(t) est une fonction de distorsion et N(t) est un bruit additif. Le bruit peut être cohérent ou incohérent. La fonction de distorsion qui caractérise des réverbérations est donnée par l'équation (2) o gR (t) = ()+ ôfr)k Jet - kT) k=i o H R(t) est l'opérateur de réverbération, r est le coefficient de réflexion au fond de l'eau, 8 est la fonction 6 de Kronecker, T est le temps de parcours d'aller et retour dans l'eau pour des ondes acoustiques. Dans le domaine de la transformée en z, l'opérateur de déréverbération peut être désigné par (3) HD (Z) I (+ rZ) Dans les cas o le fond de l'eau est relativement plat, la distorsion de réverbération descendante au point de tir, et la distorsion de réverbération montante au détecteur sont approximativement les mêmes. Par conséquent, pour un fond de l'eau plat, l'opérateur de déréverbération total peut être représenté par la convolution de l'opérateur en deux points au point de tir avec l'opérateur en deux points au détecteur. Ceci donne l'opérateur dit de Backus (4) B(z)=(1 rZ)2 (l+ rZ)' Le brevet des Etats-Unis d'Amérique N 3 689 874 de Foster et al. décrit un procédé pour la déréverbération du signal enregistré, utilisant l'autocorrélation des traces enregistrées. Le procédé suppose que l'autocorrélation du coefficient de réflexion et des séquences de bruit peut être caractérisée comme étant un bruit blanc. Le brevet des Etats-Unis d'Amérique N 4 146 871 de Ruehle décrit un procédé pour la déréverbération de signaux enregistrés, dans lequel l'opération de déréverbération est déterminée par l'essai de différentes valeurs pour le coefficient de réflexion r du fond de l'eau et le temps de parcours en aller et retour dans l'eau T et, après que la meilleure valeur a été déterminée d'après une inspection visuelle des traces traitées, les canons à air de la source sismique sont déclenchés séquentiellement pour simuler l'opérateur de déréverbération. Le déploiement de capteurs sur le fond océanique est souvent utilisé pour obtenir des données sismiques en trois dimensions dans une zone présentant des obstacles qui s'opposent à l'utilisation de rubans remorqués. Une telle situation apparaît, par exemple, lors de l'acquisition de données à proximité de plates-formes de production. Avec des capteurs sur le fond océanique, le "fantôme" affectant les données de capteurs enregistrées dû à une réflexion depuis la surface de l'eau apparaît davantage qu'avec des données de rubans, les rubans étant habituellement à des profondeurs inférieures à 10 mètres tandis que les capteurs du fond océanique peuvent être à des profondeurs de plusieurs dizaines ou centaines de mètres d'eau. Les dispositifs de capteurs multiples comprennent des capteurs de pression (hydrophones) et des capteurs de composantes verticales. Ces derniers sont habituellement des capteurs de vitesse (géophones). Si l'on connaît le coefficient de réflexion du fond de l'eau, les signaux de pression et de vitesse peuvent être combinés pour donner le champ arrivant en montant. Le brevet des Etats-Unis d'Amérique N 5 524 100 délivré à Paffenholz décrit un agencement dans lequel des dispositifs de capteurs multiples sont déployés sur le fond de la mer. Des signaux sismiques de pression et de vitesse sont combinés, le signal combiné est transformé en le domaine de fréquence et multiplié par l'opérateur de Backus inverse, ou bien le signal combiné est soumis à une convolution avec l'opérateur de Backus inverse, et un algorithme d'optimisation est utilisé pour résoudre le coefficient de réflexion du fond de l'eau. Les signaux sismiques de pression et de vitesse sont combinés, et le signal combiné est multiplié par l'opérateur de Backus inverse contenant le coefficient de réflexion du fond de l'eau pour éliminer des multiples en baïonnette du premier ordre. L'un des problèmes associés aux procédés de l'art antérieur est que la détermination de l'opérateur de déréverbération peut être troublée par des séquences de coefficient de réflexion géologiques qui ont un temps de retard égal au temps d'aller et retour dans l'eau. Les procédés de traitement sont incapables de faire la distinction entre des encoches spectrales provoquées par des réverbérations et des encoches spectrales provoquées par des réflexions souterraines. Par conséquent, l'opérateur de déréverbération supprime également des séquences de réflexion qui ont le même temps de retard que l'opérateur de Backus et les données sismiques traitées résultantes sont incorrectes par le fait que certaines réflexions associées à la géologie souterraine sont supprimées, ce qui conduit à des mauvaises interprétations possibles des données. Il serait souhaitable de disposer d'un procédé pour la déréverbération de données sismiques marines qui soit robuste en présence d'une telle interférence géologique. L'invention satisfait à ce besoin. L'invention concerne un procédé de détermination du coefficient de réflexion du fond de l'eau, qui est insensible à une interférence géologique, définie ici sous la forme de réflecteurs souterrains ayant des temps de retard proches de ceux dans la réverbération du fond de l'eau. Les temps de parcours de fond de l'eau sont obtenus par les procédés classiques. Des données sismiques de déports proches sont ensuite soumises à une déréverbération pour une série de coefficients de réflexion possible du fond de l'eau afin de donner une série de traces filtrées. Le coefficient de réflexion du fond de l'eau qui procure la valeur minimale sur la série de traces filtrées est déterminé pour donner un ensemble de coefficient de réflexion choisis pour chacun des points d'échantillonnage dans une fenêtre de temps. Une pondération est déterminée pour chaque temps dans la fenêtre de temps, la pondération étant liée à la variation de l'amplitude de la trace filtrée en ce point d'échantillonnage sur la série de coefficients de réflexion du fond de l'eau. Les coefficients de réflexion choisis sont pondérés par les pondérations et soumis à une moyenne pour réaliser une détermination du coefficient de réflexion du fond de l'eau. Ce coefficient de réflexion déterminée du fond de l'eau est utilisé pour le traitement de la trace sismique enregistrée. L'invention sera décrite plus en détail en regard des dessins annexés à titre d'exemple nullement limitatif et sur lesquels: la figure 1 illustre schématiquement un exemple de système d'acquisition de données sismiques destiné à être utilisé dans un environnement marin; la figure 2 est un diagramme illustrant des réflexions multiples du fond de l'eau; les figures 3 et 4 sont des diagrammes illustrant des réverbérations dans l'eau; les figures 5 à 8 sont des diagrammes illustrant des opérateurs inverses dans le domaine temporel à deux points et trois points; et les figures 9A et 9B illustrent des étapes de la présente invention pour la détermination du coefficient de réflexion du fond de l'eau. En référence à présent à la figure 1, celle-ci illustre un exemple de parties d'un système d'acquisition de données sismiques marines. Un navire 10 flottant sur une masse d'eau 15 s'étendant au-dessus de la terre 16 a déployé derrière lui un dispositif 20 de sources sismiques et un câble 25 de ruban. Le dispositif 20 de sources sismiques est habituellement constitué de canons à air individuels 20a, 20b...20n qui sont déclenchés sous la commande d'un dispositif de commande (non représenté) embarqué à bord du navire 10. Des impulsions sismiques se propagent dans le sol et sont réfléchies par un réflecteur 22 s'y trouvant. Pour simplifier l'illustration, un seul réflecteur est représenté; en réalité, il peut y avoir de nombreux réflecteurs, donnant chacun naissance à une impulsion réfléchie. Après la réflexion, ces impulsions reviennent à la surface o elles sont enregistrées par des détecteurs (hydrophones) 30a, 30b,...30n dans le câble de ruban. La profondeur du dispositif de sources et du câble de ruban est commandée par des dispositifs auxiliaires (non représentés). La figure 2 montre des réflexions multiples du fond de l'eau se déplaçant dans la masse d'eau, qui seraient également détectées par les détecteurs 30a, 30b,...30n. Les figures 3 et 4 montrent des réverbérations d'impulsions revenant par réflexion du sol, provoquées par des réflexions multiples dans la couche d'eau. La trace sismique qui est générée par une impulsion de tir peut être décrite mathématiquement par l'équation (1). La fonction de distorsion introduite par une réflexion unique dans la couche d'eau est illustrée par les figures 5 et 6. L'impulsion initiale à un temps 0 est suivie d'une impulsion d'amplitude r à un temps Tw, o Tw est le temps d'aller et retour dans l'eau pour une impulsion sismique et r est le coefficient de réflexion du fond de l'eau. Ce coefficient de réflexion du fond de l'eau peut être positif ou négatif, comme montré sur les figures 5 et 6. Ce coefficient de réflexion du fond de l'eau combiné avec des réverbérations multiples dans la couche d'eau a pour effet d'introduire une fonction de réverbération. Dans des conditions dans lesquelles le fond de l'eau est relativement plat, la distorsion de réverbération dans le sens descendant, à la source, et la distorsion de réverbération, dans le sens montant, au détecteur sont approximativement les mêmes. Par conséquent, pour un fond d'eau plat, un opérateur de déréverbération total peut être représenté par la convolution de l'opérateur à deux points à la source avec l'opérateur à deux points au détecteur. Ceci donne l'opérateur dit de Backus indiqué par l'équation (4). Le temps d'aller et retour dans l'eau peut être aisément obtenu par des procédés antérieurs. L'un des procédés les plus précis consiste à le calculer par bathymétrie, car la vitesse du son dans l'eau est bien connue. Un autre procédé consiste à le déterminer d'après le temps d'arrivée de la réflexion du fond de l'eau au premier détecteur 30a et à le calculer à partir de valeurs connues de la profondeur du détecteur et de la profondeur de la source au-dessous de la surface 18 de l'eau. Cependant, la détermination du coefficient de réflexion du fond de l'eau est problématique. Comme indiqué précédemment, les procédés basés sur les spectres ou sur l'autocorrélation de la trace sismique enregistrée sont sujets à une erreur due à des séquences géologiques avec la même période de temps que le temps d'aller et retour dans l'eau. Une détermination directe à partir de l'amplitude de la réflexion du fond de l'eau est difficile à réaliser car, en eau peu profonde, la réflexion du fond de l'eau peut arriver dans la même fenêtre de temps que celle d'une arrivée directe à travers l'eau, de la source au détecteur. Les figures 9A et 9B auxquelles on se réfère à présent illustrent un organigramme de la présente invention pour la détermination du coefficient de réflexion du fond de l'eau. Le procédé commence avec un temps d'aller et retour Tw dans l'eau 100 et une trace sismique enregistrée 102. La trace enregistrée comprend des impulsions de réflexion provenant de réflecteurs situés au-dessous du fond de l'eau et comprend aussi des réverbérations indésirables qui doivent être atténuées. Ce temps d'aller et retour est double du temps de parcours de la surface de l'eau jusqu'au fond de l'eau. Une plage de valeurs possibles pour le coefficient de réflexion r du fond de l'eau est choisie. Dans une forme appréciée de réalisation de l'invention, la plage de coefficient de réflexion possible du fond de l'eau est de -0,9 à +0,9. Une valeur initiale du coefficient de réflexion est choisie en 104 à la limite inférieure de la plage des coefficients de réflexion. En utilisant cette valeur initiale du coefficient de réflexion et le temps d'aller et retour dans l'eau, on applique un opérateur de Backus donné par l'équation 4 à la trace enregistrée, en 106. Une trace filtrée est sauvegardée en 112 en tant que partie d'une série de trace filtrée. La valeur possible du coefficient de réflexion est augmentée d'une quantité Ar 108 et un contrôle est réalisé pour déterminer si la valeur possible dépasse la limite supérieure de la plage de coefficients de réflexion possibles 110. Dans une forme appréciée de réalisation de l'invention, la valeur de Ar est de 0,01. Si le test effectué en 110 est négatif le processus revient à l'étape 104 avec la nouvelle valeur du coefficient de réflexion possible et les étapes 106 et 110 sont répétées, en sauvegardant à chaque fois une nouvelle trace après déréverbération en tant que partie d'une série de traces 112 ayant subi une déréverbération. Si le test effectué en 110 est positif la boucle d'itération allant de 106 à 110 est achevée et la série de traces filtrées en 114 est utilisée dans des étapes suivantes illustrées sur la figure 9B. Une fenêtre de temps pour l'analyse des données est sélectionnée. Comme cela est connu des spécialistes de la technique, le point de départ de la fenêtre ne doit pas contenir de bruits initiaux portés dans l'eau et d'autres artéfacts. Au temps de départ de la fenêtre 116, la série entière de traces ayant subi une déréverbération est examinée pour trouver les amplitudes minimale et maximale de la série de traces filtrées. Les minimum et maximum sont délivrés en sortie à une table mémorisée de valeurs en même temps que le coefficient de réflexion possible correspondant à l'amplitude minimale à l'absolu. C'est la valeur du coefficient de réflexion qui, lorsqu'elle est utilisée dans un opérateur de déréverbération, présente la réduction la plus importante d'amplitude de la trace. Le temps pour l'analyse est incrémenté d'une quantité At en 120. Il est vérifié si le temps incrémenté dépasse le temps maximal de la fenêtre d'analyse en 122. Ce temps maximal de la fenêtre d'analyse peut être choisi sur la base de l'un de plusieurs critères. Dans une forme de réalisation de l'invention, le temps maximal de la fenêtre d'analyse est basé sur un rapport signal/bruit approprié; il est évident qu'il est indésirable de sélectionner un opérateur de déréverbération provenant de bruits. En variante, le temps maximal d'analyse est basé sur des considérations géologiques; dans la plupart des zones d'exploration auxquelles on s'intéresse, on connaît à priori un temps de parcours approximatif pour l'intervalle géologique auquel on s'intéresse. L'incrément du temps d'analyse At est choisi de façon à être un multiple de l'intervalle de temps d'échantillon At utilisé pour l'enregistrement des données. Dans une forme appréciée de réalisation de l'invention, ils sont égaux, c'est-à-dire que l'incrément de temps d'analyse est égal à l'incrément de temps d'échantillonnage des données sismiques. Si le test effectué en 122 est négatif, le temps d'analyse est incrémenté et les étapes 116 à 122 sont répétées. A chaque temps d'analyse, les amplitudes absolues maximale et minimale de la série de traces ayant subi une déréverbération et le coefficient de réflexion possible correspondant à l'amplitude absolue minimale des traces ayant subi une déréverbération sont sauvegardés. Une fois que le temps d'analyse maximal est dépassé, comme indiqué par une réponse positive en 122, on procède à une analyse des valeurs mémorisées, en 124. L'analyse des valeurs mémorisées est basée sur l'analyse raisonnée suivante: chacune des valeurs mémorisées du coefficient de réflexion possible (une pour chaque temps d'analyse) est la valeur qui réduit à un minimum l'amplitude de la trace ayant été soumise à une déréverbération. C'est ce qu'un opérateur de déréverbération devrait faire s'il n'y avait aucune impulsion sismique réfléchie au temps de l'analyse. Cependant, s'il y a une impulsion de réflexion sismique au temps de l'analyse, une réduction à une valeur minimale ne convient pas. L'établissement d'une moyenne sur une plage de temps d'analyse est une moyenne établie sur différentes impulsions géologiques réfléchies et réduirait donc l'effet de la géologie sur le coefficient de réflexion déterminé. Cependant, au lieu d'un simple calcul de moyenne, la présente invention pondère le coefficient de réflexion choisi par un facteur dépendant de la plage de valeurs de coefficient de réflexion de la série de traces ayant étant soumises à une déréverbération: ceci signifie qu'une pondération plus grande est donnée à un coefficient de réflexion choisi ayant un effet plus grand. Les spécialistes de la technique reconnaîtront que d'autres fonctions de pondération pourraient également être utilisées. Ces variantes sont connues des spécialistes de la technique et entendent être incluses dans le cadre de l'invention. Le brevet N 5 524 100 précité, dont le contenu est incorporé ici totalement à titre de référence, décrit une invention dans laquelle un dispositif de capteurs doubles déployés sur le fond de l'océan est utilisé pour déterminer le coefficient de réflexion du fond de l'eau. Le dispositif de capteurs doubles comporte un géophone pour composante verticale et un hydrophone. Le coefficient de réflexion est déterminé en minimisant la puissance dans une bande spectrale. Le coefficient de réflexion déterminé est utilisé pour déterminer un champ d'onde arrivant en montant qui est exempt de réverbérations. Comme indiqué précédemment, la géologie contamine le coefficient de réflexion déterminé. Le coefficient de réflexion du fond de l'eau déterminé dans la présente invention et décrit ci-dessus en référence aux figures 9A et 9B peut être utilisé conjointement avec la divulgation du brevet N 5 524 100 précité pour déterminer un champ d'onde arrivant en montant, ayant été soumis à une déréverbération, qui n'est pas sensible à une contamination géologique de la détermination du coefficient de réflexion. 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é de détermination d'un coefficient de réflexion du fond (17) d'une masse d'eau (15), caractérisé en ce qu'il comprend: (a) le déclenchement d'au moins une source sismique (20) dans la masse d'eau (15) pour générer des impulsions sismiques qui se propagent dans la masse d'eau et dans le sol (16) situé au-dessous d'elle; (b) l'enregistrement d'au moins une trace sismique représentative d'impulsions sismiques réfléchies par des réflecteurs sismiques (22) dans le sol, cette trace sismique comprenant des réverbérations dans la masse d'eau des impulsions sismiques réfléchies par les réflecteurs sismiques; (c) l'obtention d'un temps de parcours d'une surface supérieure (18) de la masse d'eau jusqu'au fond de la masse d'eau; (d) pour chacun d'une première pluralité de coefficients de réflexion possibles du fond de la masse d'eau, l'obtention d'une trace filtrée à partir du temps de parcours, de la trace sismique ou des traces sismiques, et du coefficient de réflexion possible pour donner une première pluralité de traces filtrées; (e) pour chacun d'une seconde pluralité de temps d'analyse dans une fenêtre de temps, la détermination d'un facteur de qualité à partir de valeurs de la première pluralité de traces filtrées à chaque temps d'analyse pour donner une second pluralité de facteurs de qualité; et (f) la détermination du coefficient de réflexion au fond de la masse d'eau à partir de la première pluralité de traces filtrées et de la seconde pluralité de facteurs de qualité. 2. Procédé selon la revendication 1, caractérisé en ce que la ou chaque source sismique comporte en outre plusieurs sources sismiques (20a...20n) dans un dispositif (20). 3. Procédé selon la revendication 1, caractérisé en ce que l'enregistrement de la ou de chaque trace sismique est effectué en utilisant au moins l'un (i) d'un hydrophone dans l'eau, (ii) d'un hydrophone au fond de la masse d'eau et (iii) d'un géophone au fond de la masse d'eau. 4. Procédé selon la revendication 1, caractérisé en ce que chacun de la première pluralité de coefficients de réflexion possibles est compris entre -0,9 et 0,9. 5. Procédé selon la revendication 1, caractérisé en ce que la première pluralité de traces ayant été soumises à une déréverbération est obtenue à partir de la ou de chaque trace sismique par un opérateur de déréverbération de la forme 1 + 2R6(T)+2 R26(2T) o R est l'un de la première pluralité de coefficients de réflexion possibles, T est le temps de parcours obtenu et 6 est la fonction delta de Kronecker. 6. Procédé selon la revendication 1, caractérisé en ce que chaque facteur de qualité est une valeur absolue de la différence entre une valeur maximale de la première pluralité de traces filtrées au temps d'échantillonnage et une valeur absolue minimale de la première pluralité de traces filtrées au temps d'analyse. 7. Procédé selon la revendication 1, caractérisé en ce que la détermination du coefficient de réflexion comprend en outre: (i) la détermination d'un maximum des traces filtrées en chacun de la seconde pluralité de temps d'analyse pour donner une seconde pluralité de coefficients de réflexion choisis; (ii) la détermination d'une seconde pluralité de coefficients de pondération pour la seconde pluralité de coefficients de réflexion choisis, chaque coefficient de pondération étant associé à une plage de valeurs de la première pluralité de traces filtrées à chacun de la seconde pluralité de temps d'analyse; et (iii) la détermination du coefficient de réflexion sous la forme d'une moyenne de la seconde pluralité de coefficients de réflexion choisis pondérés par le coefficient de pondération correspondant de la seconde pluralité de coefficients de pondération. 8. Procédé selon la revendication 7, caractérisé en ce que chacun de la seconde pluralité de coefficients de pondération est une différence entre un maximum et un minimum de la première pluralité de traces ayant été soumises à une déréverbération en chacun de la seconde pluralité de temps d'analyse. 9. Procédé de prospection sismique, caractérisé en ce qu'il comprend: (a) le déclenchement d'au moins une source sismique (20a...20n) dans une masse d'eau (15) pour générer des impulsions sismiques qui se propagent dans la masse d'eau et dans le sol (16) situé au-dessous d'elle; (b) l'enregistrement d'au moins une trace sismique représentative d'impulsions sismiques réfléchies par des réflecteurs sismiques (22) dans le sol, la ou chaque trace sismique comprenant des réverbérations d'impulsions sismiques réfléchies par les réflecteurs sismiques; (c) l'obtention d'un temps de parcours depuis une surface supérieure (18) de la masse d'eau jusqu'à un fond (17) de la masse d'eau; (d) pour chacun d'une première pluralité de coefficients de réflexion possibles du fond de la masse d'eau, l'obtention d'une trace filtrée à partir du temps de parcours, de la ou de chaque trace sismique et du coefficient de réflexion possible pour donner une première pluralité de traces filtrées; (e) pour chacun d'une seconde pluralité de temps d'analyse dans une fenêtre de temps, la détermination d'un facteur de qualité à partir de valeurs de la première pluralité de traces filtrées à chaque temps d'analyse pour donner une seconde pluralité de facteurs de qualité; (f) la détermination d'un coefficient de réflexion au fond de la masse d'eau à partir de la première pluralité de traces filtrées et de la seconde pluralité de facteurs de qualité; et (g) le traitement de la ou de chaque trace sismique en utilisant le coefficient de réflexion déterminé pour donner une trace sismique traitée dans laquelle les réverbérations ont été réduites. 10. Procédé selon la revendication 9, caractérisé en ce que la ou chaque source sismique comprend en outre plusieurs sources sismiques (20a...20n) dans un dispositif (20). 11. Procédé selon la revendication 9, caractérisé en ce que l'enregistrement de la ou de chaque trace sismique est effectué en utilisant au moins l'un (i) d'un hydrophone dans l'eau, (ii) d'un hydrophone au fond de la masse d'eau, et (iii) d'un géophone au fond de la masse d'eau. 12. Procédé selon la revendication 9, caractérisé en ce que chacun de la seconde pluralité de coefficients de réflexion possibles est compris entre -0,9 et 0,9. 13. Procédé selon la revendication 9, caractérisé en ce que la première pluralité de traces ayant été soumises à une déréverbération est obtenue à partir de la ou chaque trace sismique par un opérateur de déréverbération de la forme 1 + 2R6(T)+2 R26(2T) o R est l'un de la première pluralité de coefficients de réflexion possibles, T est le temps de parcours obtenu et 6 est la fonction delta de Kronecker. 14. Procédé selon la revendication 9, caractérisé en ce que chacun des facteurs de qualité est une différence entre une valeur maximale absolue de la première pluralité de traces filtrées au temps d'échantillonnage et d'une valeur minimale absolue de la première pluralité de traces filtrées au temps d'analyse. 15. Procédé selon la revendication 1, caractérisé en ce que la détermination du coefficient de réflexion comprend en outre: (i) la détermination d'un maximum des traces filtrées en chacun de la seconde pluralité de temps d'analyse pour donner une seconde pluralité de coefficients de réflexion choisis; (ii) la détermination d'une seconde pluralité de coefficients de pondération pour la seconde pluralité de coefficients de réflexion choisis, chaque coefficient de pondération étant associé à une plage de valeurs de la première pluralité de traces filtrées à chacun de la seconde pluralité de temps d'analyse; et (iii) la détermination du coefficient de réflexion sous la forme d'une moyenne de la seconde pluralité de coefficients de réflexion choisis pondérés par le coefficient de pondération correspondant de la seconde pluralité de coefficients de pondération. 16. Procédé selon la revendication 15, caractérisé en ce que le coefficient de pondération est la différence entre un maximum et un minimum de la première pluralité de traces ayant été soumises à une déréverbération à chacun des temps d'échantillonnage.