FR2985831A1 - Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface - Google Patents

Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface Download PDF

Info

Publication number
FR2985831A1
FR2985831A1 FR1350277A FR1350277A FR2985831A1 FR 2985831 A1 FR2985831 A1 FR 2985831A1 FR 1350277 A FR1350277 A FR 1350277A FR 1350277 A FR1350277 A FR 1350277A FR 2985831 A1 FR2985831 A1 FR 2985831A1
Authority
FR
France
Prior art keywords
data
window
internal
windows
multiples
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
FR1350277A
Other languages
English (en)
Inventor
Barry Hung
Min Wang
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sercel SAS
Original Assignee
CGG Services SAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CGG Services SAS filed Critical CGG Services SAS
Publication of FR2985831A1 publication Critical patent/FR2985831A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/56De-ghosting; Reverberation compensation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Control For Baths (AREA)

Abstract

Un système et un procédé sont présentés pour éliminer sensiblement l'influence de multiples internes lors de la cartographie sismique de zones géographiques sous l'eau présentant un intérêt sans connaissance a priori des informations de sous-surface. Le système et le procédé localisent de manière itérative des horizons générant des multiples pour prédire des multiples internes et utilisent une relation inférieur-supérieur-inférieur entre les horizons générant des multiples. Le système et le procédé fournissent des moyens appropriés et rentables pour une atténuation des multiples internes sans informations de sous-surface.

Description

TRAITEMENT DE DONNEES SISMIQUES COMPRENANT UNE MODELISATION DE MULTIPLES INTERNES SANS INFORMATIONS DE SOUS-SURFACE INFORMATIONS DE PRIORITÉ [1] La présente demande comporte une revendication de priorité sous 35 USC § 119(e) fondée sur la demande de brevet US provisoire n° 61/585 830, déposée le 12 janvier 2012, dont le contenu entier est expressément incorporé ici par voie de référence. DOMAINE TECHNIQUE [2] La présente invention concerne généralement l'exploration sismique, et plus spécifiquement des systèmes et des procédés pour éliminer sensiblement l'influence de réflexions de multiples lors du traitement de données sismiques acquises sans connaissance a priori d'informations de sous-surface. CONTEXTE [3] Une technique largement utilisée pour rechercher du pétrole et du gaz consiste en l'exploration sismique de structures géophysiques de sous-surface. Le processus d'exploration sismique consiste en la génération d'ondes sismiques (c'est-à-dire, d'ondes sonores) dirigées vers la zone de sous-surface, le rassemblement des données concernant les réflexions des ondes sismiques générées au niveau des interfaces entre les couches de la sous-surface, et l'analyse des données pour générer un profil (une image) de la structure géophysique, c'est-à-dire, des couches de la sous-surface examinée. Ce type d'exploration sismique peut être utilisé à la fois sur la sous-surface de zones terrestres et pour explorer la sous-surface du fond océanique. 36-122/100306US [4] Il est connu des hommes du métier de l'exploration sismique qu'un choix approprié de fréquences pour commander un dispositif de production de son peut être utilisé pour générer des ondes sismiques dont les réflexions peuvent, à leur tour, être utilisées pour déterminer l'emplacement possible ou probable de gisements d'hydrocarbure, par exemple, sous le fond océanique. Le dispositif de production de son dans ces applications marines peut être appelé vibreur marin, et est généralement également appelé « source », c'est-à-dire une source des ondes sonores qui sont émises et ensuite réfléchies/réfractées hors du fond océanique et ensuite reçues par un ou plusieurs récepteurs, habituellement des douzaines. Des vibreurs marins (appelés ci-après « vibreurs », « vibreurs marins » et/ou « vibreurs sismiques ») peuvent être mis en oeuvre dans ce qui est appelé « ensembles remorqués » de la pluralité de sources et de leurs récepteurs, où chaque ensemble remorqué peut comprendre de nombreux vibreurs, de nombreux récepteurs, et peut comprendre plusieurs groupes de récepteurs ou plus, chacun sur son propre câble, avec une source correspondante, de nouveau sur son propre câble. Des systèmes et des procédés pour leur utilisation ont été produits pour des dispositifs qui peuvent maintenir ces câbles, par exemple, le long de lignes relativement droites alors qu'ils sont remorqués derrière des navires sur l'océan. Comme les hommes du métier peuvent l'apprécier, une industrie entière a été créée pour explorer les océans quant à de nouveaux gisements d'hydrocarbures, et a été appelée « sismologie de réflexion ». [5] Pour un processus de rassemblement sismique, comme montré sur la figure 1, un système d'acquisition de données 10 comprend un navire 2 remorquant plusieurs flûtes 6 qui peuvent s'étendre sur des kilomètres derrière le navire 2. Chacune des flûtes 6 peut comprendre un ou plusieurs oiseaux (« bird » en terminologie anglo-saxonne) 13 qui maintiennent la flûte 6 dans une position fixe connue par rapport aux autres flûtes 6, et les oiseaux 13 sont capables de déplacer la flûte 6 selon les souhaits en fonction de communications bidirectionnelles que les oiseaux 13 peuvent recevoir du navire 2. Un ou plusieurs réseaux de sources 4a,b peuvent -2- 36-122/100306US également être remorqués par le navire 2 ou un autre navire pour générer des ondes sismiques. Les réseaux de sources 4a,b peuvent être placés soit devant, soit derrière les récepteurs 14, ou à la fois derrière et devant les récepteurs 14. Les ondes sismiques générées par les réseaux de sources 4a,b se propagent vers le bas, sont réfléchies par les fonds océaniques et pénètrent dans ceux-ci, dans lesquels les ondes réfractées sont finalement réfléchies par une ou plusieurs structures de réflexion (non montrées sur la figure 1) vers la surface (voir la figure 2, examinée ci-dessous). Les ondes sismiques réfléchies se propagent vers le haut et sont détectées par des récepteurs 14 prévus sur les flûtes 6. Ce processus est généralement appelé « tir » vers une zone des fonds océaniques particulière, et la zone des fonds océaniques peut être appelée « cellule ». [006] La figure 2 illustre une vue latérale du système d'acquisition de données 10 de la figure 1. Le navire 2, situé sur la surface de l'océan 46, remorque une ou plusieurs flûtes 6, qui sont composées de câbles 12, et d'une pluralité de récepteurs 14. Deux flûtes de source sont montrées sur la figure 2, lesquelles comprennent des sources 4a,b attachées à des câbles 12a,b respectifs. Chaque source 4a,b est capable d'émettre une onde sonore respective, ou un signal émis 20a,b. A des fins de simplification des dessins, mais bien que cela n'enlève rien à la compréhension des principes présentés, uniquement un premier signal émis 20a sera montré (même si certaines ou la totalité des sources 4 peuvent émettre simultanément (ou non) des signaux émis 20 similaires). Le premier signal émis 20a se propage à travers l'océan 40 et arrive à un premier point de réfraction/réflexion 22a. Le premier signal réfléchi 24a du premier signal émis 20a se propage vers le haut à partir du fond océanique 42, de retour vers les récepteurs 14. Comme les hommes du métier peuvent l'apprécier, à chaque fois qu'un signal - optique ou acoustique - se propage d'un milieu avec un premier indice de réfraction n1 et rencontre un milieu différent, avec un deuxième indice de réfraction n2, une partie du signal émis est réfléchie selon un angle égal à l'angle d'incidence (selon la loi de Snell bien connue), et une deuxième partie du signal émis peut être réfractée (de nouveau selon la loi de Snell). -3- 36-122/100306US [7] Ainsi, comme montré sur la figure 2, le premier signal émis 20a génère un premier signal réfléchi 24a et un premier signal réfracté 26a. Le premier signal réfracté 26a se propage à travers une couche de sédiment 16 (qui peut être appelée génériquement première couche de sous-surface 16) au-dessous du fond océanique 42, et peut maintenant être considéré comme étant un « nouveau » signal émis, de sorte que, lorsqu'il rencontre un deuxième milieu au niveau d'un deuxième point de réfraction/réflexion 28a, un deuxième ensemble de signaux réfractés et réfléchis 32a et 30a, sont générés par la suite. En outre, comme montré sur la figure 2, il s'avère qu'il y a là un gisement d'hydrocarbure 44 important dans un troisième milieu, ou une couche de terre/roche solide 18 (qui peut être appelée génériquement deuxième couche de sous-surface 18). En conséquence, des signaux réfractés et réfléchis sont générés par le gisement d'hydrocarbure, et c'est l'objet du système d'acquisition de données 10 de générer des données qui peuvent être utilisées pour découvrir ces gisements d'hydrocarbure 44. [8] Les signaux enregistrés par les récepteurs sismiques 14 varient dans le temps, en ayant des pics d'énergie qui peuvent correspondre à des réflecteurs entre les couches. En réalité, étant donné que le fond océanique et l'air/eau sont très réfléchissants, certains des pics correspondent à des réflexions de multiples ou à des réflexions parasites qui devraient être éliminées avant que l'image de la structure géophysique puisse être correctement formée. Les ondes principales subissent une seule réflexion d'une interface entre des couches de la sous-surface (par exemple, le premier signal réfléchi 24a). Les ondes autres que les ondes principales sont connues en tant que multiples. Le signal 50a montré sur la figure 2 est un tel exemple d'un multiple, mais comme montré sur la figure 3, il existe d'autres manières de générer des multiples. [9] Comme illustré sur la figure 3, la source sismique 4 produit la première onde émise 20a qui se divise en une onde émise principale 26a (également appelée premier signal réfracté) pénétrant à l'intérieur d'une première couche de sous-surface 16 (également appelée « couche de sédiment » bien que ce ne soit pas nécessairement le cas) sous le fond océanique 42, et un signal -4- 36-122/100306US de multiples associé à une surface 50 qui se propage de retour vers la surface de l'océan 46 (ou quatrième interface). L'onde émise principale 26a est réfléchie une fois au niveau de la deuxième interface 48 entre différentes couches dans une première couche de sous-surface 16 et se propage de retour vers le récepteur 14 en tant que deuxième signal réfléchi 30a. Le signal de multiples associé à une surface 50 atteint également le récepteur 14, mais à un instant différent, après avoir été réfléchi (au moins) deux fois de plus : une première réflexion au niveau de la surface 46 et une deuxième réflexion au niveau du fond océanique 42. Ainsi, le récepteur 14 recevra au moins deux signaux différents provenant du même événement d'émission : le deuxième signal réfléchi 30a et le signal de multiples associé à une surface 50. Le signal de multiples associé à une surface 50 peut être reçu par le récepteur 14 soit avant le deuxième signal réfléchi 30a, soit en même temps, soit après, en fonction de la distance de propagation du premier signal réfracté 26a dans la première couche de sous-surface 16, et de la profondeur de l'océan 40 au point d'émission et de réflexion/réfraction. Ainsi, le récepteur 14 peut ne pas pouvoir « discerner » la nature vraie de l'environnement de sous-surface du fait du ou des signaux de multiples de surface 50. Comme examiné brièvement ci-dessus, d'autres multiples peuvent également être générés, dont certains peuvent également se propager à travers la sous-surface. Un multiple, par conséquent, est n'importe quel signal qui n'est pas un signal réfléchi principal. En outre, un type différent de signal de multiple est également montré sur la figure 3, le signal de multiple interne 51, qui subit une réflexion vers le bas par une couche souterraine, comme montré sur la figure 3. Les « multiples », comme cela est connu des hommes du métier, peuvent créer des problèmes lors de la détermination de la nature vraie de la géologie de la terre au-dessous du fond océanique. Des multiples (que ce soit des multiples associés à une surface 50 ou des multiples internes 51) peuvent être confondus par le système d'acquisition de données 10 avec des premier, deuxième ou troisième signaux réfléchis. Les multiples n'ajoutent pas d'informations utiles concernant la géologie au-dessous du fond océanique, et ainsi ils sont, par essence, du bruit, et il est souhaitable de les éliminer et/ou de sensiblement réduire et/ou éliminer leur influence dans un -5- 36-122/100306US traitement de signal des autres signaux réfléchis de manière à établir correctement la présence (ou l'absence) de gisements d'hydrocarbure souterrains/ sous l'eau. [0010] Alors que les multiples de surface 50 créent certains problèmes lors du traitement de signal pour déterminer la « vraie » nature de la géologie de sous-surface sous l'eau, les multiples internes 51 sont connus comme étant particulièrement problématiques. Les multiples internes 51 apparaissent généralement du fait d'une série de contrastes d'impédances de sous-surface. Ils sont généralement observés dans des données sismiques acquises à divers emplacements, tels que Santos Basin au Brésil. Ils sont souvent faiblement distingués des principaux (c'est-à-dire des premier, deuxième et troisième signaux réfléchis, entre autres), parce qu'ils ont un déplacement, des inclinaisons et une largeur de bande de fréquence similaires, faisant de ce fait d'une atténuation et/ou d'une élimination de multiples internes 51, un des problèmes clés pour fournir des images sismiques nettes en interprétant des zones présentant un intérêt. Avec le temps, divers procédés ont été développés pour aborder ce problème difficile et la plus grande partie de ceux-ci reposent sur la capacité d'identifier les générateurs de multiples. Des approches, telles que le modèle à rétroaction de Delft (Verschuur, D. J. et d'autres, 1996, « Removal of Inter-bed Multiples », 58th Meeting, EAGE, Expanded Abstracts, Paper B003, dont le contenu entier est incorporé ici par voie de référence), le procédé de convolution-corrélation de Jakubowicz (Pica, A. et d'autres, 2008, « Wave Equation Based Internai Multiple Modelling in 3D », 78th Meeting, SEG, Expanded Abstracts, pages 2476 à 2480, dont le contenu entier est incorporé ici par voie de référence), des procédés commandés par un modèle, et une déconvolution prédictive (Hugonnet, P. et d'autres, 2005, « 2D Deconvolution for OBC Data and for Internai Multiple Attenuation Part 1 : Theory », 67th Meeting, EAGE, Extended Abstracts, Paper A026, dont le contenu entier est incorporé ici par voie de référence) nécessitent des informations a priori concernant la sous-surface. Lorsque les informations sont disponibles, peut-être par l'intermédiaire d'enregistrements de puits ou de résultats d'interprétation, une suppression importante des -6- 36-122/100306US multiples internes peut être observée dans ces procédés. Néanmoins, dans de nombreuses situations, il n'est pas souvent facile d'identifier les générateurs de multiples, et cela fait du problème un défi. [0011] D'autres procédés pour gérer des multiples ont également été développés qui ne nécessitent pas des informations a priori de réflecteur de sous-surface. Par exemple, une méthodologie a été développée utilisant une série de dispersion inverse (ISS « inverse scattering serie » en terminologie anglo-saxonne) (Otnes, E. et d'autres, 2004, dans « Attenuation of Internai Multiples for Multicomponent and Towed Streamer Data », 74th Meeting, SEG, Extended Abstracts, pages 1297 à 1300, dont le contenu entier est incorporé ici par voie de référence), et a été appliqué à des données marines ou terrestres pour une atténuation de multiples internes. Le procédé ISS est une approche commandée par des données qui peut prédire tous les multiples internes d'un ordre donné sans informations de sous-surface. Comme les hommes du métier peuvent l'apprécier, l'ordre de multiples fait référence au nombre de rebonds vers le bas qu'un champ d'onde subit avant d'être capturé par un récepteur. Par comparaison, le modèle à rétroaction de Delft ou le procédé de Jakubowicz retire tous les ordres de multiples internes pour une interface donnée (voir Verschuur, D. J. et d'autres, 1999, « A Comparison of the Feedback and Inverse Scattering Internat Multiple Attenuation Methods », 61st Meeting, EAGE, Extended Abstracts, pages 1 à 14, dont le contenu entier est incorporé ici par voie de référence). Il existe des différences significatives entre les deux catégories de méthodologie de modélisation. Une différence consiste en ce que les multiples internes sont catalogués différemment et l'autre différence significative est le besoin des informations a priori dans la première catégorie des technologies de modélisation (voir Matson, K. H. et d'autres, 1998, « Comparing the Interface and Point Scatterer Methods for Attenuating Internai Multiples : A Study with Synthetic Data - Part II », 68th Meeting, SEG, Extended Abstracts, pages 1523 à 1526, dont le contenu entier est incorporé ici par voie de référence). Cependant, il a été affirmé -7- 36-122/100306US que, pour une atténuation de multiples associés à une surface, les deux procédés sont les mêmes en théorie (voir Levin, S. A., 2008, « Delft Inverse Scattering Surface-Related Multiple Attenuation in Three Lines », 78th Meeting, SEG, Extended Abstracts, pages 2512 à 2515, dont le contenu entier est incorporé ici par voie de référence). [0012] Une des spécifications importantes en ISS qui lui permet de prédire des multiples internes sans informations de sous-surface est la spécification pour une condition de monotonicité de pseudo profondeur (voir Nita, B. G. et d'autres, 2007, « Inverse Scattering Internai Multiple Attenuation Algorithm : An Analysis of the Pseudo Depth and Time Monotonicity Requirements », 77th Meeting, SEG, Expanded Abstracts, pages 2461 à 2464, dont le contenu entier est incorporé ici par voie de référence). La condition de monotonicité de pseudo profondeur signifie en fait que, pour un événement de multiples internes particulier, le diffuseur de points qui provoque la réflexion vers le bas est à une profondeur plus haute (en pseudo profondeur) que celle des diffuseurs de points qui provoquent les réflexions vers le haut, c'est-à-dire, satisfaisant une relation « inférieur-supérieur-inférieur » . [0013] La limitation inhérente de l'approche de l'ISS de monotonicité de pseudo profondeur limite sa capacité à être largement utilisée. Par conséquent, il serait souhaitable de proposer des procédés, des modes et des systèmes pour éliminer réellement et efficacement l'influence de multiples internes lors de la détermination d'une géologie du fond sous-océanique afin de la rendre plus facile pour déterminer la présence (ou non) d'hydrocarbures de sous-surface. RÉSUMÉ [0014] Divers modes de réalisation décrits ici résolvent sensiblement au moins un ou plusieurs des problèmes et/ou inconvénients examinés ci-dessus, et fournissent au moins un ou plusieurs -8- 36-122/100306US des avantages décrits ci-dessous. 100151 C'est par conséquent un aspect général des modes de réalisation de fournir un système et un procédé pour prédire des multiples internes sans avoir à identifier des générateurs de multiples qui préviendront ou réduiront à un minimum les problèmes du type précédemment décrit. 100161 Selon un premier aspect des modes de réalisation, un procédé est proposé pour éliminer sensiblement une influence des réflexions de multiples internes en déterminant une géographie sous la mer dans une zone géographique présentant un intérêt sans une connaissance a priori d'informations de sous-surface, le procédé comprenant (a) la génération d'une série de signaux sismiques par une pluralité d'émetteurs de source, (b) la réception de données brutes au niveau d'une pluralité de récepteurs sur la base de la série générée de signaux sismiques, (c) la création d'un ensemble de M fenêtres qui correspondent physiquement à un espace au-dessous de la pluralité de récepteurs et qui comprennent une zone géographique présentant un intérêt, (d) l'attribution des données brutes reçues aux fenêtres respectives de l'ensemble de fenêtres sur la base d'un temps de réception des données brutes, pour générer M trames de données de fenêtre, dans lequel une première trame de données de fenêtre la plus haute incorpore des données brutes qui correspondent à un champ d'onde le plus haut le plus proche de la pluralité de récepteurs, et une Mième trame de données de fenêtre incorpore des données brutes qui correspondent à un champ d'onde le plus bas le plus éloigné de la pluralité de récepteurs, (e) la génération de manière itérative d'un modèle de multiples internes, en utilisant un ensemble glissant de trois trames de données de fenêtre, dans lequel, pour chaque itération, le modèle de multiples internes comprend un premier produit de convolution de données provenant d'une première trame de données de fenêtre et d'une deuxième trame de données de fenêtre, et une corrélation de données provenant d'une troisième trame de données de fenêtre avec le premier produit, (f) la sommation de tous les modèles de multiples internes générés de manière itérative pour créer un modèle de -9- 36-122/100306US multiples internes complet, et la poursuite de la sommation jusqu'à ce que toutes les trames de données de fenêtre aient été utilisées, et (g) la soustraction du modèle de multiples internes complet des données brutes pour éliminer sensiblement l'influence de multiples internes lors de la détermination de la géographie de la zone géographique présentant un intérêt. Selon le premier aspect, le procédé comprend en outre le traitement des données brutes pour supprimer des multiples associés à une surface avant la génération du modèle de multiples internes. De plus, selon le premier aspect, l'étape de traitement pour supprimer des multiples associés à une surface comprend l'utilisation d'un procédé d'élimination de multiples associés à une surface pour supprimer les multiples de surface. [0017] Selon le premier aspect, l'étape de détermination de l'ensemble de M fenêtres est basée sur un temps de propagation de la série de signaux sismiques de la pluralité de sources à la pluralité de récepteurs, et en outre dans lequel les M trames temporelles de fenêtre ont sensiblement une durée similaire, et l'étape de génération d'un modèle de multiples internes en utilisant des données de fenêtre provenant d'un premier ensemble de trois trames de données de fenêtre comprend l'application de l'équation suivante : dans lequel dans chaque ensemble défini de manière itérative de trois trames de données de fenêtre, un champ d'onde supérieur généré par des données dans la trame de données de fenêtre la plus haute est défini par P,J, un premier champ d'onde inférieur généré par des données dans la deuxième trame de données de fenêtre est défini par Pwk, et un deuxième champ d'onde inférieur généré par -10- 36-122/100306US des données dans la troisième trame de données de fenêtre est défini par Pm, et en outre dans lequel, Pwi est un champ d'onde côté source qui représente une réflexion vers le bas d'un multiple interne réfléchi depuis la première fenêtre, Pwk est un champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la deuxième fenêtre, et Pwi est un champ d'onde côté récepteur qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la troisième fenêtre. 100181 Selon le premier aspect, chacune des M fenêtres a des composantes de longueur et de profondeur, et dans lequel la composante de longueur est inférieure ou égale à une distance entre une première source et une dernière source, et en outre dans lequel la composante de profondeur est corrélée à un premier nombre d'échantillons qui est corrélé à une première profondeur en distance, et en outre dans lequel des fenêtres adjacentes se superposent d'un deuxième nombre d'échantillons inférieur au premier nombre d'échantillons, qui correspond à une superposition en profondeur définie en tant que deuxième profondeur, et en outre dans lequel la deuxième profondeur est inférieure à la première profondeur, et en outre dans lequel, pour une valeur croissante de M, la profondeur de la fenêtre augmente. 100191 Selon le premier aspect, chacun de la pluralité d'ensembles de fenêtres satisfait à une condition de monotonicité de pseudo profondeur des fenêtres inférieure-supérieure-inférieure, dans lequel 13',.; est une fenêtre supérieure, et Pwk et Pm sont toutes deux des fenêtres inférieures. 100201 Selon un deuxième aspect des modes de réalisation, un procédé pour déterminer des multiples internes dans des données sismiques marines sans nécessiter une connaissance a priori des sous-surfaces est proposé, le procédé comprenant (a) l'obtention de données sismiques dans et concernant une zone géographique présentant un intérêt (GAI « geographical area of interest » en terminologie anglo-saxonne), (b) la suppression des multiples associés à une surface en relation avec la GAI, (c) la séparation des données sismiques obtenues en une pluralité de 36-122/100306US fenêtres sur la base d'un temps de propagation des champs d'onde générés par une ou plusieurs sources de champ d'onde, et reçus par un ou plusieurs récepteurs de champ d'onde, (d) l'exécution d'une convolution entre les données de fenêtre séparées qui sont responsables des réflexions vers le haut de multiples internes pour générer des données de convolution, (e) l'exécution d'une corrélation entre les données de fenêtre séparées qui sont responsables des réflexions vers le bas des multiples internes correspondants, (f) la détermination du produit des données corrélées et des données convoluées pour créer un modèle de multiples internes, et (g) la répétition des étapes (d) à (f) pour une pluralité d'ensembles de la pluralité de fenêtres, déterminant de ce fait un modèle de multiples internes sensiblement complet pour la GAI. [0021] Selon le deuxième aspect, l'étape d'exécution de la convolution et de la corrélation comprend la définition d'un premier ensemble de trois trames de données de fenêtre, l'exécution de la convolution et de la corrélation sur les données contenues dans l'ensemble de trois fenêtres pour déterminer un modèle de multiples internes, la définition de manière itérative de nouveaux ensembles de trois fenêtres avec toutes les trames de données de fenêtre restantes, et l'exécution de la convolution et de la corrélation sur les données provenant de chacun des nouveaux ensembles définis de manière itérative de trois trames de données de fenêtre, pour déterminer un nouveau modèle de multiples internes, qui est ajouté de manière itérative à un modèle de multiples précédent pour déterminer un modèle de multiples internes sensiblement complet, et dans lequel, pour chacun des ensembles définis de manière itérative de trois fenêtres, une spécification de pseudo monotonicité est satisfaite lorsqu'une relation inférieure-supérieureinférieure apparaît dans chacun des ensembles définis de manière itérative de trois trames de données de fenêtre. [0022] Selon le deuxième aspect, le procédé comprend en outre la soustraction du modèle de multiples internes sensiblement complet des données sismiques obtenues pour éliminer sensiblement une influence de multiples internes lors de la détermination de la géographie d'une -12- 36-122/100306US GAI, et comprend encore en outre l'élimination des multiples de surface des données sismiques avant de déterminer un modèle de multiples internes. [0023] Selon le deuxième aspect, l'étape d'exécution de la convolution et de la corrélation sur les ensembles définis de manière itérative de trames de données de fenêtre comprend l'application de l'équation suivante : dans lequel dans chacun des ensembles définis de manière itérative de trois trames de données de fenêtre, un champ d'onde le plus élevé généré par des données dans une trame de données de fenêtre la plus haute est défini en tant que 13,i, un premier champ d'onde inférieur généré par des données dans une deuxième trame de données de fenêtre est défini en tant que Pwk, et un deuxième champ d'onde inférieur généré par des données dans une troisième trame de données de fenêtre est défini en tant que Pwi, et en outre dans lequel, 13,i est un champ d'onde côté source qui représente une réflexion vers le bas d'un multiple interne réfléchi depuis la première fenêtre, Pwk est un champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la deuxième fenêtre, et Pwi est un champ d'onde côté récepteur qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la troisième fenêtre. [0024] Selon un troisième aspect des modes de réalisation, un système pour éliminer sensiblement une influence de réflexions de multiples internes pour la détermination de la géographie sous la mer dans une zone géographique présentant un intérêt sans connaissance a priori d'informations de sous-surface est fourni, le système comprenant (a) une pluralité -13- 36-122/100306US d'émetteurs de source configurés pour générer une série de signaux sismiques, (b) une pluralité de récepteurs configurés pour recevoir des données brutes basées sur la série générée de signaux sismiques, et (c) un processeur configuré pour créer un ensemble de M fenêtres qui correspondent physiquement à un espace au-dessous de la pluralité de récepteurs et qui comprennent une zone géographique présentant un intérêt, attribuer les données brutes reçues aux fenêtres respectives de l'ensemble de fenêtres sur la base d'un temps de réception des données brutes, pour générer M trames de données de fenêtre, dans lequel une première trame de données de fenêtre la plus haute incorpore des données brutes qui correspondent à un champ d'onde le plus haut le plus proche de la pluralité de récepteurs, et une Miènle trame de données de fenêtre incorpore des données brutes qui correspondent à un champ d'onde le plus bas le plus éloigné de la pluralité de récepteurs, générer de manière itérative un modèle de multiples internes, en utilisant un ensemble glissant de trois trames de données de fenêtre, dans lequel, pour chaque itération, le modèle de multiples internes comprend un premier produit de convolution de données provenant d'une première trame de données de fenêtre et d'une deuxième trame de données de fenêtre, et une corrélation de données provenant d'une troisième trame de données de fenêtre avec le premier produit, sommer tous les modèles de multiples internes générés de manière itérative pour créer un modèle de multiples internes complet, et poursuivre la sommation jusqu'à ce que toutes les trames de données de fenêtre aient été utilisées, et soustraire le modèle de multiples internes complet des données brutes pour éliminer sensiblement l'influence de multiples internes lors de la détermination de la géographie de la zone géographique présentant un intérêt. [0025] Selon le troisième aspect, le processeur est en outre configuré pour traiter les données brutes pour supprimer des multiples associés à une surface avant la génération du modèle de multiples internes, et en outre dans lequel le processeur est en outre configuré pour supprimer des multiples associés à une surface en utilisant un procédé d'élimination de multiples associés à -14- 36-122/100306US une surface. [0026] Selon le troisième aspect, le processeur est en outre configuré pour déterminer l'ensemble de M fenêtres sur la base d'un temps de propagation de la série de signaux sismiques de la pluralité de sources à la pluralité de récepteurs, et en outre dans lequel la totalité des M trames temporelles de fenêtre ont sensiblement une durée similaire. 100271 Selon le troisième aspect, le processeur est en outre configuré pour générer un modèle de multiples internes en utilisant des données de fenêtre provenant d'un premier ensemble de trois trames de données de fenêtre qui comprend l'application de l'équation suivante : dans lequel dans chaque ensemble défini de manière itérative de trois trames de données de fenêtre, un champ d'onde supérieur généré par des données dans la trame de données de fenêtre la plus haute est défini en tant que Pwj, un premier champ d'onde inférieur généré par des données dans la deuxième trame de données de fenêtre est défini en tant que Pwk, et un deuxième champ d'onde inférieur généré par des données dans la troisième trame de données de fenêtre est défini en tant que Pwi, et en outre dans lequel Pwj est un champ d'onde côté source qui représente une réflexion vers le bas d'un multiple interne réfléchi depuis la première fenêtre, Pwk est un champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la deuxième fenêtre, et Pwi est un champ d'onde côté récepteur qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la troisième fenêtre. -15- 36-122/100306US [0028] Selon le troisième aspect, chacune des M fenêtres a des composantes de longueur et de profondeur, et dans lequel la composante de longueur est inférieure ou égale à une distance entre une première source et une dernière source, et en outre dans lequel la composante de profondeur est corrélé à un premier nombre d'échantillons qui est corrélé à une première profondeur en distance, et en outre dans lequel des fenêtres adjacentes se superposent d'un deuxième nombre d'échantillons inférieur au premier nombre d'échantillons, qui correspond à une superposition en profondeur définie en tant que deuxième profondeur, et encore en outre dans lequel la deuxième profondeur est inférieure à la première profondeur, et encore en outre dans lequel, pour une valeur croissante de M, la profondeur de la fenêtre augmente. [0029] Selon le troisième aspect, chacun de la pluralité d'ensembles de fenêtres satisfait à une condition de monotonicité de pseudo profondeur des fenêtres inférieure-supérieure-inférieure, dans lequel 13,i est une fenêtre supérieure, et Pwk et Pwi sont toutes les deux des fenêtres inférieures. [0030] Selon un quatrième aspect des modes de réalisation, un système pour déterminer des multiples internes dans des données sismiques marines sans nécessiter une connaissance a priori des sous-surfaces est proposé, le système comprenant un processeur configuré pour obtenir des données sismiques dans et concernant une zone géographique présentant un intérêt (GAI), dans lequel le processeur est en outre configuré pour (a) supprimer des multiples associés à une surface en relation avec la GAI, (b) séparer les données sismiques obtenues en une pluralité de fenêtres sur la base d'un temps de propagation des champs d'onde générés par une ou plusieurs sources de champ d'onde, et reçus par un ou plusieurs récepteurs de champ d'onde, (c) effectuer une convolution entre les données de fenêtre séparées qui sont responsables des réflexions vers le haut de multiples internes pour générer des données de convolution, (d) effectuer une corrélation entre les données de fenêtre séparées qui sont responsables des réflexions vers le bas des multiples internes correspondants, (e) déterminer le produit des données corrélées et des données -16- 36-122/100306US convoluées pour créer un modèle de multiples internes, et répéter les étapes (c) à (e) pour une pluralité d'ensembles de la pluralité de fenêtres, déterminant de ce fait un modèle de multiples internes sensiblement complet pour la GAI. [0031] Selon le quatrième aspect, le processeur est en outre configuré pour effectuer la convolution et la corrélation en définissant un premier ensemble de trois trames de données de fenêtre, appliquer la convolution et la corrélation à des données contenues dans l'ensemble de trois fenêtres pour déterminer un modèle de multiples internes, définir de manière itérative des nouveaux ensembles de trois fenêtres avec toutes les trames de données de fenêtre restantes, et appliquer la convolution et la corrélation aux données provenant de chacun des nouveaux ensembles définis de manière itérative de trois trames de données de fenêtre, pour déterminer un nouveau modèle de multiples internes, qui est ajouté de manière itérative à un modèle de multiples précédent pour déterminer un modèle de multiples internes sensiblement complet, et dans lequel pour chacun des ensembles définis de manière itérative de trois fenêtres, une spécification de pseudo monotonicité est satisfaite lorsqu'une relation inférieure-supérieureinférieure apparaît dans chacun des ensembles définis de manière itérative de trois trames de données de fenêtre. [0032] Selon le quatrième aspect, le processeur est en outre configuré pour soustraire le modèle de multiples internes sensiblement complet des données sismiques obtenues pour éliminer sensiblement une influence de multiples internes lors de la détermination de la géographie d'une GAI, et le processeur est en outre configuré pour éliminer des multiples de surface des données sismiques avant de déterminer un modèle de multiples internes. [0033] Selon le quatrième aspect, le processeur est en outre configuré pour appliquer la convolution et la corrélation aux ensembles définis de manière itérative de trames de données de fenêtre en appliquant l'équation suivante : -17- 36-122/100306US dans lequel dans chacun des ensembles définis de manière itérative de trois trames de données de fenêtre, un champ d'onde le plus élevé généré par des données dans une trame de données de fenêtre la plus haute est défini en tant que Pwi, un premier champ d'onde inférieur généré par des données dans une deuxième trame de données de fenêtre est défini en tant que Pwk, et un deuxième champ d'onde inférieur généré par des données dans une troisième trame de données de fenêtre est défini en tant que P,I, et en outre dans lequel, Pwi est un champ d'onde côté source qui représente une réflexion vers le bas d'un multiple interne réfléchi depuis la première fenêtre, Pwk est un champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la deuxième fenêtre, et Pwi est un champ d'onde côté récepteur qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la troisième fenêtre. BRÈVE DESCRIPTION DES DESSINS [0034] Les caractéristiques ci-dessus du présent concept général de l'invention deviendront évidentes et seront plus facilement appréciées à partir de la description qui suit des modes de réalisation en faisant référence aux figures suivantes, sur lesquelles des numéros de référence identiques font référence à des parties identiques sur les diverses figures sauf spécification contraire, et sur lesquelles : [0035] la figure 1 illustre un système d'acquisition de données pour une utilisation dans un -18- 36-122/100306US processus de rassemblement sismique ; [0036] les figures 2 et 3 illustrent des vues latérales du système d'acquisition de données de la figure 1 ; [0037] la figure 4A illustre un tracé de données d'un processus connu pour prédire l'existence de multiples internes réfléchis vers le bas à partir d'un premier horizon ; [0038] la figure 4B illustre un tracé de données de procédé pour prédire l'existence de tous les multiples internes sans identifier des horizons générant des multiples selon un exemple de mode de réalisation ; [0039] les figures 5A et 5B illustrent des données de site réelles, et des données synthétiques générées à partir des données réelles, acquises respectivement à partir de la découverte de pétrole de Tupi près de Santos Basin, au large du Brésil, où des multiples internes importants sont évidents ; [0040] la figure 6A illustre un premier modèle de tous les multiples internes des champs pétroliers de Tupi prédits en utilisant le procédé à fenêtre glissante selon un exemple de mode de réalisation sans informations de sous-surface (c'est-à-dire, sans informations de sous-surface a priori) ; [0041] la figure 6B illustre un deuxième modèle de tous les multiples internes des champs pétroliers de Tupi prédits en utilisant un procédé à série de dispersion inverse ; [0042] la figure 7A illustre les données de la figure 6A auxquelles une soustraction adaptative a été appliquée pour supprimer des multiples internes selon un exemple de mode de réalisation ; [0043] la figure 7B illustre les données de la figure 6B (générées en utilisant le procédé à série -19- 36-122/100306US de dispersion inverse) auxquelles une soustraction adaptative a été appliquée pour supprimer des multiples internes ; [0044] la figure 8A illustre une section agrandie de la figure 4A ; [0045] la figure 8B illustre la section agrandie de la figure 4A comme sur la figure 8A à la suite d'une mise en oeuvre du procédé selon un exemple de mode de réalisation ; [0046] la figure 8C illustre la différence entre les figures 8A et 8B (c'est-à-dire, les données de la figure 8B soustraites des données de la figure 8A) ; [0047] la figure 9A illustre les mêmes données de la figure 8A à la suite d'un traitement en utilisant le procédé à série de dispersion inverse ; [0048] la figure 9B illustre la différence entre les figures 9A et 8A ; [0049] la figure 10 illustre un organigramme d'un procédé pour éliminer sensiblement l'influence de réflexions de multiples pour la détermination de la géographie sous la mer sans la connaissance a priori des informations de sous-surface selon un exemple de mode de réalisation ; [0050] la figure 11 illustre un schéma de principe d'un système qui peut mettre en oeuvre le procédé montré et décrit en faisant référence à la figure 10 selon un exemple de mode de réalisation ; et [0051] la figure 12 illustre la détermination de fenêtres pour une utilisation dans le procédé 100 décrit en relation avec la figure 10 selon un exemple de mode de réalisation. -20- 36-122/100306US DESCRIPTION DÉTAILLÉE [0052] Le concept de l'invention est décrit de manière plus complète ci-après en faisant référence aux dessins joints, sur lesquels des modes de réalisation du concept de l'invention sont montrés. Sur les dessins, la taille et les tailles relatives des couches et des régions peuvent être exagérées à des fins de clarté. Des numéros similaires font référence à des éléments similaires sur tous les dessins. Ce concept de l'invention peut, cependant, être mis en oeuvre sous de nombreuses formes différentes et ne devrait pas être interprété comme étant limité aux modes de réalisation exposés ici. Au lieu de cela, ces modes de réalisation sont fournis de sorte que cette description soit approfondie et complète, et transmette l'étendue du concept de l'invention aux hommes du métier. L'étendue de l'invention est par conséquent définie par les revendications jointes. Les modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec la terminologie et la structure d'un système de détermination de caractéristiques sismiques sous l'eau et sous le fond océanique en utilisant un système d'émetteurs et de récepteurs sismiques. Cependant, les modes de réalisation à examiner ensuite ne sont pas limités à ces systèmes, mais peuvent être appliqués à d'autres systèmes de détermination de caractéristiques sismiques qui utilisent les mêmes types ou des types similaires de systèmes de détermination sismiques. [0053] Une référence dans toute la spécification à « un mode de réalisation » signifie qu'une fonctionnalité, une structure, ou une caractéristique particulière décrite en relation avec un mode de réalisation est incluse dans au moins un mode de réalisation. Ainsi, l'apparition de l'expression « dans un mode de réalisation » à divers emplacements dans toute la spécification ne fait pas nécessairement référence au même mode de réalisation. En outre, la fonctionnalité, les structures, ou les caractéristiques particulières peuvent être combinées de n'importe quelle manière appropriée dans un ou plusieurs modes de réalisation. [0054] Plusieurs acronymes sont utilisés dans toute la spécification, dont la signification est -21- 36-122/100306US fournie comme suit : élimination de multiples associés à une surface (SRME « surface related multiple elimination » en terminologie anglo-saxonne) ; série de dispersion inverse (ISS « inverse scattering series » en terminologie anglo-saxonne) ; bidimensionnel (2D) ; mémoire à lecture seule (ROM) ; mémoire vive (RAM) ; mémoire à disque compact (CDROM) ; bus série universel (USB) ; disque compact (CD) ; disque polyvalent numérique (DVD) ; lecture/écriture (R/W) ; lecteur de disque dur (HDD) ; mémoire à lecture seule (ROM) ; électriquement effaçable (EE) ; mémoire ROM programmable (EEPROM) ; PROM effaçable par ultraviolets (UVPROM) ; mémoire vive (RAM) ; afficheurs à cristaux liquides (LCD) ; tubes à rayons cathodiques (CRT) ; entrées/sorties (E/S) ; système de positionnement mondial (GPS) ; réseau global (GAN) ; structures de réseau de portes programmables sur site (FPGA) ; éléments de circuit intégrés spécifiques à une application (ASIC) ; et système d'entrée/sortie de base (BIOS). 100551 Comme décrit dans l'article de 1998 attribué à Jakubowicz (ci-après « Jakubowicz »), un événement de multiples internes peut être construit en combinant trois composantes de champ d'onde, comme illustré sur la figure 4A, qui impliquent une convolution de deux traces suivie d'une corrélation d'une trace unique ou une certaine combinaison de celles-ci = _ dans laquelle Pi 402 est le champ d'onde principal de la surface à l'horizon j (Pj 402 est émis par la source 4b, réfléchi par l'horizon j, et arrive au niveau du récepteur 14a ; comme les hommes du métier peuvent le percevoir, le mot « horizon » est un autre terme pour une interface ou une frontière, entre des matériaux avec différents indices de réfraction n), Pk 406 représente le champ d'onde côté source réfléchi par l'horizon k et arrive au récepteur 14a (par conséquent, le récepteur 14a reçoit à la fois Pj 402 et Pk 406 comme indiqué sur la figure 4A) ; et P1404 est le -22- 36-122/100306US champ d'onde qui est généré par Pi 402 réfracté par l'horizon j, et ensuite réfléchi vers le haut et par l'horizon 1 (c'est-à-dire que Pi 404 est le signal réfracté émis par la source 4a, et réfléchi par l'horizon 1 vers le récepteur 14b). Dans l'équation (1), pi* est le conjugué complexe de P.i. A la fois pour la simplicité et la brièveté, le terme de signature de source est négligé dans l'équation (1). Les signes de sommation dans l'équation servent à garantir que, pour un milieu à n couches, tous les multiples internes possibles générés par l'horizon j en tant que réflecteur supérieur sont inclus. [0056] En utilisant l'ISS pour l'atténuation des multiples internes, la relation « inférieursupérieur-inférieur » implique que le diffuseur de points « supérieur » peut être à n'importe quelle pseudo profondeur tant qu'il est au-dessus du diffuseur de points « inférieur ». Par exemple, sur la figure 4A, l'horizon j est au-dessus de l'horizon k, lequel est au-dessus de l'horizon 1, et ainsi de suite. Selon un exemple de mode de réalisation préféré, une sommation supplémentaire peut être ajoutée dans une version modifiée de l'équation (1) pour qu'elle comprenne tous les horizons générant des multiples supérieurs possibles : (2), de sorte que le modèle de multiples M ne soit plus spécifique à un horizon. Avant d'examiner l'équation 2, cependant, la nomenclature va être abordée, afin d'aider à la compréhension des exemples de modes de réalisation décrits ici. Initialement, selon un exemple de mode de réalisation, les données sismiques reçues sont divisées en M fenêtres temporelles, ce qui est examiné plus en détail ci-dessous. Comme montré sur la figure 4A, Pk est un champ d'onde réfléchi par l'horizon k, lequel est un objet physique réel, fixe, inamovible (pratiquement) ; (par exemple, ce pourrait être une couche thermique dans l'océan, une frontière océan-sol, ou une -23- 36-122/100306US frontière entre le sable et la roche, parmi de nombreux autres types), Pwk, cependant, est défini selon un exemple de mode de réalisation en tant que champs d'onde côté source réfléchis par les horizons dans une fenêtre temporelle Wk ; ainsi, il s'agit d'une construction mathématique basée sur des données physiques, réelles reçues de multiples signaux, mais qui est déterminée pendant le traitement des données reçues. Selon un exemple de mode de réalisation, le champ d'onde Pwk est obtenu en traitant les données qui sont mémorisées dans une fenêtre Wk. Selon un exemple de mode de réalisation, le champ d'onde Pwi est obtenu en traitant les données qui sont mémorisées dans une fenêtre Wj. Et, selon un exemple de mode de réalisation, le champ d'onde Pm est obtenu en traitant les données qui sont mémorisées dans une fenêtre W1. Dans l'équation (2), les champs d'onde ne correspondent pas à un événement unique, mais à un groupe d'événements dans une fenêtre prédéterminée. Par exemple, Pwi représente le champ d'onde principal réfléchi vers le haut par tous les horizons dans une fenêtre Wi (l'horizon Wi défini arbitrairement représenté sur la figure 4B, par opposition à l'horizon j tel que représenté sur la figure 4A). Pwk représente le champ d'onde côté source réfléchi à partir d'une fenêtre Wk 410, qui, dans ce premier cas, peut être défini en tant que champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi par une fenêtre Wi. Pwi est un champ d'onde côté récepteur qui représente une réflexion vers le bas d'un multiple interne réfléchi par une fenêtre W1. [0057] Comme montré sur la figure 12, chaque trace sismique est segmentée en M parties horizontalement. Selon qu'elle est responsable d'une réflexion vers le haut ou vers le bas du multiple interne, une fenêtre appropriée sera sélectionnée pour garantir que la relation « inférieure-supérieure-inférieure » est satisfaite. Cela est illustré schématiquement sur la figure 12 en utilisant un rassemblement de tirs unidimensionnel. Après que la fenêtre supérieure, WA, a été déterminée (à cause de l'opération des conjugués complexes dans l'équation (2), les traces semblent être retournées verticalement ; ainsi, sur la figure 12, la fenêtre supérieure est WA, retournée, montrant l'opération des conjugués complexes), deux fenêtres inférieures (WB et -24- 36-122/100306US Wc) sont alors choisies pour satisfaire à l'équation (2). Comme les hommes du métier peuvent maintenant l'apprécier, la figure 12 représente simplement l'une des nombreuses combinaisons possibles de fenêtres. Pour une fenêtre supérieure spécifique, de nombreuses fenêtres « inférieures » peuvent être utilisées tant qu'elles satisfont à la relation « inférieure-supérieureinférieure ». Le processus, comme décrit davantage ci-dessous, est ensuite répété avec d'autres fenêtres supérieures. Par conséquent, pour obtenir l'équation (2), une approche à fenêtre glissante haut-bas est utilisée comme représenté sur la figure 4B. C'est-à-dire que la fenêtre Wi 408 est située à la position la plus haute, et est ensuite glissée jusqu'à un emplacement inférieur plusieurs fois, formant, successivement, des fenêtres Wk 410, W1 412, Wrn 414, et ainsi de suite, comme déterminé par l'équation (2), en prenant en considération les sommations dans la sommation qui sont incluses dans l'équation (2). La figure 4B illustre, selon un exemple de mode de réalisation, qu'après que la fenêtre W1408 a été traitée pour déterminer où les réflecteurs supérieurs pour les multiples internes peuvent être intégrés, une fenêtre Wk 304 peut maintenant être considérée comme une autre plage de profondeur dans laquelle les réflecteurs supérieurs peuvent être trouvés, et ainsi de suite pour chaque fenêtre, W1 412, ensuite Wm 414, et ainsi de suite. En répétant le glissement vers le bas des fenêtres, tous les multiples internes peuvent être prédits sans devoir identifier les horizons générant des multiples selon un exemple de mode de réalisation. 100581 Comme les hommes du métier peuvent l'apprécier, il y a plusieurs considérations à prendre en compte lors de la mise en oeuvre de l'équation (2). D'abord, la contrainte « inférieuresupérieure-inférieure » est satisfaite en comparant le temps de propagation des rassemblements de tirs communs et de récepteurs effectués en vue d'obliquité (« Normal Move-out » en terminologie anglo-saxonne). Selon un autre exemple de mode de réalisation, l'incorporation de la contrainte « inférieure-supérieure-inférieure » rend le traitement efficace du point de vue calcul et elle est appliquée aux données prémigrées. Une correction d'obliquité fait référence à -25- 36-122/100306US une fonction du temps et de décalage qui peut être utilisée dans un traitement sismique pour compenser les effets d'obliquité, ou le retard des temps d'arrivée de réflexion lorsque des géophones et des points de tir sont décalés les uns par rapport aux autres. Un décalage est défini, dans une acquisition sismique de surface, en tant que distance horizontale d'une source à un récepteur. Il est bien connu que des décalages entre des sources sismiques et leurs récepteurs correspondants créent un retard, ou une inclinaison, dans le temps d'arrivée d'une réflexion qui peut être corrigé avant un empilage et qui peut être utilisé pour déterminer une vitesse. [0059] Cependant, l'utilisation de la correction d'obliquité pour satisfaire à la contrainte « inférieure-supérieure-inférieure » n'est pas aussi précise que son utilisation dans une pseudo profondeur (voir par exemple, Nita, B. G. et d'autres, 2007, « Inverse Scattering Internai Multiple Attenuation Algorithm : An Analysis of the Pseudo Depth and Time Monotonicity Requirements », 77th Meeting, SEG, Expanded Abstracts, pages 2461 à 2464, dont le contenu entier est incorporé ici par voie de référence) et une vitesse est nécessaire dans le processus. Néanmoins, comme cela sera compris par les hommes du métier, des informations de vitesse pour une extrapolation de décalage proche est nécessaire pour une SRME, ainsi cela n'ajoute pas de spécification supplémentaire quant aux informations de sous-surface. Une SRME, comme cela peut être apprécié par les hommes du métier, utilise les données sismiques enregistrées pour prédire et soustraire de manière itérative la série de multiples. L'avantage principal d'une SRME est qu'elle ne nécessite absolument aucune information de sous-surface ; les multiples sont prédits complètement ou sensiblement complètement à partir des données. De plus, pour des structures de sous-surface relativement peu compliquées, le critère utilisant le temps de propagation peut être valide. [0060] La deuxième considération, selon un exemple de mode de réalisation, est que la séparation entre des générateurs de multiples (c'est-à-dire, des sources 4) doit être supérieure à la longueur de fenêtre. Il n'y a aucune telle contrainte correspondante en ISS. Pour des générateurs -26- 36-122/100306US de multiples internes, tels que des veines de charbon qui ont habituellement des strates étroitement regroupées, de très courtes fenêtres peuvent être nécessaires. Dans l'examen qui suit, les performances d'un système et d'un procédé selon un exemple de mode de réalisation sont analysées lorsque ces deux considérations sont prises en compte, et les résultats sont comparés au procédé ISS connu. Selon un autre exemple de mode de réalisation, la taille de fenêtre est basée sur la longueur des ondelettes de source, et de plus le nombre de fenêtres peut être déterminé en divisant la longueur des traces sismiques par la taille de la fenêtre. Alors que l'approche de Jakubowicz ne nécessite pas les deux contraintes examinées ci-dessus, le procédé de Jakubowicz est gêné par le besoin d'identifier les horizons générant des multiples, c'est-à-dire par une connaissance a priori. Comme les hommes du métier peuvent l'apprécier, la nécessité d'une telle connaissance a priori signifie que l'approche de Jakubowicz est plus coûteuse et moins efficace que les exemples de modes de réalisation décrits ici. [0061] Exemples de données synthétiques et de site [0062] Afin d'illustrer les avantages des exemples de modes de réalisation du système et du procédé examinés ici, à la fois un procédé basé sur l'ISS et le procédé selon des exemples de modes de réalisation ont été appliqués à des données synthétiques, bidimensionnelles (2D), ainsi qu'à des données de site d'origine à partir desquelles les données synthétiques ont été générées. Les figures 5A et 5B illustrent des données de site réelles et des données synthétiques générées à partir des données réelles, respectivement, acquises à partir de la découverte de pétrole de Tupi à proximité de Santos Basin, au large du Brésil, dans lesquelles des multiples internes importants sont évidents. Le site pétrolier de Tupi est un grand site pétrolier situé dans le Santos Basin, environ 250 km (160 milles) de la côte de Rio de Janeiro, au Brésil. La figure 5A montre les données de site provenant d'une ligne proche de la découverte de Tupi. Une série de contrastes d'impédances au-dessus du niveau antésalifère peut être observée, par exemple le fond marin, le dessus du niveau salifère, et les structures salifères en couches, ou des évaporites en couches -27- 36-122/100306US (entre autres raisons, les sites pétroliers de Tupi sont importants en ce que le pétrole découvert est situé au niveau antésalifère, ou au-dessous du niveau salifère). Tous ces réflecteurs peuvent être des générateurs de multiples internes. Les données synthétiques, affichées sur la figure 5B, ont été générées par une modélisation acoustique en utilisant la vitesse correspondante des données de site. Le modèle de densité a été mis à jour de manière itérative jusqu'à ce que les événements dans les données synthétiques ressemblent à ceux dans les données de site (notez la similitude des deux figures, les figures 5A et 5B). Cependant, les multiples internes ont été amplifiés dans les données synthétiques de sorte qu'ils puissent être facilement identifiés. Dans ce cas, les données synthétiques sont utilisées pour illustrer et examiner l'efficacité des deux procédés. 100631 En utilisant une longueur de fenêtre de 25 échantillons avec une longueur de superposition de 4 échantillons, un modèle de tous les multiples internes, comme représenté sur la figure 6A, a été prédit en utilisant le procédé à fenêtre glissante selon des exemples de modes de réalisation sans informations de sous-surface (c'est-à-dire, sans informations de sous-surface a priori). A titre d'exemple non limitatif seulement, un exemple de longueur d'échantillons peut être de 25 échantillons avec une fréquence d'échantillonnage de traitement de 4 millisecondes, ce qui signifie que la longueur de fenêtre dans le temps, ou temporellement, est segmentée en parties d'environ 100 millisecondes. Cependant, la « longueur » dans ce sens correspond en réalité à la profondeur de chaque fenêtre, comme examiné plus en détail ci-dessous. Précédemment, le même ensemble de données a été testé avant mais avec le besoin d'identifier les horizons générant des multiples (voir Griffiths, M. et d'autres, 2011, « Applications of Interbed Multiple Attenuation », The Leading Edge, 30, n° 8, pages 906 à 912, dont le contenu entier est incorporé ici par voie de référence). La figure 6B illustre un deuxième modèle de tous les multiples internes des sites pétroliers de Tupi prédits en utilisant le procédé ISS. On peut observer dans les deux prédictions que les multiples autour du niveau antésalifère, qui sont -28- 36-122/100306US principalement associés au fond marin, et les multiples internes générés par la partie supérieure du niveau salifère et des évaporites en couches sont bien modélisés par les deux procédés. Même si la structure de l'un des horizons générant des multiples supérieurs, dessus de niveau salifère, est plutôt compliquée, le système et le procédé selon des exemples de modes de réalisation sont encore capables de satisfaire à la relation « inférieure-supérieure-inférieure » parce que les horizons générant des multiples supérieurs sont bien séparés dans ce cas. La figure 7A illustre les données de la figure 6A auxquelles une soustraction adaptative a été appliquée pour supprimer des multiples internes selon un exemple de mode de réalisation. La figure 7B illustre les données de la figure 6B (générées en utilisant le procédé à série de dispersion inverse) auxquelles une soustraction adaptative a été appliquée pour supprimer des multiples internes. Selon un exemple de mode de réalisation, une soustraction adaptative implique une adaptation du spectre d'amplitude du modèle de multiples internes à celui de l'entrée avant l'exécution de la soustraction. L'adaptation qui apparaît signifie qu'une compensation se produit entre l'amplitude et un contenu fréquentiel entre des informations modélisées et des informations réelles. A la suite de l'adaptation, une soustraction est effectuée. Par conséquent, la figure 7A résulte d'une adaptation d'amplitude et de fréquence entre les données montrées sur les figures 5A et 6A, et ensuite les données de la figure 6A sont soustraites de celles de la figure 5A pour fournir les résultats montrés sur la figure 7A. Un processus similaire est effectué entre les figures 5B, 6B, et la figure résultante 7B. On peut voir que ces multiples internes au-dessous du fond du niveau salifère qui sont amplifiés par l'effet de focalisation qui piège une énergie de multiple dans les mini bassins (voir, Pica, A. et d'autres, 2008, « Wave Equation Based Internai Multiple Modeling in 3D », 78th Meeting, SEG, Expanded Abstracts, pages 2476 à 2480, dont le contenu entier est incorporé ici par voie de référence) sont particulièrement bien atténués dans les deux résultats. 100641 La figure 8A illustre une section agrandie de la figure 4A, la figure 8B illustre les mêmes données à la suite d'une mise en oeuvre du procédé selon un exemple de mode de réalisation, et la -29- 36-122/100306US figure 8C illustre la différence entre les figures 8A et 8B (c'est-à-dire, les données de la figure 8B soustraites des données de la figure 8A). En outre, la figure 9A illustre les mêmes données provenant de la figure 8A à la suite d'un traitement utilisant le procédé à série de dispersion inverse, et la figure 9B illustre la différence entre les figures 9A et 8A (permettant ainsi une comparaison entre les figures 8C et 9B). [0065] Les résultats de l'atténuation des multiples internes correspondants provenant des données de site (développées en utilisant le procédé selon un exemple de mode de réalisation, et en utilisant l'ISS) sont affichés sur les figures 8B et 9A, respectivement, après que les multiples de surfaces libres ont été retirés en utilisant la SRME. On peut voir à partir du tracé de différence de la figure 8C et du tracé de différence de la figure 9B que la plupart des multiples internes ont été modélisés et atténués par les deux procédés, bien qu'ils ne soient pas aussi nets que dans le cas des données synthétiques parce que les multiples sont beaucoup plus faibles dans les données de site. De plus, le bruit dans les données de site entrées dégrade les performances de prédiction et de soustraction. L'extension de l'utilisation du système et des procédés selon des exemples de modes de réalisation à une opération tridimensionnelle peut améliorer la précision de la modélisation de multiples. [0066] Selon des exemples de modes de réalisation examinés ici, une approche basée sur une localisation itérative des horizons générant des multiples a été proposée pour prédire des multiples internes dans l'exploration sismique sous la mer de gisements d'hydrocarbure. Les systèmes et les procédés examinés ici selon des exemples de modes de réalisation sont une variante appropriée et rentable pour une atténuation de multiples internes sans informations de sous-surface. [0067] La figure 10 illustre un organigramme de procédé 100 pour éliminer sensiblement l'influence de réflexions de multiples lors de la détermination de la géographie sous la mer sans -30- 36-122/100306US la connaissance a priori d'informations de sous-surface selon un exemple de mode de réalisation. Le procédé 100 commence avec l'étape 102, à laquelle des signaux sismiques sont générés par des sources 4. A l'étape 104, des données brutes sont reçues de tous les récepteurs 14 et mémorisées dans un dispositif de mémorisation approprié. Les données brutes sont traitées à l'étape 104 seulement dans la mesure où les multiples de surface sont supprimés. Selon un exemple de mode de réalisation, une manière de supprimer des multiples de surface est la SRME, examinée brièvement ci-dessus. D'autres procédés connus peuvent comprendre ceux utilisant une transformation de Radon et une modélisation de champ d'onde. A l'étape 108, M fenêtres sont définies pour la zone géographique présentant un intérêt. Les fenêtres sont définies en fonction d'une longueur 1 et d'une profondeur d. Selon un exemple de mode de réalisation, la longueur de la fenêtre est importante étant donné que la longueur de fenêtre 1 doit être inférieure à la distance totale entre les première et dernière sources 4. La profondeur de la fenêtre est déterminée sur la base de la vitesse du son dans l'eau de mer, qui peut être généralement supposée être d'environ 1500 mètres par seconde Ainsi, la taille de fenêtre a à la fois des dimensions de temps et de longueur. Si le navire 2 est environ 1500 mètres au-dessus du fond de l'océan, une onde sonore prendra environ deux secondes pour se propager de la source 4 jusqu'au fond océanique 42, et ensuite de retour jusqu'au récepteur 14. Selon un premier exemple de mode de réalisation, chaque fenêtre peut être définie pour correspondre à environ 100 millisecondes. Cela correspond à environ 150 mètres en profondeur, et il y aura environ dix fenêtres entre la surface de l'océan et le fond de l'océan (ainsi M = 10). Cependant, les hommes du métier peuvent apprécier que la fenêtre peut être définie comme étant pratiquement de n'importe quelle profondeur (c'est-à-dire, un intervalle de temps), dans lequel les facteurs limitant la décision du nombre de fenêtres à mettre en oeuvre peuvent être une taille d'échantillon, une vitesse et un temps de traitement, ainsi que des limitations de mise en mémoire. De plus, cependant, il peut y avoir également des contraintes en termes de résolution des données échantillonnées, de sorte que la profondeur de fenêtre s'approche d'une limitation pratique basée -31- 36-122/100306US sur des fréquences d'échantillonnage de convertisseur analogique-numérique, parmi d'autres facteurs. [0068] A la suite de l'étape 108, à laquelle les M fenêtres sont définies, le procédé 100 passe à l'étape 110 à laquelle les données sont attribuées à différentes fenêtres en fonction de l'instant d'arrivée. Si, par exemple, il existait 100 fenêtres, w1 à wloo, alors une centaines de champs d'onde seraient reconstruits, et appelés Pwi, Pw2, Pw3, jusqu'à Pwioo. Ultérieurement, dans le procédé, les ensembles de trois champs d'onde sont corrélés avec Pwi, Pwk, et Pwi et utilisés dans l'équation 2 pour déterminer un ensemble de multiples internes selon un exemple de mode de réalisation. Ensuite, à l'étape 112, le procédé 100 fixe n = 1, et passe à l'étape 114, à laquelle Pwi est fixé égal à Pw(n). Pwk, par conséquent, est toujours la fenêtre la plus haute, et Pwk, Pwi sont fixés pour être les deux fenêtres les plus basses, selon une sommation effectuée de manière itérative selon l'équation (2), comme examiné plus en détail selon un exemple de mode de réalisation ci-dessous. A l'étape 116, selon un exemple de mode de réalisation, man) est calculé selon une version modifiée de l'équation 2. Dans l'organigramme joint, figure 10, l'équation 2 est modifiée pour retirer la sommation la plus à gauche, de sorte que le calcul de M (où M est le modèle de multiples internes, c'est-à-dire, les données illustrant les multiples internes résultants 51) peut être montré sous un format d'organigramme ; c'est-à-dire que la sommation la plus à gauche, de = 1 à wi = wn, est représentée par la boucle itérative qui calcule les autres sommations pour chaque fenêtre définie (examinée ci-dessus), et la boucle indique que les sommations sont effectuées pour chaque fenêtre définie selon un exemple de mode de réalisation. A l'étape 116, selon l'équation (2), Pwi est fixé égal à Pw(l), et dans une première itération de l'équation (2), Pwk est fixé égal à Pw2, et Pwi est également fixé égal à Pw2 ; ensuite, le produit est calculé, et Pwi est itéré de 2 à 100 ; à la suite de ces séries de calcul de produit, Pwk est itéré de 2 à 100, et pour chaque itération de Pwk, Pw1 est également itéré de Pw2 à Pwioo ; de cette manière, la sommation dans la sommation se manifeste par des boucles itératives. Voir, par ailleurs, l'examen suivant -32- 36-122/100306US immédiatement l'introduction de l'équation (2) montrée ci-dessus. [0069] Une fois que Pwk et Pw1 ont été traités de manière itérative de n = 2 à 100 (parce que, à des fins d'exemple uniquement, M était fixé égal à 100), M est calculé (comme examiné ci-dessous) et ensuite le procédé 100 détermine, à l'étape 118, un champ de données de multiples total, M, selon l'équation m = m(in) ; cela ajoute les champs d'onde des fenêtres précédentes au nouveau, et répète le processus itératif jusqu'à ce que tous les champs d'onde des fenêtres définies aient été accumulés. [0070] A l'étape 120, le procédé 100 ajoute 1 à n. A l'étape de décision 122, n est vérifié pour déterminer si tous les champs d'onde ont été utilisés pour déterminer le modèle de multiples internes 51 M. Tous les champs d'onde auront été utilisés lorsque n atteint M-1 (alors, Pwi = Pw99, et Pwk et Pwi sont égaux à Pw100). Si ce n'est pas le cas (trajet « non » de l'étape de décision 122), le procédé 100 retourne à l'étape 114, et l'ensemble suivant des champs d'onde pour une nouvelle fenêtre est obtenu à partir de la mémoire de données de fenêtre. Selon un exemple de mode de réalisation, le champ d'onde le plus haut, Wff, est la fenêtre suivante sous (w'±i), et lorsque n = 2 de sorte que Pwi = Pw2, Pwk et Pm commence en Pw3 comme examiné précédemment. Ensuite, en présumant que toutes les données de fenêtre n'ont pas été utilisées, le processus de calcul des multiples internes M est effectué de nouveau. Si tous les champs d'onde pour toutes les fenêtres ont été pris en compte pour la détermination du modèle M de multiples internes 51, alors le procédé 100 continue à l'étape 124 (trajet « oui » à partir de l'étape de décision 122), et les données vraies sont déterminées en ajoutant M aux données brutes (en réalité une soustraction, parce que m(j') est défini comme étant un négatif de la sommation), et le résultat est une représentation réelle de la zone géographique présentant un intérêt avec des multiples réduits et/ou sensiblement éliminés des données brutes reçues. -33- 36-122/100306US [0071] Selon un exemple de mode de réalisation préféré, deux critères doivent être satisfaits afin d'utiliser le procédé 100 : d'abord, les critères inférieure-supérieure-inférieure examinés plus en détail ci-dessus doivent être supposés avoir été satisfaits, et ensuite, la longueur de fenêtre doit être inférieure à la séparation entre les générateurs de multiples, ou les sources 14. C'est-à-dire que la fenêtre a non seulement une profondeur (dans le temps ou en mètres), mais également une distance (de nouveau en mètres). [0072] Un système et un procédé qui peuvent prédire des multiples internes dans des données sismiques marines ou terrestres sans nécessiter des informations a priori concernant la sous-surface de la terre sont présentés ici. Ils sont destinés à être utilisés après la suppression des multiples associés à une surface. Le système et le procédé séparent d'abord les données sismiques en différentes fenêtres sur la base du temps de propagation du champ d'onde de la source aux récepteurs. Ensuite, une convolution est effectuée entre les traces dans les fenêtres qui sont responsables des réflexions vers le haut des multiples internes et une corrélation avec les traces dans la première fenêtre temporelle qui sont responsables des réflexions vers le bas des multiples internes correspondants. En répétant le processus qui comprend les fenêtres temporelles subséquentes, tous les ordres de tous les multiples internes peuvent être prédits sans informations de sous-surface. Selon certains modes de réalisation, les longueurs de fenêtre peuvent être différentes pour les trois fenêtres wj, wk et w1, par exemple, en pratique, wk et wi peuvent être suffisamment plus long que wi. [0073] La figure 11 illustre un système d'acquisition de données sismiques (système) 200 approprié pour une utilisation pour mettre en oeuvre le procédé 100 pour éliminer sensiblement l'influence de réflexions de multiples lors de la détermination de la géographie sous la mer sans la connaissance a priori d'informations de sous-surface selon un exemple de mode de réalisation. Le système 200 comprend, entre autres éléments, un serveur 201, une interface source/récepteur 202, un bus de communication/données interne (bus) 204, un ou des processeurs 208 (les -34- 36-122/100306US hommes du métier peuvent apprécier que, dans des systèmes de serveur modernes, un traitement en parallèle devient de plus en plus répandu et bien qu'un processeur unique ait été utilisé dans le passé pour mettre en oeuvre de nombreuses ou au moins plusieurs fonctions, il est plus commun actuellement d'avoir un processeur dédié unique pour certaines fonctions (par exemple, des processeurs de signaux numériques)), un port de bus série universel (USB) 210, un lecteur à lecture/écriture (RW) de disque compact (CD)/disque polyvalent numérique (DVD) 212, un lecteur de disquette 214 (bien que moins utilisé actuellement, de nombreux serveurs comprennent encore ce dispositif), et une unité de mémorisation de données 232. L'unité de mémorisation de données 232 elle-même peut comprendre un lecteur de disque dur (HDD) 216 (celui-ci peut comprendre un support de mémorisation magnétique classique, mais, comme cela devient de plus en plus répandu, peut comprendre des dispositifs de stockage de masse de type clé USB 224, entre autres types), un ou des dispositifs de mémoire à lecture seule (ROM) 218 (ceux-ci peuvent comprendre des dispositifs de mémoire ROM programmable effaçable électriquement (EE) (EEPROM), des dispositifs PROM effaçables par ultraviolets (UV) (UVPROM), entre autres types), et des dispositifs de mémoire vive (RAM) 220. Un dispositif de clé USB 224 peut être utilisé avec le port USB 210, et les disques CD/DVD 234 peuvent être utilisés avec un dispositif R/W CD/DVD 212. Les disquettes 237 peuvent être utilisées avec un dispositif de lecture de disquette 214. Chacun des dispositifs de mise en mémoire, ou des supports de mise en mémoire (216, 218, 220, 224, 234 et 237, entre autres types), peut contenir des parties ou des composants, ou dans son intégralité, un logiciel 236 de code de programmation logicielle exécutable (logiciel) 236 qui peut mettre en oeuvre le procédé décrit ici. En outre, le processeur 208 lui-même peut contenir un type ou différents types de dispositifs de mise en mémoire (très probablement, mais non d'une manière limitative, un support de mise en mémoire vive 220) qui peut mémoriser la totalité ou certains des composants de logiciel 236. [0074] En plus des composants décrits ci-dessus, le système 200 comprend également une -35- 36-122/100306US console d'utilisateur 234, qui peut comprendre un clavier 228, un afficheur 226 et une souris 230. La totalité de ces composants sont connus des hommes du métier, et cette description comprend toutes les variantes connues et futures de ces types de dispositifs. L'afficheur 226 peut être n'importe quel type d'afficheur ou d'écran de présentation connu, tel que des afficheurs à cristaux liquides (LCD), des afficheurs au plasma, des tubes à rayons cathodiques (CRT), entre autres. La console d'utilisateur 234 peut comprendre un ou plusieurs mécanismes d'interface utilisateur tels qu'une souris, un clavier, un microphone, un pavé tactile, un écran tactile, un système de reconnaissance vocale, entre autres dispositifs inter-communicants interactifs. [0075] La console d'utilisateur 234, et ses composants s'ils sont fournis séparément, s'interfacent avec le serveur 201 par l'intermédiaire d'une interface d'entrée/sortie (E/S) de serveur 222, qui peut être un port de communication RS232, Ethernet, USB ou d'un autre type, ou qui peut comprendre la totalité ou certains de ceux-ci, et comprend en outre n'importe quel autre type de moyens de communication, actuellement connus ou autrement développés. Le système 200 peut comprendre en outre un dispositif d'émission-réception de satellite de communications/système de positionnement mondial (GPS) 238, auquel au moins une antenne 240 (selon un exemple de mode de réalisation, il y aurait au moins une antenne de réception GPS seule, et au moins une antenne de communication bidirectionnelle de satellite séparée) est connectée électriquement. Le système 200 peut accéder à Internet 242, soit par une connexion câblée, par l'intermédiaire de l'interface d'entrée/sortie 222 directement, soit par une liaison sans fil par l'intermédiaire de l'antenne 240 et de l'émetteur-récepteur 238. [0076] Le serveur 201 peut être couplé à d'autres dispositifs informatiques, tels que ceux qui mettent en oeuvre ou commandent l'équipement du navire 2, par l'intermédiaire d'un ou de plusieurs réseaux. Le serveur 201 peut être une partie d'une configuration de réseau plus grande comme dans un réseau global (GAN) (par exemple, Internet 242), qui permet finalement une connexion à diverses lignes câblées. -36- 36-122/100306US [0077] Selon un autre exemple de mode de réalisation, le système 200, qui est ostensiblement conçu pour une utilisation dans l'exploration sismique, s'interfacera avec une ou plusieurs sources 4 et un ou plusieurs récepteurs 14. Ceux-ci, comme décrit précédemment, sont attachés à des flûtes 6, auxquelles sont également attachées des oiseaux 13a,b qui sont utiles pour maintenir le positionnement. Comme examiné en outre précédemment, les sources 4 et les récepteurs 14 peuvent communiquer avec le serveur 200 soit par un câble électrique qui fait partie de la flûte 6, soit par l'intermédiaire d'un système sans fil qui peut communiquer par l'intermédiaire de l'antenne 240 et de l'émetteur-récepteur 238 (collectivement décrits en tant que conduit de communication 246). [0078] Selon d'autres exemples de modes de réalisation, la console d'utilisateur 234 (par exemple, un clavier, des boutons, des commutateurs, un écran tactile et/ou une manette) fournit des moyens pour qu'un personnel entre des commandes et une configuration dans le système 200. Le dispositif d'affichage 226 peut être utilisé pour montrer : la position de la flûte 6, des représentations visuelles de données acquises, des informations d'état de la source 4 et du récepteur 14, des informations d'étude et d'autres informations importantes pour le processus d'acquisition de données sismiques. Une unité d'interface de sources et de récepteurs 202 peut recevoir les données sismiques d'hydrophone, du récepteur 14, par l'intermédiaire du conduit de communication de flûte 248 (examiné ci-dessus) qui peut être une partie de la flûte 6, ainsi que des informations de position de la flûte 6 provenant des oiseaux 13 ; la liaison est bidirectionnelle de sorte que les commandes peuvent également être envoyées aux oiseaux 13 pour maintenir un positionnement correct de la flûte. L'unité d'interface de sources et de récepteurs 202 peut également communiquer de manière bidirectionnelle avec les sources 4 par l'intermédiaire du conduit de communication de flûte 248 qui peut faire partie de la flûte 6. Des signaux d'excitation, des signaux de commande, des signaux de sortie et des informations d'état, associés à la source 4, peuvent être échangés par le conduit de communication de flûte 248 entre -37- 36-122/100306US le système 200 et la source 4. [0079] Le bus 204 fournit un trajet de données : pour le transfert et la mémorisation de données qui proviennent soit des capteurs de source, soit des récepteurs de flûte ; pour que le processeur 208 accède à des données mémorisées contenues dans la mémoire d'unité de mémorisation de données 232 ; pour que le processeur 208 envoie des informations pour un affichage visuel sur l'afficheur 226 ; ou pour que l'utilisateur envoie des commandes à des programmes d'exploitation de système qui résideraient soit dans le processeur 208, soit dans l'unité d'interface de sources et de récepteurs 202. [0080] Le système 200 peut être utilisé pour mettre en oeuvre le procédé 100 pour éliminer sensiblement l'influence de réflexions de multiples lors de la détermination de la géographie sous la mer sans la connaissance a priori d'informations de sous-surface selon un exemple de mode de réalisation. Un matériel, un micrologiciel, un logiciel ou une combinaison de ceux-ci peut être utilisé pour effectuer les diverses étapes et opérations décrites ici. Selon un exemple de mode de réalisation, le logiciel 236 pour exécuter les étapes examinées ci-dessus peut être mémorisé et distribué sur des dispositifs de mémorisation multimédias tels que les dispositifs 216, 218, 220, 224, 234 et/ou 237 (décrits ci-dessus). Ces supports de mémorisation peuvent être insérés dans et lus par des dispositifs tels qu'un port USB 210, un dispositif R/W CD/DVD 212, un lecteur de disquette 214, entre autres. [0081] Selon un exemple de mode de réalisation, une mise en oeuvre du procédé 100 peut être effectuée dans un processeur 208 dédié. Les hommes du métier dans le domaine des modes de réalisation peuvent apprécier que cette fonctionnalité peut être conçue dans divers types d'éléments de circuit, comprenant, mais sans y être limités, des structures de réseau de portes programmables sur site (FPGA), des éléments de circuit intégrés spécifiques à une application (ASIC), des systèmes basés sur des microprocesseurs, parmi d'autres types. Un examen détaillé -38- 36-122/100306US des divers types de mises en oeuvre de circuit physiques n'apporte sensiblement rien à la compréhension des modes de réalisation, et, à ce titre, celui-ci a été omis par souci de brièveté et de clarté. Cependant, comme cela est également connu des hommes du métier, les systèmes et les procédés examinés ici peuvent être mis en oeuvre comme examiné, et peuvent comprendre en outre des dispositifs programmables. 100821 Ces dispositifs programmables et/ou d'autres types d'éléments de circuit tels qu'examinés précédemment peuvent comprendre une unité de traitement, une mémoire de système et un bus de système qui couple divers composants de système comprenant la mémoire de système à l'unité de traitement. Le bus de système peut être n'importe lequel de plusieurs types de structures de bus comprenant un bus de mémoire ou un contrôleur de mémoire, un bus de périphériques, et un bus local utilisant n'importe laquelle d'un grand nombre d'architectures de bus. En outre, divers types de supports pouvant être lus par un ordinateur peuvent être utilisés pour mémoriser des instructions programmables. Un support pouvant être lu par un ordinateur peut être n'importe quel support disponible auquel un accès peut être effectué par l'unité de traitement. A titre d'exemple, et non de limitation, un support pouvant être lu par un ordinateur peut comprendre un support de mémorisation d'ordinateur et un support de communication. Le support de mémorisation d'ordinateur comprend un support volatile et non volatile, ainsi qu'un support amovible et non amovible, mis en oeuvre dans n'importe quel procédé ou technologie pour la mémorisation d'informations telles que des instructions pouvant être lues par un ordinateur, des structures de données, des modules de programme, ou d'autres données. Le support de mémorisation d'ordinateur comprend, mais sans y être limité, une mémoire vive, une mémoire morte, une mémoire EEPROM, une mémoire flash ou une autre technologie de mémoire, un CDROM, des DVD ou une autre mémoire à disque optique, des cassettes magnétiques, une bande magnétique, une mémoire à disque magnétique ou d'autres dispositifs de mémorisation magnétiques, ou n'importe quel autre support qui peut être utilisé pour mémoriser -39- 36-122/100306US les informations souhaitées et auquel un accès peut être effectué par l'unité de traitement. Les supports de communication peuvent mettre en oeuvre des instructions pouvant être lues par un ordinateur, des structures de données, des modules de programme ou d'autres données dans un signal de données modulé tel qu'une onde porteuse ou un autre mécanisme de transport et peuvent comprendre n'importe quel support de distribution d'informations approprié. [0083] La mémoire de système peut comprendre un support de mémorisation d'ordinateur sous la forme d'une mémoire volatile et/ou non volatile telle qu'une mémoire morte et/ou une mémoire vive. Un système d'entrée/sortie de base (BIOS), contenant les sous-programmes de base qui aident au transfert des informations entre des éléments connectés au processeur et entre ceux-ci et le processeur, par exemple pendant un démarrage, peut être stocké dans une mémoire. La mémoire peut contenir également des données et/ou des modules de programme qui sont immédiatement accessibles et/ou actuellement mis en oeuvre par l'unité de traitement. A titre d'exemple non limitatif, la mémoire peut également comprendre un système d'exploitation, des programmes d'application, d'autres modules de programme, et des données de programme. [0084] Le processeur peut également comprendre d'autres supports de mémorisation d'ordinateur amovibles/non amovibles et volatils/non volatils. Par exemple, le processeur peut accéder à un lecteur de disque dur qui lit ou écrit sur un support magnétique non amovible non volatil, un lecteur de disque magnétique qui lit ou écrit sur un disque magnétique, amovible, non volatil et/ou un lecteur de disque optique qui lit ou écrit sur un disque optique amovible, non volatil tel qu'un CDROM ou un autre support optique. D'autres supports de mémorisation d'ordinateur amovibles/non amovibles, volatils/non volatils qui peuvent être utilisés dans l'exemple d'environnement d'exploitation comprennent, mais sans y être limités, des cassettes de bande magnétique, des cartes de mémoire flash, des disques polyvalents numériques, des bandes vidéo numériques, une mémoire vive à semi-conducteurs, une mémoire morte à semi- conducteurs et similaire. Un lecteur de disque dur peut être connecté au bus de système par une -40- 36-122/100306US interface de mémoire non amovible telle qu'une interface, et un lecteur de disque magnétique ou un lecteur de disque optique peut être connecté au bus de système par une interface de mémoire amovible, telle qu'une interface. [0085] Les exemples de modes de réalisation examinés ici peuvent également être mis en oeuvre en tant que codes pouvant être lus par un ordinateur sur un support pouvant être lu par un ordinateur. Le support pouvant être lu par un ordinateur peut comprendre un support d'enregistrement pouvant être lu par un ordinateur et un support de transmission pouvant être lu par un ordinateur. Le support d'enregistrement pouvant être lu par un ordinateur est n'importe quel dispositif de stockage de données qui peut mémoriser des données qui peuvent être lues par la suite par un système informatique. Des exemples du support d'enregistrement pouvant être lu par un ordinateur comprennent une mémoire morte, une mémoire vive, des CDROM et des dispositifs de mémorisation de données généralement optiques, des bandes magnétiques, des clés USB, et des disquettes. Le support pouvant être lu par un ordinateur peut également être réparti sur des systèmes informatiques couplés à un réseau de sorte que le code pouvant être lu par un ordinateur soit mémorisé et exécuté d'une manière répartie. Le support de transmission pouvant être lu par un ordinateur peut transmettre des ondes porteuses ou des signaux (par exemple, une transmission de données câblée ou sans fil par Internet). Par ailleurs, des programmes fonctionnels, des codes, et des segments de code pour accomplir ou prendre en charge, lors d'une mise en oeuvre dans un matériel électronique approprié, certains éléments des revendications jointes peuvent être facilement interprétés par des programmeurs du métier concernés par les modes de réalisation. [0086] Les exemples de modes de réalisation décrits ci-dessus sont destinés à être illustratifs à tous les égards, plutôt que limitatifs, de la présente invention. Ainsi, la présente invention peut faire l'objet de nombreuses variantes de mises en oeuvre détaillées qui peuvent être déduites de la description contenue dans le présent document par un homme du métier. Aucun élément, acte, ou -41- 36-122/100306US instruction utilisé dans la description de la présente demande ne devrait être interprété comme critique ou essentiel pour l'invention, sauf spécification contraire. Par ailleurs, tel qu'utilisé ici, l'article « un » est destiné à comprendre un ou plusieurs éléments.

Claims (10)

  1. REVENDICATIONS1. Procédé pour éliminer sensiblement une influence de REVENDICATIONS1. Procédé pour éliminer sensiblement une influence de réflexions de multiples internes lors de la détermination d'une géographie sous la mer dans une zone géographique présentant un intérêt sans connaissance a priori des informations de sous-surface, le procédé comprenant : (a) la génération d'une série de signaux sismiques par une pluralité d'émetteurs de source ; (b) la réception de données brutes par une pluralité de récepteurs sur la base de la série générée de signaux sismiques ; (c) la création d'un ensemble de M fenêtres qui correspond physiquement à un espace au-dessous de la pluralité de récepteurs et qui comprend une zone géographique présentant un intérêt ; (d) l'attribution des données brutes reçues aux fenêtres respectives de l'ensemble de fenêtres sur la base d'un temps de réception des données brutes, pour générer M trames de données de fenêtre, dans lequel une première trame de données de fenêtre la plus haute incorpore des données brutes qui correspondent à un champ d'onde le plus haut le plus proche de la pluralité de récepteurs, et une Mième trame de données de fenêtre incorpore des données brutes qui correspondent à un champ d'onde le plus bas le plus éloigné de la pluralité de récepteurs ; (e) la génération de manière itérative d'un modèle de multiples internes, en utilisant un ensemble glissant de trois trames de données de fenêtre, dans lequel, pour chaque itération, le modèle de multiples internes comprend un premier produit de convolution de données provenant d'une première trame de données de fenêtre et d'une deuxième trame de données de fenêtre, et une corrélation de données provenant d'une troisième trame de données de fenêtre avec le premier produit ; (f) la sommation de tous les modèles de multiples internes générés de manière itérative pour créer un modèle de multiples internes complet, et la poursuite de la sommation jusqu'à ce que toutes les trames de données de fenêtre aient été utilisées ; et (g) la soustraction du modèle de multiples internes complet des données brutes pour éliminer sensiblement l'influence de multiples internes lors de la détermination de la géographie de la zone géographique présentant un intérêt.
  2. 2. Procédé selon la revendication 1, comprenant en outre : le traitement des données brutes pour supprimer des multiples associés à une surface avant la génération du modèle de multiples internes.
  3. 3. Procédé selon la revendication 2, dans lequel l'étape de traitement pour supprimer des multiples associés à une surface comprend : l'utilisation d'un procédé d'élimination de multiples associés à une surface pour supprimer les multiples de surface.
  4. 4. Procédé selon la revendication 1, dans lequel, en outre, l'étape de détermination de l'ensemble de M fenêtres est basée sur un temps de propagation de la série de signaux sismiques de la pluralité de sources jusqu'à la pluralité de récepteurs, et dans lequel, en outre, les M trames temporelles de fenêtre ont sensiblement une durée similaire.
  5. 5. Procédé selon la revendication 1, dans lequel l'étape de génération d'un modèle de multiples internes en utilisant des données de fenêtre provenant d'un premier ensemble de trois trames de données de fenêtre comprend l'application de l'équation suivante : -44-36-122/100306US dans lequel dans chaque ensemble défini de manière itérative de trois trames de données de fenêtre, un champ d'onde supérieur généré par des données dans la trame de données de fenêtre la plus haute est défini en tant que Pwi, un premier champ d'onde inférieur généré par des données dans la deuxième trame de données de fenêtre est défini en tant que Pwk, et un deuxième champ d'onde inférieur généré par des données dans la troisième trame de données de fenêtre est défini en tant que Pwb et en outre dans lequel, Pwi est un champ d'onde côté source qui représente une réflexion vers le bas d'un multiple interne réfléchi depuis la première fenêtre, Pwk est un champ d'onde côté source qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la deuxième fenêtre, et Pm est un champ d'onde côté récepteur qui représente une réflexion vers le haut d'un multiple interne réfléchi depuis la troisième fenêtre.
  6. 6. Procédé selon la revendication 5, dans lequel chacune des M fenêtres a des composantes de longueur et de profondeur, et dans lequel la composante de longueur est inférieure ou égale à une distance entre une première source et une dernière source, et en outre dans lequel la composante de profondeur est corrélée à un premier nombre d'échantillons qui est corrélé à une première profondeur en distance, et en outre dans lequel des fenêtres adjacentes se superposent d'un deuxième nombre d'échantillons inférieur au premier nombre d'échantillons, qui correspond à une superposition en profondeur définie en tant que deuxième profondeur, et encore en outre dans lequel -45-36-122/100306US la deuxième profondeur est inférieure à la première profondeur, et encore en outre dans lequel pour une valeur croissante de M, la profondeur de la fenêtre augmente.
  7. 7. Procédé selon la revendication 5, dans lequel chacun de la pluralité d'ensembles de fenêtres satisfait à une condition de monotonicité de pseudo profondeur des fenêtres inférieure-supérieure-inférieure, dans lequel Pwi est une fenêtre supérieure, et Pwk et Pwi sont toutes les deux des fenêtres inférieures.
  8. 8. Procédé pour déterminer des multiples internes dans des données sismiques marines sans nécessiter une connaissance a priori des sous-surfaces, le procédé comprenant : (a) l'obtention de données sismiques dans et concernant une zone géographique présentant un intérêt (GAI) ; (b) la suppression des multiples associés à une surface en relation avec la GAI ; (c) la séparation des données sismiques obtenues en une pluralité de fenêtres sur la base d'un temps de propagation des champs d'onde générés par une ou plusieurs sources de champ d'onde, et reçus par un ou plusieurs récepteurs de champ d'onde ; (d) l'exécution d'une convolution entre les données de fenêtre séparées qui sont responsables des réflexions vers le haut de multiples internes pour générer des données de convolution ; (e) l'exécution d'une corrélation entre les données de fenêtre séparées qui sont responsables des réflexions vers le bas des multiples internes correspondants ; (f) la détermination du produit des données corrélées et des données convoluées pour créer un modèle de multiples internes ; et (g) la répétition des étapes (d) à (f) pour une pluralité d'ensembles de la pluralité de fenêtres pour déterminer de ce fait un modèle de multiples internes sensiblement complet -46-36-122/100306US pour la GAI.
  9. 9. Système pour éliminer sensiblement une influence de réflexions de multiples internes lors de la détermination d'une géographie sous la mer dans une zone géographique présentant un intérêt sans connaissance a priori d'informations de sous-surface, le système comprenant : (a) une pluralité d'émetteurs de source configurés pour générer une série de signaux sismiques ; (b) une pluralité de récepteurs configurés pour recevoir des données brutes sur la base de la série générée de signaux sismiques ; et (c) un processeur configuré pour créer un ensemble de M fenêtres qui correspond physiquement à un espace au-dessous de la pluralité de récepteurs et qui comprend une zone géographique présentant un intérêt, attribuer des données brutes reçues fenêtres respectives de l'ensemble de fenêtres sur la base d'un temps de réception des données brutes, pour générer M trames de données de fenêtre, dans lequel une première trame de données de fenêtre la plus haute incorpore des données brutes qui correspondent à un champ d'onde le plus haut le plus proche de la pluralité de récepteurs, et une Mième trame de données de fenêtre incorpore des données brutes qui correspondent à un champ d'onde le plus bas le plus éloigné de la pluralité de récepteurs, générer de manière itérative un modèle de multiples internes, en utilisant un ensemble glissant de trois trames de données de fenêtre, dans lequel, pour chaque itération, le modèle de multiples internes comprend un premier produit de convolution de données d'une première trame de données de fenêtre et d'une deuxième trame de données de fenêtre, et une corrélation de données d'une troisième trame de données de fenêtre avec le premier produit, sommer la totalité des modèles de multiples internes générés de manière itérative -47-36-122/100306US pour créer un modèle de multiples internes complet, et poursuivre la sommation jusqu'à ce que toutes les trames de données de fenêtre aient été utilisées, et soustraire le modèle de multiples internes complet des données brutes pour éliminer sensiblement l'influence de multiples internes lors de la détermination de la géographie de la zone géographique présentant un intérêt.
  10. 10. Système pour déterminer des multiples internes dans des données sismiques marines sans nécessiter une connaissance a priori des sous-surfaces, le système comprenant : un processeur configuré pour obtenir des données sismiques dans et concernant une zone géographique présentant un intérêt (GAI), dans lequel le processeur est en outre configuré pour (a) supprimer des multiples associés à une surface en relation avec la GAI, (b) séparer les données sismiques obtenues en une pluralité de fenêtres sur la base d'un temps de propagation des champs d'onde générés par une ou plusieurs sources de champ d'onde, et reçus par un ou plusieurs récepteurs de champ d'onde, (c) effectuer une convolution entre les données de fenêtre séparées qui sont responsables des réflexions vers le haut de multiples internes pour générer des données de convolution, (d) effectuer une corrélation entre les données de fenêtre séparées qui sont responsables des réflexions vers le bas des multiples internes correspondants, (e) déterminer le produit des données corrélées et des données convoluées pour créer un modèle de multiples internes, et répéter les étapes (c) à (e) pour une pluralité d'ensembles de la pluralité de fenêtres, pour de ce fait déterminer un modèle de multiples internes sensiblement complet pour la GAI. -48-
FR1350277A 2012-01-12 2013-01-11 Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface Withdrawn FR2985831A1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US201261585830P 2012-01-12 2012-01-12

Publications (1)

Publication Number Publication Date
FR2985831A1 true FR2985831A1 (fr) 2013-07-19

Family

ID=47757848

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1350277A Withdrawn FR2985831A1 (fr) 2012-01-12 2013-01-11 Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface

Country Status (5)

Country Link
US (1) US9366775B2 (fr)
CA (1) CA2801734A1 (fr)
FR (1) FR2985831A1 (fr)
GB (1) GB2499882B (fr)
SG (1) SG192364A1 (fr)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2838857C (fr) * 2013-01-15 2021-10-26 Cgg Services Sa Traitement de donnees sismiques comprenant une attenuation multiple interne tridimensionnelle azimutale sans information souterraine
US9310504B2 (en) * 2013-02-01 2016-04-12 Pgs Geophysical As Systems and methods for detecting swell noise in a seismic gather
US10379245B2 (en) * 2013-07-03 2019-08-13 Pgs Geophysical As Method and system for efficient extrapolation of a combined source-and-receiver wavefield
CN105518489A (zh) * 2014-07-21 2016-04-20 杨顺伟 一种海上地震资料侧反射压制方法
CN114114419B (zh) * 2020-08-31 2024-04-09 中国石油化工股份有限公司 层间多次波预测、压制方法和提高地震资料成像品质方法
CN114428335A (zh) * 2020-09-28 2022-05-03 中国石油化工股份有限公司 一种省时层间多次波预测方法和压制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10295688B2 (en) * 2010-08-10 2019-05-21 Westerngeco L.L.C. Attenuating internal multiples from seismic data

Also Published As

Publication number Publication date
CA2801734A1 (fr) 2013-07-12
GB2499882B (en) 2017-12-06
SG192364A1 (en) 2013-08-30
US20130182535A1 (en) 2013-07-18
US9366775B2 (en) 2016-06-14
GB2499882A (en) 2013-09-04
GB201300498D0 (en) 2013-02-27

Similar Documents

Publication Publication Date Title
FR2991463B1 (fr) Systeme et procede de tomographie a haute definition et a haute resolution pour une utilisation pour la generation de modeles de vitesse et d'images de reflectivite
US20200183036A1 (en) Seismic data processing including true-azimuth three-dimensional internal multiple attentuation without subsurface information
FR2989788A1 (fr) Titre non renseigne.
US10310123B2 (en) Seismic reflection full waveform inversion for reflected seismic data
FR2985831A1 (fr) Traitement de donnees sismiques comprenant une modelisation de multiples internes sans informations de sous-surface
CA2733882C (fr) Procede de surveillance d'un site de stockage geologique de gaz par inversion stratigraphique de donnees sismiques
US9971051B2 (en) High-fidelity adaptive curvelet domain primary-multiple separation processing of seismic data
FR2990519A1 (fr) Procede, dispositif et algorithme de traitement pour l'elimination de multiples et de bruit de donnees sismiques marines
FR2989786A1 (fr) Dispositif et procede pour defantomer des donnees de flute a profondeur variable
FR3069929A1 (fr) Prediction en avant du trepan a l'aide de profil sismique vertical et de l'inversion globale
FR2985574A1 (fr) Dispositif et procede pour separer des donnees de tirs simultanes
CA2704141A1 (fr) Methode pour imager une zone cible du sous-sol a partir de donnees de type walkaway
FR2985040A1 (fr) Device and method for deghosting variable depth streamer data
FR2989787A1 (fr) Traitement de donnees sismiques comprenant une compensation des effets de fantome de source et de recepteur dans une migration inverse dans le temps
FR2985571A1 (fr) ESTIMATION CONJOINTE SIMULTANEE DE STATIQUES RESIDUELLES P-P et P-S
FR2904875A1 (fr) Traitement de donnees sismiques obtenues a partir de flutes superieures/inferieures et/ou de sources superieures/ inferieures.
FR2990769A1 (fr)
CA2845962A1 (fr) Systeme et procede de modelisation de multiples d'eaux peu profondes tridimensionnels au moyen de la reconstitution predictive du plancher marin
FR2992736A1 (fr) Trairement de donnees seismiques comprenant une prediction de multiples utilisant un retrait de migration inverse dans le temps
FR2916540A1 (fr) Procede d'exploration sismique permettant la supression de fantomes dus aux reflexions a la surface de l'eau, et procede de traitement de donnees sismiques pour la supression de ces fantomes
FR2794249A1 (fr) Procede de determination d'au moins un temps de parcours multiple sur un lever sismique
FR2845773A1 (fr) Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres
FR2972057A1 (fr) Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees
US20150355357A1 (en) Premigration deghosting for marine streamer data using a bootstrap approach in tau-p domain
FR2984525A1 (fr) Separation de champs d'onde pour des enregistrements sismiques repartis sur des surfaces d'enregistrement non planes

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 4

PLSC Publication of the preliminary search report

Effective date: 20160729

PLFP Fee payment

Year of fee payment: 5

ST Notification of lapse

Effective date: 20180928