FR2979711A1 - Procede de traitement de donnees sismiques marines - Google Patents
Procede de traitement de donnees sismiques marines Download PDFInfo
- Publication number
- FR2979711A1 FR2979711A1 FR1260128A FR1260128A FR2979711A1 FR 2979711 A1 FR2979711 A1 FR 2979711A1 FR 1260128 A FR1260128 A FR 1260128A FR 1260128 A FR1260128 A FR 1260128A FR 2979711 A1 FR2979711 A1 FR 2979711A1
- Authority
- FR
- France
- Prior art keywords
- image
- migration
- water
- seismic
- reflectivity
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012545 processing Methods 0.000 title description 13
- 230000008569 process Effects 0.000 title description 5
- 230000005012 migration Effects 0.000 claims abstract description 67
- 238000013508 migration Methods 0.000 claims abstract description 67
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 36
- 238000002310 reflectometry Methods 0.000 claims abstract description 24
- 230000001902 propagating effect Effects 0.000 claims abstract description 4
- 238000012546 transfer Methods 0.000 claims description 14
- 230000001364 causal effect Effects 0.000 description 10
- 230000008901 benefit Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 7
- 230000000630 rising effect Effects 0.000 description 7
- 238000003672 processing method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 230000001427 coherent effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 125000001475 halogen functional group Chemical group 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002745 absorbent Effects 0.000 description 1
- 239000002250 absorbent Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012421 spiking Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000002123 temporal effect 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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/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
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- 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
- G01V1/3817—Positioning of seismic devices
- G01V1/3826—Positioning of seismic devices dynamic steering, e.g. by paravanes or birds
-
- 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
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- 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
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)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Procédé pour supprimer des fantômes de données sismiques relatives à un sous-sol d'une masse d'eau, le procédé comprenant des étapes consistant à : - appliquer une procédure de migration à des données d'ondes enregistrées par des détecteurs qui sont remorqués par un navire, les ondes se propageant du sous-sol vers les détecteurs, pour obtenir une première image du sous-sol ; - appliquer une procédure de migration miroir aux données d'ondes enregistrées par les détecteurs pour obtenir une deuxième image du sous-sol ; - déconvoluer conjointement la première image et la deuxième image pour supprimer les fantômes d'une réflectivité du sous-sol ; et obtenir une image déghostée du sous-sol sur la base de la réflectivité de laquelle les fantômes ont été supprimés.
Description
PROCEDE DE TRAITEMENT DE DONNEES SISMIQUES MARINES DOMAINE TECHNIQUE La présente invention concerne le traitement de données sismiques marines acquises au moyen d'au moins une flûte sismique remorquée par un bateau et comportant une pluralité de récepteurs. Plus particulièrement l'invention concerne un procédé de traitement adapté au cas où certaines au moins des données sont acquises à des profondeurs d'eau différentes, ceci pouvant être réalisé au moyen de flûtes ayant au moins une portion inclinée par rapport à la surface de l'eau, ou de flûtes horizontales situées à des profondeurs différentes.
ETAT DE LA TECHNIQUE Une technique répandue de recherche de pétrole ou de gaz consiste à effectuer une prospection sismique du sous-sol. Pour imager la structure du sous-sol, le géophysicien utilise les techniques de « sismique réflexion ».
En sismique marine, la technique consiste à remorquer derrière un navire : - une ou plusieurs source(s) d'énergie pour l'émission d'une onde acoustique, et - des récepteurs sismiques disposés sur des flûtes (ou « streamers » selon la terminologie anglo-saxonne) positionnées horizontalement à une profondeur constante (de l'ordre de Dz = 7.5 mètres) pour enregistrer l'onde acoustique réfléchie par les interfaces entre couches géologiques. La source communique une onde acoustique à l'eau, en créant un champ d'ondes (ondes de compression) qui se propage de façon cohérente vers le bas (propagation descendante). Lorsque le champ d'ondes frappe les interfaces entre les formations terrestres, appelées réflecteurs, il est réfléchi à travers la terre et l'eau jusqu'aux récepteurs sismiques (propagation montante), où il est converti en signaux électriques et enregistré.
Les récepteurs sismiques sont disposés de façon telle et en nombre suffisant pour que les signaux enregistrés, appelés traces, constituent des données sismiques et permettent de reconstituer la configuration des couches géologiques.
Un problème que l'on rencontre est celui de la réverbération et peut être expliqué comme suit. Une onde sismique réfléchie par un réflecteur passe dans l'eau dans une direction généralement ascendante. Cette onde, que l'on appelle la « primaire », se propage dans l'eau et passe par le récepteur sismique qui enregistre sa présence.
Le champ d'ondes continue sa progression vers la surface de l'eau (qui a un coefficient de réflexion égal à -1), où il est réfléchi vers le bas. Ce champ d'ondes réfléchi ou « fantôme » se propage également dans l'eau et passe par les récepteurs où il est à nouveau enregistré avec une polarité inverse et un retard Dt qui vaut, pour des ondes se propageant verticalement : Dt=2 Dz/c, Avec : - Dt : la différence de temps entre les enregistrements par le récepteur de l'onde primaire et du fantôme, - Dz la distance entre la flûte et la surface de l'eau, - c la vitesse de propagation de l'onde dans l'eau (à savoir 1500m/s).. Cette réverbération du champ d'ondes sismique dans l'eau obscurcit les données sismiques, en amplifiant certaines fréquences et en en atténuant d'autres, ce qui rend difficile l'analyse des formations terrestres sous-jacentes. Le fantôme correspond, dans le domaine spectral, à un filtre de fonction de transfert : G(f) = 1 - exp(2jrrfAt), Cette fonction de transfert G(f) est nulle pour les fréquences f multiples de 1 c/2 750 fn, avec f, = -At = -Az - Ces fréquences pour lesquelles une fonction de transfert est nulle sont appelées « notch ». Les notchs sont particulièrement gênants car il n'est pas possible de les déconvoluer. L'usage est donc en sismique de placer les flûtes à une profondeur telle que le premier notch soit en dehors de la bande utile. 100 Hz étant une borne haute de la bande passante sismique, cela conduit à placer les flûtes à une profondeur d'environ 7,5 m. Mais cette profondeur relativement faible des flûtes a pour effet de pénaliser l'enregistrement des basses fréquences (pour les basses fréquences, la fonction G(f) est proportionelle à la profondeur Dz). Les documents US 4 353 121 et US 4 992 992 décrivent des procédés de traitement permettant de supprimer les fantômes à partir de données sismiques enregistrées en utilisant un dispositif d'acquisition comprenant une flûte sismique présentant un angle (de l'ordre de 2 degrés) avec la surface de l'eau. Dans cette configuration, c'est au cours de l'opération de sommation des données que la suppression des fantômes est effectuée. En effet les données acquises sont redondantes, et le procédé de traitement comporte une étape de sommation (ou « stack » selon la terminologie anglo-saxonne) pour obtenir l'image finale du sous-sol à partir des données redondantes. C'est au cours de l'étape de sommation que se fait la suppression des fantômes, car les enregistrements contribuant à cette somme, ayant été enregistrés par des récepteurs différents, présentent des notchs à des fréquences différentes, de sorte que l'information manquante du fait de la présence d'un notch sur un récepteur sismique est obtenue à partir d'un autre récepteur. Le document US 4 353 121 décrit un procédé de traitement basé sur les étapes suivantes : - mise en collection point miroir (« common depth point » selon la terminologie anglo-saxonne), - extrapolation 1D (une dimension) sur une surface horizontale (ou « datuming » selon la terminologie anglo-saxonne), - correction de type NMO (pour « Normale MoveOut » selon la terminologie anglo-saxonne), - sommation (ou « stack » selon la terminologie anglo-saxonne), Le « datuming » est un procédé de traitement où, à partir des données de N récepteurs sismiques R, de positions (x',z,) avec n=1,..N, on synthétise les données qu'auraient enregistrées les récepteurs sismiques s'ils avaient été placés aux mêmes positions horizontales xi, mais à une même profondeur de référence zo constante pour tous les récepteurs sismiques. Le « datuming » est dit 1D si l'on fait l'hypothèse que les ondes sismiques se propagent verticalement. Dans ce cas, le procédé se résume à appliquer à chaque enregistrement temporel enregistré par un récepteur sismique donné, un retard (ou « static shift » selon la terminologie anglo-saxonne) correspondant au temps de propagation vertical entre la profondeur réelle zi, du récepteur R, et la profondeur de référence zo. Par ailleurs, le brevet US 4,353,121 décrit le procédé consistant à obtenir une sommation de primaire (« primary-stack ») en utilisant une correction NMO alignant les réflexions, puis une sommation de fantôme (« ghost-stack ») en alignant les réflexions fantôme, puis à les combiner pour obtenir une image après stack au rapport signal-à-bruit renforcé. Tout comme US 4 353 121, le document US 4 992 992 propose de reconstituer les données sismiques qui auraient été enregistrées par un câble horizontal à partir des données sismiques enregistrées par le câble présentant un angle avec la surface de l'eau. Toutefois, le document US 4 992 992 vise à prendre en compte les propagations non verticales, en remplaçant l'étape de « datuming » 1D de US 15 4 353 121 par un « datuming » 2D. Ce « datuming » 2D prend en compte le fait que la propagation des ondes n'est pas forcément verticale, contrairement à ce qui est supposé dans le cas d'un « datuming » 1D tel que proposé dans US 4 353 121. Plus spécifiquement, US 4 992 992 décrit un procédé de traitement qui, à 20 partir des données enregistrées par les N récepteurs R, d'indice n=1,..N se trouvant à une distance horizontale xi, de la source et à une profondeur zn,=z0+xntga, (a correspondant à l'angle entre le câble et la surface de l'eau, et tg la fonction trigonométrique tangente), reconstruit les données sismiques qui auraient été enregistrées par des récepteurs sismiques situés aux mêmes 25 positions horizontales xi, mais à une même profondeur de référence zo constante pour tous les récepteurs sismiques. Pour ce faire, deux jeux de données sismiques reconstruites comme si elles avaient été enregistrées par une flûte horizontale sont calculés puis sommés après multiplication par un facteur. 30 Le premier jeu de données est synthétisé en supposant que les ondes sismiques se propagent vers le haut comme les ondes primaires, le second en supposant que les ondes sismiques se propagent vers le bas comme les fantômes.
La propagation vers le haut (onde montante) est définie par des angles de propagation par rapport à l'horizontale compris entre 0° et 180°, la propagation vers le bas (onde descendante) par des angles de propagation de 180° a 360°. Les procédés décrits dans US 4 353 121 et US 4 992 992 sont des procédés de traitement sismiques en une dimension (1D) et en deux dimensions (2D). Or de tels procédés ne sont pas généralisables en trois dimensions. En effet l'échantillonnage des capteurs dans la 3ème dimension est donné par l'écartement latéral entre les flûtes, de l'ordre de 150m, bien plus grand que l'échantillonnage des capteurs le long des flûtes qui est de l'ordre de 12.5 m.
Un but général de la présente invention est de fournir un procédé de traitement sismique 3D permettant d'imager la géologie du sous-sol à partir de données sismiques marines enregistrées au moyen de récepteurs sismiques dont certains au moins sont situés à des profondeurs d'eau différentes sans mettre en jeu une étape de datuming et sans restrictions quant à la direction de propagation des ondes. PRESENTATION DE L'INVENTION A cet effet, l'invention propose un procédé de traitement de données sismiques marines acquises au moyen d'au moins une flûte sismique remorquée par un bateau et comportant une pluralité de récepteurs sismiques (Rt, ..., Rn) sensibles à la pression situés à des profondeurs d'eau respectives (z1 , ...,z'), caractérisé en ce qu'il comprend une étape de migration des dites données agencée de telle sorte que les événements primaires bénéficient d'une sommation cohérente, une étape de migration miroir des dites données ayant subi une inversion de signe, agencée de telle façon que les événements fantômes bénéficient d'une sommation cohérente, et une étape de déconvolution conjointe dans laquelle on obtient à partir des images d1 et d2 fournies respectivement par la migration et la migration miroir une image r qui satisfait simultanément les relations = gn-nr, * r d2 - gmax* r gmin étant un opérateur causal et gmax un opérateur anticausal.35 PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des dessins annexés, sur lesquels : - la figure 1 est une représentation schématique d'une méthode d'acquisition sismique marine, - la figure 2 est une représentation schématique en vue de dessus de la méthode illustrée à la figure 1, la figure 3 est une représentation schématique illustrant une étape de migration miroir du procédé de traitement. - la figure 4 représente l'image fournie par la migration de données synthétiques représentative d'une méthode d'acquisition sismique marine avec récepteurs à des profondeurs différentes, les ondes fantômes résiduelles étant les halos blanchâtres en dessous de chaque réflecteur noir, - la figure 5 représente l'image fournie par la migration miroir des mêmes données, les ondes fantômes résiduelles étant dans ce cas les halos blanchâtres au dessus de chaque réflecteur, la figure 6 représente le résultat de la déconvolution conjointe des images des figures 4 et 5, et représente l'image « déghostée », c'est-à-dire exempte d'ondes fantômes résiduelles, - la figure 7 est une représentation schématique du procédé de déghosting. 25 DESCRIPTION DETAILLEE D'UN MODE DE REALISATION On va maintenant décrire le procédé de traitement de données sismiques marines en référence aux figures.
Acquisition des données sismiques marines A l'heure actuelle, les données sismiques marines sont enregistrées grâce à un dispositif d'acquisition dans lequel une même position de source donne lieu à des enregistrements de signaux sismiques au moyen d'un jeu de flûtes de sorte à couvrir une zone géographique. Dans le mode d'acquisition illustré par la figure 1, chaque flûte présente un angle avec la surface de l'eau, comme proposé dans les documents US 4 353 121 et US 4 992 992.
Cet angle est identique pour toutes les flûtes du jeu de flûtes de sorte qu'elles s'étendent globalement dans un même plan. La figure 1 représente une masse d'eau s'étendant au dessus d'un fond marin, avec à la surface de l'eau 2 un navire 3 de prospection sismique. Le navire 3 remorque une ou plusieurs source(s) sismiques 6 destinées à émettre une onde acoustique dans l'eau. La source 6 peut être un ensemble de canons à air comprimé, un vibrateur marin ou d'autres types de source connus de l'homme de l'art. Le navire 3 remorque également un ensemble de flûtes 1 présentant un angle avec la surface de l'eau 2.
Chaque flûte 1 comprend une pluralité de récepteurs sismiques 4, 5 destinés à enregistrer les signaux acoustiques émis par la source 6 après leurs réflexions successives sur les interfaces entre couches géologiques - appelées réflecteurs. Ces récepteurs sismiques 4, 5 sont de façon appropriée des capteurs sensibles à la pression, constitués dans cet exemple de réalisation par des hydrophones. Le dispositif d'acquisition fonctionne comme suit. La vibration sismique émise par la source sismique 6 descend le long de plusieurs trajectoires 11 et est réfléchie aux interfaces entre des matériaux d'impédances acoustiques différentes telles que l'interface 8 - le fond de la mer est référencé 8a.
Un champ d'ondes réfléchies 12 remonte vers la surface de l'eau 2 et est enregistré par les récepteurs sismiques 4, 5. Les hydrophones 4, 5 aux première et deuxième profondeurs z1, z2 enregistrent les ondes réfléchies du champ d'ondes montant 17.
Toutefois, et comme indiqué dans la partie « Etat de la technique », les enregistrements sont affectés par des réflexions parasites : un champ d'ondes descendant 13 dû à la réflexion des ondes à la surface de l'eau 2 se superpose avec le champ d'ondes montant 12. Les enregistrements comportent donc des pics correspondant à des réflexions de surface, ou réflexions parasites fantôme. Le procédé de traitement sismique décrit ci-après permet d'utiliser ces réflexions parasites fantôme pour imager le sous sol. Traitement des données sismiques marines Le procédé décrit ci-après utilise la migration 3D par point de tir, qui permet d'obtenir une image précise du sous-sol en tenant compte avec exactitude de la propagation des ondes dans des milieux complexes. Une telle migration consiste à synthétiser l'onde incidente à partir des informations sur la source sismique et l'onde réfléchie a partir des données enregistrées. Dans le cas d'une migration classique de type « one-way », le principe est le suivant. L'onde incidente D (i.e. l'onde émise par la source) est supposée descendante. Cette onde incidente D(x,y,z,t) est synthétisée de façon récursive en profondeur z, l'onde descendante étant initialisée à la profondeur de la source sismique zs. L'onde incidente D à toutes les profondeurs nAz est ensuite calculée de manière récursive en calculant l'onde incidente à la profondeur z+Az à partir de l'onde incidente à la profondeur z.
De façon analogue, l'onde réfléchie U(x,y,z,t) est supposée montante et est initialisée en z=zr avec les ondes enregistrées par les récepteurs sismiques (dans le cas où tous les récepteurs ont la même profondeur). L'onde réfléchie U dans tout le volume est ensuite calculée de manière récursive en calculant l'onde montante U à la profondeur z+Az à partir de l'onde monte à la profondeur z.
L'image du sous-sol est calculée par la cross-corrélation temporelle des deux volumes D(x,y,z,t) et U(x,y,z,t). L'altimétrie, c'est-à-dire le fait que la source et les récepteurs peuvent avoir des profondeurs non nulles (et différentes), peut être prise en compte en introduisant les sources et les récepteurs au cours des calculs récursifs en z : par exemple un récepteur à une profondeur zr comprise entre nAz et (n+1) Oz est introduit lors du calcul récursif de U((n+1) Oz) à partir de U(nAz). Par ailleurs, l'étape de migration est une migration miroir adaptée, ainsi nommée par analogie avec la migration miroir (ou « mirror migration » selon la terminologie anglo-saxonne) et le filtre adapté utilisé en traitement du signal (consistant à convoluer une mesure s(t) perturbé par une convolution par un filtre h(t) par h(-t) de façon à optimiser le rapport signal-à-bruit). Dans le cas d'une migration miroir, on utilise la surface de la mer comme un miroir : au lieu de « regarder» vers le fond marin, on « regarde » vers la surface de l'eau pour voir les réflecteurs situés en dessous des récepteurs sismiques. En pratique, on considère les données sismiques non pas comme si elles avaient été enregistrées au niveau des récepteurs sismiques de la flûte, mais à une altitude au dessus de la surface de l'eau égale à la profondeur de chaque récepteur comme illustré à la figure 3. Une technique d'imagerie miroir par migration miroir est décrite par exemple dans la publication : "Facilitating technologies for permanently instrumented oil fields", Dan Ebrom, Xiuyuan Li, and Dwight Sukup, The Leading Edge, Vol. 19, No 3, pp. 282-285, March 2000. Dans cette publication, cette technique est utilisée dans le cadre d'une acquisition avec des récepteurs sismiques situés au fond de la mer 8a. Dans cette acquisition, on utilise le principe de réciprocité et on considère donc fictivement des sources au fond de la mer (aux positions récepteurs) et des récepteurs en surface (aux positions sources). L'imagerie miroir consiste à utiliser le fantôme source fictive pour obtenir l'image, ce qui peut s'effectuer en plaçant les sources fictives à leur position miroir par rapport à la surface de l'eau, les positions sources (xs,ys,zs) étant changées en (xs,ys,-zs).
L'imagerie miroir permet de mieux éclairer les réflecteurs peu profonds. Dans le cas de la migration miroir adaptée, appliquée à des données acquises par des flûtes partiellement pentées, (xr,yr,zr) étant les positions des récepteurs sur les flûtes, on initialise l'onde réfléchie (supposée montante) U d'une migration avec altimétrie à une altitude -zmax, Zmax étant la profondeur maximale des récepteurs sismiques (le maximum des zr) et l'altitude 0 correspondant à la surface de l'eau. Lors de la descente récursive en z de l'onde U entre les valeurs -zmax et 0, on introduit aux positions miroirs des récepteurs par rapport au fond de l'eau, c'est-à-dire en (xr,yr,-zr), l'enregistrement du récepteur en question avec un changement de signe. Continuant la descente pour z=0 à zmax, on introduit à leur positions réelles (xr,yr,zr), l'enregistrement du récepteur en question. La suite du calcul récursif de U, la génération de l'onde incidente (supposée descendant) D ainsi que l'étape de cross-corrélation entre onde incidente et réfléchie pour obtenir l'image s'effectue de manière analogue à une migration one-way conventionnelle. On obtient de la sorte directement l'image du sous-sol à partir des acquisitions 3D par flûtes pentés, en prenant en compte les positions exactes des récepteurs et la propagation exacte en 3 dimensions des ondes.
L'étape d'introduction des enregistrements aux positions miroirs des récepteurs, dont le surcoût est négligeable, permet le renforcement du rapport signal-à-bruit par une image basé sur le fantôme récepteur, sans doubler le temps de calcul de la migration, ce qui serait le cas si on calculait deux images pour les sommer ensuite, comme proposé dans "Facilitating technologies for permanently instrumented oil fields". Toutefois une telle solution peut être mise en oeuvre dans une variante d'exécution du procédé. L'invention décrite ci-dessus permet d'obtenir directement une image du sous-sol à partir des données provenant d'une acquisition 3D, effectuée avec plusieurs flûtes pentées.
Contrairement au procédé décrits dans US 4 353 121 et US 4 992 992, le procédé de traitement décrit ci-dessus ne comprend pas d'étape de datuming consistant à reconstituer les données sismiques qui auraient été enregistrées par une flûte horizontale à partir des données sismiques enregistrées par la flûte pentée - préalablement à leur migration.
Le procédé de traitement décrit ci-dessus tient compte des angles de propagation à la fois en x et en y. Ce procédé permet également d'améliorer le rapport signal-à-bruit en utilisant les données fantômes pour renforcer les données réflexions primaires.
Lorsque la diversité des profondeurs de capteurs ne permet pas d'enlever complètement les ondes fantômes, la perturbation résultante sur les données finales est la convolution par un filtre symétrique (à phase nulle) et déconvoluable (sans notch). Cette étape de déconvolution est simplifiée par le fait que c'est une déconvolution à phase nulle.
Le descriptif de la migration miroir adapté ci-dessus concerne le cas d'une migration 3D pour point de tir « one-way ». Il existe d'autres types de migrations que l'on peut adapter en migration miroir adaptée en introduisant, dans le calcul de l'onde réfléchie, outre les enregistrements des récepteurs à leurs positions exactes, l'opposé des enregistrements à leurs positions miroirs.
Il existe en outre une migration 3D par point de tir appelée « Reverse Time Migration », qui ne fait pas l'hypothèse que l'onde incidente est une onde descendante et l'onde réfléchie une onde montante. La migration miroir adaptée peut s'effectuer dans ce cas en introduisant les récepteurs à leur position effective (xr,yr,zr) mais en utilisant à la surface de l'eau des conditions aux limites dites de surface libre, à la place des conditions aux limites de surface absorbantes habituellement utilisées. Une fois la migration effectuée, on obtient une image qui est entachée d'ondes fantômes résiduelles. Dans le cas où l'image est obtenue par migration miroir adaptée, le fantôme résiduel est symétrique. On peut considérer que l'image obtenue d(x,y,z) est égale à la réflectivité r(x,y,z) convoluée en z par la fonction de transfert g(z) du fantôme résiduel : d(x,y,z)=g(z) *r(x,y,z) La fonction de transfert du fantôme résiduel g(z) ne dépend que faiblement de la position (x,y,z) de sorte que cette relation est valide dans un volume donné.
L'opération consistant à estimer la réflectivité r(x,y,z) à partir des données d(x,y,z) s'appelle la déconvolution. Cette opération nécessite l'estimation de la fonction de transfert g(z). Deux procédés connus de l'homme de l'art permettent l'estimation de la fonction de transfert g(z) et le calcul de la réflectivité r(x,y,z) à partir des données d(x,y,z) : - la déconvolution à phase nulle et à réflectivité blanche : on suppose que g(z) est symétrique en z et que la réflectivité a une autocorrélation en z blanche (égale à une impulsion en z=0). La transformée de Fourier G(kz) de g(z) est réelle, et celle de r(x,y,z), R(x,y,kz) est complexe de module 1, on en déduit facilement que G(kz) est le module de D(x,y,kz), transformée de Fourier en z de d(x,y,z). - la déconvolution à phase minimale et à réflectivité blanche : on suppose que g(z) est à phase minimale et que la réflectivité a une autocorrélation en z blanche. Cette déconvolution à phase minimale est décrite par exemple dans les publications : "Wiener spiking deconvolution and minimum-phase wavelets: A tutorial", Jim Leinbach, The Leading Edge, Vol. 14, No 3, pp. 189-192, March 1995 "Seismic Data Processing", Ozdogan Yilmaz, ouvrage publié par Society of Exploration Geophysicists 1987, chapitre 2, pages 83 et suivantes.
La déconvolution à phase nulle permet de déconvoluer la migration miroir adaptée (puisque dans ce cas g(z) est symétrique) et la déconvolution à phase minimale permet de déconvoluer la migration standard, puisque dans ce cas la fonction de transfert du fantôme résiduel est causale et peut être considérée à phase minimale. Cependant ces deux procédés d'obtention de la réflectivité (migration suivie de déconvolution à phase minimale ou migration miroir adaptée suivie de déconvolution à phase nulle) nécessitent l'hypothèse de réflectivité blanche. Cette hypothèse était faite couramment dans le traitement sismique mais l'est de moins en moins car le traitement dit en amplitude préservé devient de plus en plus la norme : dans ce type de traitement, on s'intéresse non seulement à la position des réflecteurs mais aussi à leur amplitude, et dans ce contexte, on ne peut utiliser l'hypothèse de réflectivité blanche. Le mode de traitement décrit ci-après en référence aux figures 4 à 7 ne requiert pas l'hypothèse de réflectivité blanche et préserve l'amplitude. Dans ce mode de traitement, on part de deux images : di (x,y,z) est l'image obtenue par migration (où l'on introduit les enregistrements de chaque récepteur à leur position réelle (xr,yr,zr)) et d2(x,y,z) est l'image obtenue par migration miroir (où l'on introduit les enregistrements de chaque récepteur avec changement de signe à leur position miroir (xr,yr,-zr)).
La migration a pour effet d'aligner les événements primaires qui bénéficient alors d'une sommation cohérente, alors que les événements fantômes correspondent à des positions z plus grandes que les événements primaires correspondants. L'image d1(x,y,z) est donc entachée d'ondes fantômes résiduelles représentées par une fonction de transfert gi','(z), causale, voire à phase minimale, qui contamine de manière convolutionnelle la réflectivité r(x,y,z): di (x,y,z)=grrur,(z) *r(x,y,z). (1) une fonction à phase minimale étant une fonction causale dont l'inverse pour la convolution est également causale. La migration miroir a pour effet d'aligner les événements fantômes en changeant leurs signe pour faire correspondre leur polarité à celle des événements primaires. Ils bénéficient d'une sommation cohérente, alors que les événements primaires correspondent à des positions z plus petites que les événements fantômes correspondants. L'image d2(x,y,z), est donc entachée d'ondes fantômes résiduelles représentées par une fonction de transfert gmax(z), anticausale, voire à phase maximale, qui contamine de manière convolutionnelle la réflectivité r(x,y,z): d2(x,y,z)=gmax(z) *r(x,y,z) (2) une fonction à phase maximale étant une fonction anticausale dont l'inverse pour la convolution est également anticausale.
Après avoir calculé d1(x,y,z) et d2(x,y,z) par migration et migration miroir, on obtient la réflectivité r(x,y,z) par un procédé de « déconvolution conjointe ». On appelle déconvolution conjointe tout procédé de calcul permettant d'obtenir r(x,y,z), un opérateur causal gi','(z) et un opérateur anticausal gmax(z) à partir de d1(x,y,z) et d2(x,y,z) de manière que les relations (1) et (2) soient vérifiées, de manière exacte ou approchée à l'intérieur d'un certain volume de calcul V. Selon une définition plus spécifique, l'opérateur causal gi','(z) est un opérateur à phase minimale, et l'opérateur anticausal gmax(z) est un opérateur à phase maximale.
Un mode de réalisation préférentiel de la déconvolution conjointe consiste à définir un volume V = [Xmin,-rna,' x Lj rV V r7 7 ainsi qu'une longueur Dz dépendante de l'écart maximum entre un évènement et son fantôme, puis à calculer gmin(z), gmax(z) et r(x,y,z), r étant défini sur le volume V, gmm(z) sur l'intervalle [O,Dz] avec la normalisation gmm(z=0)=1, gmax(z) sur l'intervalle [-Dz,O] avec la normalisation gmax(z=0)=1, par minimisation de la fonction coût C: C = E(x,y,z) e v [di (x,Y,z)-gmin(z)*qx,y,z) 12 + [d2(x,y,z)-gmax(z)*r(x,y,z) ]2 La réflectivité r(x,y,z) étant calculée sur tout le volume d'intérêt en juxtaposant avec une zone de recouvrement les r(x,y,z) calculés sur des volumes V.
On peut aussi imposer le caractère à phase minimale de gmm(z) et le caractère à phase maximale de gmax(z). Un autre mode de réalisation consiste à remplacer les fonctions en z gmm(z), gmax(z) par des fonctions à trois dimensions gmm(x,y,z), gmax(x,y,z) qui sont causale en z et anticausale en z respectivement, la convolution en z notée * étant remplacée par une convolution en 3 dimensions (x,y,z). Ce mode de réalisation permet de prendre en compte la dépendance du fantôme résiduel avec les angles de propagation. Dans un autre mode de réalisation, qui permet lui aussi de prendre en compte la dépendance du fantôme résiduel avec les angles de propagation, on applique sur d1(x,y,z), d2(x,y,z) une transformée dite transformée (r,px,py) qui transforme les données d1(x,y,z) en Di(Px,py,T) et les données d2(x,y,z) en D2(Px, PY,T) (pour une définition de la transformée (r,px,py), également appelée « slant stack », on peut se reporter à l'ouvrage déjà cité plus haut « Seismic Data Processing », Ozdogan Yilmaz, Society of Exploration Geophysicists 1987, chapitre 7, pages 429 et suivantes. On calcule alors, pour toute les valeurs (px,py) des fantômes résiduels Gmin(Px,P0(T) et Gmax(Px,130(t), supposés respectivement causaux et anticausaux en T, ou à phase minimale et maximale en T, et une réflectivité R(px,py,T) tels que : C(px,py) = E [Di(Px,Py,T) - Gmin(Px, PeCY R(Px, Py:012 [D2(px, !Dy, T) Gmax(Px, py) (T)*R(px, T)]2 soit minimale pour tout (p,,py), l'opération * étant dans ce cas la convolution en T.. La migration déghostée, c'est-à-dire la réflectivité cherchée r(x,y,z) se déduit par transformée (,px,py) inverse de R(px,Py,T).
Dans la description ci-dessus, les migrations utilisées sont des migrations profondeur, où les images obtenues sont d1(x,y,z) et d2(x,y,z). On peut également utiliser le procédé de déconvolution conjointe avec des images provenant de migrations temps. Dans le cas d'une migration temps, l'image est d(x,y,T), le paramètre de type temporel T remplaçant le paramètre profondeur z, et la focalisation se fait à l'aide d'un modèle de vitesse vi(x,y,T) visant à ce qu'au moment de la sommation, les événement primaires bénéficient d'une sommation cohérente. L'équivalent de la migration miroir en traitement temps est une migration où après changement de polarité des données d'entrées, le modèle de vitesse v2(x,y,T) utilisé pour la focalisation est optimisé pour qu'au moment de la sommation, les événement fantômes bénéficient d'une sommation cohérente. On peut alors appliquer la déconvolution conjointe correspondant au modèle : di(x,y,T)=g,,,'(T)*r(x,y,T) et d2(x,y,T)=g,,,x(i)*r(x,y,T). La déconvolution conjointe peut être généralisée pour déconvoluer plus de 2 jeux de données, en particulier dans le cas où les récepteurs sont de type différents : récepteurs sensibles à la pression tels que hydrophone pour certains et récepteurs de type géophone pour d'autres. On peut obtenir séparément la migration d1(x,y,z) et la migration miroir d2(x,y,z) des récepteurs de type hydrophone et la migration d3(x,y,z) et la migration miroir d4(x,y,z) des données géophones. La migration miroir de données de type géophone s'obtient en introduisant les enregistrements de chaque récepteur à leur positions miroirs (xr,yr,-zr) mais sans changement de signe pour un récepteur géophone vertical et avec changement de signe pour un récepteur de type géophone horizontal. La déconvolution conjointe à 4 entrées s'obtient en modélisant les migrations et migrations miroirs par: di (x,y,z)=ghmir,(z) *r(x,y,z). d2(x,y,z)=ghniax(z) *r(x,y,z). d3(x, y, z)=ggn,;'(z) *c(z)*r(x, y, z). d4(x,y,z)=g9niax(z) *c(z)*r(x,y,z). et en estimant à l'aide d'une fonction coût de type moindre carrés les opérateurs causaux ou à phase minimale gh,,,'(z) et gg,,,'(z), les opérateurs anticausaux ou à phase maximale ghmax(z) et ggniax(z), l'opérateur de calibration c(z) ainsi que la réflectivité r(x,y,z) à partir des migrations et migrations miroirs d1(x,y,z), d2(x,y,z), d3(x,y,z) et d4(x,y,z). La fonction coût peut être pondérée pour prendre en compte la différence de spectre du bruit des capteurs hydrophone et des capteurs géophone. Par exemple, dans le cas de migrations temps, la fonction coût à minimiser s'écrit alors dans le domaine f, transformée de Fourier de : C = E [di(x,y,f)-gh,'''(f)r(x,y,f) ]2 /Bh(f) [d2(x,y,f)-gh,,,x(f)r(x,y,f) ]2 /Bh(f) [d3(x,y,f)-g9,,,'(f)c(f)r(x,y,f) ]2/B9(f) + [d4(x,y,f)-gg,,,x(f)c(f)(x,y,f) ]2/B9(f) où Bh(f) et B9(f) sont les estimations du spectre de puissance des bruits hydrophone et géophones respectivement.
Les procédés décrits ci-dessus ne sont pas limités au traitement de données acquises au moyen des flûtes linéaires de pente constante représentées à la figure 1. Ils sont tout aussi applicables à des données acquises au moyen de flûtes comportant chacune plusieurs sections de pentes différentes, ou de flûtes ayant une ou plusieurs sections pentées et une ou plusieurs sections horizontales, ou de flûtes horizontales situées à des profondeurs différentes.
Claims (8)
- REVENDICATIONS1. Procédé pour supprimer des fantômes de données sismiques relatives à un sous-sol d'une masse d'eau, le procédé comprenant des étapes consistant à : - appliquer une procédure de migration à des données d'ondes enregistrées par des détecteurs qui sont remorqués par un navire, les ondes se propageant du sous-sol vers les détecteurs, pour obtenir une première image (d1) du sous-sol ; - appliquer une procédure de migration miroir aux données d'ondes enregistrées par les détecteurs pour obtenir une deuxième image (d2) du sous-sol ; - déconvoluer conjointement la première image (d1) et la deuxième image (d2) pour supprimer les fantômes d'une réflectivité (r) du sous- sol ; et - obtenir une image déghostée du sous-sol sur la base de la réflectivité (r) de laquelle les fantômes ont été supprimés.
- 2. Procédé selon la revendication 1, dans lequel la suppression des fantômes est effectuée après une étape de migration et non avant.
- 3. Procédé selon la revendication 1, dans lequel aucune étape de datuming n'est appliquée aux données.
- 4. Procédé selon la revendication 1, dans lequel un angle de propagation des ondes se propageant du sous-sol vers les détecteurs ou d'une surface de l'eau vers les détecteurs n'est pas limité.
- 5. Procédé selon la revendication 1, dans lequel la déconvolution conjointe comprend : - obtenir la réflectivité r(x, y, z), une fonction de transfert à phase minimale gi','(z) et une fonction de transfert à phase maximale gmax(z) sur la base des équations : di (x, y, z) = gni,n(z) *r(x, y, z), et d2(x, y, z) = gmax(z) *r(x, y, z),(x, y, z) étant une position d'un point où z est une profondeur du point par rapport à la surface de l'eau.
- 6. Procédé selon la revendication 1, dans lequel la déconvolution conjointe comprend : - calculer une fonction coût C pour obtenir la réflectivité, telle que : C = Z(x,y,z) E V { [di (X, y, z) - gmir,(z)*r(x, y, z)]2 + [d2(x, y, z) - gmax(z)*r(x, y, z)]2}, où di(x, y, z) est la première image, d2(x, y, z) est la deuxième image, g,','(z) est une fonction de transfert à phase minimale, gmax(z) est une fonction de transfert à phase maximale, (x, y, z) une position d'un point, z étant une profondeur du point par rapport à la surface de l'eau, et V est un volume prédéterminé.
- 7 Procédé selon la revendication 1, comprenant en outre : - appliquer une transformée (r, px, py) à la première image d1(x, y, z) et à la deuxième image d2(x, y, z), pour transformer la première image d1(x, Y, z) en Di(Px, py, t) et la deuxième image d2(x, y, z) en D2(Px, py, t). 20
- 8. Procédé selon la revendication 1, dans lequel les données comprennent des enregistrements provenant d'hydrophones et de géophones, dans lequel un résultat de la procédure de migration est la première image d1(x, y, z) et un résultat de la procédure de migration miroir est la 25 deuxième image d2(x, y, z) pour des récepteurs du type hydrophone et un résultat de la procédure de migration est une migration d3(x, y, z) et un résultat de la procédure de migration miroir est une migration d4(x, y, z) pour des géophones, et le procédé comprend en outre une étape consistant à : - obtenir l'image déghostée en utilisant une déconvolution conjointe de la 30 première image d1(x, y, z), de la deuxième image d2(x, y, z), de la migration d3(x, y, z) et de la migration d4(x, y, z) et sur la base des équations suivantes : d1(x, y, z) = *r(x, y, z) ; 15d2(x, y, z) = ghmax(z)*r(x, y, z) ; d3(x, y, z) = gg,','(z)*c(z)*r(x, y, z) ; et d4(x, y, z) = ggmax(z)*c(z)*r(x, y, z), où ghmm et ggm,n sont des opérateurs à phase minimale, ghmax(z) et ggmax(z) sont des opérateurs à phase maximale, (x, y, z) une position d'un point, z étant une profondeur du point par rapport à la surface de l'eau, et c(z) est un opérateur de calibration.10
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1260128A FR2979711B1 (fr) | 2010-06-10 | 2012-10-24 | Procede de traitement de donnees sismiques marines |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1054599A FR2961316A1 (fr) | 2010-06-10 | 2010-06-10 | Procede de traitement de donnees sismiques marines |
FR1260128A FR2979711B1 (fr) | 2010-06-10 | 2012-10-24 | Procede de traitement de donnees sismiques marines |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2979711A1 true FR2979711A1 (fr) | 2013-03-08 |
FR2979711B1 FR2979711B1 (fr) | 2015-07-03 |
Family
ID=43558264
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1054599A Withdrawn FR2961316A1 (fr) | 2010-06-10 | 2010-06-10 | Procede de traitement de donnees sismiques marines |
FR1260128A Expired - Fee Related FR2979711B1 (fr) | 2010-06-10 | 2012-10-24 | Procede de traitement de donnees sismiques marines |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1054599A Withdrawn FR2961316A1 (fr) | 2010-06-10 | 2010-06-10 | Procede de traitement de donnees sismiques marines |
Country Status (10)
Country | Link |
---|---|
US (4) | US8451682B2 (fr) |
EP (1) | EP2395374B1 (fr) |
CN (1) | CN102298156B (fr) |
AU (1) | AU2011202739B2 (fr) |
CA (1) | CA2780081A1 (fr) |
FR (2) | FR2961316A1 (fr) |
MX (1) | MX2011006153A (fr) |
MY (1) | MY191761A (fr) |
SG (1) | SG186240A1 (fr) |
WO (1) | WO2011154463A1 (fr) |
Families Citing this family (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7400552B2 (en) | 2006-01-19 | 2008-07-15 | Westerngeco L.L.C. | Methods and systems for efficiently acquiring towed streamer seismic surveys |
US8681580B2 (en) | 2008-05-15 | 2014-03-25 | Westerngeco L.L.C. | Multi-vessel coil shooting acquisition |
US9857491B2 (en) | 2008-05-15 | 2018-01-02 | Westerngeco L.L.C. | Multi-vessel coil shooting acquisition |
US9052411B2 (en) | 2008-06-13 | 2015-06-09 | Westerngeco L.L.C. | Method to determine the deviation of seismic equipment from a planned curved path |
US9594181B2 (en) | 2008-06-13 | 2017-03-14 | Westerngeco L.L.C. | Filtering and presentation of heading observations for coil shooting |
FR2955396B1 (fr) | 2010-01-15 | 2013-03-01 | Cggveritas Services Sa | Dispositif de traitement de donnees sismiques marines |
FR2955397B1 (fr) | 2010-01-15 | 2012-03-02 | Cggveritas Services Sa | Procede et dispositif d'acquisition de donnees sismiques marines |
FR2961316A1 (fr) * | 2010-06-10 | 2011-12-16 | Cggveritas Services Sa | Procede de traitement de donnees sismiques marines |
AU2011232767B2 (en) | 2010-10-14 | 2014-05-08 | Cggveritas Services Sa | Method and device to acquire seismic data |
US8949030B2 (en) | 2011-07-29 | 2015-02-03 | Westerngeco L.L.C. | Attenuating sea-surface ghost wave effects in seismic data |
US20130028045A1 (en) * | 2011-07-29 | 2013-01-31 | Ralf Ferber | Seismic survey designs for attenuating sea-surface ghost wave effects in seismic data |
US10459099B2 (en) * | 2011-09-22 | 2019-10-29 | Cgg Services Sas | Device and method to determine shape of streamer |
US9103942B2 (en) * | 2011-10-28 | 2015-08-11 | Westerngeco L.L.C. | Methods and systems for survey designs |
US9103943B2 (en) | 2011-11-28 | 2015-08-11 | Fugro-Geoteam As | Acquisition and processing of multi-source broadband marine seismic data |
US9103941B2 (en) * | 2011-12-21 | 2015-08-11 | Cggveritas Services Sa | Device and method for deghosting variable depth streamer data |
US9176249B2 (en) | 2011-12-21 | 2015-11-03 | Cggveritas Services Sa | Device and method for processing variable depth streamer data |
US8634271B2 (en) * | 2012-01-11 | 2014-01-21 | Cggveritas Services Sa | Variable depth streamer SRME |
US9234977B2 (en) * | 2012-01-12 | 2016-01-12 | Westergeco L.L.C. | Processing collected survey data |
US10977396B2 (en) | 2012-01-13 | 2021-04-13 | Schlumberger Technology Corporation | Determining an elastic model for a geologic region |
US9256925B2 (en) | 2012-01-31 | 2016-02-09 | Shell Oil Company | Angle-dependent deghosting of seismic data |
US20130265849A1 (en) * | 2012-04-04 | 2013-10-10 | Westerngeco L.L.C. | Methods and devices for enhanced survey data collection |
US9435905B2 (en) * | 2012-04-19 | 2016-09-06 | Cgg Services Sa | Premigration deghosting of seismic data with a bootstrap technique |
US9541661B2 (en) | 2012-04-19 | 2017-01-10 | Cgg Services Sa | Device and method for deghosting variable depth streamer data |
US9322943B2 (en) | 2012-05-18 | 2016-04-26 | Cggveritas Services Sa | Method and apparatus for pre-stack deghosting of seismic data |
WO2013190513A2 (fr) | 2012-06-22 | 2013-12-27 | Schlumberger Technology B.V. | Traitement de données sismiques par empilement non linéaire |
US9665604B2 (en) * | 2012-07-31 | 2017-05-30 | Schlumberger Technology Corporation | Modeling and manipulation of seismic reference datum (SRD) in a collaborative petro-technical application environment |
US9453930B2 (en) * | 2012-11-01 | 2016-09-27 | Pgs Geophysical As | Systems and methods for high-resolution imaging using separated wavefields |
CN103018779B (zh) * | 2012-11-30 | 2015-10-21 | 中国石油大学(北京) | 一种海上地震勘探气枪震源子波模拟方法及系统 |
FR3000224B1 (fr) * | 2012-12-26 | 2015-01-09 | Cggveritas Services Sa | Method and device to deghost seismic data |
MX367875B (es) * | 2013-01-11 | 2019-09-10 | Cgg Services Sa | Retiro de fantasmas pre-migracion para datos de cable marino con el uso de una medida de muestreo en el dominio tau-p. |
WO2014114683A2 (fr) * | 2013-01-22 | 2014-07-31 | Cgg Services Sa | Systèmes et procédés pour débruitage directionnel multivolume |
AU2014209936B2 (en) | 2013-01-23 | 2018-02-08 | Cgg Services Sa | Low frequency emission and recording for seismic data acquisition |
US20140249757A1 (en) * | 2013-03-04 | 2014-09-04 | Bruno Gratacos | Apparatus and method for determination of far-field signature from variable-depth seismic data |
US20140249758A1 (en) * | 2013-03-04 | 2014-09-04 | Cgg Services Sa | Device and method for residual moveout picking of ghosted seismic data |
US9594173B2 (en) | 2013-03-08 | 2017-03-14 | Cgg Services Sas | Method and system for augmenting frequency range of conventional marine seismic source with low-frequency |
US10088588B2 (en) * | 2013-04-03 | 2018-10-02 | Cgg Services Sas | Device and method for stable least-squares reverse time migration |
EP2992362B1 (fr) | 2013-04-29 | 2022-01-19 | CGG Services SAS | Dispositif et procédé pour reconstruction de champ d'ondes |
US9791581B2 (en) * | 2013-06-07 | 2017-10-17 | Cgg Services Sas | Method and system for simultaneous acquisition of pressure and pressure derivative data with ghost diversity |
US9442208B2 (en) | 2013-06-07 | 2016-09-13 | Cgg Services Sa | Device and method for deghosting variable depth streamer data including particle motion data |
US10436923B2 (en) | 2013-06-11 | 2019-10-08 | Cgg Services Sas | Method and apparatus for receiver-side deghosting of seismic data |
US9678235B2 (en) | 2013-07-01 | 2017-06-13 | Pgs Geophysical As | Variable depth multicomponent sensor streamer |
US9671513B2 (en) * | 2014-01-10 | 2017-06-06 | Cgg Services Sas | Seismic image acquisition using ghost compensation in Kirchhoff migration |
WO2015104637A2 (fr) | 2014-01-13 | 2015-07-16 | Cgg Services Sa | Dispositif et procédé permettant de déparasiter des données sismiques au moyen d'une inversion tau-p éparse |
US9964655B2 (en) | 2014-04-01 | 2018-05-08 | Schlumberger Technology Corporation | Deghosting after imaging |
WO2015155595A1 (fr) | 2014-04-09 | 2015-10-15 | Cgg Services Sa | Procédé et système permettant de générer des signaux sismiques de basse fréquence avec une source à écoulement modulé |
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 |
CN104898166B (zh) * | 2015-05-19 | 2017-09-15 | 中国石油天然气集团公司 | 一种处理地震数据的方法及装置 |
EP3118655B1 (fr) | 2015-07-17 | 2022-09-07 | Sercel | Procédé et système d'acquisition simultanée de données sismiques de multiples lignes de sources |
US10634804B2 (en) | 2015-12-21 | 2020-04-28 | Chevron U.S.A. Inc. | System and method for dip-guided seismic image stacking |
US10520640B2 (en) | 2016-06-23 | 2019-12-31 | Cgg Services Sas | Method and system for adjusting seismic source pilot and drive signals using measured data |
US10732311B2 (en) | 2016-10-19 | 2020-08-04 | Pgs Geophysical As | Method for adjusting attributes of marine seismic survey data |
CN106873027B (zh) * | 2017-02-07 | 2019-01-04 | 国家海洋局第一海洋研究所 | 一种双震源宽方位角海上三维地震探测方法 |
US20180228393A1 (en) | 2017-02-15 | 2018-08-16 | Biosense Webster (Israel) Ltd. | Electrophysiologic device construction |
US11209564B2 (en) | 2018-07-18 | 2021-12-28 | Saudi Arabian Oil Company | Deghosting of seismic data through echo- deblending using coincidence filtering |
US11598894B2 (en) | 2020-04-21 | 2023-03-07 | Sercel | Method and system for seismic data acquisition with top and front sources |
Family Cites Families (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2614165A (en) | 1949-02-23 | 1952-10-14 | Stanolind Oil & Gas Co | Marine seismic surveying |
FR1050278A (fr) | 1950-05-02 | 1954-01-06 | Technical Assets | Procédé de fabrication d'éléments piézoélectriques artificiels |
FR1052576A (fr) | 1951-03-19 | 1954-01-26 | Licentia Gmbh | Interrupteur électrique de courant avec contact en parallele |
FR1054599A (fr) | 1952-04-17 | 1954-02-11 | Serrure | |
US3299399A (en) | 1965-12-02 | 1967-01-17 | Vector Cable Company | Methods and apparatus for indicating an underwater parameter in a marine seismic system |
US3491848A (en) | 1968-01-10 | 1970-01-27 | Texas Instruments Inc | Wave generator arrays for marine seismic exploration |
US3885515A (en) | 1974-06-07 | 1975-05-27 | Us Navy | Rigid line array suspension system |
US4068208A (en) | 1975-07-14 | 1978-01-10 | Texas Instruments Incorporated | Marine streamer position determination system |
US4101866A (en) | 1976-08-18 | 1978-07-18 | Mobil Oil Corporation | Marine detector spread having arrays of different lengths |
US4233677A (en) | 1978-02-21 | 1980-11-11 | Chevron Research Company | Automatic system and method for testing marine streamers of seismic exploration systems |
US4254480A (en) | 1978-09-11 | 1981-03-03 | Standard Oil Company (Indiana) | Frequency independent directionally sensitive array in seismic surveying |
US4353121A (en) * | 1980-07-24 | 1982-10-05 | Fairfield Industries, Inc. | High resolution, marine seismic stratigraphic system |
US4486865A (en) | 1980-09-02 | 1984-12-04 | Mobil Oil Corporation | Pressure and velocity detectors for seismic exploration |
US4404664A (en) | 1980-12-31 | 1983-09-13 | Mobil Oil Corporation | System for laterally positioning a towed marine cable and method of using same |
US4635236A (en) | 1981-09-29 | 1987-01-06 | Chevron Research Company | Submerged marine streamer locator |
US4532617A (en) | 1982-09-29 | 1985-07-30 | Baecker Donald Ray | System for locating a towed marine object |
US4581723A (en) | 1983-04-18 | 1986-04-08 | Western Geophysical Company Of America | Method for maintaining a constant tension on a seismic marine cable |
US4486863A (en) | 1983-08-11 | 1984-12-04 | Tensor Geophysical Service Corporation | Circular seismic acquisition system |
NO162398B (no) | 1987-07-27 | 1989-09-11 | Geco As | Marin streamer for bruk ved seismiske undersoekelser. |
US4875144A (en) | 1987-09-14 | 1989-10-17 | Wainwright Harry L | Fabric with illuminated changing display |
US4992992A (en) * | 1988-10-21 | 1991-02-12 | Western Atlas International, Inc. | Processing for seismic data from slanted cable |
US4872144A (en) | 1988-10-31 | 1989-10-03 | Exxon Production Research Company | Endfire seismic receiver and method of use |
EG19158A (en) | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
US5257241A (en) | 1991-05-08 | 1993-10-26 | Atlantic Richfield Company | Method and system for acquisition of 3-dimensional marine seismic data |
US5930731A (en) | 1997-02-03 | 1999-07-27 | Pgs Tensor, Inc. | Method and system for acquisition and processing of marine seismic data |
US5959938A (en) | 1997-08-22 | 1999-09-28 | Baker Hughes Incorporated | Tuned bubble attenuator for towed seismic source |
US6049507A (en) | 1997-09-30 | 2000-04-11 | Mobil Oil Corporation | Method and apparatus for correcting effects of ship motion in marine seismology measurements |
US6088299A (en) | 1998-12-04 | 2000-07-11 | Syntron, Inc. | Vertical hydrophone array |
US6024344A (en) | 1999-02-17 | 2000-02-15 | Western Atlas International, Inc. | Method for recording seismic data in deep water |
GB9906456D0 (en) * | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
GB9927395D0 (en) * | 1999-05-19 | 2000-01-19 | Schlumberger Holdings | Improved seismic data acquisition method |
FR2795527B1 (fr) | 1999-06-22 | 2001-09-07 | Thomson Marconi Sonar Sas | Systeme de prospection sismique sous-marine, notamment pour grands fonds |
US6317695B1 (en) | 2000-03-30 | 2001-11-13 | Nutec Sciences, Inc. | Seismic data processing method |
US6574567B2 (en) | 2001-01-23 | 2003-06-03 | Pgs Americas, Inc. | Weighted slant stack for attenuating seismic noise |
BRPI0101198B1 (pt) * | 2001-03-02 | 2015-09-01 | Unicamp | Processo concorrente para desconvolução autodidata de sinais digitais |
CN1110709C (zh) | 2001-09-13 | 2003-06-04 | 中国科学院声学研究所 | 用于测量海底微地貌的高分辨测深侧扫声纳系统和测量方法 |
GB2393513A (en) * | 2002-09-25 | 2004-03-31 | Westerngeco Seismic Holdings | Marine seismic surveying using a source not having a ghost at a non-zero frequency |
US6882938B2 (en) * | 2003-07-30 | 2005-04-19 | Pgs Americas, Inc. | Method for separating seismic signals from two or more distinct sources |
MXPA06008524A (es) | 2004-01-29 | 2007-04-02 | Westerngeco Seismic Holdings | Colocacion de cable sismico utilizando unidades de sistema de medicion de inercia acopladas. |
US7466632B1 (en) | 2004-05-04 | 2008-12-16 | Westerngeco L.L.C. | Method and apparatus for positioning a center of a seismic source |
US7372769B2 (en) | 2005-04-08 | 2008-05-13 | Western Geco L.L.C. | Method and apparatus for adaptive over/under combination of seismic data |
US7417924B2 (en) | 2005-04-26 | 2008-08-26 | Westerngeco L.L.C. | Apparatus, systems and methods for determining position of marine seismic acoustic receivers |
US7768869B2 (en) | 2005-05-05 | 2010-08-03 | Pgs Americas, Inc. | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
US7660191B2 (en) | 2005-07-12 | 2010-02-09 | Westerngeco L.L.C. | Methods and apparatus for acquisition of marine seismic data |
US20100135112A1 (en) | 2005-07-12 | 2010-06-03 | Johan Olof Anders Robertsson | Methods and Apparatus for Acquisition of Marine Seismic Data |
US7400552B2 (en) * | 2006-01-19 | 2008-07-15 | Westerngeco L.L.C. | Methods and systems for efficiently acquiring towed streamer seismic surveys |
US7599798B2 (en) | 2006-09-11 | 2009-10-06 | Westerngeco L.L.C. | Migrating composite seismic response data to produce a representation of a seismic volume |
US8437979B2 (en) | 2007-01-20 | 2013-05-07 | Kcf Technologies, Inc. | Smart tether system for underwater navigation and cable shape measurement |
US20090161488A1 (en) | 2007-12-19 | 2009-06-25 | Ralf Ferber | Method to Estimate Ray Parameter for Seismograms |
US8004930B2 (en) | 2008-03-17 | 2011-08-23 | Westerngeco, L.L.C. | Methods and systems for determining coordinates of an underwater seismic component in a reference frame |
US8456948B2 (en) | 2008-06-28 | 2013-06-04 | Westerngeco L.L.C. | System and technique to obtain streamer depth and shape and applications thereof |
US7813218B2 (en) | 2008-07-05 | 2010-10-12 | Westerngeco L.L.C. | Performing quality control with respect to positioning of survey hardware |
US20100054082A1 (en) | 2008-08-29 | 2010-03-04 | Acceleware Corp. | Reverse-time depth migration with reduced memory requirements |
FR2955397B1 (fr) | 2010-01-15 | 2012-03-02 | Cggveritas Services Sa | Procede et dispositif d'acquisition de donnees sismiques marines |
FR2955396B1 (fr) | 2010-01-15 | 2013-03-01 | Cggveritas Services Sa | Dispositif de traitement de donnees sismiques marines |
EP2352040A3 (fr) | 2010-01-28 | 2013-02-27 | PGS Geophysical AS | Procédé et système de contrôle de profondeur de jet |
FR2961316A1 (fr) * | 2010-06-10 | 2011-12-16 | Cggveritas Services Sa | Procede de traitement de donnees sismiques marines |
AU2011232767B2 (en) | 2010-10-14 | 2014-05-08 | Cggveritas Services Sa | Method and device to acquire seismic data |
US10459099B2 (en) * | 2011-09-22 | 2019-10-29 | Cgg Services Sas | Device and method to determine shape of streamer |
-
2010
- 2010-06-10 FR FR1054599A patent/FR2961316A1/fr not_active Withdrawn
-
2011
- 2011-06-08 US US13/155,778 patent/US8451682B2/en active Active
- 2011-06-08 SG SG2012090437A patent/SG186240A1/en unknown
- 2011-06-08 AU AU2011202739A patent/AU2011202739B2/en not_active Ceased
- 2011-06-08 EP EP11169161.4A patent/EP2395374B1/fr not_active Not-in-force
- 2011-06-08 CA CA2780081A patent/CA2780081A1/fr not_active Abandoned
- 2011-06-08 WO PCT/EP2011/059528 patent/WO2011154463A1/fr active Application Filing
- 2011-06-08 MY MYPI2012005310A patent/MY191761A/en unknown
- 2011-06-09 MX MX2011006153A patent/MX2011006153A/es active IP Right Grant
- 2011-06-10 CN CN201110201031.XA patent/CN102298156B/zh not_active Expired - Fee Related
-
2012
- 2012-05-04 US US13/464,126 patent/US8456951B2/en active Active
- 2012-10-24 FR FR1260128A patent/FR2979711B1/fr not_active Expired - Fee Related
-
2013
- 2013-04-26 US US13/871,326 patent/US8792296B2/en active Active
-
2014
- 2014-06-04 US US14/295,663 patent/US9784868B2/en active Active
Non-Patent Citations (5)
Title |
---|
DASH R ET AL: "Wide-area imaging from OBS multiples", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 74, no. 6, 1 November 2009 (2009-11-01), pages Q41 - Q47, XP001550514, ISSN: 0016-8033, DOI: 10.1190/1.3223623 * |
LASSE AMUNDSEN ET AL: "Multidimensional signature deconvolution and free-surface multiple elimination of marine mutlicomponent ocean-bottom seismic data", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 66, no. 5, 1 September 2001 (2001-09-01), pages 1594 - 1604, XP002534700, ISSN: 0016-8033, DOI: 10.1190/1.1486770 * |
REMCO MUIJS ET AL: "PRESTACK DEPTH MIGRATION OF PRIMARY AND SURFACE-RELATED MULTIPLE REFLECTIONS:PART II - IDENTIFICATION AND REMOVAL OF RESIDUAL MULTIPLES", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, US, vol. 72, no. 2, 1 March 2007 (2007-03-01), pages S71 - S76, XP001504969, ISSN: 0016-8033, DOI: 10.1190/1.2424544 * |
SOUBARAS R.: "Deghosting by joint deconvolution of a migration and a mirror migration", SEG, EXPANDED ABSTRACTS 29 (2010), October 2010 (2010-10-01) - November 2010 (2010-11-01), Denver, Co, USA, pages 3406 - 3410, XP002622658 * |
TRAD D.: "Interpolation and multiple attenuation with migration operators", GEOPHYSICS, SEG, vol. 68, no. 6, November 2003 (2003-11-01) - December 2003 (2003-12-01), pages 2043 - 2054, XP002622659 * |
Also Published As
Publication number | Publication date |
---|---|
EP2395374B1 (fr) | 2016-11-16 |
MY191761A (en) | 2022-07-14 |
CA2780081A1 (fr) | 2011-12-15 |
US8792296B2 (en) | 2014-07-29 |
WO2011154463A1 (fr) | 2011-12-15 |
EP2395374A1 (fr) | 2011-12-14 |
US20130242695A1 (en) | 2013-09-19 |
FR2979711B1 (fr) | 2015-07-03 |
US9784868B2 (en) | 2017-10-10 |
FR2961316A1 (fr) | 2011-12-16 |
US8451682B2 (en) | 2013-05-28 |
US20120213033A1 (en) | 2012-08-23 |
US8456951B2 (en) | 2013-06-04 |
US20140288839A1 (en) | 2014-09-25 |
AU2011202739B2 (en) | 2013-12-12 |
MX2011006153A (es) | 2011-12-13 |
SG186240A1 (en) | 2013-01-30 |
CN102298156B (zh) | 2015-08-19 |
CN102298156A (zh) | 2011-12-28 |
AU2011202739A1 (en) | 2012-01-12 |
US20110305109A1 (en) | 2011-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2979711A1 (fr) | Procede de traitement de donnees sismiques marines | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
US9435905B2 (en) | Premigration deghosting of seismic data with a bootstrap technique | |
FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
CA2786411A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
US20140078860A1 (en) | Interference noise attenuation method and apparatus | |
FR2790102A1 (fr) | Procede de traitement de jeux de donnees sismiques | |
US9594180B2 (en) | Removing ghost reflections from marine seismic data | |
FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
FR2916540A1 (fr) | Procede d'exploration sismique permettant la supression de fantomes dus aux reflexions a la surface de l'eau, et procede de traitement de donnees sismiques pour la supression de ces fantomes | |
FR2795181A1 (fr) | Procede de prospection sismique | |
FR2972057A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
Carrière et al. | Deep-water subsurface imaging using OBS interferometry | |
US10324208B2 (en) | Premigration deghosting for marine streamer data using a bootstrap approach in Tau-P domain | |
Edme et al. | Near-surface imaging using ambient-noise body waves | |
FR2966253A1 (fr) | Proceder et dispositif pour retirer des fantomes de donnees sismiques | |
FR2751427A1 (fr) | Methode d'attenuation multiple | |
AU2014201114B2 (en) | Method and apparatus for deghosting seismic data | |
Grion Ltd | From PZ Sum to Wavefield Separation, Mirror Imaging and Up Down Deconvolution: the Evolution of OBS Data Processing | |
Sandvin | J. Mispel, B. Arntsen, A. Kritski, L. Amundsen, M. Thompson, StatoilHydro, O. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 7 |
|
ST | Notification of lapse |
Effective date: 20180228 |