Brevet FR2890178A1 — publié le 20070302

METHODE DE TRAITEMENT D'UN ENREGISTREMENT DE TRACES SISMIQUES Un ou plusieurs modes de réalisation de la présente invention concernent généralement le traitement de données sismiques et plus particulièrement la correction des effets dissipatifs affectant les données sismiques. Dans le domaine de l'exploration sismique, les données sont obtenues en provoquant d'abord une perturbation artificielle le long de la terre au moyen de dynamite ou d'un explosif similaire. Les ondes acoustiques résultantes se propagent vers le bas dans la terre et sont réfléchies vers le haut par les interfaces réfléchissantes situées sous la surface. Les ondes réfléchies sont reçues par des détecteurs ou géophones situés au niveau du sol et enregistrées sous forme reproductible. Idéalement, les signaux enregistrés au niveau des détecteurs sont exactement représentatifs des caractéristiques de réflexion de la terre, hors présence de composantes indésirables, telles que bruit ou distorsion. Les signaux enregistrés au niveau des détecteurs contiennent malheureusement souvent de nombreuses composantes indésirables qui occultent la fonction de réflectivité de la terre et empêchent de trouver la zone de la terre susceptible de recéler des dépôts de pétrole ou de gaz. Plusieurs phénomènes contribuent à la distorsion des signaux enregistrés. Un de ces phénomènes est l'absorption, qui provoque la perte d'énergie sismique en la convertissant en d'autres formes d'énergie. Ce type de perte d'énergie sismique est généralement appelé atténuation intrinsèque. Un deuxième phénomène est l'interférence multiple intraassise. L'interférence multiple intra-assise redistribue l'énergie sismique entre les directions descendante et ascendante. Ce type de perte d'énergie sismique est généralement appelé atténuation apparente. L'atténuation apparente provoque une perte progressive des fréquences les plus élevées (en élargissant l'ondelette sismique) et une augmentation de la distorsion de phase avec le temps de propagation de l'ondelette sismique reçue. La combinaison de l'atténuation intrinsèque et apparente est généralement appelée filtre terrestre. Du fait du filtrage terrestre, rondelette sismique est variable dans le temps. L'existence d'une ondelette sismique variable dans le temps est contraire à une hypothèse élémentaire de la théorie de la déconvolution et compromet la capacité à utiliser la déconvolution pour déterminer les caractéristiques du filtre terrestre dans le cadre d'une méthode d'interprétation sismique. Une approche classique pour compenser l'atténuation par le filtre terrestre est présentée dans le document Q-Adaptive Deconvolution, par D. Hale, Stanford Exploration Project, Report 30, 1982. Hale présente deux procédures itératives pour la mise en oeuvre du filtrage Q inverse. Les procédures présentées par Hale reposent sur plusieurs hypothèses par lesquelles l'auteur arrive à une relation de dispersion approchée. L'utilisation de la relation de dispersion approchée a toutefois pour effet de dégrader la valeur de la compensation Q obtenue par Hale. Par ailleurs, l'application d'un algorithme Q inverse à décalage nul aux données avant pile (ou présentant des décalages non nuls), dans le domaine t-x conduit souvent à une compensation insuffisante du fait que des événements ayant des temps de propagation différents peuvent avoir subi des niveaux d'atténuation très différents sur leur trajet. Ce phénomène peut être particulièrement vrai en eaux profondes et dans le cas de décalages élevés. Dans ces situations, la réflexion par le fond, qui n'a subi aucune atténuation, peut coïncider en termes de temps d'arrivée avec des réflexions beaucoup plus profondes qui peuvent avoir subi une atténuation importante. De ce fait, un filtre Q unique peut ne pas suffire à compenser efficacement les réflexions. Ainsi, un besoin existe dans l'art antérieur pour une nouvelle méthode de correction des effets dissipatifs affectant les données sismiques. Divers modes de réalisation de l'invention portent sur une méthode 30 de traitement d'un enregistrement de traces sismiques. La méthode comprend la transformation de l'enregistrement de traces sismiques en un domaine à angle constant tel que chaque trace sismique correspond à un angle constant. On entend par angle constant un angle compris entre un trajet et un axe vertical à la surface. La méthode comprend en outre l'application d'un algorithme Q inverse à décalage nul dans le domaine à angle constant. Dans un mode de réalisation, l'enregistrement des traces sismiques présente une pluralité de décalages. Dans un autre mode de réalisation, l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comprend la conversion d'un modèle Q pour l'algorithme Q inverse à décalage nul en une fonction de tau. Dans un autre mode de réalisation, l'algorithme Q inverse à décalage nul est appliqué à l'enregistrement des traces sismiques au moyen 15 du modèle Q en fonction de tau. Dans un autre mode de réalisation, le modèle Q en fonction de tau est défini par dT/dt Q(to(T)), où dT représente un changement de tau, dt représente un changement du temps de propagation et to représente le temps vertical et T représente tau. Dans un autre mode de réalisation, l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comprend la conversion de l'enregistrement des traces sismiques en un domaine temps de propagation-p. Dans un autre mode de réalisation, l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comprend en outre l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques au moyen du modèle Q en fonction du temps de propagation. 2890178 4 Pour permettre la compréhension détaillée des caractéristiques ci-dessus mentionnées de la présente invention, une description plus particulière de l'invention résumée sommairement ci-dessus peut être faite par référence aux modes de réalisation, dont certains sont illustrés dans les dessins joints en annexe. Il convient toutefois de noter que les dessins joints en annexe représentent seulement des exemples de modes de réalisation de cette invention et ne doivent de ce fait pas être considérés comme limitant sa portée, l'invention pouvant admettre d'autres modes de réalisation également efficaces. La figure 1 représente une vue schématique de relevés sismiques pour lesquels divers modes de réalisation de l'invention peuvent être utilisés. La figure 2 représente une méthode de correction des effets dissipatifs de la terre affectant un enregistrement de traces sismiques conformément à un ou plusieurs modes de réalisation de l'invention. La figure 3 représente une autre méthode de correction des effets dissipatifs de la terre affectant un enregistrement de traces sismiques conformément à un ou plusieurs modes de réalisation de l'invention. La figure 4 représente un réseau informatique dans lequel les modes de réalisation de l'invention peuvent être mis en oeuvre. Un ou plusieurs modes de réalisation de l'invention peuvent être utilisés pour diverses explorations sismiques, telles que les explorations sismiques marines, les explorations sismiques terrestres, les explorations sismiques du fond de la mer, les explorations sismiques en trou de sonde, et d'autres explorations similaires. La figure 1 représente une vue schématique marine d'explorations sismiques 100 pour lesquelles divers modes de réalisation de l'invention peuvent être utilisés. Les formations sous-terraines à explorer, telles que 102 et 104, se trouvent sous une masse d'eau 106. Les sources d'énergie sismique 108 et les récepteurs sismiques 110 sont positionnés dans la masse d'eau 106, en général par un ou deux navires sismiques (non représentés). Une source sismique 108, un canon à air par exemple, crée des ondes sismiques dans la masse d'eau 106 et une partie des ondes sismiques se propage vers le bas à travers l'eau vers les formations sous-terraines 102 et 104 situées en dessous de la masse d'eau 106. Lorsque les ondes sismiques atteignent un réflecteur sismique, une partie des ondes sismiques est réfléchie vers le haut et une partie des ondes sismiques poursuit sa propagation vers le bas. Le réflecteur sismique peut être le fond de l'eau 112 ou une des interfaces entre la formation sousterraine, telle que l'interface 114 entre les formations 102 et 104. Lorsque les ondes réfléchies qui se propagent vers le haut atteignent l'interface eau/air à la surface de l'eau 116, une partie majoritaire des ondes est de nouveau réfléchie vers le bas. De cette manière, les ondes sismiques peuvent être réfléchies plusieurs fois entre des réflecteurs vers le haut, tels que le fond de l'eau 112 ou l'interface avec les formations 114, et le réflecteur vers le bas à la surface de l'eau 116 au-dessus. Chaque fois que les ondes réfléchies se propagent au niveau de la position d'un récepteur sismique 110, le récepteur 110 capte les ondes réfléchies et génère des signaux sismiques représentatifs. Ces signaux sismiques peuvent alors être utilisés pour fournir de précieuses informations concernant les caractéristiques géophysiques des formations sous-terraines explorées. La figure 2 présente une méthode 200 de correction des effets dissipatifs de la terre affectant un enregistrement de traces sismiques selon un ou plusieurs modes de réalisation de l'invention. A l'étape 210, l'enregistrement des traces sismiques présentant une pluralité de décalages est reçu. Dans un mode de réalisation, la pluralité de décalages facilite la conversion des traces sismiques en un domaine à angle constant, comme décrit ci-dessous par référence à l'étape 220. Dans un autre mode de réalisation, les traces sismiques enregistrées sont dans le domaine temps-décalage (t-x). Les traces sismiques présentant la pluralité de décalages peuvent être enregistrées et reçues par n'importe quel processus connu de l'homme du métier. A l'étape 220, les traces sismiques enregistrées sont transformées en un domaine dans lequel chaque trace sismique correspond à un angle constant, où on entend par angle constant l'angle entre un rayon et l'axe vertical à la surface, par exemple la surface de l'eau ou de la terre. Bien que divers modes de réalisation décrits ici soient en rapport avec les explorations marines, ces modes de réalisation peuvent également être utilisés dans le cadre des explorations terrestres. De ce fait, les traces sismiques enregistrées peuvent être transformées en un domaine à angle constant, qui peut également être appelé domaine à paramètres de trajet ou domaine à angle émergent. Dans un mode de réalisation, les traces sismiques enregistrées sont transformées en un domaine tau-p. Dans un autre mode de réalisation, le processus de transformation peut impliquer la transformation des traces sismiques en un domaine f-k. La transformation peut être réalisée par une transformée linéaire de Radon ou n'importe quel autre processus connu de l'homme du métier. Une fois que les traces sismiques enregistrées ont été transformées dans le domaine à angle constant, les traces sismiques peuvent être appelées traces p. A l'étape 230, une première trace sismique ayant une première valeur p (ou d'angle constant) est sélectionnée. A l'étape 240, un modèle de vitesse en fonction de to et un modèle Q en fonction de to sont reçus, où to peut représenter le temps vertical ou le temps à décalage nul. A l'étape 250, le modèle Q est converti en une fonction de tau pour la valeur p (ou l'angle constant) sélectionnée. Dans un mode de réalisation, cette conversion ne peut être effectuée qu'une seule fois si le modèle Q et le modèle de vitesse ne varient pas dans l'espace. Selon divers modes de réalisation de l'invention, la méthode 200 décrit le traitement d'un enregistrement avant pile unique (ou rassemblement), qui peut se trouver à un emplacement fixe de l'espace. Toutefois, le Q en fonction de tau calculé pour un enregistrement précédent peut être réutilisé pour le traitement d'enregistrements ultérieurs pour différents emplacements de l'espace. 2890178 7 A l'étape 260, un algorithme Q inverse à décalage nul (ou filtre) est appliqué à la trace sélectionnée au moyen du modèle Q en fonction de tau. A l'étape 270, on détermine si une autre trace sismique, par exemple une seconde trace sismique ayant une seconde valeur p (ou d'angle constant), issue de l'enregistrement des traces sismiques dans le domaine à angle constant doit être traitée. Si la réponse est affirmative, le procédé revient à l'étape 250 où le modèle Q est converti en une fonction de tau pour la seconde valeur de p (ou d'angle constant). Les étapes 250-270 se poursuivent jusqu'à ce que chaque trace sismique ayant cette valeur respective de p (ou d'angle constant) issue de l'enregistrement des traces sismiques dans le domaine à angle constant ait été traitée. Si la réponse est négative, le traitement se poursuit jusqu'à l'étape 280. A l'étape 280, les traces sismiques enregistrées peuvent être transformées en un domaine désiré, tel que le domaine temps-décalage (t-x) . Dans un mode de réalisation, cette transformation peut être effectuée au moyen d'une transformée linéaire inverse de Radon ou de toute autre méthode couramment connue de l'homme du métier. Dans un autre mode de réalisation, les traces sismiques enregistrées peuvent rester dans le domaine tau-p pour traitement ultérieur. La méthode 200 peut ainsi être utilisée pour éliminer (ou compenser) les effets Q ou l'atténuation causée par la terre affectant l'enregistrement des traces sismiques. Divers modes de réalisation de l'invention montrent que lorsqu'une hypothèse de modèle unidimensionnel est appropriée, c'està-dire un modèle constitué de couches sensiblement horizontales, des filtres Q inverses peuvent être appliqués aux traces sismiques avant pile dans le domaine tau-p pour compensation précise. De plus, tout algorithme Q inverse variable dans le temps, et précis pour des traces sismiques à décalage nul, peut être utilisé pour effectuer l'application de filtre décrite ici. Dans un mode de réalisation, la transformée tau-p décrite cidessus en référence à l'étape 220 peut être combinée à un autre filtrage dans le domaine tau-p (par exemple, pour prendre en compte le bruit ou les atténuations multiples) pour éviter les coûts liés à la réalisation de la transformée tau-p. La figure 3 présente une méthode 300 pour corriger les effets dissipatifs de la terre d'un enregistrement de traces sismiques conformément à un ou plusieurs modes de réalisation de l'invention. Les étapes 310-330 sont identiques aux étapes 210-230, c'est pourquoi les étapes 310-330 sont explicitées en se référant aux étapes 210-230 ci-dessus. A l'étape 340, un modèle de vitesse en fonction de to et un modèle Q en fonction de to sont reçus, où to peut représenter le temps vertical ou le temps à décalage nul. A l'étape 350, la trace sismique sélectionnée est convertie du domaine tau-p en un domaine temps de propagation-p, où le temps de propagation peut également être appelé temps de trajet. Cette conversion peut nécessiter l'étirement de la trace. A l'étape 360, le modèle Q est converti à partir d'une fonction de to en une fonction du temps de propagation, c'est-à-dire Q(t), où t représente le temps de propagation. A l'étape 370, un algorithme (ou filtre) Q inverse à décalage nul est appliqué à la trace sismique sélectionnée en utilisant le modèle Q en fonction du temps de propagation. A l'étape 380, on détermine si une autre trace sismique, par exemple une seconde trace sismique ayant une seconde valeur p constante, issue de l'enregistrement des traces sismiques dans le domaine à angle constant doit être traitée. Si la réponse est affirmative, le traitement revient à l'étape 350. Si la réponse est négative, le traitement se poursuit jusqu'à l'étape 390. A l'étape 390, les traces sismiques enregistrées peuvent être transformées en un domaine désiré, tel que le domaine temps-décalage (t-x) . Dans un mode de réalisation, cette transformation peut être effectuée au moyen d'une transformée linéaire inverse de Radon ou de tout autre méthode connue de l'homme du métier. Dans un autre mode de réalisation, les traces sismiques enregistrées peuvent rester dans le domaine tau-p pour traitement ultérieur. La méthode 300 peut ainsi être utilisée pour éliminer (ou compenser) les effets Q ou l'atténuation causée par la terre affectant l'enregistrement des traces sismiques. Les paragraphes qui suivent décrivent une dérivation mathématique pour la conversion du modèle en fonction du temps vertical to en modèle Q en fonction de tau. Si un trajet à un paramètre de trajet fixe (lenteur) p, l'angle entre le trajet et la verticale peut être défini par e(to) avec p = sin(9)/v, où to représente le temps vertical ou le temps à décalage nul, et non le temps de propagation. Si le temps vertical to est fixe, alors le temps de propagation et le décalage peuvent être exprimés comme: t = fodtl/(9) Équation (1) x = 1 vtan(0)dto' Équation (2) où v et e sont des fonctions de to' dans la fonction à intégrer. Pour un événement à to fixe, dt/dx = p. De ce fait, le temps à décalage nul peut être exprimé comme: T=t-xdt/dx=t px Équation (3). Si on incorpore t de l'équation (1) et x de l'équation (2) dans l'équation (3) on obtient une équation qui lie T et to, où T est tau. Pour un to et un p donnés, l'équation (3) donne le temps à décalage nul, T, de l'événement dans le domaine T-p. Les filtres Q inverses dépendant de t/Q, où t est le temps de propagation le long du trajet, il faut que dT/Q(T, p) = dt/Q(to(T)), où dT est l'évolution de T correspondant à dt, qui est un changement de t. De ce fait, Q(T, p) = dT/dt Q(to(T)). Pour un milieu à couches à vitesse constante, les intégrations dans les équations (1) et (2) peuvent devenir des sommations. Considérons la dème couche, pour laquelle Q = Q; et v = v;. A un p fixe, l'angle e; entre le trajet et la verticale est également connu. Si l'épaisseur verticale de la couche (temps de double trajet) est dto,;, alors le temps le long du trajet peut être représenté par dt; = dto,;/cos(0,), et le décalage horizontal dans la couche peut être représenté par dx; = dto,; v; tan(ei). L'incrément de T pour cette couche peut être représenté par: dT; = dt; dx; p = dto,; cos(6i) = dt; cos2(A;). Cette dérivation peut alors fournir un algorithme pour construire Q(T,p) couche par couche. Pour une réflexion dans la couche supérieure (i=1) au temps T dans le domaine T-p, le temps de propagation effectif le long du trajet est donné par t = T/cos2(01). De ce fait, Q(T,p) peut être représenté sous la forme: 15 Q(T, p) = T/t Q(to(T)) = cos2(e1) QI. L'épaisseur de la couche peut être exprimée sous la forme dri = dto, 1cos(e1). De même, les autres couches du milieu peuvent être incluses simplement en réduisant les valeurs d'intervalle Q par cos2(e) et les temps à décalage nul par cos(e). De plus, certaines régions du domaine T-p peuvent ne pas être physiques, c'est-à-dire que les rayons de lenteur donnée, p, ne peuvent se propager dans les couches qu'à une vitesse telle que pv La figure 4 présente un réseau informatique 400, dans lequel des modes de réalisation de l'invention peuvent être mis en oeuvre. Le réseau informatique 400 comprend un ordinateur système 430, qui peut être mis en oeuvre comme tout ordinateur personnel ou toute station de travail, comme par exemple une station de travail à base UNIX. L'ordinateur système 430 est en communication avec les disques de stockage 429, 431, et 433, qui peuvent être des disques durs externes. Considérons que les disques de stockage 429, 431, et 433 sont des disques durs classiques et qu'en tant que tels ils seront mis en oeuvre au moyen d'un réseau local ou par accès à distance. Il va de soi que bien que les disques de stockage 429, 431, et 433 soient représentés comme des dispositifs séparés, un disque de stockage unique peut être utilisé pour stocker toutes les instructions programmes, données de mesure et résultats, selon besoin. Dans un mode de réalisation, les données sismiques issues d'hydrophones sont stockées dans le disque 431. L'ordinateur système 430 peut récupérer les données appropriées sur le disque 431 pour appliquer la méthode de correction des traces sismiques en fonction des instructions du programme correspondant aux méthodes décrites ici. Les instructions du programme peuvent être écrites dans un langage de programmation tel que C++, Java et similaire. Les instructions du programme peuvent être stockées dans une mémoire lisible par un ordinateur, comme par exemple un disque de stockage 433. II va de soi que la mémoire dans laquelle sont stockées les instructions du programme peut être de tout type classique utilisé pour le stockage de programmes informatiques, y compris, disques durs, disquettes, CD-ROM et autres supports optiques, bande magnétique, et similaires. Selon le mode de réalisation de l'invention préféré, l'ordinateur système 430 présente les données de sortie principalement sur un écran graphique 427, ou bien via l'imprimante 428. L'ordinateur système 430 peut stocker les résultats des méthodes ci-dessus sur le disque de stockage 429, pour utilisation et analyse ultérieure. Le clavier 426 et le périphérique de pointage (souris, boule roulante ou similaire) 425 peut être fourni avec l'ordinateur système 430 pour permettre une exploitation interactive. L'ordinateur système 430 peut être localisé dans un centre informatique distant de la région d'exploration. L'ordinateur système 430 est en communication avec des hydrophones (soit directement, soit via un équipement d'enregistrement, non représenté), pour recevoir des signaux représentatifs de l'énergie sismique réfléchie. Ces signaux, après formatage classique et autre traitement initial, sont stockés par l'ordinateur système 430 sous forme de données numériques dans le disque de stockage 431 pour récupération et traitement ultérieur de la manière décrite ci-dessus. Bien que la figure 4 illustre le disque de stockage 431 connecté directement à l'ordinateur système 430, on peut également envisager que le disque de stockage 431 soit accessible via un réseau local ou par accès à distance. De plus, bien que les disques de stockage 429, 431 soient représentés comme des dispositifs séparés pour le stockage des données sismiques d'entrée et des résultats d'analyse, les disques de stockage 429, 431 peuvent être mis en oeuvre dans un disque unique (soit ensemble, soit séparément du disque de stockage 433), ou de toute autre manière classique comme le comprendra facilement l'homme du métier se reportant à cette spécification. Bien que ce qui précède soit destiné aux modes de réalisation de la présente invention, d'autres modes de réalisation de l'invention peuvent être imaginés sans s'écarter du domaine d'application de base de celle-ci, qui est déterminé par les revendications qui suivent. REVENDICATIONS 1. Une méthode pour le traitement d'un enregistrement de traces sismiques, comportant: la transformation (220) de l'enregistrement des traces sismiques en un domaine à angle constant tel que chaque trace sismique corresponde à un angle constant, où on entend par angle constant un angle entre un trajet et un axe vertical à la surface; et l'application d'un algorithme 0 inverse à décalage nul à 10 l'enregistrement des traces sismiques dans le domaine à angle constant. 2. Méthode selon la revendication 1, caractérisée en ce qu'elle comporte la transformation (220) de l'enregistrement des traces sismiques en un domaine temps-décalage. 3. Méthode selon la revendication 2, caractérisée en ce que l'on transforme l'enregistrement des traces sismiques dans le domaine tempsdécalage au moyen d'une transformée linéaire inverse de Radon. 4. Méthode selon la revendication 1, caractérisée en ce que l'enregistrement des traces sismiques comprend une pluralité de décalages. 5. Méthode selon la revendication 1, caractérisée en ce que la surface est une surface d'eau. 6. Méthode selon la revendication 1, caractérisée en ce que la transformation de l'enregistrement des traces sismiques dans le domaine à angle constant comporte la transformation de l'enregistrement des traces sismiques en un domaine tau-p. 7. Méthode selon la revendication 1, caractérisée en ce que la transformation de l'enregistrement des traces sismiques dans le domaine à angle constant comporte la transformation de l'enregistrement des traces sismiques en un domaine f-k. 8. Méthode selon la revendication 1, caractérisée en ce que l'on 5 transforme l'enregistrement des traces sismiques dans le domaine à angle constant au moyen de l'opération de transformée linéaire de Radon. 9. Méthode selon la revendication 1, caractérisée en ce que le domaine à angle constant est un domaine tau-p. 10. Méthode selon la revendication 1, caractérisée en ce que le domaine à angle constant est un domaine de paramètres de trajet. 11. Méthode selon la revendication 1, caractérisée en ce que l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comporte la conversion d'un modèle Q pour l'algorithme Q inverse à décalage nul en une fonction de tau. 12. Méthode selon la revendication 11, caractérisée en ce qu'on applique 20 l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques au moyen du modèle Q en fonction de tau. 13. Méthode selon la revendication 11, caractérisée en ce que le modèle Q en fonction de tau est défini par dT/dt Q(to(T)), où dT représente un changement de tau, dt représente un changement du temps de propagation et to représente le temps vertical et T représente tau. 14. Méthode selon la revendication 1, caractérisée en ce que l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comporte la conversion d'un modèle Q pour l'algorithme Q inverse à décalage nul en une fonction de tau pour chaque angle constant correspondant à chaque trace sismique. 15. Méthode selon la revendication 1, caractérisée en ce que l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comporte la conversion de l'enregistrement des traces sismiques en un domaine temps de propagation-p. 16. Méthode selon la revendication 15, caractérisée en ce que l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comporte la conversion d'un modèle Q pour l'algorithme Q inverse à décalage nul en une fonction du temps de propagation. 17. Méthode selon la revendication 16, caractérisée en ce que l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques comporte l'application de l'algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques au moyen du modèle Q en fonction du temps de propagation. 18. Méthode pour le traitement d'un enregistrement de traces sismiques, comportant: - la réception de l'enregistrement des traces sismiques dans un domaine temps-décalage; - la transformation de l'enregistrement des traces sismiques en un domaine tau-p; - l'application d'un algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques dans le domaine tau-p; et - la re-transformation de l'enregistrement des traces sismiques dans le domaine temps-décalage. 19. Système informatique, comportant: - un processeur (430) ; et -une mémoire (433) contenant les instructions programme exécutables par le processeur adapté à : É transformer un enregistrement de traces sismiques en un domaine à angle constant tel que chaque trace sismique corresponde à un angle constant, où on entend par angle constant un angle entre un trajet et un axe vertical à la surface; et É appliquer un algorithme Q inverse à décalage nul à l'enregistrement des traces sismiques dans le domaine à angle constant. 20. Système informatique selon la revendication 19, caractérisée en ce que la mémoire comporte les instructions programme exécutables par le processeur adapté à transformer l'enregistrement des traces sismiques en un domaine temps-décalage. 21. Système informatique selon la revendication 19, caractérisée en ce 15 que l'enregistrement des traces sismiques comporte une pluralité de décalages. 22. Système informatique selon la revendication 19, caractérisée en ce que le domaine à angle constant est un domaine tau-p. 23. Système informatique selon la revendication 19, caractérisée en ce que les instructions programme exécutables par le processeur adapté à transformer l'enregistrement des traces sismiques dans le domaine tempsdécalage comporte des instructions programme exécutables par le processeur pour convertir un modèle Q utilisable avec l'algorithme Q inverse à décalage nul en une fonction de tau. 24. Système informatique selon la revendication 23, caractérisée en ce que l'algorithme Q inverse à décalage nul est appliqué à l'enregistrement des 30 traces sismiques au moyen du modèle Q en fonction de tau. 2890178 17 25. Système informatique selon la revendication 19, caractérisée en ce que la surface est une surface d'eau.