FR 2984525 A1 20130621 FR 1161720 A 20111215 Séparation de champs d'onde pour des enregistreurs sismiques répartis sur des surfaces d'enregistrement non planes Arrière-plan Domaine Technique Les modes de réalisation de l'objet décrit dans le présent document concernent généralement des procédés et des systèmes et, plus particulièrement, des dispositifs et des techniques pour séparer des champs d'onde ascendants et descendants à partir de données sismiques enregistrées sous l'eau ou sous la surface de la terre par un récepteur sismique. Contexte de l'invention Le forage en mer ou sur terre est une procédure coûteuse. Par conséquent, les personnes engagées dans une telle entreprise coûteuse effectuent essentiellement des études géophysiques pour décider plus précisément de l'endroit à forer afin d'éviter un puits sans hydrocarbures ou avec des quantités non commercialisables d'hydrocarbures. L'acquisition et le traitement de données sismiques terrestres et marines génèrent une image de la structure géophysique (appelée sous-surface). Même si cette image/profil ne fournit pas un emplacement précis de réservoirs de pétrole et de gaz, elle suggère, à l'homme du métier, la présence ou l'absence de réservoirs de pétrole et/ou de gaz. Par conséquent, la fourniture d'une image à haute résolution de la sous-surface est une procédure toujours en cours pour explorer les ressources naturelles, y compris entre autres, du pétrole et/ou du gaz. Par exemple, les systèmes marins destinés à enregistrer les ondes sismiques sont basés sur des bandes continues remorquées ou sur des câbles ou des noeuds déployés sur le fond marin. Un exemple d'un système marin traditionnel pour enregistrer des ondes sismiques au niveau du fond marin est représenté sur la figure 1 et ce système est décrit dans le brevet européen numéro EP 1 217 390 dont le contenu est incorporé entièrement dans le présent document par cette référence. Dans ce document, plusieurs récepteurs sismiques 10 sont reliés de manière amovible à un piédestal 12 conjointement avec une mémoire 14. Plusieurs de ces récepteurs sont déployés au fond 16 de l'océan. Un navire source 18 remorque une source sismique 20 qui est configurée pour émettre des ondes sismiques 22 et 24. Les ondes sismiques 22 se propagent vers le bas, vers le fond de l'océan 16. Après avoir été réfléchie à partir d'une structure 26, l'onde sismique (principale) est enregistrée (en tant que trace) par le récepteur sismique 10 tandis que les ondes sismiques 24 réfléchies à la surface de l'eau 28 sont détectées par les récepteurs 10 plus tard. Etant donné que l'interface entre l'eau et l'air est bien considérée comme étant un réflecteur quasi parfait (c'est-à-dire que la surface de l'eau agit comme un miroir pour les ondes sonores ou sismiques), l'onde réfléchie 24 retourne vers le récepteur 10. Cette onde réfléchie est typiquement appelée onde fantôme parce que cette onde est due à une fausse réflexion. Les ondes fantômes sont également enregistrées par les récepteurs 10 mais avec une polarisation différente et un décalage temporel par rapport à l'onde principale 22. Au fur et à mesure que l'onde principale 22 se déplace suivant une direction ascendante vers le récepteur 10, cette onde est parfois appelée champ d'onde ascendant et au fur et à mesure que l'onde fantôme 24 se déplace suivant une direction descendante vers le récepteur 10, cette onde est parfois appelée champ d'onde descendant. La figure 1 représente également le récepteur 10 configuré de façon à pouvoir se détacher du piédestal 12 et s'élever à la surface de l'eau 28 pour être récupéré par un bateau de collecte 30. Sur la base des données récupérées par le récepteur 10, une image de la sous-surface est générée par des analyses plus approfondies. Comme discuté précédemment, chaque onde sismique marine qui arrive au niveau du récepteur 10 est accompagnée d'une réflexion fantôme. Il en va de même pour chaque onde sismique terrestre qui arrive, enregistrée par un récepteur enfoui sous terre. En d'autres termes, les arrivées fantômes suivent leurs arrivées principales et sont générées quand une onde qui se déplace de manière ascendante est enregistrée une première fois sur un équipement sous-marin avant d'être réfléchie au point de contact air/surface. Les signaux principaux et fantômes (fantôme du côté récepteur et non pas le fantôme côté source) sont également appelés de manière typique champs d'onde ascendants et descendants. L'écart temporel entre un évènement et son fantôme dépend entièrement de la profondeur du récepteur 10 et de la vitesse de l'onde dans l'eau (celle-ci peut être mesurée et est considérée comme étant approximativement égale à 1500 m/s) . Elle peut être de quelques millisecondes pour des données sur bandes continues remorquées (profondeurs inférieures à 15 mètres) ou atteindre plusieurs centaines de millisecondes pour des acquisitions à partir de câbles de fond marin (OBC) ou de noeuds de fond marin (OBN) en profondeur. L'effet dégénératif de l'arrivée fantôme sur la largeur et la résolution de la bande sismique est connu. En fait, l'interférence entre les arrivées principales et fantômes provoque des brèches ou des creux dans le contenu de fréquence et ces brèches ne peuvent pas être éliminées sans utiliser une combinaison de techniques de traitement et d'acquisition avancées. Ces techniques de traitement avancées comprennent la séparation de champ d'onde ou la décomposition de champ d'onde ou un défantômage. Ces techniques nécessitent une acquisition de données avancée, à savoir une acquisition marine multicomposant. L'acquisition marine multicomposant utilise des récepteurs pouvant mesurer au moins deux paramètres différents, la pression de l'eau (enregistrée par un hydrophone) et l'accélération ou la vitesse des particules d'eau (enregistrée par un géophone ou un accéléromètre) par exemple. Par conséquent, l'acquisition marine multicomposant fournit, en plus d'un enregistrement de pression P, au moins une composante de vitesse (ou d'accélération) de particule verticale Z. Une étape difficile dans le traitement des données des enregistrements marins multicomposant est la séparation des champs d'onde avant l'étape de pré-empilage. La séparation des champs d'onde permet la séparation des champs d'onde enregistrés en des parties individuelles : des ondes ascendantes et descendantes. Différentes techniques sont connues dans le domaine de la séparation des champs d'onde comme par exemple, Amundsen, 1993, « Wavenumber-based filtering of marine point source data », Geophysics ; ou Ball et Corrigan, 1996, « Dual-sensor summation of noisy ocean-bottom data », SEG Ann.Mtg. ; ou Schalkwijk et al., 2003, « Adaptive decomposition of multicomponent ocean-bottom seismic data into downgoing and upgoing P and S waves », Geophysics, dont l'ensemble des contenus est incorporé dans le présent document par cette référence. Indépendamment du type de séparation et des détails de l'algorithme utilisé, les algorithmes de séparation actuels présument que la surface d'enregistrement est une surface plane. Cependant, le fond marin est une surface d'acquisition non plane. En variante, la profondeur continue remorquée peut varier suivant sa longueur ou des récepteurs enfouis sous terre peuvent être déployés à des profondeurs variables. Par conséquent, pour ce type de situations, la présomption d'une surface plane est fausse et les données récupérées peuvent générer des effets erronés sur l'image finale sauf si elles sont corrigées. Par conséquent, il serait souhaitable de fournir des systèmes et des procédés pouvant éviter les problèmes et les inconvénients susmentionnés et pouvant, par exemple, prendre en considération une surface d'acquisition non plane. Résumé de l'invention Selon un mode de réalisation fourni à titre d'exemple, un procédé pour séparer des champs d'onde ascendants et descendants (U, D) dans des données sismiques associées à une sous-surface d'un corps d'eau ou à une sous-surface d'un corps de roches est proposé. Le procédé comprend une étape consistant à recevoir des données sismiques (P, Z0) enregistrées dans le domaine temps-espace à l'aide d'enregistreurs sismiques répartis sur un premier horizon fondamental, le premier horizon fondamental étant non plan ; une étape consistant à établir une relation mathématique entre des données sismiques transformées (P, Z) et les champs d'onde ascendants et descendants (U, D) sur un second horizon fondamental ; et une étape consistant à résoudre à l'aide d'une procédure d'inversion, exécutée sur un processeur, la relation mathématique afin d'obtenir les champs d'onde ascendants et descendants (U, D) pour le second horizon fondamental. Le second horizon fondamental étant différent du premier. Selon un autre mode de réalisation fourni à titre d'exemple, un dispositif informatique pour séparer des champs d'onde ascendants et descendants (U, D) dans des 10 données sismiques associées à une sous-surface d'un corps d'eau ou à une sous-surface d'un corps de roches est proposé. Le dispositif informatique comprend une interface configurée pour recevoir des données sismiques (Po, Zo) enregistrées dans le domaine temps- 15 espace à l'aide de récepteurs sismiques répartis sur un premier horizon fondamental, le premier horizon fondamental étant non plan ; et un processeur relié à l'interface. Le processeur est configuré pour recevoir une relation mathématique entre des données sismiques 20 transformées (P, Z) et les champs d'onde ascendants et descendants (U, D) sur un second horizon fondamental et pour résoudre, à l'aide d'une procédure d'inversion, la relation mathématique afin d'obtenir les champs ascendants et descendants (U, D) pour le second horizon 25 fondamental. Le second horizon fondamental est différent du premier horizon fondamental. Selon un autre mode de réalisation en exemple, un support pouvant être lu par un ordinateur est proposé, comprenant des instructions exécutables, les 30 instructions implémentant, quand elles sont exécutées par un processeur, des instructions destinées à séparer des champs d'onde ascendants et descendants (U, D) dans des données sismiques associées à une sous-surface d'un corps d'eau ou à une sous-surface d'un corps de roches. 35 Les instructions comprennent la réception de données sismiques (P, Z0) enregistrées dans le domaine temps-espace à l'aide de récepteurs sismiques répartis sur un premier horizon fondamental, le premier horizon fondamental n'étant pas plan ; l'établissement d'une relation mathématique entre des données sismiques transformées (P, Z) et les champs d'onde ascendants et descendants (U, D) sur un second horizon fondamental plan ; et la résolution, à l'aide d'une procédure inverse, exécutée sur un processeur, de la relation mathématique pour obtenir les champs d'onde ascendants et descendants (U, D) pour le second horizon fondamental. Le second horizon fondamental est différent du premier horizon fondamental. Brève description des dessins Les dessins joints, qui sont incorporés dans le document et qui constituent une partie du mémoire, représentent un ou plusieurs modes de réalisation et expliquent, conjointement avec la description, ces 20 modes de réalisation. Sur les dessins : la figure 1 est un schéma représentant un système d'acquisition de données sismiques conventionnel comprenant plusieurs récepteurs sismiques fournis dans le fond marin ; 25 la figure 2 est un schéma illustrant plusieurs récepteurs sismiques fournis sur un horizon fondamental non plan et un horizon fondamental plan au niveau desquels des champs d'onde ascendants et descendants sont calculés selon un mode de réalisation en exemple ; 30 la figure 3 est un organigramme représentant un procédé pour séparer des champs d'onde ascendants et descendants selon un mode de réalisation en exemple ; les figures 4 et 5 sont des graphiques illustrant des composants synthétiques P et Z selon un mode de 35 réalisation en exemple ; les figures 6 et 7 sont des graphiques illustrant des champs d'onde ascendants et descendants séparés selon un mode de réalisation en exemple ; la figure 8 représente l'organigramme d'un autre 5 procédé pour séparer des champs d'onde ascendants et descendants selon un mode de réalisation en exemple ; et la figure 9 est un schéma représentant un appareil configuré pour exécuter un procédé de séparation selon 10 un mode de réalisation en exemple. Description détaillée de l'invention La description suivante des modes de réalisation fournis à titre d'exemple se réfère aux dessins joints. 15 Les numéros de référence identiques sur les différents dessins indiquent les mêmes éléments ou des éléments similaires. La description détaillée suivante ne limite pas l'invention, au contraire, la portée de l'invention est définie par les revendications fournies en annexe. 20 Les modes de réalisation suivants sont discutés, pour plus de simplicité, par rapport à la terminologie et à la structure des algorithmes de défantômage destinés à séparer des champs d'onde ascendants et descendants qui sont enregistrés par plusieurs récepteurs sismiques 25 fournis sur le fond marin à des profondeurs différentes par rapport à la surface de l'eau. Cependant, les modes de réalisation discutés par la suite ne sont pas limités à des récepteurs placés dans le fond marin mais peuvent également être appliqués à des bandes continues 30 dont le récepteur est placé à des profondeurs différentes ou à des récepteurs qui sont enterrés sous la surface de la terre à des profondeurs différentes. Une référence tout au long du mémoire à « un mode de réalisation » signifie qu'une caractéristique, une 35 structure, ou un trait particulier décrit en relation avec un mode de réalisation est compris dans au moins un mode de réalisation de la présente invention. Par conséquent, la présence des phrases « dans un mode de réalisation » à différents endroits du texte tout au long du mémoire ne se réfère pas forcément au même mode de réalisation. En outre, les caractéristiques, structures ou traits particuliers peuvent être combinés de n'importe quelle façon appropriée dans un ou dans plusieurs modes de réalisation. Selon un mode de réalisation en exemple, un nouveau procédé pour séparer des composants ascendants et descendants est fourni, comprenant une première étape (i) consistant à déterminer des équations mettant en rapport les résultats de séparation souhaités sur un horizon fondamental plan aux enregistrements multicomposant disponibles. Ces équations peuvent être formulées dans le domaine f-k (nombre d'ondes par fréquence), dans le domaine tau-p ou dans d'autres domaines équivalents ou dans une combinaison de ces domaines. Ces équations impliquent des termes d'extrapolation de champ d'onde (qui seront discutés plus loin dans le texte) et constituent une fonction des propriétés de support correspondantes (par ex. la vitesse du bruit dans l'eau ou dans les couches de roches). Le nouveau procédé comprend en outre une étape (ii) consistant à inverser les équations de l'étape (i) pour trouver les résultats de séparation souhaités comme une fonction des enregistrements disponibles. 30 Cette étape d'inversion peut être effectuées en utilisant une variété d'algorithmes, par analyse par exemple, ou grâce à la méthode d'inversion des moindres carrés. Il est noté que la quantité des données sismiques qui est utilisée avec les équations et la 35 procédure d'inversion nécessitent qu'un logiciel informatique spécialisé soit implémenté sur un dispositif informatique. Etant donné que le volume des données sismiques devant être traitées pour séparer les composants ascendants et descendants est important, il n'est pas pratique, voire impossible, qu'un être humain effectue mentalement tous ces calculs. Un exemple d'une surface d'acquisition non plane et non plate est un système d'acquisition pour fond marin dans lequel les noeuds sont déployés sur un fond marin ayant une topographie irrégulière. La figure 2 représente un fond marin non plan 40, sur lequel ont été répartis plusieurs récepteurs 42 à 50. Un corps d'eau 52 est situé au-dessus des récepteurs et comprend une interface eau/air 54. Une structure 56 est enterrée sous le fond marin 40 et l'on souhaite en obtenir une image grâce au nouveau procédé. Un autre exemple d'une surface d'acquisition non plane est la situation dans laquelle des profondeurs de bandes continues remorquées varient suivant leur longueur ou dans laquelle des récepteurs sont enfouis sous la surface de la terre à des profondeurs variables. En considérant que chaque récepteur 42 à 50 est configuré pour enregistrer une pression d'eau P et une vitesse de particule Z le long d'un axe z, un exemple possible d'équations (relation mathématique) mettant en relation (i) les ondes ascendantes U et descendantes D au niveau d'un horizon fondamental 60 plan (et également plat) et (ii) les données sismiques enregistrées et transformées P et Z sur un horizon fondamental non plan 62 est donné par les formules : P(x, ' )=I(e+er°`'U(co,k)+ e- ik zàz' D(CO,k)) 1 -e+122,1', k Nk (1) k t 1 +2,d,' Z(x' , )=1 = e-E'kr°2'U(co,k)- e- jk 'Az' D(co,k))-e ' ' k coP Nk (2) . Il est noté que les données d'origine Po et Zo sont une fonction de la position x, et z, de chaque récepteur 42 à 50 ainsi que du temps t où les données sont enregistrées. Pour plus de simplicité, la composante y, est considérée comme étant égale à zéro dans l'exemple représenté sur la figure 2. Cependant, le procédé peut être appliqué pour l'ensemble des trois composantes spatiales de la position des récepteurs. En outre, il est noté que le procédé peut également être appliqué dans une situation où les données sismiques sont enregistrées uniquement avec des récepteurs à composant unique, à savoir le composant Po ou Zo ou un autre composant. Dans cette situation, une équation mettant en rapport le composant enregistré au champ d'onde ascendant sur un horizon fondamental plan peut être trouvé et inversé. Le résultat de l'inversion peut contenir, dans ce cas, du bruit à cause de la présence de brèches fantômes dans le récepteur. Néanmoins, ce bruit peut être réduit ou éliminé efficacement grâce à la procédure d'empilage ou en utilisant des techniques traditionnelles d' atténuation de bruit et d'amélioration de signal. Les données sismiques d'origine Po et Zo sont transformées (à partir du domaine temps-espace), dans cet exemple, à l'aide d'une transformation de Fourier rapide temporelle (FFT) de sorte que le composant temporel t soit désormais un composant de fréquence angulaire co. Par conséquent, les composants P et Z dans les équations (1) et (2) sont la transformée de Fourier rapide de Po et Zo. En outre, il est noté que les composants U et D dans les équations (1) et (2) sont écrits dans le domaine f-k (f étant la fréquence correspondant à la fréquence angulaire co, k étant le nombre d'ondes horizontales et kz étant le nombre d'ondes verticales) et ceux-ci sont les champs d'onde ascendants et descendants devant être calculés. Une autre transformation ou d'autres transformations peuvent être utilisées. La distance verticale Azi dans les équations (1) 5 et (2) est la différence de profondeur entre l'horizon fondamental plan 60 et l'horizon fondamental qui n'est pas plan 62 au niveau du récepteur i. La densité de l'eau est représentée par p, et Nk est un facteur de normalisation concernant le nombre de récepteurs situés 10 au fond de l'eau. Les termes ei2nkxi présents dans les deux équations (1) et (2) concernent une FFT inverse dans l'espace qui transforme les nombres d'ondes k en les coordonnées spatiales des capteurs. Des dispositifs d'extrapolation de champs d'onde 15 pour les champs d'onde ascendants et descendants sont également présents dans les équations (1) et (2). En fait, les dispositifs d'extrapolation de champ d'onde peuvent être trouvés dans les équations mettant en rapport les composants P et Z avec les composants U et 20 D indépendamment du domaine de transformation utilisé. Pour le présent mode de réalisation en exemple, les dispositifs d'extrapolation sont donnés par Les dispositifs d'extrapolation de champ d'onde comprennent des signes opposés pour les composants U et D et ils 25 dépendent du nombre d'ondes verticales et de la différence de profondeur entre l'horizon fondamental plan 60 et l'horizon fondamental non plan 62 au niveau du récepteur i. Dans cet exemple, les dispositifs d'extrapolation de champ d'onde s'appliquent à une 30 propagation acoustique ayant une vitesse constante. L'extension à la situation où la vitesse est variable est bien connue dans le domaine du traitement des données sismiques. Il est noté que l'horizon fondamental plan 60 auquel sont calculés les champs U 35 et D peut être au-dessus ou au-dessous d'un ou de plusieurs des récepteurs 42 à 62. Le mode de réalisation représenté sur la figure 2 représente l'horizon fondamental plan 60 au-dessus des récepteurs. En outre, il est possible que l'horizon fondamental plan 60 soit plat. En outre, il est possible que l'horizon fondamental plan 60 soit au-dessus mais proche des récepteurs 42 à 62. Il est noté que les équations (1) et (2) sont linéaires suivant U et D et par conséquent, les équations peuvent être inversées en utilisant une variété d'algorithmes bien connus. Les détails de ces algorithmes sont omis dans la présente description. Le nouveau procédé discuté ci-dessus peut être implémenté dans un dispositif informatique qui est équipé d'un logiciel dédié pour séparer les composants ascendants et descendants. Le dispositif informatique est discuté ci-après en relation avec la figure 9. Le nouveau procédé est maintenant décrit sur la base de l'organigramme de la figure 3. A l'étape 300, des données sismiques (au moins deux composants qui sont enregistrés, par exemple, P et Z) sont enregistrées à l'aide de capteurs sismiques qui sont fournis au fond de l'océan. Les données sismiques sont transformées dans un premier domaine souhaité à l'étape 302. Par exemple, le premier domaine peut être le domaine espace-fréquence. A l'étape 304, des équations mettant en rapport (1) les données sismiques transformées dans le premier domaine et (2) des champs d'onde ascendants et descendants dans un second domaine sont établies. Le second domaine est différent du premier domaine et peut être, par exemple, le domaine f-k. Les champs d'onde ascendants et descendants correspondent à un horizon fondamental non plan. D'autres domaines sont possibles pour le premier et le second domaine. Les équations sont inversées à l'étape 306 pour trouver les résultats de séparation comme fonction des enregistrements disponibles. Ensuite, après différentes étapes de traitement qui sont connues par l'homme du métier et qui ne seront pas répétées ici, une image de la sous-surface observée est générée à l'étape 308 sur la base des champs U et/ou D séparés. Le procédé mentionné ci-dessus est maintenant appliqué à un ensemble de données synthétiques P et Z. 10 Les données synthétiques P sont représentées sur la figure 4 et les données synthétiques Z sont représentées sur la figure 5. Les données sont générées comme étant enregistrées avec un certain décalage (distance le long de l'axe X) à partir de la source et 15 à un moment t (sur l'axe Y) à partir d'un horizon fondamental d'acquisition non plan. Les données représentées sur les figures 4 et 5 sont calculées, par exemple, par l'intermédiaire d'une modélisation acoustique sur un demi-espace et ne comprennent par 20 conséquent que des arrivées directes. A cet égard, il convient de noter que les arrivées directes représentent le type d'évènement le plus difficile à séparer parce que ses angles de propagation sont normalement plus larges que ceux des autres évènements. 25 Après avoir établi les équations notées à l'étape 304 sur la figure 3, et résolu les composants ascendants et descendants U et D, comme noté à l'étape 306, les composants U et D sont représentés sur les figures 6 et 7, respectivement. La figure 7 montre que le champ 30 d'onde descendant est symétrique alors que la figure 6 montre que le champ d'onde ascendant est complexe en raison de la réflexion au niveau de l'horizon fondamental non plan (c'est-à-dire le fond marin). Un procédé pour séparer des champs d'onde 35 ascendants et descendants (U et D) à partir de données sismiques associées à une sous-surface d'un corps d'eau ou de roches sera discuté maintenant en faisant référence à la figure 8. Le procédé comprend une étape 800 consistant à recevoir des données sismiques (P, Z0) enregistrées dans le domaine temps-espace à l'aide de récepteurs sismiques répartis sur un premier horizon fondamental, le premier horizon fondamental n'étant pas plan ; une étape 802 consistant à établir une relation mathématique entre des données sismiques transformées (P, Z) et les champs d'onde ascendants et descendants (U, D) sur un second horizon fondamental plan ; et une étape 804 consistant à résoudre, à l'aide d'une procédure d'inversion, exécutée sur un processeur, la relation mathématique afin d'obtenir les champs d'onde ascendants et descendants (U, D) pour le second horizon fondamental. Le second horizon fondamental est différent du premier horizon fondamental. Un exemple d'un système informatique représentatif pouvant exécuter les opérations selon les modes de 20 réalisation en exemple discutés ci-dessus est représenté sur la figure 9. Un matériel, un micrologiciel, un logiciel ou une combinaison de ceux-ci peut être utilisé pour exécuter les différentes étapes et opérations décrites dans le présent document. 25 Le système informatique fourni à titre d'exemple 900 approprié pour effectuer les activités décrites dans les modes de réalisation en exemple peut comprendre un serveur 901. Ce serveur 901 peut comprendre une unité centrale (UCT) 902 couplée à une 30 mémoire vive (MEV/RAM) 904 et à une mémoire morte (MEM/ROM) 906. La MEM 906 peut être remplacée par d'autres types de support de stockage pour stocker des programmes, tels qu'une MEM programmable (PROM), une MEMP effaçable (EPROM), etc. Le processeur 902 peut 35 communiquer avec d'autres composants internes et externes à travers un ensemble de circuits 908 d'entrée et de sortie (E/S) et des bus 910 pour fournir des signaux de commande et similaire. Le processeur 902 exécute une variété de fonctions connues par l'homme du métier, comme dicté par les instructions du logiciel et/ou du micrologiciel. Le serveur 901 peut également comprendre un ou plusieurs dispositifs de stockage de données, comprenant des lecteurs de disque dur 912, des lecteurs de CD-ROM 914 et d'autres matériels pouvant lire et/ou sauvegarder des informations, tel qu'un DVD, etc. Dans un mode de réalisation, le logiciel destiné à exécuter les étapes discutées ci-dessus peut être stocké et distribué sur un CD-ROM ou un DVD 916, un support amovible 918 ou une autre forme de support pouvant sauvegarder des informations de façon portable. Ces supports de stockage peuvent être insérés dans des dispositifs, tels que le lecteur de CD-ROM 914, le lecteur 912, etc. ou être lus par ceux-ci. Le serveur 901 peut être couplé à un écran 920, qui peut être n'importe quel type d'affichage ou d'écran de présentation connu, tel que les écrans LCD ou LED, les écrans plasma, à tubes cathodiques (CRT), etc. Une interface d'entrée utilisateur 922 est fournie, comprenant un ou plusieurs mécanismes d'interface utilisateur, comme une souris, un clavier, un microphone, un écran tactile, un clavier tactile, un système de reconnaissance vocale etc. Le serveur 901 peut être couplé à d'autres dispositifs informatiques par l'intermédiaire d'un réseau. Le serveur peut faire partie d'une configuration en réseau plus vaste, tel qu'un réseau global (GAN), comme Internet 928. L'homme du métier appréciera également que les 35 modes de réalisation fournis à titre d'exemple puissent être mis en application dans un dispositif de communication sans fil, dans un réseau de télécommunication, en tant que procédé ou dans un produit de programme informatique. Par conséquent, les modes de réalisation en exemple peuvent prendre la forme d' un mode de réalisation entièrement matériel ou d' un mode de réalisation combinant des aspects matériels et logiciels. En outre, les modes de réalisation en exemple peuvent prendre la forme d' un produit de programme informatique stocké sur un support de stockage pouvant être lu par un ordinateur comprenant des instructions pouvant être lues par un ordinateur mises en application dans le support. N' importe quel support approprié pouvant être lu par un ordinateur peut être utilisé, y compris des disques durs, des CD-ROM, des disques versatiles numériques (DVD) , des dispositifs de stockage optiques ou des dispositifs de stockage magnétiques comme des disquettes ou des bandes magnétiques. D' autres exemples non limitatifs de supports pouvant être lus par un ordinateurs comprennent des mémoires de type flash ou d' autres types de mémoires connus. Les modes de réalisation décrits en exemple fournissent un appareil et un procédé destinés au traitement de données sismiques. On comprendra que cette description n' est pas destinée à limiter l' invention . Au contraire, l'objectif des modes de réalisation en exemple est de couvrir des variantes, des modifications et des équivalences, qui sont comprises dans l' esprit et la portée de l' invention, comme définie par les revendications jointes. En outre, plusieurs détails spécifiques sont soulignés afin de fournir une compréhension cohérente de l' invention revendiquée. Cependant, l' homme du métier comprendra que plusieurs modes de réalisation peuvent être mis en pratique sans ces détails spécifiques. Bien que les caractéristiques et les éléments des présents modes de réalisation en exemple soient décrits dans les modes de réalisation dans des combinaisons particulières, chaque caractéristique ou élément peut être utilisé individuellement sans les autres caractéristiques ou éléments des modes de réalisation ou dans les différentes combinaisons avec ou sans les autres caractéristiques et éléments qui y sont décrits. Cette description écrite utilise des exemples de l'objet décrit afin de permettre à l'homme de métier de le mettre en pratique, y compris de créer et d'utiliser tout dispositif ou système et d'effectuer tout procédé incorporé. La portée pouvant être brevetée de l'objet de la présente invention est définie par les revendications et peut comprendre d'autres exemples évidents à l'homme du métier. Ces exemples sont compris dans la portée des revendications.20 REVENDICATIONS1. Procédé pour séparer des champs d'onde ascendants et descendants (U, D) dans des données sismiques associées à une sous-surface d'un corps d'eau ou à une sous-surface d'un corps de roches, le procédé comprenant les étapes consistant à : recevoir des données sismiques (P, Z0) enregistrées dans le domaine temps-espace avec des capteurs sismiques répartis sur un premier horizon fondamental, le premier horizon fondamental n'étant pas plan ; établir une relation mathématique entre des données sismiques transformées (P, Z) et les champs d'onde ascendants et descendants (U,D) sur un second horizon fondamental plan ; et résoudre à l'aide d'une procédure d'inversion, exécutée sur un processeur, la relation mathématique afin d'obtenir les champs d'onde ascendants et descendants (U, D) pour le second horizon fondamental, le second horizon fondamental étant différent du premier horizon fondamental. 2. Procédé selon la revendication 1, l'étape consistant à établir une relation comprenant en outre l'étape consistant à : transformer les données sismiques (Po, Z0) en des données sismiques transformées (P, Z) dans un premier domaine qui est différent du domaine temps-espace. 3. Procédé selon la revendication 2, l'étape consistant à établir une relation comprenant en outre l'étape consistant à : utiliser des opérateurs d'extrapolation pour 35 multiplier les champs d'onde ascendants et descendants(U, D) dans un second domaine, différent du premier domaine, afin d'établir la relation mathématique. 4. Procédé selon la revendication 3, le second domaine étant un domaine f-k, où f représente une fréquence et k un nombre d'ondes. 5. Procédé selon la revendication 1, les données sismiques étant enregistrées à l'aide de récepteurs 10 sismiques placés sur un fond du corps d'eau. 6. Procédé selon la revendication 1, les données sismiques étant enregistrées à l'aide de récepteurs sismiques placés sur des bandes continues à des 15 profondeurs différentes dans le corps d'eau. 7. Procédé selon la revendication 1, les données sismiques étant enregistrées à l'aide de récepteurs sismiques placés sous la surface de la terre. 20 8. Procédé selon la revendication 1, les données sismiques étant enregistrées à l'aide de récepteurs sismiques multicomposant. 25 9. Procédé selon la revendication 1, le second horizon fondamental étant plat. 10. Dispositif informatique pour séparer des champs d'onde ascendants et descendants (U, D) dans des 30 données sismiques associées à une sous-surface d'un corps d'eau ou à une sous-surface d'un corps de roches, le dispositif informatique comprenant : une interface configurée pour recevoir des données sismiques (P, Zo) enregistrées dans le domaine temps-35 espace à l'aide de récepteurs sismiques répartis sur unpremier horizon fondamental, le premier horizon fondamental n'étant pas plan ; et un processeur relié à l'interface et configuré pour : recevoir une relation mathématique entre des données sismiques transformées (P, Z) et les champs d'onde ascendants et descendants (U, D) sur un second horizon fondamental, et résoudre à l'aide d'une procédure d'inversion la 10 relation mathématique afin d'obtenir les champs d'onde ascendants et descendants (U, D) pour le second horizon fondamental, le second horizon fondamental étant différent du premier horizon fondamental. 15