FR2795181A1 - Procede de prospection sismique - Google Patents
Procede de prospection sismique Download PDFInfo
- Publication number
- FR2795181A1 FR2795181A1 FR0007855A FR0007855A FR2795181A1 FR 2795181 A1 FR2795181 A1 FR 2795181A1 FR 0007855 A FR0007855 A FR 0007855A FR 0007855 A FR0007855 A FR 0007855A FR 2795181 A1 FR2795181 A1 FR 2795181A1
- Authority
- FR
- France
- Prior art keywords
- data
- seismic
- transformed
- volume
- lines
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
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/003—Seismic data acquisition in general, 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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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)
- Geophysics And Detection Of Objects (AREA)
Abstract
L'invention concerne le domaine de l'acquisition et du traitement de données sismiques.Dans ce procédé, on active une source sismique pour que des ondes sismiques pénètrent dans le sous-sol. Les ondes sismiques réfléchies à au moins une interface sont reçues le long de plusieurs récepteurs pour fournir des données sismiques qui sont ensuite traitées pour fournir au moins une ligne de récepteurs additionnelle de données sismiques estimées proche des lignes de récepteurs.Domaine d'application : prospection sismique, etc.
Description
L'invention a trait au domaine de l'acquisition et du traitement de données sismiques. L'invention concerne plus particulièrement un procédé d'utilisation de techniques d'extrapolation dans deux dimensions pour la prédiction de lignes de données sismiques au-delà de lignes existantes de données sismiques et l'interpolation de lignes de données sismiques entre des lignes existantes de données sismiques.
Des prospections ou arpentages en sismique marine en 3D font appel au remorquage de dispositifs de capteurs sismiques allongés suivant un couloir contrôlé appelé andain. Les andains sont parcourus de façon répétée pour augmenter le degré de couverture dans une direction transversale à la direction des dispositifs de capteurs. Pour réduire le coût de l'acquisition, il est souhaitable d'éviter un chevauchement des andains successifs et, si possible, de ne pas laisser d'intervalles entre les andains successifs.
La présente invention peut être appliquée à toute forme d'opération sismique, qu'elle soit terrestre ou en mer. Cependant, pour des raisons de commodité et uniquement à titre d'exemple nullement limitatif, la description portera sur un arpentage sismique marin.
Dans des opérations marines en 3D, un navire de prospection sismique remorque le long d'un andain plusieurs câbles de flûtes sismiques parallèles suivant une ligne souhaitée d'arpentage, les câbles étant immergés de quelques mètres au-dessous de la surface de l'eau. Le nombre de câbles qui couvrent un andain dépend uniquement des capacités mécaniques et fonctionnelles du navire remorqueur. Il peur y avoir six de ces câbles, ou plus, espacés les uns des autres d'environ 50 à 100 mètres. Les câbles respectifs peuvent atteindre 8 à 12 km de longueur.
Chaque câble de flûte comprend habituellement 120 ou plus de 120 groupes de détecteurs sismiques espacés. Chaque groupe est constitué d'un ou plusieurs détecteurs individuels interconnectés, desservant chacun un canal unique de données. L'espacement des groupes sur la longueur du câble est de l'ordre de 25 à 50 mètres. Les détecteurs sismiques sont des transducteurs qui perçoivent l'activité mécanique due à des champs d'ondes acoustiques réfléchies et qui convertissent cette activité en signaux électriques ayant des caractéristiques représentatives de l'intensité, du positionnement temporel et de la polarité de l'activité acoustique, comme cela est bien connu dans la technique. Les détecteurs sont couplés fonctionnellement à des dispositifs de stockage et de données de tout type souhaité.
Une source acoustique telle qu'un dispositif de canons à air est remorquée dans l'eau par le navire à proximité de l'extrémité avant de l'andain de câbles des flûtes sismiques. Pendant que le navire parcourt la ligne d'arpentage, la source est déclenchée (activée) à des intervalles choisis dans l'espace, communément égaux à l'intervalle des groupes ou intertrace. En supposant que le navire se déplace à une vitesse constante telle que 4 à 6 noeuds, la source peut être commodément déclenchée à des intervalles de temps choisis, par exemple toutes les 5 à 13 secondes ou plus. Le champ d'onde émis par la source descend de manière à être réfléchi par des formations terrestres sous-marines d'où le champ d'onde revient par réflexion vers la surface de l'eau où -le champ d'onde réfléchi est reçu par les détecteurs et converti en signaux électriques. Les signaux électriques détectés sont transmis à tous moyens bien connus d'enregistrement et de traitement de signaux pour procurer un modèle physique du sous-sol.
La figure 1 des dessins annexés et décrits ci-après est une vue en plan d'une acquisition de prospection sismique 3D. Un premier andain 3D 13 de six câbles sismiques parallèles Al à A6 est remorqué dans une masse d'eau par un navire 14. Les câbles sont habituellement configurés de façon a être espacés d'une distance d. On doit comprendre que six câbles ne sont représentés qu'à titre d'exemple, et qu'un nombre de câbles plus grand ou plus petit pourrait être utilisé lors de l'acquisition. On doit en outre comprendre que, si l'on considère des opérations terrestres, les câbles pourraient être remorqués par un ou plusieurs camions ou pourraient être déployés par des câblières utilisant des techniques de couverture multiple d'une manière bien connue dans l'industrie sismique. Des signaux provenant des dispositifs de câbles respectifs Al à A6 sont transmis par l'intermédiaire d'un collecteur 20 de signaux de données à un processeur 22 de tout type bien connu, installé à bord du navire 14, lequel collecteur est relié fonctionnement au processeur 22 par des lignes de tête électriques 16 et 18. Une source acoustique discrète SL est remorquée par le navire 14 à proximité de l'extrémité avant de l'andain 13, sensiblement au centre de l'andain. On peut utiliser, si cela est souhaité, plus d'une source discrète telle que les sources SL' et SL <B>'</B>, décalées par rapport à l'axe central.
Un second andain 3D 13' de six dispositifs de câbles sismiques parallèles Al' à A6' peut être obtenu en remorquant les câbles dans une position qui est décalée latéralement de celle du premier andain 13. Ceci peut être réalisé par un navire 14'. On doit comprendre que le second andain pourrait être obtenu ultérieurement au premier andain, afin que le navire 14' puisse être le navire 14 lui-même, ou bien qu'il pourrait être acquis simultanément à l'acquisition du premier andain 13 en utilisant un second navire. L'invention est décrite ici à titre non limitatif comme si les différents andains étaient acquis à des temps différents en utilisant le même navire.
Les spécialistes de la technique reconnaîtront que pour minimiser le coût d'acquisition de données sismiques, la distance D entre les deux andains 13 et 13' devrait, au minimum, être égale à la distance d entre les câbles. Si tel n'est pas le cas, il en résulte une duplication inutile des données. Les spécialistes de la technique reconnaîtront aussi que si la distance D est supérieure à la distance d, les coûts d'acquisition de données pourraient être réduits si les données contenues dans l'intervalle entre les andains pouvaient être complétées par une simulation de données acquises au moyen de câbles se trouvant dans l'intervalle. Le terme "intervalle" n'est pas limité ici à des intervalles dans l'espace dans la géométrie d'acquisition et il entend inclure en particulier un intervalle dans l'irradiation du sous-sol par l'acquisition sismique. Le terme "insonification" a parfois été utilisé pour décrire l'irradiation par des ondes sonores. La qualité d'une image traitée qui est basée en partie sur de telles données simulées dépend de la façon dont la simulation est réalisée. Même si un seul andain de données est acquis, il serait souhaitable de simuler des câbles supplémentaires de données qui sont simplement extrapolés en partant d'un andain existant de données. Un autre problème qui est rencontré dans l'acquisition de données sismiques 3D est celui du repliement du spectre, dans lequel, par suite d'un échantillonnage incorrect de données, certaines fréquences et/ou certains pendages dans le sous-sol ne peuvent pas faire l'objet d'images convenables. Une interpolation de données entre des câbles existants est une solution possible au problème. L'invention aborde la nécessité de ces extrapolation et interpolation de câbles de données sismiques.
L'invention concerne un procédé pour l'extrapolation de câbles de données sismiques en utilisant de petites fenêtres en alignement dans des câbles existants. Un filtre d'erreur de prédiction 2D qui peut prédire les données dans des directions vers l'avant et vers l'arrière est conçu pour chaque tranche de fréquence d'une fenêtre d'espace. Un filtre à prédiction est obtenu à partir du filtre d'erreur de prédiction et est appliqué à un câble à un bord de la fenêtre d'espace pour prédire un premier câble manquant. En répétant ce processus en utilisant des fenêtres en ligne se chevauchant, on peut obtenir des fenêtres en ligne se chevauchant des câbles extrapolés. Grâce à une pondération convenable des fenêtres en ligne des câbles extrapolés, on extrapole une longueur de câble complète. Le processus peut être répété en utilisant le premier câble extrapolé dans la dérivation pour donner des câbles extrapolés supplémentaires. Un processus de prédiction similaire est utilisé pour l'interpolation de câbles manquants entre des andains existants de données sismiques 3D.
Les particularités nouvelles qui sont considérées comme étant caractéristiques de l'invention, à la fois en ce qui concerne l'organisation et les procédés de mise en #uvre, ainsi que des objets et avantages de celle-ci seront mieux compris à la lecture de la description détaillée suivante en regard des dessins annexés à titre nullement limitatif et sur lesquels la figure 1 est une illustration schématique d'une acquisition de données sismiques marines illustrant deux andains d'une acquisition 3D ; les figures 2A et 2B montrent l'extrapolation de câbles à partir d'un ensemble existant de câbles sismiques ; la figure 3 illustre l'utilisation de fenêtres spatiales se chevauchant pour obtenir un câble extrapolé ; la figure 4 montre l'application d'un filtre à prédiction pour obtenir des emplacements de données sur un câble extrapolé ; les figures 5A et 5B montrent la dérivation d'un filtre à prédiction 2D ; la figure 6 est un organigramme de certaines des étapes de la présente invention ; la figure 7 montre l'interpolation de câbles dans un andain existant d'un arpentage sismique ; et la figure 8 montre l'interpolation de traces individuelles à l'emplacement d'un récepteur à partir de traces situées sur un andain sismique existant. En référence à présent à la figure 2A, celle-ci montre une illustration schématique d'une partie d'emplacements sur un arpentage ou une prospection pour une acquisition sismique dans trois dimensions, ou 3D. Plusieurs lignes sismiques C1, C2, C3 ... CNY sont représentées et plusieurs emplacements de récepteurs sismiques R1, R2, R3 ... RNX sur chacune des lignes sismiques sont indiqués. Pour simplifier l'illustration, les emplacements des récepteurs sont indiqués sur un seul câble. Comme décrit dans le contexte ci-dessus, les lignes sismiques dans un arpentage ou une prospection sismique marine classique sont constituées de câbles connectant les récepteurs. Cependant, l'invention englobe des arpentages ou prospections sismiques dans lesquels la prospection peut être exécutée à terre, ou des arpentages dans lesquels les récepteurs le long de lignes sismiques transmettent par télémétrie des données reçues. Dans un souci de commodité, le terme "câble" est utilisé dans la description suivante pour englober tous ces agencements de lignes sismiques le long desquels se trouvent des récepteurs.
La pluralité de câbles C1, C2, C3 ... CNY et la pluralité de récepteursRl, R2, R3 ... RNX définissent une fenêtre spatiale d'emplacements ayant pour dimensions nx x ny dans les directions x et y. Le nombre nX d'emplacements de récepteurs dans le câble utilisés dans la fenêtre est en général inférieur au nombre de canaux de données disponibles sur le câble . ce dernier nombre est habituellement de 120 ou plus. En utilisant un procédé d'extrapolation faisant appel à un filtre à prédiction décrit ci-dessous en référence à la figure 9, cette fenêtre spatiale de nx x ny emplacements peut être utilisée pour définir une fenêtre de longueur nx d'un câble extrapolé désigné par CNY+l comme indiqué sur la figure 2B. En utilisant une version inversée conjuguée du câble à prédiction, on peut également obtenir un câble désigné par Co comme indiqué sur la figure 2B, c'est à dire que le câble extrapolé peut être sur un côté ou sur l'autre du groupe existant de câbles C1, CZ, C3 ... CNY.
En référence à présent à la figure 3, celle-ci montre un certain nombre de fenêtres spatiales<B>101,</B> 102 et 103 se chevauchant. Chacune des fenêtres est de longueur nx et de largeur nY, comme sur les figures 2A et 2B ci-dessus. En utilisant des données provenant de la fenêtre 101, on obtient une fenêtre de longueur nx sur un câble extrapolé comme décrit ci-dessus. Cette fenêtre extrapolée est désignée par la référence numérique 111. Similairement, les fenêtres 102 et 103 sur les fenêtres spatiales sont utilisées pour produire des fenêtres extrapolées 112 et <B>113.</B> Ceci est poursuivi sur toute la longueur des câbles (non représentés) pour produire des fenêtres extrapolées qui couvrent la longueur des câbles.
Comme on le voit sur la figure 3; il existe une région de chevauchement 120a entre les fenêtres extrapolées 111 et 112, et une région de chevauchement 120b entre les fenêtres extrapolées 112 et 113. En utilisant des procédés connus, on peut combiner les données des câbles extrapolés se trouvant dans les régions de chevauchement 120a et 120b, et, en répétant le processus pour toutes fenêtres supplémentaires, on peut extrapoler un câble entier de données. Dans une forme appréciée de réalisation de l'invention, une fonction de pondération arrondie est appliquée aux données associées aux fenêtres<B>111</B> et 112 dans la zone de chevauchement 120a et les données extrapolées pondérées provenant des fenêtres 111 et 112 sont sommées. Des techniques d'arrondissement communément utilisées, comme cela est connu des spécialistes de ce domaine, comprennent un arrondissement linéaire ou une fonction cosinus.
Comme indiqué précédemment, le nombre de récepteurs dans une prospection est de 120 ou plus, six récepteurs étant un nombre typique pour la pluralité de câbles nY. Dans une forme appréciée de réalisation de l'invention, l'extrapolation pour un câble additionnel est réalisée par l'application d'un filtre à prédiction aux données sismiques reçues sur les câbles nY. La fenêtre d'espace décrite ci-dessus en référence à la figure 3 a une longueur nx qui est avantageusement égale à dix ou plus de dix intervalles de groupes de récepteurs. Le filtre à prédiction est dérivé à l'intérieur de chaque fenêtre d'espace comme décrit ci-dessous en référence aux figures 4 à 6.
La figure 4 montre la manière dont des points sur un câble extrapolé sont obtenus lorsqu'un filtre à prédiction a été dérivé comme décrit ci-dessous. On représente par des points plusieurs emplacements de récepteurs sur des câbles existants 201a, 201b, 201c et 201d. Uniquement à titre illustratif, le filtre à prédiction est un filtre du type 5 x 4 pour l'extrapolation de données. Le filtre à prédiction utilise des données provenant de câbles existants dans de pavé 5 x 4 désigné en 210. Ces données comprennent des points précédemment déterminés désignés par des x le long d'un câble extrapolé 201e. L'application du filtre 5 x 4 aux données à l'intérieur du pavé 210 donne une trace de sortie en un emplacement 212 désigné par le cercle ouvert. Le processus est ensuite répété en déplaçant d'une trace le pavé 210 dans la direction x, en utilisant les traces de données désignées par un signe + sur les câbles existants et la trace à l'emplacement 212 pour donner une nouvelle trace de sortie en un emplacement 214. Ce processus est répété pour toutes les traces à l'intérieur d'une fenêtre spatiale comme défini ci-dessus dans la description de la figure 3. Les spécialistes de la technique reconnaîtront la présence d'un "effet d'extrémité" au commencement du processus d'extrapolation dans la génération des quelques premières traces sur un câble extrapolé.
Le processus décrit ci-dessus en référence à la figure 4 peut être répété en utilisant un câble extrapolé en tant que partie de l'entrée pour produire un deuxième câble extrapolé et des câbles extrapolés supplémentaires de données sismiques prédites.
L'extrapolation de données de câble sur l'autre côté des données existantes, c'est à dire sur la gauche sur la figure 4, est réalisée en utilisant un filtre qui est le conjugué inversé et complexe du filtre utilisé pour l'extrapolation vers la droite.
La conception d'un filtre 2D est décrite en référence aux figures 5A et 5B. I1 est montré une fenêtre d'espace de dimensions nx x nY. Caque point dans la fenêtre d'espace est une trace de données sismiques définissant un volume dans l'espace (x, y, t). Dans une forme appréciée de réalisation de l'invention, les données sont soumises à une transformation de Fourier par rapport au temps, donnant un volume de données dans un espace (x, y, co), où co est la pulsation 2nf. Un filtre séparé est dérivé pour chaque fréquence f, ou, de façon équivalente, co.
En référence à présent à la figure 5A, une région de dimensions lx x 1y désignée par 300a montre les dimensions du filtre. La représentation d'un filtre à erreur de prédiction à deux dimensions de dimensions lx x 1y sur une fenêtre d'espace à deux dimensions de format nx x ny est réalisée en inscrivant les équations de convolution énonçant que le filtre à prédiction en avant annihile des points de données de (1X, 1Y) à un point (nX, nY) comme indiqué par la flèche 301 sur la figure 5A. Similairement, la version inversée conjuguée du même filtre peut annihiler des données dans le sens inverse du point (nX-lX+l, nY-lY+l) au point (1, 1) comme indiqué par la flèche 301a sur la figure 5B.
En désignant le filtre à erreur de prédiction e sous la forme e =<B>[El,,, E2,1,</B> . . . , ELX,1, <B>E1,21</B> E2, z. . . . EIX, z, E1,LY, F,z,LY, . . . EJX, LY] T (1) où le []T désigne la transposée, l'équation du filtre est écrite sous la forme Xe = n où e est le filtre à erreur de prédiction, un vecteur de longueur lx lY, X est une matrice de données de format 2 (nx-lx+l) (nY-lY+l) par lx 1Y, et n est un vecteur bruit de longueur 2 (nx-lx+l) (nY-lY+l) .
Le filtre à erreur de prédiction e est résolu en utilisant le critère selon lequel le vecteur bruit et le signal ne peuvent être en corrélation qu'avec un retard 0, c'est à dire XT+ Xe = [l, 0, 0, 0, ... 0]T (2) où le signe + désigne un conjugué complexe. La matrice XT+ X est une matrice carrée de dimensions lx 1Y. L'équation (2) est résolue pour les coefficients du filtre e. Une condition nécessaire pour que l'équation (2) ait une solution unique est que le nombre d'équations données par l'équation (1) soit supérieur au nombre d'inconnues, c'est à dire 2 (nx-lx+1) (nY-lY+l) > lx 1Y.
Les spécialistes de la technique reconnaîtront que la forme de réalisation décrite ci-dessus pour la dérivation du filtre est exempte d'effets de bord, c'est à dire que le filtre est calculé sans parcourir les bords de la fenêtre d'espace. L'autocorrélation résultant de ce type de calcul est appelée une autocorrélation non encadrée. Ce procédé est également appelé procédé à covariance modifiée. Le procédé n'effectue aucune supposition concernant des données en dehors de la fenêtre d'espace utilisée dans la dérivation du filtre. Dans une variante de réalisation de l'invention, une autocorrélation dite encadrée peut être utilisée. Une telle autocorrélation encadrée donne un filtre qui est moins sensible au bruit de préblanchiment que le procédé non encadré, c'est à dire un filtre plus robuste. Cependant, une autocorrélation non encadrée est plus efficace en ce qui concerne la prédiction que l'approche encadrée. La description ci-dessus de la dérivation et l'application du filtre 2D ne sont indiquées qu'à titre d'exemple et d'autres procédés connus des spécialistes de la technique pourraient être utilisés pour la dérivation et l'application du filtre.
La figure 6 montre un organigramme illustrant les étapes principales de la présente invention. Un andain d'entrée de données sismiques 400 comportant plusieurs câbles est cloisonné en fenêtres spatiales comme décrit ci- dessus en référence à la figure 3. Une fenêtre initiale de données est sélectionnée pour un traitement en 404. Les données sont transformées du domaine (x, y, t) au domaine (x, y, w) en utilisant un procédé tel qu'une Transformation de Fourier Rapide (TFR). Ceci donne un nombre de tranches de fréquence sur une plage de fréquence 406 pour la fenêtre d'espace sélectionnée 404. Une tranche de fréquence est sélectionnée en 408 pour une dérivation du filtre comme décrit ci-dessus en référence à la figure 5. Le filtre est appliqué en 412 et une tranche de fréquence est délivrée en sortie en 420. Un contrôle est effectué pour déterminer si toutes les tranches de fréquence ont été traitées en 414. Si la réponse est non, on passe en 416 à la tranche de fréquence suivante et on revient en 408. Si le contrôle effectué en 414 indique que toutes les tranches de fréquence ont été traitées, les tranches filtrées de sortie 420 sont soumises à une transformation de Fourier inverse en 422 pour donner en sortie une fenêtre sur une trace extrapolée 430. Un contrôle est réalisé pour déterminer si toutes les fenêtres d'espace ont été traitées en 424. Si la réponse est non, le traitement passe à la fenêtre suivante 426 et revient à l'étape 404. Si toutes les fenêtres ont été traitées, les fenêtres extrapolées en sortie sont alors combinées comme décrit ci-dessus en référence à la figure 3. On doit comprendre que le traitement décrit en regard de la figure 6 peut également être appliqué à l'extrapolation de câbles dans chaque direction, comme décrit précédemment.
La présente invention peut également être utilisée, avec une légère modification, pour le problème de l'interpolation de câbles en des emplacements entrelacés à l'intérieur d'un andain de câbles existants. Dans une forme appréciée de réalisation de l'invention, on utilise le procédé décrit dans le brevet des États Unis d'Amérique n 5 677 892 délivré aux mêmes inventeurs et ayant la même Cessionnaire. Ce brevet, dont le contenu est incorporé ici complètement à titre de référence, traite du problème de l'interpolation de trace spatiale sans repliement dans le domaine F-K d'un volume de données dans trois dimensions. Les traces le long de la pluralité des câbles sismiques sont utilisées pour définir un volume de données dans trois dimensions dans lequel la première valeur de coordonnée dans l'espace est la position du récepteur et la deuxième valeur de coordonnée dans l'espace est le numéro du câble. Pour obtenir l'ensemble de données sans repliement (interpolées), des termes L-1 zéro sont insérés dans chaque direction (x et y) et un filtre à interpolation qui comble ces traces mortes est recherché. Un filtre est calculé à partir des (1/L)ièmes fréquences temporelles inférieures des versions "à atténuation nulle" et "à atténuation nulle, masquage nul" des données d'origine. Dans le contexte de la présente invention, le volume de données 3D d'entrée comprend plusieurs câbles (la valeur de coordonnée y du volume de données du brevet précité), chacun des câbles ayant plusieurs traces (la valeur de coordonnée x du volume de données du brevet précité). Le résultat final de l'application du procédé du brevet précité est de produire des câbles interpolés dans lesquels des récepteurs sont également interpolés. La présente invention laisse à l'écart les positions des récepteurs interpolés le long des câbles d'origine et interpolés pour produire un ensemble de câbles interpolés ayant le même espacement de récepteur que dans les câbles d'origine.
Une variante de réalisation de la présente invention utilise une modification du procédé du brevet des États Unis d'Amérique n 5 617 372 délivré aux mêmes inventeurs que ceux de la présente demande et ayant la même Cessionnaire, qui enseigne un procédé d'interpolation de traces sismiques dans un câble unique. Le contenu du brevet n 5 617 372 précité est entièrement incorporé ici à titre de référence et son procédé est utilisé pour l'interpolation de câbles entre des câbles existants de données sismiques.
Ceci est illustré sur la figure 7 sur laquelle six câbles de données 501a, 501b, 501c, 501d, 501e et 501f sont représentés. Dans une variante de réalisation de l'invention, des câbles additionnels de données 503a, 503b, 503c, 503d et 503e sont obtenus par un processus d'interpolation de données à partir des câbles existants.
Le procédé d'obtention d'une trace unique sur un câble interpolé est décrit en référence à une station 511 de récepteur sur l'andain de données. En référence à présent à la figure 8, les traces sismiques dans une position 511 sur les six câbles 501a ... 501f sont désignées par 511a ... 511f. L'axe vertical de la figure 8 est l'axe des temps. En utilisant le procédé du brevet n 5 617 372 précité, on interpole des traces additionnelles en des emplacements situés entre les emplacements sur le câble existant pour donner des traces 513a ... 513e. Le procédé de ce dernier brevet comprend une transformation des données à partir des coordonnées d'espace-temps (y et t dans le présent contexte) en coordonnées de fréquence- nombre d'onde (ky et ti) dans le présent contexte), l'interpolation étant réalisée dans le domaine transformé. Ceci est suivi par une transformation en retour vers le domaine espace-temps. Ce processus d'interpolation est ensuite répété pour d'autres emplacements de station sur l'andain existant, ce qui permet de produire des câbles entiers de données interpolées 503a ... 503e.
Il va de soi que de nombreuses modifications peuvent être apportées au procédé décrit et représenté sans sortir du cadre de l'invention.
Claims (15)
1. Procédé de prospection sismique, caractérisé en ce qu'il comprend (a) l'activation d'une source sismique (SL) pour amener des ondes sismiques à se propager sous la surface du sol ; (b) la réception de réflexions des ondes sismiques se propageant, provenant d'au moins une interface réfléchissante dans le sous-sol suivant plusieurs lignes de récepteurs, chaque ligne de récepteurs comprenant une pluralité de récepteurs sismiques (R1 ... RNY), pour fournir des données sismiques reçues le long de la pluralité de lignes de récepteurs ; et (c) le traitement des données sismiques reçues pour produire au moins une ligne de récepteurs additionnelle (CNY+1) de données sismiques estimées, proche de la pluralité de lignes de récepteurs.
2. Procédé selon la revendication 1, caractérisé en ce que la pluralité de récepteurs sur chaque ligne de récepteurs communique des données à un processeur au moyen d'au moins l'un de (i) un câble connectant la pluralité de récepteurs sur chaque ligne de récepteurs, et (ii) une télémétrie.
3. Procédé selon la revendication 1, caractérisé en ce que la ou chaque ligne de récepteurs additionnelle comporte des données extrapolées et le traitement des données reçues comprend en outre (i) la sélection d'au moins un sous-ensemble de récepteurs sismiques le long de chacune de la pluralité de lignes sismiques pour définir au moins une fenêtre d'espace ; et (ii) la dérivation d'au moins un filtre à prédiction en utilisant les données sismiques reçues à l'intérieur d'au moins une fenêtre d'espace pour prédire les données sismiques reçues à l'intérieur de la fenêtre d'espace, le filtre à prédiction ayant une dimension inférieure à celle de la ou de chaque fenêtre d'espace.
4. Procédé selon la revendication 3, caractérisé en ce qu'il comprend en outre l'application du ou de chaque filtre à prédiction à au moins une partie des données reçues à l'intérieur de la ou de chaque fenêtre d'espace pour produire au moins une fenêtre de données sur la ou chaque ligne de récepteurs additionnelle.
5. Procédé selon la revendication 3, caractérisé en ce qu'il comprend en outre l'application du ou de chaque filtre à prédiction à au moins une partie des données reçues à l'intérieur de la ou de chaque fenêtre d'espace pour produire au moins une fenêtre de données sur la ou chaque ligne de récepteurs additionnelle.
6. Procédé selon la revendication 3, dans lequel la ou chaque fenêtre d'espace comprend au moins deux fenêtres d'espace, le procédé étant caractérisé en outre en ce qu'il comprend une combinaison des données estimées sur la ou chaque ligne additionnelle en utilisant une pondération arrondie.
7. Procédé selon la revendication 3, dans lequel la ou chaque ligne de récepteurs additionnelle comporte en outre au moins une ligne de récepteurs additionnelle, le procédé étant caractérisé en ce qu'il comprend en outre le traitement des données sismiques reçues et des données estimées sur la ou chaque ligne de récepteurs pour fournir des données estimées sur la ou chaque ligne de récepteurs additionnelle.
8. Procédé selon la revendication 3, caractérisé en ce que le traitement des données sismiques reçues comprend en outre l'application d'une Transformation de Fourier Rapide aux données sismiques reçues pour fournir plusieurs tranches de fréquence des données sismiques reçues.
9. Procédé selon la revendication 8, caractérisé en ce que le ou chaque filtre à prédiction comprend plusieurs filtres à prédiction, chacun des filtres à prédiction correspondant à l'une des tranches de fréquence.
10. Procédé selon la revendication 3, caractérisé en ce que la dérivation du ou de chaque filtre à prédiction comprend en outre l'utilisation d'un procédé à covariance modifiée.
11. Procédé selon la revendication 1, caractérisé en ce que la ou chaque ligne de récepteurs additionnelle s'étend entre au moins deux de la pluralité de lignes de récepteurs et les données sismiques estimées comprennent des données sismiques interpolées.
12. Procédé selon la revendication 11, caractérisé en ce que le traitement des données comprend en outre (i) la définition d'un volume de données en trois dimensions dans lequel une première dimension spatiale est la pluralité de câbles, une deuxième dimension spatiale est la pluralité de récepteurs sismiques et la troisième dimension est le temps ; (ii) une atténuation nulle des données dans le volume de données en trois dimensions dans les première et deuxième dimensions spatiales ainsi que dans le temps pour produire un volume de données à atténuation nulle ; (iii) l'exécution d'une opération de masquage sur le volume de données à atténuation nulle pour produire un volume de données à atténuation nulle, masquage nul ; (iv) la définition d'un domaine transformé dans lequel une première dimension transformée est un nombre d'onde dans la première dimension spatiale, la deuxième dimension transformée est un nombre d'onde dans la deuxième dimension spatiale et la troisième dimension transformée est une fréquence ; (v) la transformation du volume de données à atténuation nulle pour produire un premier volume de données transformées dans le domaine transformé ; (vi) la transformation de l'ensemble de données à atténuation nulle, masquage nul pour produire un deuxième volume de données transformées dans le domaine transformé ; (vii) la définition d'un opérateur d'interpolation dans le domaine transformé en tant que rapport du premier volume de données transformées et du deuxième volume de données transformées et en ne retenant que les composantes de fréquence inférieures (d'un facteur 1/L) ; (viii) l'insertion de traces L-1 zéro dans chaque dimension spatiale du volume de données dans trois dimensions pour former un volume de données à insertion zéro ; (ix) la transformation du volume de données à insertion zéro en le domaine transformé pour fournir un volume de données à insertion zéro, transformé ; (x) l'application de l'opérateur d'interpolation de volume de données à insertion zéro, transformé pour fournir un volume de données interpolées, transformé ; (xi) la transformation inverse du volume de données interpolées, transformé pour fournir un volume de données interpolées ; et (xii) le rejet de traces du volume de données interpolées dans la deuxième dimension spatiale.
13. Procédé selon la revendication 12, caractérisé en ce que L = 2.
14. Procédé selon la revendication 11, caractérisé en ce que le traitement des données comprend en outre . (i) le tri des données pour obtenir des données triées en au moins un emplacement commun d'une station de récepteur ; (ii) la transformation des données triées en l'emplacement ou chaque emplacement commun de station de récepteur en un domaine de fréquence-nombre d'onde pour fournir des données transformées ; (iii) le filtrage des données transformées pour obtenir des données transformées filtrées ; et (iv) la transformation des données transformées filtrées en un domaine espace-temps.
15. Procédé selon la revendication 12, caractérisé en ce que l'application de l'opérateur d'interpolation comprend en outre une limitation de bande le long des axes des nombres d'onde des récepteurs.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US33711699A | 1999-06-21 | 1999-06-21 |
Publications (1)
Publication Number | Publication Date |
---|---|
FR2795181A1 true FR2795181A1 (fr) | 2000-12-22 |
Family
ID=23319197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0007855A Withdrawn FR2795181A1 (fr) | 1999-06-21 | 2000-06-20 | Procede de prospection sismique |
Country Status (5)
Country | Link |
---|---|
US (1) | US6292755B2 (fr) |
AU (1) | AU777935B2 (fr) |
CA (1) | CA2311822A1 (fr) |
FR (1) | FR2795181A1 (fr) |
GB (1) | GB2353096B (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008117121A2 (fr) * | 2006-11-06 | 2008-10-02 | Magnitude Spas | Système, procédé et produit de programme d'ordinateur pour accumuler des données de bruit sismique pour analyser des événements sismiques |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6510390B1 (en) * | 1999-10-07 | 2003-01-21 | Westerngeco, L.L.C. | 3-D seismic trace extrapolation and interpolation |
US6956084B2 (en) * | 2001-10-04 | 2005-10-18 | Bridgestone Corporation | Nano-particle preparation and applications |
GB2394045B (en) * | 2002-10-11 | 2006-07-26 | Westerngeco Seismic Holdings | Method and apparatus for positioning of seismic sensing cables |
US6745129B1 (en) | 2002-10-29 | 2004-06-01 | The University Of Tulsa | Wavelet-based analysis of singularities in seismic data |
WO2004100542A1 (fr) * | 2003-05-05 | 2004-11-18 | Thomson Licensing S.A. | Procede d'enregistrement maitrise dans un dispositif interconnecte conforme a la norme ieee-1394 |
GB2414299B (en) * | 2004-05-21 | 2006-08-09 | Westerngeco Ltd | Interpolation and extrapolation method for seismic recordings |
US7739051B2 (en) * | 2004-07-14 | 2010-06-15 | Compagnie Generale De Geophysique | Method for determination of diffractor locations at sea bottom for the purpose of attenuating such energy |
US7523003B2 (en) * | 2006-07-12 | 2009-04-21 | Westerngeco L.L.C. | Time lapse marine seismic surveying |
WO2009066057A2 (fr) * | 2007-11-19 | 2009-05-28 | Geco Technology B.V. | Procédé de traitement de données sismiques pour une atténuation multiple de surface |
US7646672B2 (en) * | 2008-01-18 | 2010-01-12 | Pgs Geophysical As | Method for wavefield separation in 3D dual sensor towed streamer data with aliased energy in cross-streamer direction |
US8861306B2 (en) * | 2008-07-05 | 2014-10-14 | Westerngeco L.L.C. | Interpolating seismic data |
GB2471456B (en) * | 2009-06-29 | 2012-06-20 | Geco Technology Bv | Interpolation and/or extrapolation of seismic data |
US8612157B2 (en) * | 2010-07-08 | 2013-12-17 | Westerngeco L.L.C. | Method to attenuate strong marine seismic noise |
CN103941290B (zh) * | 2013-01-22 | 2017-04-05 | 中国石油集团东方地球物理勘探有限责任公司 | 一种海底电缆运动轨迹模拟方法及*** |
US11598893B2 (en) | 2018-01-10 | 2023-03-07 | Landmark Graphics Corporation | Seismic rock property prediction in forward time based on 4D seismic analysis |
WO2020247503A1 (fr) * | 2019-06-03 | 2020-12-10 | Ion Geophysical Corporation | Acquisition de données sismiques éparses |
CN111551988B (zh) * | 2020-04-23 | 2021-06-25 | 中国地质大学(武汉) | 一种深度学习与预测滤波结合的地震数据反假频插值方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4594693A (en) | 1983-11-04 | 1986-06-10 | Mobil Oil Corporation | Seismic trace interpolation using f-k filtering |
US4866659A (en) * | 1984-04-06 | 1989-09-12 | Pennzoil Company | Method for selection of mining and drilling sites using synthesized three dimensional seismic data |
US4672545A (en) * | 1984-04-06 | 1987-06-09 | Pennzoil Company | Method and apparatus for synthesizing three dimensional seismic data |
US4887244A (en) | 1988-06-28 | 1989-12-12 | Mobil Oil Corporation | Method for seismic trace interpolation using a forward and backward application of wave equation datuming |
US5235556A (en) * | 1992-01-10 | 1993-08-10 | Halliburton Geophysical Services Inc. | Interpolation of aliased seismic traces |
US5596548A (en) | 1994-05-12 | 1997-01-21 | Exxon Production Research Company | Seismic imaging using wave equation extrapolation |
US5617372A (en) * | 1996-08-14 | 1997-04-01 | Western Atlas International, Inc. | Unaliased spatial trace interpolation in the f-k domain |
US5677892A (en) * | 1996-08-14 | 1997-10-14 | Western Atlas International, Inc. | Unaliased spatial trace interpolation in the f-k domain |
US5850622A (en) * | 1996-11-08 | 1998-12-15 | Amoco Corporation | Time-frequency processing and analysis of seismic data using very short-time fourier transforms |
US6021379A (en) * | 1997-07-29 | 2000-02-01 | Exxon Production Research Company | Method for reconstructing seismic wavefields |
-
2000
- 2000-06-15 CA CA002311822A patent/CA2311822A1/fr not_active Abandoned
- 2000-06-15 AU AU40870/00A patent/AU777935B2/en not_active Ceased
- 2000-06-20 FR FR0007855A patent/FR2795181A1/fr not_active Withdrawn
- 2000-06-20 GB GB0015125A patent/GB2353096B/en not_active Expired - Fee Related
-
2001
- 2001-01-03 US US09/754,135 patent/US6292755B2/en not_active Expired - Lifetime
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008117121A2 (fr) * | 2006-11-06 | 2008-10-02 | Magnitude Spas | Système, procédé et produit de programme d'ordinateur pour accumuler des données de bruit sismique pour analyser des événements sismiques |
WO2008117121A3 (fr) * | 2006-11-06 | 2009-02-05 | Magnitude Spas | Système, procédé et produit de programme d'ordinateur pour accumuler des données de bruit sismique pour analyser des événements sismiques |
GB2457614A (en) * | 2006-11-06 | 2009-08-26 | Magnitude Spas | System, method and computer program product for stacking seismic noise data to analyze seismic events |
GB2457614B (en) * | 2006-11-06 | 2011-05-25 | Magnitude Spas | System, method and computer program product for stacking seismic noise data to analyze seismic events |
US8867308B2 (en) | 2006-11-06 | 2014-10-21 | Magnitude Spas | System, method and computer program product for stacking seismic noise data to analyze seismic events |
US9678232B2 (en) | 2006-11-06 | 2017-06-13 | Baker Hughes Incorporated | System, method and computer program product for stacking seismic noise data to analyze seismic events |
Also Published As
Publication number | Publication date |
---|---|
AU777935B2 (en) | 2004-11-04 |
GB2353096A (en) | 2001-02-14 |
US6292755B2 (en) | 2001-09-18 |
GB2353096B (en) | 2003-07-16 |
GB0015125D0 (en) | 2000-08-09 |
AU4087000A (en) | 2001-01-04 |
US20010005813A1 (en) | 2001-06-28 |
CA2311822A1 (fr) | 2000-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2795181A1 (fr) | Procede de prospection sismique | |
FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
CA2786411A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
EP0313459B1 (fr) | Procédé d'acquisition de données sismiques et dispositif de mise en oeuvre du procédé | |
EP3339910B1 (fr) | Dispositif et procédé de deblending basé sur un modèle | |
EP0122824B1 (fr) | Procédé et dispositif pour l'optimisation des données sismiques | |
FR2790102A1 (fr) | Procede de traitement de jeux de donnees sismiques | |
US20060190181A1 (en) | Method for combining seismic data sets | |
FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
FR2737309A1 (fr) | Methode de reduction du bruit dans les signaux sismiques par filtrage adaptatif d'une reference du bruit | |
FR2885414A1 (fr) | Procede de deconvolution de signature de source | |
FR2989788A1 (fr) | Titre non renseigne. | |
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 | |
FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
FR2794249A1 (fr) | Procede de determination d'au moins un temps de parcours multiple sur un lever sismique | |
FR2874270A1 (fr) | Procede d'exploration sismique | |
EP0555148B1 (fr) | Procédé pour déterminer une signature en champ lointain d'une pluralité de sources sismiques | |
FR2974636A1 (fr) | Dispositif et procede d'extrapolation angulaire adaptative pour des rassemblements angulaires tridimensionnels de migration inverse dans le temps | |
FR2779532A1 (fr) | Procede d'attenuation des reverberations de colonne d'eau dans un signal sismique de capteurs doubles | |
FR2985832A1 (fr) | Dispositif et procede pour le retrait de multiples d'epoques de donnees sismiques | |
US6510390B1 (en) | 3-D seismic trace extrapolation and interpolation | |
EP0797780B1 (fr) | Methode de traitement de traces sismiques reflexion enregistrees pour des deports variables | |
FR2743896A1 (fr) | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |