FR2797501A1 - Procede de determination d'un coefficient de reflexion du fond d'une masse d'eau et procede de prospection sismique - Google Patents
Procede de determination d'un coefficient de reflexion du fond d'une masse d'eau et procede de prospection sismique Download PDFInfo
- Publication number
- FR2797501A1 FR2797501A1 FR0005710A FR0005710A FR2797501A1 FR 2797501 A1 FR2797501 A1 FR 2797501A1 FR 0005710 A FR0005710 A FR 0005710A FR 0005710 A FR0005710 A FR 0005710A FR 2797501 A1 FR2797501 A1 FR 2797501A1
- Authority
- FR
- France
- Prior art keywords
- water
- seismic
- traces
- reflection
- filtered
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 30
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 239000002689 soil Substances 0.000 claims 2
- 238000007796 conventional method Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010304 firing Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000011179 visual inspection Methods 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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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é de prospection sismique.Des temps de parcours jusqu'au fond de l'eau sont obtenus par des procédés classiques et des données sismiques de déport proche sont ensuite soumises à une déréverbération en utilisant une série de coefficients de réflexion possibles du fond (17) de la masse d'eau (15) pour donner une série de traces filtrées. Le coefficient de réflexion du fond de l'eau qui donne la valeur minimale sur la série de traces traitées est déterminé pour donner un ensemble de coefficients de réflexion choisis en chaque point d'analyse dans une fenêtre de temps. Une pondération est déterminée pour chaque temps d'analyse, cette pondération étant associée à la variation de l'amplitude de la trace filtrée en ce point d'échantillonnage.Domaine d'application : prospection sismique marine, etc.
Description
L'invention concerne l'exploration sismique marine, et plus
particulièrement la détermination des coefficients de réflexion du fond de l'eau à partir de sismogrammes enregistrés pour une utilisation dans la déréverbération de données sismiques sans endommager de façon nuisible des séquences de réflexion qui ont une périodicité proche de celle du temps d'aller et retour dans l'eau d'ondes sismiques. Dans l'exploration sismique marine, de l'énergie sismique est générée par des sources remorquées derrière un navire. Les ondes de pression acoustique résultantes sont réfléchies depuis les formations souterraines et enregistrées sous forme de sismogrammes. Habituellement, les sismogrammes enregistrés contiennent des bruits qui
masquent les caractéristiques des formations souterraines.
Ces bruits comprennent des réflexions multiples du fond de l'eau comme illustré sur la figure 2 des dessins annexés et décrits ci-après et des réverbérations comme illustré sur
les figures 3 et 4 de ces mêmes dessins.
Pendant que l'impulsion de tir initial se déplace vers le fond de l'eau, elle est partiellement transmise et partiellement réfléchie par le fond de l'eau. La partie qui est transmise est réfléchie par des discontinuités situées au-dessous de l'eau pour constituer des réflexions primaires qui sont importantes dans la formation d'images du sous-sol. La partie réfléchie revient par rebond vers la surface o elle est de nouveau réfléchie vers le fond de l'eau. Une partie de cette impulsion passe de nouveau à travers la limite et une partie est réfléchie. Le processus se poursuit indéfiniment. A chaque fois qu'une impulsion rebondit entre la surface et le fond de l'eau, une impulsion émerge de la limite située au fond pour suivre l'impulsion initiale dans la section. Ces impulsions multiples sont retardées par des intervalles de temps égaux les unes par rapport aux autres, et leurs amplitudes sont diminuées successivement par le coefficient de réflexion au fond de l'eau. Un processus similaire se produit avec les
réflexions primaires qui rebondissent dans la couche d'eau.
Le résultat net en est la production d'un sismogramme constitué d'une réflexion primaire unique suivie d'un train d'impulsions multiples ou de réverbération espacées de
façon égale.
La trace sismique qui est générée par une impulsion de tir peut être décrite mathématiquement sous la forme
(1)
X(t) = S(t) * R(t) * H(t) + N(t) o X(t) est le signal enregistré, S(t) est le signal de la source, R(t) est une séquence de coefficients de réflexion caractérisant le sous-sol, H(t) est une fonction de distorsion et N(t) est un bruit additif. Le bruit peut être
cohérent ou incohérent.
La fonction de distorsion qui caractérise des réverbérations est donnée par l'équation (2) o gR (t) = ()+ ôfr)k Jet - kT) k=i o H R(t) est l'opérateur de réverbération, r est le coefficient de réflexion au fond de l'eau, 8 est la fonction 6 de Kronecker, T est le temps de parcours d'aller et retour dans l'eau pour des ondes acoustiques. Dans le domaine de la transformée en z, l'opérateur de déréverbération peut être désigné par (3)
HD (Z) I
(+ rZ) Dans les cas o le fond de l'eau est relativement plat, la distorsion de réverbération descendante au point de tir, et la distorsion de réverbération montante au détecteur sont approximativement les mêmes. Par conséquent, pour un fond de l'eau plat, l'opérateur de déréverbération total peut être représenté par la convolution de l'opérateur en deux points au point de tir avec l'opérateur en deux points au détecteur. Ceci donne l'opérateur dit de Backus
(4)
B(z)=(1 rZ)2 (l+ rZ)' Le brevet des Etats-Unis d'Amérique N 3 689 874 de Foster et al. décrit un procédé pour la déréverbération du signal enregistré, utilisant l'autocorrélation des traces enregistrées. Le procédé suppose que l'autocorrélation du coefficient de réflexion et des séquences de bruit peut
être caractérisée comme étant un bruit blanc.
Le brevet des Etats-Unis d'Amérique N 4 146 871 de Ruehle décrit un procédé pour la déréverbération de signaux enregistrés, dans lequel l'opération de déréverbération est déterminée par l'essai de différentes valeurs pour le coefficient de réflexion r du fond de l'eau et le temps de parcours en aller et retour dans l'eau T et, après que la meilleure valeur a été déterminée d'après une inspection visuelle des traces traitées, les canons à air de la source sismique sont déclenchés séquentiellement pour simuler
l'opérateur de déréverbération.
Le déploiement de capteurs sur le fond océanique est souvent utilisé pour obtenir des données sismiques en trois dimensions dans une zone présentant des obstacles qui s'opposent à l'utilisation de rubans remorqués. Une telle situation apparaît, par exemple, lors de l'acquisition de données à proximité de plates-formes de production. Avec des capteurs sur le fond océanique, le "fantôme" affectant les données de capteurs enregistrées dû à une réflexion depuis la surface de l'eau apparaît davantage qu'avec des données de rubans, les rubans étant habituellement à des profondeurs inférieures à 10 mètres tandis que les capteurs du fond océanique peuvent être à des profondeurs de plusieurs dizaines ou centaines de mètres d'eau. Les dispositifs de capteurs multiples comprennent des capteurs de pression (hydrophones) et des capteurs de composantes verticales. Ces derniers sont habituellement des capteurs de vitesse (géophones). Si l'on connaît le coefficient de réflexion du fond de l'eau, les signaux de pression et de vitesse peuvent être combinés pour donner le champ arrivant en montant. Le brevet des Etats-Unis d'Amérique N 5 524 100 délivré à Paffenholz décrit un agencement dans lequel des dispositifs de capteurs multiples sont déployés sur le fond de la mer. Des signaux sismiques de pression et de vitesse sont combinés, le signal combiné est transformé en le domaine de fréquence et multiplié par l'opérateur de Backus inverse, ou bien le signal combiné est soumis à une convolution avec l'opérateur de Backus inverse, et un algorithme d'optimisation est utilisé pour résoudre le coefficient de réflexion du fond de l'eau. Les signaux sismiques de pression et de vitesse sont combinés, et le signal combiné est multiplié par l'opérateur de Backus inverse contenant le coefficient de réflexion du fond de l'eau pour éliminer des multiples en baïonnette du premier ordre. L'un des problèmes associés aux procédés de l'art antérieur est que la détermination de l'opérateur de déréverbération peut être troublée par des séquences de coefficient de réflexion géologiques qui ont un temps de retard égal au temps d'aller et retour dans l'eau. Les procédés de traitement sont incapables de faire la distinction entre des encoches spectrales provoquées par des réverbérations et des encoches spectrales provoquées par des réflexions souterraines. Par conséquent, l'opérateur de déréverbération supprime également des séquences de réflexion qui ont le même temps de retard que l'opérateur de Backus et les données sismiques traitées résultantes sont incorrectes par le fait que certaines réflexions associées à la géologie souterraine sont supprimées, ce qui conduit à des mauvaises interprétations possibles des données. Il serait souhaitable de disposer d'un procédé pour la déréverbération de données sismiques marines qui soit robuste en présence d'une telle
interférence géologique. L'invention satisfait à ce besoin.
L'invention concerne un procédé de détermination du coefficient de réflexion du fond de l'eau, qui est insensible à une interférence géologique, définie ici sous la forme de réflecteurs souterrains ayant des temps de retard proches de ceux dans la réverbération du fond de l'eau. Les temps de parcours de fond de l'eau sont obtenus par les procédés classiques. Des données sismiques de déports proches sont ensuite soumises à une déréverbération pour une série de coefficients de réflexion possible du
fond de l'eau afin de donner une série de traces filtrées.
Le coefficient de réflexion du fond de l'eau qui procure la valeur minimale sur la série de traces filtrées est déterminé pour donner un ensemble de coefficient de réflexion choisis pour chacun des points d'échantillonnage dans une fenêtre de temps. Une pondération est déterminée pour chaque temps dans la fenêtre de temps, la pondération étant liée à la variation de l'amplitude de la trace filtrée en ce point d'échantillonnage sur la série de coefficients de réflexion du fond de l'eau. Les coefficients de réflexion choisis sont pondérés par les pondérations et soumis à une moyenne pour réaliser une
détermination du coefficient de réflexion du fond de l'eau.
Ce coefficient de réflexion déterminée du fond de l'eau est utilisé pour le traitement de la trace sismique
enregistrée.
L'invention sera décrite plus en détail en regard des dessins annexés à titre d'exemple nullement limitatif et sur lesquels: la figure 1 illustre schématiquement un exemple de système d'acquisition de données sismiques destiné à être utilisé dans un environnement marin; la figure 2 est un diagramme illustrant des réflexions multiples du fond de l'eau; les figures 3 et 4 sont des diagrammes illustrant des réverbérations dans l'eau; les figures 5 à 8 sont des diagrammes illustrant des opérateurs inverses dans le domaine temporel à deux points et trois points; et les figures 9A et 9B illustrent des étapes de la présente invention pour la détermination du coefficient de
réflexion du fond de l'eau.
En référence à présent à la figure 1, celle-ci illustre un exemple de parties d'un système d'acquisition de données sismiques marines. Un navire 10 flottant sur une masse d'eau 15 s'étendant au-dessus de la terre 16 a déployé derrière lui un dispositif 20 de sources sismiques et un câble 25 de ruban. Le dispositif 20 de sources sismiques est habituellement constitué de canons à air individuels 20a, 20b...20n qui sont déclenchés sous la commande d'un dispositif de commande (non représenté) embarqué à bord du navire 10. Des impulsions sismiques se propagent dans le sol et sont réfléchies par un réflecteur 22 s'y trouvant. Pour simplifier l'illustration, un seul réflecteur est représenté; en réalité, il peut y avoir de nombreux réflecteurs, donnant chacun naissance à une impulsion réfléchie. Après la réflexion, ces impulsions reviennent à la surface o elles sont enregistrées par des détecteurs (hydrophones) 30a, 30b,...30n dans le câble de ruban. La profondeur du dispositif de sources et du câble de ruban est commandée par des dispositifs auxiliaires (non représentés). La figure 2 montre des réflexions multiples du fond de l'eau se déplaçant dans la masse d'eau, qui seraient également détectées par les détecteurs 30a, 30b,...30n. Les figures 3 et 4 montrent des réverbérations d'impulsions revenant par réflexion du sol, provoquées par des
réflexions multiples dans la couche d'eau.
La trace sismique qui est générée par une impulsion de
tir peut être décrite mathématiquement par l'équation (1).
La fonction de distorsion introduite par une réflexion unique dans la couche d'eau est illustrée par les figures 5 et 6. L'impulsion initiale à un temps 0 est suivie d'une impulsion d'amplitude r à un temps Tw, o Tw est le temps d'aller et retour dans l'eau pour une impulsion sismique et r est le coefficient de réflexion du fond de l'eau. Ce coefficient de réflexion du fond de l'eau peut être positif
ou négatif, comme montré sur les figures 5 et 6.
Ce coefficient de réflexion du fond de l'eau combiné avec des réverbérations multiples dans la couche d'eau a pour effet d'introduire une fonction de réverbération. Dans des conditions dans lesquelles le fond de l'eau est relativement plat, la distorsion de réverbération dans le sens descendant, à la source, et la distorsion de réverbération, dans le sens montant, au détecteur sont approximativement les mêmes. Par conséquent, pour un fond d'eau plat, un opérateur de déréverbération total peut être représenté par la convolution de l'opérateur à deux points
à la source avec l'opérateur à deux points au détecteur.
Ceci donne l'opérateur dit de Backus indiqué par l'équation (4). Le temps d'aller et retour dans l'eau peut être aisément obtenu par des procédés antérieurs. L'un des procédés les plus précis consiste à le calculer par bathymétrie, car la vitesse du son dans l'eau est bien connue. Un autre procédé consiste à le déterminer d'après le temps d'arrivée de la réflexion du fond de l'eau au premier détecteur 30a et à le calculer à partir de valeurs connues de la profondeur du détecteur et de la profondeur
de la source au-dessous de la surface 18 de l'eau.
Cependant, la détermination du coefficient de réflexion du fond de l'eau est problématique. Comme indiqué précédemment, les procédés basés sur les spectres ou sur l'autocorrélation de la trace sismique enregistrée sont sujets à une erreur due à des séquences géologiques avec la même période de temps que le temps d'aller et retour dans l'eau. Une détermination directe à partir de l'amplitude de la réflexion du fond de l'eau est difficile à réaliser car, en eau peu profonde, la réflexion du fond de l'eau peut arriver dans la même fenêtre de temps que celle d'une
arrivée directe à travers l'eau, de la source au détecteur.
Les figures 9A et 9B auxquelles on se réfère à présent illustrent un organigramme de la présente invention pour la
détermination du coefficient de réflexion du fond de l'eau.
Le procédé commence avec un temps d'aller et retour Tw dans l'eau 100 et une trace sismique enregistrée 102. La trace enregistrée comprend des impulsions de réflexion provenant de réflecteurs situés au-dessous du fond de l'eau et comprend aussi des réverbérations indésirables qui doivent être atténuées. Ce temps d'aller et retour est double du temps de parcours de la surface de l'eau jusqu'au fond de l'eau. Une plage de valeurs possibles pour le coefficient de réflexion r du fond de l'eau est choisie. Dans une forme appréciée de réalisation de l'invention, la plage de coefficient de réflexion possible du fond de l'eau est de -0,9 à +0,9. Une valeur initiale du coefficient de réflexion est choisie en 104 à la limite inférieure de la plage des coefficients de réflexion. En utilisant cette valeur initiale du coefficient de réflexion et le temps d'aller et retour dans l'eau, on applique un opérateur de Backus donné par l'équation 4 à la trace enregistrée, en 106. Une trace filtrée est sauvegardée en 112 en tant que partie d'une série de trace filtrée. La valeur possible du coefficient de réflexion est augmentée d'une quantité Ar 108 et un contrôle est réalisé pour déterminer si la valeur possible dépasse la limite supérieure de la plage de coefficients de réflexion possibles 110. Dans une forme appréciée de réalisation de l'invention, la valeur de Ar est de 0,01. Si le test effectué en 110 est négatif le processus revient à l'étape 104 avec la nouvelle valeur du coefficient de réflexion possible et les étapes 106 et 110 sont répétées, en sauvegardant à chaque fois une nouvelle trace après déréverbération en tant que partie d'une série de traces 112 ayant subi une déréverbération. Si le test effectué en 110 est positif la boucle d'itération allant de 106 à 110 est achevée et la série de traces filtrées en 114 est utilisée dans des étapes suivantes illustrées sur la
figure 9B.
Une fenêtre de temps pour l'analyse des données est sélectionnée. Comme cela est connu des spécialistes de la technique, le point de départ de la fenêtre ne doit pas contenir de bruits initiaux portés dans l'eau et d'autres artéfacts. Au temps de départ de la fenêtre 116, la série entière de traces ayant subi une déréverbération est examinée pour trouver les amplitudes minimale et maximale de la série de traces filtrées. Les minimum et maximum sont délivrés en sortie à une table mémorisée de valeurs en même temps que le coefficient de réflexion possible correspondant à l'amplitude minimale à l'absolu. C'est la valeur du coefficient de réflexion qui, lorsqu'elle est utilisée dans un opérateur de déréverbération, présente la réduction la
plus importante d'amplitude de la trace.
Le temps pour l'analyse est incrémenté d'une quantité At en 120. Il est vérifié si le temps incrémenté dépasse le temps maximal de la fenêtre d'analyse en 122. Ce temps maximal de la fenêtre d'analyse peut être choisi sur la base de l'un de plusieurs critères. Dans une forme de réalisation de l'invention, le temps maximal de la fenêtre d'analyse est basé sur un rapport signal/bruit approprié; il est évident qu'il est indésirable de sélectionner un opérateur de déréverbération provenant de bruits. En variante, le temps maximal d'analyse est basé sur des considérations géologiques; dans la plupart des zones d'exploration auxquelles on s'intéresse, on connaît à priori un temps de parcours approximatif pour l'intervalle géologique auquel on s'intéresse. L'incrément du temps d'analyse At est choisi de façon à être un multiple de l'intervalle de temps d'échantillon At utilisé pour l'enregistrement des données. Dans une forme appréciée de réalisation de l'invention, ils sont égaux, c'est-à-dire que l'incrément de temps d'analyse est égal à l'incrément
de temps d'échantillonnage des données sismiques.
Si le test effectué en 122 est négatif, le temps d'analyse est incrémenté et les étapes 116 à 122 sont répétées. A chaque temps d'analyse, les amplitudes absolues maximale et minimale de la série de traces ayant subi une déréverbération et le coefficient de réflexion possible correspondant à l'amplitude absolue minimale des traces ayant subi une déréverbération sont sauvegardés. Une fois que le temps d'analyse maximal est dépassé, comme indiqué par une réponse positive en 122, on procède à une analyse
des valeurs mémorisées, en 124.
L'analyse des valeurs mémorisées est basée sur l'analyse raisonnée suivante: chacune des valeurs mémorisées du coefficient de réflexion possible (une pour chaque temps d'analyse) est la valeur qui réduit à un minimum l'amplitude de la trace ayant été soumise à une déréverbération. C'est ce qu'un opérateur de déréverbération devrait faire s'il n'y avait aucune impulsion sismique réfléchie au temps de l'analyse. Cependant, s'il y a une impulsion de réflexion sismique au temps de l'analyse, une
réduction à une valeur minimale ne convient pas.
L'établissement d'une moyenne sur une plage de temps d'analyse est une moyenne établie sur différentes impulsions géologiques réfléchies et réduirait donc l'effet
de la géologie sur le coefficient de réflexion déterminé.
Cependant, au lieu d'un simple calcul de moyenne, la présente invention pondère le coefficient de réflexion choisi par un facteur dépendant de la plage de valeurs de coefficient de réflexion de la série de traces ayant étant soumises à une déréverbération: ceci signifie qu'une pondération plus grande est donnée à un coefficient de réflexion choisi
ayant un effet plus grand.
Les spécialistes de la technique reconnaîtront que d'autres fonctions de pondération pourraient également être utilisées. Ces variantes sont connues des spécialistes de la technique et entendent être incluses dans le cadre de l'invention. Le brevet N 5 524 100 précité, dont le contenu est incorporé ici totalement à titre de référence, décrit une invention dans laquelle un dispositif de capteurs doubles déployés sur le fond de l'océan est utilisé pour déterminer le coefficient de réflexion du fond de l'eau. Le dispositif de capteurs doubles comporte un géophone pour composante verticale et un hydrophone. Le coefficient de réflexion est déterminé en minimisant la puissance dans une bande spectrale. Le coefficient de réflexion déterminé est utilisé pour déterminer un champ d'onde arrivant en montant qui est exempt de réverbérations. Comme indiqué précédemment, la
géologie contamine le coefficient de réflexion déterminé.
Le coefficient de réflexion du fond de l'eau déterminé dans la présente invention et décrit ci-dessus en référence aux figures 9A et 9B peut être utilisé conjointement avec la divulgation du brevet N 5 524 100 précité pour déterminer un champ d'onde arrivant en montant, ayant été soumis à une déréverbération, qui n'est pas sensible à une contamination géologique de la détermination du coefficient
de réflexion.
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 (16)
1. Procédé de détermination d'un coefficient de réflexion du fond (17) d'une masse d'eau (15), caractérisé en ce qu'il comprend: (a) le déclenchement d'au moins une source sismique (20) dans la masse d'eau (15) pour générer des impulsions sismiques qui se propagent dans la masse d'eau et dans le sol (16) situé au-dessous d'elle; (b) l'enregistrement d'au moins une trace sismique représentative d'impulsions sismiques réfléchies par des réflecteurs sismiques (22) dans le sol, cette trace sismique comprenant des réverbérations dans la masse d'eau des impulsions sismiques réfléchies par les réflecteurs sismiques; (c) l'obtention d'un temps de parcours d'une surface supérieure (18) de la masse d'eau jusqu'au fond de la masse d'eau; (d) pour chacun d'une première pluralité de coefficients de réflexion possibles du fond de la masse d'eau, l'obtention d'une trace filtrée à partir du temps de parcours, de la trace sismique ou des traces sismiques, et du coefficient de réflexion possible pour donner une première pluralité de traces filtrées; (e) pour chacun d'une seconde pluralité de temps d'analyse dans une fenêtre de temps, la détermination d'un facteur de qualité à partir de valeurs de la première pluralité de traces filtrées à chaque temps d'analyse pour donner une second pluralité de facteurs de qualité; et (f) la détermination du coefficient de réflexion au fond de la masse d'eau à partir de la première pluralité de traces filtrées et de la seconde pluralité de facteurs
de qualité.
2. Procédé selon la revendication 1, caractérisé en ce que la ou chaque source sismique comporte en outre plusieurs sources sismiques (20a...20n) dans un dispositif (20).
3. Procédé selon la revendication 1, caractérisé en ce que l'enregistrement de la ou de chaque trace sismique est effectué en utilisant au moins l'un (i) d'un hydrophone dans l'eau, (ii) d'un hydrophone au fond de la masse d'eau
et (iii) d'un géophone au fond de la masse d'eau.
4. Procédé selon la revendication 1, caractérisé en ce que chacun de la première pluralité de coefficients de
réflexion possibles est compris entre -0,9 et 0,9.
5. Procédé selon la revendication 1, caractérisé en ce que la première pluralité de traces ayant été soumises à une déréverbération est obtenue à partir de la ou de chaque trace sismique par un opérateur de déréverbération de la forme 1 + 2R6(T)+2 R26(2T) o R est l'un de la première pluralité de coefficients de réflexion possibles, T est le temps de parcours obtenu et 6 est la fonction delta de Kronecker.
6. Procédé selon la revendication 1, caractérisé en ce que chaque facteur de qualité est une valeur absolue de la différence entre une valeur maximale de la première pluralité de traces filtrées au temps d'échantillonnage et une valeur absolue minimale de la première pluralité de
traces filtrées au temps d'analyse.
7. Procédé selon la revendication 1, caractérisé en ce que la détermination du coefficient de réflexion comprend en outre: (i) la détermination d'un maximum des traces filtrées en chacun de la seconde pluralité de temps d'analyse pour donner une seconde pluralité de coefficients de réflexion choisis; (ii) la détermination d'une seconde pluralité de coefficients de pondération pour la seconde pluralité de coefficients de réflexion choisis, chaque coefficient de pondération étant associé à une plage de valeurs de la première pluralité de traces filtrées à chacun de la seconde pluralité de temps d'analyse; et (iii) la détermination du coefficient de réflexion sous la forme d'une moyenne de la seconde pluralité de coefficients de réflexion choisis pondérés par le coefficient de pondération correspondant de la
seconde pluralité de coefficients de pondération.
8. Procédé selon la revendication 7, caractérisé en ce que chacun de la seconde pluralité de coefficients de pondération est une différence entre un maximum et un minimum de la première pluralité de traces ayant été soumises à une déréverbération en chacun de la seconde
pluralité de temps d'analyse.
9. Procédé de prospection sismique, caractérisé en ce qu'il comprend: (a) le déclenchement d'au moins une source sismique (20a...20n) dans une masse d'eau (15) pour générer des impulsions sismiques qui se propagent dans la masse d'eau et dans le sol (16) situé au-dessous d'elle; (b) l'enregistrement d'au moins une trace sismique représentative d'impulsions sismiques réfléchies par des réflecteurs sismiques (22) dans le sol, la ou chaque trace sismique comprenant des réverbérations d'impulsions sismiques réfléchies par les réflecteurs sismiques; (c) l'obtention d'un temps de parcours depuis une surface supérieure (18) de la masse d'eau jusqu'à un fond (17) de la masse d'eau; (d) pour chacun d'une première pluralité de coefficients de réflexion possibles du fond de la masse d'eau, l'obtention d'une trace filtrée à partir du temps de parcours, de la ou de chaque trace sismique et du coefficient de réflexion possible pour donner une première pluralité de traces filtrées; (e) pour chacun d'une seconde pluralité de temps d'analyse dans une fenêtre de temps, la détermination d'un facteur de qualité à partir de valeurs de la première pluralité de traces filtrées à chaque temps d'analyse pour donner une seconde pluralité de facteurs de qualité; (f) la détermination d'un coefficient de réflexion au fond de la masse d'eau à partir de la première pluralité de traces filtrées et de la seconde pluralité de facteurs de qualité; et (g) le traitement de la ou de chaque trace sismique en utilisant le coefficient de réflexion déterminé pour donner une trace sismique traitée dans laquelle les
réverbérations ont été réduites.
10. Procédé selon la revendication 9, caractérisé en ce que la ou chaque source sismique comprend en outre plusieurs sources sismiques (20a...20n) dans un dispositif (20).
11. Procédé selon la revendication 9, caractérisé en ce que l'enregistrement de la ou de chaque trace sismique est effectué en utilisant au moins l'un (i) d'un hydrophone dans l'eau, (ii) d'un hydrophone au fond de la masse d'eau,
et (iii) d'un géophone au fond de la masse d'eau.
12. Procédé selon la revendication 9, caractérisé en ce que chacun de la seconde pluralité de coefficients de
réflexion possibles est compris entre -0,9 et 0,9.
13. Procédé selon la revendication 9, caractérisé en ce que la première pluralité de traces ayant été soumises à une déréverbération est obtenue à partir de la ou chaque trace sismique par un opérateur de déréverbération de la forme 1 + 2R6(T)+2 R26(2T) o R est l'un de la première pluralité de coefficients de réflexion possibles, T est le temps de parcours obtenu et 6 est la fonction delta de Kronecker.
14. Procédé selon la revendication 9, caractérisé en ce que chacun des facteurs de qualité est une différence entre une valeur maximale absolue de la première pluralité de traces filtrées au temps d'échantillonnage et d'une valeur minimale absolue de la première pluralité de traces
filtrées au temps d'analyse.
15. Procédé selon la revendication 1, caractérisé en ce que la détermination du coefficient de réflexion comprend en outre: (i) la détermination d'un maximum des traces filtrées en chacun de la seconde pluralité de temps d'analyse pour donner une seconde pluralité de coefficients de réflexion choisis; (ii) la détermination d'une seconde pluralité de coefficients de pondération pour la seconde pluralité de coefficients de réflexion choisis, chaque coefficient de pondération étant associé à une plage de valeurs de la première pluralité de traces filtrées à chacun de la seconde pluralité de temps d'analyse; et (iii) la détermination du coefficient de réflexion sous la forme d'une moyenne de la seconde pluralité de coefficients de réflexion choisis pondérés par le coefficient de pondération correspondant de la
seconde pluralité de coefficients de pondération.
16. Procédé selon la revendication 15, caractérisé en ce que le coefficient de pondération est la différence entre un maximum et un minimum de la première pluralité de traces ayant été soumises à une déréverbération à chacun
des temps d'échantillonnage.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/306,122 US6141623A (en) | 1999-05-06 | 1999-05-06 | Method of deriving water bottom reflectivities in the presence of geologic interference |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| FR2797501A1 true FR2797501A1 (fr) | 2001-02-16 |
Family
ID=23183927
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0005710A Withdrawn FR2797501A1 (fr) | 1999-05-06 | 2000-05-04 | Procede de determination d'un coefficient de reflexion du fond d'une masse d'eau et procede de prospection sismique |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6141623A (fr) |
| AU (1) | AU765956B2 (fr) |
| CA (1) | CA2307467A1 (fr) |
| FR (1) | FR2797501A1 (fr) |
| GB (1) | GB2349697B (fr) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7359283B2 (en) * | 2004-03-03 | 2008-04-15 | Pgs Americas, Inc. | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers |
| US7336561B2 (en) * | 2004-09-07 | 2008-02-26 | Pgs Americas, Inc. | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors |
| EP2541282A1 (fr) * | 2011-06-29 | 2013-01-02 | Sercel | Procédé et dispositif pour obtenir une distance nýud/surface dans un réseau de nýuds acoustiques, produit de programme informatique correspondant et moyen de stockage |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4146871A (en) * | 1975-02-28 | 1979-03-27 | Mobil Oil Corporation | Determination of reflectivity and water depth for firing marine sources as an inverse distortion operator |
| US4234938A (en) * | 1977-10-25 | 1980-11-18 | Mobil Oil Corporation | Determination of reflectivity and water depth for marine seismic exploration |
| US5396472A (en) * | 1993-09-24 | 1995-03-07 | Western Atlas International | Method for deriving water bottom reflectivity in dual sensor seismic surveys |
| US5963507A (en) * | 1998-05-13 | 1999-10-05 | Western Atlas International, Inc. | Method for improving the accuracy of ocean bottom reflectivity estimations using the inverse backus filter |
| US5991238A (en) * | 1998-06-09 | 1999-11-23 | Western Atlas International, Inc. | Weighted backus filter method of combining dual sensor traces |
-
1999
- 1999-05-06 US US09/306,122 patent/US6141623A/en not_active Expired - Fee Related
-
2000
- 2000-04-12 AU AU27730/00A patent/AU765956B2/en not_active Ceased
- 2000-05-03 CA CA002307467A patent/CA2307467A1/fr not_active Abandoned
- 2000-05-04 FR FR0005710A patent/FR2797501A1/fr not_active Withdrawn
- 2000-05-04 GB GB0010829A patent/GB2349697B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| GB0010829D0 (en) | 2000-06-28 |
| AU765956B2 (en) | 2003-10-02 |
| US6141623A (en) | 2000-10-31 |
| GB2349697B (en) | 2003-09-24 |
| GB2349697A (en) | 2000-11-08 |
| CA2307467A1 (fr) | 2000-11-06 |
| AU2773000A (en) | 2000-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0680616B1 (fr) | Procede de determination de la reflexion d'un fond d'un corps d'eau dans des etudes sismiques realisees avec deux detecteurs | |
| US7551515B2 (en) | Source signature deconvolution method | |
| FR2737309A1 (fr) | Methode de reduction du bruit dans les signaux sismiques par filtrage adaptatif d'une reference du bruit | |
| FR2534696A1 (fr) | Procede et dispositif pour l'etude indirecte des ondes de cisaillement | |
| FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
| Dragoset et al. | Ocean-bottom cable dual-sensor scaling | |
| EP2174166B1 (fr) | Procédé et système pour estimer des paramètres de propagation physique associés à une étude sismique | |
| FR2489530A1 (fr) | Dispositif d'exploration sismique utilisant des detecteurs de pression et de vitesse | |
| FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
| FR2955397A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
| FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
| FR2779532A1 (fr) | Procede d'attenuation des reverberations de colonne d'eau dans un signal sismique de capteurs doubles | |
| JPH08501395A (ja) | 反響抑制のための改良された方法 | |
| EP0555148B1 (fr) | Procédé pour déterminer une signature en champ lointain d'une pluralité de sources sismiques | |
| DK2326972T3 (en) | DIMENSION OF SEISMIC INTERFERENCE NOISE USING A DUAL SENSOR RECORDING SYSTEM | |
| FR2978252A1 (fr) | Procede et dispositif pour la separation de champs d'ondes dans des donnees sismiques | |
| FR2690528A1 (fr) | Méthode de correction de la vitesse de propagation du son dans l'eau dans les explorations sismiques en mer. | |
| FR2874270A1 (fr) | Procede d'exploration sismique | |
| EP0786670B1 (fr) | Procédé de traitement de calibration d'une paire de capteurs hydrophone/géophone et procédé de prospection sismique mettant en oeuvre ce traitement | |
| FR2795181A1 (fr) | Procede de prospection sismique | |
| CA2195497C (fr) | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone | |
| FR2797501A1 (fr) | Procede de determination d'un coefficient de reflexion du fond d'une masse d'eau et procede de prospection sismique | |
| CA1270937A (fr) | Procede et dispositif pour la determination de la signature d'emission a grande distance d'un ensemble d'emission sismique | |
| FR3006451A1 (fr) | Procedes et dispositifs d'acquisition de donnees sismiques relatives a une zone du sous-sol situee sous la mer. | |
| FR3000224A1 (fr) | Method and device to deghost seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |