FR2989787A1 - SEISMIC DATA PROCESSING COMPRISING COMPENSATION OF SOURCE AND RECEIVER FANTOME EFFECTS IN REVERSE TIME MIGRATION - Google Patents
SEISMIC DATA PROCESSING COMPRISING COMPENSATION OF SOURCE AND RECEIVER FANTOME EFFECTS IN REVERSE TIME MIGRATION Download PDFInfo
- Publication number
- FR2989787A1 FR2989787A1 FR1353488A FR1353488A FR2989787A1 FR 2989787 A1 FR2989787 A1 FR 2989787A1 FR 1353488 A FR1353488 A FR 1353488A FR 1353488 A FR1353488 A FR 1353488A FR 2989787 A1 FR2989787 A1 FR 2989787A1
- Authority
- FR
- France
- Prior art keywords
- rtm
- seismic
- acoustic wave
- source
- ghost
- 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
Links
- 230000005012 migration Effects 0.000 title claims abstract description 19
- 238000013508 migration Methods 0.000 title claims abstract description 19
- 230000000694 effects Effects 0.000 title abstract description 15
- 238000012545 processing Methods 0.000 title description 10
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000010304 firing Methods 0.000 claims description 6
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 230000000087 stabilizing effect Effects 0.000 claims description 2
- 230000006854 communication Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 9
- 239000010410 layer Substances 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000007175 bidirectional communication Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000276489 Merlangius merlangus Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/366—Seismic filtering by correlation of seismic signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/37—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy specially adapted for seismic systems using continuous agitation of the ground, e.g. using pulse compression of frequency swept signals for enhancement of received signals
- G01V1/375—Correlating received seismic signals with the emitted source signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Oceanography (AREA)
- Noise Elimination (AREA)
Abstract
Des procédés et des systèmes pour compenser des effets des fantômes de sources et de récepteurs dans une équation de migration inverse dans le temps (RTM) sont décrits. Des conditions de limite associées aux équations d'onde acoustique de RTM pour les champs d'ondes de source et enregistrés sont modifiées. Les équations d'onde acoustique de RTM modifiées résultantes sont résolues pour générer des images sismiques modélisées à fantômes compensés. Selon un autre aspect, une condition de formation d'image est également modifiée et les équations d'onde acoustique de RTM résultantes sont résolues pour générer des images de perturbation de vitesse et d'impédance.Methods and systems for compensating for effects of source and receiver ghosts in a time-reversal migration equation (RTM) are described. Boundary conditions associated with the RTM acoustic wave equations for source and recorded wave fields are modified. The resulting modified RTM acoustic wave equations are solved to generate compensated ghost modeled seismic images. In another aspect, an image forming condition is also modified and the resulting RTM acoustic wave equations are resolved to generate rate and impedance perturbation images.
Description
Traitement de données sismiques comprenant une compensation des effets de fantôme de source et de récepteur dans une migration inverse dans le temps DEMANDE CONNEXE [0001] La présente demande concerne et comporte une revendication de priorité fondée sur la demande de brevet US provisoire n° 61/635 379, déposée le 19 avril 2012, intitulée « Compensating for Source and Receiver Ghost Effects in Reverse Time Migration », de Yu Zhang, Graham Roberts et Lian Duan, dont la description est incorporée ici par voie de référence. DOMAINE TECHNIQUE [0002] Les modes de réalisation de l'objet présenté ici concernent généralement des procédés et des systèmes pour le traitement de données sismiques et, plus particulièrement, des mécanismes et des techniques pour améliorer une migration inverse dans le temps en compensant les effets de fantôme de source et de récepteur. CONTEXTE [0003] Pour des données sismiques marines, les fantômes de source et de récepteur générés par la surface libre provoquent une distorsion de fréquence et d'amplitude fonction de l'angle. Pour une interprétation et une inversion améliorées, Docket No. 100333/0336-156 2 ces effets non souhaités sont mieux corrigés dans une image de profondeur avant-empilage (« pre-stack depth image » en terminologie anglo-saxonne). [0004] L'interprétation de caractéristiques géologiques subtiles dans un environnement marin nécessite que les données sismiques associées contiennent à la fois des basses fréquences et des hautes fréquences nécessaires pour une formation d'image à haute résolution. Généralement, les données basse fréquence à haute fidélité permettent une pénétration plus grande des strates sous-jacentes et une meilleure illumination associée de cibles profondes en plus de fournir une stabilité plus grande et des informations plus importantes pour une inversion sismique. Par conséquent, il est préférable d'obtenir une image sismique à large bande passante. [0005] L'obtention d'images de profondeur avant empilage à large bande passante à partir de données sismiques marines nécessite de surmonter les limitations de bande passante imposées par les fantômes de source et de récepteur. Seismic Data Processing Comprising Source and Receiver Ghost Effects in a Time-Reversed Migration RELATED REQUEST [0001] The present application is directed to and includes a priority claim based on Provisional US Patent Application No. 61 / 635,379, filed April 19, 2012, entitled "Compensating for Source and Receiver Ghost Effects in Reverse Time Migration" by Yu Zhang, Graham Roberts and Lian Duan, the description of which is hereby incorporated by reference. TECHNICAL FIELD The embodiments of the object presented herein generally relate to methods and systems for the processing of seismic data and, more particularly, to mechanisms and techniques for improving a reverse migration over time by offsetting the effects. source and receiver phantom. BACKGROUND [0003] For marine seismic data, the source and receiver phantoms generated by the free surface cause frequency and angle-dependent distortion. For improved interpretation and inversion, these unwanted effects are better corrected in a pre-stack depth image. [0004] The interpretation of subtle geological features in a marine environment requires that the associated seismic data contain both low frequencies and high frequencies necessary for high resolution image formation. Generally, high fidelity low frequency data allows for greater penetration of the underlying layers and better associated deep target illumination in addition to providing greater stability and more important information for seismic inversion. Therefore, it is preferable to obtain a high bandwidth seismic image. Obtaining depth images before high bandwidth stacking from marine seismic data requires overcoming the bandwidth limitations imposed by source and receiver ghosts.
Récemment, de nombreuses tentatives ont été faites pour compenser des fantômes de sources et de récepteurs telles qu'une acquisition de flûte à profondeur variable comme décrit par R. Soubaras et P. Whiting (ci-après « Soubaras/Whiting ») dans leur article de 2011 intitulé « Variable Depth Streamer - The New Broadband Acquisition System », publié dans le 81st Annual International Meeting, SEG, Expanded Abstracts, pages 4349 à 4353 et incorporé ici par voie de référence. Ils proposaient de tirer avantage de la réponse à faible bruit de flûtes pleines de nouvelle génération et de la diversité d'encoches provoquée par une variation de profondeur des récepteurs pour obtenir un spectre large bande de qualité élevée. Recently, numerous attempts have been made to compensate for ghosts of sources and receivers such as variable depth flute acquisition as described by R. Soubaras and P. Whiting (hereinafter "Soubaras / Whiting") in their article. of 2011 entitled Variable Depth Streamer - The New Broadband Acquisition System, published in the 81st Annual International Meeting, SEG, Expanded Abstracts, pages 4349-4353 and incorporated herein by reference. They proposed to take advantage of the low noise response of next-generation solid flutes and the diversity of notches caused by a variation in receiver depth to obtain a high quality broadband spectrum.
Docket No. 100333/0336-156 3 [0006] Toujours avec des spectres large bande appropriés, les fantômes de récepteurs peuvent être retirés à un stade de traitement ultérieur par des techniques de traitement et de formation d'image telles que décrites par R. Soubaras (ci-après « Soubaras ») dans son article de 2010 intitulé « Deghosting by Joint Deconvolution of a Migration and a Mirror Migration », publié dans le 81st Annual International Meeting, SEG, Expanded Abstracts, pages 3406 à 3410, incorporé ici par voie de référence et par R. Soubaras et Y. Lafet dans leur article de 2011 intitulé « Variable Depth Streamer Acquisition : Broadband Data for Imaging and Inversion », publié dans le 81st Annual International Meeting, SEG, Expanded Abstracts, pages 2364 à 10 2368 et incorporé ici par voie de référence. Les techniques susmentionnées génèrent une ondelette exceptionnellement nette et propre pour une interprétation en relation avec des fantômes de récepteurs. [0007] Cependant, le retrait des fantômes de sources reste une tâche plus difficile basée sur la supposition que les sources sont excitées à une profondeur 15 approximativement fixe, c'est-à-dire, un manque de diversité d'encoches, et avec un grand incrément de tir dans les positions latérales, c'est-à-dire, un fort repliement de données (« data aliasing » en terminologie anglo-saxonne), comme cela apparaît généralement dans des études à flûtes remorquées marines. En outre, l'existence de fantômes déforme à la fois le spectre de migration et la relation entre l'amplitude et 20 l'angle (AVA). [0008] Par conséquent, il serait souhaitable de proposer des systèmes et des procédés qui évitent les problèmes et les inconvénients décrits précédemment et de Docket No. 100333/0336-156 4 réaliser une compensation à la fois des fantômes de sources et de récepteurs dans un ensemble de données sismiques à large bande passante. Docket No. 100333 / 0336-156 [0006] Still with appropriate broadband spectra, receptor ghosts can be removed at a later stage of processing by image processing and imaging techniques as described by R. Soubaras (hereinafter "Soubaras") in his 2010 article entitled "Deghosting by Joint Deconvolution of a Migration and a Mirror Migration", published in the 81st Annual International Meeting, SEG, Expanded Abstracts, pages 3406 to 3410, incorporated here by reference route and by R. Soubaras and Y. Lafet in their 2011 article "Variable Depth Streamer Acquisition: Broadband Data for Imaging and Inversion", published in the 81st Annual International Meeting, SEG, Expanded Abstracts, pages 2364 to 10 2368 and incorporated herein by reference. The aforementioned techniques generate an exceptionally clean and clean wavelet for interpretation in relation to receptor ghosts. [0007] However, removal of source ghosts remains a more difficult task based on the assumption that the sources are excited to an approximately fixed depth, i.e., lack of slot diversity, and with a large firing increment in the lateral positions, that is to say, a strong aliasing ("data aliasing" in English terminology), as it usually appears in studies with marine towed flutes. In addition, the existence of ghosts deforms both the migration spectrum and the relationship between amplitude and angle (AVA). Therefore, it would be desirable to provide systems and methods that avoid the problems and disadvantages described above and Docket No. 100333 / 0336-156 4 compensate for both source and receiver ghosts in a set of high bandwidth seismic data.
Docket No. 100333/0336-156 RESUME [0009] Selon un exemple de mode de réalisation, un procédé, mémorisé dans une mémoire et s'exécutant sur un processeur, pour compenser des fantômes dans 5 une image sismique calculée basée sur une migration inverse dans le temps (RTM) d'un ensemble de données sismiques associé à une ou plusieurs sources sismiques et à un ou plusieurs récepteurs sismiques, comprend la modification d'une première condition de limite d'une première équation d'onde acoustique de RTM associée à un champ d'ondes de source ; la modification d'une deuxième condition de limite i.o d'une deuxième équation d'onde acoustique de RTM associée à un champ d'ondes enregistré ; la résolution de la première équation d'onde acoustique de RTM et de la deuxième équation d'onde acoustique de RTM et la génération d'une image sismique calculée à fantômes compensés ; et la sortie de l'image sismique calculée à fantômes compensés. 15 [0010] Un noeud pour compenser des fantômes dans des données d'image sismiques basées sur une migration inverse dans le temps (RTM) comprend un ensemble de données sismiques ; un ou plusieurs processeurs configurés pour exécuter des instructions d'ordinateur et une mémoire configurée pour mémoriser lesdites instructions d'ordinateur, dans lequel lesdites instructions d'ordinateur 20 comprennent en outre un composant de condition de limite pour modifier des conditions de limite associées à des équations d'onde acoustique de RTM ; un composant de moteur pour résoudre les équations d'onde acoustique de RTM configuré avec des conditions de limite modifiées et générer une image sismique Docket No. 100333 / 0336-156 SUMMARY [0009] According to an exemplary embodiment, a method, stored in a memory and running on a processor, for compensating for ghosts in a computed seismic image based on a reverse migration time (RTM) of a set of seismic data associated with one or more seismic sources and one or more seismic receivers, includes modifying a first boundary condition of a first associated RTM acoustic wave equation to a source wave field; modifying a second boundary condition i.o of a second RTM acoustic wave equation associated with a recorded wave field; resolving the first RTM acoustic wave equation and the second RTM acoustic wave equation and generating a compensated ghost calculated seismic image; and the output of the seismic image computed with compensated ghosts. [0010] A node for compensating for ghosts in seismic image data based on reverse time migration (RTM) comprises a set of seismic data; one or more processors configured to execute computer instructions and a memory configured to store said computer instructions, wherein said computer instructions further comprise a boundary condition component for modifying boundary conditions associated with RTM acoustic wave equations; a motor component for solving the RTM acoustic wave equations configured with modified boundary conditions and generating a seismic image
calculée à fantômes compensés ; et un composant de sortie pour sortir l'image sismique calculée à fantômes compensés. BREVE DESCRIPTION DES DESSINS [0011] Les dessins joints, qui sont incorporés dans la spécification et qui constituent une partie de celle-ci, illustrent un ou plusieurs modes de réalisation et, avec la description, expliquent ces modes de réalisation. Sur les dessins : [0012] les figures 1 et 2 montrent divers aspects d'un système d'étude sismique marine dans lequel des modes de réalisation pour générer une image sismique calculée à fantômes compensés peuvent être mis en oeuvre ; [0013] les figures 3 à 6 montrent divers aspects de composants ou de modules logiciels qui peuvent être utilisés pour mettre en oeuvre les modes de réalisation ; [0014] la figure 7 est un organigramme représentant un procédé de 15 génération d'une image sismique calculée à fantômes compensés ; [0015] la figure 8 est une représentation graphique d'un enregistrement de tir ; [0016] la figure 9a est une représentation graphique d'un CIG dans le domaine angulaire migré, la figure 9b est une représentation graphique des courbes AVA associées et les figures 9c et 9d montrent les équivalents à fantômes 20 compensés ; et [0017] la figure 10 illustre un exemple de dispositif ou de système de traitement de données qui peut être utilisé pour mettre en oeuvre les modes de réalisation. calculated with compensated ghosts; and an output component for outputting the compensated ghost seismic image. BRIEF DESCRIPTION OF THE DRAWINGS [0011] The accompanying drawings, which are incorporated in the specification and constitute a part thereof, illustrate one or more embodiments and, with the description, explain these embodiments. In the drawings: Figures 1 and 2 show various aspects of a marine seismic survey system in which embodiments for generating a compensated ghost calculated seismic image can be implemented; Figures 3 to 6 show various aspects of components or software modules that can be used to implement the embodiments; FIG. 7 is a flow diagram showing a method of generating a compensated ghost seismic image; Figure 8 is a graphical representation of a firing record; FIG. 9a is a graphical representation of a CIG in the migrated angular domain, FIG. 9b is a graphical representation of the associated AVA curves, and FIGS. 9c and 9d show compensated ghost equivalents; and [0017] FIG. 10 illustrates an exemplary data processing device or system that can be used to implement the embodiments.
Docket No. 100333/0336-156 7 DESCRIPTION DETAILLEE [0018] La description qui suit des modes de réalisation fait référence aux dessins joints. Les mêmes numéros de référence sur les différents dessins identifient les mêmes éléments ou des éléments similaires. La description détaillée qui suit ne limite pas l'invention. Au lieu de cela, l'étendue de l'invention est définie par les revendications jointes. Certains des modes de réalisation qui suivent sont examinés, par souci de simplicité, en relation avec la terminologie et la structure de compensation des effets des fantômes de sources et de récepteurs dans une migration inverse dans le temps. Cependant, les modes de réalisation qui seront examinés ensuite ne sont pas limités à ces configurations, mais peuvent être étendus à d'autres agencements tels qu'examinés ultérieurement. [0019] 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 15 réalisation de l'objet présenté. 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, les fonctionnalités, structures ou caractéristiques particulières peuvent être combinées de n'importe quelle manière appropriée dans un ou plusieurs modes de réalisation. 20 [0020] Afin de fournir un certain contexte pour les modes de réalisation qui suivent associés à la collecte de données sismiques, au retrait des fantômes de sources et/ou de récepteurs et à la génération de cartes d'illumination sismique, Docket No. 100333/0336-156 8 examinons d'abord un processus et un système d'acquisition de données sismiques tels que décrits maintenant en relation avec les figures 1 et 2. 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 13 (« bird » en terminologie anglo-saxonne) 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 conformément à des communications bidirectionnelles que les oiseaux 13 peuvent recevoir du navire 2. [0021] Un ou plusieurs réseaux de sources 4a, b peuvent é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) et retournent 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 les récepteurs 14 prévus sur les flûtes 6. Les ondes sismiques sont ensuite réfléchies par la surface libre, c'est-à-dire, la surface de la masse d'eau (voir la figure 2, examinée ci- dessous), en se propageant vers le bas et sont une fois encore détectées par les récepteurs 14 prévus sur les flûtes 6 en tant que fantômes de récepteurs. Ce processus est généralement appelé « tir » d'une zone particulière des fonds 2 989 787 Docket No. 100333/0336-156 9 océaniques, la zone des fonds océaniques étant appelée « cellule » et la surface de la mer étant appelée « surface libre ». [0022] 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 5 plusieurs flûtes 6, qui sont composées de câbles 12, et d'une pluralité de récepteurs 14. Deux flûtes de sources, qui comprennent les sources 4a, b attachées aux câbles 12a, b respectifs, sont montrées sur la figure 2. Chaque source 4a, b est capable d'émettre une onde sonore ou un signal émis 20a, b respectif. A des fins de simplification des dessins, mais sans nuire en aucune manière à la compréhension des 10 nombreux principes impliqués, seul 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 provenant du premier signal émis 20a se propage vers le haut à partir du fond 15 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 à partir 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 20 partie du signal émis peut être réfractée (de nouveau selon la loi de Snell). [0023] 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 Docket No. 100333/0336-156 10 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 à un deuxième point de réfraction/réflexion 28a, un deuxième ensemble de signaux réfracté et réfléchi 32a et 30a, est généré par la suite. En outre, comme montré sur la figure 2, il semble y avoir 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, les signaux réfracté et réfléchi sont générés par le i.o 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. [0024] Les signaux enregistrés par les récepteurs sismiques 14 varient dans le temps, comportant des pics d'énergie qui peuvent correspondre à des réflecteurs entre 15 des 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 à de 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 ne subissent qu'une seule réflexion à partir d'une interface entre les couches de la sous-surface (par 20 exemple, le premier signal réfléchi 24a). Les ondes autres que les ondes principales sont connues en tant que multiples (ou fantômes). Un signal de multiple de surface est un tel exemple d'un multiple, cependant, il existe d'autres manières de générer des multiples. Par exemple, des réflexions par la surface libre peuvent se propager en 2989787 Docket No. 100333/0336-156 11 arrière jusqu'aux récepteurs et être enregistrées en tant que fantômes. Les multiples (et fantômes) n'ajoutent pas d'informations utiles à la géologie au-dessous du fond océanique, et ainsi ils sont, par essence, un bruit, et il est souhaitable de les éliminer et/ou de réduire sensiblement et/ou d'éliminer leur influence dans un traitement de 5 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. [0025] Une fois collectées, les données sismiques sont traitées pour générer des images qui peuvent être utilisées par les hommes du métier pour identifier des caractéristiques de sous-surface. Le traitement de données sismiques peut comprendre l'application de nombreuses techniques différentes à l'ensemble de données. Par exemple, en relation avec des fantômes, d'abord, supposons qu'une profondeur de source soit désignée par Az, et qu'une vitesse d'onde acoustique dans l'eau soit désignée par vo . Le fantôme de source généré par la réflexion de surface libre est un effet dépendant de l'angle, qui modifie à la fois l'amplitude de l'ondelette et le spectre. En tant que simplification, supposons que la réflectivité de la surface soit -1. A l'emplacement de la source, le fantôme de source G, agissant sur le champ d'ondes sismiques avec un angle de propagation a, adopte une forme dans le domaine fréquentiel telle que : iwcosa,Az, scocosa,4z, COS (1) G, (co,a,)= e v" -e sin vo o osasAz, =-2ism vo - 2989787 Docket No. 100333/0336-156 12 En outre, dans les domaines des fréquences et des nombres d'onde, une relation existe telle que : cos as = 1- V° (k2 k2 (2) 02 x y) où kx et ky sont les nombres d'onde spatiaux horizontaux dans les directions x et y, 5 respectivement. De manière similaire, le fantôme de récepteur peut être exprimé par : wcosarAz,. ro cos ar4z, co cos arAzr = -2i sin Gr (co, = e v1, (3) -e v' vo [0026] La théorie et la mise en oeuvre de la migration inverse dans le temps (RTM, « reverse time migration » en terminologie anglo-saxonne) vraie ont été développées et présentées à la fois par Y. Zhang et J. Sun dans leur article de 2009 10 intitulé « Practical Issues of Reverse Time Migration : True-Amplitude Gathers, Noise Removal and Harmonic-Source Encoding », publié dans First Break, Vol. 26, pages 19 à 25, incorporé ici par voie de référence et par S. Xu, Y. Zhang et B. Tang (ci-après « XU et d'autres ») dans leur article de 2011 intitulé « 3D Angle Gathers from Reverse Time Migration », publié dans Geophysics, Vol. 76, n° 2, pages S77 à S92, incorporé 15 ici par voie de référence. [0027] La migration d'un enregistrement de tir Q(x, y; x' ys;t), avec le tir à (x' y ' zs = 0) et les récepteurs à (x, y, z = 0) , nécessite le calcul des champs d'ondes provenant de l'emplacement des sources et observés aux emplacements des récepteurs. Sur la base du champ d'ondes de sources s'étendant dans le temps et du 2989787 Docket No. 100333/0336-156 13 champ d'ondes enregistré calculé en arrière dans le temps, le champ d'ondes de sources est désigné par p F et le champ d'ondes enregistré est désigné par pB dans les équations : ' 1 a2 pF(ic;t;:i s) = 0, 0; t; -.Z, ) = 8(.Z - )r f (P)de , (4) v2 at2 -A { P F(X , Y, z = 5 et { [ 12 2t 2 A PB(i;t;i,)=- 0, (5) ) p B(x, y, z = 0;t;.X>s)= Q(x, y; xs, Ys; t), où y = v(.X) est la vitesse, f(t) est la signature de la source et Aa= 2x ± a y + z2 0 est l'opérateur laplacien. [0028] En outre, on peut montrer que l'équation de propagation d'onde (4) et 10 l'équation de propagation d'onde (5), lorsqu'elles sont combinées avec la condition de formation d'image tridimensionnelle : R(X> ; 0; ç9) v(ï)8(6r-0)8(ço'-ço)pBp,dtcrisdO'dço' sin 6P où 8 est l'angle de réflexion à l'emplacement de formation d'image et ço est l'angle d'azimut à l'emplacement de formation d'image, fournissent une amplitude de migration 15 respectant la relation AVA dans le domaine angulaire de sous-surface comme décrit (6) 2989787 Docket No. 100333/0336-156 14 par XU et d'autres. Il conviendrait de noter que la théorie de la RTM décrite précédemment ignore les effets des fantômes. [0029] En examinant maintenant un mode de réalisation, sur la base de la propagation d'un champ d'ondes artificiel en RTM, il est simple de compenser le 5 fantôme de source pendant la migration. L'exécution de la compensation de fantôme de source dans un mode de réalisation implique de modifier l'équation (4) en : i9, (x, Y, z = 0; co; ->., ) = g(8(.x- i., ) f (co) (7) i coG s(co, a s) où j),(co) représente la transformation de Fourier de p F(t) . De manière similaire, le même concept peut être appliqué dans un mode de réalisation pour une compensation 10 de fantôme de récepteur avec l'équation : Ys , y; xs s ; co) 19^ 8(x, y, z = 0; co;.is) = Q(x, . Gr(co, a r) Cependant, il conviendrait de noter que, dans le mode de réalisation, lorsque des données sismiques sont repliées, l'application de l'équation de compensation de fantôme de récepteur (8) peut provoquer directement une instabilité, ainsi un remède 15 de régularisation peut être appliqué pendant une propagation d'onde. [0030] Selon un mode de réalisation, la RTM peut être utilisée pour estimer des perturbations de vitesse et/ou d'impédance. Examinons l'équation acoustique avec à la fois des variations de vitesse v(z) et de densité p(5e) telle que : (8) Docket No. 100333/0336-156 a2 at2 v2 pV -V p(i;t;_i,)= 8(i> - 5c> s)(5(t).p) Pour des modèles de vitesse vo (i>) et de densité p0 (.x) initiales donnés, le champ d'ondes perturbé 8p(i'; t;.i>s) = p - po satisfait à l'équation suivante : 15 (9) 2 a 1 PoV V v(2) at 2 po j 28v a2V sp\ - V Sp(i;t;i's)= vo' at2 Po ) po (-i;t;.is), (10) où 81, = v- vo et Sp=p-po désignent des perturbations de vitesse et de densité, respectivement. En utilisant le procédé développé dans le document de Jin et d'autres (1992) et le document de Forgus et Lambare (1997), une relation basée sur les rayons entre des modèles géologiques perturbés et un champ d'ondes peut être obtenue par : 8v 2 (5p Sn' vo 2 ', COS as COS Ot - + cos 0- -= 3271- cos U rAA e-m) Q.5(81-0)cr cr cicoch9 (11) xr xs vo Po sin 0' y, ) Vo ) s r où A, (Ar) est l'amplitude de la fonction de Green de la source (le récepteur) au point d'image, r, (r,.) est le temps de propagation entre la source (le récepteur) et le point d'image. Dans le contexte de la RTM, l'équation (11) peut être reformulée en modifiant la condition de formation d'image exposée dans l'équation (6) comme suit : 2 81) 8(vp) = f III cos t9" 8(6" 6) 1)B(6°)15 F(6°) (12) sin u - + cos- sdcod 0'. sin 9' ico L'équation (12) indique que, si des rassemblements angulaires de sous-surface sont délivrés avec une condition de formation d'image correcte, des images d'angle proche prédisent la perturbation d'impédance (8(vp)I(vp)), tandis que les images d'angle 2 9 89 78 7 Docket No. 100333/0336-156 16 lointain peuvent être utilisées pour estimer la perturbation de vitesse (8v/v). Par conséquent, ces modes de réalisation séparent les effets de la vitesse et de la densité sur l'image empilée en délivrant des CIG de domaine angulaire. [0031] En examinant maintenant la figure 3, un schéma d'un noeud de mode de 5 réalisation 300 pour compenser des effets de fantômes de sources et de récepteurs dans des équations d'onde acoustique de RTM est représenté. Le noeud 300 comprend un composant de condition de limite 302, un composant de moteur 304, un composant de sortie 306 et un ensemble de données sismiques 308. En poursuivant avec le mode de réalisation, le composant de condition de limite 302 fournit la capacité 10 de modifier la condition de limite des équations de champ d'ondes acoustiques de RTM comme décrit précédemment. Il conviendrait de noter que, dans le mode de réalisation, la modification de la condition de limite est basée sur une transformation de Fourier des champs d'ondes de sources et de récepteurs. En outre, dans le mode de réalisation, il conviendrait de noter que la modification de la condition de limite des 15 champs d'ondes de récepteurs de l'équation de RTM peut introduire une instabilité dans l'équation de RTM, nécessitant une stabilisation. [0032] Ensuite, dans le mode de réalisation du noeud 300, le composant de moteur 304 fournit la capacité de résoudre les équations de RTM modifiées à condition de limite en tenant compte de l'ensemble de données sismiques 308 sélectionné. Il conviendrait de noter que, dans le mode de réalisation du noeud 300, des réponses en amplitude et spectrale fiables de RTM sont générées après la compensation des effets des fantômes. Ensuite, dans le mode de réalisation du noeud 300, le composant de Docket No. 100333/0336-156 17 sortie 306 fournit la capacité de délivrer une image de profondeur sismique à fantômes compensés à partir de la RTM pour une analyse associée à l'exploration des ressources. Il conviendrait de noter que dans le mode de réalisation avec les réponses en amplitude et spectrale à fantômes compensés de la RTM, la variation géologique basse fréquence est mieux décrite dans l'image sismique calculée de sortie. [0033] Ensuite, dans le mode de réalisation, en faisant référence à la figure 4, un schéma 400 d'un composant de moteur 304 d'un noeud 300 est représenté. Le composant de moteur 304 comprend en outre un composant de stabilisateur 402 pour stabiliser une équation d'onde acoustique de RTM associée à des données sismiques repliées. Il conviendrait de noter que, dans le mode de réalisation, la stabilisation est basée sur un remède de régularisation et est appliquée pendant la propagation d'onde. [0034] En examinant maintenant la figure 5, un schéma d'un mode de réalisation d'un noeud 500 pour compenser les effets de fantômes de sources et de récepteurs dans les équations d'onde acoustique de RTM et pour générer des images 15 de perturbation de vitesse et d'impédance est montré. Le noeud 500 comprend un composant de condition de limite 302, un composant de moteur 304, un composant de sortie 306 et un ensemble de données sismiques 308 comme décrit précédemment. Le noeud 500 comprend en outre un composant de condition de formation d'image 502 pour fournir la capacité de modifier la condition de limite de formation d'image des 20 équations de champ d'ondes acoustiques de RTM comme décrit précédemment. [0035] En poursuivant avec le mode de réalisation, sur la figure 6, un schéma 600 d'un composant de moteur 304 d'un noeud 500 est représenté. Le composant de Docket No. 100333/0336-156 18 moteur 304 comprend en outre un composant de perturbation de vitesse et d'impédance 602 pour estimer des perturbations de vitesse et d'impédance associées à l'ensemble de données sismiques 308 sur la base de la résolution des équations de champ d'ondes acoustiques de RTM. Il conviendrait de noter que, dans le mode de réalisation, les compensations de perturbations de vitesse et d'impédance et de fantômes peuvent avoir lieu simultanément. [0036] En examinant maintenant la figure 7, un mode de réalisation de procédé 700 de compensation des fantômes dans une image sismique calculée basée sur la RTM d'un ensemble de données sismiques associé à une ou plusieurs sources sismiques et à un ou plusieurs récepteurs sismiques est représenté. En commençant à l'étape 702, le mode de réalisation de procédé 700 modifie une première condition de limite d'une première équation d'onde acoustique de RTM. Selon un autre aspect de l'étape 702 du mode de réalisation de procédé 700, l'équation d'onde acoustique est associée à un champ d'ondes de sources. Ensuite, à l'étape 704, l'exemple de mode de réalisation de procédé 700 modifie une deuxième condition de limite d'une deuxième équation d'onde acoustique de RTM. Selon un autre aspect du mode de réalisation de procédé 700, à l'étape 704, l'équation d'onde acoustique est associée à un champ d'onde de récepteur. [0037] En poursuivant à l'étape 706 du mode de réalisation de procédé 700, les équations d'onde acoustique de RTM à condition de limite modifiée sont résolues. Selon un autre aspect de l'étape 706 du mode de réalisation de procédé 700, une image sismique calculée à fantômes compensés est générée. Ensuite, à l'étape 708 du Docket No. 100333/0336-156 19 mode de réalisation de procédé 700, l'image sismique calculée à fantômes compensés associée à l'ensemble de données sismiques est délivrée. [0038] Un exemple est fourni pour montrer les réponses en amplitude et spectrale fiables de la RTM après une compensation des effets des fantômes selon les s modes de réalisation précédents. Il est supposé pour cet exemple que les fantômes à la fois de sources et de récepteurs sont enregistrés. La figure 8 montre un enregistrement de tir unique bidimensionnel sur cinq réflecteurs horizontaux à différentes profondeurs, avec le tir au centre de la section et les récepteurs avec un décalage de 7500 mètres de chaque côté. Les profondeurs du tir et des récepteurs 10 sont de 10 mètres et de 15 mètres, respectivement, et la vitesse dans l'eau est de 1500 mètres par seconde. Dans la modélisation, il est supposé que la réflectivité est uniforme au niveau de tous les points de réflexion selon tous les angles de réflexion. Du fait de l'existence de fantômes, l'amplitude d'ondelette et le spectre sont tous les deux déformés avec le temps de propagation et la distance latérale, en plus de l'effet 15 de dispersion géométrique. En utilisant l'équation (4) et l'équation (5) de formulation de la RTM à amplitude vraie classique, en empilant tous les rassemblements de tirs d'images communes migrées pour générer des rassemblements décalés de sous-surface et ensuite en les convertissant en rassemblements d'images communes (CIG) de domaine angulaire de sous-surface, comme montré sur la figure 9, le résultat est 20 une distorsion dans le spectre de l'image migrée (figure 9a) et une tendance de l'AVA incorrecte (figure 9b). Docket No. 100333 / 0336-156 DETAILED DESCRIPTION [0018] The following description of the embodiments refers to the accompanying drawings. The same reference numbers on the different drawings identify the same or similar elements. The detailed description which follows does not limit the invention. Instead, the scope of the invention is defined by the appended claims. Some of the following embodiments are discussed, for the sake of simplicity, in connection with the terminology and compensation structure of the effects of source and receiver ghosts in a reverse time migration. However, the embodiments that will be discussed next are not limited to these configurations, but may be extended to other arrangements as discussed later. [0019] A reference throughout the specification to "an embodiment" means that a particular feature, structure, or feature described in connection with an embodiment is included in at least one embodiment of the object. present. Thus, the appearance of the phrase "in one embodiment" at various locations throughout the specification does not necessarily refer to the same embodiment. In addition, the particular features, structures or features may be combined in any suitable manner in one or more embodiments. In order to provide some context for the following embodiments associated with seismic data collection, removal of source and / or receiver ghosts and generation of seismic illumination maps, Docket No. 100333 / 0336-156 8 First, consider a seismic data acquisition process and system as described now in connection with Figures 1 and 2. In Figure 1, a data acquisition system 10 includes a ship 2 towing several flutes 6 which can extend for kilometers behind the ship 2. Each of the flutes 6 can comprise one or more birds 13 ("bird" in English terminology) which keep the flute 6 in a fixed position known by relative to the other flutes 6, and the birds 13 are able to move the flute 6 as desired in accordance with bi-directional communications that the birds 13 can receive from the ship 2. [0021] One or more each source network 4a, b may also be towed by the ship 2 or another vessel to generate seismic waves. The source networks 4a, b may be placed either in front of or behind the receivers 14, or both behind and in front of the receivers 14. The seismic waves generated by the source networks 4a, b propagate downwards, are reflected by the oceanic bottoms and into them, in which the refracted waves are finally reflected by one or more reflection structures (not shown in Figure 1) and return to the surface (see Figure 2, discussed below) . The reflected seismic waves propagate upwards and are detected by the receivers 14 provided on the flutes 6. The seismic waves are then reflected by the free surface, that is to say, the surface of the body of water ( see Figure 2, discussed below), propagating downward and are again detected by receivers 14 provided on flutes 6 as receiver ghosts. This process is generally referred to as "firing" from a particular area of the ocean floor, the ocean floor area being referred to as the "cell" and the surface of the sea being called "free surface". ". FIG. 2 illustrates a side view of the data acquisition system 10 of FIG. 1. The ship 2, situated on the surface of the ocean 46, tows one or more flutes 6, which are composed of cables 12, and a plurality of receivers 14. Two source streamers, which include the sources 4a, b attached to the respective cables 12a, b, are shown in FIG. 2. Each source 4a, b is capable of emitting a wave sound or a transmitted signal 20a, b respectively. For purposes of simplifying the drawings, but without in any way impairing the understanding of the many principles involved, only a first emitted signal 20a will be shown (even if some or all of the sources 4 may simultaneously transmit (or not) signals issued similar). The first transmitted signal 20a propagates through the ocean 40 and arrives at a first refraction / reflection point 22a. The first reflected signal 24a from the first transmitted signal 20a propagates upward from the ocean floor 42, back to the receivers 14. As can be appreciated by those skilled in the art, whenever an optical signal or acoustic - propagates from a medium with a first refractive index n1 and encounters a different medium, with a second refractive index n2, a part of the emitted signal is reflected at an angle equal to the angle of incidence (according to well known Snell's law), and a second part of the emitted signal can be refracted (again according to Snell's law). Thus, as shown in FIG. 2, the first emitted signal 20a generates a first reflected signal 24a and a first refracted signal 26a. The first refracted Dact No. 100333 / 0336-156 signal 26a propagates through a layer of sediment 16 (which may be generically referred to as the first subsurface layer 16) below the ocean floor 42, and can now be considered as a "new" emitted signal, so that when it encounters a second medium at a second refraction / reflection point 28a, a second set of refracted and reflected signals 32a and 30a is generated thereafter. In addition, as shown in Figure 2, there appears to be a large hydrocarbon deposit 44 in a third medium, or a solid soil / rock layer 18 (which may be generically referred to as the second subsurface layer 18). Accordingly, the refracted and reflected signals are generated by the hydrocarbon deposit, and it is the object of the data acquisition system 10 to generate data that can be used to discover these hydrocarbon deposits 44. The signals recorded by the seismic receivers 14 vary over time, including energy peaks which may correspond to reflectors between layers. In reality, since the ocean floor and air / water are highly reflective, some of the peaks correspond to multiple reflections or parasitic reflections that should be eliminated before the geophysical structure image can be properly formed. . The main waves undergo only a single reflection from an interface between the sub-surface layers (e.g., the first reflected signal 24a). Non-main waves are known as multiples (or ghosts). A surface multiple signal is such an example of a multiple, however, there are other ways to generate multiples. For example, reflections from the free surface can propagate back to the receivers and be recorded as ghosts. Multiples (and ghosts) do not add useful information to the geology below the ocean floor, and so they are, in essence, a noise, and it is desirable to eliminate them and / or substantially reduce and / or to eliminate their influence in a signal processing of other reflected signals so as to correctly establish the presence (or absence) of underground / underwater hydrocarbon deposits. Once collected, the seismic data is processed to generate images that can be used by those skilled in the art to identify sub-surface features. The seismic data processing may include the application of many different techniques to the data set. For example, in relation to ghosts, first, suppose that a source depth is designated by Az, and that an acoustic wave velocity in water is designated by vo. The source phantom generated by the free surface reflection is an angle dependent effect, which modifies both the amplitude of the wavelet and the spectrum. As a simplification, suppose the reflectivity of the surface is -1. At the location of the source, the source ghost G, acting on the seismic wave field with a propagation angle α, adopts a frequency-domain shape such that: iwcosa, Az, scocosa, 4z, COS (1 ) G, (co, a,) = ev "-e sin vo o osasAz, = -2ism vo - 2989787 Docket No. 100333 / 0336-156 12 In addition, in the fields of frequencies and wave numbers, a The relationship exists such that cos = V ° (k2 k2 (2) 02 xy) where kx and ky are the horizontal spatial wave numbers in the x and y directions, respectively. receiver can be expressed by: wcosarAz, ro cos ar4z, co cos arAzr = -2i sin Gr (co, = e v1, (3) -ev 'vo [0026] The theory and implementation of inverse migration in Real Time (RTM) were developed and presented by both Y. Zhang and J. Sun in their 2009 article 10 entitled "Practical Issues of Reverse Time Migration". : True-Amplitude Gathers, Noise Removal and Harmonic-Source Encoding ", published in First Break, Vol. 26, pages 19-25, incorporated herein by reference and by S. Xu, Y. Zhang and B. Tang (hereinafter "XU and others") in their 2011 article entitled "3D Angle Gathers from Reverse". Time Migration ", published in Geophysics, Vol. 76, No. 2, pages S77-S92, incorporated herein by reference. The migration of a firing record Q (x, y; x 'ys; t), with the firing at (x' y 'zs = 0) and the receivers at (x, y, z = 0) , requires the calculation of the wave fields from the location of the sources and observed at the locations of the receivers. On the basis of the time-based source wave field and the recorded back-calculated wave field in time, the source wave field is designated by p F and the recorded wave field is designated pB in the equations: ## EQU1 ## t; -Z) = 8 (.Z -) rf (P) of, (4) v2 at2 -A {PF (X, Y, z = 5 and {[12 2t 2 A PB (i; t; i, ) = - 0, (5)) p B (x, y, z = 0; t; .X> s) = Q (x, y; xs, Ys; t), where y = v (.X) is velocity, f (t) is the signature of the source and Aa = 2x ± ay + z2 0 is the Laplacian operator. In addition, it can be shown that the wave propagation equation (4) and the wave propagation equation (5), when combined with the three-dimensional image formation condition: R (X>; 0; ç9) v (ï) 8 (6r-0) 8 (ço'-ço) pBp, where δ is the angle of reflection at the formation site of image and δ is the azimuth angle at the imaging location, provide a migration amplitude respecting the AVA relation in the sub-surface angular domain as described (6) 2989787 Docket No. 100333/0336 -156 14 by XU and others. It should be noted that the RTM theory previously described ignores the effects of ghosts. [0029] Turning now to an embodiment, based on the propagation of an artificial RTM wave field, it is simple to compensate for the source ghost during migration. The execution of source phantom compensation in one embodiment involves modifying equation (4) at: i9, (x, Y, z = 0; co; ->.) = G (8 (. where (j), (co) represents the Fourier transform of p F (t). Similarly, the same concept can be applied in one embodiment for receiver ghost compensation with the equation: Ys, y; xs s; co) 19 ^ 8 (x, y, z = 0; co; .is) = Q (x, Gr (co, ar) However, it should be noted that in the embodiment, when seismic data is When folded, the application of the receiver ghost compensation equation (8) can directly cause instability, so a regulation remedy may be applied during wave propagation. [0030] According to one embodiment, the RTM can be used to estimate velocity and / or impedance perturbations Let's consider the acoustic equation with both v (z) and p (5e) density variations such that: (8) Docket No. 100333 / 0336-156 a2 at2 v2 pV -V p (i; t; _i,) = 8 (i> - 5c> s) (5 (t) .p) For velocity models vo (i>) and of given density p0 (.x), the disturbed wave field 8p (i '; t; i> s) = p - po satisfies the following equation: (9) 2a 1 PoV V v (2) ) at 2 in j 28v a2V sp \ - V Sp (i; t; i's) = vo 'at2 Po) po (-i; t; .is), (10) where 81, = v-vo and Sp = p -po designate ent perturbation of speed and density, respectively. Using the method developed in Jin et al. (1992) and Forgus and Lambare (1997), a ray-based relationship between disturbed geologic models and a wavefield can be obtained by: 8v 2 (5p Sn 'vo 2', COS as COS Ot - + cos 0- - = 3271- cos U rAA em) Q.5 (81-0) cr cr cicoch9 (11) xr xs vo Po sin 0 'y ,) Where A, (Ar) is the magnitude of the Green function of the source (the receiver) at the image point, r, (r ,.) is the propagation time between the source (the receiver) and the image point. In the context of the RTM, equation (11) can be reformulated by modifying the imaging condition set forth in equation (6) as follows: 2 81) 8 (vp) = f III cos t9 "8 (6 "6) 1) B (6 °) 15 F (6 °) (12) sin u - + cos-sdcod 0 '. sin 9 'ico Equation (12) indicates that, if sub-surface angular gatherings are delivered with a proper imaging condition, near-angle images predict the impedance disturbance (8 (vp) I (vp)), while far-angle images can be used to estimate velocity disturbance (8v / v). Therefore, these embodiments separate the effects of velocity and density on the stacked image by delivering angular domain CIGs. [0031] Referring now to FIG. 3, a diagram of embodiment node 300 for compensating for ghost effects of sources and receivers in RTM acoustic wave equations is shown. Node 300 includes a boundary condition component 302, a motor component 304, an output component 306 and a seismic data set 308. Continuing with the embodiment, the boundary condition component 302 provides the capacitance 10. to modify the boundary condition of the acoustic wave field equations of RTM as previously described. It should be noted that, in the embodiment, the modification of the boundary condition is based on a Fourier transformation of source and receiver wave fields. Further, in the embodiment, it should be noted that changing the boundary condition of the receiver wavefields of the RTM equation can introduce instability into the RTM equation, necessitating stabilization. Next, in the embodiment of the node 300, the engine component 304 provides the ability to resolve the modified RTM equations with boundary condition taking into account the selected seismic data set 308. It should be noted that in the embodiment of the node 300, reliable RTM amplitude and spectral responses are generated after ghost effect compensation. Then, in the embodiment of the node 300, the Docket component No. 100333 / 0336-156 output 306 provides the ability to output a compensated ghost seismic depth image from the RTM for analysis associated with the resource exploration. It should be noted that in the embodiment with compensated ghost amplitude and spectral responses of the RTM, the low frequency geological variation is best described in the calculated output seismic image. Then, in the embodiment, with reference to Figure 4, a diagram 400 of a motor component 304 of a node 300 is shown. The engine component 304 further comprises a stabilizer component 402 for stabilizing an RTM acoustic wave equation associated with folded seismic data. It should be noted that, in the embodiment, the stabilization is based on a regularization remedy and is applied during the wave propagation. [0034] By now examining FIG. 5, a diagram of an embodiment of a node 500 for compensating for the ghost effects of sources and receivers in the RTM acoustic wave equations and for generating images of disturbance of speed and impedance is shown. Node 500 includes a boundary condition component 302, a motor component 304, an output component 306, and a seismic data set 308 as previously described. Node 500 further comprises an image forming condition component 502 for providing the ability to modify the image forming boundary condition of the RTM acoustic wave field equations as previously described. Continuing with the embodiment, in Figure 6, a diagram 600 of a motor component 304 of a node 500 is shown. The Docket component No. 100333 / 0336-156 engine 304 further comprises a speed and impedance disturbance component 602 for estimating rate and impedance disturbances associated with the seismic data set 308 based on of the resolution of the acoustic wave field equations of RTM. It should be noted that in the embodiment, speed and impedance and ghost disturbance compensations can occur simultaneously. Referring now to FIG. 7, a ghost compensation method embodiment 700 in a computed seismic image based on the RTM of a seismic data set associated with one or more seismic sources and one or more receivers. seismic is represented. Beginning at step 702, the method embodiment 700 modifies a first boundary condition of a first RTM acoustic wave equation. In another aspect of step 702 of method embodiment 700, the acoustic wave equation is associated with a source wave field. Then, in step 704, the method embodiment example 700 modifies a second limit condition of a second RTM acoustic wave equation. According to another aspect of the method embodiment 700, in step 704, the acoustic wave equation is associated with a receiver wave field. Continuing with step 706 of the method embodiment 700, the RTM acoustic wave equations with the modified limit condition are resolved. According to another aspect of step 706 of the method embodiment 700, a compensated ghost calculated seismic image is generated. Next, at step 708 of Docket No. 100333 / 0336-156 19 method embodiment 700, the compensated ghost calculated seismic image associated with the seismic data set is outputted. An example is provided to show the reliable amplitude and spectral responses of the RTM after compensation for ghost effects according to the previous embodiments. It is assumed for this example that ghosts of both sources and receivers are recorded. Figure 8 shows a two-dimensional single shot record on five horizontal reflectors at different depths, with shooting at the center of the section and receivers with an offset of 7500 meters on each side. The depths of the shot and the receivers 10 are 10 meters and 15 meters, respectively, and the speed in the water is 1500 meters per second. In modeling, it is assumed that reflectivity is uniform at all reflection points at all angles of reflection. Because of the existence of ghosts, the wavelet amplitude and the spectrum are both distorted with propagation time and lateral distance, in addition to the geometric dispersion effect. Using Equation (4) and Formula Equation (5) for Formulation of True True Amplitude RTM, stacking all migrated common image shootings to generate shifted sub-surface shifts and then converting to subsurface angular domain common image (CIG) gatherings, as shown in FIG. 9, the result is a distortion in the spectrum of the migrated image (FIG. 9a) and a trend of the AVA incorrect (Figure 9b).
Docket No. 100333/0336-156 20 [0039] Après la compensation des effets des fantômes à la fois du côté des sources et du côté des récepteurs, les ondelettes sur le rassemblement angulaire migré ont une largeur de bande de fréquence plus large et plus équilibrée et apparaissent beaucoup plus nettes (figure 9c). Par ailleurs, les amplitudes crêtes normalisées le long des réflecteurs dans le domaine angulaire de réflexion convergent bien, ce qui indique que la réflectivité est bien récupérée et que la relation AVA est plus fiable (figure 9d). [0040] Le dispositif informatique ou les dispositifs informatiques ou d'autres noeuds de réseau impliqués dans une prédiction d'image sismique modélisée à fantômes compensés telle qu'exposée dans les modes de réalisation décrits ci-dessus peuvent être n'importe quel type de dispositif informatique capable de traiter et de communiquer les données sismiques associées à une étude sismique. Un exemple d'un système informatique représentatif capable d'exécuter des opérations selon ces modes de réalisation est illustré sur la figure 10. Le système 200 comprend, entre autres éléments, un serveur 201, une interface sources/récepteurs 202, un bus de données/communication interne (bus) 204, un ou des processeurs 208 (les hommes du métier peuvent apprécier que, dans les systèmes de serveur modernes, un traitement en parallèle est de plus en plus répandu et, alors qu'un processeur unique aurait é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) et, par conséquent, il pourrait y avoir plusieurs processeurs, agissant en série et/ou en parallèle, selon les besoins de l'application spécifique), un port de bus Docket No. 100333/0336-156 21 série universel (USB) 210, un lecteur à lecture/écriture (RM/) de disque compact (CD)/de disque vidéo 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. [0041] 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 est devenu de plus en plus répandu, peut comprendre des dispositifs de stockage de masse de type clé USB 224, entre autres types), un dispositif ou des dispositifs de mémoire morte 218 (ceux-ci peuvent comprendre des dispositifs de mémoire ROM programmable électriquement effaçable (EE) (EEPROM), des dispositifs PROM effaçables par ultraviolets (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 des disques CD/DVD 234 (qui peuvent être à la fois à lecture et à écriture) peuvent être utilisés avec un dispositif R/W CD/DVD 212. Les disquettes 237 peuvent être utilisées avec un dispositif de lecteur 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 l'intégralité, d'un code de programmation de logiciel exécutable (logiciel) 236 qui peut mettre en oeuvre une partie ou la totalité des parties du 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 pas d'une manière limitative, un support de mise en mémoire RAM 220) qui peut mémoriser la totalité ou certains des composants du logiciel 236. Docket No. 100333 / 0336-156 [0039] After compensating for the effects of ghosts on both the source and the receiver side, the wavelets on the migrated angular congregation have a wider and wider frequency bandwidth. balanced and appear much clearer (Figure 9c). On the other hand, normalized peak amplitudes along reflectors in the angular reflection range converge well, indicating that the reflectivity is well recovered and that the AVA relationship is more reliable (Fig. 9d). The computing device or the computing devices or other network nodes involved in a predicted ghost modeled seismic image prediction as set forth in the embodiments described above may be any type of computer device capable of processing and communicating the seismic data associated with a seismic survey. An example of a representative computer system capable of executing operations according to these embodiments is illustrated in FIG. 10. The system 200 comprises, among other elements, a server 201, a source / receiver interface 202, a data bus / internal communication (bus) 204, one or more processors 208 (those skilled in the art can appreciate that in modern server systems, parallel processing is becoming more common and, while a single processor would have been used in the past to implement many or at least several functions, it is more common now to have a single dedicated processor for certain functions (eg, digital signal processors) and, therefore, there could be several processors, acting in series and / or in parallel, according to the needs of the specific application), a bus port Docket No. 100333 / 0336-156 21 universal serial (USB) 210, a Read / write (RM /) compact disc (CD) / digital video disc (DVD) 212, a floppy disk drive 214 (although less used today, many servers still include this device), and a 232. The data storage unit 232 itself may comprise a hard disk drive (HDD) 216 (this may include a conventional magnetic storage medium, but as has become moreover in addition, may include USB type mass storage devices 224, among other types), a device or ROM devices 218 (these may include electrically erasable programmable (EE) ROM devices ( EEPROM), ultraviolet erasable PROM (UVPROM) devices, among other types), and RAM devices 220. USB memory device 224 may be used with USB 210, and CD / DVD discs 234 ( which can be both read and write) can be used with an R / W CD / DVD 212. The floppy disks 237 can be used with a floppy disk device 214. Each of the memory devices, or storage media (216, 218, 220, 224, 234 and 237, among other types), may contain parts or components, or all, of an executable (software) software programming code which can implement some or all of the parts of the process described herein. In addition, the processor 208 itself may contain one type or different types of storage devices (most likely, but not limitatively, a RAM storage medium 220) that can store all or some of the memory devices. software components 236.
Docket No. 100333/0336-156 22 [0042] En plus des composants décrits ci-dessus, le système 200 comprend également une 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 à diodes électroluminescentes (DEL), des afficheurs à plasma, des tubes à rayons cathodiques (CRT), entre autres. La console d'utilisateur 235 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 d'intercommunication interactifs. [0043] La console d'utilisateur 235, et ses composants s'ils sont prévus 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 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 développés ultérieurement. 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 est connectée électriquement au moins une antenne 240 (selon un exemple de mode de réalisation, ce serait une antenne GPS de réception uniquement, et au moins une antenne de communication bidirectionnelle de satellite séparée). Le système 200 peut accéder à Internet 242, soit par l'intermédiaire d'une connexion câblée, soit par l'intermédiaire de l'interface Docket No. 100333/0336-156 23 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. [0044] 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 faire partie d'une configuration de réseau plus grande comme dans un réseau mondial (GAN) (par exemple, Internet 242), qui permet finalement une connexion à diverses lignes terrestres. [0045] Selon un autre exemple de mode de réalisation, le système 200, conçu pour une utilisation dans l'exploration sismique, s'interfacera avec une ou plusieurs sources 4a, b et un ou plusieurs récepteurs 14. Ceux-ci, comme décrit précédemment, sont attachés aux flûtes 6a, b, auxquelles sont également attachés les oiseaux 13a, b qui sont utiles pour maintenir un positionnement. Comme examiné davantage précédemment, les sources 4 et les récepteurs 14 peuvent communiquer avec le serveur 201 soit par l'intermédiaire d'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 (décrits collectivement en tant que conduit de communications 246). [0046] Selon d'autres exemples de modes de réalisation, la console d'utilisateur 235 fournit des moyens pour qu'un personnel entre des commandes et une configuration dans le système 200 (par exemple, par l'intermédiaire d'un clavier, de boutons, de commutateurs, d'un écran tactile et/ou d'une manette). Le dispositif Docket No. 100333/0336-156 24 d'affichage 226 peut être utilisé pour montrer : la position de la flûte 6 ; des représentations visuelles des 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. L'unité d'interface des sources et des récepteurs 202 peut recevoir des données sismiques d'hydrophone du récepteur 14, par l'intermédiaire du conduit de communication de flûte 246 (examiné ci-dessus) qui peut faire 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 des commandes peuvent également être envoyées aux oiseaux 13 pour maintenir un positionnement correct de la flûte. L'unité d'interface des sources et des 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 246 entre le système 200 et la source 4. [0047] Le bus 204 fournit un trajet de données pour des éléments tels que : 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 à un afficheur 226 ; ou pour que l'utilisateur envoie des commandes à des programmes/logiciel de système d'exploitation 236 qui pourraient résider soit dans le processeur 208, soit dans l'unité d'interface des sources et des récepteurs 202. Docket No. 100333 / 0336-156 [0042] In addition to the components described above, the system 200 also includes a user console 234, which may include a keypad 228, a display 226 and a mouse 230. The entire of these components are known to those skilled in the art, and this description includes all known and future variants of these types of devices. The display 226 may be any type of known display or presentation screen, such as liquid crystal displays (LCDs), light-emitting diode (LED) displays, plasma displays, test tubes. cathode rays (CRT), among others. The user console 235 may include one or more user interface mechanisms such as a mouse, keyboard, microphone, touchpad, touch screen, voice recognition system, among other interactive communication devices. The user console 235, and its components if provided separately, interface with the server 201 through an input / output interface (I / O) server 222, which may be an RS232, Ethernet, USB or other communication port, or may include all or some of these, and further includes any other type of communication means, currently known or developed later . The system 200 may further include a communications satellite / global positioning system (GPS) 238 transmission-reception device, to which at least one antenna 240 is electrically connected (according to an exemplary embodiment, it would be an antenna Receiving GPS only, and at least one separate satellite bidirectional communication antenna). The system 200 can access the Internet 242, either via a wired connection, or via the input / output DIT interface No. 100333 / 0336-156 23 directly, or by a link wirelessly via antenna 240 and transceiver 238. Server 201 may be coupled to other computing devices, such as those that implement or control the ship's equipment. 2, via one or more networks. The server 201 may be part of a larger network configuration as in a global network (GAN) (e.g., Internet 242), which ultimately allows connection to various landlines. According to another exemplary embodiment, the system 200, designed for use in seismic exploration, will interface with one or more sources 4a, b and one or more receivers 14. These, as described previously, are attached to the flutes 6a, b, which are also attached to the birds 13a, b which are useful for maintaining a positioning. As discussed further above, sources 4 and receivers 14 may communicate with server 201 either via an electrical cable that is part of flute 6, or through a wireless system that can communicate via antenna 240 and transceiver 238 (collectively described as communications path 246). According to other exemplary embodiments, the user console 235 provides means for a personnel to enter commands and a configuration in the system 200 (for example, via a keyboard, buttons, switches, a touch screen and / or a joystick). Display device Docket No. 100333 / 0336-156 24 can be used to show: the position of flute 6; visual representations of acquired data; state information of the source 4 and the receiver 14; study information; and other important information for the seismic data acquisition process. The source and receiver interface unit 202 may receive hydrophone seismic data from the receiver 14, via the flute communication channel 246 (discussed above) which may be part of the flute 6, as well as position information of the flute 6 from the birds 13; the link is bidirectional so that commands can also be sent to the birds 13 to maintain proper flute positioning. The source and receiver interface unit 202 may also bi-directionally communicate with the sources 4 via the flute communication conduit 248 which may be part of the flute 6. Excitation signals, signals control, output signals and status information associated with the source 4 may be exchanged by the flute communication path 246 between the system 200 and the source 4. [0047] The bus 204 provides a data path for elements such as: the transfer and storage of data that originate either from the source sensors or from the flute receivers; for the processor 208 to access stored data contained in the data storage unit memory 232; for the processor 208 to send information for a visual display to a display 226; or for the user to send commands to operating system programs / software 236 that could reside either in the processor 208 or in the source and receiver interface unit 202.
Docket No. 100333/0336-156 25 [0048] Le système 200 peut être utilisé pour mettre en oeuvre les procédés décrits ci-dessus associés à une prédiction d'image sismique modélisée à fantômes compensés 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 exécuter 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 à multiples supports tels que les dispositifs 216, 218, 220, 224, 234 et/ou 237 (décrits ci-dessus) ou d'autres formes de supports capables de mémoriser de manière portable des informations (par exemple, la clé USB 426). Ces supports de mémorisation peuvent être insérés dans des dispositifs tels que le lecteur de CD-ROM 414, le lecteur de disque 412, entre autres types de dispositifs de mémorisation de logiciel, et lus par ceux-ci. [0049] Les exemples de modes de réalisation présentés fournissent un noeud de serveur et un procédé pour une prédiction d'image sismique calculée à fantômes compensés associée à des ensembles de données sismiques. On devrait comprendre que cette description n'est pas destinée à limiter l'invention. Au contraire, les exemples de modes de réalisation sont destinés à couvrir les variantes, les modifications et les équivalents qui sont inclus dans l'esprit et l'étendue de l'invention. En outre, dans la description détaillée des exemples de modes de réalisation, de nombreux détails spécifiques sont exposés afin de fournir une compréhension détaillée de l'invention. Cependant, un homme du métier comprendrait que divers modes de réalisation peuvent être mis en pratique sans ces détails spécifiques. Docket No. 100333 / 0336-156 [0048] The system 200 may be used to implement the methods described above associated with predicted ghost modeled seismic image prediction according to an exemplary embodiment. Hardware, firmware, software, or a combination of these can be used to perform the various steps and operations described here. According to an exemplary embodiment, the software 236 for performing the steps discussed above may be stored and distributed on multi-media storage devices such as devices 216, 218, 220, 224, 234 and / or 237 ( described above) or other forms of media capable of wearable storage of information (e.g., USB key 426). These storage media can be inserted into devices such as the CD-ROM drive 414, the disk drive 412, among other types of software storage devices, and read by them. The exemplary embodiments presented provide a server node and a method for a computed ghosted seismic image prediction associated with seismic data sets. It should be understood that this description is not intended to limit the invention. In contrast, the exemplary embodiments are intended to cover variants, modifications, and equivalents that are included in the spirit and scope of the invention. In addition, in the detailed description of the exemplary embodiments, many specific details are set forth in order to provide a detailed understanding of the invention. However, one skilled in the art would understand that various embodiments can be practiced without these specific details.
Docket No. 100333/0336-156 26 [0050] Bien que les caractéristiques et les éléments des présents exemples de modes de réalisation aient été décrits dans les modes de réalisation dans des combinaisons particulières, chaque caractéristique ou élément peut être utilisé seul sans les autres caractéristiques et éléments des modes de réalisation ou dans diverses combinaisons avec ou sans autres caractéristiques et éléments présentés ici. Les procédés ou les organigrammes fournis dans la présente demande peuvent être mis en oeuvre dans un programme, un logiciel, ou un micrologiciel d'ordinateur mis en oeuvre de manière tangible sur un support de mémorisation pouvant être lu par un ordinateur pour une exécution par un ordinateur d'usage général ou un processeur. [0051] Cette description écrite utilise des exemples de l'objet présenté pour permettre à n'importe quel homme du métier de mettre en pratique le susdit, y compris en fabriquant et en utilisant n'importe quels dispositifs ou systèmes et en exécutant n'importe quels procédés incorporés. L'étendue brevetable de l'objet est définie par les revendications et peut comprendre d'autres exemples qui apparaissent aux hommes du métier. Ces autres exemples sont destinés à être dans l'étendue des revendications. Docket No. 100333 / 0336-156 [0050] Although the features and elements of the present exemplary embodiments have been described in the embodiments in particular combinations, each feature or element may be used alone without the others. features and elements of the embodiments or in various combinations with or without other features and elements presented herein. The methods or flowcharts provided in the present application may be implemented in a computer program, software, or firmware implemented in a tangible manner on a computer-readable storage medium for execution by a computer. general purpose computer or processor. This written description uses examples of the object presented to enable any person skilled in the art to practice the aforesaid, including manufacturing and using any devices or systems and executing n ' any incorporated processes. The patentable scope of the object is defined by the claims and may include other examples that occur to those skilled in the art. These other examples are intended to be within the scope of the claims.
Claims (10)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261635379P | 2012-04-19 | 2012-04-19 |
Publications (1)
Publication Number | Publication Date |
---|---|
FR2989787A1 true FR2989787A1 (en) | 2013-10-25 |
Family
ID=48537512
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1353488A Withdrawn FR2989787A1 (en) | 2012-04-19 | 2013-04-17 | SEISMIC DATA PROCESSING COMPRISING COMPENSATION OF SOURCE AND RECEIVER FANTOME EFFECTS IN REVERSE TIME MIGRATION |
Country Status (4)
Country | Link |
---|---|
US (1) | US9465125B2 (en) |
FR (1) | FR2989787A1 (en) |
GB (1) | GB2503541A (en) |
SG (1) | SG194307A1 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140379266A1 (en) * | 2013-06-25 | 2014-12-25 | Westerngeco L.L.C. | Processing survey data containing ghost data |
US10379245B2 (en) * | 2013-07-03 | 2019-08-13 | Pgs Geophysical As | Method and system for efficient extrapolation of a combined source-and-receiver wavefield |
US20150276956A1 (en) * | 2014-03-28 | 2015-10-01 | Cgg Services Sa | Wave-equation based processing and analysis of imaged seismic data |
US10339234B2 (en) | 2014-11-25 | 2019-07-02 | Cgg Services Sas | Method for elastic model perturbation estimation from reverse time migration |
US10670751B2 (en) * | 2015-03-27 | 2020-06-02 | Cgg Services Sas | Full waveform inversion method for seismic data processing using preserved amplitude reverse time migration |
EP3387465B1 (en) * | 2015-12-11 | 2021-08-04 | ION Geophysical Corporation | System and method for reconstructed wavefield inversion |
EP3559707B1 (en) | 2016-12-20 | 2021-02-17 | ION Geophysical Corporation | System and method for reconstructed wavefield imaging |
US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100302906A1 (en) * | 2009-05-28 | 2010-12-02 | Chevron U.S.A. Inc. | Method for wavefield-based data processing including utilizing multiples to determine subsurface characteristics of a suburface region |
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5999488A (en) | 1998-04-27 | 1999-12-07 | Phillips Petroleum Company | Method and apparatus for migration by finite differences |
US8472674B2 (en) * | 2007-02-06 | 2013-06-25 | Naum Marmal Yevskyy | Method of surface seismic imaging using both reflected and transmitted waves |
US20100118651A1 (en) | 2008-11-10 | 2010-05-13 | Chevron U.S.A. Inc. | Method for generation of images related to a subsurface region of interest |
FR2955396B1 (en) | 2010-01-15 | 2013-03-01 | Cggveritas Services Sa | DEVICE FOR PROCESSING SEISMIC MARINE DATA |
FR2955397B1 (en) | 2010-01-15 | 2012-03-02 | Cggveritas Services Sa | METHOD AND DEVICE FOR ACQUIRING MARINE SEISMIC DATA |
US9658353B2 (en) | 2010-06-17 | 2017-05-23 | Westerngeco L.L.C. | Regulating coherent boundary reflections during generation of a modeled wavefield |
US8385151B2 (en) | 2010-06-24 | 2013-02-26 | Chevron U.S.A. Inc. | Reverse time migration with absorbing and random boundaries |
-
2013
- 2013-04-10 US US13/859,807 patent/US9465125B2/en active Active
- 2013-04-16 SG SG2013028766A patent/SG194307A1/en unknown
- 2013-04-17 FR FR1353488A patent/FR2989787A1/en not_active Withdrawn
- 2013-04-19 GB GB1307126.1A patent/GB2503541A/en not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100302906A1 (en) * | 2009-05-28 | 2010-12-02 | Chevron U.S.A. Inc. | Method for wavefield-based data processing including utilizing multiples to determine subsurface characteristics of a suburface region |
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
Non-Patent Citations (2)
Title |
---|
SHENG XU ET AL: "Antileakage Fourier transform for seismic data regularization", GEOPHYSICS, vol. 70, no. 4, July 2005 (2005-07-01), pages V87 - V95, XP001231046 * |
YIKE LIU ET AL: "Reverse time migration of multiples for subsalt imaging", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 76, no. 5, September 2011 (2011-09-01), pages WB209 - WB216, XP001571023, ISSN: 0016-8033, [retrieved on 20111121], DOI: 10.1190/GEO2010-0312.1 * |
Also Published As
Publication number | Publication date |
---|---|
GB2503541A (en) | 2014-01-01 |
GB201307126D0 (en) | 2013-05-29 |
SG194307A1 (en) | 2013-11-29 |
US20130343154A1 (en) | 2013-12-26 |
US9465125B2 (en) | 2016-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2989788A1 (en) | TITLE NOT RENTED. | |
FR2989787A1 (en) | SEISMIC DATA PROCESSING COMPRISING COMPENSATION OF SOURCE AND RECEIVER FANTOME EFFECTS IN REVERSE TIME MIGRATION | |
FR2991463B1 (en) | SYSTEM AND METHOD FOR HIGH DEFINITION AND HIGH RESOLUTION TOMOGRAPHY FOR USE IN GENERATING SPEED MODELS AND REFLECTIVITY IMAGES | |
US9971051B2 (en) | High-fidelity adaptive curvelet domain primary-multiple separation processing of seismic data | |
FR2974636B1 (en) | DEVICE AND METHOD FOR ADAPTIVE ANGULAR EXTRAPOLATION FOR THREE-DIMENSIONAL REVERSE MIGRATION ANGULAR GAPS IN TIME | |
FR2985040A1 (en) | DEVICE AND METHOD FOR DEGHOSTING VARIABLE DEPTH STREAMER DATA | |
FR2961316A1 (en) | PROCESS FOR PROCESSING SEISMIC MARINE DATA | |
FR2955396A1 (en) | DEVICE FOR PROCESSING SEISMIC MARINE DATA | |
US9594180B2 (en) | Removing ghost reflections from marine seismic data | |
FR2990519A1 (en) | METHOD, DEVICE AND PROCESS ALGORITHM FOR REMOVING MULTIPLE AND NOISE FROM SEISMIC MARINE DATA | |
FR2916540A1 (en) | SEISMIC EXPLORATION METHOD FOR THE REMOVAL OF GHOSTS DUE TO WATER SURFACE REFLECTIONS, AND METHOD OF PROCESSING SEISMIC DATA FOR THE SUPRESSION OF THESE GHOSTS | |
FR2972056A1 (en) | DEVICE AND METHOD FOR MULTI-DIMENSIONAL DATA COHERENCE CONTROLLED NOISE REMOVAL | |
FR2971059A1 (en) | DEVICE AND METHOD FOR DETERMINING WAVE ATTENUATION IN A CONDITION NEAR THE SURFACE | |
FR2990769A1 (en) | ||
CA2704141A1 (en) | Method for generating an image of an underground target area using walkaway data | |
FR2992736A1 (en) | SEISMIC DATA TRAINING COMPRISING MULTIPLE PREDICTION USING TIME REVERSE MIGRATION WITHDRAWAL | |
US20140172307A1 (en) | Target-oriented 4d binning in common reflection point | |
FR3001301A1 (en) | ||
FR2972057A1 (en) | DEVICE AND METHOD FOR MULTI-DIMENSIONAL DATA COHERENCE CONTROLLED NOISE REMOVAL | |
FR2985831A1 (en) | SEISMIC DATA PROCESSING COMPRISING MULTIPLE INTERNAL MODELING WITHOUT SUB-SURFACE INFORMATION | |
FR2971859A1 (en) | METHOD AND DEVICE FOR SMOOTHING WITH PRESERVATION OF THE TRAVEL TIME | |
US10324208B2 (en) | Premigration deghosting for marine streamer data using a bootstrap approach in Tau-P domain | |
US9671513B2 (en) | Seismic image acquisition using ghost compensation in Kirchhoff migration | |
FR2966253A1 (en) | Method for deghosting marine seismic data related to subsurface of body of slanted streamer, involves generating final image of subsurface in processor based on deghosted reflectivity of joint deconvoluting images | |
US20140198612A1 (en) | Ghost compensation in beam migration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20171222 |
|
ST | Notification of lapse |
Effective date: 20171229 |