FR2916540A1 - 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 - Google Patents
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 Download PDFInfo
- Publication number
- FR2916540A1 FR2916540A1 FR0755279A FR0755279A FR2916540A1 FR 2916540 A1 FR2916540 A1 FR 2916540A1 FR 0755279 A FR0755279 A FR 0755279A FR 0755279 A FR0755279 A FR 0755279A FR 2916540 A1 FR2916540 A1 FR 2916540A1
- Authority
- FR
- France
- Prior art keywords
- depth
- hydrophones
- seismic data
- acquisition
- seismic
- 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 47
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims description 29
- 229910001868 water Inorganic materials 0.000 title claims description 29
- 238000012545 processing Methods 0.000 title claims description 9
- 230000000630 rising effect Effects 0.000 claims description 11
- 238000003672 processing method Methods 0.000 claims description 9
- 230000003071 parasitic effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 description 12
- 238000004590 computer program Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002087 whitening effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 244000309466 calf Species 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001629 suppression 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/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
- 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
-
- 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)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
L'invention concerne un procédé d'exploration sismique marine à au moins deux profondeurs différentes au moyen de flûtes comprenant chacune un ensemble d'hydrophones, l'acquisition de données aux différentes profondeurs étant réalisée dans une même région d'exploration, et la position des flûtes et des hydrophones lors de l'acquisition des données à une première profondeur est indépendante de la position des flûtes et des hydrophones lors de l'acquisition des données à une deuxième profondeur.
Description
PROCEDE D'EXPLORATION SISMIQUE PERMETTANT LA SUPPRESSION DE FANTOMES DUS
AUX REFLEXIONS A LA SURFACE DE L'EAU, ET PROCEDE DE TRAITEMENT DE DONNEES SISMIQUES POUR LA SUPPRESSION DE CES FANTOMES La présente invention concerne l'exploration sismique marine. Elle concerne notamment un procédé d'exploration sismique marine permettant la suppression de "fantômes" et un procédé de traitement de données sismiques pour la suppression de fantômes.
Un fantôme (ou ghost selon la terminologie anglo-saxonne) est un signal parasite dû aux réflexions à la surface de l'eau.
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) pour enregistrer l'onde acoustique réfléchie par les interfaces entre couches géologiques. La source communique une onde acoustique à l'eau. Les ondes de compression créées se propagent 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 d'entrée, constituent des données sismiques et permettent de 30 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 direction30 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, où il est réfléchi vers le bas. Ce champ d'ondes réfléchi ou fantôme (ghost en anglais) se propage également dans l'eau et passe par les récepteurs où il est à nouveau enregistré. 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. Il a déjà été proposé, pour supprimer les fantômes, d'acquérir des données sismiques à deux profondeurs différentes. Cette méthode est appelée over-under selon la terminologie anglo-saxonne (pour over cable - under cable ). Le document US 2 757 356 décrit un procédé d'enregistrement de données sismiques par réflexion dans lequel deux récepteurs sismographiques sont disposés à deux profondeurs distinctes dans l'eau de façon à ce que les réverbérations dans la colonne d'eau qu'ils reçoivent soient déphasées de 180 degrés. En combinant la sortie des récepteurs, le document suggère que les réverbérations s'annulent. Le document US 2006/0176775 propose un dispositif permettant de relier une flûte située à une première profondeur à une flûte située à une deuxième profondeur pour assurer l'alignement vertical des flûtes aux première et deuxième profondeurs.
Toutefois, une telle solution est coûteuse et difficile à mettre en oeuvre puisqu'elle impose la mise en place de dispositifs annexes sur les flûtes. Un but de l'invention est de proposer un procédé d'exploration sismique marine plus économique et plus facile à mettre en oeuvre que le procédé décrit dans le document US 2006/0176775.
Un autre but de la présente invention est de fournir un procédé de traitement sismique permettant d'éliminer les fantômes avec un résultat au moins équivalent à celui obtenu avec le système décrit dans le document US 2006/0176775.
PRESENTATION DE L'INVENTION Selon un premier aspect, l'invention a pour objet un procédé d'exploration sismique marine à au moins deux profondeurs différentes au moyen de flûtes comprenant chacune un ensemble d'hydrophones, l'acquisition de données aux différentes profondeurs étant réalisée dans une même région d'exploration, procédé selon lequel, dans ladite région d'exploration, la position des flûtes et des hydrophones lors de l'acquisition des données sismiques à une première profondeur est indépendante de la position des flûtes et des hydrophones lors de l'acquisition des données sismiques à une deuxième profondeur. On entend par région d'exploration , l'aire délimitée par les flûtes lors de l'acquisition des signaux sismiques à l'une des profondeurs, à laquelle on ajoute une marge d' erreur. La marge d'erreur correspond à une surface entourant l'aire délimitée par les flûtes lors de l'acquisition des signaux sismiques à l'une des profondeurs. Cette marge d'erreur dépend notamment d'un facteur de qualité désiré par l'utilisateur, comme il sera expliqué plus en détail dans la suite de la description. Le procédé d'exploration sismique proposé ci-dessus présente l'avantage d'être plus économique que le procédé décrit dans le document US 2006/0176775.
En effet, en se dégageant des contraintes sur l'alignement des flûtes aux première et deuxième profondeurs, le procédé d'exploration décrit ci-dessus permet de s'affranchir de l'utilisation des dispositifs annexes proposés dans le document US 2006/0176775. Dans une variante, l'acquisition des données sismiques est effectuée en deux passages décalés dans le temps. En d'autres termes, dans le cas d'une acquisition de données sismiques à deux profondeurs, l'acquisition des données sismiques à la première profondeur est réalisée lors d'un premier passage, et l'acquisition des données sismiques à la deuxième profondeur est réalisée lors d'un deuxième passage. Cette variante permet de couvrir une zone d'exploration plus étendue et/ou d'améliorer la qualité des traces de sortie obtenues par rapport au procédé d'exploration sismique décrit dans le document US 2006/0176775. En effet, le nombre maximum de flûtes utilisées pour une prospection sismique dépend de la puissance de tractage du bateau utilisé pour effectuer celle-ci. Ainsi, dans le cas d'une exploration sismique, de type over-under , telle que décrite dans le document US 2006/0176775, le nombre de flûtes utilisées à chaque profondeur est égal à la moitié du nombre de flûtes maximum que peut tracter le bateau.
A l'inverse, en effectuant l'acquisition en deux passages tel que proposé, il est possible d'utiliser, pour chacune des deux profondeurs, un nombre de flûtes égale au nombre de flûtes maximum que peut tracter le bateau. Selon un autre aspect, l'invention a pour objet un procédé de traitement de données sismiques marines acquises à au moins deux profondeurs différentes au moyen de flûtes comprenant chacune un ensemble d'hydrophones, l'acquisition de données aux différentes profondeurs étant réalisée dans une même région d'exploration, de manière à obtenir des traces de sortie dans lesquelles des réflexions de surface parasites ont été supprimées, procédé dans lequel les traces de sorties sont obtenues par des combinaisons linéaires des données sismiques dans le domaine fréquentiel, les coefficients de ladite combinaison linéaire étant calculés, pour chaque fréquence, à partir des positions et profondeurs d'acquisition des données sismiques, des positions et de la profondeur auxquelles on souhaite obtenir les traces de sortie, et de l'équation de propagation des ondes dans la couche d'eau située au dessus des hydrophones. Ce procédé de traitement sismique est adapté aux données sismiques acquises à deux profondeurs et dont les positions d'enregistrement aux première et deuxième profondeurs sont indépendantes. On entend par position , des coordonnées en deux dimensions (x, y) dans un plan perpendiculaire à la ligne verticale sur laquelle le mesurage de la profondeur est effectué.
PRESENTATION DES FIGURES
D'autres caractéristiques et avantages des procédés d'exploration et de traitement sismique 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 : - les figures 1 et 2 sont des illustrations schématiques de l'acquisition de données sismiques (ou traces d'entrée), - la figure 3 est une vue de dessus d'hydrophones disposés à deux profondeurs, - la figure 4 est une représentation schématique d'un facteur de qualité pour différentes fréquences en fonction de la position pour laquelle un champ d'ondes montant Uo a été calculé par rapport aux grilles, - la figure 5 illustre des données sismiques acquises à des profondeurs de 9 mètres (Hom) et 18 mètres (Hi8m), ainsi qu'un champ d'ondes montant Uo calculé à partir de ces données sismiques (traces d'entrée), - la figure 6 est une représentation graphique illustrant le champ d'ondes réel Uo, le champ calculé êo ainsi que la différence entre le champ d'ondes réel et le champ d'ondes estimé êo - Uo, - la figure 7 est un organigramme du procédé de traitement sismique.
DESCRIPTION On décrira dans la suite plus en détail les procédés d'exploration sismique et de traitement présentés ci-dessus. Le procédé d'exploration comprend l'acquisition de données sismiques à deux profondeurs. Ces données sismiques 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 de deux jeux de flûtes placés à des profondeurs différentes. Pour chaque profondeur, on peut enregistrer les données sismiques à des instants différents et avec des sources différentes, mais la région d'exploration (zone couverte par un jeu de flûtes au moment de l'acquisition à laquelle on ajoute une marge d'erreur dont la valeur dépend de la fréquence étudiée) doit être la même aux première et deuxième profondeurs. A l'intérieur de cette région d'exploration, les positions des flûtes et des hydrophones aux première et deuxième profondeurs peuvent être différentes. Ainsi, la position des flûtes et des hydrophones lors de l'acquisition des données à la première profondeur est indépendante de la position des flûtes et des hydrophones lors de l'acquisition des données à la deuxième profondeur. Une fois l'acquisition des données sismiques au première et deuxième profondeurs effectuée, un procédé de traitement décrit ci- après assure la suppression des fantômes. La figure 1 illustre schématiquement une étape d'acquisition des données sismiques marines. On a en particulier représenté une masse d'eau 1 s'étendant au dessus d'un fond marin, avec à la surface de l'eau 10 un navire 13 de prospection sismique.
Le navire 13 remorque une ou plusieurs source(s) sismiques 14 destinées à émettre une onde acoustique dans l'eau. La source 14 peut être un canon à air comprimé, un vibrateur marin ou d'autres types de source connus de l'homme de l'art.
Le navire 13 remorque également un premier ensemble de flûtes situé à une première profondeur zl (par exemple 9 mètres) et un deuxième ensemble de flûtes situé à une deuxième profondeur z2 (par exemple 18 mètres). Chaque flûte comprend une pluralité de récepteurs sismiques ou hydrophones 11, 12 destinés à enregistrer les signaux acoustiques émis par la source après leurs réflexions successives sur les interfaces entre couches géologiques - appelées réflecteurs. Le dispositif d'acquisition fonctionne comme suit (cf. figure 2). La vibration sismique émise par la source sismique 14 descend le long de plusieurs trajectoires 15 et est réfléchie aux interfaces entre des matériaux d'impédances acoustiques différentes telles que l'interface 16 - le fond de la mer est référencé 16a. Un champ d'ondes réfléchies 17 remonte vers la surface de l'eau 10 et est enregistré par les hydrophones 12 à la deuxième profondeur, et par les hydrophones 11 à la première profondeur. Les hydrophones 11, 12 aux première et deuxième profondeurs zl, 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 18 dû à la réflexion des ondes à la surface de l'eau 10 se superpose avec le champ d'ondes montant 17 porteur de l'information sur la configuration des couches géologiques.
Les enregistrements comportent donc des pics correspondant à des réflexions de surface, ou réflexions parasites. Le procédé de traitement sismique décrit ci-après permet d'éliminer ces réflexions parasites.
Théorie relative au procédé de traitement Dans ce qui suit, t représente le temps et f les fréquences. On considère un ensemble {E} de n traces d'entrée hydrophones correspondant à un ou plusieurs points de tirs à une même position (Xs, Ys), les coordonnées (x, y, z) de chaque hydrophone étant connue, ainsi que 6(f) correspondant à la signature de la source convoluée par la signature de l'hydrophone et des électroniques d'enregistrement.
On souhaite donc reconstituer les traces de sortie Uo(x, y, t) du champ d'ondes montant pour les m positions de surface d'un ensemble {R}, sachant que les transformées de Fourier FKXKy du champ d'ondes à une profondeur z sont liées à celle de Uo par la relation suivante dérivée de l'équation de propagation des ondes dans un milieu à vitesse constante c (i.e. de l'eau) : Hz(f,kx,ky) =2tsin(27c z.\/(.f/c)2ùkx2ùky2)Uo(f,kx, = Gz(f,k,ky)(LU,k,ky) On a défini ainsi le fantôme dans le domaine des fréquences et des nombres d'ondes par Gz(f,kx,ky) = 2tsin( 27c z.\/(f/c)2ùkx2 ùky2) . Pour une fréquence f donnée, le champ d'ondes montant Uo extrait des données enregistrées est donné par la relation :
1 A = UoH( H+ CN)
n Uo = A Hz avec : A : matrice des coefficients d'une combinaison linéaire calculés pour chaque fréquence HZ : ensemble des données sismiques acquises quelle que soit la profondeur. La matrice de covariance de Uo , dont les éléments diagonaux représentent le niveau de bruit présent dans les données reconstituées est CUo = CUo ùA C* UoH avec : C U H : matrice (m,n) de covariance de Uo aux positions de reconstruction avec les 0 données aux positions enregistrées, CH : matrice (n, n) de covariance du "fantôme " aux positions enregistrées, C N : matrice (n, n) de covariance des niveaux de bruits enregistrées, Cu : matrice (m, m) de covariance de Uo aux positions de reconstruction, et 0 C U H : est la transposée conjuguée (ou transposée hermitienne) de la matrice 0 CUoH . 20 25 La matrice CN ne contient que des termes diagonaux (on ne cherchera pas à modéliser la propagation des bruits) qui sont la puissance estimée du bruit à la fréquence considérée pour les hydrophones. Cette matrice CN peut être obtenue soit par des mesures du bruit (enregistrement sans source d'énergie), soit par l'utilisation de modèles mathématiques de propagation en profondeur de la houle tel que décrit dans la publication de Torgeir Moan, édité dans le livre Dynamic Loading and Design of Structures , de Andreas J. Kappos, en page 180. Si on ne sait rien sur le bruit, on peut éventuellement mettre à zéro les éléments diagonaux. Il faut toutefois ajouter aux éléments diagonaux de CN un terme dit de pré-blanchiment (ou pre-whitening selon la terminologie anglo-saxonne) car à la fréquence zéro (c'est-à-dire le continu) les hydrophones ne peuvent avoir enregistré d'information et le système à résoudre devient singulier. Une fois les traces d'entrée enregistrées, un traitement permet d'éliminer les fantômes. Dans un mode de mise en oeuvre, le traitement comprend une étape 100 consistant à 15 calculer les transformées de Fourier (du domaine temporel au domaine fréquentiel) des traces d'entrée enregistrées par les hydrophones. Puis, pour toutes les fréquences f inférieures ou égale à une fréquence seuil (fm. ou Furax) : - on calcule Km, sachant que Km. = f/c (étape 110), où KTäax représente le nombre d'onde maximum possible pour la frequence f considérée - on calcule les matrices de covariance CH, CUoH, CN (étape 120), - on calcule A, êo ainsi que les éléments diagonaux de Cri. (étape 130), - on calcule le niveau de bruit (étape 140) : si ce niveau de bruit calculé est supérieur à une valeur de bruit seuil, on met à zéro la valeur du champ d'ondes êo estimé à la fréquence f (étape 150). Lorsque la fréquence f devient supérieure à la fréquence seuil (fur. ou Fm.), on calcule la transformée de Fourier inverse des transformées de Fourier des traces de sortie. L'ensemble des données est constitué de toutes les traces d'entrée d'une même position source. Les positions à reconstituer sont celles des hydrophones de la flûte supérieure.
Le théorème d'échantillonnage de Nyquist permet le calcul de la fréquence seuil. En vertu du théorème d'échantillonnage (Nyquist), il n'y a pas de perte d'information pour les longueurs d'ondes supérieures à la moitié de la distance entre deux flûtes adjacentes de la pluralité de flûtes. A partir de cette contrainte, on peut calculer la fréquence seuil Fi. à ne pas dépasser. La relation k= Veau / f entre la longueur d'onde et la fréquence permet de caractériser les longueurs d'ondes présentes à une fréquence donnée. Pour une fréquence f on a donc K.(f) = f / Veau (où Veau représente la vitesse de propagation des ondes dans l'eau). On peut réduire encore les angles d'incidence de Uo en notant que les ondes réfléchies ne peuvent avoir une incidence supérieure à l'angle critique défini par sin(e ) = V eau Vs01 Avantageusement, on peut choisir une fréquence seuil furax fixée par l'utilisateur et inférieure à la fréquence seuil Finax de Nyquist. En effet, on ne cherchera, en général, qu'à compléter les hydrophones des flûtes supérieures dans les basses fréquences. Il suffit donc d'interrompre le procédé à une certaine fréquence seuil donnée par les caractéristiques d'acquisition. Dans ce cas, pour toute fréquence f supérieure à fur , on calcule les traces de sortie à partir des traces d'entrée acquises par les hydrophones des flûtes supérieures (i.e. les plus proches de la surface de l'eau) auxquelles on applique un filtre inverse (désignature). Ceci permet d'éviter de limiter la mise en oeuvre du procédé décrit ci-dessus aux fréquences désirées par l'utilisateur. En effet, pour les hautes fréquences, les procédés classiques de filtrage inverse à partir de données sismiques acquises à une seule profondeur permettent d'obtenir de bon résultats (pour les hautes fréquences) et sont généralement moins coûteux en ressources informatiques et en temps de traitement que les procédés over-under.
Il se peut que le système à résoudre pour un point de tir complet pose des problèmes de nature informatique. On peut avantageusement le découper en plusieurs systèmes plus petits, en prenant la précaution d'inclure tous les hydrophones distants de moins de Z tan(e ) max d'un voisinage considéré : c'est-à-dire en incluant tous les trajets contribuant au voisinage. L'étape de calcul des matrices de covariance se déroule comme suit, avec = j(f /c)2 ùkx2 Soient i, j deux indices correspondants à deux positions de {R} (donc sur la surface). On note r(i,j) la distance (horizontale) entre ces deux points. L'élément (i, j) de la matrice de covariance Cu0 est donnée par (Jo et J1 sont les fonctions de Bessel usuelles): = E(U U ) =(2n)2JJe2tn(kx(xi-xi)+kx(yi-Yj))d1 dky K 27t j)) K r a 47t 2r(i, j) Jl ( Soient p, q deux indices correspondant à deux positions de {H} (donc sous la surface). On note r(p,q) la distance horizontale entre ces deux points. L'élément (p, q) de la matrice de covariance CH est donnée par : CH (p, q) = E(HpHq) = 1 ff6 p6 q GzG e2t~t (kx(xp xg>+ky(Yp Yg d k d k (27c )2 p zq x y CU.(i'j) 1 - K _ (~ v p qJo (cos(21t (zp ù zq)) ù cos(2it (z p + zq) )) Jo (2itkr(p, q))kdk Soient i, p deux indices correspondant à une position de {R} et une position de {H}, on 15 note r(i,p) la distance horizontale entre ces deux points. L'élément (i, p) de la matrice de covariance CUoH est donné par : CUo (i' p) = E(U H p) 1 ff 2tn (kx(xi-xp)+ky(yi-y (2702 JJ 6 pGz pe d kx dky _ ùt ~çù fK v 9J0sin(2it Zp ) Jo (2itkr(i, p))kdk 20 où 6 représente la signature de la source sismique, G représente le fantôme dû aux réflexions parasites de surface. Les deux dernières intégrales n'ont pas de solution analytique mais peuvent être calculées numériquement par un algorithme de type quadrature de Gauss, ou bien en remplaçant les fonctions trigonométriques par un simple développement limité et en résolvant analytiquement ces deux intégrales. La figure 5 est un exemple de résultat présentant le champ d'ondes montant Uo 20 calculé à partir : - de traces d'entrée 21 acquises à une profondeurs de neuf mètres (H9m) et - de traces d'entrée 22 acquises à une profondeur de dix-huit mètres (H18m)• Les données synthétiques illustrées à la figure 5 correspondent à un point diffractant situé à 400 mètres sous la surface de l'eau et enregistré sur deux grilles d'hydrophones 23, 24 situées respectivement à 9 mètres et 18 mètres.
On entend, dans le cadre de la présente, par grille d'hydrophones , l'ensemble des positions des hydrophones à une profondeur. Chaque grille d'hydrophones 23, 24 présente une longueur de 100 mètres et une largeur de 50 mètres. Les deux grilles d'hydrophones 23, 24 sont entrelacées, comme illustré à la figure 3, où les hydrophones 23 de la grille située à 9 mètres de profondeurs sont représentés par des étoiles, et où les hydrophones 24 de la grille située à 18 mètres de profondeurs sont représentés par des losanges. On peut qualifier la reconstitution du champ d'ondes montant Uo par le facteur de qualité tr(CuoùC'U0) P trac o) La figure 4 est une représentation schématique de ce facteur de qualité pour des 20 fréquences de 3 Hz, 6 Hz, 9 Hz et 12 Hz en fonction de la position pour laquelle le champ d'ondes montant Uo a été calculé par rapport aux grilles. Les lignes horizontales 40, 41, 42, 43 permettent de comparer la dégradation de l'extrapolation spatiale du champ d'ondes par rapport à la longueur d'onde de la fréquence considérée. 25 On constate que le facteur qualité est proche de 1 lorsque l'on reste à l'intérieur des deux grilles, puis décroît d'autant plus rapidement que l'on s'éloigne des deux grilles, et que les fréquences considérées augmentent. La marge d'erreur sera déterminée en fonction de la fréquence considérée et du facteur de qualité désiré par l'utilisateur. Ainsi, si l'utilisateur désire un facteur de qualité égale à 0,8 (ligne horizontale 40), alors la largeur de la marge d'erreur sera de 20 mètres à la fréquence f=12Hz, et de 200 mètres à la fréquence f=3Hz.
A l'intérieur des grilles, le champ d'ondes est parfaitement reconstitué, comme illustré à la figure 6 qui représente le champ d'ondes réel 29, le champ d'ondes calculées 30 et la différence 31 entre le champ d'ondes calculé et le champ d'ondes réel.
Le procédé décrit ci-dessus peut être adapté à un traitement en deux dimensions (x, z), c'est-à-dire en estimant que la géologie ne varie pas latéralement. On considère alors les câbles par paires verticales. Il suffit de poser k,, = 0 dans les formules des matrices de covariance
x(f) décrites précédemment, puis de calculer les espérances par une intégration f ...dk. La xmax(f) technique du développement limité du facteur intégré permet de trouver une approximation analytique des covariances.
De préférence, le procédé de traitement décrit ci-dessus sera mis en oeuvre au moyen d'un dispositif de traitement comprenant des moyens pour la mise en oeuvre des étapes du procédé de traitement, tel qu'un ordinateur de type PC comportant une mémoire, une unité de
traitement sur lequel est exécuté un programme d'ordinateur.
Ce programme d'ordinateur comprendra notamment un ou plusieurs algorithme permettant d'exécuter les étapes du procédé décrit précédemment. Bien évidemment, un dernier objet concerne un programme d'ordinateur en tant que tel enregistré sur un support.
Les procédés décrits ci-dessus permettent de limiter les contraintes sur l'alignement des
flûtes aux première et deuxième profondeurs dans le cas d'une acquisition de type overunder .
En particulier :
- il n'est pas nécessaire d'aligner avec précision les flûtes aux première et deuxième profondeurs;
- l'acquisition aux première et deuxième profondeurs peut être effectuée en deux passes ; en d'autres termes, on peut effectuer une acquisition à la première profondeur lors d'un premier passage, et une acquisition à la deuxième profondeur lors d'un deuxième passage.30
Claims (16)
1. Procédé d'exploration sismique marine à au moins deux profondeurs différentes (zl, z2) au moyen de flûtes comprenant chacune un ensemble d'hydrophones (11, 12), l'acquisition de données aux différentes profondeurs étant réalisée dans une même région d'exploration, caractérisé en ce que, dans ladite région d'exploration, la position des flûtes et des hydrophones (11) lors de l'acquisition des données à une première profondeur (zl) est indépendante de la position des flûtes et des hydrophones (12) lors de l'acquisition des données à une deuxième profondeur (z2).
2. Procédé d'exploration selon la revendication 1, dans lequel l'acquisition des données à la première profondeur (zl) est réalisée lors d'un premier passage, et l'acquisition des données à la deuxième profondeur (z2) est réalisée lors d'un deuxième passage.
3. Procédé d'exploration selon l'une des revendications 1 ou 2, comprenant en outre le traitement des données sismiques acquises par les hydrophones (11, 12) aux deux profondeurs (zl, z2) en vue d'obtenir des traces de sortie dans lesquelles des réflexions de surface parasites ont été supprimées, les traces de sorties étant obtenues par des combinaisons linéaires des données sismiques dans le domaine fréquentiel, les coefficients de ladite combinaison linéaire étant calculés, pour chaque fréquence, à partir des positions et profondeurs d'acquisition des données sismiques, de positions et d'une profondeur auxquelles on souhaite obtenir les traces de sortie, et d'une équation de propagation des ondes dans une couche d'eau située au dessus des hydrophones.
4. Procédé d'exploration selon la revendication 3, dans lequel les traces de sortie correspondent à un champ d'ondes montant calculé pour une profondeur égale à la surface de l'eau.
5. Procédé d'exploration selon la revendication 3, dans lequel les positions auxquelles on souhaite obtenir les traces de sortie correspondent aux positions d'acquisition des données sismiques à la première profondeur.
6. Procédé d'exploration selon l'une des revendications 3 à 5, comprenant en outre une étape de calcul d'un niveau de bruit dans les traces de sortie.
7. Procédé selon la revendication 6, dans lequel si le niveau de bruit calculé d'une trace de sortie est supérieur à un niveau de bruit seuil, on affecte à zéro la trace de sortie.
8. Procédé selon la revendication 3 à 7, dans lequel on minimise le niveau de bruit en tenant compte, dans l'étape de calcul des coefficients de la combinaison linéaire, de signatures de la source sismique utilisée lors de l'acquisition des données sismiques, et de caractéristiques des hydrophones.
9. Procédé selon l'une des revendications 3 à 8, dans lequel pour toute fréquence f supérieure à une fréquence seuil fm. prédéterminée par un utilisateur, les traces de sorties sont obtenues à partir des données sismiques, acquises par les hydrophones situés à la profondeur la plus proche de la surface de l'eau, auxquelles est appliqué un filtre inverse.
10. Procédé de traitement de données sismiques marines acquises à au moins deux profondeurs différentes (zl, z2) au moyen de flûtes comprenant chacune un ensemble d'hydrophones (11, 12), l'acquisition de données aux différentes profondeurs (zl, z2) étant réalisée dans une même région d'exploration, en vue d'obtenir des traces de sortie dans lesquelles des réflexions de surface parasites ont été supprimées, caractérisé en ce que les traces de sorties sont obtenues par des combinaisons linéaires des données sismiques dans le domaine fréquentiel, les coefficients de ladite combinaison linéaire étant calculés, pour chaque fréquence, à partir des positions et profondeurs d'acquisition des données sismiques, de positions et d'une profondeur auxquelles on souhaite obtenir les traces de sortie, et d'une équation de propagation des ondes dans une couche d'eau située au dessus des hydrophones.
11. Procédé de traitement selon la revendication 10, dans lequel les traces de sortie correspondent à un champ d'ondes montant calculé pour une profondeur égale à la surface de l'eau.
12. Procédé de traitement selon l'une des revendications 10 et 11, dans lequel les positions auxquelles on souhaite obtenir les traces de sortie correspondent aux positions d'acquisition des données sismiques à la première profondeur.
13. Procédé de traitement selon l'une des revendications 10 à 12, comprenant en outre une étape de calcul d'un niveau de bruit dans les traces de sortie.
14. Procédé de traitement selon la revendication 13, dans lequel si le niveau de bruit calculé d'une trace de sortie est supérieur à un niveau de bruit seuil, on affecte à zéro la trace de sortie.
15. Procédé de traitement selon la revendication 10 à 14, dans lequel on minimise le niveau de bruit en tenant compte des signatures de la source sismique utilisée lors de l'acquisition des données sismiques, et des caractéristiques des hydrophones dans l'étape de calcul des coefficients de la combinaison linéaire.
16. Procédé de traitement selon l'une des revendications 11 à 15, dans lequel, pour toute fréquence f supérieure à une fréquence seuil fm prédéterminée par un utilisateur, les traces de sorties sont obtenues à partir des données sismiques, acquises par les hydrophones situés à la profondeur la plus proche de la surface de l'eau, auxquelles est appliqué un filtre inverse.
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0755279A FR2916540B1 (fr) | 2007-05-25 | 2007-05-25 | 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 |
MX2009012741A MX2009012741A (es) | 2007-05-25 | 2008-05-26 | Procedimiento de exploracion sismica que permite la supresion de fantasmas debido a reflexiones en la superficie del agua y procedimiento para procesar datos sismicos para suprimir estos fantasmas. |
AU2008257573A AU2008257573B2 (en) | 2007-05-25 | 2008-05-26 | Seismic exploration process enabling the suppression of ghosts due to reflections at the water surface, and process for processing seismic data in order to suppress these ghosts |
PCT/EP2008/056409 WO2008145630A1 (fr) | 2007-05-25 | 2008-05-26 | Procédé de prospection sismique permettant la suppression de fantômes dus à des réflexions au niveau de la surface de l'eau, et procédé de traitement des données sismiques afin de supprimer ces fantômes |
BRPI0812104-4A2A BRPI0812104A2 (pt) | 2007-05-25 | 2008-05-26 | Processo de exploração sísmica que possibilita a ocultação de imagens duplas devido a reflexões na superfície da água e processo para processar sísmicos a fim de ocultar essas imagens duplas. |
GB0809583A GB2449564B (en) | 2007-05-25 | 2008-05-27 | Seismic exploration process enabling the suppression of ghostsdue to reflections at the water surface, and process for processing seismic data |
US12/154,767 US8553491B2 (en) | 2007-05-25 | 2008-05-27 | Seismic exploration process enabling the suppression of ghosts due to reflections at the water surface, and process for processing seismic data in order to suppress these ghosts |
NO20093407A NO20093407L (no) | 2007-05-25 | 2009-11-24 | Fremgangsmate for seismiske undersokelser som muliggjor undertrykkelse av "spokelsesrefleksjoner" fra vannoverflaten |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0755279A FR2916540B1 (fr) | 2007-05-25 | 2007-05-25 | 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 |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2916540A1 true FR2916540A1 (fr) | 2008-11-28 |
FR2916540B1 FR2916540B1 (fr) | 2009-08-28 |
Family
ID=38896113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0755279A Expired - Fee Related FR2916540B1 (fr) | 2007-05-25 | 2007-05-25 | 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 |
Country Status (8)
Country | Link |
---|---|
US (1) | US8553491B2 (fr) |
AU (1) | AU2008257573B2 (fr) |
BR (1) | BRPI0812104A2 (fr) |
FR (1) | FR2916540B1 (fr) |
GB (1) | GB2449564B (fr) |
MX (1) | MX2009012741A (fr) |
NO (1) | NO20093407L (fr) |
WO (1) | WO2008145630A1 (fr) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2387729A2 (fr) * | 2009-01-16 | 2011-11-23 | Geco Technology B.V. | Traitement de données sismiques |
US9110191B2 (en) | 2009-03-30 | 2015-08-18 | Westerngeco L.L.C. | Multiple attenuation for ocean-bottom seismic data |
US9733375B2 (en) * | 2011-03-01 | 2017-08-15 | CGGVeritas Services (U.S.) Inc. | Method and device for alternating depths marine seismic acquisition |
US8862408B2 (en) * | 2011-09-28 | 2014-10-14 | Westerngeco L.L.C. | Determining one or more target positions in an acquisition domain for processing survey data |
US9442209B2 (en) * | 2012-07-10 | 2016-09-13 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
US9835745B2 (en) | 2013-01-23 | 2017-12-05 | Cgg Services Sas | Low frequency emission and recording for seismic data acquisition |
EP3004941A1 (fr) * | 2013-05-29 | 2016-04-13 | CGG Services SA | Traitement de données de flûte sismique à capteurs multiples |
CN104536045A (zh) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | 一种基于子波处理的鬼波压制方法 |
CN108020864B (zh) * | 2017-12-29 | 2023-09-12 | 中国电建集团华东勘测设计研究院有限公司 | 一种海上浅层地震勘探的竖直测线***及应用 |
CN112444846B (zh) * | 2019-08-27 | 2023-11-28 | 中国石油天然气集团有限公司 | 地震资料中浅水层影响优化方法及装置 |
CN114167501B (zh) * | 2021-11-04 | 2023-06-30 | 中国海洋大学 | 压制强能量海面鬼波的ftx处理方法、存储介质、处理终端 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6028817A (en) * | 1997-12-30 | 2000-02-22 | Western Atlas International, Inc. | Marine seismic system with independently powered tow vehicles |
US6493636B1 (en) * | 1998-11-05 | 2002-12-10 | Shell Oil Company | Method of marine seismic exploration utilizing vertically and horizontally offset streamers |
US20060250890A1 (en) * | 2005-05-05 | 2006-11-09 | Van Den Berg Peter M | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2757356A (en) * | 1954-01-08 | 1956-07-31 | Texas Instruments Inc | Method and apparatus for canceling reverberations in water layers |
US3689874A (en) * | 1964-11-12 | 1972-09-05 | Manus R Foster | Processing of geophysical data |
US6590831B1 (en) * | 1997-12-30 | 2003-07-08 | Westerngeco L.L.C. | Method and apparatus for controlling and optimizing seismic data acquisition |
GB9800741D0 (en) | 1998-01-15 | 1998-03-11 | Geco As | Multiple attenuation of multi-component sea-bottom data |
GB2401180B (en) * | 2000-12-01 | 2005-05-18 | Pgs Americas Inc | Method and system for evaluating quality of deghosted seismic data |
US6477470B2 (en) * | 2000-12-01 | 2002-11-05 | Pgs Americas, Inc. | Method and system for deghosting |
US6691038B2 (en) * | 2001-06-15 | 2004-02-10 | Westerngeco L.L.C. | Active separation tracking and positioning system for towed seismic arrays |
US6704244B1 (en) * | 2002-10-08 | 2004-03-09 | Pgs Americas, Inc. | Method for acquiring and processing seismic survey data using ocean bottom cables and streamers |
GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
US7800977B2 (en) * | 2004-06-01 | 2010-09-21 | Westerngeco L.L.C. | Pre-stack combining of over/under seismic data |
FR2874270B1 (fr) | 2004-08-12 | 2006-11-17 | Geophysique Cie Gle | Procede d'exploration sismique |
US7499373B2 (en) * | 2005-02-10 | 2009-03-03 | Westerngeco L.L.C. | Apparatus and methods for seismic streamer positioning |
US7372769B2 (en) * | 2005-04-08 | 2008-05-13 | Western Geco L.L.C. | Method and apparatus for adaptive over/under combination of seismic data |
US7489590B2 (en) * | 2005-04-08 | 2009-02-10 | Westerngeco L.L.C. | Method and apparatus for source and receiver side wave field separation |
GB2446825B (en) * | 2007-02-24 | 2009-08-05 | Westerngeco Seismic Holdings | Method for seismic surveying using data collected at different depths |
-
2007
- 2007-05-25 FR FR0755279A patent/FR2916540B1/fr not_active Expired - Fee Related
-
2008
- 2008-05-26 BR BRPI0812104-4A2A patent/BRPI0812104A2/pt not_active IP Right Cessation
- 2008-05-26 MX MX2009012741A patent/MX2009012741A/es active IP Right Grant
- 2008-05-26 WO PCT/EP2008/056409 patent/WO2008145630A1/fr active Application Filing
- 2008-05-26 AU AU2008257573A patent/AU2008257573B2/en not_active Ceased
- 2008-05-27 US US12/154,767 patent/US8553491B2/en not_active Expired - Fee Related
- 2008-05-27 GB GB0809583A patent/GB2449564B/en not_active Expired - Fee Related
-
2009
- 2009-11-24 NO NO20093407A patent/NO20093407L/no not_active Application Discontinuation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6028817A (en) * | 1997-12-30 | 2000-02-22 | Western Atlas International, Inc. | Marine seismic system with independently powered tow vehicles |
US6493636B1 (en) * | 1998-11-05 | 2002-12-10 | Shell Oil Company | Method of marine seismic exploration utilizing vertically and horizontally offset streamers |
US20060250890A1 (en) * | 2005-05-05 | 2006-11-09 | Van Den Berg Peter M | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
Also Published As
Publication number | Publication date |
---|---|
FR2916540B1 (fr) | 2009-08-28 |
US20090016158A1 (en) | 2009-01-15 |
GB2449564A (en) | 2008-11-26 |
AU2008257573A1 (en) | 2008-12-04 |
AU2008257573B2 (en) | 2013-01-17 |
GB0809583D0 (en) | 2008-07-02 |
WO2008145630A1 (fr) | 2008-12-04 |
US8553491B2 (en) | 2013-10-08 |
BRPI0812104A2 (pt) | 2014-11-25 |
MX2009012741A (es) | 2009-12-11 |
NO20093407L (no) | 2009-12-18 |
GB2449564B (en) | 2011-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
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 | |
FR2989788A1 (fr) | Titre non renseigne. | |
FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
CA2733882C (fr) | Procede de surveillance d'un site de stockage geologique de gaz par inversion stratigraphique de donnees sismiques | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
CA2786411A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
FR2974636B1 (fr) | Dispositif et procede d'extrapolation angulaire adaptative pour des rassemblements angulaires tridimensionnels de migration inverse dans le temps | |
FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
FR3016970A1 (fr) | ||
FR2985040A1 (fr) | Device and method for deghosting variable depth streamer data | |
US9594180B2 (en) | Removing ghost reflections from marine seismic data | |
FR2885414A1 (fr) | Procede de deconvolution de signature de source | |
US20140027122A1 (en) | Two Component Source Seismic Acquisition and Source De-Ghosting | |
FR2831961A1 (fr) | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue | |
FR2978252A1 (fr) | Procede et dispositif pour la separation de champs d'ondes dans des donnees sismiques | |
CA2704141A1 (fr) | Methode pour imager une zone cible du sous-sol a partir de donnees de type walkaway | |
FR2990769A1 (fr) | ||
FR2989787A1 (fr) | Traitement de donnees sismiques comprenant une compensation des effets de fantome de source et de recepteur dans une migration inverse dans le temps | |
FR2794249A1 (fr) | Procede de determination d'au moins un temps de parcours multiple sur un lever sismique | |
FR2795181A1 (fr) | Procede de prospection sismique | |
FR2972057A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
EP0752114B1 (fr) | Methode d'exploration d'un milieu a tectonique complexe | |
US9671513B2 (en) | Seismic image acquisition using ghost compensation in Kirchhoff migration | |
US10324208B2 (en) | Premigration deghosting for marine streamer data using a bootstrap approach in Tau-P domain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 10 |
|
ST | Notification of lapse |
Effective date: 20180131 |